24小时热门版块排行榜    

查看: 950  |  回复: 0

烟雨笑清风

金虫 (小有名气)

[求助] 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;                                       
       }
}
回复此楼
爱拼才会赢
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 烟雨笑清风 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博发邮件 +5 Lmengk 2026-09-15 5/250 2026-09-17 17:55 by abcd..,
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 5BDX0d0WFp7t 2026-09-15 4/200 2026-09-17 16:29 by 0ngsZVSCTS7z
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 5BDX0d0WFp7t 2026-09-15 3/150 2026-09-17 16:27 by 0ngsZVSCTS7z
[基金申请] 国社科系统bug了,是不是要放榜了? +3 kynobel 2026-09-16 3/150 2026-09-17 16:11 by 接好运婷
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 14:27 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 14:27 by ObVzqQQh1rrE
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 Aj1rhIDL5ixY 2026-09-14 6/300 2026-09-17 14:17 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 40ms4Wlo4umh 2026-09-14 6/300 2026-09-17 13:53 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:28 by k06BKNblrGRH
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-17 02:02 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 6/300 2026-09-17 01:50 by Tql5LQhh5rLK
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-16 21:03 by MHf40yA0Uhyb
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 16:00 by i3KBRD19ERkG
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:48 by i3KBRD19ERkG
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:50 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
信息提示
请填处理意见