PDE 边界问题

asaki

PDE很重要 但是自己不懂  看了几天晕乎乎的 似懂非懂
看了一个文章 conference的  给出了PDE方程 但是这该怎么写入comsol还是不是很懂

忘大人能指教 各个方程如何处理 coefficient模式还是 ？

antonylau

写成标准格式

和femlab里面的方程对应一下 （u，v）是二维，Uxx这些都有对应

asaki

恩 是的 第一个方程 eps作为系数吗？如果作为洗漱  那么da项不是有参数u吗？
这怎么理解呢？

asaki

顶上去

asaki

顶上去

第一个方程中的da*diff（eps, t）这个难道作为常数项吗？

kinggun

第一个方程有点怪，没看明白。

asaki

我也就是第一个方程 diff（eps，t）这个没参量U啊

kinggun

我主要是没明白正负eps的含义。

asaki

本来是2项  分别是diff（eps，t）和 diff(sigma,t)

同时有定义 eps+sigma=1

所以diff（sigma，t）=-diff（eps，t）

六个方程应该是求解6个未知数
ug, vg, us, vs, p, and ε,
分别是这些

我初步是想用
da*diff(u,t)+∇*Γ=F general form来写
不过一个方程不该是只有一个未知数 其他作为已知吗？
所以我不太懂如何搞

asaki

对于第一项da=1/-1
Γ=0
F=eps*(diff(ug,x)+diff(vg,y)）+ug*diff(eps,x)+vg*diff(eps,y)

知道这样行吗？
eps作为未知数  ug vg作为已知？
还望kinggun版主解答

asaki

kinggun 版主 帮忙看下啊

edodo

本来是2项  分别是diff（eps，t）和 diff(sigma,t)

同时有定义 eps+sigma=1
所以diff（sigma，t）=-diff（eps，t）


还是写成两个方程然后耦合在一起的好
其实可以认为一个方程但两个变量然后写成张量形式就可以，最后在表达式里定义 eps=1-sigma就可以了

附件中的模型希望对你有帮助

asaki

具体是怎么写呢？
比如第一个
参数为epson
da=-1
F=epson*(diff(ug,x)+diff(vg,y))+ug*diff(epson,x)+vg*diff(epson,y)
第二个
参数为epson？
da=1
F=sigma*(diff(us,x)+diff(vs,y))-us*diff(epson,x)-vs*diff(epson,y)

下面几个方程也类似
只是多了Γ一项 分为x,y两个分量

不知道这么写对与错呢？

还有我加入6个pde方程发现边界条件自己不知道怎么设置

这应该是速度入口，压力出口的边界 不知道如何搞 请教大家

edodo

>还有我加入6个pde方程发现边界条件自己不知道怎么设置
>这应该是速度入口，压力出口的边界 不知道如何搞 请教大家
最简单的办法，点开一个N-S方程或者Dacy方程，设置一条边为进口速度，一条为出口压力，然后再Physics>Equation>Boundary settings找到对应的方程设置就可以了。

asaki

>还有我加入6个pde方程发现边界条件自己不知道怎么设置
>这应该是速度入口，压力出口的边界 不知道如何搞 请教大家
最简单的办法，点开一个N-S方程或者Dacy方程，设置一条边为进口速度，一条为出口压力，然后再Phy ...
edodo 发表于 2009-2-26 22:41

谢谢edodo大侠的解答
还有问题是我那样写pde方程对吗？
包括以下几个方程只要是微分都写成diff形式

或者直接写成ugx  ugy形式？
还望解答！

asaki

顶上去 版主帮忙看看!

asaki

版主帮看看

edodo

希望下面这个图对你有帮助

asaki

我再试试

asaki

我再试试