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