版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (3807)
>导师招生 (642)
>虫友互识 (533)
>文献求助 (353)
>考研 (286)
>招聘信息布告栏 (185)
>硕博家园 (179)
>休闲灌水 (159)
>论文投稿 (102)
>考博 (89)
>博后之家 (79)
>教师之家 (43)
>SciFinder/Reaxys (38)
>公派出国 (37)
>找工作 (31)
>基金申请 (29)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

enigmafwj

捐助贵宾 (小有名气)

应助: 1 (幼儿园)
金币: 0.4
散金: 21
帖子: 106
在线: 35小时
虫号: 1917147
注册: 2012-07-30
专业: 传热传质学

[求助] 【求助】FLUENT UDF 气液两相流汽蚀汽化潜热(问题详见内容)

为液相、气相质量方程和能量方程添加源项，UDF程序如下，汽化、液化由压力控制（压力条件判断），各物性参数用温度拟合。问题是：编译时，提示有语法错误(parse)，并且fv为未声明的变量。请问如何解决？另：程序中有其他不妥之处，烦请指教，谢谢！
#include "udf.h"
#include "sg_mphase.h"

DEFINE_SOURCE(liq_m,c,l_t,dS,eqn) /*液相质量输运源项*/
{
Thread *m_t,*v_t;                /*计算区域指针*/
  m_t=THREAD_SUPER_THREAD(l_t);    /*混合区液相指针*/
  v_t=THREAD_SUB_THREAD(m_t,1);    /*单相区气相指针*/
real m_dot_l;                      /*液相质量输运量*/
real stc;                         /*液氧表面张力系数*/
real fv;                         /*气相质量分数*/
  fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));
real pres_vap,pres_sat,pres_turb; /*液氧饱和蒸汽压(考虑湍流影响)*/
  if(C_T(c,m_t)>154.58)
{
   stc=0;
   pres_sat=5.041e6;
}
  else if(C_T(c,m_t)>54.36)
{
   stc=0.038066*pow(1-C_T(c,l_t)/154.58,1.2136);
   pres_sat=133.322*pow(10,20.6695-526.97/C_T(c,l_t)-6.7062*log10(C_T(c,l_t))+0.012926*C_T(c,l_t)-C_T(c,l_t)*C_T(c,l_t)*(9.8832e-13));
}
  else
{
   stc=0.0225;
   pres_sat=146.041;
}
  pres_turb=0.39*C_R(c,l_t)*C_K(c,m_t);
  pres_vap=pres_sat+0.5*pres_turb;
if(C_P(c,m_t)<=pres_vap)          /*液相汽化*/
{
m_dot_l=-0.02*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,v_t)*sqrt((2*(pres_vap-C_P(c,m_t)))/(3*C_R(c,l_t)))*(1-fv);
dS[eqn]=-0.02*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,v_t)*sqrt((2*(pres_vap-C_P(c,m_t)))/(3*C_R(c,l_t)));
}
else                               /*气相液化*/
{
m_dot_l=0.01*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,l_t)*sqrt((2*(C_P(c,m_t)-pres_vap))/(3*C_R(c,l_t)))*fv;
dS[eqn]=0;
}
return m_dot_l;
}

DEFINE_SOURCE(vap_m,c,v_t,dS,eqn) /*气相质量输运源项*/
{
Thread *m_t,*l_t;                /*计算区域指针*/
  m_t=THREAD_SUPER_THREAD(v_t);    /*混合区气相指针*/
  l_t=THREAD_SUB_THREAD(m_t,0);    /*单相区液相指针*/
real m_dot_v;                      /*气相质量输运量*/
real stc;                         /*液氧表面张力系数*/
real fv;                         /*气相质量分数*/
  fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));
real pres_vap,pres_sat,pres_turb; /*液氧饱和蒸汽压(考虑湍流影响)*/
  if(C_T(c,m_t)>154.58)
{
   stc=0;
   pres_sat=5.041e6;
}
  else if(C_T(c,m_t)>54.36)
{
   stc=0.038066*pow(1-C_T(c,l_t)/154.58,1.2136);
   pres_sat=133.322*pow(10,20.6695-526.97/C_T(c,l_t)-6.7062*log10(C_T(c,l_t))+0.012926*C_T(c,l_t)-C_T(c,l_t)*C_T(c,l_t)*(9.8832e-13));
}
  else
{
   stc=0.0225;
   pres_sat=146.041;
}
  pres_turb=0.39*C_R(c,l_t)*C_K(c,m_t);
  pres_vap=pres_sat+0.5*pres_turb;
if(C_P(c,m_t)<=pres_vap)          /*液相汽化*/
{
m_dot_v=0.02*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,v_t)*sqrt((2*(pres_vap-C_P(c,m_t)))/(3*C_R(c,l_t)))*(1-fv);
dS[eqn]=0;
}
else                               /*气相液化*/
{
m_dot_l=-0.01*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,l_t)*sqrt((2*(C_P(c,m_t)-pres_vap))/(3*C_R(c,l_t)))*fv;
dS[eqn]=-0.01*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,l_t)*sqrt((2*(C_P(c,m_t)-pres_vap))/(3*C_R(c,l_t)));
}
return m_dot_v;
}

DEFINE_SOURCE(m_enrg,c,m_t,dS,eqn) /*混合相能量方程输运源项*/
{
Thread *l_t,*v_t;                /*计算区域指针*/
  l_t=THREAD_SUB_THREAD(m_t,0);    /*混合区液相指针*/
  v_t=THREAD_SUB_THREAD(m_t,1);    /*混合区气相指针*/
real m_dot;                      /*相变质量*/
real hc;                         /*汽化潜热*/
real source;                      /*能量方程源项*/
real stc;                         /*液氧表面张力系数*/
real fv;                         /*气相质量分数*/
  fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));
real pres_vap,pres_sat,pres_turb; /*液氧饱和蒸汽压(考虑湍流影响)*/
  if(C_T(c,m_t)>154.58)
{
   stc=0;
   pres_sat=5.041e6;
   hc=0;
}
  else if(C_T(c,m_t)>54.36)
{
   stc=0.038066*pow(1-C_T(c,l_t)/154.58,1.2136);
   pres_sat=133.322*pow(10,20.6695-526.97/C_T(c,l_t)-6.7062*log10(C_T(c,l_t))+0.012926*C_T(c,l_t)-C_T(c,l_t)*C_T(c,l_t)*(9.8832e-13));
   hc=8040*pow(1-C_T(c,l_t)/154.58,0.201);
}
  else
{
   stc=0.0225;
   pres_sat=146.041;
   hc=2523;
}
  pres_turb=0.39*C_R(c,l_t)*C_K(c,m_t);
  pres_vap=pres_sat+0.5*pres_turb;
if(C_P(c,m_t)<=pres_vap)          /*液相汽化，温度升高，源项为正*/
{
m_dot=0.02*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,v_t)*sqrt((2*(pres_vap-C_P(c,m_t)))/(3*C_R(c,l_t)))*(1-fv);
source=m_dot*hc;
dS[eqn]=0;
}
else                               /*气相液化，温度降低，源项为负*/
{
m_dot=-0.01*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,l_t)*sqrt((2*(C_P(c,m_t)-pres_vap))/(3*C_R(c,l_t)))*fv;
source=-m_dot*hc;
dS[eqn]=0;
}
return source;
}

回复此楼

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

多相流UDF

传热传质

» 猜你喜欢

0854求调剂已经有6人回复
求调剂推荐已经有6人回复
材料299专硕求调剂已经有17人回复
一志愿沪9，生物学326求调剂已经有5人回复
294求调剂已经有4人回复
297，工科调剂?河南农业大学本科已经有11人回复
22专硕求调剂已经有11人回复
284求调剂已经有20人回复
291求调剂已经有7人回复
296求调剂已经有5人回复

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

求助！！！我想要在fluent中针对一个面进行积分，那么UDF应该用什么宏函数？怎么用？已经有5人回复
FLUENT求帮助啊啊啊啊啊已经有8人回复
好心人帮帮忙吧，为什么这个UDF导入fluent以后，在边界条件处没有显示~~~ 已经有17人回复
fluent中用udf编写三维速度边界条件出现错误，大家来指导一下啊已经有9人回复
FLUENT源项UDF非稳态迭代的一个疑问请教已经有8人回复
有关fluent的udf中的hook 已经有6人回复
【求助】气液分离两相流动模拟fluent gambit出入口的设置已经有11人回复
【求助】单机多核fluent并行udf 已经有13人回复
【求助】fluent计算中，udf编译通过，初始化出错已经有11人回复
【求助】如何改变Fluent中的曳力模型（UDF）已经有12人回复

1楼 2012-07-30 17:51:36

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

enigmafwj

捐助贵宾 (小有名气)

应助: 1 (幼儿园)
金币: 0.4
散金: 21
帖子: 106
在线: 35小时
虫号: 1917147
注册: 2012-07-30
专业: 传热传质学

引用回帖:

9楼: Originally posted by kylafree at 2012-08-02 22:00:06
那恭喜啊，有时间给我讲解下我的疑问啊？

...

我刚学这些个东西，只是按自己的理解用它，你觉得这两函数是否是这样用的？

赞一下

回复此楼

10楼2012-08-03 16:01:53

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

查看全部 22 个回答

kylafree

至尊木虫 (知名作家)

应助: 166 (高中生)
金币: 17562
散金: 883
红花: 18
帖子: 6227
在线: 1161.8小时
虫号: 1061079
注册: 2010-07-20
性别: GG
专业: 流体力学

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢参与！ 2012-07-30 21:50:07
enigmafwj: 金币+2 2012-12-18 20:11:11

real fv;                         /*气相质量分数*/
fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));
real pres_vap,pres_sat,pres_turb; /*液氧饱和蒸汽压(考虑湍流影响)*/

改为：
real pres_vap,pres_sat,pres_turb; /*液氧饱和蒸汽压(考虑湍流影响)*/
real fv;                         /*气相质量分数*/
fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));
定义变量要在赋值语句开始之前全部定义完毕。

另外向兄台请教一下，这两句的thread具体指什么thread，感激不尽啊。
m_t=THREAD_SUPER_THREAD(v_t);    /*混合区气相指针*/
l_t=THREAD_SUB_THREAD(m_t,0);    /*单相区液相指针*/
为啥第一句是“混合区气相指针”，而第二句跑到“单相区液相指针”了呢

赞一下(1人)

回复此楼

2楼2012-07-30 21:19:25

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

enigmafwj

捐助贵宾 (小有名气)

应助: 1 (幼儿园)
金币: 0.4
散金: 21
帖子: 106
在线: 35小时
虫号: 1917147
注册: 2012-07-30
专业: 传热传质学

引用回帖:

2楼: Originally posted by kylafree at 2012-07-30 21:19:25
real fv; /*气相质量分数*/
fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));
real pres_vap,pres_sat,pres_turb; /*液氧饱和蒸汽压(考虑湍流影 ...

因为你所问的这两句，是在气相质量方程源项里的，DEFINE_SOURCE(vap_m,c,v_t,dS,eqn)里的thread：v_t在该宏应用于气相质量方程时指向了气相区。因为之后要用到混合区气相的参数，所以用THREAD_SUPER_THREAD(v_t)把指针恢复到混合区气相并传给m_t；
THREAD_SUB_THREAD(m_t,0)的意思是在m_t得到混合区气相指针后，m_t作为最高区域(混合区)的指针，可以用这个函数把单相区液相指针传给l_t，才能用l_t调用到液相区的参数，0指的是主相(液相)，1指的是次相(气相)。

赞一下

回复此楼

3楼2012-07-31 15:49:14

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

kylafree

至尊木虫 (知名作家)

应助: 166 (高中生)
金币: 17562
散金: 883
红花: 18
帖子: 6227
在线: 1161.8小时
虫号: 1061079
注册: 2010-07-20
性别: GG
专业: 流体力学

★ ★ ★
xiegangmai: 金币+3, 鼓励讨论交流！ 2012-08-02 23:05:49

引用回帖:

3楼: Originally posted by enigmafwj at 2012-07-31 15:49:14
因为你所问的这两句，是在气相质量方程源项里的，DEFINE_SOURCE(vap_m,c,v_t,dS,eqn)里的thread：v_t在该宏应用于气相质量方程时指向了气相区。因为之后要用到混合区气相的参数，所以用THREAD_SUPER_THREAD(v_t)把 ...

感谢兄台，受教了。我本来没做多相流，但对这个有点兴趣。
再问几个小白的问题可以不。比如我按下面这样理解对头不：
气相区----就是指这个区域的CELL只有气相；
液相区----就是指这个区域的CELL只有液相；
混合区----就是指这个区域的CELL既有气相也有液相；
混全区thread----是说混合区既然有混合thread，又有气相thread，还有液相thead。气相thread液相thread混合thread怎么联系起来啊，糊涂了

赞一下

回复此楼

4楼2012-07-31 16:22:41

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

查看全部 22 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 22专硕求调剂 +9	haoyun上岸 2026-04-11	11/550	2026-04-15 14:47 by elainzaizai
[考研] 310求调剂 +16	666真好 2026-04-11	18/900	2026-04-15 13:28 by 黑科技矿业
[考研] 调剂求收留 +34	果然有我 2026-04-10	35/1750	2026-04-15 13:05 by 西北望—风沙
[考研] 本科211，报考085601-310分 +16	ararak 2026-04-13	16/800	2026-04-14 14:55 by Delta2012
[考研] 求调剂，985材料与化工348分 +9	涵竹刘 2026-04-11	14/700	2026-04-13 22:26 by 涵竹刘
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11	不会吃肉 2026-04-12	11/550	2026-04-13 21:59 by bljnqdcc
[考研] 344 材料专业求调剂211 无地域要求 +8	hualkop 2026-04-11	8/400	2026-04-12 22:24 by fqwang
[考研] 材料与化工300求调剂 +39	肖开文 2026-04-09	43/2150	2026-04-12 01:30 by 秋豆菜芽
[考研] 22408调剂315分 +3	zhuangyan123 2026-04-09	3/150	2026-04-12 00:25 by 蓝云思雨
[考研] 0860004 求调剂 309分 +9	Yin DY 2026-04-08	9/450	2026-04-11 22:55 by dongdian1
[考研] 调剂 +10	只叙离别辞 2026-04-09	12/600	2026-04-11 20:57 by 逆水乘风
[考研] 270求调剂 +14	杨乐369 2026-04-11	14/700	2026-04-11 20:16 by 蓝云思雨
[考研] 359求调剂 +5	胃痉挛累了 2026-04-11	5/250	2026-04-11 19:55 by lbsjt
[考研] 283求调剂 086004考英二数二 +17	那个噜子 2026-04-10	18/900	2026-04-11 16:27 by 明月此时有
[考研] 农业管理302分求调剂 +3	xuening1 2026-04-10	3/150	2026-04-11 10:18 by zhq0425
[考研] 求调剂 +5	不会飞的鱼@ 2026-04-10	5/250	2026-04-10 19:07 by chemisry
[考研] 314求调剂 +18	xhhdjdjsjks 2026-04-09	19/950	2026-04-10 18:53 by HPUCZ
[考研] 初试261 +3	Asht少 2026-04-10	6/300	2026-04-10 16:38 by Asht少
[考研] 本科西工大 0856 324求调剂 +10	wysyjs25 2026-04-09	11/550	2026-04-10 08:37 by 5268321
[考研] 070300化学求调剂 +13	73372112 2026-04-08	13/650	2026-04-09 20:22 by maddjdld

24小时热门版块排行榜

enigmafwj

[求助] 【求助】FLUENT UDF 气液两相流 汽蚀 汽化潜热(问题详见内容)

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

» 猜你喜欢

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

enigmafwj

kylafree

【答案】应助回帖

enigmafwj

kylafree

[求助] 【求助】FLUENT UDF 气液两相流汽蚀汽化潜热(问题详见内容)