1 2 3 4 5 6 7 8 9 10

CHAPTER 5: PROBLEMS WITH PRACTICAL PROJECT APPLICATIONS

Shortcuts:

TITLE Problem 5.01 - AXISYMMETRIC SWIRLER FLOW >>> C60A

TITLE Problem 5.02 - AXISYMMETRIC SWIRLER FLOW >>> C60B

TITLE CASE 5.03: Ramjet Combustion: Default 4 Step 8 Species Reaction

TITLE CASE 5.04: Ramjet Combustion: Default 2 Step 8 Species Reaction

TITLE CASE 5.05: Ramjet Combustion: User Specified 2 Step 6 Species

-Back To Top

************************************************************************
TITLE Problem 5.02 - AXISYMMETRIC SWIRLER FLOW >>> C60B
************************************************************************
/
GRID NODEs 62 by 18 by 7
COOR X Y Z 'ANSD502.GRD'
/----------------------------------------------------------------------/
SELEct (1,13,1) (21,999,999)
BLOCKAGE at SELEcted zone
/----------------------------------------------------------------------/
/ Initial and Boundary Conditions
/
WALL at undefined outer boundaries
SELEct (1,1,1) (1,12,999)
INLE X- at SELEcted subdomain
OUTL X+
SYMM Y-
SET U 37.0
SET V SINE THETA -44.0 1.0 0.0 0.0
SET W COSI THETA 44.0 1.0 0.0 0.0
SET K 24.3215
SET L 2.54E-3
PERIodic in Z THETA 5.729578
/----------------------------------------------------------------------/
/ Fluid Properties and Constants
DENSity 3.48
VISC 2.679E-5
/----------------------------------------------------------------------/
/ SOLUTION OPTIONS
/DIFF SECO SKEW U V W P
/
DIAGnostic output U V W P K E at (30,10,4) every 50 steps
DEBUG GEOMERTY OFF
FLUX DEFAult output OFF
SELEct from (1,1,4) to (99,99,4) interval (6,2,1)
OUTPut U V W P K L in XY planes SELEcted zone
SAVE OFF U V W P K E L MU on file '502.SAV'
/----------------------------------------------------------------------/
/ OPERATIONAL CONTROL
// ::::::> akr_note: Only partial convergence with input below
CONV GLOB 1.E-06 1 !!! Change to 1.E-10 for full run
SOLVE STEADY 200 min 200 !!! Change to 2000 steps for full run
/
END
/

-Back To Top

***********************************************************************
TITLE CASE 5.03: Ramjet Combustion: Default 4 Step 8 Species Reaction
***********************************************************************
/
GRID NODEs 32 by 22
/
COORDINATE X and R from 'ANSD504.GRD' CORNers
/
////////////////////////////////////////////////////////////////////////
* GEOMETRY SPECIFICATION
////////////////////////////////////////////////////////////////////////
/
WALL at undefined outer boundaries
SYMMetry at Y- boundary
/
LOCATE (01,01) to (01,09) with ID=INLEt
INLET X- at SELEcted boundary
OUTLET X+
/
/
////////////////////////////////////////////////////////////////////////
* Initial and Boundary Conditions
////////////////////////////////////////////////////////////////////////
/
REACTION
/
SET U = 190. at ID=INLEt at X- boundary only
SET T = 600 at ID=INLEt at X- boundary only
SET FU = 0.05 at ID=INLEt at X- boundary only
SET value of K = 1000 everywhere
SET value of L = 0.06 everywhere
LOCAte (02,02) to (15,11) FIELd nodes only
SET T 2000 SELEcted subdomian !!! Start combustion
/
ADIABATIC WALLS
LIMIT T to min = 300, max = 2500
/
////////////////////////////////////////////////////////////////////////
* OUTPUT CONTROL
////////////////////////////////////////////////////////////////////////
/
/
DIAGNOSTIC node (08,08) every 50 steps
DEBUG GEOMERTY OFF
FLUX DEFault output OFF
SELEct from (1,1) to (999,999) interval (3,2)
OUTPut U,V,P,T,K,L,FF,FU,H2O SELEcted zone
/
////////////////////////////////////////////////////////////////////////
* SOLUTION OPTIONS
////////////////////////////////////////////////////////////////////////
/
RELAx U = 0.3 V=0.3 W=0.3 RHO=0.3 P = 0.7 FF=1.0
CONVERGENCE Max mass source = 1.0E-8
CONVERGENCE PP 1.0E-9 max of 10 Iterations on P-Prime Eqn
/
SOLVE 1000 STEPS in STEADY MODE
SAVE OFF ON '504.SAV'
/
END
/

-Back To Top

***********************************************************************
TITLE CASE 5.04: Ramjet Combustion: Default 2 Step 8 Species Reaction
***********************************************************************
/
GRID NODEs 32 by 22
/
COORDINATE X and R from 'ANSD504.GRD' CORNers
/
////////////////////////////////////////////////////////////////////////
* GEOMETRY SPECIFICATION
////////////////////////////////////////////////////////////////////////
/
WALL at undefined outer boundaries
SYMMetry at Y- boundary
/
LOCATE (01,01) to (01,09) with ID=INLEt
INLET X- at SELEcted boundary
OUTLET X+
/
SET U = 190. at ID=INLEt at X- boundary only
SET T = 600 at ID=INLEt at X- boundary only
SET FU = 0.05 at ID=INLEt at X- boundary only
SET value of K = 1000 everywhere
SET value of L = 0.06 everywhere
/
////////////////////////////////////////////////////////////////////////
* Initial and Boundary Conditions
////////////////////////////////////////////////////////////////////////
/
REACTION TWO STEP
/
LOCAte (02,02) to (15,11) FIELd nodes only
SET T 2000 SELEcted subdomain !!! Start combustion
/
ADIABATIC WALLS
LIMIT T to min = 300, max = 2500
/
////////////////////////////////////////////////////////////////////////
* OUTPUT CONTROL
////////////////////////////////////////////////////////////////////////
/
DIAGNOSTIC node (08,08) every 50 steps
DEBUG GEOMERTY OFF
FLUX DEFault output OFF
SELEct from (1,1) to (999,999) interval (3,2)
OUTPut U,V,P,T,K,L,FF,FU,H2O SELEcted zone
/
////////////////////////////////////////////////////////////////////////
* SOLUTION OPTIONS
////////////////////////////////////////////////////////////////////////
/
RELAx U = 0.3 V=0.3 W=0.3 RHO=0.3 P = 0.7 FF=1.0
CONVERGENCE Max mass source = 1.0E-8
CONVERGENCE PP 1.0E-9 max of 10 Iterations on P-Prime Eqn
/
SOLVE 1000 STEPS in STEADY MODE
SAVE OFF ON '505.SAV'
/
END
/
/

-Back To Top

***********************************************************************
TITLE CASE 5.05: Ramjet Combustion: User Specified 2 Step 6 Species
***********************************************************************
/
GRID NODEs 32 by 22
/
COORDINATE X and R from 'ANSD504.GRD' CORNers
/
////////////////////////////////////////////////////////////////////////
* GEOMETRY SPECIFICATION
////////////////////////////////////////////////////////////////////////
/
WALL at undefined outer boundaries
SYMMetry at Y- boundary
/
LOCATE (01,01) to (01,09) with ID=INLEt
INLET X- at SELEcted boundary
OUTLET X+
/
SET U = 190. at ID=INLEt at X- boundary only
SET T = 600 at ID=INLEt at X- boundary only
SET FU = 0.05 at ID=INLEt at X- boundary only
SET FF = 0.05 at ID=INLEt at X- boundary only
SET value of K = 1000 everywhere
SET value of L = 0.06 everywhere
/
////////////////////////////////////////////////////////////////////////
* Initial and Boundary Conditions
////////////////////////////////////////////////////////////////////////
/
/ Set up combustion from scratch
REACTION is for a HYDROCARBON fuel
/
LOCAte (02,02) to (15,11) FIELd nodes only
SET T 2000 SELEcted subdomian !!! Start combustion
/
ADIABATIC WALLS
LIMIT T to min = 300, max = 2500
/
SPECIFIC heat Gordon-McBride relations: 'cphtDFLT.acr'
/
/ Define molecular weights of the species (FU = C12H22.5)
GAS FU =139.26293 CO = 28.01055 H2O= 18.01534
CO2= 44.00995 O2 = 31.99880 N2 = 28.01340
/
/ Allocate space for reaction rates, called R1, R2
ALLOCATE R1
ALLOCATE R2
/
/ Define reaction rates
/
/ FU + (10/2 + 19/4) O2 => 10 CO + 19/2 H2O
REACtion R1 FU O2 3.3E+14 2.7E+04 0.5 1.0 0.0 MASS TEMP EBU
/
/ CO + 0.5 O2 => CO2
REACtion R2 CO O2 6.0E+08 1.25E+04 1.0 1.0 0.0 MASS TEMP EBU
/
/ EBU coefficients for each reaction
/
EBU 3.0 10 !! start control after 10 steps
EBU R1 FU 3.0
EBU R2 CO 3.0 !! Assume only CO control; default option
/EBU R2 CO O2 4.0 1.75 !! 2 * MW_CO /MW_O2; CO & O2 control
/
/ Source terms for FU and CO transport equations
/
/ FU source term
SOURce REACtion type FU LINEar -1.0 * R1
/
/ CO source term ( 2.018 = 12*CO/FU ) <<< C12H22.5
/ CO source term ( 2.01134 = 10*CO/FU ) <<< C10H19
SOURce REACtion type CO LINEar +2.01134 * R1 -1.0 * R2
/
/ Algebraic equations for total H-atom (FH), total O-atom (FO),
/ CO2, H2O, O2, and N2 by atom balance. FF is total C-atom.
/
SET ALWAys LINEar FH = 0.0 +1.0 * FF
SET ALWAys LINEar FO = 0.2331 -0.2331 * FF
SET ALWAys LINEar CO2 = 3.1602 * FF -3.1602 * FU -1.5712 * CO + 0.0
SET ALWAys LINEar H2O = 1.2289 * FH -1.2289 * FU +0.0
SET ALWAys LINEar O2 = 1.0 * FO -0.5714 * CO -0.7273 * CO2
-0.8889 * H2O + 0.0
SET ALWAys LINEar N2 = -1.0 * FU -1.0 * CO -1.0 * CO2 -1.0 * H2O
-1.0 * O2 + 1.0
/
////////////////////////////////////////////////////////////////////////
* OUTPUT CONTROL
////////////////////////////////////////////////////////////////////////
/
/
DIAGNOSTIC node (08,08) every 50 steps
DEBUG GEOMERTY OFF
FLUX DEFault output OFF
SELEct from (1,1) to (999,999) interval (3,2)
OUTPut U,V,P,T,K,L,FF,FU,H2O SELEcted zone
/
////////////////////////////////////////////////////////////////////////
* SOLUTION OPTIONS
////////////////////////////////////////////////////////////////////////
/
RELAx U = 0.3 V=0.3 W=0.3 RHO=0.3 P = 0.7 FF=1.0
CONVERGENCE Max mass source = 1.0E-8
CONVERGENCE PP 1.0E-9 max of 10 Iterations on P-Prime Eqn
/
SOLVE 1000 STEPS in STEADY MODE
SAVE OFF ON '506.SAV'
/
END
/

-Back To Top