TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2012, Date 19/10/12
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
A direct-fired furnace for heat treating solid materials is a
rectangular 3D box in shape. The side walls and roof are
naturally cooled by ambient air with an estimated heat transfer
coefficient. The have an emmitance of 0.9. The furnace floor is
refractory and covered by three beams of load made from brick
material with an emmitance of 0.9.
The preheated methane gas premixed with air enters the furnace
with an angle to the floor through closely spaced ports at the
side wall. It ignites and steady kinetically-controled combustion
takes place producing hot non-scattering furnace gases with an
absorption coefficient equal to 0.8. The exhaust opening is at
the opposite side wall.
The task is to calculate the operating temperatures of the
load and furnace gases along with all related field
distributions.
ENDDIS
VRV USE
SLICE Z
* View and up directions
VIEW 0.000000E+00 0.000000E+00 -1.000000E+00
UP 0.000000E+00 1.000000E+00 0.000000E+00
* View centre
VIEW CENTRE 5.000000E-01 1.600000E-01 5.000000E-01
* View size
VIEW SIZE 6.524845E-01
* View perspective
VIEW DEPTH 3.000000E+00;VIEW TILT 0.8
VARIABLE TEMPERATURE
VECTOR OFF
CONTOUR ON
SURFACE OFF
PAUSE
VARIABLE T3
PAUSE
VARIABLE QRY
PAUSE
VARIABLE QRX
ENDUSE
PHOTON USE
p; ; ; ; ;
msg Temperature contours (Tmp1)
con tmp1 z 3 fi;0.001 ; pause;cl
msg Radiation temperature contours (T3)
con t3 z 3 fi;0.001 ; pause;cl
msg x-direction heat flux (qrx) contours
con qrx z 3 fi;0.001 ; pause;cl
msg y-direction heat flux (qry) contours
con qry z 3 fi;0.001 ; pause;cl
MSG Solution coefficients for H1
con #an1 z 3 fi;0.001 ; pause
con #ae1 z 3 fi;0.001 ; pause
con #ap1 z 3 fi;0.001 ; pause
con #rs1 z 3 fi;0.001 ; pause
con #cr1 z 3 fi;0.001 ; pause;cl
MSG Solution coefficients for T3
con #an3 z 3 fi;0.001 ; pause
con #ae3 z 3 fi;0.001 ; pause
con #ap3 z 3 fi;0.001 ; pause
con #rs3 z 3 fi;0.001 ; pause
con #cr3 z 3 fi;0.001 ; pause
ENDUSE
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(Heat treatment in a direct-fired furnace)
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,20,10,5)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(132)=CP1 ;NAME(133)=SCAT
NAME(134)=EMIS ;NAME(135)=EPKE
NAME(136)=ENUT ;NAME(137)=GEN1
NAME(138)=WDIS ;NAME(139)=WGAP
NAME(140)=PROD ;NAME(141)=OXID
NAME(142)=PRPS ;NAME(143)=QRZ
NAME(144)=QRY ;NAME(145)=QRX
NAME(146)=DEN1 ;NAME(147)=TMP1
NAME(148)=FUEL ;NAME(149)=T3
NAME(150) =MIXF
* Solved variables list
SOLVE(P1,U1,V1,W1,H1,FUEL,T3,MIXF)
* Stored variables list
STORE(TMP1,DEN1,QRX,QRY,QRZ,PRPS,OXID,PROD)
STORE(WGAP,WDIS,GEN1,ENUT,EPKE,EMIS,SCAT,CP1)
TURMOD(KEMODL)
************************************************************
Group 8. Terms & Devices
TERMS (H1 ,N,Y,Y,N,Y,N)
TERMS (T3 ,N,N,Y,N,Y,N)
NEWRH1 = T
************************************************************
Group 9. Properties
PRESS0 =1.0E+05
* Domain material index is 30 signifying:
* Ideal gas Law for mixture of 3 gases (ENUL=3.6E-5/GRHO)
SETPRPS(1, 30)
TMP1 = GRND8
TMP1A =1500. ;TMP1B =1500.
TMP1C =1500.
TMP2A =0.054825 ;TMP2B =4.9E+07
TMP2C =0.
DISWAL
DRH1DP =3.41E-03
DVO1DT =3.33E-03
PRNDTL(H1)=-0.0263 ;PRNDTL(T3)=-0.0263
PRT(EP)=1.314
EMISS =0.8 ;SCATT =0.
RADIA =0.8 ;RADIB =0.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(U1)=2. ;FIINIT(KE)=0.45
FIINIT(EP)=0.57404 ;FIINIT(H1)=1.1595E+06
FIINIT(EMIS)=0.8 ;FIINIT(EPKE)=1.
FIINIT(ENUT)=1.001E-10 ;FIINIT(GEN1)=1.001E-10
FIINIT(WDIS)=0.032 ;FIINIT(PROD)=1.001E-10
FIINIT(OXID)=1.001E-10 ;FIINIT(PRPS)=30.
FIINIT(DEN1)=0.3 ;FIINIT(TMP1)=300.
FIINIT(FUEL)=0.1 ;FIINIT(MIXF)=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
PATCH(CHSO, PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(CHSO, FUEL, GRND9, GRND9)
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1000
RESFAC =1.0E-06
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,0.9 )
RELAX(KE ,LINRLX,0.5 )
RELAX(EP ,LINRLX,0.5 )
RELAX(H1 ,LINRLX,0.5 )
RELAX(T3 ,LINRLX,0.5 )
RELAX(DEN1,LINRLX,0.9 )
RELAX(MIXF,LINRLX,0.9 )
KELIN = 3
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-5.0E-06
VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=300.
************************************************************
Group 19. EARTH Calls To GROUND Station
GENK = T
CONWIZ = T
RELAX(DEN1,LINRLX,0.1)
CHSOA =0.054825 ;CHSOB =1.
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
OUTPUT(W1 ,Y,Y,Y,N,Y,Y)
OUTPUT(WDIS,Y,N,N,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 10 ;IYMON = 1 ;IZMON = 3
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.,0.,-1.)
GVIEW(UP,0.,1.,0.)
> DOM, SIZE, 1.000000E+00, 3.200000E-01, 1.000000E+00
> DOM, MONIT, 4.750000E-01, 1.600000E-02, 5.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID, RSET_X_1, 20, 1.000000E+00
> GRID, RSET_Y_1, 10, 1.000000E+00
> GRID, RSET_Z_1, 5, 1.000000E+00
> OBJ, NAME, IMSWL1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube11
> OBJ, TYPE, PLATE
> OBJ, EMISSIVITY, 1.
> OBJ, NAME, IMSWL2
> OBJ, POSITION, 0.000000E+00, 3.200000E-01, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube13
> OBJ, TYPE, PLATE
> OBJ, SURF_HEAT, 1. ,0.
> OBJ, PARA_HEAT, 1. ,1. ,0.
> OBJ, T_EXT, 300.
> OBJ, EMISSIVITY, 0.9
> OBJ, NAME, B5
> OBJ, POSITION, 2.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E-01, 9.600000E-02, 6.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube
> OBJ, GRID, N,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 102,BRICK at 20 deg c
> OBJ, EMISSIVITY, 1.
> OBJ, NAME, B6
> OBJ, POSITION, 1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 3.200000E-01, 1.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube13
> OBJ, TYPE, PLATE
> OBJ, SURF_HEAT, 0. ,0.
> OBJ, PARA_HEAT, 1. ,1. ,0.
> OBJ, T_EXT, 300.
> OBJ, EMISSIVITY, 0.9
> OBJ, NAME, B7
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 3.200000E-01, 1.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube13
> OBJ, TYPE, PLATE
> OBJ, SURF_HEAT, 0. ,0.
> OBJ, PARA_HEAT, 1. ,1. ,0.
> OBJ, T_EXT, 300.
> OBJ, EMISSIVITY, 0.9
> OBJ, NAME, OUTLET
> OBJ, POSITION, 1.000000E+00, 2.400000E-01, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 8.000000E-02, 1.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube12t
> OBJ, GRID, N,N,N
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.
> OBJ, TEMPERATURE, 0.
> OBJ, COEFFICIENT, 1000.
> OBJ, TURBULENCE, SAME , SAME
> OBJ, T_EXT, 3.000000E+02
> OBJ, NAME, INLET
> OBJ, POSITION, 0.000000E+00, 1.200000E-01, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 8.000000E-02, 1.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube3t
> OBJ, GRID, N,N,N
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 3.606
> OBJ, VELOCITY, 10. ,-10. ,0.
> OBJ, TEMPERATURE, 700.
> OBJ, INLET_FUEL, 0.1
> OBJ, INLET_MIXF, 0.1
> OBJ, KE_IN, 1.000000E-02
> OBJ, EP_IN, 9.400417E-02
> OBJ, T_EXT, 3.000000E+02
> OBJ, NAME, B8
> OBJ, POSITION, 5.000001E-01, 0.000000E+00, 4.000000E-01
> OBJ, SIZE, 1.000000E-01, 9.600000E-02, 6.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube14
> OBJ, GRID, N,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 102,BRICK at 20 deg c
> OBJ, EMISSIVITY, 1.
> OBJ, NAME, B9
> OBJ, POSITION, 7.999998E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E-01, 9.600000E-02, 6.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube14
> OBJ, GRID, N,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 102,BRICK at 20 deg c
> OBJ, EMISSIVITY, 1.
> OBJ, NAME, B19
> OBJ, POSITION, 0.000000E+00, 1.200000E-01, 2.800000E-01
> OBJ, SIZE, 1.999998E-02, 7.200000E-02, 1.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, port_vca
> OBJ, ROTATION24, 9
> OBJ, TYPE, CELLTYPE
> OBJ, NAME, B11
> OBJ, POSITION, 0.000000E+00, 1.200000E-01, 0.000000E+00
> OBJ, SIZE, 1.999998E-02, 7.200000E-02, 1.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, port_vca
> OBJ, ROTATION24, 9
> OBJ, TYPE, CELLTYPE
> OBJ, NAME, B12
> OBJ, POSITION, 0.000000E+00, 1.200000E-01, 6.000000E-01
> OBJ, SIZE, 1.999998E-02, 7.200000E-02, 1.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, port_vca
> OBJ, ROTATION24, 9
> OBJ, TYPE, CELLTYPE
> OBJ, NAME, B13
> OBJ, POSITION, 0.000000E+00, 1.200000E-01, 1.400000E-01
> OBJ, SIZE, 1.999998E-02, 7.200000E-02, 1.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, port_vca
> OBJ, ROTATION24, 9
> OBJ, TYPE, CELLTYPE
> OBJ, NAME, B14
> OBJ, POSITION, 0.000000E+00, 1.200000E-01, 7.799999E-01
> OBJ, SIZE, 1.999998E-02, 7.200000E-02, 1.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, port_vca
> OBJ, ROTATION24, 9
> OBJ, TYPE, CELLTYPE
> OBJ, NAME, B15
> OBJ, POSITION, 0.000000E+00, 1.200000E-01, 4.400001E-01
> OBJ, SIZE, 1.999998E-02, 7.200000E-02, 1.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, port_vca
> OBJ, ROTATION24, 9
> OBJ, TYPE, CELLTYPE
> OBJ, NAME, B16
> OBJ, POSITION, 9.880000E-01, 2.400000E-01, 2.000000E-01
> OBJ, SIZE, 1.200000E-02, 8.000000E-02, 6.000000E-01
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, port_vca
> OBJ, ROTATION24, 21
> OBJ, TYPE, CELLTYPE
> OBJ, NAME, B17
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ, SIZE, 1.000000E+00, 3.200000E-01, 0.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube13
> OBJ, TYPE, PLATE
> OBJ, SURF_HEAT, 0. ,0.
> OBJ, PARA_HEAT, 1. ,1. ,0.
> OBJ, T_EXT, 300.
> OBJ, EMISSIVITY, 0.9
> OBJ, NAME, B18
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 3.200000E-01, 0.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube13
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, SURF_HEAT, 0. ,0.
> OBJ, PARA_HEAT, 1. ,1. ,0.
> OBJ, T_EXT, 300.
> OBJ, EMISSIVITY, 0.9
DISTIL=T
STOP
inform7begin
(stored var #an1 is anco(H1))
(stored var #ae1 is aeco(H1))
(stored var #ap1 is apco(H1))
(stored var #cr1 is corr(H1))
(stored var #rs1 is resi(H1))
(stored var #an3 is anco(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