TALK=T;RUN( 1, 1)
DISPLAY
CONVECTION AFFECTED BY RADIATION
An IMMERSOL application is presented aimed at the demonstration
of the method for the radiation-affected buoyancy-induced
convection in a rectangular enclosure with thick conductive
walls.
This case deals with conjugate heat transfer in the free-
convected flow of participating gas arising in a rectangular
cavity the thick walls of which are made from aluminium with an
emmitance of 0.9.
The outer surfaces of the bottom and top walls are supposed
to be well insulated. The temperatures of side walls are
different: one has 300K at the outer surface, another is fixed
at 600K.
The task is to compare the velocity and temperature
distributions for the different absorptive and scattering
properties of the gas.
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-direcyion heat flux (qrx) contours
con qrx z 1 fi;0.001 ; pause;cl
msg y-direcyion 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(Convection affected by radiation)
************************************************************
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(144) =PRPS ; NAME(145) =WDIS
NAME(146) =WGAP ; NAME(147) =T3
NAME(148) =QRY ; NAME(149) =QRX
NAME(150) =TEM1
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,T3 ,TEM1)
* Stored variables list
STORE(QRX ,QRY ,WGAP,WDIS,PRPS)
* 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)
************************************************************
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
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(U1 ) = 1.000000E-01 ;FIINIT(V1 ) = 1.000000E-01
FIINIT(WDIS) = 3.200000E-02 ;FIINIT(TEM1) = 3.000000E+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
PATCH (BUOYANCY,PHASEM,0,0,0,0,0,0,1,1)
COVAL (BUOYANCY,U1 , FIXFLU , GRND3 )
COVAL (BUOYANCY,V1 , FIXFLU , GRND3 )
BUOYA = 9.810000E+00 ; BUOYB = 0.000000E+00
BUOYC = 0.000000E+00
BUOYD = 3.410000E-03
BUOYE = 2.000000E+01
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 2000
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
RELAX(U1 ,FALSDT, 1.000000E-01)
RELAX(V1 ,FALSDT, 1.000000E-01)
RELAX(T3 ,LINRLX, 5.000000E-01)
RELAX(TEM1,LINRLX, 5.000000E-01)
************************************************************
Group 18. Limits
VARMAX(T3 ) = 6.000000E+02 ;VARMIN(T3 ) = 3.000000E+02
VARMAX(TEM1) = 6.000000E+02 ;VARMIN(TEM1) = 3.000000E+02
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
RADIA = 1.000000E+00 ;RADIB = 1.000000E-01
************************************************************
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 = 16 ;IYMON = 16 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
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,-1.000000E+00,0.000000E+00,0.000000E+00)
> DOM, SIZE, 3.200000E-01, 3.200000E-01, 1.500000E+00
> DOM, MONIT, 1.550000E-01, 1.550000E-01, 7.500000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_X_1, 32, 1.000000E+00
> GRID, RSET_Y_1, 32, 1.000000E+00
> DOM, RELAX, 5.000000E-01
> OBJ1, NAME, B14
> OBJ1, POSITION, 5.000000E-02, 2.700000E-01, 0.000000E+00
> OBJ1, SIZE, 2.200000E-01, 4.999998E-02, 1.500000E+00
> OBJ1, CLIPART, cube4
> OBJ1, ROTATION24, 1
> OBJ1, GRID, 2.000000E+00
> OBJ1, TYPE, BLOCKAGE
> OBJ1, MATERIAL, 100
> OBJ1, FIXED_TMP, 0.000000E+00, 6.000000E+02
> OBJ1, EMISSIVITY, 9.000000E-01
> OBJ2, NAME, B24
> OBJ2, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ2, SIZE, 3.200000E-01, 0.000000E+00, 1.500000E+00
> OBJ2, CLIPART, cube13
> OBJ2, ROTATION24, 1
> OBJ2, GRID, 2.000000E+00
> OBJ2, TYPE, PLATE
> OBJ2, SURF_TEMP, 0.000000E+00, 3.000000E+02
> OBJ2, EMISSIVITY, 1.000000E+00
> OBJ3, NAME, B25
> OBJ3, POSITION, 5.000000E-02, 5.000000E-02, 0.000000E+00
> OBJ3, SIZE, 1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ3, CLIPART, cubet
> OBJ3, ROTATION24, 1
> OBJ3, VISIBLE, -1.000000E+00
> OBJ3, TYPE, PRESSURE_RELIEF
> OBJ3, PRES_RELIEF, 1.000000E+00, 0.000000E+00
> OBJ4, NAME, B26
> OBJ4, POSITION, 5.000000E-02, 0.000000E+00, 0.000000E+00
> OBJ4, SIZE, 2.200000E-01, 5.000000E-02, 1.500000E+00
> OBJ4, CLIPART, cube14
> OBJ4, ROTATION24, 1
> OBJ4, GRID, 2.000000E+00
> OBJ4, TYPE, BLOCKAGE
> OBJ4, MATERIAL, 100
> OBJ4, ADIABATIC, 0.000000E+00, 0.000000E+00
> OBJ4, EMISSIVITY, 9.000000E-01
> OBJ5, NAME, B5
> OBJ5, POSITION, 2.700000E-01, 5.000000E-02, 0.000000E+00
> OBJ5, SIZE, 5.000000E-02, 2.200000E-01, 1.500000E+00
> OBJ5, CLIPART, cube14
> OBJ5, ROTATION24, 1
> OBJ5, GRID, 2.000000E+00
> OBJ5, TYPE, BLOCKAGE
> OBJ5, MATERIAL, 100
> OBJ5, ADIABATIC, 0.000000E+00, 0.000000E+00
> OBJ5, EMISSIVITY, 9.000000E-01
> OBJ6, NAME, B6
> OBJ6, POSITION, 0.000000E+00, 5.000000E-02, 0.000000E+00
> OBJ6, SIZE, 5.000000E-02, 2.200000E-01, 1.500000E+00
> OBJ6, CLIPART, cube14
> OBJ6, ROTATION24, 1
> OBJ6, GRID, 2.000000E+00
> OBJ6, TYPE, BLOCKAGE
> OBJ6, MATERIAL, 100
> OBJ6, ADIABATIC, 0.000000E+00, 0.000000E+00
> OBJ6, EMISSIVITY, 9.000000E-01
> OBJ7, NAME, B7
> OBJ7, POSITION, 2.700000E-01, 0.000000E+00, 0.000000E+00
> OBJ7, SIZE, 5.000000E-02, 5.000000E-02, 1.500000E+00
> OBJ7, CLIPART, cube14
> OBJ7, ROTATION24, 1
> OBJ7, GRID, 2.000000E+00
> OBJ7, TYPE, BLOCKAGE
> OBJ7, MATERIAL, 198
> OBJ8, NAME, B8
> OBJ8, POSITION, 0.000000E+00, 2.700000E-01, 0.000000E+00
> OBJ8, SIZE, 5.000000E-02, 5.000000E-02, 1.500000E+00
> OBJ8, CLIPART, cube14
> OBJ8, ROTATION24, 1
> OBJ8, GRID, 2.000000E+00
> OBJ8, TYPE, BLOCKAGE
> OBJ8, MATERIAL, 198
> OBJ9, NAME, B9
> OBJ9, POSITION, 2.700000E-01, 2.700000E-01, 0.000000E+00
> OBJ9, SIZE, 5.000000E-02, 5.000000E-02, 1.500000E+00
> OBJ9, CLIPART, cube14
> OBJ9, ROTATION24, 1
> OBJ9, GRID, 2.000000E+00
> OBJ9, TYPE, BLOCKAGE
> OBJ9, MATERIAL, 198
> OBJ10, NAME, B10
> OBJ10, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ10, SIZE, 5.000000E-02, 5.000000E-02, 1.500000E+00
> OBJ10, CLIPART, cube14
> OBJ10, ROTATION24, 1
> OBJ10, GRID, 2.000000E+00
> OBJ10, TYPE, BLOCKAGE
> OBJ10, MATERIAL, 198
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
STOP