High Order Discontinuous Galerkin Methods for Aerodynamics
|
|
- Timothy Ferdinand Jefferson
- 5 years ago
- Views:
Transcription
1 High Order Discontinuous Galerkin Methods for Aerodynamics Per-Olof Persson Massachusetts Institute of Technology Collaborators: J. Peraire, J. Bonet, N. C. Nguyen, A. Mohnot Symposium on Recent Developments and Emerging Trends in Computational Design, Analysis and Optimization WPAFB, OH January 30,
2 Outline Background Introduction to High-Order DG Cost and efficiency issues Robustness, RANS and shocks Deformable domains 2
3 Application Areas DNS/LES/DES applications Accurate RANS for engineering applications Drag prediction Rotor dynamics Hypersonics Fluid/Structure Interaction, Flapping flight Computational aeroacoustics Direct simulation using the compressible Navier-Stokes equations Accurate computation of noise sources 3
4 Numerical Schemes Three fundamental requirements: High-order/Low dispersion Unstructured meshes Stability for convervation laws 4
5 Numerical Schemes Three fundamental requirements: FVM High-order/Low dispersion Unstructured meshes Stability for convervation laws 5
6 Numerical Schemes Three fundamental requirements: FVM FDM High-order/Low dispersion Unstructured meshes Stability for convervation laws 6
7 Numerical Schemes Three fundamental requirements: FVM FDM FEM High-order/Low dispersion Unstructured meshes Stability for convervation laws 7
8 Numerical Schemes Three fundamental requirements: FVM FDM FEM DG High-order/Low dispersion Unstructured meshes Stability for convervation laws 8
9 Discontinuous Galerkin: Problems High-order geometry representation and mesh generation High CPU/memory requirements (compared to FVM or high-order FDM) Low tolerance to under-resolved features (no robustness) The challenge is to make DG methods competitive for real-world problems 9
10 Example: Aeroacoustics and K-H Instability Inspired by calculations of Munz et al (2003) (linearized) Nonlinear behavior: Large scale acoustic wave interacts with small scale flow features, leading to vorticity generation p = 7 (8th order accuracy), 100-by-20 square elements 10
11 Example: Flow around Wedge Box Simple but representative geometry Well-studied aeroacoustics test case Low Mach number flow, high order methods required to capture and propagate acoustic waves 11
12 Example: Solid Dynamics Cast governing equations as system of first order conservation laws Account for involutions by appropriate selection of approximating spaces Curl-free deformation tensor Use High-Order DG Lagrangian non-dissipative Neo- Hookean material 12
13 Outline Background Introduction to High-Order DG Cost and efficiency issues Robustness, RANS and shocks Deformable domains 13
14 The Discontinuous Galerkin Method (Reed/Hill 1973, Lesaint/Raviart 1974, Cockburn/Shu 1989-, etc) Consider non-linear hyperbolic system in conservative form: u t + F i (u) = 0 Triangulate domain into simplicial elements Use space of element-wise polynomials of degree p Solve weighted residual formulation Error O(h p+1 ) for smooth problems κ n Reduces to finite volume method for p = 0 u L u R κ 14
15 Nodal Basis and Curved Boundaries Nodal basis within each element sparse element connectivities Isoparametric treatment of curved elements (not only at boundaries) Important research topic: Isotropic mesh Meshing/Node placement for DG W1 wing Anisotropic mesh 15
16 Viscous Discretization Treat the viscous terms by rewriting as a system of first order equations: u t + F i (u) F v (u, σ) = 0 σ u = 0 Discretize using DG, various formulations and numerical fluxes such as: The BR1 scheme (Bassi/Rebay 97): Averaging of fluxes, connections between non-neighboring elements The BR2 scheme (Bassi/Rebay 98): Different lifting operator for each edge, compact connectivities The LDG scheme (Cockburn/Shu 98): Upwind/Downwind, sparse connectivities, some far connections The CDG scheme (Peraire/Persson 06): Modification of LDG for local dependence sparse and compact 16
17 Outline Background Introduction to High-Order DG Cost and efficiency issues Robustness, RANS and shocks Deformable domains 17
18 Sparsity Patterns CDG : LDG : BR2 : and and
19 Implicit DG Memory Requirements Example: Navier-Stokes + S-A in 3-D, p = 4, 10, 000 tetrahedra DOFs = , 000 = 2, 100, 000 Jacobian: matrix entries = 9.6GB Observations: For p 2, 6 tetrahedra cheaper than 1 hex (!) p = 2 and p = 3 optimal for tetrahedra Storage Memory Requirements Fixed resolution r, 3 D M CDG Simplex M CDG Block Polynomial order p 19
20 Matrix Storage and Operations Block matrix representation fundamental for high performance Up to factors of 10 difference between optimized matrix libraries and hand-coded for-loops Compromise between general purpose and application specific format The CDG scheme is well-suited for a compact block storage format With a general format (e.g. compressed column), it is hard to take advantage of the block structure 20
21 Convergence of Iterative Solvers Convergence of various solvers with block-diagonal preconditioner Block Jacobi great when it works, but unreliable Krylov solvers comparable convergence Focus on GMRES(20) Error NACA, t=10 3 Error NACA, t=10 1 Block Jacobi QMR CGS GMRES(20) GMRES # Matrix Vector products # Matrix Vector products 21
22 Preconditioners for Krylov Methods Performance of Krylov methods greatly improved by preconditioning Various standard methods: Block-diagonal: A blockdiag(a) Block-Incomplete LU: A LŨ p-multigrid (low-degree) preconditioner Highly efficient combination (Persson/Peraire 06) Use block-ilu(0) as post-smoother for coarse scale correction In particular, convergence is almost independent of Mach number 22
23 Mach Number Dependence The ILU(0)/Coarse scale preconditioner appears to make convergence almost independent of Mach number Individually, ILU(0) and coarse corrections perform poorly for low Mach Clearly, ILU(0) is an efficient smoother and reduces errors in high frequencies (within and between the elements) # GMRES Iterations NACA, t= Mach Number Block diagonal ILU(0) p0 correction ILU(0)/p0 23
24 Outline Background Introduction to High-Order DG Cost and efficiency issues Robustness, RANS and shocks Deformable domains 24
25 Artificial Viscosity Idea High-order discretization might be unstable when features are unresolved (Gibbs oscillations) Unrealistic to require the mesh to resolve all features (Persson/Peraire 06) Add a small amount of viscosity to original equations so that unresolved features can be represented by approximating space u t + F (u, u) = ( ν u) ε 2ε ε ν 25
26 Resolution SA Turbulence Model 1.0 Velocity Eddy viscosity Eddy Viscosity Profile Shear stress R θ = (from Spalart-Allmaras 1992) Benign at the wall High resolution required at boundary layer edge
27 Resolution SA Turbulence Model Eddy Viscosity Profile: hfill y/δ Coarse Mesh y/δ Fine Mesh 27
28 Resolution SA Turbulence Model Eddy Viscosity Profile: In the transition region ν ν y/ δ * ε ν y (ν + ν) 2 ν y 2 + C ν ε ν ν ε 2 ε ν!!! ( ) ν 2 y Detailed structure of ν in transition region likely to have small effect on overall solution 28
29 Artificial Viscosity vs. Limiting Low order schemes: Width δ s to capture a discontinuity: Modern Finite Volume Schemes: δ s 2h Artificial Viscosity Methods: δ s 4h High order schemes: Limiting: δ s 2h (Not obvious how to exploit increased resolution) Artificial Viscosity Methods: δ s 4h/p (Easy to extend to high order) h h/p 29
30 Example: Transonic Flow over NACA Airfoil NACA 0012, M = 0.8, p = 4 Pressure Shock stabilization by artificial viscosity with non-linear sensor (Persson/Peraire 06) Implicit steady-state solution with Newton solver Entropy Sensor 30
31 Example: Supersonic Flow over NACA Airfoil NACA 0012, M = 1.5, p = 4 Mach Shock stabilization by artificial viscosity with non-linear sensor (Persson/Peraire 06) Implicit steady-state solution with Newton solver Sensor 31
32 Proposed Approach Use Artificial Viscosity with a Non-Linear Sensor Select Parameters based on Available Resolution DG implementation Consistent Treatment of High Order Terms (CDG) Do not exploit nature of Discontinuous Approximation Artificial Viscosity Models: Shocks Viscosity Model: Laplacian form (ν 1 U) Sensor: Density Eddy Viscosity Equation Viscosity Model: Laplacian form (ν 2 ν) Sensor: ν 32
33 Sensor Determine regularity of the solution based on rate of decay of expansion coefficients Periodic Fourier case: f(x) = k= g k e ikx If f(x) has m continuous derivatives g k k (m+1) For simplices: Use orthonormal Koornwinder basis within each element 33
34 M = 0.2, Re = Example: Flat Plate Eddy Viscosity Sensor p = 4 34
35 Example: Flat Plate Velocity Skin friction Lawofthewall Exp. data p =4, grid C p =2, grid B p =1, grid A 6 x Exp. data p =4, grid C p =2, grid B p =1, grid A u + 15 C f y Re x x 10 6 p = 1, 2, 4 (constant h/p) 35
36 Example: Flat Plate Output: Drag on ( x 4.687) p =1 C D error 10 6 p = p = p =4 A B C Grid Convergence rate h p 36
37 Example: RAE2822 M = 0.675, α = 2.31 o, Re = p = 2 (constant h/p) p = 4 C L = C D = C L = C D =
38 Example: RAE2822 M = 0.675, α = 2.31 o, Re = p = 2 (constant h/p) p = 4 C L = C D = C L = C D =
39 Example: RAE2822 C p C f 1.5 p =4 p =2 x p =4 p = C p 0 C f x/c x/c 39
40 M = 0.85, Re = Example: NACA0012 p = 3 (constant h) p = 4 C D = C D =
41 M = 0.85, Re = Example: NACA0012 p = 3 (constant h) p = 4 C D = C D =
42 Example: NACA p =4 p =3 p =2 4 x p =4 p =3 p = C p C f x/c p = 2 C D = p = 3 C D = p = 4 C D = x/c 42
43 Example: NACA0012 Sensor for ν p = 3 p = 4 43
44 Example: NACA0012 Mach Number Shock Sensor p = 4 44
45 Outline Background Introduction to High-Order DG Cost and efficiency issues Robustness, RANS and shocks Deformable domains 45
46 Background and Applications Deformable domains important in many applications: Rotor-stator flows Flapping flight Fluid-structure interaction Several attempts using Arbitrary Lagrangian-Eulerian (ALE) formulations Equations discretized on a deforming grid, time-dependent metric At most third order accuracy in space and time demonstrated Alternative approach [Visbal/Gaitonde 02] in finite difference setting Map from fixed reference domain to real time-varying domain We use a mapping approach, however in a DG setting and with a conservative formulation for free-stream preservation 46
47 ALE Formulation Map from fixed reference domain V to physical deformable domain v(t) A point X in V is mapped to a point x(t) = G(X, t) in v(t) Introduce the mapping deformation gradient G and the mapping velocity v X as nda G = X G v X = G t X NdA G, g, v X x 2 v and set g = det(g) V X 2 x 1 X 1 47
48 Transformed Equations The system of conservation laws in the physical domain v(t) U x t + x F x (U x, x U x ) = 0 x can be written in the reference configuration V as where U X t + X F X (U X, X U X ) = 0 X U X = gu x, F X = gg 1 F x U X G 1 v X and x U x = X (g 1 U X )G T = (g 1 X U X U X X (g 1 ))G T Proof. See paper (only a few lines) 48
49 Geometric Conservation Law A constant solution in v(t) is not necessarily a solution of the discretized equations in V, due to inexact integration of the Jacobian g The time evolution of g is g t X (gg 1 v X ) = 0, X which in general is non-zero Visbal and Gaitonde added corrections to cancel the errors Our approach solves instead the convervative system (ḡg 1 U X ) t X F X = 0 X ḡ t X (gg 1 v X ) = 0 X 49
50 Mappings for Deformable Domains Need mapping x = G(X, t) given the boundary deformations Common approach in ALE settings: Move boundary mesh nodes, update interior nodes by mesh smoothing (Laplacian, spring analogy, etc) Would require interpolation in both space and time to obtain G, v X Stability and accuracy of scheme unclear Can also solve PDEs (usually elliptic) for x in the interior In this work, we use a blending technique which gives explicit expressions for x = G(X, t) Easy to differentiate to obtain G, v X 50
51 Blending of Meshes Our approach forms a smooth weighted combination of two meshes, displaced according to each boundary deformation Original Mesh Rigidly Moved Mesh Blended Mesh 51
52 Example: Euler Vortex Propagate an Euler vortex on a variable domain with x(ξ, η, t) = ξ sin(2πξ/20) sin(πη/7.5) sin(1.0 2πt/t 0 ) y(ξ, η, t) = η sin(2πξ/20) sin(πη/7.5) sin(2.0 2πt/t 0 ) Deformed meshes used for visualization only, everything is computed on the reference mesh 52
53 Example: Euler Vortex, Convergence Optimal order of convergence O(h p+1 ) for mapped scheme p=1 L 2 error 10 4 p=2 p= p=4 p= Element size h Mapped Unmapped 53
54 Example: Free-stream Preservation Solve with initial free-stream solution u 0 (x, y) = u L 2 -error of error u(x, y, T ) u 0 with and without ḡ equation 10 0 p=1 p=1 L 2 error 10 5 p=2 p=2 p=3 p=3 p= p=5 p=4 p= Element size h 54 With Free stream Preservation Without Free stream Preservation
55 Example: Flow around Oscillating Cylinder Oscillating unit circle, y c (t) = A sin(2πft) with A = 4/3,f = 0.1. Reynolds number 400, Mach 0.2, degree p = 4 approximation Compare rigid and blended mappings Rigid Mapping Blended Mapping 55
56 Example: Flow around Oscillating Cylinder Good agreement between different mappings of the same domain deformation Rigid Mapping Blended Mapping Drag Lift Time Drag Lift Time 56
57 Example: Heaving and Pitching Foil in Wake NACA 0012 foil heaving and pitching in wake of D-section cylinder Both oscillate y(t) = A sin(2πft), foil pitching θ = a sin(2πft + π/2) Based on experimental study [Gopalkrishnan et al 94] 57
58 Example: Heaving and Pitching Foil in Wake Thrust/Drag of combined cylinder/foil highly dependent on location/phase Future work: Parametric study to identify different modes of interaction Lift Drag 3 2 Lift (cylinder) Lift (foil) 3 2 Drag (cylinder) Drag (foil) Time Time 58
59 Example: Pitching Airfoil HT13 airfoil attached to translating and heaving point by torsional spring Fluid properties: Re = 5000, M = 0.2 Forced vertical motion r z (t) = r 0 sin ωt (at leading edge) Moment equation: I θ + Cθ S r z (t) + M aero = 0 I moment of inertia, C spring stiffness, S = mx c static unbalance M aero moment from fluid r z (t) θ 59
60 Example: Pitching Airfoil Low Amplitude r z θ M D L Time High Amplitude r z θ 1 M D L Time 60
61 Conclusions Important steps toward a practical DG solver for real-world problems Efficient discretizations and solvers Artificial viscosity approach for producing grid independent RANS solutions and sub-cell shock capturing Effect of grid distortion not fully understood (grid generation?) RANS is much harder than Euler/Laminar NS Optimal accuracy for deformable domains by mapping approach Current work: Extensions for LES/DES, parallelization for large scale 3-D problems, shock capturing for solids, etc 61
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationRANS Solutions Using High Order Discontinuous Galerkin Methods
RANS Solutions Using High Order Discontinuous Galerkin Methods Ngoc Cuong Nguyen, Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 2139, U.S.A. We present a practical
More informationDiscontinuous Galerkin Solution of the Navier-Stokes Equations on Deformable Domains
Discontinuous Galerkin Solution of the Navier-Stokes Equations on Deformable Domains P.-O. Persson a J. Bonet b and J. Peraire c, a Department of Mathematics, University of California, Berkeley, Berkeley,
More informationNEWTON-GMRES PRECONDITIONING FOR DISCONTINUOUS GALERKIN DISCRETIZATIONS OF THE NAVIER-STOKES EQUATIONS
NEWTON-GMRES PRECONDITIONING FOR DISCONTINUOUS GALERKIN DISCRETIZATIONS OF THE NAVIER-STOKES EQUATIONS P.-O. PERSSON AND J. PERAIRE Abstract. We study preconditioners for the iterative solution of the
More informationChapter 1. High-Order Discontinuous Galerkin Methods for CFD
Chapter 1 High-Order Discontinuous Galerkin Methods for CFD Jaime Peraire and Per-Olof Persson 1.1. Introduction In recent years it has become clear that the current computational methods for scientific
More informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 02139, U.S.A. We present an efficient implicit
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 informationThe hybridized DG methods for WS1, WS2, and CS2 test cases
The hybridized DG methods for WS1, WS2, and CS2 test cases P. Fernandez, N.C. Nguyen and J. Peraire Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics, MIT 5th High-Order
More informationA Linear Multigrid Preconditioner for the solution of the Navier-Stokes Equations using a Discontinuous Galerkin Discretization. Laslo Tibor Diosady
A Linear Multigrid Preconditioner for the solution of the Navier-Stokes Equations using a Discontinuous Galerkin Discretization by Laslo Tibor Diosady B.A.Sc., University of Toronto (2005) Submitted to
More informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More information1. Introduction Some Basic Concepts
1. Introduction Some Basic Concepts 1.What is a fluid? A substance that will go on deforming in the presence of a deforming force, however small 2. What Properties Do Fluids Have? Density ( ) Pressure
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationSub-Cell Shock Capturing for Discontinuous Galerkin Methods
Sub-Cell Shock Capturing for Discontinuous Galerkin Methods Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 39, U.S.A. A shock capturing strategy for higher order
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationCompressible Navier-Stokes (Euler) Solver based on Deal.II Library
Compressible Navier-Stokes (Euler) Solver based on Deal.II Library Lei Qiao Northwestern Polytechnical University Xi an, China Texas A&M University College Station, Texas Fifth deal.ii Users and Developers
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 informationAn Optimization-based Framework for Controlling Discretization Error through Anisotropic h-adaptation
An Optimization-based Framework for Controlling Discretization Error through Anisotropic h-adaptation Masayuki Yano and David Darmofal Aerospace Computational Design Laboratory Massachusetts Institute
More informationProgress in Mesh-Adaptive Discontinuous Galerkin Methods for CFD
Progress in Mesh-Adaptive Discontinuous Galerkin Methods for CFD German Aerospace Center Seminar Krzysztof Fidkowski Department of Aerospace Engineering The University of Michigan May 4, 2009 K. Fidkowski
More informationNonlinear Wave Theory for Transport Phenomena
JOSO 2016 March 9-11 2015 Nonlinear Wave Theory for Transport Phenomena ILYA PESHKOV CHLOE, University of Pau, France EVGENIY ROMENSKI Sobolev Institute of Mathematics, Novosibirsk, Russia MICHAEL DUMBSER
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 informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
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 informationMestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 2012/13
Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 212/13 Exam 2ª época, 2 February 213 Name : Time : 8: Number: Duration : 3 hours 1 st Part : No textbooks/notes allowed 2 nd Part :
More informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
HONOM 2011 in Trento DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology DLR Braunschweig 11. April 2011 1 / 35 Research
More information1st Committee Meeting
1st Committee Meeting Eric Liu Committee: David Darmofal, Qiqi Wang, Masayuki Yano Massachusetts Institute of Technology Department of Aeronautics and Astronautics Aerospace Computational Design Laboratory
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
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 informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
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 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 informationMULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544
MULTIGRID CALCULATIONS FOB. CASCADES Antony Jameson and Feng Liu Princeton University, Princeton, NJ 0544 1. Introduction Development of numerical methods for internal flows such as the flow in gas turbines
More informationMultigrid Algorithms for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations
Multigrid Algorithms for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations Khosro Shahbazi,a Dimitri J. Mavriplis b Nicholas K. Burgess b a Division of Applied
More informationPerformance tuning of Newton-GMRES methods for discontinuous Galerkin discretizations of the Navier-Stokes equations
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2685 Performance tuning of Newton-GMRES methods for discontinuous Galerkin
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 informationMasters in Mechanical Engineering Aerodynamics 1 st Semester 2015/16
Masters in Mechanical Engineering Aerodynamics st Semester 05/6 Exam st season, 8 January 06 Name : Time : 8:30 Number: Duration : 3 hours st Part : No textbooks/notes allowed nd Part : Textbooks allowed
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 informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
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 informationImplementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible Flapping Wings within Overture
10 th Symposimum on Overset Composite Grids and Solution Technology, NASA Ames Research Center Moffett Field, California, USA 1 Implementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible
More informationOutput-Based Error Estimation and Mesh Adaptation in Computational Fluid Dynamics: Overview and Recent Results
Output-Based Error Estimation and Mesh Adaptation in Computational Fluid Dynamics: Overview and Recent Results Krzysztof J. Fidkowski University of Michigan David L. Darmofal Massachusetts Institute of
More informationExtension to moving grids
Extension to moving grids P. Lafon 1, F. Crouzet 2 & F. Daude 1 1 LaMSID - UMR EDF/CNRS 2832 2 EDF R&D, AMA April 3, 2008 1 Governing equations Physical coordinates Generalized coordinates Geometrical
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationEdwin van der Weide and Magnus Svärd. I. Background information for the SBP-SAT scheme
Edwin van der Weide and Magnus Svärd I. Background information for the SBP-SAT scheme As is well-known, stability of a numerical scheme is a key property for a robust and accurate numerical solution. Proving
More informationarxiv: v1 [math.na] 25 Jan 2017
Stage-parallel fully implicit Runge-Kutta solvers for discontinuous Galerkin fluid simulations arxiv:1701.07181v1 [math.na] 25 Jan 2017 Will Pazner 1 and Per-Olof Persson 2 1 Division of Applied Mathematics,
More information[N175] Development of Combined CAA-CFD Algorithm for the Efficient Simulation of Aerodynamic Noise Generation and Propagation
The 32nd International Congress and Exposition on Noise Control Engineering Jeju International Convention Center, Seogwipo, Korea, August 25-28, 2003 [N175] Development of Combined CAA-CFD Algorithm for
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 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 informationSolving Large Nonlinear Sparse Systems
Solving Large Nonlinear Sparse Systems Fred W. Wubs and Jonas Thies Computational Mechanics & Numerical Mathematics University of Groningen, the Netherlands f.w.wubs@rug.nl Centre for Interdisciplinary
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 informationFluid-Structure Interaction Problems using SU2 and External Finite-Element Solvers
Fluid-Structure Interaction Problems using SU2 and External Finite-Element Solvers R. Sanchez 1, D. Thomas 2, R. Palacios 1, V. Terrapon 2 1 Department of Aeronautics, Imperial College London 2 Department
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationIntroduction to Turbulence AEEM Why study turbulent flows?
Introduction to Turbulence AEEM 7063-003 Dr. Peter J. Disimile UC-FEST Department of Aerospace Engineering Peter.disimile@uc.edu Intro to Turbulence: C1A Why 1 Most flows encountered in engineering and
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 informationShock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
Abstract We present a shock capturing procedure for high order Discontinuous Galerkin methods, by which shock regions are refined in sub-cells and treated by finite volume techniques Hence, our approach
More informationToward Practical Aerodynamic Design Through Numerical Optimization
Toward Practical Aerodynamic Design Through Numerical Optimization David W. Zingg, and Laura Billing, Institute for Aerospace Studies, University of Toronto 4925 Dufferin St., Toronto, Ontario M3H 5T6,
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 informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
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 informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationAPPLICATION OF SPACE-TIME MAPPING ANALYSIS METHOD TO UNSTEADY NONLINEAR GUST-AIRFOIL INTERACTION PROBLEM
AIAA 2003-3693 APPLICATION OF SPACE-TIME MAPPING ANALYSIS METHOD TO UNSTEADY NONLINEAR GUST-AIRFOIL INTERACTION PROBLEM Vladimir V. Golubev* and Axel Rohde Embry-Riddle Aeronautical University Daytona
More informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
16th International Conference on Finite Elements in Flow Problems DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology
More informationTowards Reduced Order Modeling (ROM) for Gust Simulations
Towards Reduced Order Modeling (ROM) for Gust Simulations S. Görtz, M. Ripepi DLR, Institute of Aerodynamics and Flow Technology, Braunschweig, Germany Deutscher Luft und Raumfahrtkongress 2017 5. 7. September
More informationOutput-based Mesh Adaptation for High Order Navier-Stokes Simulations on Deformable Domains
Output-based Mesh Adaptation for High Order Navier-Stokes Simulations on Deformable Domains Steven M. Kast, Krzysztof J. Fidkowski Department of Aerospace Engineering, University of Michigan, Ann Arbor,
More informationDIRECT NUMERICAL SIMULATIONS OF HIGH SPEED FLOW OVER CAVITY. Abstract
3 rd AFOSR International Conference on DNS/LES (TAICDL), August 5-9 th, 2001, Arlington, Texas. DIRECT NUMERICAL SIMULATIONS OF HIGH SPEED FLOW OVER CAVITY A. HAMED, D. BASU, A. MOHAMED AND K. DAS Department
More informationA High-Order Implicit-Explicit Fluid-Structure Interaction Method for Flapping Flight
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2690 A High-Order Implicit-Explicit Fluid-Structure Interaction Method
More informationAerodynamic Resonance in Transonic Airfoil Flow. J. Nitzsche, R. H. M. Giepman. Institute of Aeroelasticity, German Aerospace Center (DLR), Göttingen
Aerodynamic Resonance in Transonic Airfoil Flow J. Nitzsche, R. H. M. Giepman Institute of Aeroelasticity, German Aerospace Center (DLR), Göttingen Source: A. Šoda, PhD thesis, 2006 Slide 2/39 Introduction
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
More informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
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 informationHp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme Aravind Balan, Michael Woopen and Georg May AICES Graduate School, RWTH Aachen University, Germany 22nd AIAA Computational
More informationImplementation and Comparisons of Parallel Implicit Solvers for Hypersonic Flow Computations on Unstructured Meshes
20th AIAA Computational Fluid Dynamics Conference 27-30 June 2011, Honolulu, Hawaii AIAA 2011-3547 Implementation and Comparisons of Parallel Implicit Solvers for Hypersonic Flow Computations on Unstructured
More informationNonlinear Iterative Solution of the Neutron Transport Equation
Nonlinear Iterative Solution of the Neutron Transport Equation Emiliano Masiello Commissariat à l Energie Atomique de Saclay /DANS//SERMA/LTSD emiliano.masiello@cea.fr 1/37 Outline - motivations and framework
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 informationApplication of High-Order Summation-by-Parts Operators to the Steady Reynolds-Averaged Navier-Stokes Equations. Xiaoyue Shen
Application of High-Order Summation-by-Parts Operators to the Steady Reynolds-Averaged Navier-Stokes Equations Xiaoyue Shen Supervisor: Prof. David W. Zingg A thesis submitted in conformity with the requirements
More informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
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 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 informationA Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics
A Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics Allen R, PhD Candidate Peter J Attar, Assistant Professor University of Oklahoma Aerospace and Mechanical Engineering
More informationA Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization
A Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization Juan J. Alonso Department of Aeronautics & Astronautics Stanford University jjalonso@stanford.edu Lecture 19 AA200b Applied Aerodynamics
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 informationFourth-Order Implicit Runge-Kutta Time Marching Using A Newton-Krylov Algorithm. Sammy Isono
Fourth-Order Implicit Runge-Kutta Time Marching Using A Newton-Krylov Algorithm by Sammy Isono A thesis submitted in conformity with the requirements for the degree of Master of Applied Science Graduate
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationAA210A Fundamentals of Compressible Flow. Chapter 5 -The conservation equations
AA210A Fundamentals of Compressible Flow Chapter 5 -The conservation equations 1 5.1 Leibniz rule for differentiation of integrals Differentiation under the integral sign. According to the fundamental
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 informationFinite volume method on unstructured grids
Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More information41st AIAA Aerospace Sciences Meeting and Exhibit January 6 9, 2003/Reno, NV
AIAA 3 267 Calculation of Unsteady Transonic Flow by an Euler Method with Small Perturbation Boundary Conditions C. Gao, S. Luo, F. Liu Department of Mechanical and Aerospace Engineering University of
More informationAlgebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More information1. Fluid Dynamics Around Airfoils
1. Fluid Dynamics Around Airfoils Two-dimensional flow around a streamlined shape Foces on an airfoil Distribution of pressue coefficient over an airfoil The variation of the lift coefficient with the
More informationAdjoint-based aerodynamic shape optimization
IT Licentiate theses 2003-012 Adjoint-based aerodynamic shape optimization OLIVIER AMOIGNON UPPSALA UNIVERSITY Department of Information Technology Adjoint-based aerodynamic shape optimization BY OLIVIER
More informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
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 informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
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 informationApproximate tensor-product preconditioners for very high order discontinuous Galerkin methods
Approximate tensor-product preconditioners for very high order discontinuous Galerkin methods Will Pazner 1 and Per-Olof Persson 2 1 Division of Applied Mathematics, Brown University, Providence, RI, 02912
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationDerivatives for Time-Spectral Computational Fluid Dynamics using an Automatic Differentiation Adjoint
Derivatives for Time-Spectral Computational Fluid Dynamics using an Automatic Differentiation Adjoint Charles A. Mader University of Toronto Institute for Aerospace Studies Toronto, Ontario, Canada Joaquim
More information