************************************************************
Q1 created by VDI menu, Version 3.1, Date 05/11/97
CPVNAM=VDI;SPPNAM=Core
************************************************************
DISPLAY
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
************************************************************
Group 1. Run Title
TEXT(Supersonic Flow In A Duct
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,42,18,1)
* Set overall domain extent:
* xulast yvlast zwlast
name
XSI= 5.000000E+00; YSI= 1.000000E+00; ZSI= 1.000000E+00
RSET(D,CHAM )
* Set objects: x0 y0 z0
* dx dy dz
name
XPO= 1.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00
XSI= 0.000000E+00; YSI= 1.000000E+00; ZSI= 1.000000E+00
RSET(B,X1 )
XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00
XSI= 0.000000E+00; YSI= 1.000000E+00; ZSI= 1.000000E+00
RSET(B,IN )
XPO= 5.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00
XSI= 0.000000E+00; YSI= 1.000000E+00; ZSI= 1.000000E+00
RSET(B,OUTLET )
XPO= 5.600000E-01; YPO= 0.000000E+00; ZPO= 0.000000E+00
XSI= 1.130000E+00; YSI= 4.400001E-01; ZSI= 1.000000E+00
RSET(B,B4 )
XPO= 1.470000E+00; YPO= 5.599999E-01; ZPO= 0.000000E+00
XSI= 1.130000E+00; YSI= 4.400001E-01; ZSI= 1.000000E+00
RSET(B,B5 )
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* 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
RHO1 = GRND3
RHO1A = 1.000E+00 ;RHO1B = 7.140E-01 ;RHO1C = 0.000E+00
PRESS0 = 1.000E+00
CP1 = 1.005E+03
ENUL = 1.000E-05 ;ENUT = 0.000E+00
DRH1DP = GRND3
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1 ) = 0.000E+00 ;FIINIT(U1 ) = 2.000E+00
FIINIT(RHO1) = 1.000E+00
No PATCHes used for this Group
INIADD = F
************************************************************
Group 13. Boundary & Special Sources
INLET (IN ,WEST ,2,0,0,0,0,0,1,1)
VALUE (IN ,P1 , 2.000E+00)
VALUE (IN ,U1 , 2.000E+00)
VALUE (IN ,V1 , 5.000E-01)
PATCH (OUTLET ,EAST ,3,0,0,0,0,0,1,1)
COVAL (OUTLET ,P1 , FIXVAL , 0.000E+00)
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 400
SARAH = 5.000E-02
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
************************************************************
Group 18. Limits
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
ASAP = T
LSG1 = T
************************************************************
Group 22. Monitor Print-Out
IXMON = 5 ;IYMON = 14 ;IZMON = 1
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NXPRIN = 1
NYPRIN = 1
IPLTF = 1 ;IPLTL = 40
PATCH (DOMAIN ,CONTUR,1,40,1,18,1,1,1,1)
PLOT(DOMAIN ,P1 , 0.000E+00, 1.000E+01)
PLOT(DOMAIN ,U1 , 0.000E+00, 1.000E+01)
PLOT(DOMAIN ,V1 , 0.000E+00, 1.000E+01)
PLOT(DOMAIN ,RHO1, 0.000E+00, 1.000E+01)
PATCH (IYEQNY ,PROFIL,1,40,18,18,1,1,1,1)
PLOT(IYEQNY ,P1 , 0.000E+00, 0.000E+00)
PATCH (IYEQ1 ,PROFIL,1,40,1,1,1,1,1,1)
PLOT(IYEQ1 ,P1 , 0.000E+00, 0.000E+00)
************************************************************
Group 24. Dumps For Restarts
NOWIPE = T
GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 5.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, MONIT, 4.750000E-01, 7.500000E-01, 5.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> DOM, RELAX, 1.000000E+00
> OBJ1, NAME, X1
> OBJ1, POSITION, 1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ1, SIZE, 0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ1, CLIPART, wirexyz
> OBJ1, ROTATION, 1
> OBJ1, TYPE, NULL
> OBJ2, NAME, IN
> OBJ2, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ2, SIZE, 0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ2, CLIPART, cube3
> OBJ2, ROTATION, 1
> OBJ2, TYPE, USER_DEFINED
> OBJ3, NAME, OUTLET
> OBJ3, POSITION, 5.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ3, SIZE, 0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ3, CLIPART, cube3
> OBJ3, ROTATION, 1
> OBJ3, TYPE, USER_DEFINED
> OBJ4, NAME, B4
> OBJ4, POSITION, 5.600000E-01, 0.000000E+00, 0.000000E+00
> OBJ4, SIZE, 1.130000E+00, 4.400001E-01, 1.000000E+00
> OBJ4, CLIPART, wirexyz
> OBJ4, ROTATION, 1
> OBJ4, TYPE, FINE_GRID_VOL
> OBJ4, FINEFACTS, 4.000000E+00, 4.000000E+00, 1.000000E+00
> OBJ5, NAME, B5
> OBJ5, POSITION, 1.470000E+00, 5.599999E-01, 0.000000E+00
> OBJ5, SIZE, 1.130000E+00, 4.400001E-01, 1.000000E+00
> OBJ5, CLIPART, wirexyz
> OBJ5, ROTATION, 1
> OBJ5, TYPE, FINE_GRID_VOL
> OBJ5, FINEFACTS, 4.000000E+00, 4.000000E+00, 1.000000E+00
M
distil=t
EX(P1 )= 1.892E-01
EX(U1 )= 2.003E+00
EX(V1 )= 1.619E-01
EX(VABS)= 2.013E+00
EX(PTOT)= 4.948E+00
EX(MACH)= 1.707E+00
EX(RHO1)= 1.000E+00