#$r002
 
   PHOTON USE
   AUTOPLOT
   file
   PHI 5
 
   cl
   msg Temperature profile
   msg Analytical solution shows that TS=430K (middle)
   da 1 tem1
   col3 1
   le y 430
   le y 293
   le x .5
   pause
   ENDUSE
 
TEXT(1-D CARTESIAN Y-DIRECTION :   R115
  DISPLAY
   The test case considers combined radiation and conduction
   in the domain indicated below. The thin region on the left
   is solid aluminum, and the TL condition is at the left
   edge of this zone.
     ______________________________________________________
     |/|                        |/                         |/
     |/|<---TL=1000K            |/<---TS            TR=293K|/
     |/|                        |/    SOLID                |/
     |/|       AIR              |/    ALUMINUM             |/
     |/|       PRPS=0           |/    PRPS=100             |/
     |/|                        |/                         |/
     |/|<---Radiative zone 1    |/<---Radiative zone 2     |/
     |_|________________________|__________________________|
        ----------> 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 TS=430K
  ENDDIS
   ************************************************************
   *  GROUP 4.  y-direction grid specification.
NREGY=3;REGEXT(Y,5.0000E-01)
   * Region 1 for a thin solid cell
IREGY=1;GRDPWR(Y,1,0.002,1.0)
   * Region 1 for fluid
IREGY=2;GRDPWR(Y,4,5.0000E-01,-2.0)
   * Region 2 for Solid
IREGY=3;GRDPWR(Y,4,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).
SETPRPS(1,0)
PRESS0=1.0000E+05
   *  GROUP 11. Initialization of fields of variables,
   *            porosities, etc.
   * Initialize Properties Field and block-correction
   * (Meaning of integers used to represent different
   *  properties(PRPS) can be found in the Manual)
   * Define a solid cell for the linear heat transfer
CONPOR(SLD1,-1,CELL,#1,#1,#1,#1,#1,#1)
COVAL(SLD1,PRPS,0.0,100);COVAL(SLD1,TEM1,0.0,1000)
   * Define the aluminum block
CONPOR(SLD,-1,CELL,-#1,-#1,-#3,-#3,-#1,-#1)
COVAL(SLD,PRPS,0.0,100)
FIINIT(TEM1)=10
   *  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,#3,#3,#1,#1,#1,#NREGT)
COVAL(TFIXR,TEM1,1.0,293)
PATCH(@RI001,NORTH,#1,#1,#1,#1,#1,#1,#1,#NREGT)
COVAL(@RI001,TEM1,GRND1,GRND1)
PATCH(@RI002,SOUTH,#1,#1,#3,#3,#1,#1,#1,#nregt)
COVAL(@RI002,TEM1,GRND1,GRND1)
   *  GROUP 15. Termination criteria for sweeps and
   *            outer iterations.
LSWEEP=25;TSTSWP=-1
RESFAC=1E-5
IYMON=NY/2
S2SR=T
 LIBREF  =       115
spedat(set,cvd,radcvd,l,t)
STOP