DISPLAY
 
 
   Coal-combustion model; two-phase, with slip.
 
   Reactions: C (s) + 0.5 O2  > CO       (exothermic )
              CO    + 0.5 O2  > CO2      (exothermic )
              C(s)  + CO2     > 2CO      (endothermic)
              C(s)  + H2O     > CO  + H2 (endothermic)
              H2    + 0.5 O2  > H2O      (exothermic )
 
 
  ENDDIS
  Fuel composition and consequences for stoichiometry
  ** CINCL & HINCL are the mass fractions of carbon & hydrogen
     in the coal
REAL(CINCL, HINCL, NINCL)
CINCL=0.86; HINCL=0.05
     In the following formulae:
     0.232 is the mass of oxygen per unit mass of air
     0.768 is the mass of nitrogen per unit mass of air
     2.0, 12.0, 16.0, 18.0, 32.0 & 44.0 are molecular weights of
     H2,  C,    O,    H2O,  O2   & CO2  respectively
 
REAL(FS,FS2)
  ** FS is the mass of fuel per unit mass of air/fuel mixture
     to convert all carbon and oxygen to carbon monoxide.
 
FS=0.232/(0.232 + CINCL*16.0/12.0)
  ** FS2 is the mass of fuel per unit mass of air/fuel mixture
     to convert all carbon, hydrogen  and oxygen to
     carbon dioxide and water vapour.
 
FS2=0.232/(0.232 + CINCL*32.0/12.0 + HINCL*16./2.0)
 
     Thermodynamic data
 
  ** HCCO2 = heat of combustion for  C + O2      --> CO2
     HCCO  =  "    "     "        "  C + 0.5 O2  --> CO
     HHH2O =  "    "     "       "   H2 + 0.5 O2 --> H2O
 
    * the heat of reaction for    C  + O2     -> CO2, HCCO2: 3.279E7
    * the heat of reaction for    C  + 0.5*O2 -> CO , HCCO : 9.208E6
    * the heat of reaction for    H2 + 0.5*O2 -> H2O, HHH2O: 1.209E6
    * the specific heat at constant pressure,         CP   : 1.100E3
    H = CP*T + HCHX*YCHX + HCOCO2*YCO * HH2*YH2
 
REAL(HCCO2,HCCO,HHH2O,HCHX)
 HCCO2=32.792E6; HCCO=9.208E6; HHH2O=120.9E6
 HCHX=CINCL*HCCO2 + HINCL*HHH2O
 HCHX
 
REAL(HGIN,GALF,HF,HA2,HO,HSIN)
  ** take cpsolid=cpgas=1.1e3
REAL(RHOGIN); CP1=1.1E3; CP2=CP1
    Data concerning the inflows of fuel and air
REAL(FLOG,FLOS,VELO,VELG,CHATIM,LENGTH,RGIN,RSIN)
REAL(TGIN,TSIN,BURNRATE)
TGIN=500.; TSIN=350.; FLOS=1.0; VELO=1.;LENGTH=10.0
BURNRATE=10000
HGIN=CP1*TGIN
HSIN=CP2*TSIN + HCHX
HSIN
HGIN
 
PHINT(H1)=7GASES; PHINT(H2)=7GASES
 
FLOG=2.25*FLOS
mesg( Mass fraction of carbon in coal is   :cincl:
mesg( Mass fraction of hydrogen in coal is :hincl:
mesg( Mass-flow ratio air/coal is          :FLOG/FLOS:
mesg( Ratio for combustion to CO & H2 is   :1.0/FS:
mesg( Ratio for combustion to CO2 & H2O is :1.0/FS2:
 
mesg( Input data OK? (Y/n)
ans=y
readvdu(ans,char,y)
real(fac)
if(:ans:.eq.n) then
 mesg(Increase air flow by factor (default=1)
 readvdu(fac,real,1.)
 flog=flog*fac
 mesg(Increase combustion rate by factor (default=1)
 readvdu(fac,real,1.)
 burnrate=burnrate*fac
endif