Some remarks on grad-div stabilization of incompressible flow simulations
|
|
- Beverley Black
- 6 years ago
- Views:
Transcription
1 Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis for Singularly Perturbed Problems Dresden University of Technology, November 16-18, 2011 Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
2 Outline 1 Incompressible Navier-Stokes model 2 Numerical analysis of grad-div stabilized Oseen problem 3 Some recent result on limit case γ 4 A (potentially) new approach to parameter design Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
3 Incompressible Navier-Stokes model Navier-Stokes problem Incompressible Navier-Stokes model: Find velocity u, pressure p t u (2νDu) + (u )u + p = f in (0, T] Ω u = 0 in [0, T] Ω u t=0 = u 0 in Ω R d no-slip boundary conditions u = 0 (for simplicity) deformation tensor Du = 1 2 ( u + ( u)t ) viscosity ν (Reynolds number Re = UL ν ). Claude Louis Marie Henri Navier George Gabriel Stokes Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
4 Incompressible Navier-Stokes model Finite element approximation V = [H0(Ω)] 1 3, Q = L0(Ω) 2 := {q L 2 (Ω) : q dx = 0} Ω T h admissible (possibly anisotropic) mesh Ω = K Th K Conforming finite element spaces: V h V, Q h Q Basic Galerkin FE method: find (u h, p h ): [0, T] V h Q h s.t. (v h, q h ) V h Q h ( tu h, v h ) + (2νDu h, Dv h ) + b S(u h, u h, v h ) (p h, v h ) = (f, v h ) (q h, u h ) = 0 with skew-symmetric convective term b S(u, v, w) := 1 [((u )v, w) ((u )w, v)] 2 Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
5 Incompressible Navier-Stokes model Examples of inf-sup stable approximations Inf-sup stable velocity-pressure FE spaces V h Q h V Q β β(h) s.t. inf q h Q h (q h, v h) sup β > 0 v h V h q h 0 v h 0 No additional pressure stabilization required (at least for laminar flows) Taylor-Hood elements: for k N Vh TH Vh TH Q TH h = [P k+1 (T h ) H 1 0(Ω)] d [P k (T h ) C(Ω)] or Q TH h = [Q k+1 (T h ) H 1 0(Ω)] d [Q k (T h ) C(Ω)], k N problems with mass conservation with increasing order k Scott-Vogelius elements: on barycenter refined tetrahedral meshes T h V SV h Q SV h = [P k+1 (T h ) H 1 0(Ω)] d [P disc k (T h ) L 2 0(Ω)], for k d 1 with property [P k+1 (T h )] d P disc k (T h ) strong (pointwise) conservation of mass Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
6 Incompressible Navier-Stokes model Grad-div stabilized Galerkin methods Galerkin FE method with grad-div stabilization: find (u h, p h ): [0, T] V h Q h s.t. (v h, q h ) V h Q h ( tu h, v h ) + (2νDu h, Dv h ) + b S(u h, u h, v h ) (p h, v h ) + (q h, u h ) +(γ u h, v h ) = (f, v h ) classical augmented Lagrangian approach, see e.g. FORTIN/GLOWINSKI [1983] also applied to Maxwell problem introduced by HUGHES/FRANCA [1986] for equal-order interpolation with γ 0(h) Numerical analysis for inf-sup stable interpolation by GELHARD ET AL. [2005] and OLSHANSKII ET AL. [2009]: Choice γ 0(h) is not correct in general case! Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
7 Numerical analysis of grad-div stabilized Oseen problem Grad-div stabilized Oseen model Oseen model: For given b with b = 0, find velocity u, pressure p ν u + (b )u + p + σu = f u = 0 in Ω R d C.W. Oseen Grad-div stabilized Oseen problem: Find (u h, p h ) V h Q h V Q s.t. (v, q) V h Q h : a γ(u h, p h ; v h, q h ) := (ν u h, v h ) + b S(b, u h, v h ) + (σ u h, v h ) (p, v) + (q, u) + K γ K( u h, v h ) K = (f, v h ) Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
8 Numerical analysis of grad-div stabilized Oseen problem Numerical analysis of grad-div stabilized Oseen problem: M. Olshanskii, G. Lube, T. Heister, J. Löwe: Grad-div stabilization and subgrid pressure models for the incompressible Navier-Stokes equations, CMAME 198 (2009), pp Well-posedness: a γ(v h, q h ; v h, q h ) 1 2 [v h, q h ] 2 b 1 2 A-priori estimate: [u u h, p p h ] 2 b K T h h 2k K [( ( ν v h γ v h ν + γ K + h2 K b 2 L (K) ν Equilibration of error terms: leads to dynamic parameter version ( ) p H γ K max 0; k (K) ν u H k+1 (K) ) c p q h 2 0 ν + γ max + ν 1 b 2 L (Ω) ] u 2 H k+1 (K) + 1 p 2 H ν + γ k (K) K ) Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
9 Numerical analysis of grad-div stabilized Oseen problem Example 1: Vortex pairs Oseen problem on Ω = (0, 1) 2 with ν = 10 6, σ = 0 and b = u pairs of vortices strong variation of mesh Reynolds number Re K := u,k h K ν [0, h ν ] H 1 - and L 2 -errors vs. scaling parameter γ 0 of grad-div stabilization for Example 1 with σ = 0, h 1 64 Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
10 Numerical analysis of grad-div stabilized Oseen problem Example 1 Dynamic vs. constant parameter design H 1 - and L 2 -errors vs. scaling parameter γ 0 of dynamic grad-div stabilization for Example 1 H 1 - and L 2 -errors vs. scaling constant parameter γ 0 of grad-div stabilization for Example 1 with ν = 10 6, σ = 0 Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
11 Numerical analysis of grad-div stabilized Oseen problem Example 2: Vortex in boundary layer BERRONE [2001] Oseen problem on Ω = (0, 1) 2 with b = u counter-clockwise vortex in boundary layer ν-dependent solution with u 0 ν 0.35 and p 0 ν Errors in H 1 -seminorm and L 2 -norm vs. scaling parameter γ 0 of grad-div stabilization for Example 2 with ν = 10 4 Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
12 Numerical analysis of grad-div stabilized Oseen problem Example 3: Beltrami flow EITHIER ET AL. [1994] Time-dependent Navier-Stokes flow in Ω = ( 1, 1) 3 with ν = 10 6 Series of counter-rotating vortices intersecting one another at oblique angles Diagonally implicit Runge-Kutta method of order 2 with time step t = 1 64 L 2 -error vs. t [0, 1] without stabilization for different h (left) and with grad-div stabilization for fixed h for ν = 10 6 Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
13 Numerical analysis of grad-div stabilized Oseen problem Example 3 (continued) L 2 (Ω)-error (as function of t) for different values of Re = 1 for the Galerkin scheme, i.e. ν without grad-div stabilization (left) and with grad-div stabilization (right) Obvious improvement even with (time-independent) grad-div stabilization γ K = γ 0 Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
14 Numerical analysis of grad-div stabilized Oseen problem Some conclusions Constant value of grad-div parameter γ gives very often improvements of mass conservation and of other relevant norms ( ) p Dynamic design of γ K max 0; H k (K) ν is not feasible u H k+1 (K) Grad-div stabilization is not necessary in case of e.g. in shear flows! p H k (K) ν u H k+1 (K), Poiseuille-type flow: No improvement with grad-div stabilization! Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
15 Some recent result on limit case γ Stationary Navier-Stokes problem on SV-stable meshes M. Case, V. Ervin, A. Lincke, L. Rebholz: A connection between Scott-Vogelius and grad-div stabilized Taylor-Hood FE approximations of the Navier-Stokes equations, SINUM 49 (2011) 4, pp Stationary case: Let T h s.t. SV-elements are inf-sup stable, e.g. barycenter refined meshes. Theorem: CASE ET AL. Theorem mesh and barycenter refined mesh For any sequence (u h ) γi of TH-elements solutions with constant grad-div stabilization, there is a subsequence which converges to a SV-solution as the grad-div parameter γ i, i. The sequences of TH-modified pressure solutions (p h γ i u h ) γi converges to the corresponding SV pressure. Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
16 Some recent result on limit case γ Nonstationary Navier-Stokes on SV-stable meshes Time-dependent Navier-Stokes flow Remark: A similar result is valid for the time-dependent Navier-Stokes problem with Crank-Nicolson semidiscretization in time. Result from CASE ET AL. SINUM 49 (2011), Fig. 3: Flow around a 2D-cylinder with time-dependent inflow, s.t. Re [0, 100] Behavior of u n u 0 in time Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
17 Some recent result on limit case γ Stationary Navier-Stokes problem on regular meshes Consider the following steady state problem on regular meshes T h : ν( z h, v h ) + ((z h )z h, v h ) = (f, v h ) v h Vh 0 (r h, v h ) = ν( z h, v h ) + ((z h )z h, v h ) (f, v h ) v h (Vh TH ) with V 0 h := {v h V h : v h K = 0 K T h } V TH h := {v h V h : ( v h, q h ) = 0 q h Q TH h } Theorem: CASE ET AL., Theorem 4.1 For any sequence (u h, p h ) γi of TH-elements solutions with grad-div stabilization parameters γ i, a subsequence of (u h, p h γ iρ h ) γi with (ρ h,i, v h ) := ( u h,i, v h ) v h (Vh TH ) converges to (z h, r h ) as the grad-div parameter γ i, i. Remark: A similar result is valid for the time-dependent Navier-Stokes problem with Crank-Nicolson semidiscretization in time. Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
18 Some recent result on limit case γ Some further conclusions Easy implementation of grad-div stabilization into existing code with Taylor-Hood elements! Moderate grad-div stabilization may lead to considerable improvement of mass conservation without deterioration of other solution properties! Besides interesting results for constant γ in CASE ET AL. [2011] results for 3D driven cavity show a deterioration of solution properties for γ on regular meshes Optimal choice of γ:... the search of an optimal γ as a trade-off between mass conservation and energy balance in the FE system... OLSHANSKII ET AL. [2009] Potential remedy: Adaptive choice of γ K based on information from velocity gradients, see next chapter Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
19 A (potentially) new approach to parameter design Eddy-viscosity closure of Reynolds-averaged Navier-Stokes model tu i + 3 j(u iu j) + i(p + τ 3 kk 3 ) ν j 2 u i j=1 ν e := { j=1 3 ju j = 0 j=1 3 (2ν ed ij) = f i, i = 1,..., 3 j=1 c 2 B(u), if u 2 u 2 F := 3 i,j=1 ( iuj)2 0 F 0, if u F = 0 B(u) := β 11β 22 β β 11β 33 β β 22β 33 β 2 23, β ij := 3 mu i mu j m=1 Vreman, A.W.: An eddy-viscosity subgrid-scale model for turbulent shear flows: algebraic theory and application. Phys. Fluids 16 (2004), pp Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
20 A (potentially) new approach to parameter design Some observations: Design of eddy viscosity requires information from velocity gradient (similar to classical Smagorinsky model ν e = (C S ) 2 D(u) F) Data can be easily derived within time discretization and linearization Design rotationally invariant for isotropic filter width Careful algebraic considerations of 320 possible cases of local behavior of flow field u: identifies 13 types of laminar flow (shear-type flow) with zero energy transfer to subgrid scales ν e = 0 Such cases may be identified as coherent structures Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
21 A (potentially) new approach to parameter design Application to low-turbulent channel flow at Re τ = 360 Results taken from VREMAN [2004], Fig. 4: Left: Mean streamwise velocity (normalized by u τ Right: Subgrid eddy viscosity ν e (normalized by u τ H) DNS (circles), Vreman (solid), dynamic Smagorinsky (dashed), standard Smagorinsky (dotted), 0-model (triangles) Important: Vreman s eddy viscosity tends to zero in the viscous sublayer at the wall No van Driest damping in boundary layer necessary! Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
22 A (potentially) new approach to parameter design A potentially new parameter design Application to grad-div stabilization Vreman s eddy viscosity ν e models influence of unresolved velocity fluctuations u = u u Grad-div stabilization may be identified as model of unresolved pressure fluctuations p = p p Application of Vreman s eddy viscosity ν e for parameter γ K (eventually with another scaling) would deactivate grad-div stabilization for shear-type flow To be done: Idea can be transfered to variational multiscale approach with local projection stabilization!! Application of VMS framework of RÖHE/ LUBE CMAME [2010] Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
23 A (potentially) new approach to parameter design Summary. Outlook Improvement of mass conservation with grad-div stabilization of Taylor-Hood elements Further consideration of parameter choice γ Analysis and test of new non-constant parameter design THANKS FOR YOUR ATTENTION! ALL THE BEST TO YOU, DEAR MARTIN! Gert Lube (University of Göttingen) Some remarks on grad-div stabilization Dresden, November 16-18, / 27
Minimal stabilization techniques for incompressible flows
Minimal stabilization tecniques for incompressible flows G. Lube 1, L. Röe 1 and T. Knopp 2 1 Numerical and Applied Matematics Georg-August-University of Göttingen D-37083 Göttingen, Germany 2 German Aerospace
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationSuitability of LPS for Laminar and Turbulent Flow
Suitability of LPS for Laminar and Turbulent Flow Daniel Arndt Helene Dallmann Georg-August-Universität Göttingen Institute for Numerical and Applied Mathematics VMS 2015 10th International Workshop on
More informationWeierstrass Institute for Applied Analysis and Stochastics, Mohrenstr. 39, Berlin, Germany,
Volker John On the numerical simulation of population balance systems Weierstrass Institute for Applied Analysis and Stochastics, Mohrenstr. 39, 10117 Berlin, Germany, Free University of Berlin, Department
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 informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
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 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 informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
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 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 informationProjection Methods for Rotating Flow
Projection Methods for Rotating Flow Daniel Arndt Gert Lube Georg-August-Universität Göttingen Institute for Numerical and Applied Mathematics IACM - ECCOMAS 2014 Computational Modeling of Turbulent and
More informationChapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows
Chapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows Remark 7.1. Turbulent flows. The usually used model for turbulent incompressible flows are the incompressible Navier Stokes equations
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6. Mathematik Preprint Nr. 228 A Variational Multiscale Method for Turbulent Flow Simulation with Adaptive Large Scale
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 informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More 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 informationImplementation of 3D Incompressible N-S Equations. Mikhail Sekachev
Implementation of 3D Incompressible N-S Equations Mikhail Sekachev Navier-Stokes Equations The motion of a viscous incompressible fluid is governed by the Navier-Stokes equations u + t u = ( u ) 0 Quick
More informationAnisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale
More informationMultiscale Computation of Isotropic Homogeneous Turbulent Flow
Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible
More informationA Numerical Investigation of Velocity-Pressure Reduced Order Models for Incompressible Flows
A Numerical Investigation of Velocity-Pressure educed Order Models for Incompressible Flows Alfonso Caiazzo a, Traian Iliescu b,1, Volker John a,c,, Swetlana Schyschlowa a a Weierstrass Institute for Applied
More informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
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 informationLES of turbulent shear flow and pressure driven flow on shallow continental shelves.
LES of turbulent shear flow and pressure driven flow on shallow continental shelves. Guillaume Martinat,CCPO - Old Dominion University Chester Grosch, CCPO - Old Dominion University Ying Xu, Michigan State
More 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 informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationNUMERICAL INVESTIGATION OF BUOY- ANCY DRIVEN FLOWS IN TIGHT LATTICE FUEL BUNDLES
Fifth FreeFem workshop on Generic Solver for PDEs: FreeFem++ and its applications NUMERICAL INVESTIGATION OF BUOY- ANCY DRIVEN FLOWS IN TIGHT LATTICE FUEL BUNDLES Paris, December 12 th, 2013 Giuseppe Pitton
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 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 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 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 information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationA connection between grad-div stabilized FE solutions and pointwise divergence-free FE solutions on general meshes
A connection between grad-div stabilized FE solutions and pointwise divergence-free FE solutions on general meshes Sarah E. Malick Sponsor: Leo G. Rebholz Abstract We prove, for Stokes, Oseen, and Boussinesq
More informationThe Turbulent Rotational Phase Separator
The Turbulent Rotational Phase Separator J.G.M. Kuerten and B.P.M. van Esch Dept. of Mechanical Engineering, Technische Universiteit Eindhoven, The Netherlands j.g.m.kuerten@tue.nl Summary. The Rotational
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 informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationLES AND RANS STUDIES OF OSCILLATING FLOWS OVER FLAT PLATE
LES AND RANS STUDIES OF OSCILLATING FLOWS OVER FLAT PLATE By Chin-Tsau Hsu, 1 Xiyun Lu, and Man-Kim Kwan 3 ABSTRACT: Oscillatory flows over a flat plate are studied by using Large Eddy Simulation (LES)
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 informationReliability of LES in complex applications
Reliability of LES in complex applications Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) DESIDER Symposium Corfu, June 7-8, 27 Sample of complex flow
More informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
More informationVelocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step
Copyright c 2004 Tech Science Press CMC, vol.1, no.3, pp.275-288, 2004 Velocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step B. Wang 1, H.Q. Zhang 1, C.K. Chan 2 and X.L.
More informationSECONDARY MOTION IN TURBULENT FLOWS OVER SUPERHYDROPHOBIC SURFACES
SECONDARY MOTION IN TURBULENT FLOWS OVER SUPERHYDROPHOBIC SURFACES Yosuke Hasegawa Institute of Industrial Science The University of Tokyo Komaba 4-6-1, Meguro-ku, Tokyo 153-8505, Japan ysk@iis.u-tokyo.ac.jp
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 informationTowards pressure-robust mixed methods for the incompressible Navier Stokes equations
Mohrenstrasse 39 10117 Berlin Germany Tel. +49 30 20372 0 www.wias-berlin.de 2016-05-05, Pittsburgh Weierstrass Institute for Applied Analysis and Stochastics Towards pressure-robust mixed methods for
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 informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More informationOn the parameter choice in grad-div stabilization for the Stokes equations
Noname manuscript No. (will be inserted by the editor) On the parameter choice in grad-div stabilization for the Stokes equations Eleanor W. Jenkins Volker John Alexander Linke Leo G. Rebholz Abstract
More informationBOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW
Proceedings of,, BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW Yunho Jang Department of Mechanical and Industrial Engineering University of Massachusetts Amherst, MA 01002 Email:
More information1. Introduction. The flow of an incompressible fluid is governed by the incompressible Navier Stokes equations. u(0, x) =u 0 in Ω,
SIAM J. SCI. COMPUT. Vol. 6, No. 5, pp. 1485 153 c 5 Society for Industrial and Applied Mathematics A FINITE ELEMENT VARIATIONAL MULTISCALE METHOD FOR THE NAVIER STOKES EQUATIONS VOLKER JOHN AND SONGUL
More informationRegularity diagnostics applied to a turbulent boundary layer
Center for Turbulence Research Proceedings of the Summer Program 208 247 Regularity diagnostics applied to a turbulent boundary layer By H. J. Bae, J. D. Gibbon, R. M. Kerr AND A. Lozano-Durán Regularity
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
More informationBoundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity
Boundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity Tobias Knopp D 23. November 28 Reynolds averaged Navier-Stokes equations Consider the RANS equations with
More informationDatabase analysis of errors in large-eddy simulation
PHYSICS OF FLUIDS VOLUME 15, NUMBER 9 SEPTEMBER 2003 Johan Meyers a) Department of Mechanical Engineering, Katholieke Universiteit Leuven, Celestijnenlaan 300, B3001 Leuven, Belgium Bernard J. Geurts b)
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 MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 NUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES Steven
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationTurbulence and its modelling. Outline. Department of Fluid Mechanics, Budapest University of Technology and Economics.
Outline Department of Fluid Mechanics, Budapest University of Technology and Economics October 2009 Outline Outline Definition and Properties of Properties High Re number Disordered, chaotic 3D phenomena
More 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 informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW DATA SET Data provenance: J. Graham 1, M. Lee 2, N. Malaya 2, R.D. Moser 2, G. Eyink 1 & C. Meneveau 1 Database ingest and Web Services: K. Kanov
More informationWhat is Turbulence? Fabian Waleffe. Depts of Mathematics and Engineering Physics University of Wisconsin, Madison
What is Turbulence? Fabian Waleffe Depts of Mathematics and Engineering Physics University of Wisconsin, Madison it s all around,... and inside us! Leonardo da Vinci (c. 1500) River flow, pipe flow, flow
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 informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
More informationINTRODUCTION OBJECTIVES
INTRODUCTION The transport of particles in laminar and turbulent flows has numerous applications in engineering, biological and environmental systems. The deposition of aerosol particles in channels and
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 informationLattice-Boltzmann vs. Navier-Stokes simulation of particulate flows
Lattice-Boltzmann vs. Navier-Stokes simulation of particulate flows Amir Eshghinejadfard, Abouelmagd Abdelsamie, Dominique Thévenin University of Magdeburg, Germany 14th Workshop on Two-Phase Flow Predictions
More informationRegularization modeling of turbulent mixing; sweeping the scales
Regularization modeling of turbulent mixing; sweeping the scales Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) D 2 HFest, July 22-28, 2007 Turbulence
More informationPredicting natural transition using large eddy simulation
Center for Turbulence Research Annual Research Briefs 2011 97 Predicting natural transition using large eddy simulation By T. Sayadi AND P. Moin 1. Motivation and objectives Transition has a big impact
More informationTurbulence: Basic Physics and Engineering Modeling
DEPARTMENT OF ENERGETICS Turbulence: Basic Physics and Engineering Modeling Numerical Heat Transfer Pietro Asinari, PhD Spring 2007, TOP UIC Program: The Master of Science Degree of the University of Illinois
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW AT Re τ = 5200 DATA SET Data provenance: M. Lee 1 & R. D. Moser 1 Database ingest and Web Services: Z. Wu 2, G. Lemson 2, R. Burns 2, A. Szalay
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationElectronic Transactions on Numerical Analysis Volume 32, 2008
Electronic Transactions on Numerical Analysis Volume 32, 2008 Contents 1 On the role of boundary conditions for CIP stabilization of higher order finite elements. Friedhelm Schieweck. We investigate the
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 informationLagrangian acceleration in confined 2d turbulent flow
Lagrangian acceleration in confined 2d turbulent flow Kai Schneider 1 1 Benjamin Kadoch, Wouter Bos & Marie Farge 3 1 CMI, Université Aix-Marseille, France 2 LMFA, Ecole Centrale, Lyon, France 3 LMD, Ecole
More informationCOMPUTATION OF MEAN DRAG FOR BLUFF BODY PROBLEMS USING ADAPTIVE DNS/LES
COMPUTATION OF MEAN DRAG FOR BLUFF BODY PROBLEMS USING ADAPTIVE DNS/LES JOHAN HOFFMAN Abstract. We compute the time average of the drag in two benchmark bluff body problems: a surface mounted cube at Reynolds
More informationSimulations of the turbulent channel flow at Re τ = 180 with projection-based finite element variational multiscale methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Met. Fluids 7; 55:47 49 Publised online 4 Marc 7 in Wiley InterScience (www.interscience.wiley.com). DOI:./fld.46 Simulations of te
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationProblem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 1600
Problem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 6 Overview This problem is aimed at testing the accuracy and the performance of high-order methods on the direct numerical simulation
More informationFinite Element LES and VMS Methods on Tetrahedral Meshes
Finite Element LES and VMS Methods on Tetrahedral Meshes Volker John 1 Adela Kindl 2 and Carina Suciu 3 FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11 50, 66041 Saarbrücken, Germany Abstract
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 informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More informationAnalysis and modelling of subgrid-scale motions in near-wall turbulence
J. Fluid Mech. (1998), vol. 356, pp. 327 352. Printed in the United Kingdom c 1998 Cambridge University Press 327 Analysis and modelling of subgrid-scale motions in near-wall turbulence By CARLOS HÄRTEL
More informationLarge Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena
29 Review Large Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena Masahide Inagaki Abstract Computational Fluid Dynamics (CFD) has been applied
More informationNUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM
Proceedings of ALGORITMY 212 pp. 44 415 NUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM NAVEED AHMED AND GUNAR MATTHIES Abstract. In this paper, we combine
More information43rd AIAA Aerospace Sciences Meeting and Exhibit, Jan 2005, Reno, Nevada
43rd AIAA Aerospace Sciences Meeting and Exhibit, 10-13 Jan 2005, Reno, Nevada A Dynamic Procedure for the Lagrangian Averaged Navier-Stokes-α Model of Turbulent Flows Kamran Mohseni and Hongwu Zhao Aerospace
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
More informationBasic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue
11 Basic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue Yang GUO*, Chisachi KATO** and Yoshinobu YAMADE*** 1 FrontFlow/Blue 1) is a general-purpose finite element program that calculates
More informationPreliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data
35th AIAA Fluid Dynamics Conference, June 6 9, 2005/Toronto,Canada Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data Ellen M. Taylor, M. Pino Martín and Alexander
More informationSystems Theory and Shear Flow Turbulence Linear Multivariable Systems Theory is alive and well
Systems Theory and Shear Flow Turbulence Linear Multivariable Systems Theory is alive and well Dedicated to J. Boyd Pearson Bassam Bamieh Mechanical & Environmental Engineering University of California
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 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 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 informationA fundamental study of the flow past a circular cylinder using Abaqus/CFD
A fundamental study of the flow past a circular cylinder using Abaqus/CFD Masami Sato, and Takaya Kobayashi Mechanical Design & Analysis Corporation Abstract: The latest release of Abaqus version 6.10
More informationModeling of turbulence in stirred vessels using large eddy simulation
Modeling of turbulence in stirred vessels using large eddy simulation André Bakker (presenter), Kumar Dhanasekharan, Ahmad Haidari, and Sung-Eun Kim Fluent Inc. Presented at CHISA 2002 August 25-29, Prague,
More informationFEM techniques for nonlinear fluids
FEM techniques for nonlinear fluids From non-isothermal, pressure and shear dependent viscosity models to viscoelastic flow A. Ouazzi, H. Damanik, S. Turek Institute of Applied Mathematics, LS III, TU
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More information+ = + t x x x x u. The standard Smagorinsky model has been used in the work to provide the closure for the subgridscale eddy viscosity in (2):
International Conference on Methods of Aerophysical Research, ICMAR 008 LARGE EDDY SIMULATION OF TURBULENT ROUND IMPINGING JET B.B. Ilyushin, D.V. Krasinsky Kutateladze Institute of Thermophysics SB RAS
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationThe Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water Interface
Turbulence, Heat and Mass Transfer 4 K. Hanjalic, Y. Nagano and M. Tummers (Editors) 3 Begell House, Inc. The Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water
More informationAER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly
1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence
More information