Difference between revisions of "Evaluation AC3-11"

From KBwiki
Jump to navigation Jump to search
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]]
 
 
        Top              Next
+
{{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

Front Page

Description

Test Data

CFD Simulations

Evaluation

Best Practice Advice

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 (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)


Image335.gif



Appendix D

APPENDIX D: Experimental Temperature and Heat Transfer Results


Image336.gif


Image337.gif


Image338.gif


Image339.gif


Image340.gif


Image341.gif


Image342.gif


Image343.gif


Image344.gif


Image345.gif


Image346.gif


Image347.gif


Image348.gif


Image349.gif


Image350.gif


Image351.gif


Image352.gif


Image353.gif


Image354.gif


Image355.gif



Appendix E

APPENDIX E: Experimental Conditions and Heat Transfer Results


Image356.gif


Image357.gif


Image358.gif


Image359.gif


Image360.gif


Image361.gif


Image362.gif


Image363.gif


Image134.jpg




© copyright ERCOFTAC 2004


Contributors: Mike Rabbitt - British Energy

Site Design and Implementation: Atkins and UniS


Front Page

Description

Test Data

CFD Simulations

Evaluation

Best Practice Advice