| 查看: 751 | 回复: 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 ] |
回复此楼» 猜你喜欢
梦想
已经有5人回复
-
中青基了要发朋友圈吗？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
30