[交流] fluent udf 模拟出错 已有3人参与

#include "udf.h" #include "dpm.h" #include "surf.h" #include "math.h" DEFINE_DPM_HEAT_MASS(multiabs,p,Cp,hgas,hvap,cvap_surf,dydt,dzdt) { cell_t c=P_CELL(p); /* cell and thread*/这句有问题么 Thread *t=P_CELL_THREAD(p);/* where the particle is in*/这句有问题么 real Tp=P_T(p);/* particle temperature*/ real Tg=C_T(c,t); real mp=P_MASS(p);/* particle mass*/这句有问题么 real Dp=P_DIAM(p);/* particle diameter*/ real Ap=3.14*Dp*Dp;/*particle surface*/ real x_vel_g,x_vel_p,y_vel_g,y_vel_p,abs_v,slip_x,slip_y,reyp,rho_g,rho_p,mu_g,YI,pg,vol_c,D,Sc,Sh,kg,p_so2,pe_so2,Ns,mp_dot; x_vel_g=C_U(c,t)+1.e-8; y_vel_g=C_V(c,t)+1.e-8; x_vel_p=P_VEL(p)[0]; y_vel_p=P_VEL(p)[1]; slip_x=x_vel_g-x_vel_p; slip_y=y_vel_g-y_vel_p; rho_g=C_R(c,t)+1.e-8; rho_p=2.716; mu_g=C_MU_T(c,t)+1.e-8; YI=C_YI(c,t,0); pg=C_P(c,t)+1.e-8; vol_c=C_VOLUME(c,t)+1.e-8; /*compute slip*/ abs_v=sqrt(slip_x*slip_x+slip_y*slip_y); /*compute Reynold's number*/ reyp=rho_g*abs_v*Dp/mu_g; /* diffusion coefficient*/ D=1.e-7*pow(Tp,1.75)*pow(1.0/29+1.0/64,0.5)/(pg*pow(pow(41.8,1./3)+pow(20.1,1./3),2.)+1.e-8); /* Schmidt number*/ Sc=mu_g/(rho_g*D+1.e-8); Sh=2.+0.55*pow(reyp,0.5)*pow(Sc,1./3); kg=Sh*D/(Dp*8.314*Tp+1.e-8);/*solve kg*/ p_so2=rho_g/(rho_p+1.e-8)*YI*pg; pe_so2=103.9951*exp(-7.5778*Tg); Ns=kg/(Tg*8.314+1.e-8)*(p_so2-pe_so2); dydt[1]-=mp_dot;这里的1带表什么 dzdt->species[0]+=mp_dot;这里的0代表什么 } 气相组分一为so2，第二组分为空气 离散项为h2o 最后想加一句确定浓度的公式，公式为    dC/dt=-Ns*（单位体积内气液传质面积），dC为液相中，可以实现么 11.png

回复此楼