TALK=T;RUN(1,1)
PHOTON USE
p
msg the grid. press return for temperature contours
gr z 1;pause; gr off;red;gr ou z 1
msg temperature contours. press return for velocity vectors
con temp z 1 fi;0.0002;pause; con off;red;se re ve 2;vec z 1 sh
msg velocity vectors. press return for reduced-pressure contours
pause;vec off;red;con p1 z 1 fi;0.0002
ENDUSE
DISPLAY
A two-dimensional square cavity is formed between two
vertical walls, one of which is heated and the other cooled.
the top and bottom of the cavity are bounded by walls at which
there is friction but no heat transfer.
The left wall is the hot wall of a temperature Thot=1 K
The right wall is the cold wall of a temperature Tcold=0 K
ENDDIS
DISPLAY
Parameters:
-----------
1) The flow characteristics are set by way of
Length of the cavity=CAVL
Prandtl Number=Pr=enul/diffusity
Rayleigh Number=Ra=Beta*g*(Th-Tc)/(enul*alfa)
The Thermal Expansion Coefficient Beta as well as
the Thermal diffusivity Alfa remain constant.
2) This library case also enables the user to select
the numerical solution procedure, namely:
-the choice between two solvers:
STONE or CNGR (Conjugate-residuals gradient)
-the possibilty of the activation of EXPERT
the Self-Adjusting-Relaxation Algorithm.
ENDDIS
DISPLAY
Library case 30 loaded, to activate EXPERT
expsol f .... expert is not used in the solver
expdtf t .... expert is used to adjust false time step
indtf 2 ... adjust the false time step of all velocities
inres 3 ... monitor the residuals of variable 3
ifrsts 20 ... first adjustment sweep
ifrequ 10 ... adjustments will be made every ifrequ sweeps
facdec 0.5 .. decreases will be by a factor of facdec
facinc 2.0 .. increases will be by a factor of facinc
exprin t .... print adjusted false time steps in RESULT file
expend ..... end of the inputs to expert via Q1
xprtend
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Laminar Free Convection In Cavity )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 26
XULAST =1.
XFRAC(1)=0.038462 ;XFRAC(2)=0.076923
XFRAC(3)=0.115385 ;XFRAC(4)=0.153846
XFRAC(5)=0.192308 ;XFRAC(6)=0.230769
XFRAC(7)=0.269231 ;XFRAC(8)=0.307692
XFRAC(9)=0.346154 ;XFRAC(10)=0.384615
XFRAC(11)=0.423077 ;XFRAC(12)=0.461538
XFRAC(13)=0.5 ;XFRAC(14)=0.538462
XFRAC(15)=0.576923 ;XFRAC(16)=0.615385
XFRAC(17)=0.653846 ;XFRAC(18)=0.692308
XFRAC(19)=0.730769 ;XFRAC(20)=0.769231
XFRAC(21)=0.807692 ;XFRAC(22)=0.846154
XFRAC(23)=0.884615 ;XFRAC(24)=0.923077
XFRAC(25)=0.961538 ;XFRAC(26)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 26
YVLAST =1.
YFRAC(1)=0.038462 ;YFRAC(6)=0.230769
YFRAC(11)=0.423077 ;YFRAC(16)=0.615385
YFRAC(21)=0.807692 ;YFRAC(26)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =1.
ZFRAC(1)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(14)=TEMP
* 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(U1,Y,Y,Y,N,N,N)
SOLUTN(V1,Y,Y,Y,N,N,N)
SOLUTN(TEMP,Y,Y,N,N,N,Y)
************************************************************
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(TEMP,N,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.207 ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =1.0E-02 ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =0.71 ;ENUT =0.
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(TEMP)=0.71
PRT(U1)=1. ;PRT(V1)=1.
PRT(TEMP)=1.
CP1 =1. ;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(TEMP)=1.0E-10
No PATCHes yet used for this Group
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(HOT ,WWALL , 1, 1, 1, 26, 1, 1, 1, 1)
COVAL(HOT ,V1 ,1. ,0. )
COVAL(HOT ,TEMP,1. ,1. )
PATCH(COLD ,EWALL , 26, 26, 1, 26, 1, 1, 1, 1)
COVAL(COLD ,V1 ,1. ,0. )
COVAL(COLD ,TEMP,1. ,0. )
PATCH(LOWAL ,SWALL , 1, 26, 1, 1, 1, 1, 1, 1)
COVAL(LOWAL ,U1 ,1. ,0. )
PATCH(HIWAL ,NWALL , 1, 26, 26, 26, 1, 1, 1, 1)
COVAL(HIWAL ,U1 ,1. ,0. )
PATCH(BUOY ,PHASEM, 1, 26, 1, 26, 1, 1, 1, 1)
COVAL(BUOY ,V1 , FIXFLU , GRND3 )
PATCH(REFP ,CELL , 13, 13, 13, 13, 1, 1, 1, 1)
COVAL(REFP ,P1 ,1. ,0. )
COVAL(REFP ,U1 ,0. ,0. )
COVAL(REFP ,V1 ,0. ,0. )
COVAL(REFP ,TEMP,0. , SAME )
XCYCLE = F
EGWF = T
WALLCO = GRND2
BUOYA =0. ; BUOYB =-7.1E+04 ; BUOYC =0.
BUOYE =0.5
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 200 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=2
LITER(V1)=2 ;LITER(TEMP)=20
ENDIT(P1)=1.0E-06 ;ENDIT(U1)=1.0E-06
ENDIT(V1)=1.0E-06 ;ENDIT(TEMP)=1.0E-06
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,1.0E-03)
RELAX(V1,FALSDT,1.0E-03)
RELAX(TEMP,FALSDT,1.)
OVRRLX =0.
EXPERT = T ;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(TEMP)=1.0E+10 ;VARMIN(TEMP)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG62 = 1
CSG3 ='CNGR'
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(U1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(TEMP,Y,N,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 6 ;IYMON = 6 ;IZMON = 1
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 = 5 ;IYPRF = 1 ;IYPRL = 10000
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
PATCH(PROF ,PROFIL, 1, 26, 1, 26, 0, 0, 1, 1)
PLOT(PROF ,TEMP,0. ,1. )
PATCH(CONT ,CONTUR, 1, 26, 1, 26, 1, 1, 1, 1)
PLOT(CONT ,TEMP,0. ,10. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = T
NSAVE =CHAM
STOP