UDF水蒸汽冷凝初始化时出现fatal error

279514060

#include "udf.h"
#include "sg_mphase.h"
#define A1_P 58042
#define A2_P 273
#define A3_P 2.963
#define A4_P 58660
#define A5_P 0.0151
#define A6_P 4.273
#define A7_P 2.963
#define A8_P 0.0151
#define A9_P 4.237
#define C1_DOT 2.022
#define C2_DOT -0.498
#define C3_DOT 0.666
#define C4_DOT 3.471
#define C5_DOT -2.223
#define C6_DOT 0.114
#define VOF_RY 6.278e-05
#define A10_P 101325
#define A11_P 0.662
#define LAT_HT 2260

DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn)

{

  Thread *mix_th, *sec_th, *thi_th;
  real m_dot_l;
  real nd_liq;
  real tem_liq,rot_liq,rot_vap,rot_air,tem_air;
  real d,pp_liq,pp_vap;
  real vof_liq,vof_vap,vof_air;

  mix_th = THREAD_SUPER_THREAD(pri_th);
  sec_th = THREAD_SUB_THREAD(mix_th, 1);

  vof_liq = C_VOF(cell,pri_th);
  vof_vap = C_VOF(cell,sec_th);
  vof_air = C_VOF(cell,thi_th);
  tem_liq = C_T(cell,pri_th);
  rot_liq = C_R(cell,pri_th);
  rot_vap = C_R(cell,sec_th);
  rot_air = C_R(cell,thi_th);
  tem_air = C_T(cell,thi_th);


  nd_liq =(7/13)*0.005/(vof_liq+(7/13)*0.005);
  d = (rot_vap*vof_vap)/(rot_air*vof_air);

  pp_liq = A1_P+A9_P*pow((tem_liq-A2_P),A3_P)-A4_P*pow(nd_liq,A8_P)-A6_P*pow((tem_liq-A2_P),A3_P)*pow(nd_liq,A8_P);
  pp_vap = (d*A10_P)/(d+A11_P);


  if(pp_vap > pp_liq){
              m_dot_l = 0.001*3697*pow(d,C1_DOT)*pow(tem_air,C2_DOT)*pow(1.1*(rot_vap+rot_air)*0.04,C3_DOT)*pow(nd_liq,C4_DOT)*pow(tem_liq,C5_DOT)*pow(VOF_RY,C6_DOT);
              dS[eqn] = -0.001*3697*pow(d,C1_DOT)*pow(tem_air,C2_DOT)*pow(1.1*(rot_vap+rot_air)*0.04,C3_DOT)*C4_DOT*pow(nd_liq,(C4_DOT-1))*(-(7/13)*0.005/pow((vof_liq+(7/13)*0.005),2))*pow(tem_liq,C5_DOT)*pow(VOF_RY,C6_DOT);
              }
              
  else{
             m_dot_l = -0.001*3697*pow(d,C1_DOT)*pow(tem_air,C2_DOT)*pow(1.1*(rot_vap+rot_air)*0.04,C3_DOT)*pow(nd_liq,C4_DOT)*pow(tem_liq,C5_DOT)*pow(VOF_RY,C6_DOT);
             dS[eqn] = -0.001*3697*pow(d,C1_DOT)*pow(tem_air,C2_DOT)*pow(1.1*(rot_vap+rot_air)*0.04,C3_DOT)*C4_DOT*pow(nd_liq,(C4_DOT-1))*(-(7/13)*0.005/pow((vof_liq+(7/13)*0.005),2))*pow(tem_liq,C5_DOT)*pow(VOF_RY,C6_DOT);
    }
   
  return m_dot_l;

}




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

{

  Thread *mix_th, *pri_th, *thi_th;
  real m_dot_v;
  real nd_liq;
  real tem_liq,rot_liq,rot_vap,rot_air,tem_air;
  real d,pp_liq,pp_vap;
  real vof_liq,vof_vap,vof_air;

  mix_th = THREAD_SUPER_THREAD(sec_th);
  pri_th = THREAD_SUB_THREAD(mix_th, 0);

  vof_liq = C_VOF(cell,pri_th);
  vof_vap = C_VOF(cell,sec_th);
  vof_air = C_VOF(cell,thi_th);
  tem_liq = C_T(cell,pri_th);
  rot_liq = C_R(cell,pri_th);
  rot_vap = C_R(cell,sec_th);
  rot_air = C_R(cell,thi_th);
  tem_air = C_T(cell,thi_th);


  nd_liq = (7/13)*0.005/(vof_liq+(7/13)*0.005);
  d = (rot_vap*vof_vap)/(rot_air*vof_air);

  pp_liq = A1_P+A9_P*pow((tem_liq-A2_P),A3_P)-A4_P*pow(nd_liq,A8_P)-A6_P*pow((tem_liq-A2_P),A3_P)*pow(nd_liq,A8_P);
  pp_vap = (d*A10_P)/(d+A11_P);

  if(pp_vap > pp_liq){
                   m_dot_v = -0.001*3697*pow(d,C1_DOT)*pow(tem_air,C2_DOT)*pow(1.1*(rot_vap+rot_air)*0.04,C3_DOT)*pow(nd_liq,C4_DOT)*pow(tem_liq,C5_DOT)*pow(VOF_RY,C6_DOT);
                   dS[eqn] = -0.001*3697*C1_DOT*pow(d,(C1_DOT-1))*(rot_vap)/(rot_air*vof_air)*pow(tem_air,C2_DOT)*pow(1.1*(rot_vap+rot_air)*0.04,C3_DOT)*pow(nd_liq,C4_DOT)*pow(tem_liq,C5_DOT)*pow(VOF_RY,C6_DOT);

                            }
              
      else{
                   m_dot_v = 0.001*3697*pow(d,C1_DOT)*pow(tem_air,C2_DOT)*pow(1.1*(rot_vap+rot_air)*0.04,C3_DOT)*pow(nd_liq,C4_DOT)*pow(tem_liq,C5_DOT)*pow(VOF_RY,C6_DOT);
                   dS[eqn] = -0.001*3697*C1_DOT*pow(d,(C1_DOT-1))*(rot_vap)/(rot_air*vof_air)*pow(tem_air,C2_DOT)*pow(1.1*(rot_vap+rot_air)*0.04,C3_DOT)*pow(nd_liq,C4_DOT)*pow(tem_liq,C5_DOT)*pow(VOF_RY,C6_DOT);
                   }
  return m_dot_v;

}



DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn)

{
  Thread *pri_th , *sec_th, *thi_th;
  real m_dot, m_dot_e;
  real nd_liq;
  real tem_liq,rot_liq,rot_vap,rot_air,tem_air,tem_vap;
  real d,pp_liq,pp_vap;
  real vof_liq,vof_vap,vof_air;

  pri_th = THREAD_SUB_THREAD(mix_th, 0);
  sec_th = THREAD_SUB_THREAD(mix_th, 1);
  thi_th = THREAD_SUB_THREAD(mix_th, 2);


  vof_liq = C_VOF(cell,pri_th);
  vof_vap = C_VOF(cell,sec_th);
  vof_air = C_VOF(cell,thi_th);
  tem_liq = C_T(cell,pri_th);
  tem_vap = C_T(cell,sec_th);
  rot_liq = C_R(cell,pri_th);
  rot_vap = C_R(cell,sec_th);
  rot_air = C_R(cell,thi_th);
  tem_air = C_T(cell,thi_th);


  nd_liq = (7/13)*0.005/(vof_liq+(7/13)*0.005);
  d = (rot_vap*vof_vap)/(rot_air*vof_air);

  pp_liq = A1_P+A9_P*pow((tem_liq-A2_P),A3_P)-A4_P*pow(nd_liq,A8_P)-A6_P*pow((tem_liq-A2_P),A3_P)*pow(nd_liq,A8_P);
  pp_vap = (d*A10_P)/(d+A11_P);



    if(pp_vap > pp_liq){

         m_dot = 0.001*3697*pow(d,C1_DOT)*pow(tem_air,C2_DOT)*pow(1.1*(rot_vap+rot_air)*0.04,C3_DOT)*pow(nd_liq,C4_DOT)*pow(tem_liq,C5_DOT)*pow(VOF_RY,C6_DOT);
         dS[eqn] = -0.001*3697*pow(d,C1_DOT)*pow(tem_air,C2_DOT)*pow(1.1*(rot_vap+rot_air)*0.04,C3_DOT)*pow(nd_liq,C4_DOT)*C5_DOT*pow(tem_liq,(C5_DOT-1))*pow(VOF_RY,C6_DOT);

                             }
   
      else {

         m_dot = -0.001*3697*pow(d,C1_DOT)*pow(tem_air,C2_DOT)*pow(1.1*(rot_vap+rot_air)*0.04,C3_DOT)*pow(nd_liq,C4_DOT)*pow(tem_liq,C5_DOT)*pow(VOF_RY,C6_DOT);
         dS[eqn] = -0.001*3697*pow(d,C1_DOT)*C2_DOT*pow(tem_air,(C2_DOT-1))*pow(1.1*(rot_vap+rot_air)*0.04,C3_DOT)*pow(nd_liq,C4_DOT)*pow(tem_liq,C5_DOT)*pow(VOF_RY,C6_DOT);

             }
    return  m_dot_e = m_dot*LAT_HT;

}
求助：我做的是三相的传热传质，水、水蒸汽、干空气，水蒸汽发生冷凝凝结成水，水蒸汽放出汽化潜热！
我用的混合模型！可是还是在初始化的过程中出现问题了
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()
编译通过，但是初始化不行，求助：我查了好久法相网上人都说是变量御值方面出了问题，请大家给些意见和讨论，谢谢

leider2010

回复 279514060 的帖子

我也正在做相关的题目，大家可以交流一下。

不过不知道为什么是3相，不就是气液两相么？我觉得水蒸汽可以按照一个气相组分处理。

Heidiadalheid

279514060 发表于 2010-8-3 15:35

                               
登录/注册后可看大图

#include "udf.h"
#include "sg_mphase.h"
#define A1_P 58042

这个错误信息好像不是针对某一个具体的特定错误的，他不能说明任何问题
他可能是计算机配置的问题，内存不足，有的人重启计算机就好使了，也有的人换了计算机就好使了
也有可能是模型的问题，需要重新划分网格
当然也不能排除udf的使用问题，看看你的头文件是否能够正确的调用

Heidiadalheid

leider2010 发表于 2010-10-8 18:14

                               
登录/注册后可看大图

回复 279514060 的帖子

我也正在做相关的题目，大家可以交流一下。

湿空气中的水蒸气在一定条件下有相变，所以作为一个气相组分来处理不好处理

h04581803