小女子请教各位大虾们，有溶液蒸发过程的UDF算例吗？

shangyiyin

        来请各位高手帮忙。我做的模拟是往煤粉炉炉膛里喷尿素溶液，用双欧拉做烟气和尿素溶液的混合。然后要编一个UDF模拟尿素溶液的蒸发和化学反应的过程。请问有蒸发的UDF的例子吗？如果有相似的化学反应的UDF例子也一并给我。在这里先谢谢了！
我的邮箱是：shangyiyin@126.com
比较急，万分感谢！

[ 本帖最后由 shangyiyin 于 2010-5-24 10:32 编辑 ]

shangyiyin

怎么没有人回复我呀？O(∩_∩)O谢谢大家呀。

limingyiban

#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_1;
  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_1=-0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,pri_th)-T_SAT)/T_SAT;
       dS[eqn]=-0.1*C_R(cell,pri_th)*fabs(C_T(cell,pri_th)-T_SAT)/T_SAT;
    }
  else
    {  m_dot_1=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
       dS[eqn]=0.;
    }
    return m_dot_1;
}

DEFINE_SOURCE(vap_src,cell,sec_th,dS,eqn)
{
  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[eqn]=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[eqn]=-0.1*C_R(cell,sec_th)*fabs(C_T(cell,sec_th)-T_SAT)/T_SAT;
    }
    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);
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,pri_th)-T_SAT)/T_SAT;
     dS[eqn]=-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[eqn]=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)/T_SAT;
  }
return LAT_HT*m_dot;
}

shangyiyin

热心的朋友，非常感谢呀！我现在在用这个例子。但是我模拟的是多相流多组分，其中第一相有水蒸气、O2、NO、N2，第二相有尿素（液态）、水（液态）。在这里用上面的那个蒸发的udf不行。因为它只有水和水蒸气。需要怎么修改呀？在这种情况下还可以用SOURCE宏吗？
可否加我QQ:546515188 指教一二