| 查看: 1017 | 回复: 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 ] |
回复此楼» 猜你喜欢
310求调剂
已经有14人回复
-
290调剂生物0860
已经有28人回复
-
山东双非院校考核超级无底线，领导幸灾乐祸，教师遭殃恐
已经有4人回复
-
279学硕食品专业求调剂院校
已经有14人回复
-
0831生医工第一轮调剂失败求助
已经有16人回复
-
电子信息270求调剂
已经有18人回复
-
食品与营养（0955）271求调剂
已经有5人回复
-
293求调剂
已经有11人回复
-
322求调剂
已经有13人回复
-
272分材料子求调剂
已经有50人回复
5