沸腾换热的UDF程序讨论!

jielimi

  本人在用沸腾UDF程序时,出现问题,算了很多次结果和理论的相差很大,请各位高手指教!
我用的程序也是本网站发表的,可以肯定的时候里面错误不少,下面是按照我的理解修改的,只是一部分,不知道理解对不对!请大家说说看!
#include "udf.h"  
#include "sg_mphase.h"  
#define T_SAT 373   
#define LAT_HT 2.253e+6 (水潜热）
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 = -4220*C_VOF(cell, pri_th)*C_R(cell, pri_th)*   
                fabs(C_T(cell, mix_th) - T_SAT)/ LAT_HT; （4220为水的比热，LAT_HT为水潜热）

    dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
              （dS这相我不明白什么意思，但估计前面的0.1也应该换成4220）
                               }
  else { 应该是气态转换成液体，需要计算气态所含多余的热量}
    m_dot_l = 2028*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,sec_th))/ LAT_HT;
              （2028蒸汽比热）
    dS[eqn] = 0.; （不明白什么意思）   
       }
   
  return m_dot_l;
}
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);  mix-th指针代表气相
  pri_th = THREAD_SUB_THREAD(mix_th, 0);   pri-th指针代表液相
  if(C_T(cell, mix_th)>=T_SAT) 如果气相温度大于373（估计这里需要改成液相的）
{
    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; 计算转换成气态的水，公式应该不对，后面一些也是这样的问题，

jielimi

     有知道的请指教! 本人的QQ:41957527

crystalshine

Thread *mix_th, *sec_th;  
real m_dot_l;
这几个变量代表什么？

ch06

Thread *mix_th, *sec_th; 代表定义线指针
*mix_th下面给出是指向液相
*sec_th指向气相
real m_dot_l 定义液相汽化时转移的质量

pvxd

有几个问题你理解应该存在偏差吧？这个程序里面用到的沸腾公司非常简单
公式前面那个系数，0.1是指的蒸发系数和凝结系数。虽然数值上个人觉得并不正确

本人在用沸腾UDF程序时,出现问题,算了很多次结果和理论的相差很大,请各位高手指教!
我用的程序也是本网站发表的,可以肯定的时候里面错误不少,下面是按照我的理解修改的,只是一部分,不知道理解对不对!请大家说说看!
#include "udf.h"  
#include "sg_mphase.h"  
#define T_SAT 373   
#define LAT_HT 2.253e+6 (水潜热）
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 = -4220*C_VOF(cell, pri_th)*C_R(cell, pri_th)*   
               fabs(C_T(cell, mix_th) - T_SAT)/ LAT_HT; （4220为水的比热，LAT_HT为水潜热）
   dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
             （dS这相我不明白什么意思，但估计前面的0.1也应该换成4220）
                              }
else { 应该是气态转换成液体，需要计算气态所含多余的热量}
   m_dot_l = 2028*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
               fabs(T_SAT-C_T(cell,sec_th))/ LAT_HT;
             （2028蒸汽比热）
   dS[eqn] = 0.; （不明白什么意思）   
      }
   
return m_dot_l;
}
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);  mix-th指针代表气相
pri_th = THREAD_SUB_THREAD(mix_th, 0);   pri-th指针代表液相
if(C_T(cell, mix_th)>=T_SAT) 如果气相温度大于373（估计这里需要改成液相的）
{
   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; 计算转换成气态的水，公式应该不对，后面一些也是这样的问题，

ch06

这个模型好像来自
Lee, W. H. (1979). A Pressure Iteration Scheme for Two-Phase Flow Modeling (Technical Paper No. LA-UR-79-975). Los Alamos, New Mexico, USA: Los Alamos National Laboratory.
我是在一篇英文文献的参考文献中见到的，但一直找不到上面的文献，
如果有找到的希望共享一下

wf0609

这个程序中的 几个温度的相对变化没搞明白