| 查看: 869 | 回复: 1 | |||
[交流]
【求助】用species transport模型模拟蒸发UDF的疑问
|
|
我用某教材里的下面这个UDF模拟130摄氏度（或更高）的热空气加热80摄氏度的水，热空气采用组分运输模型，多项流采用VOF模型（因为有自由页面），含有氧气、氮气、水蒸气。 但是模拟的结果空气中产生不出水蒸气，而用FLUENT自带的沸腾蒸发模型可以产生。 现在初步判断可能是UDF中几个指针有问题，但是不知道怎么改，或者是源项加载的地方不对。 其中，air_src加给空气中的H2O组分，liq_src加给水，enrg_src加给混合物 麻烦各位帮助，谢谢！！！ #include "udf.h" //包括常规宏 #include "sg_mphase.h" #define T_SAT 373 //定义蒸发温度100℃ #define LAT_HT 2.3e6 //定义蒸发潜热J/Kg DEFINE_SOURCE(air_src,cell,pri_th,dS,eqn) //气相质量源项UDF { Thread * mix_th, *sec_th; real m_dot_a; mix_th = THREAD_SUPER_THREAD(pri_th); //指向混合区的第二相即气相的指针 sec_th = THREAD_SUB_THREAD(mix_th,1); //指向单相控制区的液相的指针，液相为第二相 if(C_T(cell, mix_th) >= T_SAT) //如果混合区单元的温度高于蒸发温度，液相向气相的质量质量转移 { m_dot_a = 0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn] = 0; //由于是液相向气相转移，所以气相的质量源项对来自液相的质量转移的偏导为零 } else { m_dot_a = -0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; //如果指向混合区的单元温度小于蒸发温度，气相向液相的质量转移，气相失 dS[eqn] = -0.1*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; //由于是气相向液相转移，所以气相的质量源项对自身的质量转移的偏导不为零 } return m_dot_a; } DEFINE_SOURCE(liq_src,cell,sec_th,dS,eqn) //液相质量源项UDF { Thread *mix_th, *pri_th; //定义计算区线指针 real m_dot_l; //定义液相质量转移 kg/(m2.s) mix_th = THREAD_SUPER_THREAD(sec_th); //指向混合区的主相即液相的指针 pri_th = THREAD_SUB_THREAD(mix_th,0); //指向单相控制区的气相的指针，气相为主相 if(C_T(cell, mix_th) >= T_SAT) //如果液相单元的温度高于蒸发温度，液相向气相的质量质量转移 { m_dot_l = -0.1*C_VOF(cell,sec_th)*C_R(cell, sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn] = -0.1*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; //定义源项对质量转移偏导 } else { m_dot_l = 0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; //如果指向混合区液相的单元温度小于蒸发温度，气相向液相的质量转移，液相得 dS[eqn] = 0; //由于是气相向液相转移，所以液相的质量源项对质量转移的偏导为零 } return m_dot_l; } DEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn) //混合模型能量源项UDF { Thread *pri_th, *sec_th; real m_dot; pri_th = THREAD_SUB_THREAD(mix_th,0); //指向混合区的液相的指针 sec_th = THREAD_SUB_THREAD(mix_th,1); //指向混合区的气相的指针 if(C_T(cell, mix_th)>=T_SAT) //如果混合区的单元温度高于蒸发温度。质量转移由液相向气相转移，吸热，质量转移量前有负号 { m_dot = -0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn] = -0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)/T_SAT; } else { //相反，气相向液相转移则放热 m_dot = 0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn] = 0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)/T_SAT; } return LAT_HT*m_dot; //气化潜热与质量转移率相乘得能量源项W/m3 } [ Last edited by kobe8875 on 2013-10-18 at 14:31 ] |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
VOF模型静止壁面上速度不为0
已经有10人回复
- 哪位大侠会fluent定义化学反应
已经有3人回复
- FLUENT加入VOF模型后的一些问题
已经有15人回复
- 用UDF定义了volumetric reaction rate 还怎么在fluent里设置阿累尼乌斯定律参数
已经有11人回复
- 采用VOF模型对磨削加工中的流场进行仿真，仿真结果不对，摸索中，望高手相助！！
已经有28人回复
- fluent中的species transport中如何正确输入反应参数
已经有5人回复
- ansys13.0启动fluent出错
已经有30人回复
- udf获得多相流中某一项中某一组分的含量
已经有14人回复
- Fluent，VOF模型，怎么在追踪的界面上施加反冲压
已经有7人回复
- UDF的编译通不过
已经有8人回复
- vof模型和species transport模型到底冲突不
已经有18人回复
- fluent-udf自定义标量方程边界条件
已经有5人回复
- 在气相混合中如何选择正确的物理模型
已经有7人回复
- fluent 关于prePDF模型求助
已经有13人回复
- 组分输运模型的出口组分设置？
已经有26人回复
- 【求助】通气空泡多相流（气、水蒸汽、水）空化模型的UDF如何入手啊？？？？
已经有4人回复
- 【求助】求UDF高手帮忙看看我的问题
已经有18人回复
- 【求助】用fluent模拟超重力气液流体动力模型和传质过程，中间有化学吸收反应，fluent
已经有23人回复
2楼2013-10-17 12:20:08
