#immersol
#declare
   IMMERSOL  Non-participating media (1D, T3 only)
  **************************************************************
   PHOTON USE
   AUTOPLOT
   FILE
   PHI 5
 
   CL;DA 1 T3;COL3 1
   MSG    T3 profile
   ENDUSE
  **************************************************************
TEXT(IMMERSOL 1D transparent medium     :201
#cls
TITLE
  DISPLAY 
  

  
        solids have artificially high thermal conductivity
      hotter solid       temperature = Th    emissivity = eph
   // 0.01 m thick////////////////////////////////////////////////
 
 
      non-participating medium, 1 m thick
 
 
 
      colder solid       temperature = Tc    emissivity = epc
   // 0.01 m thick ////////////////////////////////////////////// 
   
   
   Exact solution:  Qr = (Eh - Ec)/(1/eph + 1/epc - 1)
              where Eh and Ec are the emissive powers ( E = S*T**4 )
                    S is the Stefan-Boltzman constant.
   Values:  Th = 800K, Tc = 400K, eph = 0.3, epc = 0.1
   Confirmation or modification of these settings will be requested.
 
   The macro r193 is loaded to set emissivity via a 3D store, and
   to activate line-printer plots.
   
   It is necessary to store PRPS so that the LTLS solver knows which
   cells contain solids, and therefore computes the correct WGAP.
   
  ENDDIS
#pause
STORE(PRPS)
  ** start of optional data input by the user **
#direct
#temps
WSL1 = 0.01;  WSL2= 0.01;  WGAP= 1.0;
QRAD= SIGMA*(THOT**4 - TCLD**4)
MESG( Expected radiative heat flux (W/m**2):
QRAD = QRAD/(1./EMISH + 1./EMISC - 1.)
QRAD
  ** end of optional data input by the user **
#geom1
#radflux
INIADD= F;  FIINIT(T3)= 0.5*( THOT + TCLD) ; FIINIT(PRPS) = 0.0
  !!!! the practice of ascribing emissivities to materials by way 
      of SPEDAT is no longer recommended. 
      Macro r193 illustrates the better practice of setting the
      values directly by way of the 3D-stored EMIS variable !!!!
      
      Consequently the value ascribed to PRPS is no longer 
      significant, so long as it is greater than or equal to
      SOLPRP (usually 100.0).
INIT(SOL1,PRPS,0.0,111.); INIT(SOL2,PRPS,0.0,111.)
  
  SPEDAT(SET,EMISSIVITY,111,R,:EMISH:)
  SPEDAT(SET,EMISSIVITY,112,R,:EMISC:)
COVAL(HOT,T3,FIXVAL,THOT); COVAL(COLD,T3,FIXVAL,TCLD)
LSWEEP= 50; SELREF=F; RESREF(T3)= 1.E-3
IXMON=NX/2+1;IYMON=NY/2+1;IZMON=NZ/2+1;
NXPRIN=1; NYPRIN=1; NZPRIN=1
OUTPUT(LTLS,N,N,N,N,N,N); OUTPUT(WDIS,N,N,N,N,N,N)
load($r193)

VARMAX(T3)=THOT;VARMIN(T3)=TCLD

stop