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
  DISPLAY
  This case is like case 922; but the fluid is glycerine
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Free Conv.; Steel Block In Glycerine    )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 924
 ************************************************************
  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)=65.
 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,65,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