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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 网络生活区 » 有奖问答 » 学术帮助 » 想模拟单个气泡的溃灭，想用fluent编程把压力导入进去，用的vof模型。
北京石油化工学院2026年研究生招生接收调剂公告
61/1返回列表
查看: 2255  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

renzhihao

版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！

[求助] 想模拟单个气泡的溃灭，想用fluent编程把压力导入进去，用的vof模型。

我在做空化，想模拟单个气泡的溃灭，但是现在要求内部流体的压力是个随着位置变化的函数，想用fluent编程把压力导入进去，用的vof模型。

下面的为我从别的帖子里面看到的。想问下这些东西需要在fluent里面设置还是编程？
为液相、气相质量方程和能量方程添加源项，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;
}
这些东西直接在fluent里面设置行吗，还是必须要编程？
回复此楼

» 猜你喜欢
1楼 2016-04-21 17:27:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiachaolis

超级版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！

【答案】应助回帖

楼主，请把你相变的模型公式贴出来，这样就好查错误了
赞一下(1人) 回复此楼
付出就会有收获
2楼2016-04-24 20:01:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renzhihao

专家顾问

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
2楼: Originally posted by jiachaolis at 2016-04-24 20:01:44
楼主，请把你相变的模型公式贴出来，这样就好查错误了

怎么查错误
回复此楼
3楼2016-04-25 16:41:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renzhihao

版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
2楼: Originally posted by jiachaolis at 2016-04-24 20:01:44
楼主，请把你相变的模型公式贴出来，这样就好查错误了

怎么查错误
回复此楼
4楼2016-04-26 16:53:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wolfghost

新虫
楼主，长期无人应助金币可直接申请返还(点击右下角“通知版主处理”)，重新开贴以提高求助成功率！为营造互助、和谐论坛共同努力，谢谢~
赞一下 回复此楼
大刀向鬼子的头上砍去！ChopJapstoDeath!
5楼2016-05-25 13:50:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

聆听风声love

管理员

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！

【答案】应助回帖

质量源项的方程出自哪里?
赞一下 回复此楼
努力！奋斗！
6楼2016-07-03 11:41:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 renzhihao 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂申请 +7 张张张张zy 2026-03-31 7/350 2026-03-31 22:25 by piaoyi415
[电化学] 070300化学调剂 +7 山顶见α 2026-03-25 7/350 2026-03-31 22:21 by piaoyi415
[考研] 085602 307分 求调剂 +10 不知道叫什么！ 2026-03-26 10/500 2026-03-31 19:53 by Dyhoer
[考研] 070300化学354求调剂 +15 101次希望 2026-03-28 15/750 2026-03-31 17:58 by jp9609
[考研] 生医工0831调剂求推荐 +4 小熊睿睿_s 2026-03-27 6/300 2026-03-31 17:12 by 记事本2026
[考研] 物理学调剂 +4 小羊36 2026-03-30 4/200 2026-03-31 16:16 by lishahe
[基金申请] 面上5B能上会吗？ +8 redcom 2026-03-29 8/400 2026-03-31 15:53 by niuailing
[考研] 083000学硕274求调剂 +12 Li李鱼 2026-03-26 12/600 2026-03-31 10:01 by cal0306
[考研] 342求调剂 +4 加油a李zs 2026-03-26 4/200 2026-03-30 16:39 by 晶体之美
[考研] 298求调剂 +3 什么是胖头鱼 2026-03-30 5/250 2026-03-30 14:41 by 青海小西牛
[考研] 求调剂 +7 青春裁为三截 2026-03-29 7/350 2026-03-30 13:14 by laoshidan
[考研] 0703本科郑州大学求调剂 +7 nhj_ 2026-03-25 7/350 2026-03-30 12:44 by fangnagu
[考研] 总分293求调剂 +8 加一一九 2026-03-25 11/550 2026-03-29 19:53 by 无际的草原
[考研] 327求调剂 +6 汲亦昊 2026-03-29 6/300 2026-03-29 13:40 by peike
[考研] 0856求调剂 +13 zhn03 2026-03-25 14/700 2026-03-29 08:13 by fmesaito
[考研] 11408软件工程求调剂 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 21:50 by zhq0425
[考研] 本科新能源科学与工程，一志愿华理能动285求调剂 +7 AZMK 2026-03-28 11/550 2026-03-28 21:01 by xxxsssccc
[考研] 本科新能源科学与工程，一志愿华理能动285求调剂 +3 AZMK 2026-03-27 5/250 2026-03-28 16:19 by xxxsssccc
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-26 3/150 2026-03-26 19:57 by nihaoar
[考研] 296求调剂 +4 汪！？！ 2026-03-25 7/350 2026-03-25 16:41 by 汪！？！
信息提示
关闭
请填处理意见
关闭