Uncertainty quantification for RANS simulation of flow over a wavy wall
|
|
- Chrystal Cain
- 5 years ago
- Views:
Transcription
1 Uncertainty quantification for RANS simulation of flow over a wavy wall Catherine Gorlé 1,2,3, Riccardo Rossi 1,4, and Gianluca Iaccarino 1 1 Center for Turbulence Research, Stanford University, Stanford, CA, 9435, USA. gorle.catherine@gmail.com 2 Von Karman Institute for Fluid Dynamics, Sint-Genesius-Rode, Belgium. 3 EMAT, Department of Physics, University of Antwerp, Antwerp, Belgium. 4 Dipartimento di Ingegneria Energetica Nucleare e del Controllo Ambientale, Università di Bologna, Italy. Abstract The inability of Reynolds-averaged Navier-Stokes (RANS) simulations with two-equation linear eddy viscosity models to correctly predict flow separation and reattachment limits the reliability of such simulations in many wind engineering problems. A methodology to quantify the turbulence model form uncertainty in the RANS simulations would provide valuable additional information. In this study we consider this problem for the flow over a wavy wall. The methodology is based on perturbing the Reynolds stress tensor in the momentum equations, and was shown previously to correctly predict the uncertainty in the location of the reattachment point along the wavy wall for simulations with the SST k ω model Gorlé et al. (212). We first present further details on the comparison of the turbulence kinetic energy obtained from two RANS models, the SST k ω and the realizable k ɛ model, and from direct numerical simulation for the flow over the wavy wall. Subsequently the results of the uncertainty quantification are further evaluated by investigating the variation in the velocity field predicted by the different simulations. 1 Introduction The predictive capability of Computational Fluid Dynamics (CFD) for wind engineering problems is limited by uncertainties in the simulations. The incoming boundary layer characteristics are an obvious example of uncertain parameters for which the influence on the solution should be quantified if one is to base engineering decisions on the simulation outcome. An uncertainty quantification (UQ) study inevitably requires a large set of simulations to correctly capture the influence of the uncertainties and in order to keep the overall computational cost of such a study within reasonable limits, the cost of each individual simulation should be minimized. This can enforce the use of reduced-order models for physical processes that are expensive to resolve, such as turbulence and turbulent mixing. Reynoldsaveraged Navier-Stokes (RANS) models remain the most affordable technique for simulating complex turbulent flows and they are therefore well suited for performing a large set of simulations. However, the use of turbulence models introduces additional uncertainties in the simulation outcome and these epistemic, or model-form, uncertainties, need to be quantified as well. In Emory et al. (211) the idea to quantify the uncertainty related to the turbulence model through the introduction of perturbations in the Reynolds stress tensor was first introduced. The method defines the perturbations in terms of the magnitude of the tensor, i.e. the turbulence kinetic energy, and the eigenvalues and eigenvectors of the normalized anisotropy tensor, thereby being completely 1
2 6 th European and African Wind Engineering Conference 2 independent of the initial turbulence model form. In Gorlé & Iaccarino (213) this framework was extended to include UQ of turbulent scalar flux models for RANS simulations of the mixing of a jet in a supersonic crossflow. It was shown that the uncertainty in a quantity of interest related to the downstream mixing could be captured when defining the range of discrepancies considered based on a LES/RANS comparison of the decomposed Reynolds stresses. The results indicated a number of remaining questions for future research. In particular, the definition of the perturbation functions in a more affordable and general way, relying less on the availability of high-fidelity data, was identified as the main challenge. In Gorlé et al. (212), this problem was further addressed for the simulations of the flow over a wavy wall using the SST k ω model, and the uncertainty in the location of the reattachment point for the flow over the wavy wall was successfully quantified. The present paper presents further details and additional results of this uncertainty quantification study for the flow over a wavy wall. Two turbulence models are considered: the realizable k ɛ and the SST k ω models. The comparison between the turbulence kinetic energy, k, obtained by both RANS models and that computed by direct numerical simulations (DNS) for the flow over a wavy wall (Rossi, 211) is presented. Three different sets of simulations were performed. First, both the RANS flow and turbulence equations were solved. Second, the effect of the coupling to the mean flow was eliminated by freezing the flow to the time-averaged DNS flow field, solving the transport equations for the turbulence quantities only. Third, the turbulence kinetic energy production term was also frozen to the production term calculated from the DNS. Based on this analysis the definition of a perturbation function for k that will be used in the model form UQ study for the flow over the wavy wall was established. Finally, the variation in the velocity field obtained from this model form UQ study for both turbulence models is presented. The details of the configuration and the RANS simulations are presented in Section 2. The methodology used for the model form UQ is summarized in the third section. An overview of the results of the comparative analysis for the turbulence kinetic energy is presented in the Section 4. Section 5 of this paper discusses how the full UQ study was formulated and presents the resulting variation in the velocity field. 2 Flow setup, DNS data and RANS simulations The flow setup considered in the present study is identical to the one used in the DNS simulations by Rossi (211) and is shown in Figure 1.The geometry is defined as a two-dimensional channel, with a flat upper wall and a wavy surface at the lower wall defined by: ( y = A cos 2πx λ ), (1) with A the amplitude and λ the wavelength. The wave steepness is defined by 2A/λ and is equal to.1. The DNS database includes results for the flow and for scalar dispersion at Re = 685, based on the bulk velocity U b at the wave crest and the average channel height H. The comparison presented in Section 4 focuses on the magnitude of the Reynolds stress tensor, which is given by the turbulence kinetic energy, k. In Section 5 the time-averaged velocity is compared to the results obtained with the UQ study. The Reynolds stresses and the velocity obtained from the time-averaged DNS solution were averaged in the spanwise direction and over the 4 wave crests, since the flow is assumed to be fully developed and hence periodic in both directions. The RANS simulations were performed on a two-dimensional mesh of 256 x 192 cells, and the near wall resolution was sufficient (y + 1) to avoid the use of wall functions. In the streamwise and spanwise direction periodic conditions are applied. Two different two-equation turbulence models based on the turbulent viscosity hypothesis have been tested: the SST k-ω and the realizable k-ɛ model.
3 6 th European and African Wind Engineering Conference 3 Figure 1: Flow configuration (Rossi, 211). The uncertainty introduced by these models will be investigated using the methodology from Emory et al. (211), summarized in the following section. 3 Epistemic UQ methodology In the past, the influence of turbulence models has mainly been investigated by comparing different models, often belonging to the same class of two-equation linear eddy viscosity models. Such traditional sensitivity studies can not capture the full model form uncertainty, because this uncertainty is largely determined by the fundamental model assumptions, i.e. those made in the turbulent viscosity hypothesis: u i u j = 2 3 kδ ij 2ν t S ij, (2) where u i u j are the Reynolds stresses, ν t is the turbulent viscosity and S ij is the strain rate tensor. To overcome this limitation, a framework for the epistemic UQ of the Reynolds stresses independent of the initial model form was established in Emory et al. (211). The method consists in directly introducing perturbations in the Reynolds stress tensor computed by the model and used in the momentum equations. The definition of the perturbations is based on the eigenvalue decomposition of the normalized anisotropy tensor a ij : a ij = R ij 2k 1 2 δ ij = v ik Λ kl v jl, (3) where v ij is the matrix of orthonormal eigenvectors and Λ kl is the diagonal matrix of eigenvalues λ l. This decomposition allows reformulating the Reynolds stress tensor as ( ) 1 R ij = 2k 3 δ ij + v ik Λ kl v jl. (4) The above formulation does not involve any modeling assumptions, thereby presenting a general way of introducing epistemic uncertainty in the Reynolds stress tensor by writing the perturbed Reynolds stresses as: ( ) 1 Rij = 2k 3 δ ij + vikλ klvjl. (5) The perturbations are thus specified in terms of a discrepancy in the turbulence kinetic energy k = k + k, the diagonal matrix Λ kl of perturbed eigenvalues λ l, and perturbed eigenvectors v ij = q ikv kj, where q ik is an orthonormal rotation matrix. As stated in the Section 1, the main challenge in this approach is the definition of the perturbation functions in an affordable and general way. For the eigenvalues, the perturbations will be defined in
4 6 th European and African Wind Engineering Conference 4 Figure 2: Eigenvalue perturbation defined as a displacement in the Barycentric map (Banerjee et al., 27). terms of the location of the Reynolds stress tensor in the Barycentric map (Banerjee et al., 27), by specifying a displacement in the map as shown in Figure 2. This imposes clear bounds on the possible perturbations, since all realizable states of turbulence are located inside the triangular map. For the eigenvectors, the perturbations consist of rotations, which are bound between and 36. Since the flow field is two-dimensional, only one eigenvector rotation angle has to be considered. For the turbulence kinetic energy, the bounds are less easily defined, which motivated the investigation presented in the following section, focusing on the evaluation of the discrepancies in k between the RANS and DNS simulations. 4 Analysis of discrepancies in turbulence kinetic energy Figure 3 presents the comparison for the turbulence kinetic energy from the DNS and from the SST k ω and the realizable k ɛ model for the solutions obtained when freezing the flow field to the time-averaged DNS solution (case A) and when combining this with the exact production term P k = u i u j U i x j computed from the DNS (case B). The difference between both sets of simulations is that when only freezing the DNS flow field, P k is still determined from the inexact Reynolds stresses as computed using the turbulent viscosity hypothesis, while in the second case the exact Reynolds stresses are used. The result shows an overprediction of up to 1% in k when using the correct flow field (case A) with the realizable k ɛ model and an underprediction of the same order of magnitude with the SST k ω model. However, when also using the correct production term for k (case B) both models predict k with much higher accuracy. The comparison of the turbulence kinetic energy field presented in Figure 3 shows that when the exact production term is used in the transport equations for the turbulence quantities, a much more accurate solution for the turbulence kinetic energy can be obtained. Based on this result, the UQ study was performed with perturbations introduced in the eigenvalues and eigenvectors of the Reynolds stress tensor, and subsequently the turbulence kinetic energy was perturbed indirectly using the updated Reynolds stress tensor to compute the production term in the transport equations for the turbulent quantities.
5 6 th European and African Wind Engineering Conference 5.4 k DNS k.6 k k RKE case A k RKE case B k RKE case A k RKE case B k SSTKO case A k SSTKO case B k SSTKO case A k SSTKO case B Figure 3: Comparison of k from DNS, realizable k ɛ (RKE) and SST k ω (SSTKO) models for solutions to the turbulence quantity transport equations based on the DNS flow field (case A) and based on the DNS flow field and DNS turbulence kinetic energy production term (case B). 5 UQ study 5.1 Introducing the perturbations First, the introduction of uncertainties in regions of the flow where the turbulence model performs sufficiently accurate was avoided by defining a marker. This marker identifies regions that deviate from parallel shear flow as regions where perturbations should be introduced, since the models are not expected to provide accurate predictions in those regions. The exact definition and validation of the marker can be found in Gorlé et al. (212). Second, instead of exploring the entire possible space of eigenvalues and eigenvectors, only two sets of perturbations were selected. Those corresponded to the perturbations that produce the maximum and minimum turbulence kinetic energy production term integrated over the marked region, since these can be expected to provide the smallest and largest length of the separated region. In Gorlé et al. (212), it was shown that the maximum and minimum integrated production terms are found when moving the eigenvalues to the one- and three-component corners of the map, respectively, without a rotation of the eigenvectors. 5.2 Uncertainty predicted in the flow field In Gorlé et al. (212), the two perturbed simulations as defined above were shown to capture the uncertainty in the size of the separated region and give a very realistic representation of the uncertainty in the streamwise wall shear stress for the SST k ω model. The same can be demonstrated for the realizable k ɛ model. In Figures 4 and 5 we further explore this result by considering the contour plots for the horizontal and vertical velocities obtained from the DNS, from the original realizable k ɛ (RKE) and SST k ω (SSTKO) models, and from the perturbed cases for both models. As mentioned in Section 5.1
6 6 th European and African Wind Engineering Conference 6 DNS RKE RKE C1 RKE C SSTKO SSTKO C SSTKO C Figure 4: Comparison of horizontal velocity from the DNS, the original realizable k ɛ (RKE) and SST k ω (SSTKO) models, and the perturbed cases (C1 and C3). DNS RKE RKE C1 RKE C SSTKO SSTKO C SSTKO C Figure 5: Comparison of horizontal velocity from the DNS, the original realizable k ɛ (RKE) and SST k ω models (SSTKO), and the perturbed cases (C1 and C3). eigenvalue perturbations towards the one-component (C1) and three-component (C3) corners of the Barycentric map were considered. The plots show that quantitatively the uncertainty in the velocity field is well represented, especially in the region below the crest of the wave. In general, the variation obtained with the SST k ω model is larger than that obtained with the realizable k ɛ model. A more quantitative comparison of the horizontal velocity component is presented in Figure 6. The plots show the DNS profiles for the horizontal velocities, and the lines that present the maximum and minimum values obtained from the perturbed runs. The plots confirm the observations from Figure 4, with the uncertainty in the velocity field being well represented by the two perturbed runs at most locations. For the realizable k ɛ model the uncertainty in the velocity field is not fully captured at x/λ =.5, while for the SST k ω model it is only on top of the wave crest that the uncertainty in the velocity field is less well represented. 6 Conclusions The goal of the present work is to quantify turbulence model form uncertainty for RANS simulations of the flow over a wavy wall. The method presented in Emory et al. (211) is used, which is based on introducing perturbations in the Reynolds stress tensor in terms of the magnitude of the tensor, i.e. the turbulence kinetic energy, and the eigenvalues and eigenvectors of the normalized anisotropy tensor. In Gorlé et al. (212), it was shown that this method enables quantifying the uncertainty in the
7 6 th European and African Wind Engineering Conference 7 =. =.25 =.5 = =..5 1 = = = Figure 6: Comparison of horizontal velocity for the realizable k ɛ (top) SST k ω (bottom) models. Solid lines: DNS, black dotted lines: unperturbed models, blue dotted lines: minimum and maximum velocities from the perturbed runs. location of the reattachment point for the flow over the wavy wall for simulations with the SST k ω model. The present paper presents further details and additional results of this uncertainty quantification study considering two turbulence models: the realizable k ɛ and the SST k ω models. For both models a comparison between the modeled turbulence kinetic energy and that computed by direct numerical simulations shows that when the exact production term is used in the transport equations for the turbulence quantities, a much more accurate solution for the turbulence kinetic energy can be obtained. This supports the procedure of indirectly perturbing the turbulence kinetic energy, by using the Reynolds stress tensor obtained after perturbing the eigenvalues and eigenvectors to compute the production term in the transport equations for the turbulent quantities. It is shown that a good representation of the uncertainty in the velocity field can be obtained when performing simulations for the perturbations that result in the minimum and maximum integrated production terms. Acknowledgements This work was supported by the Department of Energy [NNSA] under Award Number NA The first author is supported by a Pegasus-Marie Curie fellowship from the Research Foundation Flanders (FWO), co-funded by the Marie Curie action, since October 212. References Banerjee, S., Krahl, R., Durst, F., & Zenger, C. 27. Presentation of anisotropy properties of turbulence, invariants versus eigenvalue approaches. Journal of Turbulence, 8, Emory, M., Pecnik, R., & Iaccarino, G Modeling Structural Uncertainties in Reynolds Averaged Computations of Shock/Boundary Layer Interactions. AIAA paper, Gorlé, C., & Iaccarino, G A framework for epistemic uncertainty quantification of turbulent scalar flux models for Reynolds-averaged Navier-Stokes simulations. Physics of Fluids.
8 6 th European and African Wind Engineering Conference 8 Gorlé, C., Emory, M., Larsson, J., & Iaccarino, G Epistemic uncertainty quantification of RANS modeling of the flow over a wavy wall. Annual Research Briefs 212, Center for Turbulence Research, Stanford, CA, Rossi, R A numerical study of algebraic flux models for heat and mass transport simulation in complex flows. International Journal of Heat and Mass Transfer, 53,
Componentality-based wall-blocking for RANS models
Center for Turbulence Research Annual Research Briefs 2014 193 Componentality-based wall-blocking for RANS models By M. Emory AND G. Iaccarino 1. Motivation and objectives Although there has been increased
More informationarxiv: v1 [physics.flu-dyn] 19 May 2016
Quantifying Model Form Uncertainty in RANS Simulation of Wing Body Junction Flow Jin-Long Wu, Jian-Xun Wang, Heng Xiao Department of Aerospace and Ocean Engineering, Virginia Tech, Blacksburg, VA 24060,
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationDNS, LES, and wall-modeled LES of separating flow over periodic hills
Center for Turbulence Research Proceedings of the Summer Program 4 47 DNS, LES, and wall-modeled LES of separating flow over periodic hills By P. Balakumar, G. I. Park AND B. Pierce Separating flow in
More informationSimulations for Enhancing Aerodynamic Designs
Simulations for Enhancing Aerodynamic Designs 2. Governing Equations and Turbulence Models by Dr. KANNAN B T, M.E (Aero), M.B.A (Airline & Airport), PhD (Aerospace Engg), Grad.Ae.S.I, M.I.E, M.I.A.Eng,
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationRANS simulations of rotating flows
Center for Turbulence Research Annual Research Briefs 1999 257 RANS simulations of rotating flows By G. Iaccarino, A. Ooi, B. A. Pettersson Reif AND P. Durbin 1. Motivation and objectives Numerous experimental
More informationComputation of separated turbulent flows using the ASBM model
Center for Turbulence Research Proceedings of the Summer Program 2008 365 Computation of separated turbulent flows using the ASBM model By H. Radhakrishnan, R. Pecnik, G. Iaccarino AND S. Kassinos The
More informationLES of turbulent shear flow and pressure driven flow on shallow continental shelves.
LES of turbulent shear flow and pressure driven flow on shallow continental shelves. Guillaume Martinat,CCPO - Old Dominion University Chester Grosch, CCPO - Old Dominion University Ying Xu, Michigan State
More informationAn evaluation of a conservative fourth order DNS code in turbulent channel flow
Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationNumerical Simulation of Flow Around An Elliptical Cylinder at High Reynolds Numbers
International Journal of Fluids Engineering. ISSN 0974-3138 Volume 5, Number 1 (2013), pp. 29-37 International Research Publication House http://www.irphouse.com Numerical Simulation of Flow Around An
More informationNumerical simulations of heat transfer in plane channel flow
Numerical simulations of heat transfer in plane channel flow Najla EL GHARBI 1, 3, a, Rafik ABSI 2, b and Ahmed BENZAOUI 3, c 1 Renewable Energy Development Center, BP 62 Bouzareah 163 Algiers, Algeria
More informationA dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries
Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation
More informationComparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2
Comparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2 Department of Mechanical and Materials Engineering, Military Engineering Institute,
More informationTurbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i
Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)
More informationExplicit algebraic Reynolds stress models for boundary layer flows
1. Explicit algebraic models Two explicit algebraic models are here compared in order to assess their predictive capabilities in the simulation of boundary layer flow cases. The studied models are both
More informationThe mean shear stress has both viscous and turbulent parts. In simple shear (i.e. U / y the only non-zero mean gradient):
8. TURBULENCE MODELLING 1 SPRING 2019 8.1 Eddy-viscosity models 8.2 Advanced turbulence models 8.3 Wall boundary conditions Summary References Appendix: Derivation of the turbulent kinetic energy equation
More informationAnisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale
More informationDivergence free synthetic eddy method for embedded LES inflow boundary condition
R. Poletto*, A. Revell, T. Craft, N. Jarrin for embedded LES inflow boundary condition University TSFP Ottawa 28-31/07/2011 *email: ruggero.poletto@postgrad.manchester.ac.uk 1 / 19 SLIDES OVERVIEW 1 Introduction
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationTurbulence and its modelling. Outline. Department of Fluid Mechanics, Budapest University of Technology and Economics.
Outline Department of Fluid Mechanics, Budapest University of Technology and Economics October 2009 Outline Outline Definition and Properties of Properties High Re number Disordered, chaotic 3D phenomena
More informationThere are no simple turbulent flows
Turbulence 1 There are no simple turbulent flows Turbulent boundary layer: Instantaneous velocity field (snapshot) Ref: Prof. M. Gad-el-Hak, University of Notre Dame Prediction of turbulent flows standard
More informationHybrid LES RANS Method Based on an Explicit Algebraic Reynolds Stress Model
Hybrid RANS Method Based on an Explicit Algebraic Reynolds Stress Model Benoit Jaffrézic, Michael Breuer and Antonio Delgado Institute of Fluid Mechanics, LSTM University of Nürnberg bjaffrez/breuer@lstm.uni-erlangen.de
More informationPreliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data
35th AIAA Fluid Dynamics Conference, June 6 9, 2005/Toronto,Canada Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data Ellen M. Taylor, M. Pino Martín and Alexander
More informationTurbulent flow and convective heat transfer in a wavy wall channel
Turbulent flow and convective heat transfer in a wavy wall channel A.Z. Dellil, A. Azzi, B.A. Jubran Heat and Mass Transfer 40 (2004) 793 799 DOI 10.1007/s00231-003-0474-4 Abstract This paper reports the
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationSimulation of Aeroelastic System with Aerodynamic Nonlinearity
Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School
More informationABSTRACT OF ONE-EQUATION NEAR-WALL TURBULENCE MODELS. Ricardo Heinrich Diaz, Doctor of Philosophy, 2003
ABSTRACT Title of dissertation: CRITICAL EVALUATION AND DEVELOPMENT OF ONE-EQUATION NEAR-WALL TURBULENCE MODELS Ricardo Heinrich Diaz, Doctor of Philosophy, 2003 Dissertation directed by: Professor Jewel
More informationCHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1.
1/3 η 1C 2C, axi 1/6 2C y + < 1 axi, ξ > 0 y + 7 axi, ξ < 0 log-law region iso ξ -1/6 0 1/6 1/3 Figure 11.1: The Lumley triangle on the plane of the invariants ξ and η of the Reynolds-stress anisotropy
More informationDirect numerical simulation of interfacial instability in gas-liquid flows
Direct numerical simulation of interfacial instability in gas-liquid flows Iain Bethune 1, Lennon Ó Náraigh 2, David Scott 1, Peter Spelt 3,4, Prashant Valluri 5, Zlatko Solomenko 3 1 Edinburgh Parallel
More informationRECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD
RECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD N. Jarrin 1, A. Revell 1, R. Prosser 1 and D. Laurence 1,2 1 School of MACE, the University of Manchester,
More informationOn the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows
Center for Turbulence Research Annual Research Briefs 1998 267 On the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows By Jeffrey S. Baggett 1. Motivation and objectives
More informationINVESTIGATION OF THE FLOW OVER AN OSCILLATING CYLINDER WITH THE VERY LARGE EDDY SIMULATION MODEL
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationZonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach
Zonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach S. Arvidson 1,2, L. Davidson 1, S.-H. Peng 1,3 1 Chalmers University of Technology 2 SAAB AB, Aeronautics 3 FOI, Swedish Defence
More informationColloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.
Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš
More informationNumerical simulation of scalar dispersion downstream of a square obstacle
Center for Turbulence Research Annual Research Briefs 8 87 Numerical simulation of scalar dispersion downstream of a square obstacle By R. Rossi AND G. Iaccarino. Motivation and objectives The analysis
More informationarxiv: v2 [physics.flu-dyn] 27 Feb 2017
A Physics Informed Machine Learning Approach for Reconstructing Reynolds Stress Modeling Discrepancies Based on DNS Data Jian-Xun Wang, Jin-Long Wu, Heng Xiao a Department of Aerospace and Ocean Engineering,
More informationMulti-source Aeroacoustic Noise Prediction Method
Multi-source Aeroacoustic Noise Prediction Method Jonathan SCOTT CFD Engineer 03/12/2013 1 Introduction Trend to reduce engine displacement while increasing break power by turbo charging engines Traditionally
More informationSG Turbulence models for CFD
SG2218 2012 Turbulence models for CFD Stefan Wallin Linné FLOW Centre Dept of Mechanics, KTH Dept. of Aeronautics and Systems Integration, FOI There are no simple turbulent flows Turbulent boundary layer:
More informationNumerical Studies of Supersonic Jet Impingement on a Flat Plate
Numerical Studies of Supersonic Jet Impingement on a Flat Plate Overset Grid Symposium Dayton, OH Michael R. Brown Principal Engineer, Kratos/Digital Fusion Solutions Inc., Huntsville, AL. October 18,
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN
International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay
More informationTurbulent boundary layer
Turbulent boundary layer 0. Are they so different from laminar flows? 1. Three main effects of a solid wall 2. Statistical description: equations & results 3. Mean velocity field: classical asymptotic
More informationA TURBULENT HEAT FLUX TWO EQUATION θ 2 ε θ CLOSURE BASED ON THE V 2F TURBULENCE MODEL
TASK QUARTERLY 7 No 3 (3), 375 387 A TURBULENT HEAT FLUX TWO EQUATION θ ε θ CLOSURE BASED ON THE V F TURBULENCE MODEL MICHAŁ KARCZ AND JANUSZ BADUR Institute of Fluid-Flow Machinery, Polish Academy of
More informationNONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS
June - July, 5 Melbourne, Australia 9 7B- NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS Werner M.J. Lazeroms () Linné FLOW Centre, Department of Mechanics SE-44
More informationIntroduction to Turbulence and Turbulence Modeling
Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly
More informationA Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM
206 9th International Conference on Developments in esystems Engineering A Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM Hayder Al-Jelawy, Stefan Kaczmarczyk
More informationAtmospheric Boundary Layer Studies with Unified RANS-LES and Dynamic LES Methods
5st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 7 - January 23, Grapevine (Dallas/Ft. Worth Region), Texas AIAA 23-747 Atmospheric Boundary Layer Studies with
More informationPOTENTIAL TURBULENCE MODEL PREDICTIONS OF FLOW PAST A TRIANGULAR CYLINDER USING AN UNSTRUCTURED STAGGERED MESH METHOD
POTENTIAL TURBULENCE MODEL PREDICTIONS OF FLOW PAST A TRIANGULAR CYLINDER USING AN UNSTRUCTURED STAGGERED MESH METHOD Xg Zhang Blair Perot Department of Mechanical and Industrial Engeerg, University of
More informationTURBULENCE MODEL AND VALIDATION OF AIR FLOW IN WIND TUNNEL
International Journal of Technology (2016) 8: 1362-1372 ISSN 2086-9614 IJTech 2016 TURBULENCE MODEL AND VALIDATION OF AIR FLOW IN WIND TUNNEL Gun Gun Ramdlan G. 1,2*, Ahmad Indra Siswantara 1, Budiarso
More informationTowards an Accurate and Robust Algebraic Structure Based Model
Center for Turbulence Research Proceedings of the Summer Program 22 283 Towards an Accurate and Robust Algebraic Structure Based Model By R. Pecnik, J. P. O. Sullivan, K. Panagiotou, H. Radhakrishnan,
More informationInitial Efforts to Improve Reynolds Stress Model Predictions for Separated Flows
Initial Efforts to Improve Reynolds Stress Model Predictions for Separated Flows C. L. Rumsey NASA Langley Research Center P. R. Spalart Boeing Commercial Airplanes University of Michigan/NASA Symposium
More informationarxiv: v1 [physics.flu-dyn] 11 Oct 2012
Low-Order Modelling of Blade-Induced Turbulence for RANS Actuator Disk Computations of Wind and Tidal Turbines Takafumi Nishino and Richard H. J. Willden ariv:20.373v [physics.flu-dyn] Oct 202 Abstract
More informationLarge eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions
June 30 - July 3, 2015 Melbourne, Australia 9 P-26 Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions Jungwoo Kim Department of Mechanical System Design Engineering
More informationCOMMUTATION ERRORS IN PITM SIMULATION
COMMUTATION ERRORS IN PITM SIMULATION B. Chaouat ONERA, 93 Châtillon, France Bruno.Chaouat@onera.fr Introduction Large eddy simulation is a promising route. This approach has been largely developed in
More informationAlgorithms for Uncertainty Quantification
Algorithms for Uncertainty Quantification Lecture 9: Sensitivity Analysis ST 2018 Tobias Neckel Scientific Computing in Computer Science TUM Repetition of Previous Lecture Sparse grids in Uncertainty Quantification
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationDirect and Large Eddy Simulation of stably stratified turbulent Ekman layers
Direct and Large Eddy Simulation of stably stratified turbulent Ekman layers Stimit Shah, Elie Bou-Zeid Princeton University 64 th APS DFD Baltimore, Maryland Nov 21, 211 Effect of Stability on Atmospheric
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
More informationWind Flow Modeling The Basis for Resource Assessment and Wind Power Forecasting
Wind Flow Modeling The Basis for Resource Assessment and Wind Power Forecasting Detlev Heinemann ForWind Center for Wind Energy Research Energy Meteorology Unit, Oldenburg University Contents Model Physics
More informationLarge Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena
29 Review Large Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena Masahide Inagaki Abstract Computational Fluid Dynamics (CFD) has been applied
More informationA NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS
June 30 - July 3, 2015 Melbourne, Australia 9 7B-4 A NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS C.-Y. Chang, S. Jakirlić, B. Krumbein and C. Tropea Institute of Fluid Mechanics and Aerodynamics /
More informationThe interaction of vorticity and rate-of-strain in homogeneous sheared turbulence
PHYSICS OF FLUIDS VOLUME 12, NUMBER 4 APRIL 2000 The interaction of vorticity and rate-of-strain in homogeneous sheared turbulence K. K. Nomura a) and P. J. Diamessis Department of Mechanical and Aerospace
More informationTransition to turbulence in plane Poiseuille flow
Proceedings of the 55th Israel Annual Conference on Aerospace Sciences, Tel-Aviv & Haifa, Israel, February 25-26, 2015 ThL2T5.1 Transition to turbulence in plane Poiseuille flow F. Roizner, M. Karp and
More informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
More informationOn the transient modelling of impinging jets heat transfer. A practical approach
Turbulence, Heat and Mass Transfer 7 2012 Begell House, Inc. On the transient modelling of impinging jets heat transfer. A practical approach M. Bovo 1,2 and L. Davidson 1 1 Dept. of Applied Mechanics,
More informationDIRECT NUMERICAL SIMULATIONS OF HIGH SPEED FLOW OVER CAVITY. Abstract
3 rd AFOSR International Conference on DNS/LES (TAICDL), August 5-9 th, 2001, Arlington, Texas. DIRECT NUMERICAL SIMULATIONS OF HIGH SPEED FLOW OVER CAVITY A. HAMED, D. BASU, A. MOHAMED AND K. DAS Department
More informationConvection in Three-Dimensional Separated and Attached Flow
Convection in Three-Dimensional Separated and Attached Flow B. F. Armaly Convection Heat Transfer Laboratory Department of Mechanical and Aerospace Engineering, and Engineering Mechanics University of
More informationComputation of hypersonic shock boundary layer interaction on a double wedge using a differential Reynolds Stress Model
Computation of hypersonic shock boundary layer interaction on a double wedge using a differential Reynolds Stress Model A. Bosco AICES, Templergraben 55, Aachen, 52056, Germany B. Reinartz CATS, Theaterplatz
More informationOptimizing calculation costs of tubulent flows with RANS/LES methods
Optimizing calculation costs of tubulent flows with RANS/LES methods Investigation in separated flows C. Friess, R. Manceau Dpt. Fluid Flow, Heat Transfer, Combustion Institute PPrime, CNRS University
More informationCFD analysis of convective heat transfer at the surfaces of a cube immersed in a turbulent boundary layer
Accepted for publication in International Journal of Heat and Mass Transfer CFD analysis of convective heat transfer at the surfaces of a cube immersed in a turbulent boundary layer Thijs Defraeye a, *,
More informationTURBULENT SCALAR MIXING IN A SKEWED JET IN CROSSFLOW: EXPERIMENTS AND MODELING
June 30 - July 3, 2015 Melbourne, Australia 9 4C-4 TURBULENT SCALAR MIXING IN A SKEWED JET IN CROSSFLOW: EXPERIMENTS AND MODELING Kevin J. Ryan Dept. of Mechanical Engineering Stanford University 488 Panama
More informationTURBULENT VORTEX SHEDDING FROM TRIANGLE CYLINDER USING THE TURBULENT BODY FORCE POTENTIAL MODEL
Proceedgs of ASME FEDSM ASME 2 Fluids Engeerg Division Summer Meetg June 11-15, 2 Boston, Massachusetts FEDSM2-11172 TURBULENT VORTEX SHEDDING FROM TRIANGLE CYLINDER USING THE TURBULENT BODY FORCE POTENTIAL
More informationTurbulent Rankine Vortices
Turbulent Rankine Vortices Roger Kingdon April 2008 Turbulent Rankine Vortices Overview of key results in the theory of turbulence Motivation for a fresh perspective on turbulence The Rankine vortex CFD
More informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC
More informationTurbulent eddies in the RANS/LES transition region
Turbulent eddies in the RANS/LES transition region Ugo Piomelli Senthil Radhakrishnan Giuseppe De Prisco University of Maryland College Park, MD, USA Research sponsored by the ONR and AFOSR Outline Motivation
More informationCalculations on a heated cylinder case
Calculations on a heated cylinder case J. C. Uribe and D. Laurence 1 Introduction In order to evaluate the wall functions in version 1.3 of Code Saturne, a heated cylinder case has been chosen. The case
More informationcompression corner flows with high deflection angle, for example, the method cannot predict the location
4nd AIAA Aerospace Sciences Meeting and Exhibit 5-8 January 4, Reno, Nevada Modeling the effect of shock unsteadiness in shock-wave/ turbulent boundary layer interactions AIAA 4-9 Krishnendu Sinha*, Krishnan
More informationModel Studies on Slag-Metal Entrainment in Gas Stirred Ladles
Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Anand Senguttuvan Supervisor Gordon A Irons 1 Approach to Simulate Slag Metal Entrainment using Computational Fluid Dynamics Introduction &
More informationAdvanced RANS modeling of wingtip vortex flows
Center for Turbulence Research Proceedings of the Summer Program 2006 73 Advanced RANS modeling of wingtip vortex flows By A. Revell, G. Iaccarino AND X. Wu The numerical calculation of the trailing vortex
More informationTurbulence - Theory and Modelling GROUP-STUDIES:
Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence
More informationINFLUENCE OF THE BOUNDARY CONDITIONS ON A TEMPERATURE FIELD IN THE TURBULENT FLOW NEAR THE HEATED WALL
International Conference Nuclear Energy for New Europe 22 Kranjska Gora, Slovenia, September 9-2, 22 www.drustvo-js.si/gora22 INFLUENCE OF THE BOUNDARY CONDITIONS ON A TEMPERATURE FIELD IN THE TURBULENT
More informationOn modeling pressure diusion. in non-homogeneous shear ows. By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3
Center for Turbulence Research Proceedings of the Summer Program 1996 63 On modeling pressure diusion in non-homogeneous shear ows By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3 New models
More informationA Solution Method for the Reynolds-Averaged Navier-Stokes Equation
A Solution Method for the Reynolds-Averaged Navier-Stokes Equation T.-W. Lee Mechanical and Aerospace Engineering, Arizona State University, Tempe, AZ, 8587 Abstract- Navier-Stokes equations are difficult
More informationModeling Complex Flows! Direct Numerical Simulations! Computational Fluid Dynamics!
http://www.nd.edu/~gtryggva/cfd-course/! Modeling Complex Flows! Grétar Tryggvason! Spring 2011! Direct Numerical Simulations! In direct numerical simulations the full unsteady Navier-Stokes equations
More informationWALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationInteraction(s) fluide-structure & modélisation de la turbulence
Interaction(s) fluide-structure & modélisation de la turbulence Pierre Sagaut Institut Jean Le Rond d Alembert Université Pierre et Marie Curie- Paris 6, France http://www.ida.upmc.fr/~sagaut GDR Turbulence
More informationmeters, we can re-arrange this expression to give
Turbulence When the Reynolds number becomes sufficiently large, the non-linear term (u ) u in the momentum equation inevitably becomes comparable to other important terms and the flow becomes more complicated.
More informationCOMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE
COMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE L. Velázquez-Araque 1 and J. Nožička 2 1 Division of Thermal fluids, Department of Mechanical Engineering, National University
More informationOn the validation study devoted to stratified atmospheric flow over an isolated hill
On the validation study devoted to stratified atmospheric flow over an isolated hill Sládek I. 2/, Kozel K. 1/, Jaňour Z. 2/ 1/ U1211, Faculty of Mechanical Engineering, Czech Technical University in Prague.
More informationNovel uses of DNS with turbulent separation for RANS models
Novel uses of DNS with turbulent separation for RANS models Gary Coleman*, Chris Rumsey* and Philippe Spalart** *NASA Langley Research Center **Boeing Commercial Airplanes 1 Outline q Flow: spatial separation
More informationConjugate Heat Transfer Simulation of Internally Cooled Gas Turbine Vane
Conjugate Heat Transfer Simulation of Internally Cooled Gas Turbine Vane V. Esfahanian 1, A. Shahbazi 1 and G. Ahmadi 2 1 Department of Mechanical Engineering, University of Tehran, Tehran, Iran 2 Department
More informationOpenFOAM selected solver
OpenFOAM selected solver Roberto Pieri - SCS Italy 16-18 June 2014 Introduction to Navier-Stokes equations and RANS Turbulence modelling Numeric discretization Navier-Stokes equations Convective term {}}{
More informationROLE OF THE VERTICAL PRESSURE GRADIENT IN WAVE BOUNDARY LAYERS
ROLE OF THE VERTICAL PRESSURE GRADIENT IN WAVE BOUNDARY LAYERS Karsten Lindegård Jensen 1, B. Mutlu Sumer 1, Giovanna Vittori 2 and Paolo Blondeaux 2 The pressure field in an oscillatory boundary layer
More informationAER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly
1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence
More information