Evaluation AC3-11

From KBwiki
Revision as of 16:24, 17 September 2008 by Tonyh (talk | contribs)
Jump to navigation Jump to search

Downward flow in a heated annulus

Application Challenge 3-11 © copyright ERCOFTAC 2004


Comparison of Test data and CFD

The CFD explorations have brought out several important features relevant to the prediction of vertically downward flows through an annulus with a heated core tube.

• The study of Guy et al (1999) has shown that these annular flows (with the possible exception of very low Reynolds number cases) can be treated as axisymmetric, and that for water flows it is important to include the variation of molecular properties with temperature.

• Where buoyant effects are significant the computational studies have highlighted the need for an accurate modeling of the near-wall sublayer. The usual form of wall functions, based on a universal logarithmic wall law, do not account for these buoyancy effects and hence fail to predict the enhancement of heat transfer as the buoyancy parameter is increased.

• The analytical wall function, which accounts for buoyant forces within the viscous sublayer and mimics the effects of changes in the viscous sublayer thickness, returns results in generally good agreement with the data. It thus provides an economic alternative to low-Reynolds-number models, which would otherwise be needed in order to resolve accurately the near-wall region of these annular flows

.

• The present calculations have not included a lengthscale correction term in the ε equation. Earlier computations showed such a term to have some influence, particularly at high heat loadings. Further calculations using the present scheme under such conditions, particularly when reverse flow may occur, should therefore also include a correction of this form.


References

1) Cotton, M.A., Jackson, J.D., 1987, “Calculation of Turbulent Mixed Convection in a Vertical Tube using a Low-Reynolds-Number k-ε turbulence model”, Proc. 6th Turbulent Shear Flows Symposium, Toulouse, France.

2) Gerasimov A.V., 2002, “CFD Quality and Trust: Develoment and Validation of an Analytical Wall-Function Strategy for Modelling Forced, Mixed and Natural Convection Flows”, Report PM/GNSR/5106, Department of Mechanical, Aerospace & Manufacturing Engineering, UMIST.

3) Guy A., Iacovides H., Launder B.E., 1999, “Study of Downward Water Flow in a Heated Annulus, Phase 2”, Technical Report HTH/GNSR/5032, Department of Mechanical Engineering, UMIST.

4) Huang P.G., Leschziner M.A., 1983, “An Introduction and Guide to the Computer Code TEAM”, Report TFD/83/9/(R), Thermofluids Division, Department of Mechanical Engineering, UMIST.

5) Ince, N.Z., Launder, B.E., 1989, “On the Computation of Buoyancy-Driven Turbulent Flows in Rectangular Enclosures”, Int. J. Heat and Fluid Flow, 10, 110-117.

6) Jackson, J.D., Hall, W.B., 1979, “Influences of Buoyancy on Heat Transfer to Fluids Flowing in Vertical Tubes under Turbulent Conditions”, In Turbulent Forced Convection in Channels and Bundles Theory and Applications to Heat Exchangers and Nuclear Reactors, (eds. S. Kakac & D.B. Spalding), Vol 2, pp 613-640.

7) Jackson J.D., He S., Xu Z., Wu T., 2000, “CFD Quality and Trust – Generic Studies of Thermal Convection”, Technical Report HTH/GNSR/5029, School of Engineering, Univeristy of Manchester.

8) Jones, W.P., Launder, B.E., 1972, “The Prediction of Laminarization with a Two-Equation Model of Turbulence”, Int. J. Heat Mass Transfer, 15, 301-314.

9) Launder B.E., Sharma B.I., 1974, “Application of the Energy-Dissipation Model of Turbulence to the Calculation of Flow Near a Spinning Disc”, Letters in Heat and Mass Transfer, 1, pp131-138 .

10) Launder, B.E., 1986, “Low-Reynolds-Number Turbulence near Walls”, Report TFD/86/4, Dept of Mechanical Engineering, UMIST.

11) Leonard B.P., 1979, “A stable and accurate convective modeling procedure based on quadratic upstream interpolation”, Comp. Meth. Appl. Mech. Eng., 19, p59.

12) Patankar S.V., 1980, “Numerical Heat Transfer and Fluid Flow”, Hemisphere Publishing Corporation, Taylor and Francis Group, New York.

13) Patankar S.V., Spalding D.B., 1972, “A Calculation Procedure of Heat, Mass and Momentum transfer in Three-Dimensional Parabolic Flows”, Int J. Heat and Mass Transfer, 15, p1787.

14) Spalding D.B, 1967, “Heat Transfer from Turbulent Separated Flows”, J. Fluid Mechanics, 27, p97.

15) Yap, C.R., 1987, “Turbulent heat and momentum transfer in recirculating and impinging flows”, PhD Thesis, Dept of Mech. Eng., Faculty of Technology, University of Manchester.


Figures

Image323.gif


Figure 4: Variation of Nusselt number and temperature along the heated inner wall for Re=6000, Bo=0.22.


Image325.gif


Figure 5: Variation of Nusselt number and temperature along the heated inner wall for Re=6000, Bo=0.78.


Image327.gif



Figure 6: Variation of Nusselt number and temperature along the heated inner wall for Re=4000, Bo=0.83.


Image329.gif



Figure 7: Variation of Nusselt number and temperature along the heated inner wall for Re=4000, Bo=2.89.


Image331.gif



Figure 8: Profiles of mean vertical velocity across the annular gap at a streamwise position x/deff=35.9. Upper graph: unheated case; Lower graph: buoyant case at Re=4000, Bo=2.89.


Image334.gif


Figure 9: Normalized Mean temperature and vertical velocity profiles close to the inner heated wall at x/deff=35.9 for the buoyant case at Re=4000, Bo=2.89. (Position of main grid nodes is indicated by symbols on the wall function solutions).


Acknowledgement

This work was funded under the HSE Generic Nuclear Safety Research programme and is published with the permission of the UK Nuclear Industry Management Committee (IMC). The authors gratefully acknowledge the financial assistance provided for this investigation. The Manchester University experiments were carried out under the terms of the research Contract entitled ‘CFD Quality and Trust – Generic Studies of Thermal Convection’. The UMIST computational studies were carried out under the terms of the research Contract entitled ‘CFD Quality and Trust – Model Evaluation, Refinement and Application Advice’.


Appendix A

APPENDIX A: LRN MODEL (Low Reynolds Number k-ε Model)


The model of Launder & Sharma employs transport equations for the turbulent kinetic energy and its dissipation rate, and includes effects of molecular viscosity, allowing it to be applied across the viscous sublayer to the wall.

The equations for the turbulent energy and modified dissipation rate Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbases.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \[\tilde{\varepsilon}\]} are:


Failed to parse (syntax error): {\displaystyle \[\frac{D\kappa}{D{t}} = \frac{\delta}{\delta x_{j}}\left[{(\nu + \nu_{t})} \frac{\delta\kappa}{\delta x_{j}}\right] + P_{\kappa} + G_{\kappa} - \tilde{\varepsilon} - 2\nu\left(\frac{\delta\sqrt{\kappa}}{\delta x_{j}}\right)^2\]}


Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbases.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \[\frac{D\tilde{\varepsilon}}{D{t}} = \frac{\delta}{\delta x_{j}}\left[{(\nu + \nu_{t}/\sigma_{\varepsilon})} \frac{\delta\tilde{\varepsilon}}{\delta x_{j}}\right] + c_\varepsilon {1} \frac{\tilde{\varepsilon}}{\kappa} {(P_{\kappa} + G_{\kappa})} - c_{\varepsilon 2} f_2 \frac{\tilde{\varepsilon}^2}{\kappa} + 2\nu\nu_{t}\left(\frac{\delta^2 U_{i}}{\delta x_{j}\delta x_{\kappa}}\right)^2\]}


where the production rates of turbulent kinetic energy due to shear and buoyancy are:

Failed to parse (syntax error): {\displaystyle \[P_{\kappa} = -\overline{u_{i}u_{j}} \frac{\delta U_{i}}{\delta x_{j}}\]}


and Failed to parse (syntax error): {\displaystyle \[G_{\kappa} = -\beta g_{i}\overline{u_{i}\theta}\]}

The turbulent viscosity is then calculated as:


Failed to parse (syntax error): {\displaystyle \[\nu_{t} = c_{\mu} f_{\mu}\frac{\kappa^2}{\varepsilon}\]}


The stresses and turbulent heat fluxes are modeled as:

Failed to parse (syntax error): {\displaystyle \[\overline{u_{i}u_{j}} = \frac{2}{3}\kappa\delta_{j} - \nu_{t}\left(\frac{\delta U_{i}}{\delta x_{j}} + \frac{\delta U_{j}}{\delta x_{i}}\right)\]}

and Failed to parse (syntax error): {\displaystyle \[\overline{u_{i}\theta} = -\frac{\nu_{t}}{\sigma_{t}}\frac{\delta T}{\delta x_{i}}\]}


The various constants and function appearing in the model are:

Failed to parse (syntax error): {\displaystyle \[c_{\mu} = 0.09 c_{\varepsilon 1} = 1.44 c_{\varepsilon 2} = 1.92 \sigma_{\varepsilon} = 1.3 \sigma_{t} = 0.9\]}

Failed to parse (syntax error): {\displaystyle \[f_{\mu} = exp\left[\frac{-3.4}{{(1 + Re_{t}/50)}^2}\right]f_2 = 1-0.3 exp {(-Re^2_{t})}Re_{t} = \frac{\kappa}{\nu\varepsilon}\]}


Appendix B

APPENDIX B: StWF (Standard Wall Function)


The standard wall function employs the log-laws:

where

and

9.24{(Pr/st)3/4-1}[1+0.28exp(-0.007Pr/st)] cm-1/4 9.79 0.41 0.9


In solving the k equation over the near-wall cell, the generation and dissipation terms are replaced by cell-averaged values, defined as

and the dissipation rate at the near-wall node is prescribed as

© ERCOFTAC 2004 Appendix C

APPENDIX C: AWF (Algebraic Wall Function)








© ERCOFTAC 2004 Appendix D APPENDIX D: Experimental Temperature and Heat Transfer Results







































© ERCOFTAC 2004 Appendix E APPENDIX E: Experimental Conditions and Heat Transfer Results















© copyright ERCOFTAC 2004

Contributors: Mike Rabbitt - British Energy

Site Design and Implementation: Atkins and UniS

       Top        	      Next