24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 639  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿复旦材料,英一专硕,总分357调剂 +4 1050389037 2026-04-02 5/250 2026-04-02 21:40 by dongzh2009
[考研] 070300化学求调剂 +8 起个名咋这么难 2026-03-27 9/450 2026-04-02 20:42 by dongzh2009
[考研] 求调剂 +7 Aniyaio 2026-04-02 7/350 2026-04-02 16:42 by zzsw+
[考研] 310求调剂 +17 争取九点睡 2026-03-30 17/850 2026-04-02 16:40 by guanxin1001
[考研] 324分 085600材料与化工 +20 呆鹅oor 2026-03-27 20/1000 2026-04-02 10:13 by oooqiao
[考研] 348环境工程调剂 +3 吴彦祖24k 2026-04-01 3/150 2026-04-02 09:14 by nanaliuyun
[考研] 301求调剂 +13 A_JiXing 2026-04-01 13/650 2026-04-02 09:01 by sanrepian
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
[考研] 调剂 +5 好好读书。 2026-03-28 7/350 2026-04-01 15:32 by 王亮_大连医科大
[考研] 262求调剂 +9 励志一定发文章 2026-03-31 10/500 2026-04-01 12:22 by sunshine0013
[考研] 326求调剂 +4 崽崽仔 2026-03-31 4/200 2026-04-01 09:58 by 我的船我的海
[考研] 0856材料化工调剂 总分330 +18 zhubinhao 2026-03-27 18/900 2026-04-01 09:37 by oooqiao
[考研] 254材料与化工求调剂 +3 翰冬林楠 2026-03-30 4/200 2026-03-31 17:53 by yishunmin
[考研] 08工科,295,接受跨专业调剂 +6 lmnlzy 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[考研] 22408 359分调剂 +4 Qshers 2026-03-27 8/400 2026-03-31 08:53 by Qshers
[考研] 一志愿中海洋320化学工程与技术学硕求调剂 +8 披星河 2026-03-30 8/400 2026-03-31 08:53 by lbsjt
[考研] 08工科求调剂286 +5 tgs_001 2026-03-28 5/250 2026-03-31 08:18 by 一只好果子?
[考研] 哈尔滨工业大学材料与化工专硕378求调剂 +3 塔比乌斯 2026-03-30 3/150 2026-03-30 22:55 by 无际的草原
[考研] 085602 化学工程专硕 340分求调剂 +4 qianbai11 2026-03-29 4/200 2026-03-30 11:34 by 唐沐儿
信息提示
请填处理意见