TALK=T;RUN(1,1)
PHOTON USE
p;;;;
up x; view z; Gr ou z 1
msg temperature contours
con tem1 z 1 fi;0.001; Pause;con off;red
msg velocity vectors
vec z 1 sh; pause
ENDUSE
DISPLAY
The problem considered is the 2d thermal convection of air
between isothermal horizontal walls in an enclosure of 5:1 aspect
ratio with adiabatic side walls. The lower wall is at a higher
temperature than the upper wall. This configuration frequently
occurs in technology, for example, in energy storage systems and
in many industrial heating and cooling systems. For Rayleigh
numbers Ra less than the critical value of 1708, no motion
occurs. For Ra in the range 1708 < Ra < 4E5, the fluid motion
consists of regularly spaced roll cells, while for larger Ra,
the cells break down and the motion is turbulent.
PHOENICS simulations may be performed for laminar flow (the
default) at a Rayleigh number Ra of 6.468E4 or turbulent flow at
Ra=2.35E6. The Boussinesq approximation is used and the energy
equation is solved with the temperature TEM1 as dependent
variable. For turbulent flow, the low-Re k-e model of Lam and
Bremhorst is used with buoyancy extensions and a turbulence
buoyancy coefficient of C3e=1.44.
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(2D Rayleigh-Benard Convection )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 297
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 28
XULAST =1.
XFRAC(1)=5.896926E-04 ;XFRAC(6)=5.810251E-03
XFRAC(11)=0.016416 ;XFRAC(16)=0.03083
XFRAC(21)=0.042543 ;XFRAC(26)=0.048698
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 42
YVLAST =1.
YFRAC(1)=1.902297E-03 ;YFRAC(9)=0.029051
YFRAC(17)=0.086583 ;YFRAC(25)=0.163417
YFRAC(33)=0.220949 ;YFRAC(41)=0.248098
************************************************************
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(149)=PRPS
NAME(150)=TEM1
* 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,N,N,N,Y)
SOLUTN(U1,Y,Y,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,Y)
SOLUTN(TEM1,Y,Y,N,N,N,Y)
PRPS = 149
************************************************************
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(TEM1,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.157136 ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =3.273322E-03 ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1.761505E-05 ;ENUT =0.
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(TEM1)=-0.026657
PRT(U1)=1. ;PRT(V1)=1.
PRT(TEM1)=0.9
CP1 =1008. ;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(PRPS)=-1.
FIINIT(TEM1)=305.5
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, 42, 1, 1, 1, 1)
COVAL(HOT ,V1 ,1. ,0. )
COVAL(HOT ,TEM1, GRND2 ,308. )
PATCH(COLD ,EWALL , 28, 28, 1, 42, 1, 1, 1, 1)
COVAL(COLD ,V1 ,1. ,0. )
COVAL(COLD ,TEM1, GRND2 ,303. )
PATCH(TOPWAL ,NWALL , 1, 28, 42, 42, 1, 1, 1, 1)
COVAL(TOPWAL ,U1 ,1. ,0. )
PATCH(BOTWAL ,SWALL , 1, 28, 1, 1, 1, 1, 1, 1)
COVAL(BOTWAL ,U1 ,1. ,0. )
PATCH(BUOY ,PHASEM, 1, 28, 1, 42, 1, 1, 1, 1)
COVAL(BUOY ,U1 , FIXFLU , GRND3 )
PATCH(REFP ,CELL , 14, 14, 21, 21, 1, 1, 1, 1)
COVAL(REFP ,P1 ,1. ,0. )
COVAL(REFP ,V1 ,0. ,0. )
COVAL(REFP ,TEM1,0. , SAME )
XCYCLE = F
EGWF = T
WALLCO = GRND2
BUOYA =-9.81 ; BUOYB =0. ; BUOYC =0.
BUOYD =3.273322E-03
BUOYE =305.5
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 800 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=2
LITER(V1)=2 ;LITER(TEM1)=50
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(TEM1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,0.099651)
RELAX(V1,FALSDT,0.099651)
RELAX(PRPS,LINRLX,1.)
RELAX(TEM1,FALSDT,1.0E+09)
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(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
VARMAX(TEM1)=308. ;VARMIN(TEM1)=303.
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=3.703E-04 ;EX(U1)=6.746E-03
EX(V1)=0.01186 ;EX(PRPS)=0.
EX(TEM1)=305.299988
************************************************************
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(PRPS,Y,N,Y,N,N,N)
OUTPUT(TEM1,Y,N,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 5 ;IYMON = 21 ;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 = 5 ;IXPRF = 1 ;IXPRL = 10000
NYPRIN = 8 ;IYPRF = 1 ;IYPRL = 10000
IPLTF = 1 ;IPLTL = -1 ;NPLT = 50
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