关于边界条件不连续的设置

liuyimi0732 · 发表于 2008-11-14 15:00:17

本帖最后由 TBE_Legend 于 2009-6-5 22:16 编辑

边界是一个突变边界，我想在comsol里设置u1-u2=0.3，请问这里边界条件该如何在软件的boundary setting里设置？

[ 本帖最后由 liuyimi0732 于 2008-11-18 14:52 编辑 ]

liuyimi0732 · 发表于 2008-11-26 10:29:37

或者comsol是否能够实现这种不连续边界的表达？
其他地方，比如在equation system里面的boundary settings里可不可以找到解决方法

TBE_Legend · 发表于 2008-11-26 11:07:18

1）最原始的方法，用扩展耦合变量来做这种约束。
2）identify pairs边界。很方便。

kinggun · 发表于 2008-11-26 14:18:38

如果是内部边界，可以增加一个同样的应用模式，但是不同的应用模式用来求解不同的域。

liuyimi0732 · 发表于 2008-11-27 10:22:45

谢谢3L
请问这两种方法具体该如何操作呢？
或者，帮助文档中哪块有相应介绍

liuyimi0732 · 发表于 2008-12-6 08:24:28

我试过了identity pairs，它只能实现c*gradu1-c*gradu2+qu=g的边界，并不能做到u1-u2=g的设置呀

kinggun · 发表于 2008-12-6 17:54:54

Sorry，我没说清楚，我的意思把这种内部边界通过增加应用模式进行耦合调整成外部边界，然后就可以设定不同的边界条件了。
例如两个矩形域连接，只有一个模式时，中间线是内部边界。但是如果增加一个模式，一个域内用一个模式，中间线就变成了外部边界，就可以设置两个不同的边界值了。

liuyimi0732 · 发表于 2008-12-9 09:28:09

我是这样做的，两块的subdomain分别是不同的公式描述
但是交界处的边界如图1所示，不能设置
在groups的边界设置里有设置项，但是依然是纽曼边界，无法设u2-u1=g
还望进一步指教

uracratia · 发表于 2008-12-26 09:18:53

他应该是说他只约束两个边界的差值。（实际上比较怀疑大部分情况下这样的约束不够做出解来）
如果只是在边界两边设定不同的边界条件，那是很简单的，assembly就专门做这个事情吧

原帖由 kinggun 于 2008-12-6 17:54 发表
Sorry，我没说清楚，我的意思把这种内部边界通过增加应用模式进行耦合调整成外部边界，然后就可以设定不同的边界条件了。
例如两个矩形域连接，只有一个模式时，中间线是内部边界。但是如果增加一个模式，一个域内用 ...

uracratia · 发表于 2008-12-26 12:19:56

建立pair以后假设u1是源，u2是目的
那么在domain1调用u2就是dst2src_ipn(u)
constr=if(dst2src_ipn,u-dst2src_ipn(u)+g,[...非接触部分表达式])
在domain2调用u1就是src2dst_ipn(u)
constr=if(src2dst_ipn,src2dst_ipn(u)-u+g,[...非接触部分表达式])
没试过这样能不能解，但是理论上应该如此吧

liuyimi0732 · 发表于 2008-12-29 19:37:01

谢谢
根据10楼，我在pairs里的weak里的constr作了相应设置，可算出来还是边界上的连续解
还请进一步指教

uracratia · 发表于 2008-12-29 23:08:06

画法线方向的line,在其上对u作图，看看你到底设置对了没有

kinggun · 发表于 2009-1-1 10:32:29

这样的话，看来你最好是使用COMSOL的超弱解约束条件了，参考：Model library手册的A transport problem，该模型用超弱项定义了求解域内部的变量的一条不连续线。

uracratia · 发表于 2009-1-1 21:13:27

请看我的这个例子
V在中间的边界上不连续差值为5

原帖由 liuyimi0732 于 2008-12-29 19:37 发表
谢谢
根据10楼，我在pairs里的weak里的constr作了相应设置，可算出来还是边界上的连续解
还请进一步指教

liuyimi0732 · 发表于 2009-4-24 21:52:06

本帖最后由 liuyimi0732 于 2009-4-25 14:16 编辑

谢谢uracratia，高人啊
不过我的模型不大一样，能帮我看看嘛，谢谢
只要界面处差值为0.3就好

chbo76 · 发表于 2009-4-27 21:36:09

if(src2dst_ip3,V-src2dst_ip3(V)+5,0)
和src2dst_ip3
是什么意思啊？

uracratia · 发表于 2009-6-5 22:14:07

15# liuyimi0732
回头帮你看看，不知道你是不是已经解决了

uracratia · 发表于 2009-6-5 23:39:10

15# liuyimi0732
你完全设置错了啊，
1.只需要在一个边界设置即可
2.xxxx_ip'n'的这个n是要到identity pair 里面看n=几，一维的就是1,二维就是2，三维就是3....
已经改好了

yanhongtao · 发表于 2009-6-6 00:54:28

谢谢楼主分享，谢谢!

yydjane20 · 发表于 2009-6-6 09:18:27

本帖最后由 yydjane20 于 2009-6-6 09:36 编辑

19# yanhongtao
那么if(src2dst_ip3,V-src2dst_ip3(V),0) 分别代表什么意思呢？例如括号中最后一项有时是0，可有时是-u,可不可以解释一下

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[01.PDE] 关于边界条件不连续的设置

相关帖子

评分

回复 3# TBE_Legend 的帖子

本帖子中包含更多资源

评分

本帖子中包含更多资源

评分

本帖子中包含更多资源

本帖子中包含更多资源