|
![](static/image/common/ico_lz.png)
楼主 |
发表于 2004-7-27 18:48:31
|
显示全部楼层
NUSSELT NUMBER COMPUTATION FOR TURBULENT HEAT TRANSFER
TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.5, Date 08/10/02
CPVNAM=VDI;SPPNAM=Core
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(INFORM - Nusselt number computation )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,40,40)
* Cylindrical-polar grid
CARTES=F
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(135) =HTCB ; NAME(137) =ENUL
NAME(138) =TWAL ; NAME(139) =YPLS
NAME(140) =STAN ; NAME(141) =HTCO
NAME(142) =SPH1 ; NAME(143) =DEN1
NAME(144) =ENUT ; NAME(145) =QWAL
NAME(146) =MACH ; NAME(147) =EPKE
NAME(148) =PRPS ; NAME(149) =TEM1
NAME(150) =KOND
* Solved variables list
SOLVE(P1 ,V1 ,W1 ,TEM1)
* Stored variables list
STORE(KOND,PRPS,EPKE,MACH,QWAL,ENUT,DEN1,SPH1)
STORE(HTCO,STAN,YPLS,TWAL,ENUL,HTCB)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
TURMOD(KEMODL)
************************************************************
Group 8. Terms & Devices
NEWRH1 = T
NEWENL = T
DENPCO = T
************************************************************
Group 9. Properties
PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02
SETPRPS(1, 2)
DRH1DP = GRND5
DVO1DT = 3.330000E-03
PRNDTL(TEM1) = -2.630000E-02
PRT (EP ) = 1.314000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS) = -1.000000E+00
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 (SURF ,CELL ,5,0,0,0,0,0,1,2)
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 700
SARAH = 5.000000E-02
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER (P1 ) = 50 ;LITER (TEM1) = 50
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
RELAX(V1 ,FALSDT, 1.000000E-02)
RELAX(W1 ,FALSDT, 1.000000E-02)
RELAX(KE ,LINRLX, 5.000000E-01)
RELAX(EP ,LINRLX, 5.000000E-01)
RELAX(TEM1,LINRLX, 3.000000E-01)
KELIN = 3
************************************************************
Group 18. Limits
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
GENK = T
ASAP = T
IDISPB = 1 ;IDISPC = 30
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 40 ;IZMON = 30
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NYPRIN = 1
IYPRF = 1 ;IYPRL = 40
NZPRIN = 1
IZPRF = 28 ;IZPRL = 36
YZPR = T
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 2
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
NOWIPE = T
IDISPA = 1 ;IDISPB = 1 ;IDISPC = 30
CSG1 ='N'
CSG2 ='C'
GVIEW(P,-9.004471E-01,0.000000E+00,-4.349655E-01)
GVIEW(UP,4.349655E-01,0.000000E+00,-9.004471E-01)
> DOM, SIZE, 1.000000E-01, 1.600000E-01, 2.000000E-01
> DOM, MONIT, 5.000000E-02, 1.580000E-01, 1.475000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> OBJ, NAME, OUTLETSI
> OBJ, POSITION, 0.000000E+00, 1.600000E-01, 0.000000E+00
> OBJ, SIZE, 1.000000E-01, 0.000000E+00, 1.500000E-01
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, 2.500000E+01
> OBJ, COEFFICIENT, 1.000000E+03
> OBJ, TURBULENCE, -1.026000E+04,-1.026000E+04
> OBJ, NAME, JET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E-01, 2.000000E-02, 0.000000E+00
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 4.000000E+01
> OBJ, TEMPERATURE, 5.000000E+02
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, CONCRETE
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.500000E-01
> OBJ, SIZE, 1.000000E-01, 1.600000E-01, 5.000000E-02
> OBJ, CLIPART, cube14
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 122
> OBJ, INI_TEMP, 2.500000E+01
> OBJ, NAME, OUTTOP
> OBJ, POSITION, 0.000000E+00, 2.000000E-02, 0.000000E+00
> OBJ, SIZE, 1.000000E-01, 1.400000E-01, 0.000000E+00
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, 2.500000E+01
> OBJ, COEFFICIENT, 1.000000E+03
> OBJ, TURBULENCE, -1.026000E+04,-1.026000E+04
> OBJ, NAME, SURF
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.500000E-01
> OBJ, SIZE, 1.000000E-01, 1.600000E-01, 0.000000E+00
> OBJ, CLIPART, poldef
> OBJ, ROTATION24, 1
> OBJ, TYPE, USER_DEFINED
STOP
The case considered is a steady, turbulent, axisymmetric
heated jet impinging on a block of lightweight concrete.
The task is to predict the distribution of the heat-transfer
coefficient and Nusselt number along the block. The Reynolds
number of the jet is 20,000 and it discharges at 500 deg C
into ambient at 25 deg C. The jet is located 3.75 nozzle
diameters above the concrete block, which has a thickness of
1.25 nozzle diameters. The standard k-e model is used in the
calculations. A polar computational grid of 40 radial and 40
axial cells is used. The outer boundary of the solution domain
is located 4 nozzle diameters from the jet axis. All model
settings have been made in VR-based Q1 of PHOENICS 3.5. The
IN-FORM facility is used to compute both the Nusselt number
and heat transfer coefficient distribution along the surface
of the block.
** Compute Nusselt number distribution along plate
on the last sweep
CHAR(PART1,PART2)
** +1 refers to high-z slab otherwise current z slab
PART2=((KOND[&&+1]/(0.5*DZW[&&+1]))*TEM1[&&+1]+HTCO*TEM1)
PART1=((KOND[&&+1]/(0.5*DZW[&&+1]))+HTCO)
** compute surface temperature
(STORED OF TWAL at SURF is PART2/PART1 with IF(ISWEEP.EQ.LSWEEP))
** compute surface heat flux
(STORED OF QWAL at SURF is HTCO*(TEM1-TWAL) with IF(ISWEEP.EQ.LSWEE$
P))
** compute bulk heat-transfer coefficient
CHAR(TJET,TAMB,DIAM);TJET=500.0;TAMB=25.0;DIAM=0.04
(STORED OF HTCB at SURF is ABS(QWAL/(:TJET:-TWAL)) with IF(ISWEEP.E$
Q.LSWEEP))
** compute bulk nusselt number
(STORED OF NUSS at SURF is HTCB* IAM:/KOND with IF(ISWEEP.EQ.LSWEE$
P))
|
|