24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xxkobett

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 14.8
散金: 10
帖子: 13
在线: 3.2小时
虫号: 5129008
注册: 2016-10-17
专业: 工程热物理相关交叉领域

[求助] 气液相变模型，可是每次初始化都提示Segmentation Fault，新人求解答

#include "udf.h"
#include "sg_mphase.h"
#include "mem.h"
#define T_SAT 373.15
#define LAT_HT 2.257e6

DEFINE_SOURCE(vap_src, cell, pri_th, dS, eqn)
{

  Thread *mix_th, *sec_th;
  real m_dot_v;

  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_v = -2*C_K_L(cell,sec_th)*NV_DOT(C_VOF_G(cell,sec_th),C_T_G(cell,mix_th))/LAT_HT;

dS[eqn] = 0;
}
  else
{
m_dot_v = 2*C_K_L(cell,pri_th)*NV_DOT(C_VOF_G(cell,pri_th),C_T_G(cell,mix_th))/LAT_HT;

dS[eqn] = 0;

   }

  return m_dot_v;

}

DEFINE_SOURCE(lip_src, cell, sec_th, dS, eqn)

{

  Thread * mix_th, *pri_th;
  real m_dot_l;

  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_l = 2*C_K_L(cell,sec_th)*NV_DOT(C_VOF_G(cell,sec_th),C_T_G(cell,mix_th))/LAT_HT;
dS[eqn] = 0;
                           }
  else
{
m_dot_l = -2*C_K_L(cell,pri_th)*NV_DOT(C_VOF_G(cell,pri_th),C_T_G(cell,mix_th))/LAT_HT;

dS[eqn] = 0;

   }

  return m_dot_l;

}

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 = 2*C_K_L(cell,sec_th)*NV_DOT(C_VOF_G(cell,sec_th),C_T_G(cell,mix_th))/LAT_HT;

dS[eqn] = 0;
                           }

  else {
m_dot = -2*C_K_L(cell,pri_th)*NV_DOT(C_VOF_G(cell,pri_th),C_T_G(cell,mix_th))/LAT_HT;

dS[eqn] = 0;
}

  return LAT_HT*m_dot;
}@wuming524

回复此楼