24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2309  |  回复: 7

俊采懿范

新虫 (初入文坛)

[求助] 参照论文写了湿空气结霜的udf,成功编译后初始化出错,恳请各位大神指点!!! 已有3人参与

#include "udf.h"
#include "sg_mphase.h"


DEFINE_MASS_TRANSFER(liq_gas_source,cell,thread,from_phase_index,from_species_index,to_phase_index,to_species_index)
{
        real hr;                   //饱和空气含湿量//
        real Wvs;                  //饱和水蒸气质量分数//
        real Wv;                  //湿空气水蒸气质量分数//
        real m_lg_m;                 //质量转移//
        real Pvs;                   //水蒸气饱和压力//
        real Pv;                //水蒸气饱和分压力//
        real Td;               
       
       
        Thread *tf= THREAD_SUPER_THREAD(thread,thread);
        Thread *mix = THREAD_SUB_THREAD(thread,from_phase_index);
        Thread *ice = THREAD_SUB_THREAD(thread,to_phase_index);
       
        if(C_T(cell,tf)<=273.15 && C_T(cell,tf)>173.15)
        {
                Pvs=exp(-5674.5359/C_T(cell,tf)
                        +6.3925247
                        -9.677843*pow(10,-3)*C_T(cell,tf)
                        +6.22157*pow(10,-7)*C_T(cell,tf)*C_T(cell,tf)
                        +2.0747825*pow(10,-9)*C_T(cell,tf)*C_T(cell,tf)*C_T(cell,tf)
                        +9.484024*pow(10,-13)*C_T(cell,tf)*C_T(cell,tf)*C_T(cell,tf)*C_T(cell,tf)+
                        4.1635019*log(C_T(cell,tf)));
        }
    else if
     
           (C_T(cell,tf)<=473.15 && C_T(cell,tf)>273.15)
    {
        Pvs=exp(-5.8002206*pow(10,3)/C_T(cell,tf)
                +1.3914993-4.8640239*pow(10,-2)*C_T(cell,tf)
                -4.1764768*pow(10,-5)*C_T(cell,tf)*C_T(cell,tf)
                -1.4452093*pow(10,-8)*C_T(cell,tf)*C_T(cell,tf)*C_T(cell,tf)+6.5459673*log(C_T(cell,tf)));
        }
                                                                                /*定义计算单元温度对应的水蒸气饱和压力*/
        hr=0.622*Pvs/(101325-Pvs);
        Wvs=hr/(1+hr);
        Wv=C_YI(cell,mix,0);
        Pv=(ABS_P(C_P(cell,mix),op_pres)*hr)/(hr+0.622);
       
        if
           (C_T(cell,tf)>=213.16&&C_T(cell,tf)<=273.15)
          {
            Td=(-60.45+7.0332*log(Pv)+0.37*log(Pv)*log(Pv))+273.15;
          }
        else if
            (C_T(cell,tf)>273.15&&C_T(cell,tf)<343.15)
            {
              Td=(-35.957-1.8726*log(Pv)+1.1689*log(Pv)*log(Pv))+273.15;
            }
     
        if
                (C_T(cell,tf)<=(Td<273.15?Td:273.15))
        {
                m_lg_m=C_VOF(cell,mix)*C_R(cell,mix)*Wv*(Wv-Wvs)/Wvs;   /*水蒸气向冰相质量转移的速率*/

        }
        else if (C_T(cell,tf)>(Td<273.15?Td:273.15))
        {
                m_lg_m=0;
        }
                return (m_lg_m);

}






出错方式:
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()


模型是使用欧拉+组分运输模型进行,同时自定义添加质量运输方程,恳请各位大神帮忙
@wuming524
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gyctju

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
其他的先不说,这一句  Thread *tf= THREAD_SUPER_THREAD(thread,thread); 是有逻辑错误的
2楼2017-07-28 13:20:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

俊采懿范

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by gyctju at 2017-07-28 13:20:51
其他的先不说,这一句  Thread *tf= THREAD_SUPER_THREAD(thread,thread); 是有逻辑错误的

谢谢,不过我想获取混合的温度,不是单相的,不知道在质量转移方程里面能定义吗?

发自小木虫IOS客户端
3楼2017-07-29 11:28:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果子酒

新虫 (著名写手)

楼主,我也在做结霜的课题,也是UDF来算,请问你进行的怎么样啦

发自小木虫IOS客户端
4楼2017-09-15 11:18:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果子酒

新虫 (著名写手)

引用回帖:
3楼: Originally posted by 俊采懿范 at 2017-07-29 11:28:43
谢谢,不过我想获取混合的温度,不是单相的,不知道在质量转移方程里面能定义吗?
...

楼主有没有试过用这个mass transfer算一下蒸发的算例?我用这个算不通相变啊

发自小木虫IOS客户端
5楼2017-09-15 11:45:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

metalfoamPCM

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by 果子酒 at 2017-09-15 11:45:36
楼主有没有试过用这个mass transfer算一下蒸发的算例?我用这个算不通相变啊
...

这个问题可有进展,看了一些论文成功算出来了,我也调不出,有空可以商量
6楼2017-11-16 14:57:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

muyuyou111

新虫 (初入文坛)

【答案】应助回帖

请问楼主现在有结霜的udf吗,我最近也在做结霜的模拟
7楼2020-05-08 16:56:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

koobm

铁虫 (小有名气)

【答案】应助回帖

引用回帖:
7楼: Originally posted by muyuyou111 at 2020-05-08 16:56:41
请问楼主现在有结霜的udf吗,我最近也在做结霜的模拟

可以交流下,我也在做
8楼2023-05-19 07:58:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 俊采懿范 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 论文撤稿了 +5 bjvtcliu 2026-05-24 8/400 2026-05-24 23:24 by zju2000
[考博] 博士申请 +6 星…… 2026-05-18 7/350 2026-05-24 22:45 by 预约这个秋天
[基金申请] 青B发送上会通知了吗 +5 chemBioBro 2026-05-22 8/400 2026-05-24 22:10 by Max0601
[基金申请] 评审有感 +16 popular289 2026-05-18 27/1350 2026-05-24 17:34 by hhs666
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 hvkbtfonbv 2026-05-23 4/200 2026-05-24 17:21 by 75ui6h7z2t
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:10 by 75ui6h7z2t
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 5/250 2026-05-24 16:21 by hhx1yx9evi
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 4/200 2026-05-24 16:16 by hhx1yx9evi
[基金申请] 河北省自然科学基金 +6 Peterchao 2026-05-18 9/450 2026-05-24 16:02 by 130067131
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路! +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[论文投稿] 投稿求助,期刊 +4 希冀,有书读 2026-05-20 8/400 2026-05-22 10:16 by 希冀,有书读
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
信息提示
请填处理意见