#$r002
 
   PHOTON USE
   AUTOPLOT
   file
   PHI 5
 
   cl
   msg Temperature profile
   msg Analytical solution shows that
   msg TS1=1000K (distance 0.5) and TS2=418K (distance 1.0)
   da 1 tem1
   col3 1
   le x .5
   le x 1
   le y 293
   le y 418
   le y 1000
   ENDUSE
 
   *  GROUP 1.  Run identifiers and other preliminaries.
TEXT(1-D SOLID-FLUID-SOLID SET-UP: R116
TITLE
  DISPLAY
   The test case considers combined radiation and conduction
   in the domain indicated below.
    ___________________________________________________________
   /|                  /|                   |/                |/
   /|-->qflx           /|<---TS1            |/<---TS2    TR-->|/
   /|                  /|                   |/                |/
   /|    SOLID 1       /|   FLUID           |/  SOLID 2       |/
   /|    ALUMINIUM     /|   AIR             |/  ALUMINIUM     |/
   /|    PRPS=100      /|   PRPS=0          |/  PRPS=100      |/
   /|                  /|                   |/                |/
   /|                  /|<--thermal zone 1  |/<-thermal zone 2|/
   /|__________________/|___________________|/________________|/
    --------> Y
     qflx = 55kW/m2; TR=293K; TS1 and TS2 to be computed.
     length of domain = 1.5 m, height of domain = 1.0m
 
   The PHOENICS solution can be compared with the analytical
   solution which shows that TS1=1000.K, TS2=418.K and
   the net radiative flux at zone 2 =54.97kW/m2
  ENDDIS
 
   *  User declared Real variables
   * NYG1   - the Y solid cell at zone 1 used in Group 12
   * NYG2   - the Y gas cell at zone 2 used in Group 12
   * NYG1 and NYG2 set in Group 4
INTEGER(NYG1,NYG2)
   *  GROUP 4.  y-direction grid specification.
NREGY=3;REGEXT(Y,5.0000E-01)
   * Region 1 for solid
IREGY=1;GRDPWR(Y,-4,5.0000E-01,1.9)
   * Region 2 for fluid
NYG1=4
IREGY=2;GRDPWR(Y,-4,5.0000E-01,1.9)
NYG2=8
   * Region 3 for solid
IREGY=3;GRDPWR(Y,4,5.0000E-01,1.9)
   *  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; TEMP0=0
   *  GROUP 11. Initialization of fields of variables,
   *            porosities, etc.
   * Initialize Properties Field
   * SOLID boundary condition, name SLD
CONPOR(SLD,-1,CELL,-#1,-#1,-#1,-#1,-#1,-#1)
COVAL(SLD,PRPS,0.0,100)
   * SOLID boundary condition, name SLD2
CONPOR(SLD2,-1,CELL,-#1,-#1,-#3,-#3,-#1,-#1)
COVAL(SLD2,PRPS,0.0,100)
FIINIT(TEM1)=1000
   *  GROUP 13. Boundary conditions and special sources
   * boundary condition, name TFLXL to set the heat flux
     for the conductive boundary condition
PATCH(TFLXL,SOUTH,#1,#1,#1,#1,#1,#1,#1,#NREGT)
COVAL(TFLXL,TEM1,FIXFLU,55000.)
   * 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)
   * 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. The surface temperature of this zone
     is calculated in GXS2SR
PATCH(@RI001,NORTH,#1,#1,#1,#1,#1,#1,#1,#NREGT)
COVAL(@RI001,TEM1,GRND1,GRND1)
  * Radiation zone 2. The surface temperature of this zone
    is calculated in GXS2SR
PATCH(@RI002,SOUTH,#1,#1,#3,#3,#1,#1,#1,#nregt)
COVAL(@RI002,TEM1,GRND1,GRND1)
LSWEEP=50;TSTSWP=-1
RELAX(TEM1,FALSDT,5000.)
S2SR=T
spedat(set,cvd,radcvd,l,t)
RESFAC=1E-6