TALK=T;RUN(1,1)
DISPLAY
Simulation of laminar steady flow in z direction around 2 spheres
in a uniform stream of atmospheric air, treated as incompressible.
The basic grid is cartesian and rather coarse (18*18*31), but the
flow follows the curved surface smoothly because the PARSOL (i.e.
cut-cell) technique is used.
Provision is made for:
1. Solving for only one quarter of the domain; this is allowed
by reason of symmetry and desirable for economy and accuracy.
2. Refining the grid near the spheres to increase accuracy more.
3. Varying the Reynolds number.
The variables computed are:
* p1, the pressure,
* u1, the x-direction velocity,
* v1, the y-direction velocity, and
* w1, the z-direction velocity.
This file contains a macro of commands which cause the Viewer
(when the macro button is pressed) to display contours
automatically.
[Case derived from core linrary case 805 by dbs 11.09.08]
ENDDIS
VRV USE
* Start of frame
VARIABLE Pressure; CON ON
msg Pressure contours
PAUSE
* Start of frame
VARIABLE Velocity; VEC ON
msg Pressure contours superimposed by velocity vectors
PAUSE
* Start of frame
* Setting object visibility and painting status
OBJECT SHOW TYPE BLOCKAGE
OBJECT PAINT TYPE BLOCKAGE ON
VARIABLE Velocity; VEC OFF
CONTOUR OFF
msg Velocity contours on the surface of the sphere
PAUSE
* Start of frame
* Setting object visibility and painting status
OBJECT SHOW TYPE BLOCKAGE
OBJECT PAINT TYPE BLOCKAGE ON
VARIABLE Pressure; VEC OFF
msg Pressure contours on the surface of the sphere
PAUSE
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(2 spheres Re=40. quarter=T finegrid=T )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 18
XULAST =2.
XFRAC(1)=0.0625 ;XFRAC(2)=0.125
XFRAC(3)=0.1875 ;XFRAC(4)=0.25
XFRAC(5)=0.3 ;XFRAC(6)=0.35
XFRAC(7)=0.4 ;XFRAC(8)=0.45
XFRAC(9)=0.5 ;XFRAC(10)=0.55
XFRAC(11)=0.6 ;XFRAC(12)=0.65
XFRAC(13)=0.7 ;XFRAC(14)=0.76
XFRAC(15)=0.82 ;XFRAC(16)=0.88
XFRAC(17)=0.94 ;XFRAC(18)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 18
YVLAST =2.
YFRAC(1)=0.0625 ;YFRAC(2)=0.125
YFRAC(3)=0.1875 ;YFRAC(4)=0.25
YFRAC(5)=0.3 ;YFRAC(6)=0.35
YFRAC(7)=0.4 ;YFRAC(8)=0.45
YFRAC(9)=0.5 ;YFRAC(10)=0.55
YFRAC(11)=0.6 ;YFRAC(12)=0.65
YFRAC(13)=0.7 ;YFRAC(14)=0.76
YFRAC(15)=0.82 ;YFRAC(16)=0.88
YFRAC(17)=0.94 ;YFRAC(18)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 31
ZWLAST =10.
ZFRAC(1)=0.0364 ;ZFRAC(2)=0.0728
ZFRAC(3)=0.1092 ;ZFRAC(4)=0.1456
ZFRAC(5)=0.182 ;ZFRAC(6)=0.202
ZFRAC(7)=0.222 ;ZFRAC(8)=0.255333
ZFRAC(9)=0.288667 ;ZFRAC(10)=0.322
ZFRAC(11)=0.355333 ;ZFRAC(12)=0.388667
ZFRAC(13)=0.422 ;ZFRAC(14)=0.442
ZFRAC(15)=0.462 ;ZFRAC(16)=0.4988
ZFRAC(17)=0.5356 ;ZFRAC(18)=0.5724
ZFRAC(19)=0.6092 ;ZFRAC(20)=0.646
ZFRAC(21)=0.666 ;ZFRAC(22)=0.699333
ZFRAC(23)=0.732667 ;ZFRAC(24)=0.766
ZFRAC(25)=0.786 ;ZFRAC(26)=0.821667
ZFRAC(27)=0.857333 ;ZFRAC(28)=0.893
ZFRAC(29)=0.928667 ;ZFRAC(30)=0.964333
ZFRAC(31)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(7)=W1
NAME(148)=STGP ;NAME(149)=VLSQ
NAME(150)=PRPS
* 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,Y)
SOLUTN(U1,Y,Y,Y,N,N,Y)
SOLUTN(V1,Y,Y,Y,N,N,Y)
SOLUTN(W1,Y,Y,Y,N,N,Y)
SOLUTN(STGP,Y,N,N,N,N,N)
SOLUTN(VLSQ,Y,N,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,Y)
PRPS = 150
************************************************************
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(U1,Y,Y,Y,Y,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1. ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =0.05 ;ENUT =0.
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(W1)=1.
PRT(U1)=1. ;PRT(V1)=1.
PRT(W1)=1.
CP1 =1005. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10
FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.
FIINIT(STGP)=1.0E-10 ;FIINIT(VLSQ)=1.0E-10
FIINIT(PRPS)=-1.
PATCH(^OB5 ,INIVAL, 5, 0, 0, 0, 0, 0, 1, 1)
INIT(^OB5 ,PRPS,0. ,198. )
PATCH(^OB6 ,INIVAL, 6, 0, 0, 0, 0, 0, 1, 1)
INIT(^OB6 ,PRPS,0. ,198. )
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(IMBL3&2 ,VOLUME, 3, 16, 3, 16, 3, 29, 1, 1)
COVAL(IMBL3&2 ,W1 ,0. ,0. )
PATCH(IMBL4&3 ,VOLUME, 3, 16, 3, 16, 4, 28, 1, 1)
COVAL(IMBL4&3 ,W1 ,0. ,0. )
PATCH(IMBL5&4 ,VOLUME, 3, 16, 3, 16, 5, 27, 1, 1)
COVAL(IMBL5&4 ,W1 ,0. ,0. )
PATCH(IMBL810 ,VOLUME, 3, 16, 3, 16, 8, 21, 1, 1)
COVAL(IMBL810 ,W1 ,0. ,0. )
PATCH(IMBL815 ,VOLUME, 5, 14, 5, 14, 8, 16, 1, 1)
COVAL(IMBL815 ,W1 ,0. ,0. )
PATCH(^OB1 ,CELL , 1, 0, 0, 0, 0, 0, 1, 1)
PATCH(^OB2 ,CELL , 2, 0, 0, 0, 0, 0, 1, 1)
PATCH(^OB3 ,LOW , 3, 0, 0, 0, 0, 0, 1, 1)
COVAL(^OB3 ,P1 , FIXFLU ,1. )
COVAL(^OB3 ,U1 ,0. ,0. )
COVAL(^OB3 ,V1 ,0. ,0. )
COVAL(^OB3 ,W1 ,0. ,1. )
PATCH(^OB4 ,HIGH , 4, 0, 0, 0, 0, 0, 1, 1)
COVAL(^OB4 ,P1 ,1000. ,0. )
COVAL(^OB4 ,U1 ,0. ,0. )
COVAL(^OB4 ,V1 ,0. ,0. )
COVAL(^OB4 ,W1 ,0. ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 30 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
LITER(P1)=250 ;LITER(U1)=10
LITER(V1)=10 ;LITER(W1)=10
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,1.)
RELAX(V1,FALSDT,1.)
RELAX(W1,FALSDT,1.)
RELAX(STGP,LINRLX,1.)
RELAX(VLSQ,LINRLX,1.)
RELAX(PRPS,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(STGP)=1.0E+10 ;VARMIN(STGP)=-1.0E+10
VARMAX(VLSQ)=1.0E+10 ;VARMIN(VLSQ)=-1.0E+10
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
ISG50 = 1
ISG52 = 2
ISG62 = 0
SPEDAT(SET,STORED,STGP,C,=P1+0.5*RHO1*VLSQ!SWPFIN)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,OBJNAM,^OB1,C,G1)
SPEDAT(SET,OBJTYP,^OB1,C,FINE_GRID_VOL)
SPEDAT(SET,^OB1,FINEX,I,4)
SPEDAT(SET,^OB1,FINEY,I,4)
SPEDAT(SET,^OB1,FINEZ,I,4)
SPEDAT(SET,OBJNAM,^OB2,C,G2)
SPEDAT(SET,OBJTYP,^OB2,C,FINE_GRID_VOL)
SPEDAT(SET,^OB2,FINEX,I,4)
SPEDAT(SET,^OB2,FINEY,I,4)
SPEDAT(SET,^OB2,FINEZ,I,4)
SPEDAT(SET,OBJNAM,^OB3,C,AIRIN)
SPEDAT(SET,OBJTYP,^OB3,C,INLET)
SPEDAT(SET,OBJNAM,^OB4,C,AIROUT)
SPEDAT(SET,OBJTYP,^OB4,C,OUTLET)
SPEDAT(SET,ARATIO,^OB4,R,1.)
SPEDAT(SET,OBJNAM,^OB5,C,FIRSTSPH)
SPEDAT(SET,OBJTYP,^OB5,C,BLOCKAGE)
SPEDAT(SET,FIRSTSPH,MATERIAL,R,198.)
SPEDAT(SET,OBJNAM,^OB6,C,SECONDSP)
SPEDAT(SET,OBJTYP,^OB6,C,BLOCKAGE)
SPEDAT(SET,SECONDSP,MATERIAL,R,198.)
SPEDAT(SET,FACETDAT,NUMOBJ,I,6)
SPEDAT(SET,MATERIAL,198,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=0.1298 ;EX(U1)=0.0324
EX(V1)=0.03335 ;EX(W1)=1.012
EX(STGP)=0. ;EX(VLSQ)=0.
EX(PRPS)=3.354
************************************************************
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(U1,Y,Y,Y,N,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(STGP,Y,N,Y,N,N,N)
OUTPUT(VLSQ,Y,N,Y,N,N,N)
OUTPUT(PRPS,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 9 ;IYMON = 9 ;IZMON = 17
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
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