求大神看看结冰UDF是否正确

lj675230266

#include "udf.h"
  DEFINE_MASS_TRANSFER(liq_sol_source, cell, thread, from_index,
  from_species_index, to_index, to_species_index)
  {
    real m_ls;
    real T_SAT = 273.15;/*凝固温度*/
    Thread *liq, *sol;
    liq = THREAD_SUB_THREAD(thread, from_index);
    sol = THREAD_SUB_THREAD(thread, to_index);
    m_ls = 0.0;

    if (C_T(cell, liq) < T_SAT)
      {                                               /* soliding 凝固*/
        m_ls = 0.1*C_VOF(cell,liq)*C_R(cell,liq)*
               (C_T(cell,liq)-T_SAT)/T_SAT;
      }
    else if (C_T(cell, sol) > T_SAT)
      {                                               /* thawing 融化*/
        m_ls = -0.1*C_VOF(cell,sol)*C_R(cell,sol)*
               (T_SAT-C_T(cell,sol))/T_SAT;
      }

    return (m_ls);
  }

易阳

大神好，求助模拟60度热水在敞口绝热容器里，大范围温度为零下20度时，热水结冰情况，可以有偿，可以共同研究。

巢老根

这里的m-ls是什么物理意义，为什么它可以代表物理状态的变化呢？

巢老根

这里的m-ls是什么物理意义，为什么它可以代表物理状态的变化呢？

巢老根

这里的m-ls是什么物理意义，为什么它可以代表物理状态的变化呢？

lionel123

易阳 发表于 2017-12-29 11:19
大神好，求助模拟60度热水在敞口绝热容器里，大范围温度为零下20度时，热水结冰情况，可以有偿，可以共同 ...

还在做水结冰吗?希望加Q1272137826交流