请问:下面这个蒸发过程UDF里传质公式的出处?
水－水蒸汽蒸发过程UDF#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 373
#define LAT_HT 1.e3
DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn)
{
Thread *mix_th, *sec_th;
real m_dot_l;
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_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
dS = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
else {
m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
fabs(T_SAT-C_T(cell,mix_th))/T_SAT; /*就是这个公式,请问这个公式的出处?谢谢*/
dS = 0.;
}
return m_dot_l;
}
[ 本帖最后由 quakefans 于 2009-6-15 09:57 编辑 ] 刚才发现,这个是UDF手册关于DEFINE_MASS_TRANSFER 一个例子里用到的公式,可惜也没有给出参考文献. 找到这个公式的出处了,但是文献太老了(Wen Ho Lee (1979)),没找到,哪位如果能下的到这篇文献的话,麻烦共享下.
文献的名字没有,只有个作者和年代,下面是公式的相关信息. 请把附件中的内容传全。
谢谢了。:handshake 原帖由 chbo76 于 2009-6-26 11:22 发表 http://www.chinavib.com/forum/images/common/back.gif
请把附件中的内容传全。
谢谢了。:handshake
是个UDF教程，教程里关于蒸发只有这点内容，而且回帖不支持PDF和RAR文件上传。 哎,什么事情还是自己最可靠啊,呵呵.这个公式的最原始的出处没找到,但找到了篇别人用这个公式的文章.有需要的,自己去下载吧.
Influence of surface evaporation on stratification in liquid hydrogen tanks of different aspect ratios,S. Prasanth Kumar, B.V.S.S.S. Prasad, G. Venkatarathnam, K. Ramamurthi, S. Srinivasa Murthy,International Journal of Hydrogen Energy 32 (2007) 1954-1960 本帖最后由 Heidiadalheid 于 2010-10-9 17:06 编辑
我来吧这个udf补全吧
#include "udf.h"‘包括常规宏
#include "sg_mphase.h"‘包括体积分数宏CVOF(C,T)
#define T_SAT 373 ‘定义蒸发温度100℃
#define LAT_HT 1.e3 ‘定义蒸发潜热J/Kg
DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) ‘液相质量源项UDF
{
Thread *mix_th, *sec_th;‘定义计算区域线指针
real m_dot_l; ‘定义液相质量转移 kg/(m2.s)
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_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
dS = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
‘ 定义源项对质量转移偏导，
}
else {
m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
‘如果指向混合区液相的单元温度小于蒸发温度，气相向液相的质量转移，液相得
dS = 0.;‘由于是气相向液相转移，所以液相的质量源项对质量转移的偏导为零
}
return m_dot_l;
}
DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) ‘气相质量源项UDF
{
Thread * mix_th, *pri_th;
real m_dot_v;
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_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
dS = 0.; ‘由于是液相向气相转移，所以气相的质量源项对来自液相的质量转移的偏导为零
}
else {
m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
‘如果指向混合区的单元温度小于蒸发温度，气相向液相的质量转移
dS = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
由于是气相向液相转移，所以气相的质量源项对自身的质量转移的偏导不为零
}
return m_dot_v;
}
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, pri_th)*C_R(cell, pri_th)*
fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
dS =-0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT;
}
else { 相反，气相向液相转移则放热
m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
dS = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;}
return LAT_HT*m_dot;气化潜热与质量转移率相乘得能量源项W/m3
}
请问楼上的各位，你们有没有用这个UDF做过类似的模拟，结果如何？有没有实际的情况做过对比？谢谢~！ 我做过，趋势还可以。但是具体的量值还没法验证。没有实验对比。{:{23}:} 只能说还能用吧，定量的效果一般 这个UDF能否并行运行时使用呢？应该如何改一改，能并行运行？ jimee_lzj 发表于 2010-10-12 20:06 static/image/common/back.gif
这个UDF能否并行运行时使用呢？应该如何改一改，能并行运行？
这类问题我也没处理过，不过udf帮助文件里有这一章是介绍这方面工作的 另外我这儿有一个飞昂公司的ppt，最后部分涉及到这个问题，你可以参考一下