应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 参照论文写了湿空气结霜的udf，成功编译后初始化出错，恳请各位大神指点！！！
81/1返回列表
查看: 2175  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

俊采懿范

新虫 (初入文坛)

[求助] 参照论文写了湿空气结霜的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
回复此楼

» 收录本帖的淘帖专辑推荐

多相流UDF

» 猜你喜欢
1楼 2017-07-28 10:00:26
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 俊采懿范 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭