在 fluent 的两相流模型 VOF 中,如果利用 UDF 计算物质传递(mass transfer),UDF 中给出的例子为:
/* UDF to define a simple mass transfer based on Saturation
Temperature */
DEFINE_EXCHANGE_PROPERTY(liq_gas_source, cell, thread, liq_i, gas_i)
{
real m_lg;
real T_SAT = 373.15.;
Thread *gas = THREAD_SUB_THREAD(thread, gas_i);
Thread *liq = THREAD_SUB_THREAD(thread, liq_i);
m_lg = 0.;
if (C_T(cell, liq) >= T_SAT)
{
m_lg = -0.1*C_VOF(cell,liq)*C_R(cell,liq)*
fabs(C_T(cell,liq)-T_SAT)/T_SAT;
}
if ((m_lg == 0. ) && (C_T(cell, gas) <= T_SAT))
{
m_lg = 0.1*C_VOF(cell,gas)*C_R(cell,gas)*
fabs(T_SAT-C_T(cell,gas))/T_SAT;
}
return (m_lg);
}
我想请问下,按照 help 上的说法,这里的 m_lg 是 mass transfer coenficient 的意思,但所谓的 mass transfer coenficient 又是什么意
思?
并且,m_lg 代表的是从液体到气体的蒸发还是气体到液体的凝固?
十分感谢。
2. 这是一个粗糙的相变传质公式,相传是 1975 年由华人科学家:李文和提出。显然,在这个 udf 中 m_lg 是对液相来说的:
if (C_T(cell, liq) >= T_SAT)
{
m_lg = -0.1*C_VOF(cell,liq)*C_R(cell,liq)*
fabs(C_T(cell,liq)-T_SAT)/T_SAT;
}
先判断 cell 内的液体温度是不是高于饱和温度 T_SAT,大的话就发生蒸发系数为负;
if ((m_lg == 0. ) && (C_T(cell, gas) <= T_SAT))
{
m_lg = 0.1*C_VOF(cell,gas)*C_R(cell,gas)*
fabs(T_SAT-C_T(cell,gas))/T_SAT;
}
如果液体的温度没有达到饱和,并且 cell 内气体的温度低于饱和温度,则发生凝结,液相源相为正。则 m_lg 为正。
公式中的 0.1 是为了增加计算过程的收敛性而加入的时间松弛系数。
评论0