| 查看: 673 | 回复: 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 ] |
回复此楼» 猜你喜欢
青C资助名额大幅增加！
已经有11人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有6人回复
-
有带发论文的吗
已经有3人回复
-
导师各种操作恶心咋办
已经有13人回复
-
本人最近太闲了，谁有问题可以提，每天会统一回复
已经有12人回复
-
护理论文 晋升
已经有5人回复
-
求助大佬sci投稿哪个好中
已经有4人回复
-
评审感受-评审感受-评审感受
已经有14人回复
-
这年头没有找到涵评专家，还有中面上的可能吗
已经有6人回复
-
2026博士申请求助
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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