24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 642  |  回复: 2

japus1988

新虫 (初入文坛)

[求助] 求助属性,这样写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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

b07988

金虫 (正式写手)

星期六可以休息

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
japus1988: 金币+2, ★★★很有帮助, 谢谢,我也发现了这个问题了 2012-12-11 20:05:39
xiegangmai: 金币+1, 谢谢参与 2012-12-11 20:33:28
The value of domain_id is always 1 for the mixture domain.所以#define domain_ID 应该为2或者3
后面的就不清楚正确不正确了。
先做最重要的事情,再做次重要的事情
2楼2012-12-11 00:13:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

highway1

新虫 (初入文坛)

【答案】应助回帖

楼主问题解决了吗?
3楼2018-05-02 17:04:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 japus1988 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕调剂 +18 椰椰。 2026-03-29 18/900 2026-04-03 16:45 by 玲玲0606
[考研] 调剂 +6 JLLLLLLLLLL 2026-04-03 6/300 2026-04-03 15:47 by 晴空210210
[考研] 283求调剂 +3 jiouuu 2026-04-03 4/200 2026-04-03 13:28 by jiouuu
[考研] 调剂 +7 祉岷. 2026-04-02 7/350 2026-04-03 09:11 by 花呗还欠600
[考研] 一志愿山东大学,085600,344 +7 魏子per 2026-04-02 8/400 2026-04-02 21:12 by 百灵童888
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 找调剂 +6 楚乔乔 2026-04-01 6/300 2026-04-02 12:06 by xiaoranmu
[考研] 0856初试324分求调剂 +6 想上学求调 2026-04-01 6/300 2026-04-02 11:42 by 星空星月
[考研] 372求调剂 +3 jj涌77 2026-04-02 3/150 2026-04-02 09:57 by olim
[考研] 一志愿346上海大学生物学 +3 上海大学346调剂 2026-04-01 3/150 2026-04-02 08:36 by w虫虫123
[考研] 379求调剂 +3 ?苦瓜不苦 2026-04-01 3/150 2026-04-01 20:09 by 暮云清寒
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 11:35 by wangjy2002
[考研] 求调剂:一志愿:南京大学 专业:0705 总分320 ,本科985,四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
[考研] 本2一志愿C9-333分,材料科学与工程,求调剂 +9 升升不降 2026-03-31 9/450 2026-03-31 18:01 by 无际的草原
[考研] 080500-315分复试调剂 +9 上岸3821 2026-03-31 9/450 2026-03-31 17:29 by 唐沐儿
[考研] 340求调剂 +4 希望如此i 2026-03-31 4/200 2026-03-31 16:40 by 690616278
[考研] 材料与化工272求调剂 +25 阿斯蒂芬2004 2026-03-28 25/1250 2026-03-31 16:27 by hypershenger
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
[考研] 291求调剂 +5 Y-cap 2026-03-29 6/300 2026-03-29 13:18 by mumin1990
[考研] 求调剂 +7 争取九点睡 2026-03-28 8/400 2026-03-28 21:07 by 争取九点睡
信息提示
请填处理意见