Numerical modeling of complex turbulent flows
|
|
- Calvin Stephens
- 5 years ago
- Views:
Transcription
1 ISTP-16, 5, PRAGUE 16 TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA Numerical modeling of complex turbulent flows Karel Kozel Petr Louda Jaromír Příhoda Dept. of Technical Mathematics CTU Prague, Karlovo nám. 13, CZ Prague Institute of Thermomechanics AS CR, Dolejškova 5, CZ-18 Prague 8 Corresponding author: P. Louda, louda@marian.fsik.cvut.cz Keywords: jet in cross-flow, impinging jet, artificial compressibility, model,, heat transfer Abstract The work deals with the solution of two cases of turbulent incompressible flows: impinging jet flow and round jet in cross-flow. The flows are considered turbulent and statistically steady. The mathematical model is based on Reynolds averaged Navier-Stokes equations with a twoequation turbulence model. The turbulence models are low Reynolds number eddy-viscosity one an model. The system of equations is solved by artificial compressibility method. The discretisation uses finite volume method with higher order upwind approximation for convective term and implicit discretisation in time. For both cases results are compared with experimental data. In the case of impinging jet flow, the heat transfer from impingement wall is also solved. 1 Introduction In this work, we investigate two cases of 3D turbulent flow using the Reynolds averaging approach. The first one is round jet in cross-flow, the second one impinging jet flow. Using sufficiently accurate implicit finite volume method, we study the inluence of turbulence modeling on the results. We use one eddy viscosity and explicit algebraic Reynolds stress model for turbulent momentum trabsfer, and three different turbulent heat flux clusores. Mathematical model The turbulent flow-field in both cases is assumed to be statistically steady. Then the mean flowfield is governed by Reynolds averaged Navier- Stokes (RANS) equations in Cartesian coordinates x i u i t + u ju i ui x i =, = p + (νs ij τ ij ), x i S ij = 1 ( ui + u ) j, (1) x i where u i is mean velocity vector, p kinematic pressure, ν laminar kinematic viscosity, and τ ij Reynolds stress tensor. We apply two closures for τ ij : the (extended) eddy viscosity one of turbulence model and explicit algebraic Reynolds stress model () due to Wallin, Hellsten. The Reynolds stress by model reads τ ij = ν T S ij + 3 kδ ij, ν T = a 1 k max(a 1 ω, ΩF ), a 1 =.31, () where k is turbulent kinetic energy, ω specific dissipation rate, Ω magnitude of vorticity, and F a switching function given in [6]. The 1
2 constitutive relation is given in terms of anisotropy tensor a ij [11, ] τ ij = a ij k + 3 kδ ij, (3) a ij = β 1 τ T S ij + + β 3 τ (Ω ik Ω kj II Ω δ ij /3) + + β 4 τ (S ik Ω kj Ω ik S kj ) + + β 6 τ 3 (S ik Ω kl Ω lj + Ω ik Ω kl S lj IV δ ij /3) + + β 9 τ 4 (Ω ik S kl Ω lm Ω mj Ω ik Ω kl S lm Ω mj ), where τ is turbulent time scale limited from below by Kolmogorov time scale ν/ɛ τ = max Ω ij is rotation rate tensor ( ) k ν ɛ, 6, ɛ =.9kω, (4) ɛ Ω ij = 1 ( ui u ) j, (5) x i and II Ω, IV are invariants formed from S ij, Ω ij given in [] as well as the coefficients β. Both turbulence models use k-ω model equations of similar form Dk Dt = P β kω + [ Dω Dt = γ ω k P βω + + [ (ν + σ ω ν T ) ω ] (ν + σ k ν T ) k + σ d ω ], k ω, (6) where D /Dt / t + (u j )/, and coefficients for model are given in [6]. The efective eddy-viscosity in is given by [] ν T = C µ kτ, C µ = 1 (β 1 + II Ω β 6 ). (7) For the remaining coefficients see also []. 3 Numerical method The system of RANS and turbulence model equations is solved by an artificial compressibility method, which consists of adding a pressure time derivative into the continuity equation: 1 β p t + u i x i =, (8) where the parameter β has been chosen equal to the maximum velocity in the solution domain. The steady solution of the original system is obtained for t. To discretize the equations in space, the cellcentered finite volume method is used. The grid consists of hexahedral cells and is structured and multi-block. An example of handling circular nozzle in the case of jet in cross-flow is shown in Fig. 1. The convective terms are discretized using 3rd order van Leer upwind interpolation (without a limiter). The viscous terms are approximated using nd order central scheme, with cell face derivatives computed on a dual grid of octahedrons constructed over each face of primary grid. Pressure stabilization in the form of pressure Laplacian is added to the continuity equation to prevent pressure-velocity decoupling. The time integration uses backward Euler (implicit) scheme. The Newton linearized implicit operator contains also convective terms, however has the 1st order upwind stencil. In, linear part of turbulent diffusion corresponding to eddy viscosity (7) is discretized implicitely, the remaining part explicitely. Only the negative part of source terms in the turbulence model is discretized implicitly. The resulting block 7-diagonal system of algebraic equations is solved using block relaxation method with direct tri-diagonal solver for slected family of grid lines [3]. The CFL numbers range between 3 and 1. 4 Jet in cross flow This section deals with the simulation of flowfield generated by a round jet injected into a cross-flow of same fluid. The case modeled is given by the available measurement [1], where the speeds were 4.8 m/s and 1. m/s in the jet nozzle and in the cross-flow, respectively, with
3 Fig. 1: Detail of the grid in plane (x, y) z x B H y d Uj A Ue C C/ Fig. : Computational domain for jet in cross-flow no effects of compressibility. The turbulence intensity in both inlets was very low (under.%). The rectangular computational domain approximating aerodynamic tunnel test section is shown in Fig.. The dimensions are (d, A, B, C, H) = (1, 4, 4, 19, 14.8) d, where d = 13.3 mm, and the Reynolds number for jet nozzle diameter and velocity equals 35 (air). The following boundary conditions were prescribed for velocity components u, v, w, kinetic energy of turbulence k and specific dissipation rate ω (in nondimensional form with 4.8 m/s and d as reference velocity and length): x = A (cross flow inlet): u = U e =.5, v = w = k = ω = k/(.1ν) z =, (x + y ).5 (jet nozzle): u = v =, w = 1. ( x + y /.5) 6 k = 3 1 ω = k/(.5ν) x = 4 (outlet): Neumann b.c. for u, v, w, k, ω z = (bottom wall): u = v = w = k = ω = 6ν/(β 1 y 1) 3
4 side and upper walls: inviscid (slip) wall, i.e. normal component of velocity =, Neumann b.c. for: tangential velocity, k, ω The lateral and upper walls were considered inviscid in order to reduce the number of grid points. This simplification is made possible by the fact, that the flow rate is prescribed not in terms of pressure drop in the channel but directly in terms of velocities, thus there is no need to capture friction on walls exactly, and further justified by the fact, that in the measurement measures to eliminate the displacement effect of neglected boundary layers were taken as well (slightly diverging walls). Figures 3 5 show profiles of horizontal and vertical velocity component in the mid-plane. The measured values [1] are for two inlet turbulence intensities and different runs. It is apparent that in the simulation the jet bents too slowly with respect to measurements. A test computation with slightly different heights of the domain showed no improvement. Thus we could say that this is not due to the inviscid forces but rather due to the inaccurate model for turbulent mixing. The model predicts higher vertical velocity than, which is caused by worse capturing of the secondary vortices the jet breaks in, see Fig. 6. Moreover, it seems according to measurements that the flow-field contains large scale unsteady vortices and the errors have more fundamental origin in the steady RANS approach. 5 Impinging jet flow 5.1 Isothermal case Here we consider the solution of RANS equations with and EARS turbulence models in the case of ERCOFTAC s impinging jet [1] for nozzle-to-wall distance nozzle diameters D and nozzle Reynolds number 3. In the previous works [5, 4] it has been found that eddy viscosity models highly over-predict the thickness of wall jet, not to mention the over-prediction of turbulence kinetic energy around impingement. The improvement brings the Bradshaw hypothesis used in the model as an eddy viscosity limiter. The model however has poor near wall behavior typical for k-ω models, apparent in Fig. 7, where the kink in wall shear stress at radial position r/d is missing. A k-ɛ model is capable of predicting the kink and when the -like eddy viscosity limiter is used, also the wall jet prediction is correct [4]. In the present work we consider the k-ω model only and show the influence of constitutive relation. Indeed, the improves the k-ω s behavior here. The and provide similar results in terms of radial velocity, Fig. 8, the Reynolds shear stress shown in Fig. 1 seems better in the model. The also attenuates the over-prediction of turbulence kinetic energy at impingement typical for eddy viscosity, see Fig. 9, however the measured value is still well lower. 5. Heat transfer In this section we consider the heat transfer from the heated impingement wall. Since the temperature difference of approx. 5 C as well as velocity is small, the energy equation is simplified to a passive scalar transport equation for temperature T : T t + u it = [ ] ν T u x i x i P r x i x it, (9) i where P r is laminar Prandtl number and u it turbulent heat transfer. Three closure models for u it were used. First one is the eddy diffusivity model u it = ν T P r T T x i, (1) where ν T is eddy viscosity and P r T turbulent Prandtl number set equal to.85. Then, with the anisotropy models for Reynolds stress we used two anisotropic models for heat transfer. The generalized gradient-diffusion hypothesis (GGDH) model reads u it = C t τu iu T j, C t =.3, (11) 4
5 z / D u / U j w / U j Fig. 3: Horizontal (left) and vertical (right) component of velocity, center-plane, x/d = 1 where τ is the turbulent time scale. In the fully turbulent region, τ = k/ɛ. The lower limit is given by Kolmogorov time scale ν/ɛ, whereas the realizability gives the upper limit. The final expression is [8] [ ( ) ] k ν τ = min max ɛ, 6 C R,, ɛ 6C µ II S C R =.46, C µ =.9. (1) The third model considered, which formally includes the previous two, is the one of Abe, Suga (AS) [9]. u it = c t τ(σ ij + α ij ) T, σ ij = c σ δ ij + c σ1 u iu j/k + + c σ u iu l u l u j/k, α ij = c α τω ij + + c α1 τ(ω il u l u j/k + Ω lj u iu l /k), (13) where the coefficients are.38 c t = (1 exp( R T /1)), 1/4 c σ =, c σ1 =.f b +.1f P r, c σ = 1 f b f P r, c α =, c α1 = (1 f P r ) [.5gA ξ 1 + 5ξ. ] exp[ (ξ/.) ], g A + (ξ +.) f b = (1 f P r ) exp[ (ξ/.) (g A /.3) ], f P r = (1 + (P r/.85) 1.5 ) 1, g A =.3(1 exp[ (R T /7) ]), ξ = τ S ij S ij, (14) and the turbulent time scale τ is taken same as in the, Eq. (4). 5.3 Heat transfer results The temperature of impingement wall is kept constant at T w = 313 K and the temperature of fluid at nozzle exit was T = 93 K. The Nusselt number is defined as Nu = D T T n w, (15) T w T 5
6 z / D u / U j w / U j Fig. 4: Horizontal (left) and vertical (right) component of velocity, center-plane, x/d = z / D u / U j w / U j Fig. 5: Horizontal (left) and vertical (right) component of velocity, center-plane, x/d = 6
7 Fig. 6: Secondary flow in the plane x = d, model (left) and (right) τ w / (ρ w b ) r / d Fig. 7: Wall shear stress on impingement wall where n denotes wall normal direction. The result in terms of Nusselt number is shown in Fig. 11. All models over-predict the heat transfer in the stagnation region. This allows us a statement that even a sofisticated constitutive relation for Reynolds stress as well as for turbulent heat transfer hardly guarantees quantitatively acceptable results. Instead, the turbulent time scale prediction is of higher importance. We recall here a very good results of Merci, Dick [7] achieved with eddy diffusivity heat flux model, but with a non-linear k-ɛ model, where attention has been payed both to the constitutive relation as well as to the ɛ-equation (i.e. time scale prediction). Thus we can conclude that a more complex constitutive relations remain here limited by turbulent scale prediction, here the ω-equation. As for constitutive relations, the eddy diffusivity assumption leads to even qualitatively incorrect result with no drop of heat transfer at r/d 1.5, analogous error as in wall friction. Other models predict the presence of the drop correctly, with GGDH even being in good agreement with experiment starting from this radial position. A simple mean to improve the time scale prediction here could be the well known Yap correction, which (originally in a k-ɛ model) increases 7
8 u / w b u / w b r/d =.5 r/d = u / w b u / w b r/d =.5 r/d = Fig. 8: Velocity in the impinging jet at different radial positions.5. k / w b Fig. 9: Turbulence kinetic energy k on impinging jet axis 8
9 .1.8 τ 13 / w b τ 13 / w b r/d =.5 r/d = τ 13 / w b..4.6 τ 13 / w b r/d =.5 r/d = Fig. 1: Reynolds shear stress in the impinging jet at different radial positions 35 3 ERCOFTAC + eddy diffusivity + GGDH + AS 5 Nu r / D Fig. 11: Nusselt number on the impingement wall 9
10 the source in the ɛ-equation proportionally to k and yw 3 and thus diminishing the near wall overprediction of k in the impingement region. We added the analogous term Y c in the ω-equation of (6) in the following form with the switching functions from [7] Dω Dt = γ ω k P C µ k Y cɛ, S 1.5Ω, Y cɛ =.13(1 f Ry ) k max (.4k 3/ 1, ), otherwise y ɛy S = S ij S ij, Ω = Ω ij Ω ij, f Ry = [ π sin ( { } ) ] Ry min max 5 3, 1, 1, ɛ = β kω, β C µ =.9. The influence of Yap correction in the and GGDH model on Nusselt number is shown in Fig. 1. While it is effective the improvement is insufficient. Same degree of change is observed also for the AS turbulent heat flux model. 6 Conclusions Several turbulence models were used to simulate turbulent momentum and heat transport in jet in cross-flow and in impinging jet flow. A good qualitative agreement with experimental data was achieved in both cases. While gives much better results than model for both flow cases, its accuracy remains limited by the turbulent time scale model which is very similar to that of model. This limitation is most apparent in the impinging jet heat transfer simulation. Acknowledgment This work was partially supported by grant No. 1/5/681 GA CR, and Research Plan MSM References [1] D. Cooper, D. C. Jackson, B. E. Launder, G. C. Liao, J. W. Baughn, X. Yan, and M. Masbah. Normally impinging jet from a circular nozzle. Classic ER- COFTAC Database at University of Surrey. [] Antti Hellsten. New advanced k-ω turbulence model for high-lift aerodynamics. In 4nd AIAA Aerospace Sciences Meeting. Reno, Nevada, 4. [3] K. Kozel, P. Louda, and J. Příhoda. Numerical solution of D and 3D incompressible viscous flow problems. In Internal Flows, 5th International Symposium on Experimental and Computational Aerothermodynamics of Internal Flows, Gdansk, 1. [4] K. Kozel, P. Louda, and J. Příhoda. Numerical solution of turbulent backward facing step and impinging jet flows. In T. Lajos and J. Vad, editors, Conference on modelling fluid flow, pages , Budapest, 3. [5] P. Louda. Numerické řešení dvourozměrného a třírozměrného turbulentního impaktního proudění. PhD thesis, FS ČVUT, Praha,. [6] F. R. Menter. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA Journal, 3(8): , August [7] Bart Merci and Erik Dick. Heat transfer predictions with a cubic kɛ model for axisymmetric turbulent jets impinging onto a flat plate. Int. J. of Heat nad Mass Transfer, 46:469 48, 3. [8] Masoud Rokni and Thomas B. Gatski. Predicting turbulent convective heat transfer in fully developed duct flows. Int. J. of Heat and Fluid Flow, :381 39, 1. 1
11 35 3 ERCOFTAC + GGDH + GGDH + Yap cor. 5 Nu r / D Fig. 1: Influence of Yap correction on Nusselt number [9] K. Suga and K. Abe. Nonlinear eddy viscosity modelling for turbulence and heat transfer near wall and shear-free boundaries. Int. J. of Heat and Fluid Flow, 1:37 48,. [1] V. Uruba, O. Mazur, and P. Jonáš. The instabilities in jet-cross flow interaction. In International Conference Engineering Mechanics, volume IV, pages , Svratka,. [11] S. Wallin. Engineering turbulence modeling for CFD with a focus on explicit algebraic Reynolds stress models. PhD thesis, Royal Institute of Technology, Stockholm,. 11
ON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
Conference Applications of Mathematics 212 in honor of the 6th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 212 ON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
More informationENGINEERING MECHANICS 2012 pp Svratka, Czech Republic, May 14 17, 2012 Paper #195
. 18 m 2012 th International Conference ENGINEERING MECHANICS 2012 pp. 309 315 Svratka, Czech Republic, May 14 17, 2012 Paper #195 NUMERICAL SIMULATION OF TRANSITIONAL FLOWS WITH LAMINAR KINETIC ENERGY
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 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 TRANSITIONAL FLOWS WITH LAMINAR KINETIC ENERGY
Engineering MECHANICS, Vol. 20, 2013, No. 5, p. 379 388 379 NUMERICAL SIMULATION OF TRANSITIONAL FLOWS WITH LAMINAR KINETIC ENERGY JiříFürst* The article deals with the numerical solution of transitional
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 informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
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 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 informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationSTEADY AND UNSTEADY 2D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS FLOW. Radka Keslerová, Karel Kozel
Conference Applications of Mathematics 1 in honor of the th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 1 STEADY AND UNSTEADY D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS
More informationTable of Contents. Foreword... xiii. Preface... xv
Table of Contents Foreword.... xiii Preface... xv Chapter 1. Fundamental Equations, Dimensionless Numbers... 1 1.1. Fundamental equations... 1 1.1.1. Local equations... 1 1.1.2. Integral conservation equations...
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 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 informationMODELLING TURBULENT HEAT FLUXES USING THE ELLIPTIC BLENDING APPROACH FOR NATURAL CONVECTION
MODELLING TURBULENT HEAT FLUXES USING THE ELLIPTIC BLENDING APPROACH FOR NATURAL CONVECTION F. Dehoux Fluid Mechanics, Power generation and Environment Department MFEE Dept.) EDF R&D Chatou, France frederic.dehoux@edf.fr
More informationρ t + (ρu j ) = 0 (2.1) x j +U j = 0 (2.3) ρ +ρ U j ρ
Chapter 2 Mathematical Models The following sections present the equations which are used in the numerical simulations documented in this thesis. For clarity, equations have been presented in Cartesian
More informationMATHEMATICAL MODELING AND NUMERICAL SOLUTION OF 3D ATMOSPHERIC BOUNDARY LAYER
, Vol, Pt, Special Issue Proceedings of International Conference RDAMM 585 MATHEMATICAL MODELING AND NUMERICAL SOLUTION OF D ATMOSPHERIC BOUNDARY LAYER L. Beneš, K. Kozel Department of Technical Mathematics,
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 informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
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 informationSTRESS TRANSPORT MODELLING 2
STRESS TRANSPORT MODELLING 2 T.J. Craft Department of Mechanical, Aerospace & Manufacturing Engineering UMIST, Manchester, UK STRESS TRANSPORT MODELLING 2 p.1 Introduction In the previous lecture we introduced
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 informationfluid mechanics as a prominent discipline of application for numerical
1. fluid mechanics as a prominent discipline of application for numerical simulations: experimental fluid mechanics: wind tunnel studies, laser Doppler anemometry, hot wire techniques,... theoretical fluid
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 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 informationChapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
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 informationModelling of turbulent flows: RANS and LES
Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationFlow characteristics of curved ducts
Applied and Computational Mechanics 1 (007) 55-64 Flow characteristics of curved ducts P. Rudolf a *, M. Desová a a Faculty of Mechanical Engineering, Brno University of Technology,Technická,616 69 Brno,
More informationFLOW OVER A GIVEN PROFILE IN A CHANNEL WITH DYNAMICAL EFFECTS
Proceedings of Czech Japanese Seminar in Applied Mathematics 2004 August 4-7, 2004, Czech Technical University in Prague http://geraldine.fjfi.cvut.cz pp. 63 72 FLOW OVER A GIVEN PROFILE IN A CHANNEL WITH
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 informationHeat Transfer from An Impingement Jet onto A Heated Half-Prolate Spheroid Attached to A Heated Flat Plate
1 nd International Conference on Environment and Industrial Innovation IPCBEE vol.35 (1) (1) IACSIT Press, Singapore Heat Transfer from An Impingement Jet onto A Heated Half-Prolate Spheroid Attached to
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 informationCFD ANALYSIS OF TURBULENT THERMAL MIXING OF HOT AND COLD AIR IN AUTOMOBILE HVAC UNIT
ISTP-6, 005, PRAGUE 6 TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA CFD ANALYSIS OF TURBULENT THERMAL MIING OF HOT AND COLD AIR IN AUTOMOBILE HVAC UNIT Hideo Asano ((, Kazuhiko Suga (3, Masafumi Hirota
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
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 informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationNumerical Simulation of Newtonian and Non-Newtonian Flows in Bypass
Numerical Simulation of Newtonian and Non-Newtonian Flows in Bypass Vladimír Prokop, Karel Kozel Czech Technical University Faculty of Mechanical Engineering Department of Technical Mathematics Vladimír
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 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 informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationFLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes
SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.) INL 1 (3 cr.) 3 sets of home work problems (for 10 p. on written exam) 1 laboration TEN1 (4.5 cr.) 1 written exam
More informationThe current issue and full text archive of this journal is available at
The current issue and full text archive of this journal is available at www.emeraldinsight.com/0961-5539.htm HFF 16,6 660 Received February 2005 Revised December 2005 Accepted December 2005 3D unsteady
More informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More informationComputational Investigations of High-Speed Dual-Stream Jets
9th AIAA/CEAS Aeroacoustics Conference and Exhibit -4 May 3, Hilton Head, South Carolina AIAA 3-33 Computational Investigations of High-Speed Dual-Stream Jets Nicholas J. Georgiadis * National Aeronautics
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 informationROBUST EDDY VISCOSITY TURBULENCE MODELING WITH ELLIPTIC RELAXATION FOR EXTERNAL BUILDING FLOW ANALYSIS
August 11 13, ROBUST EDDY VISCOSITY TURBULENCE MODELING WITH ELLIPTIC RELAXATION FOR EXTERNAL BUILDING FLOW ANALYSIS Mirza Popovac AIT Austrian Institute of Technology Österreichisches Forschungs- und
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 informationA Discussion of Low Reynolds Number Flow for the Two-Dimensional Benchmark Test Case
A Discussion of Low Reynolds Number Flow for the Two-Dimensional Benchmark Test Case M. Weng, P. V. Nielsen and L. Liu Aalborg University Introduction. The use of CFD in ventilation research has arrived
More information2. Conservation Equations for Turbulent Flows
2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationNumerical Simulation of the Transitional Flow on Airfoil
Numerical Simulation of the Transitional Flow on Airfoil Ing. Miroslav Ďuriš Supervisor: Prof. Ing. František Maršík, DrSc. Abstract This paper considers to design and to validate the transitional method
More informationHEAT TRANSFER BY CONVECTION. Dr. Şaziye Balku 1
HEAT TRANSFER BY CONVECTION Dr. Şaziye Balku 1 CONDUCTION Mechanism of heat transfer through a solid or fluid in the absence any fluid motion. CONVECTION Mechanism of heat transfer through a fluid in the
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 4 HEAT TRANSFER IN CHANNEL FLOW BASIC CONCEPTS BASIC CONCEPTS Laminar
More informationThe Simulation of Wraparound Fins Aerodynamic Characteristics
The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:
More informationDetailed Outline, M E 320 Fluid Flow, Spring Semester 2015
Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 I. Introduction (Chapters 1 and 2) A. What is Fluid Mechanics? 1. What is a fluid? 2. What is mechanics? B. Classification of Fluid Flows 1. Viscous
More informationUncertainty quantification for RANS simulation of flow over a wavy wall
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,
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 19 Turbulent Flows Fausto Arpino f.arpino@unicas.it Introduction All the flows encountered in the engineering practice become unstable
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationPhysics of turbulent flow
ECL-MOD 3A & MSc. Physics of turbulent flow Christophe Bailly Université de Lyon, Ecole Centrale de Lyon & LMFA - UMR CNRS 5509 http://acoustique.ec-lyon.fr Outline of the course A short introduction to
More informationOn the Turbulence Modelling for an Air Cavity Interface
1 Introduction On the Turbulence Modelling for an Air Cavity Interface Gem Rotte, Maarten Kerkvliet, and Tom van Terwisga TU Delft, Delft/The Netherlands MARIN, Wageningen/The Netherlands G.M.Rotte@tudelft.nl
More informationComparison of two equations closure turbulence models for the prediction of heat and mass transfer in a mechanically ventilated enclosure
Proceedings of 4 th ICCHMT May 17-0, 005, Paris-Cachan, FRANCE 381 Comparison of two equations closure turbulence models for the prediction of heat and mass transfer in a mechanically ventilated enclosure
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
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 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 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 informationNUMERICAL INVESTIGATION ON THE FLOW CHARACTERISTICS OF A SUPERSONIC JET IMPINGING ON AN AXI-SYMMETRIC DEFLECTOR
ICAS 2002 CONGRESS NUMERICAL INVESTIGATION ON THE FLOW CHARACTERISTICS OF A SUPERSONIC JET IMPINGING ON AN AXI-SYMMETRIC DEFLECTOR S.Sankaran, M.Rajeswara Rao, T.N.V.Satyanarayana, N.Satyanarayana K.Visvanathan
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 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 informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationTurbulence modelling. Sørensen, Niels N. Publication date: Link back to DTU Orbit
Downloaded from orbit.dtu.dk on: Dec 19, 2017 Turbulence modelling Sørensen, Niels N. Publication date: 2010 Link back to DTU Orbit Citation (APA): Sørensen, N. N. (2010). Turbulence modelling. Paper presented
More informationExplicit algebraic Reynolds stress models for internal flows
5. Double Circular Arc (DCA) cascade blade flow, problem statement The second test case deals with a DCA compressor cascade, which is considered a severe challenge for the CFD codes, due to the presence
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationExam in Fluid Mechanics 5C1214
Eam in Fluid Mechanics 5C1214 Final eam in course 5C1214 13/01 2004 09-13 in Q24 Eaminer: Prof. Dan Henningson The point value of each question is given in parenthesis and you need more than 20 points
More informationLaminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition
Sādhanā Vol. 40, Part 2, April 2015, pp. 467 485. c Indian Academy of Sciences Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition RAMBIR BHADOURIYA,
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationCurvature correction and application of the v 2 f turbulence model to tip vortex flows
Center for Turbulence Research Annual Research Briefs 5 157 Curvature correction and application of the v f turbulence model to tip vortex flows By K. Duraisamy AND G. Iaccarino 1. Motivation and objectives
More informationTURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS
HEFAT2014 10 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 16 July 2014 Orlando, Florida TURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS Everts, M.,
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 informationTurbulence Models for Flows with Free Surfaces and Interfaces
AIAA JOURNAL Vol. 44, No. 7, July 2006 Turbulence Models for Flows with Free Surfaces and Interfaces Ebrahim Shirani, Ali Jafari, and Nasser Ashgriz University of Toronto, Toronto, Ontario M5S 3G8, Canada
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
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 informationNumerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel
Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel Arunanshu Chakravarty 1* 1 CTU in Prague, Faculty of Mechanical Engineering, Department of Process Engineering,Technická
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
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 informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationPeriodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain.
Flow Characterization of Inclined Jet in Cross Flow for Thin Film Cooling via Large Eddy Simulation Naqavi, I.Z. 1, Savory, E. 2 and Martinuzzi, R. J. 3 1,2 The Univ. of Western Ontario, Dept. of Mech.
More informationLaminar Flow. Chapter ZERO PRESSURE GRADIENT
Chapter 2 Laminar Flow 2.1 ZERO PRESSRE GRADIENT Problem 2.1.1 Consider a uniform flow of velocity over a flat plate of length L of a fluid of kinematic viscosity ν. Assume that the fluid is incompressible
More informationINTERNATIONAL JOURNAL OF SCIENTIFIC & TECHNOLOGY RESEARCH VOLUME 5, ISSUE 09, SEPTEMBER 2016 ISSN
Numerical Analysis Of Heat Transfer And Fluid Flow Characteristics In Different V-Shaped Roughness Elements On The Absorber Plate Of Solar Air Heater Duct Jitesh Rana, Anshuman Silori, Rohan Ramola Abstract:
More informationComputation of turbulent Prandtl number for mixed convection around a heated cylinder
Center for Turbulence Research Annual Research Briefs 2010 295 Computation of turbulent Prandtl number for mixed convection around a heated cylinder By S. Kang AND G. Iaccarino 1. Motivation and objectives
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 informationPressure-velocity correction method Finite Volume solution of Navier-Stokes equations Exercise: Finish solving the Navier Stokes equations
Today's Lecture 2D grid colocated arrangement staggered arrangement Exercise: Make a Fortran program which solves a system of linear equations using an iterative method SIMPLE algorithm Pressure-velocity
More informationAalto University School of Science and Technology CFD-group/ Department of Applied Mechanics. MEMO No CFD/MECHA DATE: March 15, 2012
Aalto University School of Science and Technology CFD-group/ Department of Applied Mechanics MEMO No CFD/MECHA-2-212 DATE: March 15, 212 TITLE The Effect of Free-Stream Turbulence Parameters on the SST
More informationFLUID MECHANICS. ! Atmosphere, Ocean. ! Aerodynamics. ! Energy conversion. ! Transport of heat/other. ! Numerous industrial processes
SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.)! INL 1 (3 cr.)! 3 sets of home work problems (for 10 p. on written exam)! 1 laboration! TEN1 (4.5 cr.)! 1 written
More information