TALK=T;RUN(1,1)
DISPLAY
The problem considered is an axisymmetric, supersonic, rocket
exhaust plume discharging into a subsonic moving air stream.
The exhaust gas has a specific-heat ratio of 1.21, a Mach number
of 2.39, a total pressure of 5.786MPa, and a total temperature
of 3462 K.
The free-stream has a Mach number of 0.029, a static temperature
of 288K, and a static pressure of 1.01325 bar.
The plume is underexpanded with a nozzle-to-ambient static
pressure ratio of 3.8.
In addition to the continuity and momentum equations,
conservation equations are solved for the total enthalpy equation
and the mass fraction of the exhaust gas.
The turbulence is represented by means of the k-e turbulence mode
The results of this elliptic-solver case may be compared with
those produced for library case 161 by the much more economical
parabolic solver.
ENDDIS
PHOTON USE
P
0.20443E+04 0.15633E+04 CR
CON MACH X 1 FI;.5
msg Mach number contours
msg Press to continue
pau;cl
vec x 1 sh;
pau;cl
con tmp1 x 1 fi;.5
pau;cl
con rho1 x 1 fi;.5
pau;cl
con c1 x 1 fi;.5
pause
ENDUSE
AUTOPLOT USE
FILE
PHI 5
D 1 MACH Y 1;PLOT;REDR;LEVEL Y 1
msg Mach number distribution along flow axis
msg Press to continue
pause
cl
D 1 P1 Y 1;PLOT
shift y
1.01325e5 1
scale;redr
msg static pressure distribution along flow axis
msg Press to continue
pause
cl
D 1 RHO1 Y 1;PLOT;REDR
msg static density distribution along flow axis
msg Press to continue
pause
cl
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(2D Supersonic Rocket Exhaust Plume )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 14
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 1
XULAST =0.1
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 40
YVLAST =0.181625
YFRAC(1)=1.0E-02 ;YFRAC(2)=0.02
YFRAC(3)=0.03 ;YFRAC(4)=0.04
YFRAC(5)=0.05 ;YFRAC(6)=0.06
YFRAC(7)=0.07 ;YFRAC(8)=0.08
YFRAC(9)=0.09 ;YFRAC(10)=0.1
YFRAC(11)=0.11 ;YFRAC(12)=0.12
YFRAC(13)=0.13 ;YFRAC(14)=0.14
YFRAC(15)=0.15 ;YFRAC(16)=0.16
YFRAC(17)=0.17 ;YFRAC(18)=0.18
YFRAC(19)=0.19 ;YFRAC(20)=0.2
YFRAC(21)=0.212068 ;YFRAC(22)=0.231849
YFRAC(23)=0.256185 ;YFRAC(24)=0.284049
YFRAC(25)=0.31487 ;YFRAC(26)=0.348272
YFRAC(27)=0.383986 ;YFRAC(28)=0.421806
YFRAC(29)=0.46157 ;YFRAC(30)=0.503143
YFRAC(31)=0.546416 ;YFRAC(32)=0.591293
YFRAC(33)=0.637692 ;YFRAC(34)=0.685542
YFRAC(35)=0.734781 ;YFRAC(36)=0.78535
YFRAC(37)=0.837201 ;YFRAC(38)=0.890287
YFRAC(39)=0.944566 ;YFRAC(40)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 80
ZWLAST =0.5812
ZFRAC(1)=0.0125 ;ZFRAC(2)=0.025
ZFRAC(3)=0.0375 ;ZFRAC(4)=0.05
ZFRAC(5)=0.0625 ;ZFRAC(6)=0.075
ZFRAC(7)=0.0875 ;ZFRAC(8)=0.1
ZFRAC(9)=0.1125 ;ZFRAC(10)=0.125
ZFRAC(11)=0.1375 ;ZFRAC(12)=0.15
ZFRAC(13)=0.1625 ;ZFRAC(14)=0.175
ZFRAC(15)=0.1875 ;ZFRAC(16)=0.2
ZFRAC(17)=0.2125 ;ZFRAC(18)=0.225
ZFRAC(19)=0.2375 ;ZFRAC(20)=0.25
ZFRAC(21)=0.2625 ;ZFRAC(22)=0.275
ZFRAC(23)=0.2875 ;ZFRAC(24)=0.3
ZFRAC(25)=0.3125 ;ZFRAC(26)=0.325
ZFRAC(27)=0.3375 ;ZFRAC(28)=0.35
ZFRAC(29)=0.3625 ;ZFRAC(30)=0.375
ZFRAC(31)=0.3875 ;ZFRAC(32)=0.4
ZFRAC(33)=0.4125 ;ZFRAC(34)=0.425
ZFRAC(35)=0.4375 ;ZFRAC(36)=0.45
ZFRAC(37)=0.4625 ;ZFRAC(38)=0.475
ZFRAC(39)=0.4875 ;ZFRAC(40)=0.5
ZFRAC(41)=0.5125 ;ZFRAC(42)=0.525
ZFRAC(43)=0.5375 ;ZFRAC(44)=0.55
ZFRAC(45)=0.5625 ;ZFRAC(46)=0.575
ZFRAC(47)=0.5875 ;ZFRAC(48)=0.6
ZFRAC(49)=0.6125 ;ZFRAC(50)=0.625
ZFRAC(51)=0.6375 ;ZFRAC(52)=0.65
ZFRAC(53)=0.6625 ;ZFRAC(54)=0.675
ZFRAC(55)=0.6875 ;ZFRAC(56)=0.7
ZFRAC(57)=0.7125 ;ZFRAC(58)=0.725
ZFRAC(59)=0.7375 ;ZFRAC(60)=0.75
ZFRAC(61)=0.7625 ;ZFRAC(62)=0.775
ZFRAC(63)=0.7875 ;ZFRAC(64)=0.8
ZFRAC(65)=0.8125 ;ZFRAC(66)=0.825
ZFRAC(67)=0.8375 ;ZFRAC(68)=0.85
ZFRAC(69)=0.8625 ;ZFRAC(70)=0.875
ZFRAC(71)=0.8875 ;ZFRAC(72)=0.9
ZFRAC(73)=0.9125 ;ZFRAC(74)=0.925
ZFRAC(75)=0.9375 ;ZFRAC(76)=0.95
ZFRAC(77)=0.9625 ;ZFRAC(78)=0.975
ZFRAC(79)=0.9875 ;ZFRAC(80)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(5)=V1
NAME(7)=W1 ;NAME(12)=KE
NAME(13)=EP ;NAME(14)=H1
NAME(16)=C1 ;NAME(144)=EPKE
NAME(145)=CP1 ;NAME(146)=TMP1
NAME(147)=LEN1 ;NAME(148)=ENUT
NAME(149)=MACH ;NAME(150)=RHO1
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,Y,N,N,N)
SOLUTN(V1,Y,Y,N,N,N,N)
SOLUTN(W1,Y,Y,N,N,N,N)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(H1,Y,Y,N,N,N,Y)
SOLUTN(C1,Y,Y,N,N,N,Y)
SOLUTN(EPKE,Y,N,N,N,N,N)
SOLUTN(CP1,Y,N,N,N,N,Y)
SOLUTN(TMP1,Y,N,N,N,N,Y)
SOLUTN(LEN1,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
SOLUTN(MACH,Y,N,N,N,N,Y)
SOLUTN(RHO1,Y,N,N,N,N,Y)
DEN1 = 150
VIST = 148
TEMP1 = 146
LEN1 = 147
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(P1,Y,Y,Y,N,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
TERMS(H1,N,Y,Y,Y,Y,Y)
TERMS(C1,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWRH1 = T
NEWENT = T
UCONV = T
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
DENPCO = T
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 = GRND5 ;TMP1 = GRND6 ;EL1 = GRND4
TSURR =0. ;TEMP0 =0. ;PRESS0 =1.01325E+05
DVO1DT =0. ;DRH1DP = GRND5
RHO1A =289.096985 ;RHO1B =337.162598 ;RHO1C =0.714286
TMP1A =0. ;TMP1B =0. ;TMP1C =0.
TMP2A =0. ;TMP2B =0. ;TMP2C =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL =1.0E-05 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(V1)=1. ;PRNDTL(W1)=1.
PRNDTL(KE)=1. ;PRNDTL(EP)=1.
PRNDTL(H1)=1. ;PRNDTL(C1)=1.
PRT(V1)=1. ;PRT(W1)=1.
PRT(KE)=1. ;PRT(EP)=1.314
PRT(H1)=1. ;PRT(C1)=1.
CP1 = GRND7 ;CP2 =1.
CP1A =1011.839478 ;CP1B =1942.698486 ;CP1C =0. ;CP1D =0.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=2.83775E+05 ;FIINIT(V1)=0.
FIINIT(W1)=2246.340088 ;FIINIT(KE)=5.046044E+04
FIINIT(EP)=5.126943E+08 ;FIINIT(H1)=6.725079E+06
FIINIT(C1)=1. ;FIINIT(EPKE)=1.0E-10
FIINIT(CP1)=1.0E-10 ;FIINIT(TMP1)=2163.
FIINIT(LEN1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
FIINIT(MACH)=1.0E-10 ;FIINIT(RHO1)=0.528053
PATCH(INITA ,INIVAL, 1, 1, 21, 40, 1, 80, 1, 1)
INIT(INITA ,P1 ,0. ,0. )
INIT(INITA ,W1 ,0. ,10. )
INIT(INITA ,KE ,0. ,1.0E-10 )
INIT(INITA ,EP ,0. ,9.000001E-15 )
INIT(INITA ,H1 ,0. ,2.914598E+05 )
INIT(INITA ,C1 ,0. ,0. )
INIT(INITA ,TMP1,0. ,288. )
INIT(INITA ,RHO1,0. ,1.216972 )
INIADD = F
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
PATCH(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(IN ,LOW , 1, 1, 1, 20, 1, 1, 1, 1)
COVAL(IN ,P1 , FIXFLU ,1186.18689 )
COVAL(IN ,V1 ,0. ,0. )
COVAL(IN ,W1 ,0. ,2246.340088 )
COVAL(IN ,KE ,0. ,5.046044E+04 )
COVAL(IN ,EP ,0. ,5.126943E+08 )
COVAL(IN ,H1 ,0. ,6.725079E+06 )
COVAL(IN ,C1 ,0. ,1. )
PATCH(LB ,LOW , 1, 1, 21, 40, 1, 1, 1, 1)
COVAL(LB ,P1 , FIXFLU ,12.16972 )
COVAL(LB ,W1 ,0. ,10. )
COVAL(LB ,KE ,0. ,1.0E-10 )
COVAL(LB ,EP ,0. ,9.000001E-15 )
COVAL(LB ,H1 ,0. ,2.914598E+05 )
PATCH(NB ,NORTH , 1, 1, 40, 40, 1, 80, 1, 1)
COVAL(NB ,P1 ,1000. ,0. )
COVAL(NB ,W1 ,0. ,10. )
COVAL(NB ,KE ,0. ,1.0E-10 )
COVAL(NB ,EP ,0. ,9.000001E-15 )
COVAL(NB ,H1 ,0. ,2.914598E+05 )
COVAL(NB ,C1 ,0. ,0. )
PATCH(OUT ,HIGH , 1, 1, 1, 40, 80, 80, 1, 1)
COVAL(OUT ,P1 ,1000. ,0. )
COVAL(OUT ,KE ,0. , SAME )
COVAL(OUT ,EP ,0. , SAME )
COVAL(OUT ,H1 ,0. , SAME )
COVAL(OUT ,C1 ,0. , SAME )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 700 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = F
RESREF(P1)=7.825901E-14 ;RESREF(V1)=1.757964E-10
RESREF(W1)=1.757964E-10 ;RESREF(KE)=3.948984E-09
RESREF(EP)=4.012295E-05 ;RESREF(H1)=5.26298E-07
RESREF(C1)=7.825901E-14
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(V1)=10
LITER(W1)=10 ;LITER(KE)=20
LITER(EP)=20 ;LITER(H1)=20
LITER(C1)=20
ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
ENDIT(EP)=1.0E-03 ;ENDIT(H1)=1.0E-03
ENDIT(C1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(V1,FALSDT,4.56098E-06)
RELAX(W1,FALSDT,4.56098E-06)
RELAX(KE,FALSDT,4.56098E-06)
RELAX(EP,FALSDT,4.56098E-06)
RELAX(H1,FALSDT,4.56098E-06)
RELAX(C1,LINRLX,1.)
RELAX(EPKE,LINRLX,1.)
RELAX(CP1,LINRLX,1.)
RELAX(TMP1,LINRLX,1.)
RELAX(LEN1,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
RELAX(MACH,LINRLX,1.)
RELAX(RHO1,LINRLX,1.)
KELIN = 3
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
VARMAX(C1)=1. ;VARMIN(C1)=1.0E-10
VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
VARMAX(CP1)=1.0E+10 ;VARMIN(CP1)=-1.0E+10
VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=200.
VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
VARMAX(MACH)=1.0E+10 ;VARMIN(MACH)=-1.0E+10
VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=1.0E-05
************************************************************
Group 19. Data transmitted to GROUND
NAMGRD =CONV
GENK = T
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(P1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(KE,Y,N,Y,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(H1,Y,N,Y,Y,Y,Y)
OUTPUT(C1,Y,N,Y,Y,Y,Y)
OUTPUT(EPKE,Y,N,Y,N,N,N)
OUTPUT(CP1,Y,N,Y,N,N,N)
OUTPUT(TMP1,Y,N,Y,N,N,N)
OUTPUT(LEN1,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
OUTPUT(MACH,Y,N,Y,N,Y,N)
OUTPUT(RHO1,Y,N,Y,N,Y,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 20 ;IZMON = 40
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 700 ;NUMCLS = 5
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = 20
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP