TALK=T;RUN(1,1)
  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
  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
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(FREEZING WATER IN A LID-DRIVEN CAVITY:Z1)
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 106
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 15
 XULAST =0.1
 XFRAC(1)=0.066667 ;XFRAC(2)=0.133333
 XFRAC(3)=0.2 ;XFRAC(4)=0.266667
 XFRAC(5)=0.333333 ;XFRAC(6)=0.4
 XFRAC(7)=0.466667 ;XFRAC(8)=0.533333
 XFRAC(9)=0.6 ;XFRAC(10)=0.666667
 XFRAC(11)=0.733333 ;XFRAC(12)=0.8
 XFRAC(13)=0.866667 ;XFRAC(14)=0.933333
 XFRAC(15)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 15
 YVLAST =0.1
 YFRAC(1)=0.066667 ;YFRAC(2)=0.133333
 YFRAC(3)=0.2 ;YFRAC(4)=0.266667
 YFRAC(5)=0.333333 ;YFRAC(6)=0.4
 YFRAC(7)=0.466667 ;YFRAC(8)=0.533333
 YFRAC(9)=0.6 ;YFRAC(10)=0.666667
 YFRAC(11)=0.733333 ;YFRAC(12)=0.8
 YFRAC(13)=0.866667 ;YFRAC(14)=0.933333
 YFRAC(15)=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(148)=VISL
 NAME(149)=FLIQ ;NAME(150)=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,N,N,N,Y)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(VISL,Y,N,N,N,N,Y)
 SOLUTN(FLIQ,Y,N,N,N,N,Y)
 SOLUTN(TEMP,Y,Y,N,N,N,Y)
 VISL = 148
 ************************************************************
  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 =1000. ;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-03 ;ENUT =0.
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(TEMP)=5.
 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(VISL)=1.0E-10
 FIINIT(FLIQ)=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(RELIEF ,CELL , 7, 7, 7, 7, 1, 1, 1, 1)
 COVAL(RELIEF ,P1 ,1000. ,0. )
 COVAL(RELIEF ,U1 ,0. ,0. )
 COVAL(RELIEF ,V1 ,0. ,0. )
 COVAL(RELIEF ,TEMP,0. ,0. )
 
 PATCH(SOUTH ,SWALL , 1, 15, 1, 1, 1, 1, 1, 1)
 COVAL(SOUTH ,U1 ,1. ,0. )
 
 PATCH(MOVING ,NWALL , 1, 15, 15, 15, 1, 1, 1, 1)
 COVAL(MOVING ,U1 ,1. ,1. )
 
 PATCH(HOT ,WWALL , 1, 1, 1, 15, 1, 1, 1, 1)
 COVAL(HOT ,V1 ,1. ,0. )
 COVAL(HOT ,TEMP,1. ,283. )
 
 PATCH(COLD ,EWALL , 15, 15, 1, 15, 1, 1, 1, 1)
 COVAL(COLD ,V1 ,1. ,0. )
 COVAL(COLD ,TEMP,1. ,263. )
 
 PATCH(MUSHY ,VOLUME, 1, 15, 1, 15, 1, 1, 1, 1)
 COVAL(MUSHY ,U1 , GRND ,0. )
 COVAL(MUSHY ,V1 , GRND ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 500 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = F
 RESREF(P1)=1.0E-08 ;RESREF(U1)=1.0E-08
 RESREF(V1)=1.0E-08 ;RESREF(TEMP)=1.0E-08
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(TEMP)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(TEMP)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.2)
 RELAX(U1,FALSDT,2.0E-04)
 RELAX(V1,FALSDT,2.0E-04)
 RELAX(VISL,LINRLX,1.)
 RELAX(FLIQ,LINRLX,1.)
 RELAX(TEMP,LINRLX,1.)
 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(VISL)=1.0E+10 ;VARMIN(VISL)=-1.0E+10
 VARMAX(FLIQ)=1.0E+10 ;VARMIN(FLIQ)=-1.0E+10
 VARMAX(TEMP)=1.0E+10 ;VARMIN(TEMP)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMSAT =MOSG
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 RG( 1) =271.
 RG( 2) =275.
 RG( 3) =100.
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=28.780001 ;EX(U1)=0.0887
 EX(V1)=0.04645 ;EX(VISL)=1.0E-03
 EX(FLIQ)=0.6455 ;EX(TEMP)=273.600006
 ************************************************************
  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(VISL,Y,N,Y,N,N,N)
 OUTPUT(FLIQ,Y,N,Y,N,N,N)
 OUTPUT(TEMP,Y,N,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 7 ;IYMON = 14 ;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
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP