蒸发模拟,不是沸腾!!!
RT,诸位大神,我要模拟常温下的蒸发情况,不是沸腾情况,fluent自带的蒸发模型是沸腾的。
网上看了下,说蒸发是液体表面发生汽化,而沸腾是液体内部和表面都发生。
肯定是要编写UDF的,那么问题是该如何让只在表面发生相变?
然后我看网上流传的蒸发UDF,里面质量传递单位是 kg/m3·s,也就是说单位体积液体相变的质量,而我的情况是只在表面相变。
大神们,有没有什么建议,该从何处下手?
你的蒸发做的怎样了,我最近也在做着个,很是茫然
同问 楼主,您好,请问这个案例您最后是怎么完成的?我研究的是膜蒸馏过程中潜热蒸发仿真 楼主最后有结果了么? 同问,我也在迷茫 #include "udf.h"
#include "sg_mphase.h"
#define T_SAT 377 //饱和温度
#define VOF_SAT 0.5 //定义边界
#define LAT_HT 2.257e6//相变潜热1000j/kg
#define M_W 0.018 //水的摩尔质量18g/mol
#define g_c 8.314 //气体分子常数8.314j/(mol*k)
#define v_md 0.58979 //水蒸气的密度
//液体首相的传质
DEFINE_SOURCE(liq_src,cell,pri_th,dS,eqn)
{
Thread *mix_th, *sec_th;
real m_dot_l;
//mix_th
mix_th =THREAD_SUPER_THREAD(pri_th);
sec_th =THREAD_SUB_THREAD(mix_th,1);
//混合相温度导致相变
m_dot_l= -5*(C_T(cell,pri_th)-T_SAT)*v_md*LAT_HT*pow((1-C_VOF(cell,pri_th))*C_VOF(cell,pri_th),2)/(T_SAT*sqrt(2*3.14*g_c*T_SAT));
//对体积分数求导
dS= 0;
return m_dot_l;
}
//定义汽相的传质
DEFINE_SOURCE(vap_src,cell,sec_th,dS,eqn)
{
Thread *mix_th, *pri_th;
real m_dot_v;
//mix_th
mix_th =THREAD_SUPER_THREAD(sec_th);
pri_th =THREAD_SUB_THREAD(mix_th,0);
//混合相温度导致相变
m_dot_v= 5*(C_T(cell,pri_th)-T_SAT)*v_md*LAT_HT*pow((1-C_VOF(cell,pri_th))*C_VOF(cell,pri_th),2)/(T_SAT*sqrt(2*3.14*g_c*T_SAT));
//对体积分数求导
dS= 0;
return m_dot_v;
}
//能量传递方程
DEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn)
{
Thread *pri_th, *sec_th;
real m_dot;
pri_th =THREAD_SUB_THREAD(mix_th,0);
sec_th=THREAD_SUB_THREAD(mix_th,1);
//混合相温度导致相变
m_dot= -5*(C_T(cell,pri_th)-T_SAT)*v_md*LAT_HT*pow((1-C_VOF(cell,pri_th))*C_VOF(cell,pri_th),2)/(T_SAT*sqrt(2*3.14*g_c*T_SAT));
//对体积分数求导
dS= 0;
//热量的传递
return LAT_HT*m_dot;
} 我目前是这样做的,不知道有没有效果 名刀马 发表于 2023-3-2 18:39
#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 377 //饱和温度
您好您现在做的怎么样?可以交流吗?qq 2529693979 duyuxuan1997 发表于 2022-5-8 15:27
同问,我也在迷茫
您好 我也是做这个的 可以一起交流吗 qq 2529693979 名刀马 发表于 2023-3-2 18:39
#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 377 //饱和温度
您好 您这个UDF 中 公式是参考哪里的?可以请教一下嘛
页:
[1]