TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2007, Date 06/02/07
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
PHOTON USE
p
phi
20 1
gr ou x 1
msg fuel concentration; vertical (radial) size enlarged ten-fold
con fuel x 1 fi;0.001
pause
msg oxidant concentration; flow is from right to left
con oxid x 1 fi;0.001
pause
msg product concentration
con prod x 1 fi;0.001
pause
msg turbulence energy
con ke x 1 fi;0.001
pause
msg temperature
con tmp1 x 1 fi;0.001
pause
msg density
con rho1 x 1 fi;0.001
pause; con off; red
msg velocity vectors
vec x 1 sh
ENDUSE
DISPLAY
The problem considered is the two-dimensional, axisymmetric,
turbulent combustion of methane and air in a cylindrical
furnace. The fuel and air enter separately, and the EBU
(i.e. eddy-break-up) combustion model is used.
This involves the solution of conservation equations for the
fuel mass fraction, the mixture fraction, the turbulence energy,
ke, and the disipation rate, ep.
The enthalpy equation is also solved, and the temperature is
deduced from it and the concentrations.
An eddy-break-up constant ( CHSOB) of 0.5 is specified.
Other constants are:
stoichiometric mixture fraction (TMP2A) = 5.482456E-02
heat of combustion (TMP2B) = 4.900000E+07
The results show that 87% of the fuel is consumed within the
furnace.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 129
************************************************************
Group 1. Run Title
TEXT(2-D methane-air combustion, EBU model )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,24,35,1.000000E-05)
* Cylindrical-polar grid
CARTES=F
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(143) =EPKE ; NAME(144) =ENUT
NAME(145) =FUEL ; NAME(146) =PROD
NAME(147) =OXID ; NAME(148) =RHO1
NAME(149) =TMP1 ; NAME(150) =MIXF
* Solved variables list
SOLVE(P1 ,V1 ,W1 ,H1 ,FUEL,MIXF)
* Stored variables list
STORE(TMP1,RHO1,OXID,PROD,ENUT,EPKE)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
TURMOD(KEMODL)
************************************************************
Group 8. Terms & Devices
TERMS (H1 ,Y,Y,Y,N,Y,N)
NEWRH1 = T
DENPCO = T
************************************************************
Group 9. Properties
PRESS0 = 1.000000E+05
RHO1 = GRND6
RHO1A = 1.604303E+01 ;RHO1B = 2.894200E+01
RHO1C = 2.918912E+01
ENUL = 1.000000E-05
CP1 = GRND10
CP1A = 1.500000E+03 ;CP1B = 1.500000E+03
CP1C = 1.500000E+03 ;CP1D = 0.000000E+00
TMP1 = GRND8
TMP1A = 1.500000E+03 ;TMP1B = 1.500000E+03
TMP1C = 1.500000E+03
TMP2A = 5.482456E-02 ;TMP2B = 4.900000E+07
TMP2C = 0.000000E+00
PRT (EP ) = 1.314000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(W1 ) = 6.000000E-01 ;FIINIT(KE ) = 9.000000E-04
FIINIT(EP ) = 8.872200E-03 ;FIINIT(H1 ) = 1.215000E+06
FIINIT(ENUT) = 1.001000E-10 ;FIINIT(FUEL) = 1.001000E-10
FIINIT(PROD) = 1.001000E-10 ;FIINIT(OXID) = 1.001000E-10
FIINIT(RHO1) = 1.001000E-10 ;FIINIT(TMP1) = 1.001000E-10
FIINIT(MIXF) = 1.001000E-10
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 (CHSO ,PHASEM,0,0,0,0,0,0,1,1)
COVAL (CHSO ,FUEL, GRND9 , GRND9 )
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1000
RESREF(P1 ) = 3.998963E-17 ;RESREF(V1 ) = 2.250403E-16
RESREF(W1 ) = 2.250403E-16 ;RESREF(KE ) = 3.599066E-20
RESREF(EP ) = 3.547960E-19 ;RESREF(H1 ) = 1.451732E-10
RESREF(MIXF) = 2.003087E-18
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER (P1 ) = 40 ;LITER (H1 ) = 40
LITER (FUEL) = 1 ;LITER (PROD) = 1
LITER (OXID) = 1 ;LITER (MIXF) = 1
ENDIT (H1 ) = 1.000000E-06 ;ENDIT (FUEL) = 1.000000E-06
ENDIT (MIXF) = 1.000000E-06
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
RELAX(V1 ,FALSDT, 3.000000E-03)
RELAX(W1 ,FALSDT, 3.000000E-03)
RELAX(KE ,LINRLX, 3.000000E-01)
RELAX(EP ,LINRLX, 3.000000E-01)
RELAX(H1 ,FALSDT, 1.000000E-03)
RELAX(FUEL,FALSDT, 1.000000E-03)
RELAX(RHO1,LINRLX, 1.000000E-02)
RELAX(MIXF,LINRLX, 3.000000E-01)
KELIN = 3
************************************************************
Group 18. Limits
VARMAX(P1 ) = 1.000000E+10 ;VARMIN(P1 ) =-1.000000E+04
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06
VARMAX(FUEL) = 1.000000E+00 ;VARMIN(FUEL) = 0.000000E+00
VARMAX(PROD) = 1.000000E+00 ;VARMIN(PROD) = 0.000000E+00
VARMAX(OXID) = 1.000000E+00 ;VARMIN(OXID) = 0.000000E+00
VARMAX(RHO1) = 5.000000E+01 ;VARMIN(RHO1) = 1.000000E-03
VARMAX(TMP1) = 4.000000E+03 ;VARMIN(TMP1) = 2.890000E+01
VARMAX(MIXF) = 1.000000E+00 ;VARMIN(MIXF) = 0.000000E+00
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
GENK = T
ASAP = T
PARSOL = T
ISG52 = 2
CHSOA = 5.482456E-02 ;CHSOB = 5.000000E-01
************************************************************
Group 20. Preliminary Printout
ECHO = T
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL = 1.000000E-02
EX (P1 ) = 1.983000E+02 ;EX (V1 ) = 6.237000E-02
EX (W1 ) = 1.968000E+01 ;EX (KE ) = 2.238000E+00
EX (EP ) = 1.306000E+03 ;EX (H1 ) = 6.346000E+06
EX (EPKE) = 3.381000E+02 ;EX (ENUT) = 6.504000E-04
EX (FUEL) = 1.021000E-01 ;EX (PROD) = 3.332000E-01
EX (OXID) = 5.648000E-01 ;EX (RHO1) = 4.834000E-01
EX (TMP1) = 8.967000E+02 ;EX (MIXF) = 1.203000E-01
************************************************************
Group 21. Print-out of Variables
OUTPUT(PROD,Y,N,Y,N,Y,Y)
OUTPUT(OXID,Y,N,Y,N,Y,Y)
OUTPUT(RHO1,Y,N,Y,N,Y,Y)
OUTPUT(TMP1,Y,N,Y,N,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 9 ;IZMON = 32
NPRMON = 100
NPRMNT = 100
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NYPRIN = 1
NZPRIN = 1
YZPR = T
IPLTF = 1 ;IPLTL = 10000
NPLT = 25
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,-9.457793E-01,3.009828E-01,1.221105E-01)
GVIEW(UP,3.054609E-01,9.520084E-01,1.933038E-02)
> DOM, SIZE, 1.000000E-01, 1.000000E-02, 6.000000E-01
> DOM, MONIT, 5.000000E-02, 4.206370E-03, 4.851390E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E-01
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_X_1, 1, 1.000000E+00
> GRID, RSET_Y_1, 12,-1.500000E+00
> GRID, RSET_Y_2, 12, 1.500000E+00
> GRID, RSET_Z_1, 35, 1.050000E+00,G
> OBJ, NAME, INLETF
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E-01, 5.000000E-03, 0.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, polcu5t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 3.000000E+05
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 6.000000E-01
> OBJ, TEMPERATURE, 3.000000E+02
> OBJ, INLET_FUEL, 1.000000E+00
> OBJ, INLET_MIXF, 1.000000E+00
> OBJ, KE_IN, 9.000000E-04
> OBJ, EP_IN, 8.872200E-03
> OBJ, NAME, INLETA
> OBJ, POSITION, 0.000000E+00, 5.000000E-03, 0.000000E+00
> OBJ, SIZE, 1.000000E-01, 5.000000E-03, 0.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, polcu5t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 5.000000E+04
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 5.892573E+00
> OBJ, TEMPERATURE, 3.000000E+02
> OBJ, KE_IN, 8.680604E-02
> OBJ, EP_IN, 8.404124E+00
> OBJ, NAME, OUTLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 6.000000E-01
> OBJ, SIZE, 1.000000E-01, 1.000000E-02, 0.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, polcubet
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+03
> OBJ, TURBULENCE, SAME , SAME
> OBJ, NAME, WFNN
> OBJ, POSITION, 0.000000E+00, 1.000000E-02, 0.000000E+00
> OBJ, SIZE, 1.000000E-01, 0.000000E+00, 6.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, polcu10
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
STOP