公式一:
Aahaρcyc(x,t)e−EaRT(x,t)ybO2(x,t)h+Aohoρfe−EoRT(x,t)ybO2(x,t)m(−ya(x,t)−yc(x,t)+1)f+Aphpρfe−EpRT(x,t)(−ya(x,t)−yc(x,t)+1)g+∂T(x,t)∂t((Cc−Cf)ya(x,t)+(Cc−Cf)yc(x,t)+Cf)((ρc−ρf)ya(x,t)+(ρc−ρf)yc(x,t)+ρf)=Akg(Tg(x,t)−T(x,t))⎛⎝1110710−−−√3(dfMpatmug(x,t)μgRT(x,t))3/5+2⎞⎠df+∂T(x,t)∂x(16dfσ∂T(x,t)∂xT(x,t)2+(kc−kf)(∂ya(x,t)∂x+∂yc(x,t)∂x))+∂2T(x,t)∂x2(163dfσT(x,t)3+kc(ya(x,t)+yc(x,t))−kf(ya(x,t)+yc(x,t)−1))(1)