A note on the numerical treatment of the k-epsilon turbulence model Λ
|
|
- Silvester Harrell
- 6 years ago
- Views:
Transcription
1 A note on the numerical treatment of the k-epsilon turbulence model Λ Adrián J. Lew y, Gustavo C. Buscaglia z and Pablo M. Carrica Centro Atómico Bariloche and Instituto Balseiro, 8400 Bariloche, Argentina. To appear in International Journal of Computational Fluid Dynamics Abstract Numerical solution of the equations arising from the» ffl turbulence model has difficulties inherent to nonlinear convection-reaction-diffusion equations with strong reaction terms, resulting in that numerical schemes easily become unstable. We present a formulation that stresses on the robustness of the solution method, tackling common problems that produce instability. The main contribution concerns the loss of positivity of» and ffl, which is addressed by acting on the coefficients of the reaction and diffusion terms rather than on the turbulent variables themselves. In addition, a linearized implicit, non-iterative, treatment of the wall law is proposed. Keywords: turbulence, numerical methods, k-epsilon model, finite elements 1 Introduction The numerical simulation of turbulent flows modeled by the» ffl model has been the subject of much research in the last years. This article is mainly focused on the stability of numerical approximations for this model. Though mathematical results exist ensuring the well-posedness of the equations [7], the strong nonlinearities may interact with discretization errors in such a way as to instabilize computations. This indeed happens if, as frequently occurs, no physically meaningful initial condition is at hand. A typical behavior of unstable computations involves the loss of positivity of» or ffl. This changes the sign of several terms in the equations, with disastrous effects. To avoid such loss of positivity is not an easy task. One may limit» and ffl from below (e.g., [5]) but the solution strongly depends on the limitation strategy and numerical results are frequently unacceptable. The failure of this simple approach motivates the use of numerical techniques that preserve positivity [3, 7], but in most unstructured meshes obtuse angles exist that preclude these algorithms from guaranteing positive results. Though elaborate iterative techniques provide further improvement [3], the difficulty is not removed. Another alternative is to solve for the logarithms of» and ffl [4], but one then deals with a modified set of equations involving exponentials of the unknowns. The numerical trick we propose here is very simple and can be implemented in most» ffl solvers. The idea is to look at the linearized equations for» and ffl as convection-diffusion-reaction equations, and to notice that their coefficients are always positive for the exact solution. Since a negative coefficient must be due to discretization errors, it can be put to zero remorselessly. One never touches the values of» or ffl, just the coefficients of the linearized equations. This trick, together with an improved treatment of the wall-law, results in a very robust scheme. It was combined with a recent equal-order method introduced by Codina and Blasco [2], using linear interpolation for all the unknowns. The following sections describe the proposed method, with numerical tests assessing both its stability and accuracy. 2 Problem Formulation Let ~ U, p, f, ρ and μ` be the Reynolds-averaged velocity, pressure, volumetric force, density and viscosity of the fluid, respectively, and also let ν` = μ`=ρ. We start from the Reynolds Averaged Navier-Stokes (RANS) equations for ΛDedicated to F. Basombrío on his 60th birthday y Present address: Mail Code , Graduate Aeronautical Laboratories, California Institute of Technology, Pasadena, CA 91125, USA z Corresponding author, gustavo@cab.cnea.gov.ar
2 incompressible flow ~ + ρ~ U r ~ U + rp μ`4 ~ U + ρr R = f (1) r ~U = 0 (2) Here R is the Reynolds Stress Tensor for which the» ffl turbulence model proposes R = 2 3»I + ν T r U ~ + r U ~ T» 2 ; ν T = c μ (3) ffl where c μ =0:09,» is the turbulent kinetic energy and ffl is the dissipation of». Replacing R from (3) in the RANS equations we obtain equations identical to the laminar case but with an effective viscosity ν tot = ν` + ν T and an effective pressure P = p + 2». 3 The equations for» and ffl are of the convection-diffusion-reaction type, with non-linear coefficients, as given below in a form suitable for later + U ~ r» r (D» r») +fl»» = + U ~ rffl r (D ffl rffl) +fl ffl ffl = F ffl (5) where the diffusion coefficients are given by D» = ν T ff» + ν`; D ffl = ν T ff ffl + ν` (6) the reaction coefficients by fl» = ffl» ; fl ffl ffl = c 2 (7)» and the source terms by F» = ν T 2 jr~ U + r U ~ T j 2 ; F ffl = c 1» 2 jr~ U + r U ~ T j 2 (8) with c 1 =0:126, c 2 =1:92, ff» =1:0 and ff ffl =1:3. Notice that the diffusion and reaction coefficients, and the source terms, are non-negative for physically admissible solutions (non-negative» and ffl). The previous equations are defined on a domain Ω, that we assume bounded The portion corresponding to solid walls is denoted by w, along which a layer of thickness ffi has been removed and replaced by a wall law. This amounts to the imposition of the no-penetration condition U ~ n = 0 (n is the normal), of a tangential stress fi w (opposite to local velocity) obeying ψ sfi j U ~ j fififi fi! rfi fifi fi 1 fi w χ ln E ffi fi w fi fi ν` ρ fi =0 (9) ρ fi fi and of Dirichlet-type values for» and ffl» = uλ2 p cμ ffl = uλ3 χffi (10) In (9)-(10), E =9for smooth walls, χ =0:41, u Λ = p jfi w j=ρ, and defining ffi + = u Λ ffi=ν`, the inequalities 20» ffi +» 100 must hold. It is customary to specify the value for ffi over w, so that in order to find fi w one must solve the nonlinear equation (9). A second possibility has been successfully applied in this work: To specify the value of ffi + over the boundary with wall law condition. The main advantage of this choice is that given ffi + we can calculate fi w in a direct fashion (no iterations that could fail to converge). As a balancing drawback, now ffi depends on the solution, and can take different values at different points on the boundary. This is however not an issue when ffi is negligible compared to the dimensions of the domain.
3 3 Numerical Scheme The implemented numerical scheme consists of four fractional steps: (1) Calculate the projection of the pressure gradient, (2) Solve continuity and momentum equations, (3) Solve equation for», and (4) Solve equation for ffl. The first two will not be discussed here (see [2]), except for the treatment of the wall law that is deferred to paragraph 3.2. Let K h and E h be the interpolation spaces for» and ffl, respectively, consisting of continuous, piecewise linear functions. The implemented variational formulations are: Fractional step 3: Let R» =»n+1» n + U ~ n+1 r» n+1 + fl 4t»» n+1 F», find» n+1 2 K h such that Fractional step 4: Let R ffl = ffln+1 ffl n 4t (R» ;s h )+ D» r» n+1 ; rs h + X K2T h (R» ;fi» P» (s h )) K = 0 8s h 2 _ K h (11) + ~ U n+1 rffl n+1 + fl ffl ffl n+1 F ffl, find ffl n+1 2 E h such that (R ffl ;t h )+ D ffl rffl n+1 ; rt h + X K2T h (R ffl ;fi ffl P ffl (t h )) K = 0 8t h 2 _ E h (12) R where (u; v) = uv dω and spaces with a dot above indicate homogeneous Dirichlet data. The first two terms in Eq. Ω (11) (and analogously for Eq. (12)) come from the Galerkin formulation. The third one is a stabilization term, involving an algorithmic parameter fi» and a perturbation function P», which introduces upwinding. Specific designs for this term are provided by the SUPG or SGS methods [1]. 3.1 Positivity of the coefficients The numerical scheme as explained so far performs poorly due to its inability to guarantee that» and ffl remain positive. Undershoots develop in boundary and internal layers, and negative values of either» or ffl result in the immediate instabilization of the scheme. However, by inspection of the formulation it is realized that it is not the negative values of» and ffl that really matter, but instead the appearance of negative reaction or diffusion coefficients that cause the exponential growth of the solution. We propose to limit ν T, fl», fl ffl and F ffl from below without touching the nodal values of» or ffl. The rationale for this is that we can accept local regions with negative values of these fields as long as these problematic regions remain localized and do not lead to global instability. The modifications introduced in the numerical scheme are: ffl ν T is limited from below to a small fraction of ν`, and evaluated at the previous time step. ffl The rest of the coefficients are limited from below by zero, 3.2 Implementation of the wall law From (9), fi w can be written as (» n ) 2 ν T = maxfc μ ; 10 4 ν`g (13) ffl fl» = maxf c μ» n ; 0g (14) ν T fl ffl = maxf c 2ffl n ; 0g (15)» n+1 fi F ffl = maxf c 1» n+1 fi fi fir U ~ + r U fifi ~ T 2 ; 0g (16) 2 ~U fi w = g j U ~ j j U ~ j where the function g simply equals ρ [χ j ~ Uj= ln(e ffi + ))] 2 when ffi + is specified, but is an implicit expression (iterations needed) if ffi is specified. All of the following can be done in either case, the latter requiring differentiation of implicit functions. (17)
4 Evaluating fi w at the previous time step leads to instability. We use instead a linearized implicit scheme, i.e., fii w ~U n+1 ο= fi w i U n+1 j U ~ j n ~U n j fi fififi~u = ~ U n where summation over repeated indices holds. Using (17) and some algebra we get fii w ~U n+1 ψ! = g j U ~ n Ui n + U n+1 i j j U ~ U i n U j n U j U ~ n j j n j j U ~ n j U ~ n j ψ U n i U n j U n+1 j j ~ U n j 2 U n i! (18) (19) This formulation turns out to be stable. Furthermore, the terms multiplying ~ U n+1 in (19) lead to a positive contribution to the system matrix. 4 Numerical results for two test problems 4.1 Fully developed flow in a channel This test consists of the flow between parallel plates due to a constant pressure gradient. The following flow conditions were = 0:52;ν` = 10 4 ;H = 1, resulting in a Reynolds number Re H = 120; 000, based on H, the channel height, and U 0, the maximum velocity. Wall-law boundary conditions were imposed, fixing ffi =0:01196H. Several numerical tests were performed, those chosen as representative are listed in Table 1. Initial condition A corresponds to U,» and ffl set to zero at t =0. This initial condition is physically meaningfull. Conversely, initial condition B corresponds to U=100,»=0, ffl=0, which is non-physical and generates a strong initial transient. Case 1 addresses the accuracy of the method, while cases 2 to 4 were used to test its robustness. In cases 1, 2 and 4 a steady state was reached, and case 3 oscillates. The solutions obtained for the steady state, not shown here for brevity, agree with the experimental results by Laufer and the numerical results by Utnes [9]. Also, the solutions obtained with different grids and time steps (cases 1, 2 and 4) practically coincide. The usefulness of the proposed tricks is realized looking at Figs. 1 and 2, where the maxima and minima of» and ffl during the transient are shown. Wild oscillations appear (locally!) in cases 2, 3 and 4. However, in cases 2 and 4 the steady state is reached without problem. Moreover, persistant oscillations such as those in case 3 can be avoided by refining either the temporal or spatial discretizations (leading to cases 2 and 4, respectively). 4.2 Backward facing step (BFS) The geometry of this problem (2D) is shown in Fig.3. The dimensions are H =1, H 1 =2, L 1 =14and L 2 =30. Wall laws (imposing ffi + =30) hold along top and bottom boundaries. Inflow conditions correspond to fully-developed flow. The case discussed here corresponds to Re =44479(based on step height H and maximum inlet velocity U 0 ), with ν` =2:88:10 4 and ρ =1:0, which lead to U 0 =12:81. The initial condition is uniform over Ω, U =1:0, P =0:0,» =0:01, ffl =0:01. Notice that this does not respect mass conservation and is far away from the steady state solution. The grid consists of 6823 nodes and linear triangles (see Fig. 4). With a quite large time step ( t =0:3125) a quasi-steady state was obtained without difficulty in 240 steps. Nevertheless, at a few nodes near corner a the turbulent dissipation ffl kept oscillating. This local phenomenon did not pollute the solution away from this corner, and could be damped by further reduction of the time step (100 additional steps with t =0:003 proved sufficient). Though not presented here, comparison of velocity and Reynolds stress profiles with previous experimental and numerical results [10] is satisfactory (see [6] for details). The obtained reattachment length, frequently used for comparison, was L r =6:59, in good agreement with both an experimental value L r ß 7 and several works reporting that the» ffl model underestimates L r by 10% to 25% [8]. Coming back to the central topics of this article, there are regions of the flow in which» or ffl are negative, even at steady state. It must be recalled that our method allows for the appearance of such regions, and that our speculation is that in physically meaningfull transients or at the steady state they will be sharply localized and will not affect the overall solution. For the present example at steady state (things are much worse during the initial transient) such regions are shown in Fig. 5. It is clear that» or ffl are negative at just a few nodes. Furthermore, this happens where crossflow gradients are steep, and could be avoided by refining the grid or by introducing a discontinuity-capturing operator in the formulation. The minimum values of» and ffl (see Table 2) are indeed small (in modulus) as compared to the maximum values.
5 5 Concluding remarks A robust algorithm to solve the» ffl equations has been presented. Its design is based on the manipulation of the different terms so as to add as many positive terms as possible to the system matrix and to keep source terms also positive. Then, as discretization errors may induce non-physical negative values of diffusion or reaction coefficients, these are limited from below by zero. This approach works much better than the frequently used limiters on the turbulent variables. Along the same design line is the proposed treatment of the wall boundary condition by semi-implicit linearization, an additional trick being the imposition of ffi + instead of ffi, avoiding iterations in the wall-law algorithm that may fail to converge. Remarkably, the proposed algorithm can easily be implemented in most flow solvers, not just finite element ones. Also worth mentioning is the good behavior of the recent stabilized method of Codina and Blasco [2], for which few reported experiments (none of them turbulent) exist. ACKNOWLEDGMENTS: Financial support from Agencia Nacional de Promoción Científica y Tecnológica through grant PICT is acknowledged. GCB also belongs to Consejo Nacional de Investigaciones Científicas y Técnicas. References [1] R. Codina, Comparison of some finite element methods for solving the diffusion-convection-reaction equation, Comp. Meth. Appl. Mech. Eng., 156, (1998). [2] R. Codina and J. Blasco, Stabilized finite element method for the transient Navier-Stokes equations based on a pressure gradient projection, Comp. Meth. Appl. Mech. Eng., 182, (2000). [3] R. Codina and O. Soto, Finite element implementation of two-equation and algebraic stress turbulence models for steady incompressible flows, Int. J. Numer. Meth. in Fluids, 30, (1999). [4] F. Ilinca and D. Pelletier, Positivity preservation and adaptive solution for the» ffl model of turbulence, AIAA Paper No , 35th Aerospace Sciences Meeting & Exhibit, January 1997, Reno, USA. [5] F. Ilinca, D. Pelletier and F. Arnoux-Guisse, An adaptive finite element scheme for turbulent free shear flows, Int. J. Comput. Fluid Dyn., 8, (1997). [6] A. Lew, The Finite Element Method in High Performance Computing Environments, Master s Thesis in Nuclear Engineering, Instituto Balseiro, Argentina. In Spanish. [7] B. Mohammadi and O. Pironneau, Analysis of the» ffl turbulence model, J. Wiley & Sons (1994). [8] S. Tangham and C. Speziale, Turbulent separated flow past a backward-facing step: A critical evaluation of twoequation turbulence model, Rep , Inst. for Computer Appl. in Sci. and Eng., NASA Langley Res. Ctr. (1991). [9] T. Utnes, Two equations (»,ffl) turbulence computations by the use of a finite element model, Int. J. Numer. Meth. in Fluids, 8, (1988). [10] V. Yakhot, S. Tangham, T. Gatski, S. Orszag and C. Speziale, Development of turbulence models for shear flows by a double expansion technique, Rep , Inst. for Computer Appl. in Sci. and Eng., NASA Langley Res. Ctr. (1991).
6 Case Number of nodes (N) t initial condition A B B B Table 1: Details of the numerical tests for the channel flow problem. The grids are nonuniform, refined near the walls. max» min» max ffl min ffl Table 2: Maximum and minimum values of» and ffl at the steady state of the BFS test problem.
7 100 Case 2 Max(K) 10 Case 4 Case 3 K 1 Case time 4 Min(K) Case 2 a) 2 Case 4 0 Case 1-2 Case 3 K time b) Figure 1: Channel flow problem. Evolution of the maximum (a) and minimum (b) of» in the initial transient for each of the cases.
8 1e+06 Max(E) Case Case Case 4 E Case time a) 0 Case 1 Case 2 Case 4 Min(E) -20 Case 3-40 E time b) Figure 2: Channel flow problem. Evolution of the maximum (a) and minimum (b) of ffl in the initial transient for each of the cases.
9 e f H1 L1 Y a H b X L2 H2 d c Figure 3: BFS problem. Geometry and dimensions of the computational domain. Figure 4: BFS problem. Central portion of the grid
10 K< 0 a) e < 0 b) Figure 5: BFS problem. Regions with negative values of» (a) and ffl (b).
AN 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 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 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 informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
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 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 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 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 informationAerodynamic force analysis in high Reynolds number flows by Lamb vector integration
Aerodynamic force analysis in high Reynolds number flows by Lamb vector integration Claudio Marongiu, Renato Tognaccini 2 CIRA, Italian Center for Aerospace Research, Capua (CE), Italy E-mail: c.marongiu@cira.it
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationAn example of the Rvachev function method
arxiv:1603.00320v1 [physics.flu-dyn] 1 Mar 2016 An example of the Rvachev function method Alexander V. Proskurin Altai State University, Altai State Technical University, k210@list.ru Anatoly M. Sagalakov
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 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 informationThe two-dimensional streamline upwind scheme for the convection reaction equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2001; 35: 575 591 The two-dimensional streamline upwind scheme for the convection reaction equation Tony W. H. Sheu*,1
More informationEigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations
NASA/CR-1998-208741 ICASE Report No. 98-51 Eigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations David L. Darmofal Massachusetts Institute of Technology, Cambridge,
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 informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationOn discontinuity capturing methods for convection diffusion equations
On discontinuity capturing methods for convection diffusion equations Volker John 1 and Petr Knobloch 2 1 Universität des Saarlandes, Fachbereich 6.1 Mathematik, Postfach 15 11 50, 66041 Saarbrücken, Germany,
More informationNUMERICAL AND EXPERIMENTAL INVESTIGATIONS OF AIR FLOW AND TEMPERATURE PATTERNS OF A LOW VELOCITY DIFFUSER
NUMERICAL AND EXPERIMENTAL INVESTIGATIONS OF AIR FLOW AND TEMPERATURE PATTERNS OF A LOW VELOCITY DIFFUSER M Cehlin and B Moshfegh Division of Energy and Mechanical Engineering, Department of Technology,
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 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 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 informationMostafa Momen. Project Report Numerical Investigation of Turbulence Models. 2.29: Numerical Fluid Mechanics
2.29: Numerical Fluid Mechanics Project Report Numerical Investigation of Turbulence Models Mostafa Momen May 2015 Massachusetts Institute of Technology 1 Numerical Investigation of Turbulence Models Term
More informationApplication of the immersed boundary method to simulate flows inside and outside the nozzles
Application of the immersed boundary method to simulate flows inside and outside the nozzles E. Noël, A. Berlemont, J. Cousin 1, T. Ménard UMR 6614 - CORIA, Université et INSA de Rouen, France emeline.noel@coria.fr,
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 informationNUMERICAL STUDY OF CONVECTION DIFFUSION PROBLEM IN TWO- DIMENSIONAL SPACE
IJRRAS 5 () November NUMERICAL STUDY OF CONVECTION DIFFUSION PROBLEM IN TWO- DIMENSIONAL SPACE. K. Sharath Babu,* & N. Srinivasacharyulu Research Scholar, Department of Mathematics, National Institute
More informationUtilising high-order direct numerical simulation for transient aeronautics problems
Utilising high-order direct numerical simulation for transient aeronautics problems D. Moxey, J.-E. Lombard, J. Peiró, S. Sherwin Department of Aeronautics, Imperial College London! WCCM 2014, Barcelona,
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 informationA Stabilized Formulation for Incompressible Elasticity Using Linear Displacement and Pressure Interpolations
A Stabilized Formulation for Incompressible Elasticity Using Linear Displacement and Pressure Interpolations M. Chiumenti y, Q. Valverde Λ, C. Agelet de Saracibar y and M. Cervera y ydept. de Resistència
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 informationCalculation of the separation streamlines of barchans and transverse dunes
Calculation of the separation streamlines of barchans and transverse dunes H. J. Herrmann a,b, J. S. Andrade Jr. b, V. Schatz a, G. Sauermann a and E. J. R. Parteli a a ICP, University of Stuttgart, Pfaffenwaldring
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 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 informationVERTICAL TURBULENT BUOYANT HELIUM JET CFD MODELING AND VALIDATION
VERTICAL TURBULENT BUOYANT HELIUM JET CFD MODELING AND VALIDATION Cheng Z, Agranat V.M. and Tchouvelev A.V. A.V.Tchouvelev & Associates, Inc., 659 Spinnaker Circle, Mississauga, Ontario, Canada L5W R Hydrogenics
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 informationLarge Scale Fluid-Structure Interaction by coupling OpenFOAM with external codes
Large Scale Fluid-Structure Interaction by coupling OpenFOAM with external codes Thomas Gallinger * Alexander Kupzok Roland Wüchner Kai-Uwe Bletzinger Lehrstuhl für Statik Technische Universität München
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
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 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 informationProject 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions
Project 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions R. S. Sellers MAE 5440, Computational Fluid Dynamics Utah State University, Department of Mechanical and Aerospace Engineering
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 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 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 informationDirect Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme
19th AIAA Computational Fluid Dynamics - 5 June 009, San Antonio, Texas AIAA 009-3797 19th AIAA Computational Fluid Dynamics Conference, - 5 June 009, San Antonio, Texas Direct Numerical Simulations of
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 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 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 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 informationNumerical investigation of swirl flow inside a supersonic nozzle
Advances in Fluid Mechanics IX 131 Numerical investigation of swirl flow inside a supersonic nozzle E. Eslamian, H. Shirvani & A. Shirvani Faculty of Science and Technology, Anglia Ruskin University, UK
More information^Laboratorio Nacional de Computaqao Cientifica, Rua Lauro Muller 455, Rio de Janeiro, RJ , Brasil
A space-timefiniteelement formulation for the shallow water equations F.L.B. Ribeiro,' A.C. Galeao," L. Landau* "Programa de Engenharia Civil, COPPE / Universidade Federal do Rio de Janeiro, Caixa Postal
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 informationSYMMETRY BREAKING PHENOMENA OF PURELY VISCOUS SHEAR-THINNING FLUID FLOW IN A LOCALLY CONSTRICTED CHANNEL
ISSN 1726-4529 Int j simul model 7 (2008) 4, 186-197 Original scientific paper SYMMETRY BREAKING PHENOMENA OF PURELY VISCOUS SHEAR-THINNING FLUID FLOW IN A LOCALLY CONSTRICTED CHANNEL Ternik, P. University
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 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 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 informationBEM for turbulentfluidflow. L. Skerget & M. Hribersek Faculty of Mechanical Engineering, University ofmaribor, Slovenia.
BEM for turbulentfluidflow L. Skerget & M. Hribersek Faculty of Mechanical Engineering, University ofmaribor, Slovenia. Abstract Boundary-Domain Integral Method is applied for the solution of incompressible
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 informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
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 SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
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 informationA combined application of the integral wall model and the rough wall rescaling-recycling method
AIAA 25-299 A combined application of the integral wall model and the rough wall rescaling-recycling method X.I.A. Yang J. Sadique R. Mittal C. Meneveau Johns Hopkins University, Baltimore, MD, 228, USA
More informationLARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL
LARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL 09/11/2017 Paolo Casco Stephie Edwige Philippe Gilotte Iraj Mortazavi LES and flow control over a 25 ramp model : context 2 Context Validation
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
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 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 informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
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 informationNUMERICAL SIMULATION OF LDI COMBUSTOR WITH DISCRETE-JET SWIRLERS USING RE-STRESS MODEL IN THE KIVA CODE
NUMERICAL SIMULATION OF LDI COMBUSTOR WITH DISCRETE-JET SWIRLERS USING RE-STRESS MODEL IN THE KIVA CODE S. L. Yang, C. Y. Teo, and Y. K. Siow Department of Mechanical Engineering Engineering Mechanics
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More information5.1 2D example 59 Figure 5.1: Parabolic velocity field in a straight two-dimensional pipe. Figure 5.2: Concentration on the input boundary of the pipe. The vertical axis corresponds to r 2 -coordinate,
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationLarge-eddy simulations for wind turbine blade: rotational augmentation and dynamic stall
Large-eddy simulations for wind turbine blade: rotational augmentation and dynamic stall Y. Kim, I.P. Castro, and Z.T. Xie Introduction Wind turbines operate in the atmospheric boundary layer and their
More informationThe effect of geometric parameters on the head loss factor in headers
Fluid Structure Interaction V 355 The effect of geometric parameters on the head loss factor in headers A. Mansourpour & S. Shayamehr Mechanical Engineering Department, Azad University of Karaj, Iran Abstract
More informationDesign and simulation of Open Circuit Blowdown type Wind Tunnel
Design and simulation of Open Circuit Blowdown type Wind Tunnel Sanjeev Kumar Gupta a, V.K.Dwivedi b, Jitendra Kumar Chauhan c, and Rahul Goswami c a Assistant Professor, Department of Mechanical Engineering,
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 informationA Finite-Element based Navier-Stokes Solver for LES
A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,
More information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
More informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationThe Effect of Undulated Wall on Natural Convection Heat Transfer in Vertical Channels
Proceedings of the 3 rd International Conference on Fluid Flow, Heat and Mass Transfer (FFHMT 16) Ottawa, Canada May 2 3, 2016 Paper No. 123 The Effect of Undulated Wall on Natural Convection Heat Transfer
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 informationHEAT TRANSFER IN A RECIRCULATION ZONE AT STEADY-STATE AND OSCILLATING CONDITIONS - THE BACK FACING STEP TEST CASE
HEAT TRANSFER IN A RECIRCULATION ZONE AT STEADY-STATE AND OSCILLATING CONDITIONS - THE BACK FACING STEP TEST CASE A.K. Pozarlik 1, D. Panara, J.B.W. Kok 1, T.H. van der Meer 1 1 Laboratory of Thermal Engineering,
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 informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
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 informationNumerical modeling of complex turbulent flows
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
More informationME 144: Heat Transfer Introduction to Convection. J. M. Meyers
ME 144: Heat Transfer Introduction to Convection Introductory Remarks Convection heat transfer differs from diffusion heat transfer in that a bulk fluid motion is present which augments the overall heat
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 informationSUPG STABILIZATION PARAMETERS CALCULATED FROM THE QUADRATURE-POINT COMPONENTS OF THE ELEMENT-LEVEL MATRICES
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 24 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) W. Rodi and P. Le Quéré (assoc. eds.) Jyväskylä,
More information2.3 The Turbulent Flat Plate Boundary Layer
Canonical Turbulent Flows 19 2.3 The Turbulent Flat Plate Boundary Layer The turbulent flat plate boundary layer (BL) is a particular case of the general class of flows known as boundary layer flows. The
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 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 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 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 informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
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 informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
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 information