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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 【求助】FLUENT UDF 气液两相流 汽蚀 汽化潜热(问题详见内容)
53/1返回列表
查看: 5065  |  回复: 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;
}
回复此楼

» 猜你喜欢
1楼2012-07-30 17:51:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梅傥ojps

银虫 (正式写手)
我想问下各位这个蒸发冷凝公式是哪位提出来的，在哪偏文章里有的，急用
赞一下 回复此楼
22楼2018-03-21 16:33:35
 已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[高分子] 寻找聚酯反应釜 +3 茕茕恭煮 2024-06-15 5/250 2024-06-18 11:44 by 茕茕恭煮
[基金申请] 有人中过人文社科类的博后特助吗？ +3 outsider1986 2024-06-16 5/250 2024-06-18 11:10 by 袁天未然
[找工作] 药学硕士找不到工作，打算去做科研助理了 +6 pom戴墨镜 2024-06-14 17/850 2024-06-18 09:54 by lialee22
[精细化工] 试剂价格 +3 旭必上岸 2024-06-17 3/150 2024-06-18 09:13 by kzx_33
[考博] 34岁读博士晚吗 +32 emitdne 2024-06-13 32/1600 2024-06-18 08:44 by yuan0806
[基金申请] F口信息学部拿面上，大概需要什么样的成果 +7 _奋黎_ 2024-06-16 12/600 2024-06-18 07:41 by ld_cauchy
[论文投稿] ACS Nano投稿后分配到副编辑手里12天了，能确定送审了吗？ +6 潇洒怡惜 2024-06-12 11/550 2024-06-18 01:48 by 玩耍绝缘体
[论文投稿] 审稿问题：为什么荧光激发波长和紫外吸收波长差的大？ 10+5 sdawege 2024-06-14 10/500 2024-06-17 18:54 by HH-探针
[教师之家] 请问事业编制和年薪制冲突吗？ +10 ZHONGWU_U 2024-06-14 10/500 2024-06-17 14:44 by 周周520
[找工作] 江西双非一本和四川双一流高校如何选择？ 5+8 寒山敲钟 2024-06-12 25/1250 2024-06-16 22:05 by zhaojiang427
[论文投稿] 投稿被一个审稿人恶意评审了怎么样？ +5 1chen 2024-06-14 7/350 2024-06-15 23:15 by xy66xy
[论文投稿] 投稿时忘记修改一作 +7 gll123456 2024-06-13 11/550 2024-06-15 11:49 by gll123456
[基金申请] 博后基金，以往的结果点不开，怎么回事呢？最后一次机会了，两次都没中前面。 +7 kyukitu 2024-06-14 13/650 2024-06-15 06:46 by 我是王小帅
[论文投稿] 投了一篇4区的SCI，审稿人一个拒稿，一个小修，编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[考研] 物理化学一对一辅导 +3 林大diao 2024-06-12 5/250 2024-06-14 20:57 by 林大diao
[基金申请] 博士后基金需要结题吗？ +8 zhouchuck 2024-06-13 8/400 2024-06-14 17:27 by liuyupu132
[基金申请] 工材E口JQ有消息了吗 +4 babyduck 2024-06-11 4/200 2024-06-14 17:23 by firepick
[论文投稿] 投稿后发现其他作者的邮箱填错了该怎么办呀 10+4 在飞的猪 2024-06-13 6/300 2024-06-14 04:45 by 小虫子咔咔
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
信息提示
关闭
请填处理意见
关闭