24小时热门版块排行榜    

查看: 2502  |  回复: 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
回复此楼

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

多相流UDF

» 猜你喜欢

已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 俊采懿范 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 什么时候开奖? +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见