| 查看: 1119 | 回复: 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 ] |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
200