初来乍到，关于UDF有一事不明

mack1983

[这个贴子最后由mack1983在 2005/09/25 04:18pm 第 1 次编辑]

在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代表的是从液体到气体的蒸发还是气体到液体的凝固？

十分感谢。

Joywe

这是一个粗糙的相变传质公式，相传是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是为了增加计算过程的收敛性而加入的时间松弛系数。
最后，为什么说m_lg是系数呢？你仔细看看它的单位就会发现它仅仅是个系数，源相的话还要乘以其质量。

沃尔玛69

有没有关于这个公式出处的参考文献？谢谢！