DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
|
|
- Louisa Lang
- 5 years ago
- Views:
Transcription
1 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 / 35
2 Research group working on discontinuous Galerkin methods for aerodynamic flows at DLR The current group members are: Dr. Ralf Hartmann Tobias Leicht (PhD student) Stefan Schoenawa (PhD student) Marcel Wallraff (PhD student) former group member were: Dr. Joachim Held Florian Prill (PhD student) Numerical results are based on: The DLR-PADGE code which is based on a modified version of deal.ii. 2 / 35
3 Overview Higher-order discontinuous Galerkin methods Error estimation and adaptive mesh refinement for force coefficients Residual-based mesh refinement Numerical results for aerodynamic test cases considered in the EU-project ADIGMA turbulent flow around the 3-element L1T2 high-lift configuration turbulent flow around the DLR-F6 wing-body configuration considered in the EU-project IDIHOM subsonic turbulent flow around the VFE-2 delta wing configuration transonic turbulent flow around the VFE-2 delta wing configuration 3 / 35
4 DG discretization of the RANS-kω equations RANS and Wilcox k-ω turbulence model equations: (F c (u) F v (u, u) ) = S(u, u) Discontinuous Galerkin discretization of order p + 1: Find u h V p h such that R(u h, v h ) R(u h ) v h dx + r(u h ) v + h + ρ(u h) : v + h ds κ T h with the element residual, Ω κ\γ + r Γ (u h ) v + h + ρ (u Γ h) : v + h ds = 0 v h V p h, Γ R(u h ) = S(u h, u h ) F c (u h ) + F v (u h, u h ), and face and boundary residuals r(u h ), ρ(u h ) and r Γ (u h ), ρ Γ (u h ). 4 / 35
5 Error estimation with respect to target quantities Target quantities J(u) of interest are the drag, lift and moment coefficients pressure induced and viscous stress induced parts of the force coefficients 5 / 35
6 Error estimation with respect to target quantities Target quantities J(u) of interest are the drag, lift and moment coefficients pressure induced and viscous stress induced parts of the force coefficients We want to quantity the error of the discrete function u h in terms of a target quantity J( ), i.e. we want to quantity the error Here, J(u) J(u h ) J(u h ) is the computed force coefficient, and J(u) is the exact (but unknown) value of the force coefficient 5 / 35
7 Error estimation for single target quantities Given a discretization: find u h V h,p such that N (u h, v h ) = 0 v h V h,p. and a target quantity J. Using a duality argument we obtain an error representation wrt. J( ): J(u) J(u h ) = R(u h, z) := N (u h, z) R(u h, z h ) = κ η κ. where z h is the solution to the discrete adjoint problem: find z h V h,p such that N [u h ](w h, z h ) = J [u h ](w h ) w h V h,p, and η κ are adjoint-based indicators which are particularly suited for the accurate and efficient approximation of the target quantity J(u). 6 / 35
8 Residual-based mesh refinement The DG discretization: Find u h V p h such that R(u h, v h ) R(u h ) v h dx + r(u h ) v + h + ρ(u h) : v + h ds κ T h Error representation: Residual-based indicators: Ω κ\γ + r Γ (u h ) v + h + ρ (u Γ h) : v + h ds = 0 v h V p h, Γ J(u) J(u h ) = R(u h, z) ( ) 1/2 J(u) J(u h ) (ηκ res ) 2 κ T h ηκ res = h κ R(u h ) κ + hκ 1/2 r(u h ) κ\γ + hκ 1/2 ρ(u h ) κ\γ + hκ 1/2 r Γ (u h ) κ Γ + hκ 1/2 ρ Γ (u h ) κ Γ 7 / 35
9 hp-refinement with anisotropic element subdivision hp-refinement: After having selected an element for refinement, e.g. by residual-based or adjoint-based refinement indicators, decide whether to split the element in subelements, i.e. use h-refinement, when the solution (or the adjoint solution) is smooth/regular increase the polynomial degree, i.e. use p-refinement, when the solution is non-smooth (shocks, sharp trailing edges,... ) The decision is based on the decay of the Legendre series coefficients. 8 / 35
10 hp-refinement with anisotropic element subdivision hp-refinement: After having selected an element for refinement, e.g. by residual-based or adjoint-based refinement indicators, decide whether to split the element in subelements, i.e. use h-refinement, when the solution (or the adjoint solution) is smooth/regular increase the polynomial degree, i.e. use p-refinement, when the solution is non-smooth (shocks, sharp trailing edges,... ) The decision is based on the decay of the Legendre series coefficients. Anisotropic element subdivision: After having selected an element for h-refinement decide upon the specific refinement case based on anisotropic error estimation or on an anisotropic jump indicator: the jump of the discrete solution over element faces is associated with the approximation quality orthogonal to the face 8 / 35
11 The L1T2 high lift configuration full mesh zoom of mesh Coarse mesh of 4740 elements. Grid lines are given by polynomials of degree 4. 9 / 35
12 Turbulent flow around the L1T2 high lift configuration Freestream conditions: M = 0.197, α = and Re = Mach number and streamlines turbulent intensity 10 / 35
13 Turbulent flow around the L1T2 high lift configuration Freestream conditions: M = 0.197, α = and Re = / 35
14 Turbulent flow around the L1T2 high lift configuration Freestream conditions: M = 0.197, α = and Re = / 35
15 Turbulent flow around the L1T2 high lift configuration Freestream conditions: M = 0.197, α = and Re = / 35
16 Turbulent flow around the L1T2 high lift configuration 1.00 global DG(1) DG(1) hp(σ 0 = 1.0) CL C L ref e+05 1e+06 degrees of freedom lift convergence hp-adaptive mesh Tobias Leicht 14 / 35
17 Turbulent flow around the L1T2 high lift configuration 1.00 global DG(1) DG(1) aniso DG(1) hp(σ 0 = 1.0) aniso hp(σ 0 = 1.0) CL C L ref e+05 1e+06 degrees of freedom lift convergence hp-adaptive mesh Tobias Leicht 14 / 35
18 Turbulent flow around the L1T2 high lift configuration 1.00 DG(1) aniso DG(1) hp(σ 0 = 1.0) aniso hp(σ 0 = 1.0) CL C L ref CPU time lift convergence hp-adaptive mesh Tobias Leicht 14 / 35
19 Turbulent flow around the L1T2 high lift configuration 1 isotropic 1e-3 1e-6 1 1e-3 1e-6 1 1e-3 hp isotropic hp anisotropic A suitable solution-adaptive mesh can improve the solver behavior. 1e convergence: residual vs. nonlinear iterations Tobias Leicht 15 / 35
20 The DLR-F6 wing-body configuration without fairing The original mesh of elements has been agglomerated twice. The elements of the coarse mesh of elements are curved based on additional points taken from the original mesh geometry curved mesh with lines given by polynomials of degree 4 16 / 35
21 Subsonic turbulent flow around the DLR-F6 wing-body Modification of the DPW III test case: M = 0.5 (instead of M = 0.75) α = (instead of target lift C l = 0.5) Re = coarse mesh 2 nd order solution DG solutions on coarse mesh of curved elements. 3 rd order solution 4 th order solution 17 / 35
22 Subsonic turbulent flow around the DLR-F6 wing-body TAU on original grid 3 rd order solution after one refinement of coarse mesh 18 / 35
23 Subsonic turbulent flow around the DLR-F6 wing-body TAU on original grid 3 rd order solution after one refinement of coarse mesh 19 / 35
24 Subsonic turbulent flow around the DLR-F6 wing-body Adjoint-based refinement for C d : Mesh after 2 adjoint-based refinement steps Density adjoint 20 / 35
25 Subsonic turbulent flow around the DLR-F6 wing-body Adjoint-based refinement for C d : Convergence of C d (global mesh refinement): p=1, global p=2, global p=3, global TAU 0.04 Mesh after 2 adjoint-based refinement steps Density adjoint e+06 1e+07 degrees of freedom 20 / 35
26 Subsonic turbulent flow around the DLR-F6 wing-body Adjoint-based refinement for C d : Mesh after 2 adjoint-based refinement steps Convergence of C d (global & anisotropic h-refinement): p=1, global p=2, global p=3, global residual-based adjoint-based(cd) adj(cd) + est. TAU Density adjoint e+06 1e+07 degrees of freedom 20 / 35
27 Subsonic turbulent flow around the DLR-F6 wing-body Adjoint-based refinement for C d : Mesh after 2 adjoint-based refinement steps Convergence of C d (global, anisotropic h- & hp-refinement): p=1, global p=1, adjoint-based(cd) p=1, adj(cd) + est. hp, adjoint-based(cd) hp, adj(cd) + est. TAU Density adjoint e+06 1e+07 degrees of freedom 20 / 35
28 The VFE-2 delta wing with medium rounded leading edge The original mesh of elements has been agglomerated twice. The elements of the coarse mesh of elements are curved based on additional points taken from the original mesh geometry original mesh curved coarse mesh with lines with straight lines given by polynomials of degree 4 21 / 35
29 Fully turbulent flow around the VFE-2 delta wing configuration Underlying flow case U.1 in the EU-project IDIHOM The VFE-2 delta wing with medium rounded leading edge at two different flow conditions: U.1b: RANS-kω, subsonic flow at M = 0.4, α = 13.3 and Re = U.1c: RANS-kω, transonic flow at M = 0.8, α = 20.5 and Re = / 35
30 Subsonic flow around the VFE-2 delta wing U.1b: Fully turbulent flow at M = 0.4, α = 13.3 and Re = coarse mesh with curved elements c p distribution 5 th order solution vs. experiment (PSP) 23 / 35
31 Subsonic flow around the VFE-2 delta wing U.1b: Fully turbulent flow at M = 0.4, α = 13.3 and Re = residual-based refined mesh with curved elements c p distribution 4 th order solution vs. experiment (PSP) 24 / 35
32 Subsonic flow around the VFE-2 delta wing U.1b: Fully turbulent flow at M = 0.4, α = 13.3 and Re = th -order solution on residual-based refined mesh with curved elements 25 / 35
33 Subsonic flow around the VFE-2 delta wing U.1b: Fully turbulent flow at M = 0.4, α = 13.3 and Re = th -order solution on residual-based refined mesh with curved elements 2 nd -order solution on residual-based refined mesh with curved elements 25 / 35
34 Fully turbulent flow around the VFE-2 delta wing configuration Underlying flow case U.1 in the EU-project IDIHOM The VFE-2 delta wing with medium rounded leading edge at two different flow conditions: U.1b: RANS-kω, subsonic flow at M = 0.4, α = 13.3 and Re = U.1c: RANS-kω, transonic flow at M = 0.8, α = 20.5 and Re = requires shock capturing 26 / 35
35 Shock-capturing based on artificial viscosity (1) N sc (u h, v) ε(u h ) u h : vdx ε klm (u h ) xl uh m xk v m dx, κ κ κ κ For the compressible Navier-Stokes equations (2nd order DG discretization), 1 ε klm (u h ) = C ε δ kl h 2 β k R m (u h ), k, l = 1,..., d, m = 1,..., n, R m (u h ) = n R q (u h ), m = 1,..., n, q=1 where R(u h ) = (R q (u h ), q = 1,..., n) is the residual of the PDE given by R(u h ) = (F c (u h ) F v (u h, u h )). 1 R. Hartmann. Adaptive discontinuous Galerkin methods with shock-capturing for the compressible Navier-Stokes equations. Int. J. Numer. Meth. Fluids, 51(9 10): , / 35
36 Shock-capturing based on artificial viscosity (2) N sc (u h, v) ε(u h ) u h : vdx ε klm (u h ) xl uh m xk v m dx, κ κ κ κ For the RANS-kω equations (2nd and higher order discretization), 2 ε klm (u h ) = C ε b k b l hκ 2 f p(u h ) Rp(u h) + s p(u + h, u h ) p, b = p p + ε R p(u h ) = Xd+2 m=1 p Xd+2 R m(u h ), s p(u + h, u h ) = u m m=1 p u m s m(u + h, u h ), with R(u h ) = F c (u h ), Z Z s m(u + h, u + h ) v h dx = `H(u h, u h, n ) F c (u + h ) n + m v h ds κ κ 2 F. Bassi et. al. Very high-order accurate Discontinuous Galerkin Computation of transonic turbulent flows on Aeronautical configurations, ADIGMA, NNFMMD 113, / 35
37 Shock-capturing based on artificial viscosity (combines 1 and 2) N sc (u h, v) κ κ ε(u h ) u h : vdx κ κ ε klm (u h ) xl u m h xk v m dx, For the RANS-kω equations (2nd and higher order discretization) ε klm (u h ) = C ε δ kl h2 k f p(u h ) Rp(u h), k, l = 1,..., d, m = 1,..., d + 2, p R p(u h ) = Xd+2 m=1 p u m R m(u h ), R(u h ) = S(u h, u h ) F c (u h ) + F v (u h, u h ), h i = h i /(degree + 1), where h i is the dimension of the element κ in the x i -coordinate direction 29 / 35
38 Transonic flow around the VFE-2 delta wing U.1c: Fully turbulent flow at M = 0.8, α = 20.5 and Re = refined mesh with curved elements c p distribution 4 th order solution vs. experiment (PSP) 30 / 35
39 Transonic flow around the VFE-2 delta wing U.1c: Fully turbulent flow at M = 0.8, α = 20.5 and Re = th -order solution on residual-based refined mesh with curved elements 31 / 35
40 Transonic flow around the VFE-2 delta wing U.1c: Fully turbulent flow at M = 0.8, α = 20.5 and Re = th -order solution on residual-based refined mesh with curved elements 31 / 35
41 Transonic flow around the VFE-2 delta wing U.1c: Fully turbulent flow at M = 0.8, α = 20.5 and Re = th -order solution on residual-based refined mesh with curved elements 31 / 35
42 Transonic flow around the VFE-2 delta wing U.1c: Fully turbulent flow at M = 0.8, α = 20.5 and Re = th -order solution on residual-based refined mesh with curved elements 31 / 35
43 Transonic flow around the VFE-2 delta wing U.1c: Fully turbulent flow at M = 0.8, α = 20.5 and Re = c p 32 / 35
44 Transonic flow around the VFE-2 delta wing U.1c: Fully turbulent flow at M = 0.8, α = 20.5 and Re = c p and M = 1 isosurface 33 / 35
45 Transonic flow around the VFE-2 delta wing U.1c: Fully turbulent flow at M = 0.8, α = 20.5 and Re = c p and M = 1 isosurface zoom of M = 1 isosurface 34 / 35
46 Summary Higher-order discontinuous Galerkin methods Error estimation and adaptive mesh refinement for force coefficients Residual-based mesh refinement Numerical results for aerodynamic flows around the 3-element L1T2 high-lift configuration the DLR-F6 wing-body configuration the VFE-2 delta wing configuration (subsonic and transonic) Computations have been performed with the DLR-PADGE code Thank you 35 / 35
DG 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 informationCurved grid generation and DG computation for the DLR-F11 high lift configuration
ECCOMAS Congress 2016, Crete Island, Greece, 5-10 June 2016 Curved grid generation and DG computation for the DLR-F11 high lift configuration Ralf Hartmann 1, Harlan McMorris 2, Tobias Leicht 1 1 DLR (German
More informationAnisotropic mesh refinement for discontinuous Galerkin methods in two-dimensional aerodynamic flow simulations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2008; 56:2111 2138 Published online 4 September 2007 in Wiley InterScience (www.interscience.wiley.com)..1608 Anisotropic
More informationImplicit 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 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 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 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 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 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 informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationEntropy-Based Drag Error Estimation and Mesh. Adaptation in Two Dimensions
Entropy-Based Drag Error Estimation and Mesh Adaptation in Two Dimensions Krzysztof J. Fidkowski, Marco A. Ceze, and Philip L. Roe University of Michigan, Ann Arbor, MI 48109 This paper presents a method
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 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 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 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 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 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 informationOn Dual-Weighted Residual Error Estimates for p-dependent Discretizations
On Dual-Weighted Residual Error Estimates for p-dependent Discretizations Aerospace Computational Design Laboratory Report TR-11-1 Masayuki Yano and David L. Darmofal September 21, 2011 Abstract This report
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 informationCFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing. Lakhdar Remaki
CFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing Lakhdar Remaki Outline What we doing in CFD? CFD in Industry Shock-filter model for mesh adaptation
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
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 informationCOMPLETE CONFIGURATION AERO-STRUCTURAL OPTIMIZATION USING A COUPLED SENSITIVITY ANALYSIS METHOD
COMPLETE CONFIGURATION AERO-STRUCTURAL OPTIMIZATION USING A COUPLED SENSITIVITY ANALYSIS METHOD Joaquim R. R. A. Martins Juan J. Alonso James J. Reuther Department of Aeronautics and Astronautics Stanford
More informationHigh Order Discontinuous Galerkin Methods for Aerodynamics
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
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 informationSYMMETRIC INTERIOR PENALTY DG METHODS FOR THE COMPRESSIBLE NAVIER STOKES EQUATIONS II: GOAL ORIENTED A POSTERIORI ERROR ESTIMATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 3, Number, Pages 4 62 c 26 Institute for Scientific Computing and Information SYMMETRIC INTERIOR PENALTY DG METHODS FOR THE COMPRESSIBLE
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 informationHigh order methods have the potential of delivering higher accuracy with less CPU time than lower order
21st AIAA Computational Fluid Dynamics Conference June 24-27, 2013, San Diego, CA AIAA 2013-2869 Adjoint Based Anisotropic Mesh Adaptation for the CPR Method Lei Shi and Z.J. Wang Department of Aerospace
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 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 informationFinal abstract for ONERA Taylor-Green DG participation
1st International Workshop On High-Order CFD Methods January 7-8, 2012 at the 50th AIAA Aerospace Sciences Meeting, Nashville, Tennessee Final abstract for ONERA Taylor-Green DG participation JB Chapelier,
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 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 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 informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
More informationDiscontinuous Galerkin methods for compressible flows: higher order accuracy, error estimation and adaptivity
Discontinuous Galerkin methods for compressible flows: higher order accuracy, error estimation and adaptivity Ralf Hartmann Institute of Aerodynamics and Flow Technology German Aerospace Center (DLR) Lilienthalplatz
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 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 informationEfficient A Posteriori Error Estimation for Finite Volume Methods
Efficient A Posteriori Error Estimation for Finite Volume Methods Richard P. Dwight Institute of Aerodynamics and Flow Technology German Aerospace Center (DLR) Braunschweig, D-38108 GERMANY Tel: +49 531
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 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 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 COUPLED-ADJOINT METHOD FOR HIGH-FIDELITY AERO-STRUCTURAL OPTIMIZATION
A COUPLED-ADJOINT METHOD FOR HIGH-FIDELITY AERO-STRUCTURAL OPTIMIZATION Joaquim Rafael Rost Ávila Martins Department of Aeronautics and Astronautics Stanford University Ph.D. Oral Examination, Stanford
More informationAERODYNAMIC SHAPING OF PAYLOAD FAIRING FOR A LAUNCH VEHICLE Irish Angelin S* 1, Senthilkumar S 2
e-issn 2277-2685, p-issn 2320-976 IJESR/May 2014/ Vol-4/Issue-5/295-299 Irish Angelin S et al./ International Journal of Engineering & Science Research AERODYNAMIC SHAPING OF PAYLOAD FAIRING FOR A LAUNCH
More informationSlip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics
Slip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics J.J.W. van der Vegt and H. van der Ven Nationaal Lucht- en Ruimtevaartlaboratorium National
More informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationThe Simulation of Wraparound Fins Aerodynamic Characteristics
The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:
More informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
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 informationNear-wall Reynolds stress modelling for RANS and hybrid RANS/LES methods
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Near-wall Reynolds stress modelling for RANS and hybrid RANS/LES methods Axel Probst (now at: C 2 A 2 S 2 E, DLR Göttingen) René Cécora,
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 informationICES REPORT A DPG method for steady viscous compressible flow
ICES REPORT 13-3 August 013 A DPG method for steady viscous compressible flow by Jesse Chan, Leszek Demkowicz, and Robert Moser The Institute for Computational Engineering and Sciences The University of
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 informationModel reduction of parametrized aerodynamic flows: stability, error control, and empirical quadrature. Masayuki Yano
Model reduction of parametrized aerodynamic flows: stability, error control, and empirical quadrature Masayuki Yano University of Toronto Institute for Aerospace Studies Acknowledgments: Anthony Patera;
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
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 informationAdjoint code development and optimization using automatic differentiation (AD)
Adjoint code development and optimization using automatic differentiation (AD) Praveen. C Computational and Theoretical Fluid Dynamics Division National Aerospace Laboratories Bangalore - 560 037 CTFD
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 informationMultigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations. Todd A.
Multigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations by Todd A. Oliver B.S., Massachusetts Institute of Technology (22) Submitted to the
More informationApplications of adjoint based shape optimization to the design of low drag airplane wings, including wings to support natural laminar flow
Applications of adjoint based shape optimization to the design of low drag airplane wings, including wings to support natural laminar flow Antony Jameson and Kui Ou Aeronautics & Astronautics Department,
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationDetailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition
ILASS-Americas 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 23 Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and
More informationUncertainty Quantification in Viscous Hypersonic Flows using Gradient Information and Surrogate Modeling
Uncertainty Quantification in Viscous Hypersonic Flows using Gradient Information and Surrogate Modeling Brian A. Lockwood, Markus P. Rumpfkeil, Wataru Yamazaki and Dimitri J. Mavriplis Mechanical Engineering
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationLocal Preconditioning for Low Mach-number Aerodynamic Flow Simulation with a Steady Compressible Navier-Stokes Flow Solver
Local Preconditioning for Low Mach-number Aerodynamic Flow Simulation with a Steady Compressible Navier-Stokes Flow Solver Kurt Sermeus and Eric Laurendeau Advanced Aerodynamics Department Bombardier Aerospace
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 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 informationApplication of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Matthew McMullen and Antony Jameson and Juan J. Alonso Dept. of Aeronautics & Astronautics Stanford University
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 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 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 informationA STUDY OF MULTIGRID SMOOTHERS USED IN COMPRESSIBLE CFD BASED ON THE CONVECTION DIFFUSION EQUATION
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More 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 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 informationChapter 9 Flow over Immersed Bodies
57:00 Mechanics of Fluids and Transport Processes Chapter 9 Professor Fred Stern Fall 009 1 Chapter 9 Flow over Immersed Bodies Fluid flows are broadly categorized: 1. Internal flows such as ducts/pipes,
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 informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationCOMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL
COMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL Kengo Asada 1 and Kozo Fujii 2 ABSTRACT The effects of body force distribution on the
More informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More informationRoyal Aeronautical Society 2016 Applied Aerodynamics Conference Tuesday 19 th Thursday 21 st July Science Centre, Bristol, UK
Assessment and validation of aerodynamic performance results for a Natural Laminar Flow transonic wing tested in cryogenic conditions via simulation of turbulent wedges in CFD. Royal Aeronautical Society
More informationLINEARISED FREQUENCY DOMAIN GUST ANALYSIS OF LARGE CIVIL AIRCRAFT
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationADJOINT AND DEFECT ERROR BOUNDING AND CORRECTION FOR FUNCTIONAL ESTIMATES
ADJOINT AND DEFECT ERROR BOUNDING AND CORRECTION FOR FUNCTIONAL ESTIMATES Niles A. Pierce and Michael B. Giles Applied & Computational Mathematics, California Institute of Technology Computing Laboratory,
More informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
More informationOptimizing Runge-Kutta smoothers for unsteady flow problems
Optimizing Runge-Kutta smoothers for unsteady flow problems Philipp Birken 1 November 24, 2011 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany. email:
More informationAnisotropic meshes for PDEs: a posteriori error analysis and mesh adaptivity
ICAM 2013, Heraklion, 17 September 2013 Anisotropic meshes for PDEs: a posteriori error analysis and mesh adaptivity Simona Perotto MOX-Modeling and Scientific Computing Department of Mathematics, Politecnico
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 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 informationLocally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials
Locally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials Harald Klimach, Michael Gaida, Sabine Roller harald.klimach@uni-siegen.de 26th WSSP 2017 Motivation Fluid-Dynamic
More informationAdjoint based multi-objective shape optimization of a transonic airfoil under uncertainties
EngOpt 2016-5 th International Conference on Engineering Optimization Iguassu Falls, Brazil, 19-23 June 2016. Adjoint based multi-objective shape optimization of a transonic airfoil under uncertainties
More informationNumerical Investigation of the Transonic Base Flow of A Generic Rocket Configuration
1 Numerical Investigation of the Transonic Base Flow of A Generic Rocket Configuration A. Henze, C. Glatzer, M. Meinke, W. Schröder Institute of Aerodynamics, RWTH Aachen University, Germany March 21,
More informationAn Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond
45 th Aerospace Sciences Meeting and Exhibit, January 8 11, 2007, Reno, Nevada An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond Antony Jameson Department of Aeronautics
More informationA Study of Transonic Flow and Airfoils. Presented by: Huiliang Lui 30 th April 2007
A Study of Transonic Flow and Airfoils Presented by: Huiliang Lui 3 th April 7 Contents Background Aims Theory Conservation Laws Irrotational Flow Self-Similarity Characteristics Numerical Modeling Conclusion
More informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationDrag (2) Induced Drag Friction Drag Form Drag Wave Drag
Drag () Induced Drag Friction Drag Form Drag Wave Drag Outline Nomenclature and Concepts Farfield Drag Analysis Induced Drag Multiple Lifting Surfaces Zero Lift Drag :Friction and Form Drag Supersonic
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 informationComputation of hypersonic shock boundary layer interaction on a double wedge using a differential Reynolds Stress Model
Computation of hypersonic shock boundary layer interaction on a double wedge using a differential Reynolds Stress Model A. Bosco AICES, Templergraben 55, Aachen, 52056, Germany B. Reinartz CATS, Theaterplatz
More 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 informationNumerical Simulation of Flows in Highly Heterogeneous Porous Media
Numerical Simulation of Flows in Highly Heterogeneous Porous Media R. Lazarov, Y. Efendiev, J. Galvis, K. Shi, J. Willems The Second International Conference on Engineering and Computational Mathematics
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 information