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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 【求助】FLUENT UDF 气液两相流 汽蚀 汽化潜热(问题详见内容)
51/1返回列表
查看: 5069  |  回复: 21
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

enigmafwj

捐助贵宾 (小有名气)

[求助] 【求助】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 传热传质

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-07-30 17:51:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enigmafwj

捐助贵宾 (小有名气)
引用回帖:
4楼: Originally posted by kylafree at 2012-07-31 16:22:41
感谢兄台，受教了。我本来没做多相流，但对这个有点兴趣。
再问几个小白的问题可以不。比如我按下面这样理解对头不：
气相区----就是指这个区域的CELL只有气相；
液相区----就是指这个区域的CELL只有液相；
混 ...

你看下UDF帮助的第3章第11节的区域结构和第6章THREAD函数。
赞一下 回复此楼
5楼2012-08-02 09:54:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 22 个回答

kylafree

至尊木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +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

捐助贵宾 (小有名气)
引用回帖:
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

至尊木虫 (知名作家)
★ ★ ★
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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 博后网上办公系统维护 +14 取名字烦人 2024-06-25 22/1100 2024-06-26 22:15 by tuanbao2023
[找工作] offer选择 +3 羡渔樵 2024-06-26 3/150 2024-06-26 21:57 by 陈大胖
[硕博家园] 夏至，要不要硕博联谊 +6 我是王小帅 2024-06-21 8/400 2024-06-26 21:05 by zxw3578
[有机交流] 大佬们，打的核磁氢谱与chemdraw预测的有偏差 +5 来了个 2024-06-26 5/250 2024-06-26 17:28 by Jasonlin7758
[论文投稿] 我在写关于多维空间的论文，希望能与大家交流，准备在science上发表 +8 lgf519 2024-06-24 10/500 2024-06-26 16:25 by 梦渺岚烟
[基金申请] 青年基金E02口青基去年几个函评专家？ +6 他山攻玉之石 2024-06-25 9/450 2024-06-26 15:09 by 他山攻玉之石
[基金申请] 省基金对国自然申请是有利还是有弊 +6 zycumt04 2024-06-26 7/350 2024-06-26 15:02 by 潇湘之迷
[硕博家园] 【45岁以上博士】柔性人才引进项目，有补贴 +9 Dreamsummit 2024-06-24 18/900 2024-06-26 12:28 by hujm159
[基金申请] 基金申请书名称有变化 +5 xuel2011 2024-06-25 7/350 2024-06-26 00:07 by 老虎当猫养
[基金申请] 博后面上到底什么时候出结果？？？ +7 爱学lsy 2024-06-24 7/350 2024-06-25 19:24 by 暴走的蒜泥
[基金申请] 2024安徽省哲社立项名单公示 +3 robin_work 2024-06-24 6/300 2024-06-25 14:37 by Pickfoot
[基金申请] 演员杨幂在核心期刊发表论文，中国广播电视学期刊 +9 babu2015 2024-06-22 12/600 2024-06-25 11:05 by 我是勇敢牛牛
[有机交流] 高温酯化反应喷料 20+4 hl24678 2024-06-21 7/350 2024-06-25 09:45 by ghj-hebust
[考博] 申请25博士，可以提前进组做科研助理 +6 逐梦途中w 2024-06-22 8/400 2024-06-24 21:49 by keyaner23
[基金申请] 说博后基金7月出的真打电话了吗？ +12 antonysole 2024-06-24 14/700 2024-06-24 13:39 by sizhouyi
[基金申请] 青年和面上，哪个上会难度更大 +12 今晚推荐22 2024-06-21 18/900 2024-06-24 11:08 by 半简体
[催化] 求助原位红外 +3 jfdhj 2024-06-21 7/350 2024-06-24 10:52 by qvhm2609
[基金申请] 国自然青年基金，1A4B能上会吗？青年和面上的上会标准是一样的吗？ +19 今晚推荐22 2024-06-20 32/1600 2024-06-23 23:17 by andywei1028
[有机交流] 生成亚胺的反应怎么能进行完全 +3 1369836 2024-06-23 3/150 2024-06-23 18:44 by hwqMSE
[基金申请] 听大佬说今年信息口本子数量大幅增加？ +8 wutzxt 2024-06-21 9/450 2024-06-21 19:58 by wutzxt
信息提示
关闭
请填处理意见
关闭