HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
|
|
- Candice Miles
- 6 years ago
- Views:
Transcription
1 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 Politécnica de Madrid jesus.garicano.mena@upm.es III/04/2015
2 2 / 36 Outline 1 Introduction 2 Goal of this presentation 3 Shocked flows Perfect Gas flow over Adiabatic Flat Plate Inviscid Perfect Gas flow around Cylinder 4 Conclusions and Future Work
3 3 / 36 Outline Introduction 1 Introduction 2 Goal of this presentation 3 Shocked flows Perfect Gas flow over Adiabatic Flat Plate Inviscid Perfect Gas flow around Cylinder 4 Conclusions and Future Work
4 4 / 36 Objective Introduction Can Discontinuous Galerkin Spectral Element Methods predict accurately hypersonic stagnation heating?
5 5 / 36 Objective Introduction Can Discontinuous Galerkin Spectral Element Methods predict accurately hypersonic stagnation heating?
6 Introduction Stagnation heating 6 / 36 The stagnation heating problem Requirements Robust Conservative Monotonicity respecting Accurate (boundary layer) V > 5a
7 Introduction Stagnation heating 7 / 36 Governing Equations Conservation form Conserved variables U t + F c = F d U = [ ρ, ρu j, ρe ] t Continuous regime Kn 1 Perfect Gas model F c = Convective flux tensor ρ u j ρ u i u j + pδ i,j ρ H u j = F c x j 1 t x j F d = Diffusive flux tensor 0 τ i,j τ i,j u j q i = F d x j 1 t x j
8 8 / 36 Introduction Discontinuous Galerkin Spectral Element Methods Discontinuous Galerkin Spectral Element Method [1] Galerkin: weak formulation ( U ) t, ϕ + ( F c F d ), ϕ = 0, ϕ V }{{} F L 2 projection (u, v) = Ω u v dω [1] Implementing Spectral Methods for Partial Differential Equations, D.A. Kopriva, 2009.
9 Introduction Discontinuous Galerkin Spectral Element Methods 9 / 36 Discontinuous Galerkin Spectral Element Method Domain Ω and a possible tesselation Ω h. n Elem Ω Ω h = Ω k Ω i Ω j = if i j. k=1 Physical and computational domains: iso-parametric mapping
10 Introduction Discontinuous Galerkin Spectral Element Methods Discontinuous Galerkin Spectral Element Method Lagrange polynomial expansion - 1D n U (x, t) U DoF N (x, t) = U J (t) l J (x), l α(x β ) = δ α,β J=1 Spectral magic Choose as nodes the zeros of orthogonal polynomials 1.0 Approximation is optimal Exact quadrature for polynomials up to degree 2N T 3 T 4 T L 2 L 3 L 4 L 5 T k (x) 0 L k (x) T x Chebyshev polynomials x Legendre Polynomials 10 / 36
11 11 / 36 Introduction Discontinuous Galerkin Spectral Element Methods Discontinuous Galerkin Spectral Element Method Discrete form ( U ) ( t, ϕn + F, ) ϕ N = 0, ϕ N V h 0 Ωk Approximate solution and fluxes U U N N (ξ, η) = U µ,νφ N µ,ν ; µ,ν=0 N F F N (ξ, η) = µ,ν=0 F N µ,νφ µ,ν Basis functions 2D - tensor product φ µ,ν = l µ (ξ)l ν (η) P N P N
12 12 / 36 Introduction Discontinuous Galerkin Spectral Element Methods Discontinuous Galerkin Spectral Element Method ( ϕ N = φ N = Test function N µ,ν=0 ϕ N µ,νφ µ,ν Linearity of projection operator U ) N ), φ i,j + ( F N, φ i,j = 0, i, j = 0,..., N t
13 13 / 36 Introduction Discontinuous Galerkin Spectral Element Methods Discontinuous Galerkin Spectral Element Method ( Departing point U ) N ), φ i,j + ( F N, φ i,j = 0, i, j = 0,..., N t ( Integration by parts U ) N ), φ i,j + F N φ i,j 1 ext dσ ( F N, φ i,j = 0, i, j = 0,..., N t δω k Integration Numerical Quadrature 1 1 v(ξ, η)dξdη = 0 0 n,m=0 N v(ξ n, η m)w nw m v P 2N+1,2N+1
14 Introduction Discontinuous Galerkin Spectral Element Methods Discontinuous Galerkin Spectral Element Method ( Time derivative term U ) N d, φ i,j = w i w j U t dt i,j N Space operator - edge (surface) contribution F N φ i,j 1 ext dσ = F N (1, η j )φ i,j (1, η j )w j F N (0, η j )φ i,j (0, η j )w j δω k + F N (ξ i, 1)φ i,j (ξ i, 1)w i F N (ξ i, 0)φ i,j (ξ i, 0)w i Space operator - area (volume) contribution ) ( F N, φ i,j = N i,j=0 F N i,j φ i,j w i w j 14 / 36
15 15 / 36 Introduction Discontinuous Galerkin Spectral Element Methods Discontinuous Galerkin Spectral Element Method d w i w j U dt i,j N + δω k Assembled discrete ode F N, φ i,j 1 ext dσ N i,j=0 w i w j F N i,j φ i,j = 0, i, j = 1,..., N; k = 1,..., n Elem R N ( U) = Discrete partial differential operator δω k F 1 ext φ i,j d σ }{{} Information propagation N i,j=0 w i w j F N i,j φ i,j i, j = 0,..., N
16 Introduction Discontinuous Galerkin Spectral Element Methods Discontinuous Galerkin Spectral Element Method d w i w j U dt i,j N + δω k Assembled discrete ode F N, φ i,j 1 ext dσ N i,j=0 w i w j F N i,j φ i,j = 0, i, j = 1,..., N; k = 1,..., n Elem Boundary Conditions [1] Weakly enforced, exploiting characteristic information [1] Implementing Spectral Methods for Partial Differential Equations, D.A. Kopriva, / 36
17 17 / 36 Outline Goal of this presentation 1 Introduction 2 Goal of this presentation 3 Shocked flows Perfect Gas flow over Adiabatic Flat Plate Inviscid Perfect Gas flow around Cylinder 4 Conclusions and Future Work
18 18 / 36 Goal of this presentation Goal of this presentation Evaluate DGSEM performance for Perfect Gas (PG) shocked flows Two configurations: Flat plate Straight cylinder
19 19 / 36 Goal of this presentation Handling shock waves Shock capturing term (Euler) [1] U t + F c = F d + F ad F ad = ε µ F d [1] Sub-cell Shock Capturing for DG Methods, P.-O. Persson and J. Peraire, AIAA
20 Goal of this presentation Handling shock waves Shock capturing term (Euler) U t + F c = F d + F ad F ad = ε µ F d Artificial viscosity 0 if ( s e s 0 κ ) ε ε = sin( π(se s 0) ) 2κ ε 0 if s 0 + κ s e ε 0 h N, s 0 if s 0 κ s e s 0 + κ 1, κ empirically chosen N4 s e = lg S e, with shock sensor S e = ( p N p N 1, p N p N 1) (p N, p N ) 20 / 36
21 21 / 36 Outline Shocked flows 1 Introduction 2 Goal of this presentation 3 Shocked flows Perfect Gas flow over Adiabatic Flat Plate Inviscid Perfect Gas flow around Cylinder 4 Conclusions and Future Work
22 22 / 36 Flat plate in PG flow Shocked flows Perfect Gas flow over Adiabatic Flat Plate Ma Re Ld Flat plate. Definition of the domain. Surface y = 0, x (0, L d ): no-slip adiabatic wall Surface x = L d is a supersonic outlet Surfaces x = L u and y = H: supersonic inlets, with the upstream state U Surface y = 0, x ( L u, 0): slip wall
23 23 / 36 Shocked flows Perfect Gas flow over Adiabatic Flat Plate Viscous Flow over Flat Plate - Ma = 2.0, Re = 10 4 N = 2 N = 3 Non-dimensional pressure field. O2 vs O3.
24 24 / 36 Shocked flows Perfect Gas flow over Adiabatic Flat Plate Viscous Flow over Flat Plate - Ma = 2.0, Re = 10 4 N = 2 N = 3 Non-dimensional temperature field. O2 vs O3.
25 25 / 36 Shocked flows Perfect Gas flow over Adiabatic Flat Plate Viscous Flow over Flat Plate - Ma = 2.0, Re = 10 4 N = 2 N = 3 Ma field. O2 vs O3.
26 Shocked flows Perfect Gas flow over Adiabatic Flat Plate Viscous Flow over Flat Plate - Ma = 2.0, Re = 10 4 Wall quantities - C p. Reference value from [1]. [1] A second-order weak interaction expansion for moderately hypersonic flow past a flat plate, Ko and Kubota, AIAAJ, Vol. 5, No. 10 (1967). 26 / 36
27 Shocked flows Perfect Gas flow over Adiabatic Flat Plate Viscous Flow over Flat Plate - Ma = 2.0, Re = 10 4 Wall quantities - C f. Reference value from [1]. [1] A second-order weak interaction expansion for moderately hypersonic flow past a flat plate, Ko and Kubota, AIAAJ, Vol. 5, No. 10 (1967). 27 / 36
28 28 / 36 Shocked flows Inviscid Perfect Gas flow around Cylinder Inviscid Flow around Cylinder - Ma = 1.2 N = 3 Non-dimensional pressure field.
29 29 / 36 Shocked flows Inviscid Perfect Gas flow around Cylinder Inviscid Flow around Cylinder - Ma = 1.2 N = 3 Non-dimensional temperature field.
30 30 / 36 Shocked flows Inviscid Perfect Gas flow around Cylinder Inviscid Flow around Cylinder - Ma = 1.2 N = 3 Ma field.
31 31 / 36 Outline Conclusions and Future Work 1 Introduction 2 Goal of this presentation 3 Shocked flows Perfect Gas flow over Adiabatic Flat Plate Inviscid Perfect Gas flow around Cylinder 4 Conclusions and Future Work
32 32 / 36 Conclusions and Future Work Conclusions and Future Work Can Discontinuous Galerkin Spectral Element Methods predict accurately hypersonic stagnation heating? Not yet...
33 Conclusions and Future Work 33 / 36 Conclusions and Future Work Can Discontinuous Galerkin Spectral Element Methods predict accurately hypersonic stagnation heating? Not yet... but now shortcomings have been identified cures are known (from RDS, FV, SUPG)
34 34 / 36 Conclusions and Future Work Conclusions and Future Work Can Discontinuous Galerkin Spectral Element Methods predict accurately hypersonic stagnation heating? To do list Improve robustness Turn implicit [1,2] Shock handling Explictly enthalpy-preserving [3] PDE-based artificial viscosity [4] [1] Jacobian-free Newton-Krylov methods, Knoll and Keyes, JCP, Vol.193, [2] Ph.D. dissertation, Universidad Politécnica de Madrid (in preparation), C. Redondo Calle. [3] Ph.D. dissertation, University of Texas at Austin, Kirk, [4] Shock capturing with PDE-based artificial viscosity for DGFEM, Barter and Darmofal, JCP, Vol.225, 2010.
35 35 / 36 Acknowledgements Conclusions and Future Work This research has been funded by Universidad Politécnica de Madrid
36 Conclusions and Future Work 36 / 36 Thanks for your attention Questions?
Sub-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 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 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 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 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 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 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 information58:160 Intermediate Fluid Mechanics Bluff Body Professor Fred Stern Fall 2014
Professor Fred Stern Fall 04 Chapter 7 Bluff Body Fluid flows are broadly categorized:. Internal flows such as ducts/pipes, turbomachinery, open channel/river, which are bounded by walls or fluid interfaces:
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 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 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 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 informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu and Roy H. Stogner roystgnr@cfdlab.ae.utexas.edu Univ. of Texas at Austin September 9, 2008 1 Introduction 2 Weighted Residuals
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
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 informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationGeneralised Summation-by-Parts Operators and Variable Coefficients
Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation
More informationInterior penalty tensor-product preconditioners for high-order discontinuous Galerkin discretizations
Interior penalty tensor-product preconditioners for high-order discontinuous Galerkin discretizations Will Pazner Brown University, 8 George St., Providence, RI, 9, U.S.A. Per-Olof Persson University of
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 informationChapter 7: Natural Convection
7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection
More informationEffect of Applied Magnetic Field on Shock Boundary Layer Interaction
University of Kentucky UKnowledge Mechanical Engineering Faculty Publications Mechanical Engineering 1-2012 Effect of Applied Magnetic Field on Shock Boundary Layer Interaction Ovais U. Khan University
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 informationEfficient and Flexible Solution Strategies for Large- Scale, Strongly Coupled Multi-Physics Analysis and Optimization Problems
University of Colorado, Boulder CU Scholar Aerospace Engineering Sciences Graduate Theses & Dissertations Aerospace Engineering Sciences Spring 1-1-2016 Efficient and Flexible Solution Strategies for Large-
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 informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
More informationControl Theory Approach to Aero Shape Optimization
Control Theory Approach to Aero Shape Optimization. Overview. Finite Volume Method 3. Formulation of Optimal Design Problem 4. Continuous Adjoint Approach 5. Discrete Adjoint Approach Overview of Adjoint
More informationEntry Aerodynamics MARYLAND U N I V E R S I T Y O F. Entry Aerodynamics. ENAE Launch and Entry Vehicle Design
Atmospheric Regimes on Entry Basic fluid parameters Definition of Mean Free Path Rarified gas Newtonian flow Continuum Newtonian flow (hypersonics) 2014 David L. Akin - All rights reserved http://spacecraft.ssl.umd.edu
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 informationExcerpt from the Proceedings of the COMSOL Users Conference 2006 Boston
Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu
More informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationWALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID
WALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID Zia Ghiasi sghias@uic.edu Dongru Li dli@uic.edu Jonathan Komperda jonk@uic.edu Farzad
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
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 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 informationNumerical Simulation of Hypersonic Boundary-Layer Instability in a Real Gas with Two-Dimensional Surface Roughness
AIAA Aviation 22-26 June 2015, Dallas, TX 45th AIAA Fluid Dynamics Conference AIAA 2015-3077 Numerical Simulation of Hypersonic Boundary-Layer Instability in a Real Gas with Two-Dimensional Surface Roughness
More informationHydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010
Hydrodynamics Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 What is Hydrodynamics? Describes the evolution of physical systems (classical or quantum particles, fluids or fields) close to thermal
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 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 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 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 informationDevelpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow
Develpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow F. Piscaglia, A. Montorfano Dipartimento di Energia, POLITECNICO DI MILANO Content Introduction
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationarxiv: v1 [physics.comp-ph] 21 Nov 2018
Implicit High-Order Flux Reconstruction Solver for High-Speed Compressible Flows Ray Vandenhoeck a,, Andrea Lani b a Department of Mechanical Engineering, KU Leuven, Celestijnenlaan 300, 3001, Leuven,
More informationA Hybridized Discontinuous Galerkin Method for Time-Dependent Compressible Flows
A Hybridized Discontinuous Galerkin Method for Time-Dependent Compressible Flows Master s Thesis by Alexander Jaust Institut für Geometrie und Praktische Mathematik Rheinisch-Westfälische Technische Hochschule
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 informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationMultiscale Hydrodynamic Phenomena
M2, Fluid mechanics 2014/2015 Friday, December 5th, 2014 Multiscale Hydrodynamic Phenomena Part I. : 90 minutes, NO documents 1. Quick Questions In few words : 1.1 What is dominant balance? 1.2 What is
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 information3 Generation and diffusion of vorticity
Version date: March 22, 21 1 3 Generation and diffusion of vorticity 3.1 The vorticity equation We start from Navier Stokes: u t + u u = 1 ρ p + ν 2 u 1) where we have not included a term describing a
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 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 informationMHD Linear Stability Analysis Using a Full Wave Code
US-Japan JIFT Workshop on Progress of Extended MHD Models NIFS, Toki,Japan 2007/03/27 MHD Linear Stability Analysis Using a Full Wave Code T. Akutsu and A. Fukuyama Department of Nuclear Engineering, Kyoto
More informationA Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations
A Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations Institute for Mathematics Applied to the Geosciences (IMAGe) National Center for Atmospheric Research (NCAR) Boulder CO
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 informationNonlinear Evolution of a Vortex Ring
Nonlinear Evolution of a Vortex Ring Yuji Hattori Kyushu Institute of Technology, JAPAN Yasuhide Fukumoto Kyushu University, JAPAN EUROMECH Colloquium 491 Vortex dynamics from quantum to geophysical scales
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
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 informationRobust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations
Robust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with
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. 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 informationModeling Shock Waves Using Exponential Interpolation Functions with the Least-Squares Finite Element Method
Old Dominion University ODU Digital Commons Mechanical & Aerospace Engineering Theses & Dissertations Mechanical & Aerospace Engineering Spring 26 Modeling Shock Waves Using Exponential Interpolation Functions
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
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 informationFINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS
FINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS Lucia Catabriga Department of Computer Science, Federal University of Espírito Santo (UFES) Av. Fernando Ferrari,
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 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 informationSimulation of Condensing Compressible Flows
Simulation of Condensing Compressible Flows Maximilian Wendenburg Outline Physical Aspects Transonic Flows and Experiments Condensation Fundamentals Practical Effects Modeling and Simulation Equations,
More informationIntroduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)
Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationA numerical study of heat transfer and fluid flow over an in-line tube bank
Fluid Structure Interaction VI 295 A numerical study of heat transfer and fluid flow over an in-line tube bank Z. S. Abdel-Rehim Mechanical Engineering Department, National Research Center, Egypt Abstract
More informationPerformance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain
March 4-5, 2015 Performance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain M. Bonnasse-Gahot 1,2, H. Calandra 3, J. Diaz 1 and S. Lanteri 2
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 informationComputation of Incompressible Flows: SIMPLE and related Algorithms
Computation of Incompressible Flows: SIMPLE and related Algorithms Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de SIMPLE-Algorithm I - - - Consider
More informationNumerical Study of Natural Unsteadiness Using Wall-Distance-Free Turbulence Models
Numerical Study of Natural Unsteadiness Using Wall-Distance-Free urbulence Models Yi-Lung Yang* and Gwo-Lung Wang Department of Mechanical Engineering, Chung Hua University No. 707, Sec 2, Wufu Road, Hsin
More informationIMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS
Approved for public release; distribution is unlimited IMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS Sean L. Gibbons Captain, United
More informationStructural LES modeling with high-order spectral difference schemes
Center for Turbulence Research Annual Research Briefs 211 123 Structural LES modeling with high-order spectral difference schemes By G. Lodato, P. Castonguay AND A. Jameson 1. Motivation and objectives
More informationADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS
ECCOMAS Computational Fluid Dynamics Conference 2001 Swansea, Wales, UK, 4-7 September 2001 c ECCOMAS ADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS Tayfun E. Tezduyar Team for Advanced
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationarxiv: v2 [math.na] 9 Apr 2018
oname manuscript o. will be inserted by the editor) The BR1 Scheme is Stable for the Compressible avier-stokes Equations Gregor J. Gassner Andrew R. Winters Florian J. Hindenlang David A. Kopriva arxiv:1704.03646v
More informationFundamentals of compressible and viscous flow analysis - Part II
Fundamentals of compressible and viscous flow analysis - Part II Lectures 3, 4, 5 Instantaneous and averaged temperature contours in a shock-boundary layer interaction. Taken from (Pasquariello et al.,
More informationSpectral element schemes for the. Korteweg-de Vries and Saint-Venant equations
Spectral element schemes for the Korteweg-de Vries and Saint-Venant equations R. PASQUETTI a a. Université Côte d Azur, CNRS, Inria, LJAD, France, richard.pasquetti@unice.fr... Résumé : Les sytèmes hyperboliques
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
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 informationDevelopment of Two-Dimensional Convergent-Divergent Nozzle Performance Rapid Analysis Project
International Forum on Energy, Environment Science and Materials (IFEESM 015) Development of Two-Dimensional Convergent-Divergent Nozzle Performance Rapid Analysis Project Yaxiong YANG 1,a *, Eri Qitai,b,
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 informationA STABLE PENALTY METHOD FOR THE COMPRESSIBLE NAVIER STOKES EQUATIONS: III. MULTIDIMENSIONAL DOMAIN DECOMPOSITION SCHEMES
SIAM J. SCI. COMPUT. Vol. 20, No., pp. 62 93 c 998 Society for Industrial and Applied Mathematics A STABLE PENALTY METHOD FOR THE COMPRESSIBLE NAVIER STOKES EQUATIONS: III. MULTIDIMENSIONAL DOMAIN DECOMPOSITION
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 information. α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω. Friday April 1 ± ǁ
. α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω Friday April 1 ± ǁ 1 Chapter 5. Photons: Covariant Theory 5.1. The classical fields 5.2. Covariant
More informationPREDICTIVE SIMULATION OF UNDERWATER IMPLOSION: Coupling Multi-Material Compressible Fluids with Cracking Structures
PREDICTIVE SIMULATION OF UNDERWATER IMPLOSION: Coupling Multi-Material Compressible Fluids with Cracking Structures Kevin G. Wang Virginia Tech Patrick Lea, Alex Main, Charbel Farhat Stanford University
More informationarxiv: v1 [math.na] 23 Jun 2017
On the use of kinetic energy preserving DG-schemes for large eddy simulation David Flad 1 arxiv:1706.07601v1 [math.na] 23 Jun 2017 Institute for Aerodynamics and Gas Dynamics, University of Stuttgart,
More informationFinite Element Solver for Flux-Source Equations
Finite Element Solver for Flux-Source Equations Weston B. Lowrie A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics Astronautics University
More informationON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY ON INTERIOR NODAL SETS
6th European Conference on Computational Mechanics (ECCM 6 7th European Conference on Computational Fluid Dynamics (ECFD 7 11 15 June 018, Glasgow, UK ON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY
More informationExam in Fluid Mechanics SG2214
Exam in Fluid Mecanics G2214 Final exam for te course G2214 23/10 2008 Examiner: Anders Dalkild Te point value of eac question is given in parentesis and you need more tan 20 points to pass te course including
More informationMultiscale hp-adaptive quiescent reactors
Multiscale hp-adaptive quiescent reactors C.E.Michoski Institute for Computational Engineering and Sciences (ICES) University of Texas (UT) at Austin Group Meeting Nov. 7, 2011 Outline 1 2 Equation (or
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationModule-5: Hypersonic Boundary Layer theory. Lecture-20: Hypersonic boundary equation
Module-5: Hypersonic Boundary Layer theory Lecture-0: Hypersonic boundary equation 0.1 Governing Equations for Viscous Flows The Navier-Stokes (NS) equaadtions are the governing equations for the viscous
More information