1st International Workshop on High-Order CFD Methods
|
|
- Kathleen Francis
- 5 years ago
- Views:
Transcription
1 1st International Worksop on Hig-Order CFD Metods Test Case C1.3: Steady Flow over te NACA0012 Airfoil Georg May, Micael Woopen Graduate Scool AICES, RWTH Aacen January 26, 2012 Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
2 Test Case 1.3: Summary Trends in numerical discretization and relaxation procedures Summary of results Tis report: Not too eavy on statistics (small sample size) 1.3a (subsonic, inviscid): 6 data sets 1.3b (transonic, inviscid): 3 data sets 1.3c (subsonic, viscous): 5 data sets Note: Not all contributions represented (missing data) Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
3 Discretization Metods Metod of coice seems to be DG Tis includes many DG-flavors (e.g. CPR-DG) Non-DG submissions: Finite-Volume (1) SBP-SAT (1) Most use fairly standard DG scemes Winner of popularity contest: DG wit BR2 (used by a tird of submissions) Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
4 Meses and Fidelity of Results Meses Very many different meses found use provided quad meses self-generated meses (quad, tri, and mixed) Correlation between mes type and fidelity of results not attempted ere Submissions include adaptive metods (goal-oriented) General Remark on Results For subsonic cases: Convergence studies reveal superiority of ig-order approximation Fewer DOFS for same error Lower CPU time for same error Work needs to be done for transonic flow Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
5 Solution Metods Tis is a steady test case! Look for efficient relaxation metods Metod of coice seems to be (damped) Newton / Krylov Almost all use GMRES (one BICGSTAB) No consensus on preconditioners ILU-n Gauss-Seidel Multigrid Important note: We do not compare and evaluate efficiency insufficient number of samples We cannot assume all codes are fully optimized Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
6 Convergence Summary Test Case C1.3a (Subsonic inviscid) Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
7 Convergence Towards Trut Solution Do we agree on te trut solution? Trut lift and drag coefficients Average Standard Deviation Drag E E-06 Lift 2.865E E-05 Compare to p=3 solutions Average Standard Deviation Drag E E-05 Lift 2.865E E-05 Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
8 Convergence towards Trut Solution order Mean lift and drag coefficients for eac order on te finest mes approacing te mean reference values Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
9 Convergence Test Case C1.3a (Subsonic Inviscid) 2nd Order 2nd Order (p=1) 2nd Order Gooc Darmofal 2nd Order nd Order Darmofal Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
10 Convergence Test Case C1.3a (Subsonic Inviscid) 3rd Order p=2 3rd Order Darmofal Gooc Darmofal Darmofal 3nd Order rd Order Darmofal Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
11 Convergence Test Case C1.3a (Subsonic Inviscid) 4t Order p=3 4t Order Gooc Svaerd Darmofal Svaerd Svaerd 4t Order t Order Darmofal Svaerd Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
12 Convergence Test Case C1.3a (Subsonic Inviscid) 5t Order p=4 5t Order Svaerd Svaerd 5t Order t Order Svaerd Svaerd Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
13 Convergence Test Case C1.3a (Subsonic Inviscid) 6t Order p=5 6t Order Svaerd Svaerd 6t Order t Order Svaerd Svaerd Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
14 Convergence Summary Test Case C1.3b (Transonic Inviscid) Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
15 Convergence Test Case C1.3b (Transonic Inviscid) 2nd Order 2nd order (p=1) 2nd Order Darmofal Darmofal 2nd Order nd Order ( uses p=1..4. Tis curve is te same for all values of p) Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
16 Convergence Test Case C1.3b (Transonic Inviscid) 3rd Order p=2 3rd Order Darmofal Darmofal 3rd Order rd Order ( uses p=1..4. Tis curve is te same for all values of p) Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
17 Convergence Summary Test Case C1.3c (Subsonic Viscous) Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
18 Convergence Towards Trut Solution order Mean lift and drag coefficients for eac order on te finest mes approacing te mean reference values Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
19 Convergence Test Case C1.3c (Subsonic Viscous) st Order 1st order (p=0) 1st Order 1st Order st Order c l Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
20 Convergence Test Case C1.3c (Subsonic Viscous) 2nd Order p=1 2nd Order 10 1 Darmofal 10 0 Mavriplis Mavriplis Mavriplis 2nd Order nd Order Mavriplis ( uses p=1..5. Tis curve is te same for all values of p) Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
21 Convergence Test Case C1.3c (Subsonic Viscous) 3nd Order p=2 3rd Order 3rd Order rd Order ( uses p=1..5. Tis curve is te same for all values of p) Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
22 Convergence Test Case C1.3c (Subsonic Viscous) 4t Order p=3 4t Order 4t Order t Order ( uses p=1..5. Tis curve is te same for all values of p) Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
23 Convergence Test Case C1.3c (Subsonic Viscous) t Order p=4 5t Order 5t Order t Order ( uses p=1..5. Tis curve is te same for all values of p) Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
24 Convergence Test Case C1.3c (Subsonic Viscous) t Order p=5 6t Order 6t Order t Order ( uses p=1..5. Tis curve is te same for all values of p) Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
25 Summary of Conclusions Participants use: Predominantly DG Discretization metods Similar Solution metods for te steady problem Higer order pays in terms of work units versus error i.e. for te subsonic cases! Convergence in lift and drag sows considerable scatter Adaptive mes refinement demonstrates advantages Very difficult to establis clear trends Georg May, Micael Woopen 1st International Worksop on Hig-Order CFD Metods
Hp-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 informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
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 information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
More informationA Polynomial Adaptive LCP Scheme for Viscous Compressible Flows
A Polynomial Adaptive LCP Sceme for Viscous Compressible Flows J.S. Cagnone, B.C. Vermeire, and S.. Nadaraja Department of Mecanical Engineering, McGill University, Montreal, Canada, H3A 2S6 Email: jean-sebastien.cagnone@mail.mcgill.ca
More informationCompressible Potential Flow: The Full Potential Equation. Copyright 2009 Narayanan Komerath
Compressible Potential Flow: The Full Potential Equation 1 Introduction Recall that for incompressible flow conditions, velocity is not large enough to cause density changes, so density is known. Thus
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 informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 27, No. 4, pp. 47 492 c 26 Society for Industrial and Applied Matematics A NOVEL MULTIGRID BASED PRECONDITIONER FOR HETEROGENEOUS HELMHOLTZ PROBLEMS Y. A. ERLANGGA, C. W. OOSTERLEE,
More information(4.2) -Richardson Extrapolation
(.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as
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 informationEdwin van der Weide and Magnus Svärd. I. Background information for the SBP-SAT scheme
Edwin van der Weide and Magnus Svärd I. Background information for the SBP-SAT scheme As is well-known, stability of a numerical scheme is a key property for a robust and accurate numerical solution. Proving
More 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 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 informationHow Many Steps are Required to Solve the Euler Equations of Steady, Compressible Flow: In Search of a Fast Solution Algorithm
AIAA 2 2673 How Many Steps are Required to Solve the Euler Equations of Steady, Compressible Flow: In Search of a Fast Solution Algorithm Antony Jameson Stanford University Stanford, CA 9435 D. A. Caughey
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationA Weak Galerkin Method with an Over-Relaxed Stabilization for Low Regularity Elliptic Problems
J Sci Comput (07 7:95 8 DOI 0.007/s095-06-096-4 A Weak Galerkin Metod wit an Over-Relaxed Stabilization for Low Regularity Elliptic Problems Lunji Song, Kaifang Liu San Zao Received: April 06 / Revised:
More informationMultigrid Methods for Discretized PDE Problems
Towards Metods for Discretized PDE Problems Institute for Applied Matematics University of Heidelberg Feb 1-5, 2010 Towards Outline A model problem Solution of very large linear systems Iterative Metods
More informationIntroduction to Multigrid Method
Introduction to Multigrid Metod Presented by: Bogojeska Jasmina /08/005 JASS, 005, St. Petersburg 1 Te ultimate upsot of MLAT Te amount of computational work sould be proportional to te amount of real
More informationThe Derivative as a Function
Section 2.2 Te Derivative as a Function 200 Kiryl Tsiscanka Te Derivative as a Function DEFINITION: Te derivative of a function f at a number a, denoted by f (a), is if tis limit exists. f (a) f(a + )
More informationLimit Cycle Oscillations of a Typical Airfoil in Transonic Flow
Limit Cycle Oscillations of a Typical Airfoil in Transonic Flow Denis B. Kholodar, United States Air Force Academy, Colorado Springs, CO 88 Earl H. Dowell, Jeffrey P. Thomas, and Kenneth C. Hall Duke University,
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationA Hybridizable Discontinuous Galerkin Method for the Compressible Euler and Navier-Stokes Equations
A Hybridizable Discontinuous Galerkin Metod for te Compressible Euler and Navier-Stokes Equations J. Peraire and N. C. Nguyen Massacusetts Institute of Tecnology, Cambridge, MA 02139, USA B. Cockburn University
More informationExplicit Hyperbolic Reconstructed Discontinuous Galerkin Methods for Time-Dependent Problems
AIAA AVIATION Forum June 25-29 218 Atlanta Georgia 218 Fluid Dynamics Conference 1.2514/6.218-427 Explicit Hyperbolic Reconstructed Discontinuous Galerkin Metods for Time-Dependent Problems Jialin Lou
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More information, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim
Mat 132 Differentiation Formulas Stewart 2.3 So far, we ave seen ow various real-world problems rate of cange and geometric problems tangent lines lead to derivatives. In tis section, we will see ow to
More informationSmoothed aggregation for Helmholtz problems
NUMERICL LINER LGEBR WITH PPLICTIONS Numer. Linear lgebra ppl. 2010; 17:361 386 Publised online 2 February 2010 in Wiley InterScience (www.interscience.wiley.com)..686 Smooted aggregation for Helmoltz
More informationScalable parallel approach for high-fidelity steady-state aeroelastic analysis and adjoint derivative computations
This is a preprint of the following article, which is available from http://mdolab.engin.umich.edu Gaetan K. W. Kenway, Graeme J. Kennedy, and Joaquim R. R. A. Martins, Scalable Parallel Approach for High-
More informationSimulations of the turbulent channel flow at Re τ = 180 with projection-based finite element variational multiscale methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Met. Fluids 7; 55:47 49 Publised online 4 Marc 7 in Wiley InterScience (www.interscience.wiley.com). DOI:./fld.46 Simulations of te
More informationNumerical Simulation and Aerodynamic Energy Analysis of Limit Cycle Flutter of a Bridge Deck
Proceedings of te 018 World Transport Convention Beijing, Cina, June 18-1, 018 Numerical Simulation and Aerodynamic Energy Analysis of Limit Cycle Flutter of a Bridge Deck Xuyong Ying State Key Laboratory
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 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 informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
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 informationOptimal iterative solvers for linear nonsymmetric systems and nonlinear systems with PDE origins: Balanced black-box stopping tests
Optimal iterative solvers for linear nonsymmetric systems and nonlinear systems wit PDE origins: Balanced black-box stopping tests Pranjal, Prasad and Silvester, David J. 2018 MIMS EPrint: 2018.13 Mancester
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationAlgebraic multigrid within defect correction for the linearized Euler equations
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2010; 17:307 324 Published online 10 December 2009 in Wiley InterScience (www.interscience.wiley.com)..687 Algebraic multigrid within
More informationIterative Methods and Multigrid
Iterative Methods and Multigrid Part 3: Preconditioning 2 Eric de Sturler Preconditioning The general idea behind preconditioning is that convergence of some method for the linear system Ax = b can be
More informationConvergence and Descent Properties for a Class of Multilevel Optimization Algorithms
Convergence and Descent Properties for a Class of Multilevel Optimization Algoritms Stepen G. Nas April 28, 2010 Abstract I present a multilevel optimization approac (termed MG/Opt) for te solution of
More informationMAT Calculus for Engineers I EXAM #1
MAT 65 - Calculus for Engineers I EXAM # Instructor: Liu, Hao Honor Statement By signing below you conrm tat you ave neiter given nor received any unautorized assistance on tis eam. Tis includes any use
More informationarxiv: v1 [math.na] 3 Nov 2011
arxiv:.983v [mat.na] 3 Nov 2 A Finite Difference Gost-cell Multigrid approac for Poisson Equation wit mixed Boundary Conditions in Arbitrary Domain Armando Coco, Giovanni Russo November 7, 2 Abstract In
More informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
More informationGRID CONVERGENCE ERROR ANALYSIS FOR MIXED-ORDER NUMERICAL SCHEMES
GRID CONVERGENCE ERROR ANALYSIS FOR MIXED-ORDER NUMERICAL SCHEMES Cristoper J. Roy Sandia National Laboratories* P. O. Box 5800, MS 085 Albuquerque, NM 8785-085 AIAA Paper 00-606 Abstract New developments
More informationRecent Advances in Time-Domain Maxwell s Equations in Metamaterials
Recent Advances in Time-Domain Maxwell s Equations in Metamaterials Yunqing Huang 1, and Jicun Li, 1 Hunan Key Laboratory for Computation and Simulation in Science and Engineering, Xiangtan University,
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 informationAN EFFICIENT AND ROBUST METHOD FOR SIMULATING TWO-PHASE GEL DYNAMICS
AN EFFICIENT AND ROBUST METHOD FOR SIMULATING TWO-PHASE GEL DYNAMICS GRADY B. WRIGHT, ROBERT D. GUY, AND AARON L. FOGELSON Abstract. We develop a computational metod for simulating models of gel dynamics
More informationLyapunov spectrum of scale-resolving turbulent simulations. Application to chaotic adjoints
Lyapunov spectrum of scale-resolving turbulent simulations. Application to caotic adjoints Te MIT Faculty as made tis article openly available. Please sare ow tis access benefits you. Your story matters.
More informationThe use of high-order methods to compute turbulent flows governed by the Reynolds-averaged Navier-
AIAA SciTech 5-9 January 2015, Kissimmee, Florida 53rd AIAA Aerospace Sciences Meeting AIAA 2015-1286 CPR High-order Discretization of the RANS Equations with the SA Model Cheng Zhou and Z.J. Wang The
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 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 informationarxiv: v1 [math.na] 7 Mar 2019
Local Fourier analysis for mixed finite-element metods for te Stokes equations Yunui He a,, Scott P. MacLaclan a a Department of Matematics and Statistics, Memorial University of Newfoundland, St. Jon
More informationHP-MULTIGRID AS SMOOTHER ALGORITHM FOR HIGHER ORDER DISCONTINUOUS GALERKIN DISCRETIZATIONS OF ADVECTION DOMINATED FLOWS PART I. MULTILEVEL ANALYSIS
HP-MULTIGRID AS SMOOTHER ALGORITHM FOR HIGHER ORDER DISCONTINUOUS GALERKIN DISCRETIZATIONS OF ADVECTION DOMINATED FLOWS PART I. MULTILEVEL ANALYSIS J.J.W. VAN DER VEGT AND S. RHEBERGEN AMS subect classifications.
More informationUNIT #6 EXPONENTS, EXPONENTS, AND MORE EXPONENTS REVIEW QUESTIONS
Answer Key Name: Date: UNIT # EXPONENTS, EXPONENTS, AND MORE EXPONENTS REVIEW QUESTIONS Part I Questions. Te epression 0 can be simpliied to () () 0 0. Wic o te ollowing is equivalent to () () 8 8? 8.
More informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
More informationCalculus I Practice Exam 1A
Calculus I Practice Exam A Calculus I Practice Exam A Tis practice exam empasizes conceptual connections and understanding to a greater degree tan te exams tat are usually administered in introductory
More informationThe Advection-Diffusion equation!
ttp://www.d.edu/~gtryggva/cf-course/! Te Advectio-iffusio equatio! Grétar Tryggvaso! Sprig 3! Navier-Stokes equatios! Summary! u t + u u x + v u y = P ρ x + µ u + u ρ y Hyperbolic part! u x + v y = Elliptic
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationA fully-coupled discontinuous Galerkin method for two-phase flow in porous media with discontinuous capillary pressure
arxiv:39.7555v [pysics.flu-dyn] 29 Sep 23 A fully-coupled discontinuous Galerkin metod for two-pase flow in porous media wit discontinuous capillary pressure Peter Bastian Heidelberg University, Interdisciplinary
More informationDirect Numerical Solution of the Steady 1D Compressible Euler Equations for Transonic Flow Profiles with Shocks
Direct Numerical Solution of the Steady 1D Compressible Euler Equations for Transonic Flow Profiles with Shocks Hans De Sterck, Scott Rostrup Department of Applied Mathematics, University of Waterloo,
More informationSection 2.1 The Definition of the Derivative. We are interested in finding the slope of the tangent line at a specific point.
Popper 6: Review of skills: Find tis difference quotient. f ( x ) f ( x) if f ( x) x Answer coices given in audio on te video. Section.1 Te Definition of te Derivative We are interested in finding te slope
More informationMultigrid Methods for Obstacle Problems
Multigrid Metods for Obstacle Problems by Cunxiao Wu A Researc Paper presented to te University of Waterloo in partial fulfillment of te requirement for te degree of Master of Matematics in Computational
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationHanieh Mirzaee Jennifer K. Ryan Robert M. Kirby
J Sci Comput (2010) 45: 447 470 DOI 10.1007/s10915-009-9342-9 Quantification of Errors Introduced in te Numerical Approximation and Implementation of Smootness-Increasing Accuracy Conserving (SIAC) Filtering
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 informationCONVERGENCE ANALYSIS OF YEE SCHEMES FOR MAXWELL S EQUATIONS IN DEBYE AND LORENTZ DISPERSIVE MEDIA
INTRNATIONAL JOURNAL OF NUMRICAL ANALYSIS AND MODLING Volume XX Number 0 ages 45 c 03 Institute for Scientific Computing and Information CONVRGNC ANALYSIS OF Y SCHMS FOR MAXWLL S QUATIONS IN DBY AND LORNTZ
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationEffect of Small-Scale Output Unsteadiness on Adjoint-Based Sensitivity
AIAA JOURNAL Vol. 48, No. 11, November 21 Effect of Small-Scale Output Unsteadiness on Adjoint-Based Sensitivity Joshua A. Krakos and David L. Darmofal Massachusetts Institute of Technology, Cambridge,
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationCFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow Neale, A.; Derome, D.; Blocken, B.; Carmeliet, J.E.
CFD calculation of convective eat transfer coefficients and validation Part I: Laminar flow Neale, A.; Derome, D.; Blocken, B.; Carmeliet, J.E. Publised in: IEA Annex 41 working meeting, Kyoto, Japan Publised:
More informationOptimizing Runge-Kutta smoothers for unsteady flow problems
Optimizing Runge-Kutta smoothers for unsteady flow problems Philipp Birken 1 November 24, 2011 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany. email:
More informationSupersonic Aerodynamics. Methods and Applications
Supersonic Aerodynamics Methods and Applications Outline Introduction to Supersonic Flow Governing Equations Numerical Methods Aerodynamic Design Applications Introduction to Supersonic Flow What does
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
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 informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
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 informationSMAI-JCM SMAI Journal of Computational Mathematics
SMAI-JCM SMAI Journal of Computational Matematics Compatible Maxwell solvers wit particles II: conforming and non-conforming 2D scemes wit a strong Faraday law Martin Campos Pinto & Eric Sonnendrücker
More informationc 2003 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 4, No. 4, pp. 388 45 c 23 Society for Industrial and Applied Mathematics CONVERGENCE ANALYSIS OF THE GAUSS SEIDEL PRECONDITIONER FOR DISCRETIZED ONE DIMENSIONAL EULER EQUATIONS
More informationTHE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein
Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationJian-Guo Liu 1 and Chi-Wang Shu 2
Journal of Computational Pysics 60, 577 596 (000) doi:0.006/jcp.000.6475, available online at ttp://www.idealibrary.com on Jian-Guo Liu and Ci-Wang Su Institute for Pysical Science and Tecnology and Department
More informationA parallel Newton-Krylov-Schur flow solver for the Reynolds-Averaged Navier-Stokes equations
A parallel Newton-Krylov-Schur flow solver for the Reynolds-Averaged Navier-Stokes equations Michal Osusky and David W. Zingg Institute for Aerospace Studies, University of Toronto, Toronto, Ontario, M3H
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationA New Fifth Order Finite Difference WENO Scheme for Hamilton-Jacobi Equations
A New Fift Order Finite Difference WENO Sceme for Hamilton-Jacobi Equations Jun Zu, 1 Jianxian Qiu 1 College of Science, Nanjing University of Aeronautics and Astronautics, Nanjing, Jiangsu 10016, People
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 informationModule 4 Linear Boundary value problems
Module 4 Linear Boundary value problems Lecture Content Hours 1 Finite Difference Metods: Diriclet type boundary 1 condition Finite Difference Metods: Mixed boundary condition 1 3 Sooting Metod 1 4 Sooting
More informationInf sup testing of upwind methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Met. Engng 000; 48:745 760 Inf sup testing of upwind metods Klaus-Jurgen Bate 1; ;, Dena Hendriana 1, Franco Brezzi and Giancarlo
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationRP 2.4. SEG/Houston 2005 Annual Meeting 1513
P 2.4 Measurement of sear ave velocity of eavy oil De-ua Han, Jiajin Liu, University of Houston Micael Batzle, Colorado Scool of Mines Introduction It is ell knon tat te fluids ave no sear modulus and
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 informationEFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS
Statistica Sinica 24 2014, 395-414 doi:ttp://dx.doi.org/10.5705/ss.2012.064 EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS Jun Sao 1,2 and Seng Wang 3 1 East Cina Normal University,
More informationArbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations
Arbitrary order exactly divergence-free central discontinuous Galerkin metods for ideal MHD equations Fengyan Li, Liwei Xu Department of Matematical Sciences, Rensselaer Polytecnic Institute, Troy, NY
More informationThird order Approximation on Icosahedral Great Circle Grids on the Sphere. J. Steppeler, P. Ripodas DWD Langen 2006
Tird order Approximation on Icosaedral Great Circle Grids on te Spere J. Steppeler, P. Ripodas DWD Langen 2006 Deasirable features of discretisation metods on te spere Great circle divisions of te spere:
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
More informationHigh-Order Energy and Linear Momentum Conserving Methods for the Klein-Gordon Equation
matematics Article Hig-Order Energy and Linear Momentum Conserving Metods for te Klein-Gordon Equation He Yang Department of Matematics, Augusta University, Augusta, GA 39, USA; yang@augusta.edu; Tel.:
More informationNumerical Study To Assess Sulfur Hexafluoride as a Medium for Testing Multielement Airfoils
NASA Technical Paper 3496 Numerical Study To Assess Sulfur Hexafluoride as a Medium for Testing Multielement foils Daryl L. Bonhaus and W. Kyle Anderson Langley Research Center Hampton, Virginia Dimitri
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 informationDownloaded 01/08/14 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 3, No. 3, pp. 506 56 c 00 Society for Industrial and Applied Matematics NESTED ITERATION AND FIRST-ORDER SYSTEM LEAST SQUARES FOR INCOMPRESSIBLE, RESISTIVE MAGNETOHYDRODYNAMICS
More information27. A Non-Overlapping Optimized Schwarz Method which Converges with Arbitrarily Weak Dependence on h
Fourteent International Conference on Domain Decomposition Metods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 003 DDM.org 7. A Non-Overlapping Optimized Scwarz Metod wic Converges
More information