TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2019, Date 25/01/19
CPVNAM=VDI; SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
SUPERSONIC FLOW OVER DIAMOND-SHAPE BODY
The supersonic flow over a 2D diamond-shaped body is solved
here on Cartesian grid by "cut-off" technique of PARSOL.
The case is aimed to demonstrate the PARSOL performance for
compressible flow along with the relevant treatments and settings.
The demonstration case considers the supersonic air flow over a
diamond-shaped profile placed in a rectangular channel.
In this case the pressure distribution and accompanying velocity
field have to be calculated along with auxiliary field of density
to get the Mach number contours.
Symmetry of the flow allows to consider the half of the actual
flow pattern.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(SUPERSONIC FLOW OVER DIAMOND-SHAPE BODY )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,60,30,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(146)=PRPS ;NAME(147)=VABS
NAME(148)=MACH ;NAME(149)=PTOT
NAME(150) =RHO1
* Solved variables list
SOLVE(P1,U1,V1)
* Stored variables list
STORE(RHO1,PTOT,MACH,VABS,PRPS)
************************************************************
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(PRPS)=-1. ;FIINIT(RHO1)=1.
No PATCHes used for this Group
INIADD = T
************************************************************
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 = 500
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,0.3 )
RELAX(U1 ,FALSDT,5. )
RELAX(V1 ,FALSDT,5. )
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
ISG62 = 0
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=0.4267 ;EX(U1)=3.659
EX(V1)=0.3438 ;EX(PRPS)=0.9767
EX(VABS)=3.698 ;EX(MACH)=3.119
EX(PTOT)=52.59 ;EX(RHO1)=1.08
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 4 ;IYMON = 21 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 4
NYPRIN = 2
IPLTF = 1 ;IPLTL = 40
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.,0.,1.)
GVIEW(UP,0.,1.,0.)
GVIEW(DEPTH,2.73014E+04)
GVIEW(VDIS,1.457694)
GVIEW(CENTRE,1.307784,0.214941,0.5)
> DOM, SIZE, 5.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, MONIT, 2.955557E-01, 6.833333E-01, 5.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID, RSET_X_1, 9, 1.000000E+00
> GRID, RSET_X_2, 7, 1.000000E+00
> GRID, RSET_X_3, 7, 1.000000E+00
> GRID, RSET_X_4, 37, 1.000000E+00
> GRID, RSET_Y_1, 30, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> 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, GEOMETRY, cube3t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.
> OBJ, VELOCITY, 4. ,0. ,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, GEOMETRY, cube12t
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.
> OBJ, COEFFICIENT, 1000.
> OBJ, NAME, WEDGE1
> OBJ, POSITION, 7.600000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 5.733000E-01, 2.400000E-01, 1.000000E+00
> OBJ, GEOMETRY, wedge2
> OBJ, ROTATION24, 18
> OBJ, GRID, Y,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 199,Solid allowing fluid-slip at walls
> OBJ, NAME, WEDGE2
> OBJ, POSITION, 1.333300E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 5.700000E-01, 2.400000E-01, 1.000000E+00
> OBJ, GEOMETRY, wedge2
> OBJ, ROTATION24, 9
> OBJ, GRID, Y,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 199,Solid allowing fluid-slip at walls
STOP