**** PIL fragment for IMMERSOL method
REAL(SIGMA); SIGMA = 5.6697E-8
STORE(EMIS,SCAT)
EMISS=0.0;SCATT=0.0;FIINIT(EMIS)=EMISS;FIINIT(SCAT)=SCATT
SOLVE(T3);RELAX(T3,LINRLX,0.5);STORE(WGAP)
DISWAL
TERMS(T3,N,N,Y,N,Y,N)
char(temps,direct,declare);temps=$r198;direct=$r199;declare=$r200
char(geom1,geom2,radflux);geom1=$r195;geom2=$r196;radflux=$r197
char(geom3);geom3=$r194
MESG(IMMERSOL settings have been loaded
** LOAD(x 1) from the x Input Library
** LOAD(x 1) from the x Input Library
** declarations needed for many immersol examples **
REAL(QRAD,WSL1,WSL2,WGAP,THOT,TCLD,EMISH,EMISC,EMISB)
REAL(OPTHI,KROSS,TGHT,TGCL,KGAS,NN1,NN2,AA1,AA2)
REAL(HHOT,HCLD,CP111,CP112)
REAL(SIZX,SIZZ,SIZY,UIN)
REAL(QCNV,AEMISH,AEMISC,NNH,NNC,TGH,TGC)
REAL(LENG1,LENG2,GAP,WPLT)
REAL(EMISM1,EMISM2,EMISW1,EMISW2,REYNO)
REAL(SPH0,TOXY,TFUEL,TWALL,HAIRIN,CPAIR,CPPR,DTF)
REAL(CPFU,WAIR,WFU,WPR,GRHO,STOIC,FSTOI,HFU,TKEI1,EMISI1,FINF,HINF)
** LOAD(x200) from the x Input Library
** LOAD(x200) from the x Input Library
IMMERSOL Emitting and absorbing gray media
(1D, TEM1)
TEXT(IMMERSOL 1D Radiative equilibrium :202
TITLE
**************************************************************
PHOTON USE
AUTOPLOT
FILE
PHI 5
CL;DA 1 TEM1;DA 1 T3;COL3 1;BLB4 2;REDR
MSG Temperature (curve) and T3 (crosses) profiles
ENDUSE
**************************************************************
#cls
DISPLAY
The problem is as for case r202; but
1. the solid walls at the boundaries are both thicker and divided
into five cells each;
2. The heat conduction due to true (ie not T3) temperature
gradients, which was switched off in r202, is swiched back on
again;
3. In-Form "stored" statements are used for eliciting print-out
of coefficients, which is useful for checking correctness.
4. Also printed is the variable #3-1, which is the value of
t3-tem1. This should ideally be zero in the solid regions;
but differences of a few degrees will be noted.
ENDDIS
** sequence used in several immersol examples
INTEGER(IDIR,I1,I2,J1,J2,K1,K3)
MESG( The problem may be stated and solved in any direction
MESG( ENTER 1 for X; 2 (or RETURN) for Y, 3 for Z.
READVDU(IDIR,INT,2)
MESG( ENTER Temperature of hot solid (default 800K):
READVDU(THOT,REAL,800.)
MESG( ENTER Temperature of cold solid (default 400K):
READVDU(TCLD,REAL,400.)
MESG( ENTER Surface emissivity of hot solid (default 0.3):
READVDU(EMISH,REAL,0.3)
MESG( ENTER Surface emissivity of cold solid (default 0.1):
READVDU(EMISC,REAL,0.1)
** LOAD(x198) from the x Input Library
** LOAD(x198) from the x Input Library
MESG( ENTER optical thickness Kr*Wg (default 1.):
READVDU(OPTHI,REAL,1.0)
settings of wall thicknesses and the gap between them
WSL1 = 0.01; WSL2= 0.01; WGAP= 1.0;
conductivity of medium swet to zero
KGAS= 0.0
QRAD first used to store EH-EC (see above)
QRAD= SIGMA*(THOT**4 - TCLD**4)
KROSS deduced from above-set optical thickness; scattering
coefficient set to zero; then medium emissivity deduced
KROSS= OPTHI/WGAP; SCATT= 0.0; EMISS= KROSS-SCATT
MESG( Expected radiative heat flux (W/m**2):
QRAD now computed, based on above settings.
Note that KGAS has been set to zero, so that NN1 and NN2 = 0
NN1 = KGAS*KROSS/4./SIGMA/THOT**3
NN2 = KGAS*KROSS/4./SIGMA/TCLD**3
QRAD= QRAD + KGAS/WGAP*(THOT-TCLD)*0.75*OPTHI
AA1 = (1./EMISH-0.5)/(1.+0.75*NN1)
AA2 = (1./EMISC-0.5)/(1.+0.75*NN2)
QRAD= QRAD/(0.75*OPTHI + AA1 + AA2)
QRAD
MESG( Expected gas temperatures near walls are:
TGCL = ((SIGMA*TCLD**4 + QRAD*AA2)/SIGMA)**0.25
TGHT = ((SIGMA*THOT**4 - QRAD*AA1)/SIGMA)**0.25
TGCL
TGHT
#pause
SOLVE(TEM1)
STORE(PRPS)
** geom1 **
IF(IDIR.EQ.1) THEN
NREGX= 3
IREGX= 1; GRDPWR(X,5,0.5*wgap,1.0)
IREGX= 2; GRDPWR(X,10,WGAP,1.0)
IREGX= 3; GRDPWR(X,5,0.5*wgap,1.0)
PATCH(SOL1,INIVAL, #1,#1,1,NY,1,NZ,1,LSTEP)
PATCH(SOL2,INIVAL,#3,#3,1,NY,1,NZ,1,LSTEP)
PATCH(HOT,CELL,1, 1,1,NY,1,NZ,1,LSTEP)
PATCH(COLD,CELL,nx,nx,1,NY,1,NZ,1,LSTEP)
ENDIF
IF(IDIR.EQ.2) THEN
NREGY= 3; IREGY= 1; GRDPWR(Y,5,0.5*wgap,1.0)
IREGY= 2; GRDPWR(Y,10,WGAP,1.0);
IREGY= 3; GRDPWR(Y,5,0.5*wgap,1.0)
PATCH(SOL1,INIVAL,1,NX, #1,#1,1,NZ,1,LSTEP)
PATCH(SOL2,INIVAL,1,NX,#3,#3,1,NZ,1,LSTEP)
PATCH( HOT,CELL,1,NX, 1,1,1,NZ,1,LSTEP)
PATCH(COLD,CELL,1,NX,ny,ny,1,NZ,1,LSTEP)
ENDIF
IF(IDIR.EQ.3) THEN
NREGZ= 3; IREGZ= 1; GRDPWR(Z,5,0.5*wgap,1.0)
IREGZ= 2; GRDPWR(Z,10,WGAP,1.0)
IREGZ= 3; GRDPWR(Z,5,0.5*wgap,1.0)
PATCH(SOL1,INIVAL,1,NX,1,NY, #1, #1,1,LSTEP)
PATCH(SOL2,INIVAL,1,NX,1,NY,#3,#3,1,LSTEP)
PATCH( HOT,CELL,1,NX,1,NY,1,1,1,LSTEP)
PATCH(COLD,CELL,1,NX,1,NY,nz,nz,1,LSTEP)
ENDIF
** LOAD(x195) from the x Input Library
** LOAD(x195) from the x Input Library
** radflux **
if(nx.gt.1) then
store(qrx)
endif
if(ny.gt.1) then
store(qry)
endif
if(nz.gt.1) then
store(qrz)
endif
** LOAD(x197) from the x Input Library
** LOAD(x197) from the x Input Library
TERMS(TEM1,N,N,y,N,Y,N); PRNDTL(TEM1)=1.e10
INIADD= F; FIINIT(TEM1)= (THOT+TCLD)/2.; FIINIT(PRPS)= 0.
FIINIT(T3)=FIINIT(TEM1)
COVAL(HOT,TEM1,FIXVAL,THOT); COVAL(COLD,TEM1,FIXVAL,TCLD)
COVAL(HOT,TEM1,1.e5,THOT); COVAL(COLD,TEM1,1.e5,TCLD)
INIT(SOL1,PRPS,0.0,111.); INIT(SOL2,PRPS,0.0,112.)
LSWEEP= 500; 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 macro for storing and setting emissivity and scattering
coefficient, using EMISS, SCATT, EMISH and EMISC
use of 3D store for emissivity,
over-writing what may be supplied by SPEDAT
STORE(EMIS,SCAT)
FIINIT(EMIS)=EMISS;FIINIT(SCAT)=SCATT
IF(IDIR.EQ.1) THEN
PATCH( HOTINI,INIVAL, #1, #1,1,NY,1,NZ,1,LSTEP)
PATCH(COLDINI,INIVAL,#3,#3,1,NY,1,NZ,1,LSTEP)
line-printer output
PATCH(PROFILES,PROFIL,1,NX,1,1,1,1,1,1)
COVAL(PROFILES,QRX,0.0,0.0)
ENDIF
IF(IDIR.EQ.2) THEN
PATCH( HOTINI,INIVAL, 1, NX,#1,#1,1,NZ,1,LSTEP)
PATCH(COLDINI,INIVAL,1,NX,#3,#3,1,NZ,1,LSTEP)
line-printer output
PATCH(PROFILES,PROFIL,1,1,1,NY,1,1,1,1)
COVAL(PROFILES,QRY,0.0,0.0)
ENDIF
IF(IDIR.EQ.3) THEN
PATCH( HOTINI,INIVAL, 1, NX,1,NY,#1,#1,1,LSTEP)
PATCH(COLDINI,INIVAL,1,NX,1,NY,#3,#3,1,LSTEP)
line-printer output
PATCH(PROFILES,PROFIL,1,1,1,1,1,NZ,1,1)
COVAL(PROFILES,QRZ,0.0,0.0)
ENDIF
setting absorptivity/emissivity of solids
COVAL(HOTINI,EMIS,0.0,EMISH)
COVAL(COLDINI,EMIS,0.0,EMISC)
COVAL(PROFILES,T3,0.0,0.0)
IF(LIBREF.NE.201.AND.LIBREF.NE.204) THEN
COVAL(PROFILES,TEM1,0.0,0.0)
ENDIF
** LOAD(x193) from the x Input Library
** LOAD(x193) from the x Input Library
VARMAX(T3)=THOT;VARMIN(T3)=TCLD
VARMAX(TEM1)=THOT;VARMIN(TEM1)=TCLD
inform7begin
(stored var #3-1 is t3-tem1)
(stored var #rs1 is resi(tem1))
(stored var #rs3 is resi(t3))
(stored var #ap1 is apco(tem1))
(stored var #ap3 is apco(t3))
if(idir.eq.1) then
(stored var #ae1 is aeco(tem1))
(stored var #ae3 is aeco(t3))
endif
if(idir.eq.2) then
(stored var #an1 is anco(tem1))
(stored var #an3 is anco(t3))
endif
if(idir.eq.3) then
(stored var #ah1 is ahco(tem1))
(stored var #ah3 is ahco(t3))
endif
inform7end