TALK=T;RUN(1,1)
DISPLAY
In the early years of PHOENICS, the energy equation was solved
with enthalpy as the dependent variable; temperature was
treated as an auxiliary variable, the values of which were to
be deduced from those of enthalpy by way of the specific heat
capacity of the fluid.
However, this was troublesome when so-called 'conjugate heat
transfer' problems had to be solved with heat-conducting
solid materials immersed in fluids; for each material could
have a different specific heat.
This library case was introduced at the time at which solving
for temperature directly was introduced, this being achieved
by use of the material-marker variable PRPS, wht value of
which at each point indicated what material, and therefore
what specific heat, was to be used there.
The steady flow in a cavity with one moving wall is modelled.
The stationary walls are at one temperature, and the wall
moving with a constant velocity is adiabatic
Within the cavity is a heated block of steel. The fluid is
water.
ENDDIS
PHOTON USE
p;;;
gr ou z 1
msg contours of PRPS
set prop off
con prps z 1 fi;0.001
pause;con off;red
msg Vector field
vec z 1 sh
msg Press RETURN to continue
pause;vec off;red
con tem1 z 1 fi;0.001
con tem1 z 1; int 50
msg Temperature contours
ENDUSE
DISPLAY
The geometry is as for case 921, but the north wall no longer
moves. Instead, the motion is induced by allowing for the
effect of buoyancy, by way of the Boussinesq approximation.
The PHOTON USE macro of case 921 can be activated by issuing the
command 'use Q1EAR'; but not 'use Q1'. The reason is that the
date of case 921 have been brought in by way of the compact
#$921 command above.
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Free Conv.; Steel Block In Water )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 922
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 24
XULAST =0.1
XFRAC(1)=0.041667 ;XFRAC(2)=0.083333
XFRAC(3)=0.125 ;XFRAC(4)=0.166667
XFRAC(5)=0.208333 ;XFRAC(6)=0.25
XFRAC(7)=0.291667 ;XFRAC(8)=0.333333
XFRAC(9)=0.375 ;XFRAC(10)=0.416667
XFRAC(11)=0.458333 ;XFRAC(12)=0.5
XFRAC(13)=0.541667 ;XFRAC(14)=0.583333
XFRAC(15)=0.625 ;XFRAC(16)=0.666667
XFRAC(17)=0.708333 ;XFRAC(18)=0.75
XFRAC(19)=0.791667 ;XFRAC(20)=0.833333
XFRAC(21)=0.875 ;XFRAC(22)=0.916667
XFRAC(23)=0.958333 ;XFRAC(24)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 24
YVLAST =0.1
YFRAC(1)=0.041667 ;YFRAC(2)=0.083333
YFRAC(3)=0.125 ;YFRAC(4)=0.166667
YFRAC(5)=0.208333 ;YFRAC(6)=0.25
YFRAC(7)=0.291667 ;YFRAC(8)=0.333333
YFRAC(9)=0.375 ;YFRAC(10)=0.416667
YFRAC(11)=0.458333 ;YFRAC(12)=0.5
YFRAC(13)=0.541667 ;YFRAC(14)=0.583333
YFRAC(15)=0.625 ;YFRAC(16)=0.666667
YFRAC(17)=0.708333 ;YFRAC(18)=0.75
YFRAC(19)=0.791667 ;YFRAC(20)=0.833333
YFRAC(21)=0.875 ;YFRAC(22)=0.916667
YFRAC(23)=0.958333 ;YFRAC(24)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =0.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,Y,N,N,Y)
SOLUTN(U1,Y,Y,Y,N,N,Y)
SOLUTN(V1,Y,Y,Y,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,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. ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1.0E-05 ;ENUT =0.
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(TEM1)=1.
PRT(U1)=1. ;PRT(V1)=1.
PRT(TEM1)=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(PRPS)=67.
FIINIT(TEM1)=20.
PATCH(SOLID ,INIVAL, 7, 18, 7, 18, 1, 1, 1, 1)
INIT(SOLID ,PRPS,0. ,111. )
INIT(SOLID ,TEM1,0. ,25. )
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(NORTHW ,NWALL , 1, 24, 24, 24, 1, 1, 1, 1)
COVAL(NORTHW ,U1 ,1. ,0. )
COVAL(NORTHW ,TEM1,1. ,20. )
PATCH(WESTW ,WWALL , 1, 1, 1, 24, 1, 1, 1, 1)
COVAL(WESTW ,V1 ,1. ,0. )
COVAL(WESTW ,TEM1,1. ,20. )
PATCH(EASTW ,EWALL , 24, 24, 1, 24, 1, 1, 1, 1)
COVAL(EASTW ,V1 ,1. ,0. )
COVAL(EASTW ,TEM1,1. ,20. )
PATCH(RELIEF ,CELL , 1, 1, 1, 1, 1, 1, 1, 1)
COVAL(RELIEF ,P1 ,1. ,0. )
PATCH(HEATEDBL,VOLUME, 7, 18, 7, 18, 1, 1, 1, 1)
COVAL(HEATEDBL,TEM1, FIXFLU ,1.0E+04 )
PATCH(BUOYANCY,PHASEM, 1, 24, 1, 24, 1, 1, 1, 1)
COVAL(BUOYANCY,V1 , FIXFLU , GRND3 )
XCYCLE = F
EGWF = T
WALLCO = GRND2
BUOYA =0. ; BUOYB =-9.81 ; BUOYC =0.
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1000 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=10
LITER(V1)=10 ;LITER(TEM1)=20
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,1.)
RELAX(V1,FALSDT,1.)
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)=1.0E+10 ;VARMIN(TEM1)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
CONWIZ = T
ISG50 = 1
ISG52 = 1
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,MATERIAL,67,L,T)
SPEDAT(SET,MATERIAL,111,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,Y,Y,Y,Y,Y)
OUTPUT(U1,Y,Y,Y,Y,Y,Y)
OUTPUT(V1,Y,Y,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 = 13 ;IYMON = 2 ;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 = -1 ;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(MIDDLE ,PROFIL, 13, 13, 1, 24, 1, 1, 1, 1)
PLOT(MIDDLE ,U1 ,-0.1 ,0.1 )
PLOT(MIDDLE ,TEM1,0. ,0. )
PATCH(MAP ,CONTUR, 1, 24, 1, 24, 1, 1, 1, 1)
PLOT(MAP ,U1 ,0. ,10. )
PLOT(MAP ,V1 ,0. ,10. )
PLOT(MAP ,TEM1,0. ,10. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP