1st Committee Meeting
|
|
- Annabelle Hamilton
- 6 years ago
- Views:
Transcription
1 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
2 Outline 1 Introduction and Motivation 2 Existing Work 3 Expected Contributions 4 Logistics E. Liu (MIT) 1st Committee Meeting 2 / 41
3 Outline 1 Introduction and Motivation 2 Existing Work 3 Expected Contributions 4 Logistics E. Liu (MIT) 1st Committee Meeting 3 / 41
4 Motivation: Next Generation CFD Constantly-increasing computational power allows for higher-fidelity simulations of multiscale problems; e.g., unsteady RANS LES/DES/DNS aeroacoustics Most of these problems are unsteady Need for automatic error control and methods that can utilize DOFs efficiently (a) Axial velocity at combustor/turbine interface (Medic 2006) (b) Pressure field of SSLV (Mavriplis 2007) E. Liu (MIT) 1st Committee Meeting 4 / 41
5 Motivation: Space-Time Anisotropy Many unsteady problems are anisotropic in space-time e.g., wave-propagation: wave support propagation distance Appropriate mesh anisotropy efficiency gains (Yano et al. 2011) Recognized in the space-time literature: Key to an effective space-time adaptation strategy is attributing the error indicator to spatial versus temporal resolution using a measure of space-time anisotropy. (Luo, Fidkowski, 2011) Simplex meshes can match arbitrary solution anisotropy Consider a simple 1D + 1 example... E. Liu (MIT) 1st Committee Meeting 5 / 41
6 Motivation: DOF Reduction s=1 Ideal 1D "Realistic" 1D (a) 1D Meshes t x (b) Ideal Simplex Consider a 1D, stationary shock problem (Left) Ideal 1D mesh with 2 elements is shown Shock-tracking unreasonable; more realistic, 5 element mesh also shown Now consider a 1D, moving shock problem (Right) Ideal, 3 element simplex mesh is shown Shock path highlighted in red (s = 1) E. Liu (MIT) 1st Committee Meeting 6 / 41
7 Motivation: DOF Reduction t t x (a) Uniform Quad (N uni elem) x (b) Adapted Quad (O( N uni log N uni) elem) Assumed current state-of-the-art in ST adaptation: e.g., Bangerth 1999, Hartmann 2001, Fidkowski 2011 Sample mesh shown in (b) Hierarchical, (usually) tensor-product mesh structure 1 hanging node per edge E. Liu (MIT) 1st Committee Meeting 7 / 41
8 Motivation: DOF Reduction t (a) Anisotropic Simplex x Simplex-adaptive ST potential: N 2D = C 1 N 1D elements! Best-case scenario N 1D = C 2 + O(log N uni ), C 2 = number of elements spanning shock C 1 is a constant dependent on: maximum element anisotropy, discretization dissipation, dispersion properties Result: dimension of the problem is reduced E. Liu (MIT) 1st Committee Meeting 8 / 41
9 Sample Works: Adaptive Space-Time Methods Bar-Yoseph (1989) First application of DG to space-time problems 1D + 1 and 2D + 1 scalar advection Hartmann (2001, 2002) DWR-driven (isotropic) mesh adaptation 1D + 1 Burgers and Euler (a) Hartmann (2002) Luo, Fidkowski (2011, 2012) Applies anisotropy measures of Leicht, Hartmann (2010) Linear solve cost nearly reduced to spatial linear solve (Richter 2009) Navier-Stokes Note: All prior works (all examples here) use hierarchical meshes. E. Liu (MIT) 1st Committee Meeting 9 / 41
10 Objective Development of an autonomous, robust, and practical solution strategy for (certain) unsteady problems using an adaptive, space-time DG method. Approaches High-order DG space-time solver DWR-based error estimation Mesh-optimization adaptation scheme (Yano 2012) Metric-based mesh adaptation (BAMG, EPIC/MADCAP) E. Liu (MIT) 1st Committee Meeting 10 / 41
11 General Approach Higher-order solutions with the Discontinuous Galerkin (DG) method p-th order polynomial basis on each element; discontinuous between elements Nearest-neighbor stencil y x T H u(x,y) Anisotropic simplex meshes for complex geometries Output-based mesh adaptation: Error estimation on integrated outputs of interest (e.g., lift, drag) Element-wise localization of error estimates E. Liu (MIT) 1st Committee Meeting 11 / 41
12 Output Adaptive Solution Strategy Problem: Find output J(u) where u satisfies governing equations DGFEM: Find discrete output, J h,p, with u h,p V h,p s.t. R h,p (u h,p, v h,p ) = 0 v h,p V h,p. Output-based adaptation: estimate and control the error: ε J(u) J h,p (u h,p ). E. Liu (MIT) 1st Committee Meeting 12 / 41
13 Error Estimation Dual-Weighted Residual (DWR) of Becker and Rannacher (2001): ε = R h,p (u h,p, ψ) = R h,p[uu h,p ](u u h,p, ψ ψ h,p ) R h,p (u h,p, ψ h,p ) Localized error indicator: η K R h,p (u h,p, ψ h,p K ), where p > p. Primal solution features (left) and dual solution features (right). E. Liu (MIT) 1st Committee Meeting 13 / 41
14 Mesh Optimization via Error Sampling and Synthesis (MOESS) Driven by DWR (a posteriori) error estimator Resulting mesh resolves primal and dual solution features No a priori assumptions about solution irregularity Robust method can start from almost any mesh (a) Mesh-Metric Duality E. Liu (MIT) 1st Committee Meeting 14 / 41
15 Mesh Optimization via Error Sampling and Synthesis (MOESS) Mesh optimization (relaxed): M = arg inf M E(M) s.t. C(M) = DOF target 1 Sample error response to local metric modifications 2 Synthesize {M Ki, η Ki } pairs into local approximations of (roughly) E M 3 Solve optimization problem via steepest descent E. Liu (MIT) 1st Committee Meeting 15 / 41
16 Outline 1 Introduction and Motivation 2 Existing Work 3 Expected Contributions 4 Logistics E. Liu (MIT) 1st Committee Meeting 16 / 41
17 Overview Yano (2012) produced a series of preliminary adaptive space-time results. Results Included... 1D + 1/2D + 1 Wave Equation (not shown) 1D + 1/2D + 1 Euler Equation Yano s results are now presented to further motivate space-time schemes. Validation We will show the entropy mode (1D) and shear mode (2D) We will also show a 1D + 1 Riemann Problem (modified Sod s) E. Liu (MIT) 1st Committee Meeting 17 / 41
18 Euler Equation Validation (1D + 1) (Yano 2012) (a) entropy (b) mesh (c) convergence Gaussian perturbation of density at constant pressure Output: entropy perturbation at final time Adaptation reduces DOF usage substantially Mesh anisotropy is nearly perfect, with 3 elements spanning time Hence the effective dimension is reduced E. Liu (MIT) 1st Committee Meeting 18 / 41
19 Euler Equation Validation (2D + 1) (Yano 2012) (a) density Excitation of shear mode using bivariate Gaussian with the procedure of Wang (2006) High level of anisotropy observed in the mesh (b) mesh E. Liu (MIT) 1st Committee Meeting 19 / 41
20 Euler Equation Validation (2D + 1) (Yano 2012) Output: momentum perturbation at the final time Expected convergence rate: h 2p ; correlates to N p DOF Adapted cases superconverge Obtains 2D rates even though problem is 3D Appropriate anisotropy dimensionality reduction in 2D, N2p/3 DOF in 3D E. Liu (MIT) 1st Committee Meeting 20 / 41
21 (Modified) Sod s Problem: Inviscid Shock (1D + 1) (Yano 2012) Initial Conditions: x < 0: (ρ L, v L, p L ) = (1.0, 0.0, 1.0) x > 0: (ρ R, v R, p R ) = (0.4, 0.0, 0.4) (a) density (b) mesh Output: density perturbation at final time Space-time adaptation appears effective for this problem Shock-capturing (Barter 2008) not robust for stronger shocks E. Liu (MIT) 1st Committee Meeting 21 / 41
22 Outline 1 Introduction and Motivation 2 Existing Work 3 Expected Contributions 4 Logistics E. Liu (MIT) 1st Committee Meeting 22 / 41
23 Overview Objective Development of an autonomous, robust, and practical solution strategy for (certain) unsteady problems using an adaptive, space-time DG method. Challenges Performance: space-time is expensive without special handling of time Robustness: current shock-capturing scheme lacks robustness Expected Contributions Preconditioning, solver efficiency improvements Shock-capturing robustness and error estimation improvements Demonstration of capability on complex 2D + 1 problems (Time permitting) Exploration of adaptive space-time schemes applied to near-turbulent flow E. Liu (MIT) 1st Committee Meeting 23 / 41
24 Preconditioning: Long Journey Ahead Currently, we run steady 3D cases at... p = 2 or p = 3 with 1M DOFs (e.g., 100K p = 2 elements) Insufficient for RANS problems Potentially workable for Euler shock problems Memory Usage for p = 2 primal (p = 3 dual) (Jacobian Only) N elem p = 2, 3D p = 3, 3D p = 2, 4D p = 3, 4D 100K 9.31GB 37.3GB 402GB 2.19TB 1M 93.1GB 373GB 4.02TB 21.9TB 10M 931GB 3.73TB 40.2TB 219TB 100M 9.31TB 37.3TB 402TB 2.19PB For reference, the upcoming ACDL cluster will have 40 32GB = 1.28TB or about 50,000 4D elements. E. Liu (MIT) 1st Committee Meeting 24 / 41
25 Preconditioning Ignoring causality of time = high cost Time-direction has strong hyperbolic structure Potentially reduce cost with appropriate preconditioning Care needed so that causality is not violated May also reduce cost with solver modifications/optimizations Some Ideas Parallel decomposition cut purely in the time direction? Time-slab structure upper block triangular Jacobian Insert cut planes and use cut cell techniques? Use mesher capable of internal, material boundaries? Line-preconditioning: order elements along strongly causal lines and perform Gauss-Seidel-like sweep E. Liu (MIT) 1st Committee Meeting 25 / 41
26 Reducing Jacobian Storage t t (a) Time slabs via cut cells x (b) Meshing with time slabs x Time-slabs can help with the storage issue: Time-slabs allow serialization in time Store Jacobian for only one slab at a time Traditional domain decomposition (e.g., METIS) in space to parallelize Try avoiding Jacobian storage for p + 1 truth surrogate: Strong solvers (e.g., GMRES) implemented matrix-free Weak smoothers (e.g., block Jacobi) that do not require Jacobian storage E. Liu (MIT) 1st Committee Meeting 26 / 41
27 Shock-Capturing Currently use a scheme that governs artificial viscosity with a PDE Pros PDE adds artificial viscosity smoothly across elements Reduces oscillations/solution pollution downstream compared to element-wise constant viscosity Cons Lack of robustness for strong shocks (steady and space-time settings) Increased cost due to extra PDE coupled to flow equations How to handle error estimation? Don t necessarily want to adapt to features of the shock PDE Cons may be from poor tuning or incorrect extension to space-time Try p = 0 subdivision technique of Huerta, Casoni, Peraire (2011)? Error estimation questions remain with or without the shock PDE E. Liu (MIT) 1st Committee Meeting 27 / 41
28 Complex 2D + 1 Problems 3D + 1 problems are probably too cost-prohibitive (memory). Instead, consider more complex 2D + 1 problems, for example... Flow through rotor-stator rows: many complex interactions Stationary/rotating blade rows Interactions between wakes and downstream blades Shock-boundary layer interaction Complex heat-transfer problem Large, unsteady force variations Other problems possible; may try 3D + 1 depending on progress (c) Entropy, transonic turbine (Paniagua 2008) (d) Density, 3 stator-rotor rows (Cizmas 1998) E. Liu (MIT) 1st Committee Meeting 28 / 41
29 Complex 2D + 1 Problems: Lead-up Cannot start solving rotor-stator problems tomorrow. Work toward goal with simpler problems: 1D + 1 Sod s Problem, other traditional shock-tube problems Shock-tube with moving walls 2D + 1 Euler, Navier-Stokes, then RANS-SA with subsonic, then transonic flow Shock-tube, extruded Isolated airfoil with periodic BCs, time-varying BCs Euler-only: isolated airfoil with specified vortex-wake inflow Single blade rotor Multiple airfoils (stationary, then rotating) (Hopefully) Comparison with industry code from Pratt E. Liu (MIT) 1st Committee Meeting 29 / 41
30 Broader Classes of Problems Have noted that space-time makes sense for problems propagating few, complex phenomena with limited feature width (shocks, waves, vortices, etc.) What about flows with less distinct/obvious features? Near-turbulent or onset-turbulent flow Rotorcraft Or PDE problems from other fields? Solid Mechanics (e.g., elasticity) Maxwell s Equations, Magnetohydrodynamics Shallow water (e.g., climate problems) E. Liu (MIT) 1st Committee Meeting 30 / 41
31 Outline 1 Introduction and Motivation 2 Existing Work 3 Expected Contributions 4 Logistics E. Liu (MIT) 1st Committee Meeting 31 / 41
32 Proposed Timeline Time Task Now - 09/12 Develop, defend thesis proposal 08/12-12/12 Improve preconditioning 1/13-4/13 Improve shock capturing implementation (robustness, error estimation) 5/13-10/13 Apply method to complex 2D + 1 problems (time allowing) Investigate method effectiveness on nearturbulent flow 11/13-4/14 Write, defend thesis E. Liu (MIT) 1st Committee Meeting 32 / 41
33 Completed Coursework Major: Aerospace Computational Engineering 2.25 Advanced Fluid Mechanics A Optimization Methods A * Flight Vehicle Aerodynamics A Aerodynamics of Viscous Fluids A Computational Mechanics of Materials B * Numerical Methods for PDEs A Adv. Topics in Numerical Methods for PDEs A Introduction to Numerical Methods A Numerical Methods for PDEs A+ Minor: Theoretical Computer Science Geometric Folding Algorithms A Advanced Data Structures A Advanced Algorithms A Randomized Algorithms A Topics in Theoretical CS A * Indicates that the course was taken as an undergraduate. E. Liu (MIT) 1st Committee Meeting 33 / 41
34 Backup Slides E. Liu (MIT) 1st Committee Meeting 34 / 41
35 Motivation: DOF Reduction Ignoring the dissipation issue... Consider feature-propagation in 1D. Let: L propagation distance, l feature width (=0 for shocks) x element size (1D), N elem,1d = l x Then the number of elements needed scales like: L l Current (isotropic) methods: Potential (anisotropic) method: Dimension Reduction x x = L l N2 elem,1d L l AR x x = 1 AR L l N2 elem,1d = Ω(N elem,1d) 1 L Generalizes to higher-dimensions: AR l N2 elem,nd+1 = Ω(N elem,nd) Dimension reduction saves DOFs by exploiting higher AR to reduce the effect of L x E. Liu (MIT) 1st Committee Meeting 35 / 41
36 Euler Equation Implementation The Euler Equations can be written: F(U(x)) = F ij x j = 0, x Ω where Ω is a space-time domain and = ( t, x 1, x 2, x 3 ). F : R 5 R 5 4 is defined as: ρ ρu 1 ρu 2 ρu 3 ρu 1 ρu 2 1 F = + p ρu 1u 2 ρu 1 u 3 ρu 2 ρu 1 u 2 ρu p ρu 2u 3 ρu 3 ρu 1 u 3 ρu 2 u 3 ρu p ρe (ρe + p)u 1 (ρe + p)u 2 (ρe + p)u 3 Roe Flux Requirements Reduces to full-upwinding for faces with normal (t, 0, 0, 0) Reduces to standard Roe for faces with normal (0, x 1, x 2, x 3 ) E. Liu (MIT) 1st Committee Meeting 36 / 41
37 Euler Equation Implementation Standard Roe F Roe = 1 2 (F n,r + F n,l ) 1 2 A Roe (U R U L ) with eigenvalues v n ± a and v n. Space-Time Roe Apply Roe s method to Un 0 + F spatial n s Roe-averaging unchanged (for normal velocities, etc. use n s, the spatial normal, which need not have unit length) Eigenvalues become v n ± a n s + n t and v n + n t Validation Verified that Riemann invariants were preserved along characteristic lines and propagated at the correct speed. We will show the entropy mode (1D) and shear mode (2D) E. Liu (MIT) 1st Committee Meeting 37 / 41
38 Preconditioning: Long Journey Ahead The cost issue is severe. Let s look at some (rough) numbers: Number of Basis Functions 2D 3D 4D 1 2 (p + 1)(p + 2) 1 6 (p + 1)(p + 2)(p + 3) 1 24 (p + 1)... (p + 4) (p = 0) (p = 1) (p = 2) (p = 3) (p = 4) (p = 5) E. Liu (MIT) 1st Committee Meeting 38 / 41
39 Preconditioning: Long Journey Ahead Currently, we run 3D cases at... p = 2 or p = 3 100K p = 2 elements or 1M DOFs Insufficient for RANS problems Potentially workable for Euler shock problems Consider memory usage for p = 2 primal (p = 3 dual) solves: Memory Usage for 3D and 4D (Jacobian+GMRES) N elem p = 2, 3D p = 3, 3D p = 2, 4D p = 3, 4D 10K 1.6GB 5.2GB 44.7GB 229GB 30K 5.0GB 15.6GB 134GB 688GB 100K 16.8GB 52.2GB 447GB 2.30TB 300K 50.3GB 156GB 1.3TB 6.9TB 1M 167GB 522GB 4.5TB 23TB For reference, the upcoming ACDL cluster will have 40 32GB = 1.28TB GB or about 50,000 4D elements. E. Liu (MIT) 1st Committee Meeting 39 / 41
40 3D + 1 Capability Interesting engineering problems are essentially all fully 3D. 3D + 1 adaptive space-time simulations pose meshing challenges: Anisotropic, metric-based simplex meshing in 4D Ability to handle internal (material) boundaries Visualization Impossible to even look at the entire mesh Even in 3D, looking at 2D mesh slices is minimally informative Very challenging problem! Simplify by considering simple geometry only. Two current ideas for meshing: Generalize method of Loseille (edge-split primitives) to 4D 4D Anisotropic Delaunay (disadvantage: leaves no path toward more complex geometries) Depending on time, may stick with solving more complex 2D problems; e.g., stator-rotor flows E. Liu (MIT) 1st Committee Meeting 40 / 41
41 Disk Space Requirements: Space-Time vs. Method of Lines Memory Usage for 3D and 4D Solution Vector N elem p = 2, 3D p = 3, 3D p = 2, 4D p = 3, 4D 1M 0.447GB 0.894GB 2.68GB 6.26GB 1.8M 0.804GB 1.61GB 4.83GB 11.2GB 10M 4.47GB 8.94GB 26.8GB 62.6GB 100M 44.7GB 89.4GB 268GB 626GB New ACDL cluster has enough memory for at most p = 3 adjoint on 1.8M elements in 3D. Require checkpoints for unsteady adjoint solves 80GB to 800GB storage for p = 2, 3D, 1.8M elem Checkpoints are subset of 10 4 to 10 5 time-steps Space-time solution for the same problem size requires 500GB Assumes 100 time-spanning elements (= 100M total) But ACDL cluster may be unable to solve such a large problem E. Liu (MIT) 1st Committee Meeting 41 / 41
Is 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 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 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 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 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 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 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 informationEigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations
NASA/CR-1998-208741 ICASE Report No. 98-51 Eigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations David L. Darmofal Massachusetts Institute of Technology, Cambridge,
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 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 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 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 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 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 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 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 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 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 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 informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
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 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 informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
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 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 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 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 informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
More informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
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 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 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 informationAPPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS
TASK QUARTERLY Vol. 17, Nos 3 4, 2013, pp. 145 154 APPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS SŁAWOMIR DYKAS, WŁODZIMIERZ WRÓBLEWSKI AND SEBASTIAN RULIK
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
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 informationMixing-Plane Method for Flutter Computation in Multi-stage Turbomachines
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 2009, Orlando, Florida AIAA 2009-862 Mixing-Plane Method for Flutter Computation in Multi-stage
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 informationLARGE EDDY SIMULATION OF FLOW OVER NOZZLE GUIDE VANE OF A TRANSONIC HIGH PRESSURE TURBINE
20 th Annual CFD Symposium, August 09-10, 2018, Bangalore LARGE EDDY SIMULATION OF FLOW OVER NOZZLE GUIDE VANE OF A TRANSONIC HIGH PRESSURE TURBINE Bharathan R D, Manigandan P, Vishal Tandon, Sharad Kapil,
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 informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
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 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 informationParallel Computations of Unsteady Three-Dimensional Flows in a High Pressure Turbine
Parallel Computations of Unsteady Three-Dimensional Flows in a High Pressure Turbine Dongil Chang and Stavros Tavoularis Department of Mechanical Engineering, University of Ottawa, Ottawa, ON Canada Stavros.Tavoularis@uottawa.ca
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 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 informationChallenges and Complexity of Aerodynamic Wing Design
Chapter 1 Challenges and Complexity of Aerodynamic Wing Design Kasidit Leoviriyakit and Antony Jameson Department of Aeronautics and Astronautics Stanford University, Stanford CA kasidit@stanford.edu and
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 informationModeling Unsteady Flow in Turbomachinery Using a Harmonic Balance Technique
Modeling Unsteady Flow in Turbomachinery Using a Harmonic Balance Technique Torsten Palenschat 317220 30.04.2014 Master seminar CES WS 2013/2014 Center for Computational Engineering Sciences RWTH Aachen
More informationAssessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow
Assessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow Simone Colonia, René Steijl and George N. Barakos CFD Laboratory - School of Engineering - University
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 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 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 informationNonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations p.1/50
Nonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations Matthew McMullen Advisor: Antony Jameson Co-advisor: Juan Alonso Sponsor: Accelerated Strategic Computing Initiative
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 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 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 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 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 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 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 informationTransition Modeling Activities at AS-C²A²S²E (DLR)
Transition Modeling Activities at AS-C²A²S²E (DLR) Andreas Krumbein German Aerospace Center (DLR) Institute of Aerodynamics and Flow Technology (AS) C²A²S²E - Center for Computer Applications in AeroSpace
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 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 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 informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationA Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws
A Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws S. Jaisankar, S. V. Raghurama Rao CFD Centre, Department of Aerospace Engineering, Indian Institute of Science, Bangalore, India-56 Abstract
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationLarge Scale Simulations of Turbulent Flows for Industrial Applications. Lakhdar Remaki BCAM- Basque Centre for Applied Mathematics
Large Scale Simulations o Turbulent Flows or Industrial Alications Lakhdar Remaki BCAM- Basque Centre or Alied Mathematics Outline Flow Motion Simulation Physical Model Finite-Volume Numerical Method CFD
More informationNEWTON-GMRES PRECONDITIONING FOR DISCONTINUOUS GALERKIN DISCRETIZATIONS OF THE NAVIER-STOKES EQUATIONS
NEWTON-GMRES PRECONDITIONING FOR DISCONTINUOUS GALERKIN DISCRETIZATIONS OF THE NAVIER-STOKES EQUATIONS P.-O. PERSSON AND J. PERAIRE Abstract. We study preconditioners for the iterative solution of the
More informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
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 informationProspects for High-Speed Flow Simulations
Prospects for High-Speed Flow Simulations Graham V. Candler Aerospace Engineering & Mechanics University of Minnesota Support from AFOSR and ASDR&E Future Directions in CFD Research: A Modeling & Simulation
More informationA Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes
A Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory SIAM Annual Meeting
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 informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
More informationSimulation of unsteady muzzle flow of a small-caliber gun
Advances in Fluid Mechanics VI 165 Simulation of unsteady muzzle flow of a small-caliber gun Y. Dayan & D. Touati Department of Computational Mechanics & Ballistics, IMI, Ammunition Group, Israel Abstract
More informationTurbomachinery Applications with the Time Spectral Method
35th AIAA Fluid Dynamics Conference and Exhibit, June 6-9,005, Toronto,Ontario Turbomachinery Applications with the Time Spectral Method Edwin van der Weide, Arathi K Gopinath and Antony Jameson Stanford
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 informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
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 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 informationGT UNSTEADY SIMULATION OF A TWO-STAGE COOLED HIGH PRESSURE TURBINE USING AN EFFICIENT NON-LINEAR HARMONIC BALANCE METHOD
Proceedings of ASME Turbo Expo 213: Turbine Technical Conference and Exposition GT213 June 3-7, 213, San Antonio, Texas, USA GT213-94574 UNSTEADY SIMULATION OF A TWO-STAGE COOLED HIGH PRESSURE TURBINE
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 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 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 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 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 informationAffordable, entropy-consistent, Euler flux functions
Affordable, entropy-consistent, Euler flux functions (with application to the carbuncle phenomenon) Phil Roe Aerospace Engineering University 0f Michigan Ann Arbor Presented at HYP 2006 1 Entropy pairs
More informationCourant Institute of Mathematical Sciences, New York University, New York, NY 10012,
EFFICIENT VARIABLE-COEFFICIENT FINITE-VOLUME STOKES SOLVERS MINGCHAO CAI, ANDY NONAKA, JOHN B. BELL, BOYCE E. GRIFFITH, AND ALEKSANDAR DONEV Abstract. We investigate several robust preconditioners 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 informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationModel reduction of parametrized aerodynamic flows: discontinuous Galerkin RB empirical quadrature procedure. Masayuki Yano
Model reduction of parametrized aerodynamic flows: discontinuous Galerkin RB empirical quadrature procedure Masayuki Yano University of Toronto Institute for Aerospace Studies Acknowledgments: Anthony
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 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 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 informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More information