#$r002
 
   AUTOPLOT USE
   file
   phi 5
 
   cl
   msg S-S Radiation for combined radiation and conduction case
   msg for distance < 0.5 --- AIR
   msg for distance > 0.5 --- Aluminum
   msg Analytical solution shows that TS=418K (distance 0.5)
   da 1 tem1
   col3 1
   le y 418
   le y 293
   le x .5
   ENDUSE
 
   *  GROUP 1.  Run identifiers and other preliminaries.
TEXT(1-D FIXED SURFACE TEMPERATURE : R111
TITLE
  DISPLAY
   The test case considers combined radiation and conduction in the
   domain indicated below.
        ___________________________________________________
      /|                        |/                         |/
      /|<---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=418.K, and the net radiative
   flux at zone 2 =54.97kW/m2
  ENDDIS
   ************************************************************
   *  GROUP 4.  y-direction grid specification.
   *  Total number of REGIONS in Y direction 2
   *  Total number of cells in Y direction 2
   *  Total length in Y direction   1 m
NREGY=2;REGEXT(Y,5.0000E-01)
   * Region 1 for fluid
IREGY=1;GRDPWR(Y,4,5.0000E-01,-2.0)
   * Region 2 for Solid
IREGY=2;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).
   *  Fluid/medium is Air-ideal_gas
SETPRPS(1,0)
PRESS0=1.0000E+05;TEMP0=0
   ************************************************************
   *  GROUP 11. Initialization of fields of variables,
   *            porosities, etc.
INIADD=F
   * Initialize Properties Field and block-correction
   * SOLID boundary condition, name SLD
CONPOR(SLD,-1,CELL,-#1,-#1,-#2,-#2,-#1,-#1)
COVAL(SLD,PRPS,0.0,100)
FIINIT(TEM1)=800
   ************************************************************
   *  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. The surface temperature of this zone
     is calculated in GXS2SR
PATCH(@RI002,SOUTH,#1,#1,#2,#2,#1,#1,#1,#nregt)
COVAL(@RI002,TEM1,GRND1,GRND1)
   ************************************************************
   *  GROUP 15. Termination criteria for sweeps and
   *            outer iterations.
LSWEEP=10;TSTSWP=-1
RESFAC=1E-4
   ************************************************************
   *  GROUP 16. Termination criteria for inner iterations.
S2SR=T
   ************************************************************
   *  GROUP 17. Under-relaxation and related devices.
   ************************************************************
   *  GROUP 23. Variable-by-variable field printout and plot
   *            and/or tabulation of spot-values and residuals.
ITABL=3;NPLT=1
   ************************************************************
   *  GROUP 24. Preparation for continuation runs.
 LIBREF  =       111
 spedat(set,cvd,radcvd,l,t)
STOP