************************************************************
Q1 created by VDI menu, Version 3.4, Date 04/09/01
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
PHOTON USE
p
use patgeo
rot z ang -10
con tfav x 1 fil;.001
con tfav z 10 fil;.001
con tfav x 5 fil;.001
vec x m; vec z m; vec y m
msg( Fluid-average temperature
pause;con cl;vec cl;red
con tmp1 x 1 fil;.001
con tmp1 z 10 fil;.001
con tmp1 x 5 fil;.001
msg( Single-fluid temperature
pause;con cl;red
con crav x 1 fil;.001
con crav z 10 fil;.001
con crav x 5 fil;.001
msg( Fluid-average radiation flux
pause;con cl;red
con crad x 1 fil;.001
con crad z 10 fil;.001
con crad x 5 fil;.001
msg( Single-fluid radiation flux
pause;con cl;red
con cr1 x 1 fil;.001
con cr1 z 10 fil;.001
con cr1 x 5 fil;.001
msg( 1st-fluid radiation
pause;con cl;red
con cr2 x 1 fil;.001
con cr2 z 10 fil;.001
con cr2 x 5 fil;.001
msg( 2nd-fluid radiation
pause;con cl;red
con cr3 x 1 fil;.001
con cr3 z 10 fil;.001
con cr3 x 5 fil;.001
msg( 3rd-fluid radiation
pause;con cl;red
con cr4 x 1 fil;.001
con cr4 z 10 fil;.001
con cr4 x 5 fil;.001
msg( 4th-fluid radiation
pause;con cl;red
con cr5 x 1 fil;.001
con cr5 z 10 fil;.001
con cr5 x 5 fil;.001
msg( 5th-fluid radiation
pause;con cl;red
con f1 x 1 fil;.001
con f1 z 10 fil;.001
con f1 x 5 fil;.001
msg( 1st-fluid 'presence probability'
pause;con cl;red
con f2 x 1 fil;.001
con f2 z 10 fil;.001
con f2 x 5 fil;.001
msg( 2nd-fluid 'presence probability'
pause;con cl;red
con f3 x 1 fil;.001
con f3 z 10 fil;.001
con f3 x 5 fil;.001
msg( 3th-fluid 'presence probability'
pause;con cl;red
con f4 x 1 fil;.001
con f4 z 10 fil;.001
con f4 x 5 fil;.001
msg( 4th-fluid 'presence probability'
pause;con cl;red
con f5 x 1 fil;.001
con f5 z 10 fil;.001
con f5 x 5 fil;.001
msg( 5th-fluid 'presence probability'
ENDUSE
DISPLAY
MFM: 5 fluid population
Reaction: CO + 0.5 (O2+0.86N2) -> CO2 + 0.43N2
i.e. FUEL + 0.5 OXID -> PROD
Stoichiometric ratio: S=1.
Rate: Rfu = 4*rho*EP/KE*min(FUEL,OXID)
Radiation: P-1 model
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 302
************************************************************
Group 1. Run Title
TEXT(P-1 radiation model for 5-fluid MFM )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,5,5,20)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME( 70) =RFU5 ; NAME( 71) =RFU4
NAME( 72) =RFU3 ; NAME( 73) =RFU2
NAME( 74) =RFU1 ; NAME( 75) =SUMA
NAME( 76) =SUM5 ; NAME( 77) =SUM4
NAME( 78) =SUM3 ; NAME( 79) =SUM2
NAME( 80) =SUM1 ; NAME( 81) =SUMS
NAME( 82) =PR5 ; NAME( 83) =PR4
NAME( 84) =PR3 ; NAME( 85) =PR2
NAME( 86) =PR1 ; NAME( 87) =PRAV
NAME( 88) =VPR4 ; NAME( 89) =VPR3
NAME( 90) =VPR2 ; NAME( 91) =OXAV
NAME( 92) =VOX4 ; NAME( 93) =VOX3
NAME( 94) =VOX2 ; NAME( 95) =OX5
NAME( 96) =OX4 ; NAME( 97) =OX3
NAME( 98) =OX2 ; NAME( 99) =OX1
NAME(100) =VFU5 ; NAME(101) =VFU4
NAME(102) =VFU3 ; NAME(103) =VFU2
NAME(104) =FUAV ; NAME(105) =FU5
NAME(106) =FU4 ; NAME(107) =FU3
NAME(108) =FU2 ; NAME(109) =FU1
NAME(110) =VHF4 ; NAME(111) =VHF3
NAME(112) =VHF2 ; NAME(113) =COF4
NAME(114) =COF3 ; NAME(115) =COF2
NAME(116) =TF5 ; NAME(117) =TF4
NAME(118) =TF3 ; NAME(119) =TF2
NAME(120) =TF1 ; NAME(121) =TFAV
NAME(122) =HFAV ; NAME(123) =HF5
NAME(124) =HF4 ; NAME(125) =HF3
NAME(126) =HF2 ; NAME(127) =HF1
NAME(128) =CRAV ; NAME(129) =CR5
NAME(130) =CR4 ; NAME(131) =CR3
NAME(132) =CR2 ; NAME(133) =CR1
NAME(134) =FAV ; NAME(135) =F5
NAME(136) =F4 ; NAME(137) =F3
NAME(138) =F2 ; NAME(139) =F1
NAME(140) =ENUT ; NAME(141) =EL1
NAME(142) =EPKE ; NAME(143) =DEN1
NAME(144) =FSUM ; NAME(145) =FUEL
NAME(146) =TMP1 ; NAME(147) =PROD
NAME(148) =OXID ; NAME(149) =CRAD
NAME(150) =MIXF
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,W1 ,H1 ,PR5 ,PR4 ,PR3 )
SOLVE(PR2 ,PR1 ,OX5 ,OX4 ,OX3 ,OX2 ,OX1 ,FU5 )
SOLVE(FU4 ,FU3 ,FU2 ,FU1 ,HF5 ,HF4 ,HF3 ,HF2 )
SOLVE(HF1 ,CR5 ,CR4 ,CR3 ,CR2 ,CR1 ,F5 ,F4 )
SOLVE(F3 ,F2 ,F1 ,FUEL,CRAD,MIXF)
* Stored variables list
STORE(OXID,PROD,TMP1,FSUM,DEN1,EPKE,EL1 ,ENUT)
STORE(FAV ,CRAV,HFAV,TFAV,TF1 ,TF2 ,TF3 ,TF4 )
STORE(TF5 ,COF2,COF3,COF4,VHF2,VHF3,VHF4,FUAV)
STORE(VFU2,VFU3,VFU4,VFU5,VOX2,VOX3,VOX4,OXAV)
STORE(VPR2,VPR3,VPR4,PRAV,SUMS,SUM1,SUM2,SUM3)
STORE(SUM4,SUM5,SUMA,RFU1,RFU2,RFU3,RFU4,RFU5)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
TURMOD(KEMODL)
************************************************************
Group 8. Terms & Devices
TERMS (H1 ,N,Y,Y,Y,Y,Y)
TERMS (CR5 ,N,N,Y,Y,N,Y)
TERMS (CR4 ,N,N,Y,Y,Y,Y)
TERMS (CR3 ,N,N,Y,Y,N,Y)
TERMS (CR2 ,N,N,Y,Y,Y,Y)
TERMS (CR1 ,N,N,Y,Y,N,Y)
TERMS (CRAD,N,N,Y,Y,Y,N)
NEWRH1 = T
************************************************************
Group 9. Properties
RHO1 = GRND6
RHO1A = 2.800000E+01 ;RHO1B = 5.600000E+01
RHO1C = 5.600000E+01
PRESS0 = 1.000000E+05
TMP1 = GRND8
TMP1A = 1.500000E+03 ;TMP1B = 1.500000E+03
TMP1C = 1.500000E+03
TMP2A = 5.000000E-01 ;TMP2B = 7.450000E+06
TMP2C = 0.000000E+00
CP1 = 1.500000E+03
ENUL = 1.000000E-05
DVO1DT = 3.330000E-03
PRNDTL(CR5 ) = GRND ;PRNDTL(CR4 ) = GRND
PRNDTL(CR3 ) = GRND ;PRNDTL(CR2 ) = GRND
PRNDTL(CR1 ) = GRND ;PRNDTL(CRAD) = GRND
PRT (EP ) = 1.314000E+00 ;PRT (CR5 ) = 1.000000E+09
PRT (CR4 ) = 1.000000E+09 ;PRT (CR3 ) = 1.000000E+09
PRT (CR2 ) = 1.000000E+09 ;PRT (CR1 ) = 1.000000E+09
PRT (CRAD) = 1.000000E+09
EMISS = 4.000000E-01 ;SCATT = 1.000000E-01
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Echo PLANT settings
PLANTBEGIN
* Exchange coefficients for incident radiations
LAMPR(CRAD)=:ENUL:*DEN1*3.*(:EMISS:+:SCATT:)
LAMPR(CR1) =:ENUL:*DEN1*3.*(:EMISS:+:SCATT:)
LAMPR(CR2) =:ENUL:*DEN1*3.*(:EMISS:+:SCATT:)
LAMPR(CR3) =:ENUL:*DEN1*3.*(:EMISS:+:SCATT:)
LAMPR(CR4) =:ENUL:*DEN1*3.*(:EMISS:+:SCATT:)
LAMPR(CR5) =:ENUL:*DEN1*3.*(:EMISS:+:SCATT:)
* Single-fluid sources due to radiation
PATCH(RADSOR,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
VAL=:EMISS:*(CRAD-4.*5.68e-8*TMP1**4)
COVAL(RADSOR,H1,FIXFLU,GRND)
VAL=4.*5.68e-8*TEMP1**4
COVAL(RADSOR,CRAD,EMISS,GRND)
* Within-fluid sources due to radiation
PATCH(RADFLUID,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
** 1st fluid
VAL=:EMISS:*(CR1-4.*5.68e-8*TF1**4)*F1
COVAL(RADFLUID,HF1,FIXFLU,GRND)
VAL=4.*5.68e-8*TF1**4*F1
COVAL(RADFLUID,CR1,EMISS,GRND)
** 2nd fluid
VAL=:EMISS:*(CR2-4.*5.68e-8*TF2**4)*F2
COVAL(RADFLUID,HF2,FIXFLU,GRND)
VAL=4.*5.68e-8*TF2**4*F2
COVAL(RADFLUID,CR2,EMISS,GRND)
** 3rd fluid
VAL=:EMISS:*(CR3-4.*5.68e-8*TF3**4)*F3
COVAL(RADFLUID,HF3,FIXFLU,GRND)
VAL=4.*5.68e-8*TF3**4*F3
COVAL(RADFLUID,CR3,EMISS,GRND)
** 4th fluid
VAL=:EMISS:*(CR4-4.*5.68e-8*TF4**4)*F4
COVAL(RADFLUID,HF4,FIXFLU,GRND)
VAL=4.*5.68e-8*TF4**4*F4
COVAL(RADFLUID,CR4,EMISS,GRND)
** 5th fluid
VAL=:EMISS:*(CR5-4.*5.68e-8*TF5**4)*F5
COVAL(RADFLUID,HF5,FIXFLU,GRND)
VAL=4.*5.68e-8*TF5**4*F5
COVAL(RADFLUID,CR5,EMISS,GRND)
* Within-fluid fuel sources
PATCH(FUELFLUI,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
VAL=-Rfu1*F1
COVAL(FUELFLUI,FU1,FIXFLU,GRND)
VAL=-Rfu2*F2
COVAL(FUELFLUI,FU2,FIXFLU,GRND)
VAL=-Rfu3*F3
COVAL(FUELFLUI,FU3,FIXFLU,GRND)
VAL=-Rfu4*F4
COVAL(FUELFLUI,FU4,FIXFLU,GRND)
VAL=-Rfu5*F5
COVAL(FUELFLUI,FU5,FIXFLU,GRND)
* Within-fluid oxidant sources
PATCH(OXIDFLUI,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
VAL=-Rfu1*F1
COVAL(OXIDFLUI,OX1,FIXFLU,GRND)
VAL=-Rfu2*F2
COVAL(OXIDFLUI,OX2,FIXFLU,GRND)
VAL=-Rfu3*F3
COVAL(OXIDFLUI,OX3,FIXFLU,GRND)
VAL=-Rfu4*F4
COVAL(OXIDFLUI,OX4,FIXFLU,GRND)
VAL=-Rfu5*F5
COVAL(OXIDFLUI,OX5,FIXFLU,GRND)
* Within-fluid product sources
PATCH(PRODFLUI,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
VAL=2.*Rfu1*F1
COVAL(PRODFLUI,PR1,FIXFLU,GRND)
VAL=2.*Rfu2*F2
COVAL(PRODFLUI,PR2,FIXFLU,GRND)
VAL=2.*Rfu3*F3
COVAL(PRODFLUI,PR3,FIXFLU,GRND)
VAL=2.*Rfu4*F4
COVAL(PRODFLUI,PR4,FIXFLU,GRND)
VAL=2.*Rfu5*F5
COVAL(PRODFLUI,PR5,FIXFLU,GRND)
* Micromixing contributions to CVA transport*
PATCH(MICRCVA,PHASEM,1,NX,1,NY,1,NZ,1,1)
** Within-fluid enthalpy
CO = 5.*EPKE*COF2
VAL= VHF2
COVAL(MICRCVA,HF2,GRND,GRND)
CO = 5.*EPKE*COF3
VAL= VHF3
COVAL(MICRCVA,HF3,GRND,GRND)
CO = 5.*EPKE*COF4
VAL= VHF4
COVAL(MICRCVA,HF4,GRND,GRND)
** Within-fluid fuel
CO = 5.*EPKE*COF2
VAL= VFU2
COVAL(MICRCVA,FU2,GRND,GRND)
CO = 5.*EPKE*COF3
VAL= VFU3
COVAL(MICRCVA,FU3,GRND,GRND)
CO = 5.*EPKE*COF4
VAL= VFU4
COVAL(MICRCVA,FU4,GRND,GRND)
** Within-fluid oxidant
CO = 5.*EPKE*COF2
VAL= VOX2
COVAL(MICRCVA,OX2,GRND,GRND)
CO = 5.*EPKE*COF3
VAL= VOX3
COVAL(MICRCVA,OX3,GRND,GRND)
CO = 5.*EPKE*COF4
VAL= VOX4
COVAL(MICRCVA,OX4,GRND,GRND)
** Within-fluid products
CO = 5.*EPKE*COF2
VAL= VPR2
COVAL(MICRCVA,PR2,GRND,GRND)
CO = 5.*EPKE*COF3
VAL= VPR3
COVAL(MICRCVA,PR3,GRND,GRND)
CO = 5.*EPKE*COF4
VAL= VPR4
COVAL(MICRCVA,PR4,GRND,GRND)
* COefficients for child fluids
COF2=F1*F3+F1*F4/2.+F1*F5/3.+tiny
COF3=F2*F4+F1*F4/2.+F2*F5/2.+F1*F5/3.+tiny
COF4=F3*F5+F2*F5/2.+F1*F5/3.+tiny
* VALues for child fluids
** Enthalpy
VHF2=HF1/2.+(F1*F3/2.*HF3+F1*F4/4.*HF4+F1*F5/6.*HF5)/COF2
VHF3=(F1*F4/4.+F1*F5/6.)*HF1+(F2*F4/2.+F2*F5/4.)*HF2
VHF3=VHF3+(F1*F4/4.+F2*F4/2.)*HF4+(F1*F5/6.+F2*F5/4.)*HF5
VHF3=VHF3/COF3
VHF4=HF5/2.+(F3*F5/2.*HF3+F2*F5/4.*HF2+F1*F5/6.*HF1)/COF4
** Fuel
VFU2=FU1/2.+(F1*F3/2.*FU3+F1*F4/4.*FU4+F1*F5/6.*FU5)/COF2
VFU3=(F1*F4/4.+F1*F5/6.)*FU1+(F2*F4/2.+F2*F5/4.)*FU2
VFU3=VFU3+(F1*F4/4.+F2*F4/2.)*FU4+(F1*F5/6.+F2*F5/4.)*FU5
VFU3=VFU3/COF3
VFU4=FU5/2.+(F3*F5/2.*FU3+F2*F5/4.*FU2+F1*F5/6.*FU1)/COF4
** Oxidant
VOX2=OX1/2.+(F1*F3/2.*OX3+F1*F4/4.*OX4+F1*F5/6.*OX5)/COF2
VOX3=(F1*F4/4.+F1*F5/6.)*OX1+(F2*F4/2.+F2*F5/4.)*OX2
VOX3=VOX3+(F1*F4/4.+F2*F4/2.)*OX4+(F1*F5/6.+F2*F5/4.)*OX5
VOX3=VOX3/COF3
VOX4=OX5/2.+(F3*F5/2.*OX3+F2*F5/4.*OX2+F1*F5/6.*OX1)/COF4
** Products
VPR2=PR1/2.+(F1*F3/2.*PR3+F1*F4/4.*PR4+F1*F5/6.*PR5)/COF2
VPR3=(F1*F4/4.+F1*F5/6.)*PR1+(F2*F4/2.+F2*F5/4.)*PR2
VPR3=VPR3+(F1*F4/4.+F2*F4/2.)*PR4+(F1*F5/6.+F2*F5/4.)*PR5
VPR3=VPR3/COF3
VPR4=PR5/2.+(F3*F5/2.*PR3+F2*F5/4.*PR2+F1*F5/6.*PR1)/COF4
** Within-fluid temperatures
* TMP2B is the neat of combustion
* TMP1A is specific heat
TF1=(HF1-:TMP2B:*FU1)/:TMP1A:
TF2=(HF2-:TMP2B:*FU2)/:TMP1A:
TF3=(HF3-:TMP2B:*FU3)/:TMP1A:
TF4=(HF4-:TMP2B:*FU4)/:TMP1A:
TF5=(HF5-:TMP2B:*FU5)/:TMP1A:
** Within-fluid reaction rates
Rfu1=4.*EPKE*AMIN1(FU1,OX1)
Rfu2=4.*EPKE*AMIN1(FU2,OX2)
Rfu3=4.*EPKE*AMIN1(FU3,OX3)
Rfu4=4.*EPKE*AMIN1(FU4,OX4)
Rfu5=4.*EPKE*AMIN1(FU5,OX5)
** Population averaged properties
* Mixture fraction
FAV =F1*0.0+F2*0.25+F3*0.5+F4*0.75+F5*1.
* Radiation
CRAV=F1*CR1+F2*CR2+F3*CR3+F4*CR4+F5*CR5
* Enthalpy
HFAV=F1*HF1+F2*HF2+F3*HF3+F4*HF4+F5*HF5
* Temperature
TFAV=F1*TF1+F2*TF2+F3*TF3+F4*TF4+F5*TF5
* Fuel mass fraction
FUAV=F1*FU1+F2*FU2+F3*FU3+F4*FU4+F5*FU5
* Oxygen mass fraction
OXAV=F1*OX1+F2*OX2+F3*OX3+F4*OX4+F5*OX5
* Product mass fraction
PRAV=F1*PR1+F2*PR2+F3*PR3+F4*PR4+F5*PR5
** Checks for continuity
* Sum of fluid mass fractions
FSUM=F1+F2+F3+F4+F5
* Sum of single-fluid components
SUMS=OXID+FUEL+PROD
* Sum of 1st-fluid components
SUM1=OX1+FU1+PR1
* Sum of 2nd-fluid components
SUM2=OX2+FU2+PR2
* Sum of 3rd-fluid components
SUM3=OX3+FU3+PR3
* Sum of 4th-fluid components
SUM4=OX4+FU4+PR4
* Sum of 5th-fluid components
SUM5=OX5+FU5+PR5
* Sum of fluid-averged components
SUMA=OXAV+FUAV+PRAV
PLANTEND
************************************************************
Group 11.Initialise Var/Porosity Fields
No PATCHes used for this Group
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
PATCH (MICROMIX,PHASEM,1,5,1,5,1,20,1,1)
COVAL (MICROMIX,F5 , GRND , SAME )
COVAL (MICROMIX,F4 , GRND , SAME )
COVAL (MICROMIX,F3 , GRND , SAME )
COVAL (MICROMIX,F2 , GRND , SAME )
COVAL (MICROMIX,F1 , GRND , SAME )
PATCH (CHSO ,PHASEM,0,0,0,0,0,0,1,1)
COVAL (CHSO ,FUEL, GRND9 , GRND9 )
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 2000
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 3.000000E-01)
RELAX(U1 ,FALSDT, 1.000000E-01)
RELAX(V1 ,FALSDT, 1.000000E-01)
RELAX(W1 ,FALSDT, 1.000000E-01)
RELAX(KE ,FALSDT, 1.000000E-01)
RELAX(EP ,FALSDT, 1.000000E-01)
RELAX(H1 ,FALSDT, 1.000000E-01)
RELAX(PR5 ,FALSDT, 1.000000E-01)
RELAX(PR4 ,FALSDT, 1.000000E-01)
RELAX(PR3 ,FALSDT, 1.000000E-01)
RELAX(PR2 ,FALSDT, 1.000000E-01)
RELAX(PR1 ,FALSDT, 1.000000E-01)
RELAX(OX5 ,FALSDT, 1.000000E-01)
RELAX(OX4 ,FALSDT, 1.000000E-01)
RELAX(OX3 ,FALSDT, 1.000000E-01)
RELAX(OX2 ,FALSDT, 1.000000E-01)
RELAX(OX1 ,FALSDT, 1.000000E-01)
RELAX(FU5 ,FALSDT, 1.000000E-01)
RELAX(FU4 ,FALSDT, 1.000000E-01)
RELAX(FU3 ,FALSDT, 1.000000E-01)
RELAX(FU2 ,FALSDT, 1.000000E-01)
RELAX(FU1 ,FALSDT, 1.000000E-01)
RELAX(HF5 ,FALSDT, 1.000000E-01)
RELAX(HF4 ,FALSDT, 1.000000E-01)
RELAX(HF3 ,FALSDT, 1.000000E-01)
RELAX(HF2 ,FALSDT, 1.000000E-01)
RELAX(HF1 ,FALSDT, 1.000000E-01)
RELAX(CR5 ,FALSDT, 1.000000E-01)
RELAX(CR4 ,FALSDT, 1.000000E-01)
RELAX(CR3 ,FALSDT, 1.000000E-01)
RELAX(CR2 ,FALSDT, 1.000000E-01)
RELAX(CR1 ,FALSDT, 1.000000E-01)
RELAX(F5 ,FALSDT, 1.000000E-01)
RELAX(F4 ,FALSDT, 1.000000E-01)
RELAX(F3 ,FALSDT, 1.000000E-01)
RELAX(F2 ,FALSDT, 1.000000E-01)
RELAX(F1 ,FALSDT, 1.000000E-01)
RELAX(DEN1,LINRLX, 1.000000E-01)
RELAX(FUEL,FALSDT, 1.000000E-01)
RELAX(CRAD,FALSDT, 1.000000E+01)
RELAX(MIXF,FALSDT, 1.000000E-01)
KELIN = 3
************************************************************
Group 18. Limits
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06
VARMAX(H1 ) = 1.000000E+10 ;VARMIN(H1 ) = 3.000000E+05
VARMAX(PR5 ) = 1.000000E+00 ;VARMIN(PR5 ) = 0.000000E+00
VARMAX(PR4 ) = 1.000000E+00 ;VARMIN(PR4 ) = 0.000000E+00
VARMAX(PR3 ) = 1.000000E+00 ;VARMIN(PR3 ) = 0.000000E+00
VARMAX(PR2 ) = 1.000000E+00 ;VARMIN(PR2 ) = 0.000000E+00
VARMAX(PR1 ) = 1.000000E+00 ;VARMIN(PR1 ) = 0.000000E+00
VARMAX(OX5 ) = 1.000000E+00 ;VARMIN(OX5 ) = 0.000000E+00
VARMAX(OX4 ) = 1.000000E+00 ;VARMIN(OX4 ) = 0.000000E+00
VARMAX(OX3 ) = 1.000000E+00 ;VARMIN(OX3 ) = 0.000000E+00
VARMAX(OX2 ) = 1.000000E+00 ;VARMIN(OX2 ) = 0.000000E+00
VARMAX(OX1 ) = 1.000000E+00 ;VARMIN(OX1 ) = 0.000000E+00
VARMAX(FU4 ) = 1.000000E+00 ;VARMIN(FU4 ) = 0.000000E+00
VARMAX(FU3 ) = 1.000000E+00 ;VARMIN(FU3 ) = 0.000000E+00
VARMAX(FU2 ) = 1.000000E+00 ;VARMIN(FU2 ) = 0.000000E+00
VARMAX(FU1 ) = 1.000000E+00 ;VARMIN(FU1 ) = 0.000000E+00
VARMAX(HF5 ) = 1.000000E+10 ;VARMIN(HF5 ) = 3.000000E+05
VARMAX(HF4 ) = 1.000000E+10 ;VARMIN(HF4 ) = 3.000000E+05
VARMAX(HF3 ) = 1.000000E+10 ;VARMIN(HF3 ) = 3.000000E+05
VARMAX(HF2 ) = 1.000000E+10 ;VARMIN(HF2 ) = 3.000000E+05
VARMAX(HF1 ) = 1.000000E+10 ;VARMIN(HF1 ) = 3.000000E+05
VARMAX(F5 ) = 1.000000E+00 ;VARMIN(F5 ) = 0.000000E+00
VARMAX(F4 ) = 1.000000E+00 ;VARMIN(F4 ) = 0.000000E+00
VARMAX(F3 ) = 1.000000E+00 ;VARMIN(F3 ) = 0.000000E+00
VARMAX(F2 ) = 1.000000E+00 ;VARMIN(F2 ) = 0.000000E+00
VARMAX(F1 ) = 1.000000E+00 ;VARMIN(F1 ) = 0.000000E+00
VARMAX(DEN1) = 2.000000E+00 ;VARMIN(DEN1) = 1.000000E-03
VARMAX(FUEL) = 1.000000E+00 ;VARMIN(FUEL) = 0.000000E+00
VARMAX(TMP1) = 1.000000E+04 ;VARMIN(TMP1) = 3.000000E+02
VARMAX(PROD) = 1.000000E+00 ;VARMIN(PROD) = 0.000000E+00
VARMAX(OXID) = 1.000000E+00 ;VARMIN(OXID) = 0.000000E+00
VARMAX(MIXF) = 1.000000E+00 ;VARMIN(MIXF) = 0.000000E+00
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
NAMSAT =MOSG
GENK = T
ASAP = T
CHSOA = 5.000000E-01 ;CHSOB = 4.000000E+00
SPEDAT(SET,MFM,MFMMOD,C,MFM)
SPEDAT(SET,MFM,NFLUIDS,I,5)
SPEDAT(SET,MFM,NFLR,I,1)
SPEDAT(SET,MFM,NFLF,I,5)
SPEDAT(SET,MFM,CONREA,R,0.00000E+00)
SPEDAT(SET,MFM,CONMIX,R,5.00000E+00)
SPEDAT(SET,MFM,VISCON,R,0.00000E+00)
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
OUTPUT(P1 ,Y,Y,Y,Y,Y,Y)
OUTPUT(U1 ,Y,Y,Y,N,Y,Y)
OUTPUT(V1 ,Y,Y,Y,N,Y,Y)
OUTPUT(W1 ,Y,Y,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 3 ;IYMON = 5 ;IZMON = 20
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NZPRIN = 1
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
NOWIPE = T
GVIEW(P,-5.953568E-01,6.788709E-01,-4.297495E-01)
GVIEW(UP,6.290849E-01,7.265853E-01,2.762717E-01)
> DOM, SIZE, 1.000000E+00, 1.000000E+00, 1.000000E+01
> DOM, MONIT, 5.000000E-01, 9.000000E-01, 9.750000E+00
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_Y_1, 5, 1.000000E+00
> DOM, RELAX, 3.000000E-01
> OBJ, NAME, FUELIN
> OBJ, POSITION, 6.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 4.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, 0.000000E+00, 2.000000E+00, 0.000000E+00
> OBJ, TEMPERATURE, 3.000000E+02
> OBJ, INLET_FU5, 1.000000E+00
> OBJ, INLET_FU4, 1.000000E+00
> OBJ, INLET_FU3, 1.000000E+00
> OBJ, INLET_FU2, 1.000000E+00
> OBJ, INLET_FU1, 1.000000E+00
> OBJ, INLET_HF5, 7.900000E+06
> OBJ, INLET_HF4, 7.900000E+06
> OBJ, INLET_HF3, 7.900000E+06
> OBJ, INLET_HF2, 7.900000E+06
> OBJ, INLET_HF1, 7.900000E+06
> OBJ, INLET_F5, 1.000000E+00
> OBJ, INLET_FUEL, 1.000000E+00
> OBJ, INLET_MIXF, 1.000000E+00
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, OUTLET
> OBJ, POSITION, 0.000000E+00, 1.000000E+00, 9.000000E+00
> OBJ, SIZE, 6.000000E-01, 0.000000E+00, 1.000000E+00
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, -1.026000E+04
> OBJ, COEFFICIENT, 1.000000E+03
> OBJ, TURBULENCE, -1.026000E+04,-1.026000E+04
> OBJ, NAME, AIRIN1
> OBJ, POSITION, 6.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, SIZE, 4.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube10
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, 0.000000E+00, 2.000000E+00, 0.000000E+00
> OBJ, TEMPERATURE, 3.000000E+02
> OBJ, INLET_OX5, 1.000000E+00
> OBJ, INLET_OX4, 1.000000E+00
> OBJ, INLET_OX3, 1.000000E+00
> OBJ, INLET_OX2, 1.000000E+00
> OBJ, INLET_OX1, 1.000000E+00
> OBJ, INLET_HF5, 4.500000E+05
> OBJ, INLET_HF4, 4.500000E+05
> OBJ, INLET_HF3, 4.500000E+05
> OBJ, INLET_HF2, 4.500000E+05
> OBJ, INLET_HF1, 4.500000E+05
> OBJ, INLET_F1, 1.000000E+00
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, AIRIN2
> OBJ, POSITION, 4.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube10
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, 0.000000E+00, 2.000000E+00, 0.000000E+00
> OBJ, TEMPERATURE, 3.000000E+02
> OBJ, INLET_OX5, 1.000000E+00
> OBJ, INLET_OX4, 1.000000E+00
> OBJ, INLET_OX3, 1.000000E+00
> OBJ, INLET_OX2, 1.000000E+00
> OBJ, INLET_OX1, 1.000000E+00
> OBJ, INLET_HF5, 4.500000E+05
> OBJ, INLET_HF4, 4.500000E+05
> OBJ, INLET_HF3, 4.500000E+05
> OBJ, INLET_HF2, 4.500000E+05
> OBJ, INLET_HF1, 4.500000E+05
> OBJ, INLET_F1, 1.000000E+00
> OBJ, TURB-INTENS, 5.000000E+00
STOP