FLUENT-UDF代码-颗粒(煤或生物质)床层多相流燃烧过程数值模拟-【第2部分-完结】.doc

PAGE1

Fluent-UDF_Coal-Combustion-Multiphase-Flow-Processes

Codesforvariouskindsofheterogeneousreactions

【第二部分-Part2】

DEFINE_HET_RXN_RATE(SteamGasif_Rev_MGAS,c,t,hr,mw,yi,rr,rr_t)

{

Thread**pt=THREAD_SUB_THREADS(t);

Thread*tp=pt[0];/*gasphase*/

intindex_phase=Get_Phase_Index(hr);

Thread*ts=pt[index_phase];/*solidphase*/

*rr=0;

doubledirection=0.0,mol_weight,y_carbon,rr_turb=1e+20;

/*

CO+H2H2O+1/25Soot

Thereversesteamgasificationreaction,whichisCO+H21/25Soot+H2Oiswrittenasthereaction

shownabove.SoanegativeratemeansCOandH2isconsumedandH2OandSootaregenerated.Notethatno

C(s)isgenerated,thestoichcoeffforC(s)intheabovereationiszero.

Setthephaseandspeciesindices.Ashspeciesindexisinitializedtozero,withallotherindices.

AshspeciesindexisusedasaflagtoexecuteSetSpeciesIndexonlyonce.Thisisdonebythefirst

reaction,definedintheheterogeneousreactionpanelinFLUENTGUI.

*/

if(IS_ASH==0)

SetSpeciesIndex();

if(MGAS_Gasif)

{

doubleRoRT=C_R(c,tp)*UNIVERSAL_GAS_CONSTANT*C_T(c,tp);

doublep_h2o=RoRT*yi[IP_H2O][IS_H2O]/mw[IP_H2O][IS_H2O]/101325.;

doublep_co=RoRT*yi[IP_CO][IS_CO]/mw[IP_CO][IS_CO]/101325.;

doublep_h2=RoRT*yi[IP_H2][IS_H2]/mw[IP_H2][IS_H2]/101325.;

y_carbon=yi[IP_SOOT][IS_SOOT];

mol_weight=mw[IP_SOOT][IS_SOOT];

if(rp_ke)

rr_turb=Turbulent_rr(c,t,hr,yi);

if(C_VOF(c,ts)=eps_s_small)

{

*rr=rr_steam_gasif(c,t,ts,tp,p_h2o,p_co,p_h2,y_carbon,mol_weight,direction);/*mol/(cm^3.s)*/

if(direction0.0)/*positivevalueimpliesC(s)+H2OCO+H2*/

*rr=0.0;

else/*negativevalueimpliesCO+H2H2O+1/25Soot*/

{