| 查看: 575 | 回复: 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 ] |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有7人回复
-
请教限项目规定
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有15人回复
-
Materials Today Chemistry审稿周期
已经有6人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
推荐一本书
已经有13人回复
-
国自然申请面上模板最新2026版出了吗？
已经有17人回复
-
纳米粒子粒径的测量
已经有8人回复
-
疑惑？
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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