| 查看: 1090 | 回复: 0 | |||
llongtan木虫 (小有名气)
老木虫
|
[求助]
Define_DPM_Heat_Mass的UDF中 solver_par.molWeight[gas_index]的意义
|
|
这两天看UDF的英文帮助中，看到DEFINE_DPM_HEAT_MASS的一个例子，对于solver_par.molWeight的含义不太清楚，请哪位大侠能够解释下。 谢谢了。 DEFINE_DPM_HEAT_MASS(multivap,p,Cp,hgas,hvap,cvap_surf,dydt,dzdt) { int ns; int nc = TP_N_COMPONENTS( p ); /* number of particle components */ cell_t c0 = RP_CELL(&(p->cCell)); /* cell and thread */ Thread *t0 = RP_THREAD( &(p->cCell) ); /* where the particle is in */ Material *gas_mix = THREAD_MATERIAL( t0 ); /* gas mixture material */ Material *cond_mix = p->injection->material;/* particle mixture material */ cphase_state_t *c = &(p->cphase); /* cell info of particle location */ real molwt[MAX_SPE_EQNS]; /* molecular weight of gas species */ real Tp = P_T(p); /* particle temperature */ real mp = P_MASS(p); /* particle mass */ real molwt_bulk = 0.; /* average molecular weight in bulk gas */ real Dp = DPM_DIAM_FROM_VOL( mp / P_RHO(p) ); /* particle diameter */ real Ap = DPM_AREA(Dp); /* particle surface */ real Pr = c->sHeat * c->mu / c->tCond; /* Prandtl number */ real Nu = 2.0 + 0.6 * sqrt( p->Re ) * pow( Pr, 1./3. ); /* Nusselt number */ real h = Nu * c->tCond / Dp; /* Heat transfer coefficient */ real dh_dt = h * ( c->temp - Tp ) * Ap; /* heat source term */ dydt[0] += dh_dt / ( mp * Cp ); dzdt->energy -= dh_dt; { Material *sp; mixture_species_loop(gas_mix,sp,ns) { molwt[ns] = MATERIAL_PROP(sp,PROP_mwi); /* molecular weight of gas species */ molwt_bulk += C_YI(c0,t0,ns) / molwt[ns]; /* average molecular weight */ } } /* prevent division by zero */ molwt_bulk = MAX(molwt_bulk,DPM_SMALL); for( ns = 0; ns < nc; ns++ ) { /* gas species index of vaporization */ int gas_index = TP_COMPONENT_INDEX_I(p,ns); if( gas_index >= 0 ) { /* condensed material */ Material * cond_c = MIXTURE_COMPONENT(cond_mix, ns ); /* vaporization temperature */ real vap_temp = MATERIAL_PROP(cond_c,PROP_vap_temp); /* diffusion coefficient */ real D = MATERIAL_PROP_POLYNOMIAL( cond_c, PROP_binary_diffusivity, c->temp); /* Schmidt number */ real Sc = c->mu / ( c->rho * D ); /* mass transfer coefficient */ real k = ( 2. + 0.6 * sqrt(p->Re) * pow( Sc, 1./3. ) ) * D / Dp; /* bulk gas concentration */ real cvap_bulk = c->pressure / UNIVERSAL_GAS_CONSTANT / c->temp * c->yi[gas_index] / molwt_bulk / solver_par.molWeight[gas_index]; /* vaporization rate */ real vap_rate = k * molwt[gas_index] * Ap * ( cvap_surf[ns] - cvap_bulk ); /* only condensation below vaporization temperature */ if( 0. < vap_rate && Tp < vap_temp ) vap_rate = 0.; dydt[1+ns] -= vap_rate; dzdt->species[gas_index] += vap_rate; /* dT/dt = dh/dt / (m Cp)*/ dydt[0] -= hvap[gas_index] * vap_rate / ( mp * Cp ); /* gas enthalpy source term */ dzdt->energy += hgas[gas_index] * vap_rate; } } } [ Last edited by llongtan on 2012-11-12 at 08:49 ] |
回复此楼» 猜你喜欢
filecode=后面第一个是大写字母
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有81人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
感觉是下周放榜了
已经有7人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
25