| 查看: 2239 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
想模拟单个气泡的溃灭，想用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里面设置行吗，还是必须要编程？ |
回复此楼» 猜你喜欢
271材料工程求调剂
已经有5人回复
-
281求调剂（0805）
已经有16人回复
-
304求调剂
已经有6人回复
-
材料工程专硕调剂
已经有6人回复
-
一志愿天大材料与化工（085600）总分338
已经有4人回复
-
085700资源与环境308求调剂
已经有3人回复
-
求材料调剂
已经有8人回复
-
294求调剂材料与化工专硕
已经有5人回复
-
一志愿华中科技大学，080502，354分求调剂
已经有4人回复
-
一志愿吉林大学材料学硕321求调剂
已经有6人回复
3楼2016-04-25 16:41:57
jiachaolis
铜虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 184.1
- 红花: 2
- 帖子: 134
- 在线: 51.5小时
- 虫号: 3307685
- 注册: 2014-07-04
- 性别: GG
- 专业: 动力学与控制
2楼2016-04-24 20:01:44
4楼2016-04-26 16:53:41
wolfghost
版主 (知名作家)
- 博学EPI: 432
- 应助: 101 (高中生)
- 贵宾: 7.987
- 金币: 36642.2
- 散金: 11283
- 红花: 198
- 沙发: 8
- 帖子: 7788
- 在线: 1298小时
- 虫号: 573180
- 注册: 2008-06-13
- 性别: GG
- 专业: 电化学
- 管辖: 有奖问答
5楼2016-05-25 13:50:12
