| 查看: 2059 | 回复: 4 | |||
[交流]
UDF加载问题，加载在哪，如何设置 已有2人参与
|
|
这个UDF在complied完之后，怎么加载进去，大神如何设置，求指导 #include "udf.h" 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; /* 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 += dh_dt / ( mp * Cp ); dzdt->energy -= dh_dt; { Material *sp; mixture_species_loop(gas_mix,sp,ns) { molwt = MATERIAL_PROP(sp,PROP_mwi); /* molecular weight of gas species */ molwt_bulk += C_YI(c0,t0,ns) / molwt; /* 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 / molwt_bulk / solver_par.molWeight; /* vaporization rate */ real vap_rate = k * molwt * Ap * ( cvap_surf - cvap_bulk ); /* only condensation below vaporization temperature */ if( 0. < vap_rate && Tp < vap_temp ) vap_rate = 0.; dydt -= vap_rate; dzdt->species += vap_rate; /* dT/dt = dh/dt / (m Cp)*/ dydt -= hvap * vap_rate / ( mp * Cp ); /* gas enthalpy source term */ dzdt->energy += hgas * vap_rate; } } } |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然评审意见
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有13人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
-
系统查不到
已经有11人回复
2楼2016-10-22 20:24:27
zhouxman
金虫 (小有名气)
- 应助: 51 (初中生)
- 金币: 1363.4
- 散金: 32
- 红花: 5
- 帖子: 299
- 在线: 115.5小时
- 虫号: 1419015
- 注册: 2011-09-27
- 专业: 零件成形制造
3楼2016-10-23 18:25:25
4楼2016-10-24 15:52:19
5楼2016-10-24 15:53:15
