TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 3.6, Date 11/05/06 CPVNAM=VDI;SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY BU$ OYANCY-INDUCED COOLING OF RADIATING ROD An IMMERSOL application is presented aimed at the demonstration of the method for the situation when the thermal radiation is coupled with natural convection to cool down the heated rod. This case deals with the free-convection flow of non-absorbing air arising in an annular cavity the inner aluminium rod of which is heated by the total heat flux of 150W. The emittance of the rod is 0.5 and the outer surface of black cylindrical wall is maintained at 20C. ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 402 ************************************************************ Group 1. Run Title TEXT(Buoyancy-induced cooling of a radiating ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,20,20,1) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T * Non-default variable names NAME(142) =SCAT ; NAME(143) =EMIS 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,EMIS,SCAT) * 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 PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02 * Domain material index is 0 signifying: * Air at 20 deg C, 1 atm, treated as incompressible SETPRPS(1, 0) RHO1 = 1.189000E+00 RHO1A = 0.000000E+00 ;RHO1B = 3.496503E-03 RHO1C = 7.142857E-01 ENUL = 1.544000E-05 ENULA = 1.460000E-06 ;ENULB = 1.100000E+02 ENULC = 0.000000E+00 DISWAL 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(SCAT) = 0.000000E+00 ;FIINIT(EMIS) = 0.000000E+00 FIINIT(WDIS) = 3.200000E-03 fiinit(tem1)=20; fiinit(t3)=20 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 = 0.000000E+00 ; BUOYB =-9.810000E+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 = 1000 RESFAC = 1.000000E-03 ************************************************************ Group 16. Terminate Iterations LITER (T3 ) = 200 ;LITER (TEM1) = 200 ENDIT (T3 ) = 1.000000E-05 ;ENDIT (TEM1) = 1.000000E-05 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 1.000000E+00) RELAX(T3 ,LINRLX, 7.500000E-01) RELAX(TEM1,LINRLX, 7.500000E-01) ************************************************************ Group 18. Limits VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06 VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06 VARMAX(TEM1) = 1.000000E+10 ;VARMIN(TEM1) = 2.000000E+01 ************************************************************ Group 19. EARTH Calls To GROUND Station USEGRD = T ;USEGRX = T ASAP = T PARSOL = T CONWIZ = T SPEDAT(SET,MAXINC,T3,R,1.00000E+02) SPEDAT(SET,MAXINC,TEM1,R,1.00000E+02) ************************************************************ 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 = 11 ;IYMON = 17 ;IZMON = 1 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,5.099395E-03,6.618215E-06,-9.999870E-01) GVIEW(UP,-6.182039E-04,9.999998E-01,3.469538E-06) > DOM, SIZE, 3.200000E-02, 3.200000E-02, 1.400000E+00 > DOM, MONIT, 1.680000E-02, 2.640000E-02, 7.000000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 3.000000E-02 > DOM, SNAPSIZE, 1.000000E-02 > GRID, RSET_X_1, 20, 1.000000E+00 > GRID, RSET_Y_1, 20, 1.000000E+00 > OBJ, NAME, B1 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 3.200000E-02, 3.200000E-02, 1.400000E+00 > OBJ, ROT-MODE, OLD > OBJ, GEOMETRY, wallpipe > OBJ, ROTATION24, 1 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c > OBJ, EMISSIVITY, 5.000000E-01 > OBJ, NAME, B2 > OBJ, POSITION, 1.600000E-02, 1.600000E-02, 0.000000E+00 > OBJ, SIZE, 1.600000E-02, 1.600000E-02, 1.400000E+00 > OBJ, ROT-MODE, OLD > OBJ, GEOMETRY, hh1 > OBJ, ROTATION24, 1 > OBJ, GRID, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c > OBJ, LINR_HEVT, 1.000000E+06, 2.000000E+01 > OBJ, EMISSIVITY, 1.000000E+00 > OBJ, NAME, B3 > OBJ, POSITION, 0.000000E+00, 1.600000E-02, 0.000000E+00 > OBJ, SIZE, 1.600000E-02, 1.600000E-02, 1.400000E+00 > OBJ, ROT-MODE, OLD > OBJ, GEOMETRY, hh1 > OBJ, ROTATION24, 2 > OBJ, GRID, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c > OBJ, LINR_HEVT, 1.000000E+06, 2.000000E+01 > OBJ, EMISSIVITY, 1.000000E+00 > OBJ, NAME, B4 > OBJ, POSITION, 1.600000E-02, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.600000E-02, 1.600000E-02, 1.400000E+00 > OBJ, ROT-MODE, OLD > OBJ, GEOMETRY, hh1 > OBJ, ROTATION24, 4 > OBJ, GRID, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c > OBJ, LINR_HEVT, 1.000000E+06, 2.000000E+01 > OBJ, EMISSIVITY, 1.000000E+00 > OBJ, NAME, B5 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.600000E-02, 1.600000E-02, 1.400000E+00 > OBJ, ROT-MODE, OLD > OBJ, GEOMETRY, hh1 > OBJ, ROTATION24, 3 > OBJ, GRID, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c > OBJ, LINR_HEVT, 1.000000E+06, 2.000000E+01 > OBJ, EMISSIVITY, 1.000000E+00 > OBJ, NAME, B6 > OBJ, POSITION, 1.300000E-02, 1.300000E-02, 0.000000E+00 > OBJ, SIZE, 6.000000E-03, 6.000000E-03, 1.400000E+00 > OBJ, ROT-MODE, OLD > OBJ, GEOMETRY, cylinder > OBJ, ROTATION24, 1 > OBJ, GRID, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c > OBJ, HEAT_FLUX, 0.000000E+00, 1.500000E+02 > OBJ, EMISSIVITY, 5.000000E-01 > OBJ, NAME, B7 > OBJ, POSITION, 3.200000E-03, 1.440000E-02, 0.000000E+00 > OBJ, SIZE, 1.600000E-03, 1.600000E-03, 1.400000E+00 > OBJ, ROT-MODE, OLD > OBJ, GEOMETRY, cubet > OBJ, ROTATION24, 1 > OBJ, TYPE, PRESSURE_RELIEF > OBJ, PRES_RELIEF, 1.000000E+03, 0.000000E+00 STOP