24小时热门版块排行榜

返回列表

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

俊采懿范

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 120.5
帖子: 10
在线: 9.1小时
虫号: 4402942
注册: 2016-02-15
专业: 工程热物理相关交叉领域

[求助] 参照论文写了湿空气结霜的udf，成功编译后初始化出错，恳请各位大神指点！！！已有3人参与

#include "udf.h"
#include "sg_mphase.h"

DEFINE_MASS_TRANSFER(liq_gas_source,cell,thread,from_phase_index,from_species_index,to_phase_index,to_species_index)
{
real hr;                //饱和空气含湿量//
real Wvs;                //饱和水蒸气质量分数//
real Wv;                //湿空气水蒸气质量分数//
real m_lg_m;                //质量转移//
real Pvs;                //水蒸气饱和压力//
      real Pv;             //水蒸气饱和分压力//
real Td;

Thread *tf= THREAD_SUPER_THREAD(thread,thread);
Thread *mix = THREAD_SUB_THREAD(thread,from_phase_index);
Thread *ice = THREAD_SUB_THREAD(thread,to_phase_index);

if(C_T(cell,tf)<=273.15 && C_T(cell,tf)>173.15)
{
Pvs=exp(-5674.5359/C_T(cell,tf)
+6.3925247
-9.677843*pow(10,-3)*C_T(cell,tf)
+6.22157*pow(10,-7)*C_T(cell,tf)*C_T(cell,tf)
+2.0747825*pow(10,-9)*C_T(cell,tf)*C_T(cell,tf)*C_T(cell,tf)
+9.484024*pow(10,-13)*C_T(cell,tf)*C_T(cell,tf)*C_T(cell,tf)*C_T(cell,tf)+
4.1635019*log(C_T(cell,tf)));
}
else if

   (C_T(cell,tf)<=473.15 && C_T(cell,tf)>273.15)
{
Pvs=exp(-5.8002206*pow(10,3)/C_T(cell,tf)
+1.3914993-4.8640239*pow(10,-2)*C_T(cell,tf)
-4.1764768*pow(10,-5)*C_T(cell,tf)*C_T(cell,tf)
-1.4452093*pow(10,-8)*C_T(cell,tf)*C_T(cell,tf)*C_T(cell,tf)+6.5459673*log(C_T(cell,tf)));
}
                                                                        /*定义计算单元温度对应的水蒸气饱和压力*/
hr=0.622*Pvs/(101325-Pvs);
Wvs=hr/(1+hr);
Wv=C_YI(cell,mix,0);
      Pv=(ABS_P(C_P(cell,mix),op_pres)*hr)/(hr+0.622);

      if
         (C_T(cell,tf)>=213.16&&C_T(cell,tf)<=273.15)
      {
         Td=(-60.45+7.0332*log(Pv)+0.37*log(Pv)*log(Pv))+273.15;
      }
      else if
         (C_T(cell,tf)>273.15&&C_T(cell,tf)<343.15)
         {
            Td=(-35.957-1.8726*log(Pv)+1.1689*log(Pv)*log(Pv))+273.15;
         }

if
(C_T(cell,tf)<=(Td<273.15?Td:273.15))
{
m_lg_m=C_VOF(cell,mix)*C_R(cell,mix)*Wv*(Wv-Wvs)/Wvs; /*水蒸气向冰相质量转移的速率*/

}
else if (C_T(cell,tf)>(Td<273.15?Td:273.15))
{
m_lg_m=0;
}
return (m_lg_m);

}

出错方式：
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: ()

模型是使用欧拉+组分运输模型进行，同时自定义添加质量运输方程，恳请各位大神帮忙