DISPLAY
FREEZING WATER IN A LID-DRIVEN CAVITY FLOW
The simplified porosity method is introduced to simulate
how the water flow is affected by solid-liquid phase change.
It is assumed that the phase change is occured over a finite
range of temperatures defined by RG(1) and RG(2) which
are used to calculate a volumetric liquid fraction, FLIQ.
An additional source term appears in the momentum equation
to reduce gradually the velocities in the solidifying zone.
No account is taken for neither the latent heat effect nor
differences in phase physical properties.
A full account of enthalpy-porosity method for phase-change
problem and its PHOENICS implementation can be found in
V.R. Voller and C. Prakash 'A fixed grid numerical modelling
methodology for convection-diffusion mushy region phase-change
problems', Int. Journal of Heat and Mass Transfer, v.30, 1987,
pp. 1709-1720.
ENDDIS
PHOTON USE
p
gr ou z 1
con fliq z 1 fi;0.001
vec z 1
msg Velocity vectors and liquid volume fraction.
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(FREEZING WATER IN A LID-DRIVEN CAVITY:Z106 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 15
XULAST = 1.000000E-01
XFRAC ( 1) = 6.666667E-02 ;XFRAC ( 2) = 1.333333E-01
XFRAC ( 3) = 2.000000E-01 ;XFRAC ( 4) = 2.666667E-01
XFRAC ( 5) = 3.333333E-01 ;XFRAC ( 6) = 4.000000E-01
XFRAC ( 7) = 4.666667E-01 ;XFRAC ( 8) = 5.333334E-01
XFRAC ( 9) = 6.000000E-01 ;XFRAC (10) = 6.666667E-01
XFRAC (11) = 7.333333E-01 ;XFRAC (12) = 8.000000E-01
XFRAC (13) = 8.666667E-01 ;XFRAC (14) = 9.333333E-01
XFRAC ( 15) = 1.000000E+00
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 15
YVLAST = 1.000000E-01
YFRAC ( 1) = 6.666667E-02 ;YFRAC ( 2) = 1.333333E-01
YFRAC ( 3) = 2.000000E-01 ;YFRAC ( 4) = 2.666667E-01
YFRAC ( 5) = 3.333333E-01 ;YFRAC ( 6) = 4.000000E-01
YFRAC ( 7) = 4.666667E-01 ;YFRAC ( 8) = 5.333334E-01
YFRAC ( 9) = 6.000000E-01 ;YFRAC (10) = 6.666667E-01
YFRAC (11) = 7.333333E-01 ;YFRAC (12) = 8.000000E-01
YFRAC (13) = 8.666667E-01 ;YFRAC (14) = 9.333333E-01
YFRAC ( 15) = 1.000000E+00
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST = 1.000000E+00
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
SOLVE(P1,U1,V1,TEMP)
STORE(FLIQ,VISL)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
RHO1 = 1000.
ENUL = 1.000000E-03
PRNDTL(TEMP) = 5.0
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
************************************************************
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
Parent VR object for this patch is: RELIEF
PATCH(RELIEF ,CELL , 7, 7, 7, 7, 1, 1, 1, 1)
COVAL(RELIEF ,P1 , 1.000000E+03, 0.000000E+00)
COVAL(RELIEF ,U1 , 0.000000E+00, 0.000000E+00)
COVAL(RELIEF ,V1 , 0.000000E+00, 0.000000E+00)
COVAL(RELIEF ,TEMP, 0.000000E+00, 0.000000E+00)
Parent VR object for this patch is: SOUTH
PATCH(SOUTH ,SWALL , 1, 15, 1, 1, 1, 1, 1, 1)
COVAL(SOUTH ,U1 , 1.000000E+00, 0.000000E+00)
Parent VR object for this patch is: MOVING
PATCH(MOVING ,NWALL , 1, 15, 15, 15, 1, 1, 1, 1)
COVAL(MOVING ,U1 , 1.000000E+00, 1.000000E+00)
Parent VR object for this patch is: WEST
PATCH(HOT ,WWALL , 1, 1, 1, 15, 1, 1, 1, 1)
COVAL(HOT ,V1 , 1.000000E+00, 0.000000E+00)
COVAL(HOT ,TEMP, 1.000000E+00, 2.830000E+02)
Parent VR object for this patch is: EAST
PATCH(COLD ,EWALL , 15, 15, 1, 15, 1, 1, 1, 1)
COVAL(COLD ,V1 , 1.000000E+00, 0.000000E+00)
COVAL(COLD ,TEMP, 1.000000E+00, 2.630000E+02)
PATCH(MUSHY ,VOLUME, 1, 15, 1, 15, 1, 1, 1, 1)
COVAL(MUSHY ,U1 , GRND , 0.000000E+00)
COVAL(MUSHY ,V1 , GRND , 0.000000E+00)
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
SELREF=F
LSWEEP = 500
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 2.000000E-01)
RELAX(U1 ,FALSDT, 2.000000E-04)
RELAX(V1 ,FALSDT, 2.000000E-04)
RELAX(FLIQ,LINRLX, 1.000000E+00)
RELAX(VISL,LINRLX, 1.000000E+00)
RELAX(TEMP,LINRLX, 1.000000E+00)
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
NAMSAT =MOSG
RG( 1) = 2.710000E+02
RG( 2) = 2.750000E+02
RG( 3) = 1.000000E+02
************************************************************
Group 20. Preliminary Printout
************************************************************
tstswp=-1
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 7 ;IYMON = 14 ;IZMON = 1
************************************************************
Group 23.Field Print-Out & Plot Control
************************************************************
Group 24. Dumps For Restarts
PLANTBEGIN
PATCH (MUSHY,VOLUME,1,NX,1,NY,1,NZ,1,1)
CO=RG(3)*(1-FLIQ)/FLIQ
COVAL (MUSHY,U1 , GRND , 0.0 )
CO=RG(3)*(1-FLIQ)/FLIQ
COVAL (MUSHY,V1 , GRND , 0.0 )
STORE(FLIQ)
FLIQ=AMAX1(1.e-5,AMIN1(1.,(TEMP-RG(1))/(RG(2)-RG(1))))
PLANTEND
dmpstk=t
DISTIL=T
EX(P1)=2.878E+01; EX(U1)=8.870E-02; EX(V1)=4.645E-02
EX(VISL)=1.000E-03; EX(FLIQ)=6.455E-01; EX(TEMP)=2.736E+02
LIBREF=106
STOP