| 查看: 1012 | 回复: 4 | ||
| 【悬赏金币】回答本帖问题，作者SH1201将赠送您 10 个金币 | ||
[求助]
fluent加载udf后，初始化就闪退
|
||
|
写了一个简单的蒸发冷凝udf，不加载udf时可以正常运行，加入udf编辑的源项后，点击初始化后就闪退了，udf可以通过编译不知道为什么，求大神帮助 代码： #include "udf.h" #include "math.h" #include "sg_mphase.h" #define operating_p 101325 #define la_heat 2.257e6 DEFINE_SOURCE(liq,c,t,dS,eqn) { real sat_tem , tem_pressure,pressure,p_wv,m_dot_l; Thread *t_air = THREAD_SUB_THREAD(t,0); Thread *t_wl = THREAD_SUB_THREAD(t,1); Thread *t_wv = THREAD_SUB_THREAD(t_air,0); pressure = C_P(c,t_air)+operating_p; //pressure of wet_air p_wv = pressure*C_VOF(c,t_wv); //pressure of water vapor sat_tem = 1657.459/(7.074056-log10(p_wv+101325))+46.13; //define sat_tem if (C_T(c,t_wl)>=sat_tem) { m_dot_l = -0.1*C_VOF(c,t_wl)*C_R(c,t_wl)*fabs(C_T(c,t_wl)-sat_tem)/sat_tem; //liquid vaporation dS[eqn] = -0.1*C_R(c,t_wl)*fabs(C_T(c,t_wl)-sat_tem)/sat_tem; } else { m_dot_l = 0.1*C_VOF(c,t_wl)*C_R(c,t_wl)*fabs(C_T(c,t_wl)-sat_tem)/sat_tem; dS[eqn] = 0; } return m_dot_l; } DEFINE_SOURCE(vap,c,t,dS,eqn) { real sat_tem , tem_pressure,pressure,p_wv,m_dot_v; Thread *t_air = THREAD_SUB_THREAD(t,0); Thread *t_wl = THREAD_SUB_THREAD(t,1); Thread *t_wv = THREAD_SUB_THREAD(t_air,0); pressure = C_P(c,t_air)+operating_p; //pressure of wet_air p_wv = pressure*C_VOF(c,t_wv); //pressure of water vapor sat_tem = 1657.459/(7.074056-log10(p_wv+101325))+46.13; //define sat_tem if (C_T(c,t_wl)>=sat_tem) { m_dot_v = 0.1*C_VOF(c,t_wl)*C_R(c,t_wl)*fabs(C_T(c,t_wl)-sat_tem)/sat_tem; //liquid vaporation dS[eqn] = 0; } else { m_dot_v = -0.1*C_VOF(c,t_wl)*C_R(c,t_wl)*fabs(C_T(c,t_wl)-sat_tem)/sat_tem; dS[eqn] = 0.1*C_R(c,t_wl)*fabs(C_T(c,t_wl)-sat_tem)/sat_tem; } return m_dot_v; } DEFINE_SOURCE(energy,c,t,dS,eqn) { real sat_tem , tem_pressure,pressure,p_wv,m_dot; Thread *t_air = THREAD_SUB_THREAD(t,0); Thread *t_wl = THREAD_SUB_THREAD(t,1); Thread *t_wv = THREAD_SUB_THREAD(t_air,0); pressure = C_P(c,t_air)+operating_p; //pressure of wet_air p_wv = pressure*C_VOF(c,t_wv); //pressure of water vapor sat_tem = 1657.459/(7.074056-log10(p_wv+101325))+46.13; //define sat_tem if (C_T(c,t_wl)>=sat_tem) { m_dot = -0.1*C_VOF(c,t_wl)*C_R(c,t_wl)*fabs(C_T(c,t_wl)-sat_tem)/sat_tem; dS[eqn] = -0.1*C_VOF(c,t_wl)*C_R(c,t_wl)/sat_tem; } else { m_dot = 0.1*C_VOF(c,t_wl)*C_R(c,t_wl)*fabs(C_T(c,t_wl)-sat_tem)/sat_tem; dS[eqn] = 0.1*C_VOF(c,t_wl)*C_R(c,t_wl)/sat_tem; } return la_heat*m_dot; |
回复此楼» 猜你喜欢
博士找工作真难
已经有16人回复
-
竟然失眠了
已经有4人回复
-
迷茫的毕业季
已经有3人回复
-
这样的课题组值得去吗
已经有5人回复
-
一般很空的审稿意见大家一般是怎么回复的
已经有6人回复
-
蹲一个6月7号江苏省青基公示
已经有25人回复
-
间接经费用于绩效支持比例人员等问题
已经有4人回复
-
席夫碱反应ph敏感 控制ph在多少呢？
已经有5人回复
-
难上的课安排给我，杂事安排给我。评优时我就没份儿了，即便我做得不比绝大多数人差
已经有11人回复
-
看到大家一直在说代表作NC，不明觉厉
已经有22人回复
2楼2023-02-01 21:20:56
3楼2023-10-31 18:51:23
4楼2023-10-31 18:53:31
5楼2023-11-07 21:44:52
30