TALK=T;RUN(1,1)
R400
DISPLAY
"RADIANT TRANSFER IN A FURNACE"
Comparison is made of the heat flux predicted by IMMERSOL with
that predicted by classical methods for a problem presented in the
textbook "Heat and Mass Transfer" by A.F. Mills, published by
Irwin in 1995.
A 3m x 3m furnace is in question, long enough for variations in
the length direction to be neglected.The top wall is held at
1400 deg K,the side walls at 1700 deg K, and ths floor at
600 deg K.
All surfaces have an emissivity of 0.5.
The task is to predict the heat flux to the floor, which Mills
calculates as 159 kW/m**2
ENDDIS
PHOTON USE
p
con tem1 y 1 fi;0.001 ; pause
con t3 y 1 fi;0.001 ; pause
con qrx y 1 fi;0.001 ; pause
con qrz y 1 fi;0.001 ; pause
con #ah1 y 1 fi;0.001 ; pause
con #ae1 y 1 fi;0.001 ; pause
con #ap1 y 1 fi;0.001 ; pause
con #rs1 y 1 fi;0.001 ; pause
con #cr1 y 1 fi;0.001 ; pause
con #ah3 y 1 fi;0.001 ; pause
con #ae3 y 1 fi;0.001 ; pause
con #ap3 y 1 fi;0.001 ; pause
con #rs3 y 1 fi;0.001 ; pause
con #cr3 y 1 fi;0.001 ; pause
ENDUSE
************************************************************
Q1 created by VDI menu, Version 3.3, Date 24/07/00
CPVNAM=VDI;SPPNAM=Core
************************************************************
IRUNN = 1 ;LIBREF = 400
*************************************************************
Group 1. Run Title
TEXT(IMMERSOL for heat-treatment furnace
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,30,1,30)
* Set overall domain extent:
* xulast yvlast zwlast
name
XSI= 3.000000E+00; YSI= 1.000000E+00; ZSI= 3.000000E+00
RSET(D,CHAM )
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(145) =WDIS ; NAME(146) =WGAP
NAME(147) =T3 ; NAME(148) =QRZ
NAME(149) =QRX ; NAME(150) =TEM1
* Solved variables list
SOLVE(T3 ,TEM1)
STORE(EMIS,SCAT)
* Stored variables list
STORE(QRX ,QRZ ,WGAP,WDIS)
* Additional solver options
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
TERMS (T3 ,N,N,Y,N,Y,N)
************************************************************
Group 9. Properties
SETPRPS(1, 0)
RHO1 = 1.189000E+00
PRESS0 = 1.000000E+05
CP1 = 1.005000E+03
DISWAL
ENUL = 1.544000E-05 ;ENUT = 0.000000E+00
DVO1DT = 3.410000E-03
PRNDTL(T3 ) = -2.580000E-02 ;PRNDTL(TEM1) = -2.580000E-02
EMISS = 1.000000E-05 ;SCATT = 1.000000E-05
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(WDIS) = 1.000000E-01
fiinit(emis)=emiss;fiinit(scat)=scatt
No PATCHes used for this Group
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
No PATCHes used for this Group
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1000
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(TEM1,LINRLX,0.5)
RELAX(T3,LINRLX,0.5)
************************************************************
Group 18. Limits
VARMIN(T3)=600. ;VARMAX(T3)=1700.
VARMIN(TEM1)=600. ;VARMAX(TEM1)=1700.
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
OUTPUT(WDIS,Y,N,N,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 5 ;IYMON = 1 ;IZMON = 5
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
NOWIPE = T
GVIEW(P,0.000000E+00,-1.000000E+00,0.000000E+00)
GVIEW(UP,0.000000E+00,0.000000E+00,1.000000E+00)
> DOM, SIZE, 3.000000E+00, 1.000000E+00, 3.000000E+00
> DOM, MONIT, 4.500000E-01, 5.000000E-01, 4.500000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> DOM, RELAX, 5.000000E-01
> OBJ1, NAME, B1
> OBJ1, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ1, SIZE, 1.000000E-01, 1.000000E+00, 3.000000E+00
> OBJ1, CLIPART, cube4
> OBJ1, ROTATION24, 1
> OBJ1, TYPE, BLOCKAGE
> OBJ1, MATERIAL, 100
> OBJ1, FIXED_TMP, 0.000000E+00, 1.700000E+03
> OBJ1, EMISSIVITY, 5.000000E-01
> OBJ2, NAME, B2
> OBJ2, POSITION, 2.900000E+00, 0.000000E+00, 0.000000E+00
> OBJ2, SIZE, 1.000000E-01, 1.000000E+00, 3.000000E+00
> OBJ2, CLIPART, cube4
> OBJ2, ROTATION24, 1
> OBJ2, TYPE, BLOCKAGE
> OBJ2, MATERIAL, 100
> OBJ2, FIXED_TMP, 0.000000E+00, 1.700000E+03
> OBJ2, EMISSIVITY, 5.000000E-01
> OBJ3, NAME, B3
> OBJ3, POSITION, 0.000000E+00, 0.000000E+00, 2.900000E+00
> OBJ3, SIZE, 3.000000E+00, 1.000000E+00, 1.000000E-01
> OBJ3, CLIPART, cube4
> OBJ3, ROTATION24, 1
> OBJ3, TYPE, BLOCKAGE
> OBJ3, MATERIAL, 100
> OBJ3, FIXED_TMP, 0.000000E+00, 1.400000E+03
> OBJ3, EMISSIVITY, 5.000000E-01
> OBJ4, NAME, B4
> OBJ4, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ4, SIZE, 3.000000E+00, 1.000000E+00, 1.000000E-01
> OBJ4, CLIPART, cube4
> OBJ4, ROTATION24, 1
> OBJ4, TYPE, BLOCKAGE
> OBJ4, MATERIAL, 100
> OBJ4, FIXED_TMP, 0.000000E+00, 6.000000E+02
> OBJ4, EMISSIVITY, 5.000000E-01
inform7begin
(stored var #ah1 is ahco(tem1))
(stored var #ae1 is aeco(tem1))
(stored var #ap1 is apco(tem1))
(stored var #cr1 is corr(tem1))
(stored var #rs1 is resi(tem1))
(stored var #ah3 is ahco(t3))
(stored var #ae3 is aeco(t3))
(stored var #ap3 is apco(t3))
(stored var #cr3 is corr(t3))
(stored var #rs3 is resi(t3))
inform7end
STOP