| 查看: 570 | 回复: 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 ] |
回复此楼» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有21人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- fluent能量源项中温度T对时间的t的偏导数，即dT/dt如何用udf编写能量源项
已经有6人回复
- 求助：如何求已得曲线某点二阶导数
已经有5人回复
- 【求助】各位大侠救我，请教一个简单的UDF程序怎么写，有金币送！
已经有10人回复
- 【求助】物料属性可以用UDF修改吗
已经有8人回复
- 【求助】二阶导数 UV
已经有3人回复
- 【求助】怎样在origin里把一个曲线求二阶导数？
已经有2人回复