TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2011, Date 23/03/11
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
Library case Y101:
SUPERSONIC FLOW IN A DUCT: FGEM
This run illustrates what happens in a supersonic flow through
a rectangular duct. The use of GREX3 facilities is illustrated
via the activation of the isentropic gas law.
Interesting variations include changes to: the gas-law constants,
and to the inlet velocity.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 302
************************************************************
Group 1. Run Title
TEXT(SUPERSONIC FLOW IN DUCT : Y101 )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,42,18,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(147)=VABS ;NAME(148)=PTOT
NAME(149)=MACH ;NAME(150)=RHO1
* Solved variables list
SOLVE(P1,U1,V1)
* Stored variables list
STORE(RHO1,MACH,PTOT,VABS)
************************************************************
Group 8. Terms & Devices
TERMS (U1 ,Y,Y,N,Y,Y,Y)
TERMS (V1 ,Y,Y,N,Y,Y,Y)
NEWRH1 = T
************************************************************
Group 9. Properties
PRESS0 =1.
RHO1 = GRND3
RHO1A =1. ;RHO1B =0.714
RHO1C =0.
ENUL =1.0E-05
CP1 =1005.
ENUT =0.
DRH1DP = GRND3
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=0. ;FIINIT(U1)=2.
FIINIT(RHO1)=1.
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
No PATCHes used for this Group
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 200
SARAH =0.05
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 5 ;IYMON = 14 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 1
NYPRIN = 1
IPLTF = 1 ;IPLTL = 40
ISWPRF = 1 ;ISWPRL = 100000
PATCH(DOMAIN, CONTUR, 1, 40, 1, 18, 1, 1, 1, 1)
PLOT(DOMAIN, P1, 0., 10.)
PLOT(DOMAIN, U1, 0., 10.)
PLOT(DOMAIN, V1, 0., 10.)
PLOT(DOMAIN, RHO1, 0., 10.)
PATCH(IYEQNY, PROFIL, 1, 40, 18, 18, 1, 1, 1, 1)
PLOT(IYEQNY, P1, 0., 0.)
PATCH(IYEQ1, PROFIL, 1, 40, 1, 1, 1, 1, 1, 1)
PLOT(IYEQ1, P1, 0., 0.)
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,-0.069532,0.104168,0.992126)
GVIEW(UP,-3.951097E-03,0.994497,-0.104694)
> DOM, SIZE, 5.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, MONIT, 5.040000E-01, 7.524999E-01, 5.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID, RSET_X_1, 5, 1.000000E+00
> GRID, RSET_X_2, 4, 1.000000E+00
> GRID, RSET_X_3, 4, 1.000000E+00
> GRID, RSET_X_4, 2, 1.000000E+00
> GRID, RSET_X_5, 8, 1.000000E+00
> GRID, RSET_X_6, 19, 1.000000E+00
> GRID, RSET_Y_1, 8, 1.000000E+00
> GRID, RSET_Y_2, 2, 1.000000E+00
> GRID, RSET_Y_3, 8, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> OBJ, NAME, X1
> OBJ, POSITION, 1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, wirexyz
> OBJ, TYPE, NULL
> OBJ, NAME, IN
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube3t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VELOCITY, 2. ,0.5 ,0.
> OBJ, NAME, OUTLET
> OBJ, POSITION, 5.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube12t
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, COEFFICIENT, 1000.
> OBJ, NAME, B4
> OBJ, POSITION, 5.600000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.130000E+00, 4.400001E-01, 1.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, fine
> OBJ, TYPE, FINE_GRID_VOL
> OBJ, FINEFACTS, 3, 3, 1
> OBJ, NAME, B5
> OBJ, POSITION, 1.470000E+00, 5.599999E-01, 0.000000E+00
> OBJ, SIZE, 1.130000E+00, 4.400001E-01, 1.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, fine
> OBJ, TYPE, FINE_GRID_VOL
> OBJ, FINEFACTS, 3, 3, 1
STOP