| 查看: 790 | 回复: 2 | |||
[求助]
求助属性，这样写udf可以用来求二阶导数嘛？ 已有1人参与
|
|
我想利用下面的代码来求C_VOF的二阶导数，然后添加一个动量源项，可以编译成功，但是为什么一开始计算就会出现divergence detected in AMG solver: x-momentum 的错误啊 #define domain_ID 1 DEFINE_ADJUST(adjust_gradient,domain) { Thread *t; cell_t c; face_t f; domain=Get_Domain(domain_ID); thread_loop_c(t,domain) { begin_c_loop(c,t) { C_UDSI(c,t,0)=C_VOF(c,t); } end_c_loop(c,t) } } DEFINE_ADJUST(adjust_gradient,domain) { Thread *t; cell_t c; face_t f; domain=Get_Domain(domain_ID); thread_loop_c(t,domain) { begin_c_loop(c,t) { C_UDSI(c,t,1)=C_UDSI(c,t,0)[0]; } end_c_loop(c,t) } } DEFINE_ADJUST(adjust_gradient,domain) { Thread *t; cell_t c; face_t f; domain=Get_Domain(domain_ID); thread_loop_c(t,domain) { begin_c_loop(c,t) { C_UDSI(c,t,2)=C_UDSI(c,t,0)[1]; } end_c_loop(c,t) } } 这样的话C_UDSI_G（c,t,1）[0]和C_UDSI_G(c,t,2)[1]是不是就是二阶导数了？ [ Last edited by japus1988 on 2012-12-10 at 20:10 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- fluent能量源项中温度T对时间的t的偏导数，即dT/dt如何用udf编写能量源项
已经有6人回复
- 求助：如何求已得曲线某点二阶导数
已经有5人回复
- 【求助】各位大侠救我，请教一个简单的UDF程序怎么写，有金币送！
已经有10人回复
- 【求助】物料属性可以用UDF修改吗
已经有8人回复
- 【求助】二阶导数 UV
已经有3人回复
- 【求助】怎样在origin里把一个曲线求二阶导数？
已经有2人回复
2楼2012-12-11 00:13:41
3楼2018-05-02 17:04:37
