# Difference between revisions of "Evaluation AC3-11"

David.Fowler (talk | contribs) |
|||

Line 1: | Line 1: | ||

+ | {{AC|front=AC 3-11|description=Description_AC3-11|testdata=Test Data_AC3-11|cfdsimulations=CFD Simulations_AC3-11|evaluation=Evaluation_AC3-11|qualityreview=Quality Review_AC3-11|bestpractice=Best Practice Advice_AC3-11|relatedUFRs=Related UFRs_AC3-11}} | ||

+ | |||

='''Downward flow in a heated annulus'''= | ='''Downward flow in a heated annulus'''= | ||

Line 331: | Line 333: | ||

Site Design and Implementation: [[Atkins]] and [[UniS]] | Site Design and Implementation: [[Atkins]] and [[UniS]] | ||

− | + | {{AC|front=AC 3-11|description=Description_AC3-11|testdata=Test Data_AC3-11|cfdsimulations=CFD Simulations_AC3-11|evaluation=Evaluation_AC3-11|qualityreview=Quality Review_AC3-11|bestpractice=Best Practice Advice_AC3-11|relatedUFRs=Related UFRs_AC3-11}} |

## Revision as of 10:28, 19 March 2009

## Contents

**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**

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

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

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

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

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.

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 (syntax error): {\displaystyle \[\tilde{\varepsilon}\]}**
are:

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

**Failed to parse (syntax error): {\displaystyle \[\frac{D\tilde{\varepsilon}}{D{t}} = \frac{\partial}{\partial x_{j}}\left[{(\nu + \nu_{t}/\sigma_{\varepsilon})} \frac{\partial\tilde{\varepsilon}}{\partial 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{\partial^2 U_{i}}{\partial x_{j}\partial 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{\partial U_{i}}{\partial 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 (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \[\overline{u_{i}u_{j}} = \frac{2}{3}\kappa\partial_{j} - \nu_{t}\left(\frac{\partial U_{i}}{\partial x_{j}} + \frac{\partial U_{j}}{\partial x_{i}}\right)\]}**

and **Failed to parse (syntax error): {\displaystyle \[\overline{u_{i}\theta} = -\frac{\nu_{t}}{\sigma_{t}}\frac{\partial T}{\partial 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:

**Failed to parse (syntax error): {\displaystyle \[U\cdot = \frac{1}{\kappa}log {(E\cdot Y\cdot)} T\cdot = \sigma_{t}{(U\cdot + P\cdot)}\]}**

where

**Failed to parse (syntax error): {\displaystyle \[U\cdot = U\kappa^{1/2}/\nu T\cdot = \sigma C_{p}\kappa^{1/2} {(T_{w}-T)}/q_{w}Y\cdot = y\kappa^{1/2}/\nu\]}**

and

**Failed to parse (syntax error): {\displaystyle \[P\cdot\]}**
= 9.24{(Pr/st)3/4-1}[1+0.28exp(-0.007Pr/st)] cm-1/4 **Failed to parse (syntax error): {\displaystyle \[E\cdot\]}**
= 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

**Failed to parse (syntax error): {\displaystyle \[\overline{P_{\kappa}} = \frac{s^2_{w}}{\kappa{c}^{1/2}_{\mu} {d}\kappa^{1/2}_{p} Y_{n}} log{(Y_{n}/Y_{\nu})}\]}**

**Failed to parse (syntax error): {\displaystyle \[-_{\varepsilon} = \frac{1}{y_{n}}\left[\frac{2{\mu\kappa}_{p}}{y_{\nu}} + \frac{\kappa^{3/2}_{p}}{c_{l}} log {(y_{n}/y{\nu})}\right]\]}**

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

**Failed to parse (syntax error): {\displaystyle \[\varepsilon_{p} = \kappa^{3/2}_{p}/{(c_{l}y_{p})}\]}**

**Appendix C**

**APPENDIX C: AWF (Algebraic Wall Function)**

**Appendix D**

**APPENDIX D: Experimental Temperature and Heat Transfer Results**

**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