#$r002
   PHOTON USE
   AUTOPLOT
   file
   PHI 5
 
   cl
   msg Temperature profile
   msg Analytical solution shows that TSL=842K (distance 0.5)
   da 1 tem1
   col3 1
   ENDUSE
 
   *  GROUP 1.  Run identifiers and other preliminaries.
TEXT(1-D THIN PLATE : R118
TITLE
  DISPLAY
   The test case considers internal black-body radiative heat
   transfer between the various surfaces, and conductive heat
   transfer through the air and a thin plate located mid-way
   between the two bounding surfaces of the solution domain.
        ___________________________________________________
      /|                        |/                         |/
      /|<---TL=1000K    TSL---->|/<---TSR           TR=293K|/
      /|                        |/                         |/
      /|       AIR              |/    AIR                  |/
      /|       PRPS=0           |/    PRPS=0               |/
      /|                        |/                         |/
      /|<---Radiative zone 1    |<----Radiative zone 2     |/
      /|    Radiative zone 3--->|/    Radiative zone 4---> |/
       |________________________|__________________________|
        ----------> Y
        length of domain = 1.0 m, height of domain = 1.0m
 
   The PHOENICS solution can be compared with the analytical
   solution which shows that TSL=842K=TSR, and  the net radiative
   flux at zone 3 =28.kW/m2
  ENDDIS
   ************************************************************
   *  User declared Real variables
   * NYG - the Y gas cell at the solid/fluid interface used in
           Group 12
   * NYG is set in Group 4
   * GCOND - the laminar thermal conductivity of the plate (W/m/K)
   * GPLTH - the thickness of the plate (m)
INTEGER(NYG)
REAL(GCOND,GPLTH);GCOND=220.;GPLTH=0.005
   ************************************************************
   *  GROUP 4.  y-direction grid specification.
NREGY=2;REGEXT(Y,5.0000E-01)
   * Region 1 for fluid
IREGY=1;GRDPWR(Y,5,5.0000E-01,-2.0)
NYG=5
   * Region 2 for Solid
IREGY=2;GRDPWR(Y,5,5.0000E-01,2.0)
   ************************************************************
   *  GROUP 7.  Variables (including porosities) named,
   *            stored & solved.
SOLVE(TEM1);SOLUTN(TEM1,Y,Y,Y,N,N,Y)
STORE(PRPS)
   ************************************************************
   *  GROUP 9.  Properties of the medium (or media).
   *  Fluid/medium is Air-ideal_gas
   *  Reference pressure (N/m^2)
SETPRPS(1,0)
PRESS0=1.0000E+05
   ************************************************************
   *  GROUP 11. Initialization of fields of variables,
   *            porosities, etc.
   * Thin plate
CONPOR(PLT,0.0,NORTH,#1,#1,-NYG,-NYG,#1,#1)
FIINIT(TEM1)=800.00
   ************************************************************
   *  GROUP 13. Boundary conditions and special sources
   * WALL boundary condition, name TFIXL to set the fixed TL
     for the conductive boundary condition
PATCH(TFIXL,SWALL,#1,#1,#1,#1,#1,#1,#1,#NREGT)
COVAL(TFIXL,TEM1,1.0,1000)
   * WALL boundary condition, name TFIXR to set the fixed TR
     for the conductive boundary condition
PATCH(TFIXR,NWALL,#1,#1,#2,#2,#1,#1,#1,#NREGT)
COVAL(TFIXR,TEM1,1.0,293)
   * Radiative zones have PATCH name of the form RI###, ###
     are digits indicating the zone number. RI indicates that
     the PATCH is for an internal radiating zone.
   * Radiation zone 1. This zone is located at the left
     wall and has a fixed temperature
PATCH(@RI001,SOUTH,#1,#1,#1,#1,#1,#1,#1,#NREGT)
COVAL(@RI001,TEM1,0.0,1000.)
   * Radiation zone 2 and 3 are surfaces of the plate.
     The surface temperature of these zones  are calculated
     in GXS2SR
     The COefficient is GRND3 and the VALue is
     k/d for the plate
PATCH(@RI002,SOUTH,#1,#1,NYG+1,NYG+1,#1,#1,#1,#nregt)
COVAL(@RI002,TEM1,GRND3,GCOND/GPLTH)
PATCH(@RI003,NORTH,#1,#1,NYG,NYG,#1,#1,#1,#1)
COVAL(@RI003,TEM1,GRND3,GCOND/GPLTH)
     * Fixed temperature radiative thermal zone 4
PATCH(@RI004,NORTH,#1,#1,#2,#2,#1,#1,#1,#1)
COVAL(@RI004,TEM1,0.0,293.)
   ************************************************************
   *  GROUP 15. Termination criteria for sweeps and
   *            outer iterations.
LSWEEP=2000;TSTSWP=-1
   ************************************************************
   *  GROUP 16. Termination criteria for inner iterations.
S2SR=T
   ************************************************************
   *  GROUP 19. Data communicated by SATELLITE to GROUND
IYMON=NY/2
   ************************************************************
   *  GROUP 24. Preparation for continuation runs.
 LIBREF  =       118
spedat(set,cvd,radcvd,l,t)
STOP