DISPLAY
 
  This run restarts from the data produced by case 481, and
  performs two further steps, namely:-
 
  1. by calling GXCREK, it computes the equilibrium values of
     O, OH, and H corresponding to the combustion reaction; and
 
  2. by solving the source and transport equations for N and NO
     in accordance with the Zeldovich mechanism, it computes the
     concentrations of these substances.
 
  The coding is to be found in file GXNOX.F
  press RETURN to continue
readvdu(NPHI,INT,50)
  ENDDIS
  noxdbegin
  ns 7
  nr 3
  noxdend
load(481)
  noxdbegin
  ns 7
  nr 3
  noxdend
  PRT0BEGIN
  PRINT0 0.0
  PRT0END
TEXT(NOX Post-Processing Of Case 481   
TITLE
    NOXCAL = T  activates solution for nox
BOOLEAN(NOXCAL);NOXCAL=T
 
STORE(P1,GAS,FUE,SHAD,H1,H2,MIXF,C2)
SOLUTN(5,Y,N,N,N,N,N);SOLUTN(5,Y,N,N,N,N,N)
SOLVE(C3,C5);STORE(C4,C6);NAME(C3)=XN;NAME(C5)=XNO
STORE(XO,XO2,XH,XOH,CRKT,NOSR)
CINT(XNO)=0.0;CINT(XN)=0.0;CINT(C4)=0.0;CINT(C6)=0.0
 
RESTRT(P1,GVEL,CVEL,GAS,FUE,SHAD,H1,H2,MIXF,C2,TMP1,TMP2)
RESTRT(RHO1,RMIX,YH2O,YH2,YN2,YCO2,YO2,YCO,CFIP,MDOT)
 
  ** NO & N reaction-rate source terms
PATCH(NOXSR,FREEVL,1,NX,1,NY,1,NZ,1,1)
COVAL(NOXSR,XNO,FIXFLU,GRND);COVAL(NOXSR,XN,GRND,GRND)
 
PATCH(FUELNOX,SOUTH,1,NX,1,1,1,NZ,1,1)
COVAL(FUELNOX,XNO,FLOG,0.0);COVAL(FUELNOX,XN,FLOG,0.0)
 
PATCH(OUTLET,NORTH,1,NX,NY,NY,1,NZ,1,1);COVAL(OUTLET,P1,OUTCO1,0.0)
 
INLETS=SKIP;INLETG=SKIP
TSTSWP=10;LSWEEP=30
RELAX(XNO,FALSDT,CHATIM); RELAX(XN,FALSDT,CHATIM)
VARMIN(XNO)=0.; VARMAX(XNO)=1.; VARMIN(XN)=0.; VARMAX(XN)=1.
 
LSG62=NOXCAL
 
OUTPUT(C4,N,N,N,N,N,N); OUTPUT(C6,N,N,N,N,N,N)
OUTPUT(H1,N,N,N,N,N,N); OUTPUT(H2,N,N,N,N,N,N)
OUTPUT(7,N,N,N,N,N,N); OUTPUT(11,N,N,N,N,N,N)
OUTPUT(FUE,N,N,N,N,N,N); OUTPUT(TMP2,N,N,N,N,N,N)
OUTPUT(CFIP,N,N,N,N,N,N); OUTPUT(MDOT,N,N,N,N,N,N)
 
SAVE=F