Curved grid generation and DG computation for the DLR-F11 high lift configuration
|
|
- Shon Sims
- 5 years ago
- Views:
Transcription
1 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 Aerospace Center), Institute of Aerodynamics and Flow Technology, Braunschweig 2 CentaurSoft, Austin, TX 10. June / 22
2 HioCFD-4 Computational/Meshing Challenge DLR-F11 Meshing challenge: Quadratic hybrid/mixed-element meshes by the CENTAUR grid generator Computational challenge: Stabilized DG discretization and implicit solver in the DLR-PADGE code DLR-F11 high-lift configuration with slat tracks and flap track fairings hybrid/mixed-element 3 rd -order grid 3 rd -order flow solution 2 / 22
3 Linear and quadratic mesh generation process (CENTAUR) (1) Generate a coarse, linear hybrid grid. (2) Insert an additional midpoint for every grid edge and quadrilateral face. (3) Use the CAD information to map each new boundary point onto the underlying CAD surface. (4) Adjust the position of the interior points based on the mapped position of the boundary edge midpoints, in order to prevent self-intersecting grid elements and to ensure grid validity. (5) Quadratic interpolation of points and interior points for a curved element representation. (1) (2) (3) (4) (5) 3 / 22
4 Linear and quadratic mesh generation process (CENTAUR) For the grid generation of a particularly coarse linear mesh a lower analytic curvature clustering, and a larger maximum element size are used than for linear meshes typically generated for 2nd order schemes. After curving of the mesh the grid quality is evaluated using the volumes, and the the volume ratios of the complete elements and sub-elements. 4 / 22
5 Linear and quadratic mesh generation process: Examples linear mesh quadratic mesh linear mesh quadratic mesh 5 / 22
6 Linear and quadratic mesh generation: Application to DLR-F11 Coarse mesh with prisms pyramids tetrahedra elements Outboard flap track fairing surface mesh: prism layers around the slat and leading edge of the main linear mesh quadratic mesh 6 / 22
7 Linear and quadratic mesh generation: Application to DLR-F11 Coarse mesh with prisms pyramids tetrahedra elements Hybrid mesh cuts in the slat region: linear mesh quadratic mesh 7 / 22
8 The DG discretization of the compressible Euler equations The problem: u F c (u) = 0 in Ω R 2 + u h h, κ with u = (ϱ, ϱv 1, ϱv 2, ϱe) T. The DG(p) discretization: Find u h in V p h such that N h (u h, v h ) { κ T h κ } F c (u h ) : v h dx + ĥ(u + h, u h, n) v+ h ds κ\γ + ĥ Γ (u + h, n) v+ h ds = 0 v h V p h, Γ Numerical flux function ĥ: Roe flux, local Lax-Friedrichs flux,... Numerical flux function ĥ Γ at the boundary: ĥ Γ (u + h, n) = n F c (u Γ (u + h )) allows for an adjoint consistent discretization ĥ Γ (u + h, n) = ĥ(u+ h, u Γ (u+ h ), n) is, to our experience, significantly more stable 8 / 22
9 Numerical flux function at the boundary Use the same numerical flux on the boundary like on interior faces 1 2 ĥ Γ,h = ĥγ(u + h, n) = ĥ(u+ h, u Γ (u+ h ), n), where ( the boundary exterior state u Γ (u+ h ) is computed from u + h + u Γ (u+ h )) = u Γ (u + h ) with n v Γ = 0, the discretization is given by F c (u h ) : h v h dx+ ĥ h v + h ds+ ĥ(u + h, u Γ (u+ h ), n) v+ h ds = 0. Ω κ T κ\γ Γ h This discretization is adjoint inconsistent in combination with both J h (u h ) = p ( u Γ (u + h )) n ψ ds = J(u Γ (u + h )), Γ W and J(u h ) = p(u h ) n ψ ds. Γ W 9 / 22
10 Adjoint consistent discretization of force coefficients Consider the discretization: find u h in V p h such that F c (u h ) : v h dx + ĥ h v + h ds + ĥ Γ,h v + h ds = 0 Ω κ T κ\γ Γ h v h V p h (1) Consider the target quantity and its discretization J(u) = p(u) n ψ ds, J h (u h ) = ĥ Γ,h ψ ds, (2) Γ W Γ W with ψ = (0, ψ 1, ψ 2, 0) for ψ = (ψ 1, ψ 2 ), and ĥ Γ,h = ĥ Γ (u + h, n) = ĥ(u+ h, u Γ (u+ h ), n), It can be shown 1, that (1) in combination with (2) is adjoint consistent. - 1 R. Hartmann and T. Leicht. Generalized adjoint consistent treatment of wall boundary conditions for compressible flows, J. Comput. Phys., 300: , / 22
11 Example: Inviscid flow around the NACA0012 airfoil M = 0.5, α = 2, numerical flux at the wall boundary: ĥ Γ,h = ĥ(u + h, u Γ (u+ h ), n) z 1 isolines of the discrete adjoint solution z h for C dp : y 0.02 y x J h (u h ) = J(u Γ (u + h )) adjoint inconsistent x J h (u h ) = Γ W ĥ Γ,h ψ ds adjoint consistent 11 / 22
12 Extension to viscous flows Ω ( F c (u h ) + F v (u h, h u h )) : h v h dx (ĥγ,h ˆσ Γ,h n) v h ds is adjoint consistent in combination with J(u) = (p n τ n) ψ ds, J h (u h ) = Γ W with ψ = (0, ψ 1, ψ 2, 0) for ψ = (ψ 1, ψ 2 ). Γ W Γ (ĥγ,h ˆσ Γ,h n) ψ ds, 12 / 22
13 Adjoint consistent treatment of integral and local quantities Discretization of flow equations: ( F c (u h ) + F v (u h, h u h )) : h v h dx (ĥγ,h ˆσ Γ,h n) v h ds Ω Adjoint consistent discretization of integral quantities (drag, lift coefficients): J(u) = (p n τ n) ψ ds, J h (u h ) = (ĥγ,h ˆσ Γ,h n) ψ ds. Γ W Adjoint consistent discretization of local quantities (surface pressure, skin friction): c p (u) = p(u) p 1 2 ρ v 2 Γ W, c p,h (u h ) = ĥγ,h ñ p 1 2 ρ v 2 c f (u, u) = τ W (u, u) 1 2 ρ, c v 2 f,h (u h, u h ) = with ñ = (0, n 1, n 2, 0) for the normal vector n = (n 1, n 2 ), and t = (0, t 1, t 2, 0) for the tangential vector t = (t 1, t 2 ). Γ (ˆσΓ,h n ) t 1 2 ρ, v 2, 13 / 22
14 Example: L1T2 high-lift configuration RANS & Wilcox k-ω model, M = 0.197, α = 20.18, Re = ĥ Γ,h = ĥγ(u + h, n) = ĥ(u+ h, u Γ (u+ h ), n). c p- and c f -distributions of the p = 1 flow solution 14 / 22
15 Regularity check of a curved mesh Discrete regularity check: Consider the mapping σ from the reference element ˆκ to the element κ in physical space. The determinant of its Jacobian matrix dσ dx is required in numerical quadrature: κ = det(j) dx det(j(ˆx q ))w q ˆκ q Check whether the Jacobian determinants det(j(ˆx q )) are positive in each of the quadrature points. Note, that this discrete regularity check checks necessary conditions only, whereas the check based on Bézier functions 1 would give a necessary and sufficient condition for validity. - 1 A. Johnen, J.-F. Remacle, and C. Geuzaine. Geometrical validity of curvilinear finite elements. J. Comput. Phys., 233: , / 22
16 Regularity check of the quadratic DLR-F11 mesh κ = det(j) dx ˆκ q det(j(x q ))w q Discrete regularity check: det(j(x q )) > 0 for all quadrature points x q? Discrete regularity check of the quadratic mesh around the DLR-F11 mesh with curved elements (prisms, pyramids and tetrahedras): degree # DoFs/eqn # quadrature # irregular worst # neg. points/elem. elements Jac./element / / /64 16 / 22
17 Application to the DLR-F11 (Config 4): The solution process RANS & Wilcox k-ω model, flow conditions: M = 0.175, Re = , α = 7 Flow solver: Fully implicit (Backward-Euler) solver Convergence criterion: nonlinear residual below Linear systems solved with GMRes and ILU per process The CFL number is increased as the nonlinear residual decreases Diverging steps are recomputed with CFL/2. Numerical boundary flux: ĥγ = ĥ(u+ h, u Γ (u+ h ), n). degree # DoFs/eqn res tol start CFL end CFL steps lost 0 3.5e / e / / e / / 22
18 Application to the DLR-F11 (Config 4): The solution process RANS & Wilcox k-ω model, flow conditions: M = 0.175, Re = , α = 7 Flow solver: Fully implicit (Backward-Euler) solver Convergence criterion: nonlinear residual below Linear systems solved with GMRes and ILU per process The CFL number is increased as the nonlinear residual decreases legend numerical flux recompute diverging steps with CFL/2 p=1 (1) ĥ(u + h, u Γ (u+ h ), n) yes p=1 (2) n F c (u Γ (u + h )) yes p=1 (3) n F c (u Γ (u + h )) no 18 / 22
19 DLR-F11 high lift configuration: The flow solution RANS & Wilcox k-ω model, flow conditions: M = 0.175, Re = , α = 7. 3 rd -order solution on the quadratic (3 rd -order) mesh: 19 / 22
20 DLR-F11 high lift configuration: The flow solution RANS & Wilcox k-ω model, flow conditions: M = 0.175, Re = , α = 7. 3 rd -order solution on the quadratic (3 rd -order) mesh: separation above flap track consider pressure tab locations 01, 04, 06, 10 fairings and at wing tip 20 / 22
21 DLR-F11 high lift configuration: The flow solution RANS & Wilcox k-ω model M = Re = α = 7 PS01(η=0.150) PS04(η=0.449) 3 rd -order solution on 3 rd -order mesh PS06(η=0.681) PS10(η=0.891) C L C D C M Exp TAU % +4.1% -4.8% DG-O % +2.1% -5.8% 21 / 22
22 Computational/Meshing Challenge DLR-F11: Summary Quadratic hybrid/mixed-element mesh by the CENTAUR grid generator. Stabilized DG discretization and solver in the DLR-PADGE code: For stability use interior flux also on the boundary. For adjoint consistency use the same flux also for the evaluation of integral quantities (force coefficients), and for an according evaluation of local quantities (c p- and c f -distributions) Recompute diverging solver steps with CFL/2. DLR-F11 high-lift configuration with slat tracks and flap track fairings hybrid/mixed-element 3 rd -order grid 3 rd -order flow solution 22 / 22
DG 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 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 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 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 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 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 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 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 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 informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
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 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 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 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 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 informationTotal Pressure Losses Minimization in Turbomachinery Cascades, using a New Continuous Adjoint Formulation
Total Pressure Losses Minimization in Turbomachinery Cascades, using a New Continuous Adjoint Formulation D.I. Papadimitriou, K.C. Giannakoglou Laboratory of Thermal Turbomachines, National Technical University
More informationContinuous adjoint based error estimation and r-refinement for the active-flux method
Continuous adjoint based error estimation and r-refinement for the active-flux method Kaihua Ding, Krzysztof J. Fidkowski and Philip L. Roe Department of Aerospace Engineering, University of Michigan,
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 informationFLOW OVER A GIVEN PROFILE IN A CHANNEL WITH DYNAMICAL EFFECTS
Proceedings of Czech Japanese Seminar in Applied Mathematics 2004 August 4-7, 2004, Czech Technical University in Prague http://geraldine.fjfi.cvut.cz pp. 63 72 FLOW OVER A GIVEN PROFILE IN A CHANNEL WITH
More 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 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 informationDefinitions. Temperature: Property of the atmosphere (τ). Function of altitude. Pressure: Property of the atmosphere (p). Function of altitude.
Definitions Chapter 3 Standard atmosphere: A model of the atmosphere based on the aerostatic equation, the perfect gas law, an assumed temperature distribution, and standard sea level conditions. Temperature:
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 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 informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE Kevin G. Wang (1), Patrick Lea (2), and Charbel Farhat (3) (1) Department of Aerospace, California Institute of Technology
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 informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
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 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 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 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 informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
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 informationCost target result Temperature
Adjoint Techniques for Aerodynamic Shape Optimization Problems: New Mathematical Formulations and Software Assessment for Well Posed Problems with Acceptable Solutions I. Kampolis, Student D. Papadimitriou,
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 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 informationAIRFRAME NOISE MODELING APPROPRIATE FOR MULTIDISCIPLINARY DESIGN AND OPTIMIZATION
AIRFRAME NOISE MODELING APPROPRIATE FOR MULTIDISCIPLINARY DESIGN AND OPTIMIZATION AIAA-2004-0689 Serhat Hosder, Joseph A. Schetz, Bernard Grossman and William H. Mason Virginia Tech Work sponsored by NASA
More informationGerman Aerospace Center (DLR)
German Aerospace Center (DLR) AEROGUST M30 Progress Meeting 23-24 November 2017, Bordeaux Presented by P. Bekemeryer / J. Nitzsche With contributions of C. Kaiser 1, S. Görtz 2, R. Heinrich 2, J. Nitzsche
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 informationDesign of a Droopnose Configuration for a Coanda Active Flap Application. Marco Burnazzi and Rolf Radespiel
Design of a Droopnose Configuration for a Coanda Active Flap Application Marco Burnazzi and Rolf Radespiel Institute of Fluid Mechanics, Technische Universität Braunschweig, 38108 Braunschweig, Germany
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 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 informationSection 26.1: Reporting Conventions. Section 26.2: Fluxes Through Boundaries. Section 26.3: Forces on Boundaries
Chapter 26. Alphanumeric Reporting FLUENT provides tools for computing and reporting integral quantities at surfaces and boundaries. These tools enable you to find the mass flow rate and heat transfer
More informationFOR many computational fluid dynamics (CFD) applications, engineers require information about multiple solution
AIAA Aviation 13-17 June 2016, Washington, D.C. 46th AIAA Fluid Dynamics Conference AIAA 2016-3809 Improved Functional-Based Error Estimation and Adaptation without Adjoints William C. Tyson *, Katarzyna
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 informationPhysic-based Preconditioning and B-Splines finite elements method for Tokamak MHD
Physic-based Preconditioning and B-Splines finite elements method for Tokamak MHD E. Franck 1, M. Gaja 2, M. Mazza 2, A. Ratnani 2, S. Serra Capizzano 3, E. Sonnendrücker 2 ECCOMAS Congress 2016, 5-10
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 informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
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 informationAnderson acceleration for time-dependent problems
for time-dependent s PROBLEM SETTNG : 1D FLOW FGURE : One-dimensional flow in a tube. incompressible inviscid gravity neglected conservative form PROBLEM SETTNG : 1D FLOW gu t + gu2 x ( ) g gu t + + 1
More informationImproved Seventh-Order WENO Scheme
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Eposition 4-7 January 2, Orlando, Florida AIAA 2-45 Improved Seventh-Order WENO Scheme Yiqing Shen Gecheng Zha Dept.
More informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
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 informationActive Flutter Control using an Adjoint Method
44th AIAA Aerospace Sciences Meeting and Exhibit 9-12 January 26, Reno, Nevada AIAA 26-844 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 9 12 Jan, 26. Active Flutter Control using an
More informationEfficient Hessian Calculation using Automatic Differentiation
25th AIAA Applied Aerodynamics Conference, 25-28 June 2007, Miami Efficient Hessian Calculation using Automatic Differentiation Devendra P. Ghate and Michael B. Giles Oxford University Computing Laboratory,
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 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 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 informationAerodynamics. High-Lift Devices
High-Lift Devices Devices to increase the lift coefficient by geometry changes (camber and/or chord) and/or boundary-layer control (avoid flow separation - Flaps, trailing edge devices - Slats, leading
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 informationNumerical Solution of Partial Differential Equations governing compressible flows
Numerical Solution of Partial Differential Equations governing compressible flows Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
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 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 informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationQUELQUES APPLICATIONS D UN SCHEMA SCHEMA MIXTE-ELEMENT-VOLUME A LA LES, A L ACOUSTIQUE, AUX INTERFACES
1 QUELQUES APPLICATIONS D UN SCHEMA SCHEMA MIXTE-ELEMENT-VOLUME A LA LES, A L ACOUSTIQUE, AUX INTERFACES O. ALLAIN(*), A. DERVIEUX(**), I. ABALAKIN(***), S. CAMARRI(****), H. GUILLARD(**), B. KOOBUS(*****),
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 informationGiven the water behaves as shown above, which direction will the cylinder rotate?
water stream fixed but free to rotate Given the water behaves as shown above, which direction will the cylinder rotate? ) Clockwise 2) Counter-clockwise 3) Not enough information F y U 0 U F x V=0 V=0
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 informationHIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume LIII, Number 1, March 2008 HIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IOAN TELEAGA AND JENS
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 informationThe Discontinuous Galerkin Method on Cartesian Grids with Embedded Geometries: Spectrum Analysis and Implementation for Euler Equations
The Discontinuous Galerkin Method on Cartesian Grids with Embedded Geometries: Spectrum Analysis and Implementation for Euler Equations by Ruibin Qin A thesis presented to the University of Waterloo in
More informationNUMERICAL INVESTIGATION ON THE POTENTIAL OF STEAM COOLING FOR THE SKYLON SPACEPLANE IN HYPERSONIC FLOW
28 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES NUMERICAL INVESTIGATION ON THE POTENTIAL OF STEAM COOLING FOR THE SKYLON SPACEPLANE IN HYPERSONIC FLOW Th. Eggers* *German Aerospace Center (DLR)
More informationApplication of turbulence. models to high-lift airfoils. By Georgi Kalitzin
Center for Turbulence Research Annual Research Briefs 997 65 Application of turbulence models to high-lift airfoils By Georgi Kalitzin. Motivation and objectives Accurate prediction of the ow over an airfoil
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 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 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 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 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 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 information2010 Mathematics Subject Classification. 49Q10, 49Q12, 65K10, 76Nxx, 76M30. Introduction
SHAPE DERIVATIVES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS IN VARIATIONAL FORM arxiv:32.586v [math.oc] 20 Dec 203 Matthias Sonntag, Stephan Schmidt 2 and Nicolas R. Gauger Abstract. Shape optimization
More informationAn Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond
An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond Antony Jameson Department of Aeronautics and Astronautics AIAA Aerospace Sciences Meeting, Reno, NV AIAA Paper 2007-0037
More informationSIMULATION OF GAS FLOW OVER MICRO-SCALE AIRFOILS USING A HYBRID CONTINUUM-PARTICLE APPROACH
33rd AIAA Fluid Dynamics Conference and Exhibit 3-6 June 3, Orlando, Florida AIAA 3-44 33 rd AIAA Fluid Dynamics Conference and Exhibit / Orlando, Florida / 3-6 Jun 3 SIMULATION OF GAS FLOW OVER MICRO-SCALE
More informationSimulation of low Mach number flows
Simulation of low Mach number flows Improvement of accuracy and convergence of the TAU code using time derivative preconditioning Ralf Heinrich, Braunschweig, 9th February 008 Simulation of low Mach number
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 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 informationTime-resolved adaptive FEM simulation of the DLR-F11 aircraft model at high Reynolds number
Time-resolved adaptive FEM ulation of the DLR-F11 aircraft model at high Reynolds number Johan Hoffman Computational Technology Laboratory, School of Computer Science and Communication, KTH Johan Jansson
More information( ) A i,j. Appendices. A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in
Appendices A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in Eq.(3.11), can be found in the literature [9,172,173] and are therefore
More informationAerodynamic Inverse Design and Shape Optimization via Control Theory
Aerodynamic Inverse Design and Shape Optimization via Control Theory 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics & Astronautics Stanford University SciTech 2015 January 5, 2015
More informationDiscrete and continuous adjoint method for compressible CFD J. Peter ONERA
Discrete and continuous adjoint method for compressible CFD J. Peter ONERA J. Peter 1 1 ONERA DMFN October 2, 2014 J. Peter (ONERA DMFN) Adjoint method for compressible CFD October 2, 2014 1 / 60 Outline
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 informationLab Reports Due on Monday, 11/24/2014
AE 3610 Aerodynamics I Wind Tunnel Laboratory: Lab 4 - Pressure distribution on the surface of a rotating circular cylinder Lab Reports Due on Monday, 11/24/2014 Objective In this lab, students will be
More informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
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 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 informationReduced order model of three-dimensional Euler equations using proper orthogonal decomposition basis
Journal of Mechanical Science and Technology 24 (2) (2010) 601~608 www.springerlink.com/content/1738-494x DOI 10.1007/s12206-010-0106-0 Reduced order model of three-dimensional Euler equations using proper
More information1082. Computational study on aerodynamic characteristics of a flying wing MAV
1082. Computational study on aerodynamic characteristics of a flying wing MAV Hongming Cai 1, Zhilin Wu 2, Tianhang Xiao 3, Shuanghou Deng 4 1, 2 School of Mechanical Engineering, Nanjing University of
More information