24小时热门版块排行榜    

查看: 582  |  回复: 2

gjmfreeboy

新虫 (初入文坛)

[求助] UDF水蒸气体积分数传递问题 已有1人参与

目的在于在mixture模型下,二维模型,将某个面的水蒸气体积分数传递至两一个面,但不知道为什么编译没问题,一运行就出问题?求解答,谢谢
DEFINE_ADJUST(humidity_avg_temp, mixture_domain)
{      
        real sum_T_B=0.0,sum_B=0.0;
        cell_t c;
        Domain *subdomain;
        int phase_domain_index;
        Thread*c_thread;
        c_thread = Lookup_Thread(mixture_domain, 33);
      
        sub_domain_loop(subdomain,mixture_domain, phase_domain_index)

         {
           if(DOMAIN_ID(subdomain==1))
          thread_loop_c(c_thread,subdomain)
         {
         begin_c_loop_all(c,c_thread)
        {
            
                sum_B+=1;                  
                sum_T_B+=C_VOF(c, c_thread);
              
        }
        end_c_loop_all(c, c_thread)
           
}
}
avg_temp7=sum_T_B/sum_B;  
}   

DEFINE_PROFILE(humidity_Y,t,position)
{
        face_t f;
        begin_f_loop(f,t)
        {
                F_PROFILE(f,t,position) = avg_temp7;
  }
        end_f_loop(f,t)
}
同时,我尝试获取质量分数,同样出现错误
DEFINE_ADJUST(humidity_avg_temp, domain)
{      
        real sum_T_B=0.0,sum_B=0.0;
        cell_t c;
        int i=2;
        Thread*c_thread;
        c_thread = Lookup_Thread(domain, 33);
        begin_c_loop(c, c_thread)
        {
            
                sum_B+=1;                  
                sum_T_B+=C_YI(c, c_thread,i);
              
        }
        end_c_loop(c, c_thread)
         avg_temp7=sum_T_B/sum_B;   
}
      

DEFINE_PROFILE(humidity_Y,t,position)
{
        face_t f;
        begin_f_loop(f,t)
        {
                F_PROFILE(f,t,position) = avg_temp7;
  }
        end_f_loop(f,t)
}
回复此楼

» 猜你喜欢

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

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

小七工作室1

新虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
检查下程序是否有问题
2楼2015-05-07 12:02:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gjmfreeboy

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 小七工作室1 at 2015-05-07 12:02:55
检查下程序是否有问题

检查了好几遍,都没发现问题
3楼2015-05-07 12:41:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gjmfreeboy 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 20:50 by MHf40yA0Uhyb
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 18:37 by rRvxEj3zorCy
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-16 18:00 by i3KBRD19ERkG
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:48 by i3KBRD19ERkG
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-16 15:24 by i3KBRD19ERkG
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 14:39 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-16 07:12 by LH5NkK5Ud4bh
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-16 06:51 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:27 by LH5NkK5Ud4bh
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 05:03 by LH5NkK5Ud4bh
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 04:51 by LH5NkK5Ud4bh
[考博] 申博发邮件 +4 Lmengk 2026-09-15 4/200 2026-09-16 00:05 by 木心想继续深造
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见