IMMERSOL  Non-participating media (1D, T3 only)
   PHI 5
   CL;DA 1 T3;COL3 1
   MSG    T3 profile
TEXT(IMMERSOL 1D transparent medium     :201

        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