TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2025, Date 27/01/25
CPVNAM=VDI; SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
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
************************************************************
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)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(141)=T3SV ;NAME(143)=SCAT
NAME(144)=EMIS ;NAME(145)=WDIS
NAME(146)=WGAP ;NAME(147)=T3
NAME(148)=QRZ ;NAME(149)=QRX
NAME(150) =TEM1
* Solved variables list
SOLVE(T3,TEM1)
* Stored variables list
STORE(QRX,QRZ,WGAP,WDIS,EMIS,SCAT,T3SV)
* Additional solver options
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
************************************************************
Echo save-block settings for Group 7
SAVE7BEGIN
(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))
(make1 t31s is 0)
(store1 t31s is sum(volu*t3sv) with imat<100)
(print of radiative_source is t31s)
SAVE7END
************************************************************
Group 8. Terms & Devices
TERMS (T3 ,N,N,Y,N,Y,N)
************************************************************
Group 9. Properties
PRESS0 =1.0E+05
* Domain material index is 0 signifying:
* Air at 20 deg C, 1 atm, treated as incompressible
SETPRPS(1, 0)
DISWAL
ENUT =0.
DVO1DT =3.41E-03
PRNDTL(T3)=-0.0258 ;PRNDTL(TEM1)=-0.0258
EMISS =1.0E-05 ;SCATT =1.0E-05
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(SCAT)=1.0E-05 ;FIINIT(EMIS)=1.0E-05
FIINIT(WDIS)=0.1
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.0E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(T3 ,LINRLX,0.5 )
RELAX(TEM1,LINRLX,0.5 )
************************************************************
Group 18. Limits
VARMAX(T3)=1700. ;VARMIN(T3)=600.
VARMAX(TEM1)=1700. ;VARMIN(TEM1)=600.
************************************************************
Group 19. EARTH Calls To GROUND Station
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
************************************************************
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
GVIEW(P,0.,-1.,0.)
GVIEW(UP,0.,0.,1.)
GVIEW(DEPTH,1.17391E+05)
GVIEW(VDIS,2.092272)
GVIEW(CENTRE,1.5,0.5,1.5)
> DOM, SIZE, 3.000000E+00, 1.000000E+00, 3.000000E+00
> DOM, MONIT, 4.499999E-01, 5.000000E-01, 4.499999E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID, RSET_X_1, 1, 1.000000E+00
> GRID, RSET_X_2, 28, 1.000000E+00
> GRID, RSET_X_3, 1, 1.000000E+00
> GRID, RSET_Y_1, 1, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> GRID, RSET_Z_2, 28, 1.000000E+00
> GRID, RSET_Z_3, 1, 1.000000E+00
> OBJ, NAME, B1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E-01, 1.000000E+00, 3.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube4
> OBJ, TYPE, BLOCKAGE
> OBJ, COLOR-MODE, FROM_FIL
> OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c
> OBJ, FIXED_TMP, 0. ,1700.
> OBJ, EMISSIVITY, 0.5
> OBJ, NAME, B2
> OBJ, POSITION, 2.900000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E-01, 1.000000E+00, 3.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube4
> OBJ, TYPE, BLOCKAGE
> OBJ, COLOR-MODE, FROM_FIL
> OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c
> OBJ, FIXED_TMP, 0. ,1700.
> OBJ, EMISSIVITY, 0.5
> OBJ, NAME, B3
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.900000E+00
> OBJ, SIZE, 3.000000E+00, 1.000000E+00, 1.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube4
> OBJ, TYPE, BLOCKAGE
> OBJ, COLOR-MODE, FROM_FIL
> OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c
> OBJ, FIXED_TMP, 0. ,1400.
> OBJ, EMISSIVITY, 0.5
> OBJ, NAME, B4
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 3.000000E+00, 1.000000E+00, 1.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube4
> OBJ, TYPE, BLOCKAGE
> OBJ, COLOR-MODE, FROM_FIL
> OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c
> OBJ, FIXED_TMP, 0. ,600.
> OBJ, EMISSIVITY, 0.5
STOP