24小时热门版块排行榜

>论坛更新日志 (3067)
>虫友互识 (493)
>导师招生 (172)
>文献求助 (154)
>休闲灌水 (82)
>硕博家园 (80)
>考博 (37)
>博后之家 (35)
>基金申请 (29)
>论文投稿 (28)
>找工作 (25)
>教师之家 (22)
>公派出国 (22)
>考研 (21)
>论文道贺祈福 (17)
>外文书籍求助 (17)

返回列表

烟雨笑清风

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 510.2
散金: 96
红花: 2
帖子: 172
在线: 100.8小时
虫号: 1396938
注册: 2011-09-09
专业: 化工系统工程

[求助] FLUENT received fatal signal (ACCESS_VIOLATION) 初始化出错

我进行颗粒的湿度模拟，导入udf发现不能计算，初始化后就出现以下错误，然后我初始化后加载udf，可以初始化但是计算又出现以下错误。请各位大神给解决一下，在此先谢谢了。
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.
以下是我的udf
#include "udf.h"
#define PSAT_A 0.01
#define PSAT_TP 338.15
#define C_LOOP 8
#define H2O_PC 22.089E6
#define H20_TC 647.286
/*用户定义变量*/
#define MAX_SPE_EQNS_PRIM 2 /*主相的组分数*/
#define MAX_SPE_EQNS_SEC 2 /*次相的组分数*/
#define  index_evap_primary 0/*主相蒸发组分*/
#define  index_evap_sec 1/*次相蒸发组分*/
#define  prim_index 0/*主相*/
#define  P_OPER 101325/*操作压力为常压*/

/*结束输入*/

double psat_h2o(double tsat)
/*计算水的饱和蒸汽压*/
/*作为温度的函数*/
/*返回值单位为PASCALS*/
{
      int i;
      double  var1,sum1,ans1,psat;
      double  constants[8]={-7.4192420,2.97221E-1,-1.155286E-1,8.68563E-3,1.094098E-3,-4.39993E-3,2.520658E-3,-5.218684E-4};

/* var1用于循环*/

  var1=PSAT_A*(tsat-PSAT_TP);
/*计算总循环*/
i=0;
sum1=0.0;
while (i<C_LOOP)
{sum1+=constants*pow(var1,i);++i;}
ans1=sum1*(H20_TC/tsat-1.0);
/*psat的单位是Pascals*/
psat=H2O_PC*exp(ans1);
return psat;
}

DEFINE_HET_RXN_RATE(user_evap_condens_react,c,t,hr,mw,yi,rr,rr_t)
{
Thread **pt=THREAD_SUB_THREADS(t);
Thread *tp=pt [0];
Thread *ts=pt[1];
int i;
real concentration_evap_primary,accum=0.,mole_frac_evap_prim,concentration_sat;
real  T_prim=C_T(c,tp);/*气相温度*/
real  T_sec=C_T(c,ts);/*颗粒相温度*/
real  diam=C_PHASE_DIAMETER(c,ts);/*颗粒相直径*/
real  D_evap_prim=C_DIFF_EFF(c,tp,index_evap_primary)-0.7*C_MU_T(c,tp)/C_R(c,tp);
      /*主相的湍流扩散率*/
real Re,Sc,Nu,urel,urelx,urely,urelz=0.,mass_coeff,area_density,flux_evap;
if(Data_Valid_P())/*判断是否获得数据*/
{
   urelx = C_U(c,tp)-C_U(c,ts);
   urely = C_V(c,tp)-C_V(c,ts);
  #if RP_3D
   urelz = C_W(c,tp)-C_W(c,ts);

#endif

   urel = sqrt(urelx*urelx+urely*urely+urelz*urelz);
      /*relative velocity*/

      Re = urel*diam*C_R(c,tp)/C_MU_T(c,tp);

      Sc = C_MU_T(c,tp)/C_R(c,tp)/D_evap_prim ;
      Nu = 2.+0.6*pow(Re,0.5)*pow(Sc,0.333);
      mass_coeff = Nu*D_evap_prim/diam;  /*传质系数*/
   for (i=0;i<MAX_SPE_EQNS_PRIM;i++)
      {
         accum = accum+C_YI(c,tp,i)/mw[prim_index];;
      }
   mole_frac_evap_prim = C_YI(c,tp,index_evap_primary )/mw[index_evap_primary][prim_index]/accum;
      concentration_evap_primary = mole_frac_evap_prim*P_OPER/UNIVERSAL_GAS_CONSTANT/T_prim;
      concentration_sat = psat_h2o(T_sec)/UNIVERSAL_GAS_CONSTANT/T_sec;
      area_density = 6.*C_VOF(c,ts)/diam ;
      flux_evap = mass_coeff *(concentration_sat - concentration_evap_primary );
if(0.002< yi[1][1]<0.125)
   {*rr=(yi[1][1]-0.125)/(0.125-0.002)*area_density*flux_evap;}
else
*rr = area_density*flux_evap;
   }
}