TALK=T;RUN( 1, 1)
DISPLAY
LID-DRIVEN FLOW IN A RADIATING CAVITY
An IMMERSOL application is presented aimed at the demonstration
of the method for the situation when both thermal radiation
and convection affect the heat transfer.
This simulation analyses the flow in a chamber, the floor of
which moves with a constant velocity, 1m/s. The stationary
walls are adiabatic, the moving wall is at a 300K.
A radiant heater panel is placed at the top of the chamber.
It consists of a row of cylindrical electrical heating
elements 1 cm in diameter, 150 cm long, spaced at a 3 cm pitch,
and backed by a well insulated wall. The panel has dimensions
30x150 cm and is located 30 cm above a workpiece floor, which
is also 30x150 cm. The heater elements are rated at 5 kW each.
The emittances of the elements and back wall are 0.9 and 0.8,
respectively; the workpiece floor is assumed black at 300K.
The black side walls are supposed to be made from refractory
(adiabatic) bricks.
The task is to calculate the operating temperatures of the heater
elements, side wall temperatures from the velocity and temperature
fields and radiative heat fluxes distributions.
ENDDIS
PHOTON USE
p; ; ; ; ;
msg Temperature contours (TEM1)
con tEM1 z 1 fi;0.001 ; pause;cl
msg Radiation temperature contours (T3)
con t3 z 1 fi;0.001 ; pause;cl
msg x-direction heat flux (qrx) contours
con qrx z 1 fi;0.001 ; pause;cl
msg y-direction heat flux (qry) contours
con qry z 1 fi;0.001 ; pause;cl
MSG Solution coefficients for Tem1
con #an1 z 1 fi;0.001 ; pause
con #ae1 z 1 fi;0.001 ; pause
con #ap1 z 1 fi;0.001 ; pause
con #rs1 z 1 fi;0.001 ; pause
con #cr1 z 1 fi;0.001 ; pause;cl
MSG Solution coefficients for T3
con #an3 z 1 fi;0.001 ; pause
con #ae3 z 1 fi;0.001 ; pause
con #ap3 z 1 fi;0.001 ; pause
con #rs3 z 1 fi;0.001 ; pause
con #cr3 z 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 = 197
************************************************************
Group 1. Run Title
TEXT(Lid-driven flow in a radiating cavity )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,32,32,1)
* Set overall domain extent:
* xulast yvlast zwlast
name
XSI= 3.200000E-01; YSI= 3.200000E-01; ZSI= 1.500000E+00
RSET(D,CHAM )
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(143) =ENUL ; NAME(144) =QRY
NAME(145) =QRX ; NAME(147) =WDIS
NAME(148) =WGAP ; NAME(149) =T3
NAME(150) =TEM1
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,T3 ,TEM1)
* Stored variables list
STORE(WGAP,WDIS,QRX ,QRY ,ENUL)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
TERMS (T3 ,N,N,Y,N,Y,N)
NEWENL = T
************************************************************
Group 9. Properties
SETPRPS(1, 0)
RHO1 = 1.189000E+00
PRESS0 = 1.000000E+05
CP1 = 1.005000E+03
DISWAL
ENUL = GRND6 ;ENUT = 0.000000E+00
ENULA = 1.460000E-06 ;ENULB = 1.100000E+02
ENULC = 0.000000E+00
DVO1DT = 3.410000E-03
PRNDTL(T3 ) = -2.580000E-02 ;PRNDTL(TEM1) = -2.580000E-02
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(WDIS) = 3.200000E-02
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 = 1200
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
************************************************************
Group 18. Limits
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
VARMAX(T3 ) = 2.000000E+03 ;VARMIN(T3 ) = 3.000000E+02
VARMAX(TEM1) = 2.000000E+03 ;VARMIN(TEM1) = 3.000000E+02
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
RADIA = 1.000000E-09 ;RADIB = 1.000000E-09
CONWIZ = T
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
OUTPUT(P1 ,Y,Y,Y,N,Y,Y)
OUTPUT(U1 ,Y,Y,Y,N,Y,Y)
OUTPUT(V1 ,Y,Y,Y,N,Y,Y)
OUTPUT(WDIS,Y,N,N,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 18 ;IYMON = 18 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 1
NYPRIN = 1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
NOWIPE = T
GVIEW(P,0.000000E+00,0.000000E+00,-1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 3.200000E-01, 3.200000E-01, 1.500000E+00
> DOM, MONIT, 1.750000E-01, 1.750000E-01, 7.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, 3.100000E-01, 0.000000E+00, 0.000000E+00
> OBJ1, SIZE, 1.000000E-02, 3.200000E-01, 1.500000E+00
> OBJ1, CLIPART, cube14
> OBJ1, ROTATION24, 1
> OBJ1, GRID, 2.000000E+00
> OBJ1, TYPE, BLOCKAGE
> OBJ1, MATERIAL, 102
> OBJ1, ADIABATIC, 0.000000E+00, 0.000000E+00
> OBJ1, EMISSIVITY, 1.000000E+00
> OBJ2, NAME, B2
> OBJ2, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ2, SIZE, 1.000000E-02, 3.200000E-01, 1.500000E+00
> OBJ2, CLIPART, cube14
> OBJ2, ROTATION24, 1
> OBJ2, GRID, 2.000000E+00
> OBJ2, TYPE, BLOCKAGE
> OBJ2, MATERIAL, 102
> OBJ2, ADIABATIC, 0.000000E+00, 0.000000E+00
> OBJ2, EMISSIVITY, 1.000000E+00
> OBJ3, NAME, B3
> OBJ3, POSITION, 0.000000E+00, 3.100000E-01, 0.000000E+00
> OBJ3, SIZE, 3.200000E-01, 1.000000E-02, 1.500000E+00
> OBJ3, CLIPART, cube14
> OBJ3, ROTATION24, 1
> OBJ3, GRID, 2.000000E+00
> OBJ3, TYPE, BLOCKAGE
> OBJ3, MATERIAL, 198
> OBJ4, NAME, B5
> OBJ4, POSITION, 0.000000E+00, 3.100000E-01, 0.000000E+00
> OBJ4, SIZE, 3.200000E-01, 0.000000E+00, 1.500000E+00
> OBJ4, CLIPART, cube11
> OBJ4, ROTATION24, 1
> OBJ4, GRID, 2.000000E+00
> OBJ4, TYPE, PLATE
> OBJ4, POROSITY, 0.000000E+00
> OBJ4, SIDE, BOTH
> OBJ4, ADIABATIC_L, 0.000000E+00, 0.000000E+00
> OBJ4, EMI_LOW, 1.000000E+00
> OBJ4, ADIABATIC_H, 0.000000E+00, 0.000000E+00
> OBJ4, EMI_HIGH, 1.000000E+00
> OBJ5, NAME, B15
> OBJ5, POSITION, 5.000000E-02, 2.800000E-01, 0.000000E+00
> OBJ5, SIZE, 1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ5, CLIPART, cube
> OBJ5, ROTATION24, 1
> OBJ5, GRID, 2.000000E+00
> OBJ5, TYPE, BLOCKAGE
> OBJ5, MATERIAL, 100
> OBJ5, HEAT_FLUX, 0.000000E+00, 5.000000E+03
> OBJ5, EMISSIVITY, 9.000000E-01
> OBJ6, NAME, B16
> OBJ6, POSITION, 8.000002E-02, 2.800000E-01, 0.000000E+00
> OBJ6, SIZE, 1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ6, CLIPART, cube
> OBJ6, ROTATION24, 1
> OBJ6, GRID, 2.000000E+00
> OBJ6, TYPE, BLOCKAGE
> OBJ6, MATERIAL, 100
> OBJ6, HEAT_FLUX, 0.000000E+00, 5.000000E+03
> OBJ6, EMISSIVITY, 9.000000E-01
> OBJ7, NAME, B17
> OBJ7, POSITION, 1.100000E-01, 2.800000E-01, 0.000000E+00
> OBJ7, SIZE, 1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ7, CLIPART, cube
> OBJ7, ROTATION24, 1
> OBJ7, GRID, 2.000000E+00
> OBJ7, TYPE, BLOCKAGE
> OBJ7, MATERIAL, 100
> OBJ7, HEAT_FLUX, 0.000000E+00, 5.000000E+03
> OBJ7, EMISSIVITY, 9.000000E-01
> OBJ8, NAME, B18
> OBJ8, POSITION, 1.400000E-01, 2.800000E-01, 0.000000E+00
> OBJ8, SIZE, 1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ8, CLIPART, cube
> OBJ8, ROTATION24, 1
> OBJ8, GRID, 2.000000E+00
> OBJ8, TYPE, BLOCKAGE
> OBJ8, MATERIAL, 100
> OBJ8, HEAT_FLUX, 0.000000E+00, 5.000000E+03
> OBJ8, EMISSIVITY, 9.000000E-01
> OBJ9, NAME, B19
> OBJ9, POSITION, 1.700000E-01, 2.800000E-01, 0.000000E+00
> OBJ9, SIZE, 1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ9, CLIPART, cube
> OBJ9, ROTATION24, 1
> OBJ9, GRID, 2.000000E+00
> OBJ9, TYPE, BLOCKAGE
> OBJ9, MATERIAL, 100
> OBJ9, HEAT_FLUX, 0.000000E+00, 5.000000E+03
> OBJ9, EMISSIVITY, 9.000000E-01
> OBJ10, NAME, B20
> OBJ10, POSITION, 2.000000E-01, 2.800000E-01, 0.000000E+00
> OBJ10, SIZE, 1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ10, CLIPART, cube
> OBJ10, ROTATION24, 1
> OBJ10, GRID, 2.000000E+00
> OBJ10, TYPE, BLOCKAGE
> OBJ10, MATERIAL, 100
> OBJ10, HEAT_FLUX, 0.000000E+00, 5.000000E+03
> OBJ10, EMISSIVITY, 9.000000E-01
> OBJ11, NAME, B21
> OBJ11, POSITION, 2.300000E-01, 2.800000E-01, 0.000000E+00
> OBJ11, SIZE, 1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ11, CLIPART, cube
> OBJ11, ROTATION24, 1
> OBJ11, GRID, 2.000000E+00
> OBJ11, TYPE, BLOCKAGE
> OBJ11, MATERIAL, 100
> OBJ11, HEAT_FLUX, 0.000000E+00, 5.000000E+03
> OBJ11, EMISSIVITY, 9.000000E-01
> OBJ12, NAME, B22
> OBJ12, POSITION, 2.600000E-01, 2.800000E-01, 0.000000E+00
> OBJ12, SIZE, 1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ12, CLIPART, cube
> OBJ12, ROTATION24, 1
> OBJ12, GRID, 2.000000E+00
> OBJ12, TYPE, BLOCKAGE
> OBJ12, MATERIAL, 100
> OBJ12, HEAT_FLUX, 0.000000E+00, 5.000000E+03
> OBJ12, EMISSIVITY, 9.000000E-01
> OBJ13, NAME, B14
> OBJ13, POSITION, 2.000000E-02, 2.800000E-01, 0.000000E+00
> OBJ13, SIZE, 1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ13, CLIPART, cube
> OBJ13, ROTATION24, 1
> OBJ13, GRID, 2.000000E+00
> OBJ13, TYPE, BLOCKAGE
> OBJ13, MATERIAL, 100
> OBJ13, HEAT_FLUX, 0.000000E+00, 5.000000E+03
> OBJ13, EMISSIVITY, 9.000000E-01
> OBJ14, NAME, B23
> OBJ14, POSITION, 2.900002E-01, 2.800000E-01, 0.000000E+00
> OBJ14, SIZE, 1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ14, CLIPART, cube
> OBJ14, ROTATION24, 1
> OBJ14, GRID, 2.000000E+00
> OBJ14, TYPE, BLOCKAGE
> OBJ14, MATERIAL, 100
> OBJ14, HEAT_FLUX, 0.000000E+00, 5.000000E+03
> OBJ14, EMISSIVITY, 9.000000E-01
> OBJ15, NAME, B24
> OBJ15, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ15, SIZE, 3.200000E-01, 0.000000E+00, 1.500000E+00
> OBJ15, CLIPART, cube13
> OBJ15, ROTATION24, 1
> OBJ15, GRID, 2.000000E+00
> OBJ15, TYPE, PLATE
> OBJ15, SURF_TEMP, 0.000000E+00, 3.000000E+02
> OBJ15, VELOCITY, 1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ15, EMISSIVITY, 1.000000E+00
> OBJ16, NAME, B26
> OBJ16, POSITION, 1.000000E-02, 0.000000E+00, 0.000000E+00
> OBJ16, SIZE, 1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ16, CLIPART, cubet
> OBJ16, ROTATION24, 1
> OBJ16, TYPE, PRESSURE_RELIEF
> OBJ16, PRES_RELIEF, 1.000000E+03, 0.000000E+00
STOP
inform7begin
(stored var #an1 is anco(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 #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