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
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Square Cavity; Steel Block; Water       )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 921
 ************************************************************
  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(MOVING ,SWALL , 1, 24, 1, 1, 1, 1, 1, 1)
 COVAL(MOVING ,U1 ,1. ,-0.1 )
 
 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 )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  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