沸腾换热

slz361

求助：水的沸腾换热的udf

tingyue2008

二、传热
#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_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 = -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_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,mix_th))/T_SAT;

tingyue2008

我也是在别的地方找的  希望能帮到你

slz361

高手，不胜感激。

tingyue2008

里面水蒸气的潜热和饱和温度根据你的实际情况定
我现在也搞不懂这个UDf  里面单位都不一致  不知道能用不呢

还有0.1 怎么来的  液体温度减去饱和温度的绝对值/饱和温度这个又是什么意思
知道的请教一下？？？

浪子轻狂

这个udf还是相当粗糙，练练手还可以。

0.1是蒸发速度，是一个相当复杂的系数，和很多因素有关，在相变过程中，实际上是会发生变化的

lsq

我的毕设任务好像就是研究这个0.1的问题

浪子轻狂

够你研究的，呵呵