First, Second, and Third Order Finite-Volume Schemes for Diffusion
|
|
- Neil Norris
- 6 years ago
- Views:
Transcription
1 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.
2 All Created Equal all men are created equal, - Thomas Jefferson, The Declaration of Independence Are all PDEs created equal?
3 Hyperbolic vs Parabolic Hyperbolic (inviscid): - Principle of upwinding (dissipation) led to many useful schemes. - Robust 1st-order schemes - a home to come back. - A variety of unstructured, high-order schemes. Parabolic (viscous): - Lack of universal guiding principles. - Lack of robust 1st-order schemes (rely on inconsistent scheme). - Much less variety for unstructured, high-order schemes. - Degraded accuracy of derivatives on irregular grids. They don t seem created equal...
4 Who Created PDEs? all men are created equal, that they are endowed by their Creator with certain unalienable Rights, - Thomas Jefferson, The Declaration of Independence We created PDEs. Then, we can recreate them equal.
5 Recreate Them Hyperbolic First-Order Hyperbolic System Method JCP2007, 2010, 2012, AIAA2009, 2011, 2013 Methods for hyperbolic systems apply to all PDEs. Dramatic simplification/improvements to numerical methods
6 Turn Every Food into a Burger Simple, Efficient, Accurate. Sushi Burger! It looks eccentric, but the taste is the same, or even better.
7 Nishikawa, JCP2007 Hyperbolic Diffusion System t u = ν ( xx u + yy u) Sushi Burger for Diffusion steady t u = ν ( x p + y q) t p = ( x u p)/t r t q = ( y u q)/t r 0 = ν ( x p + y q), p = x u, q = y u, 0=ν ( xx u + yy u), System is equivalent to diffusion in the steady state for any Tr : T r = L2 r ν, L r = 1 2π Unsteady computation possible by dual-time formulation (implicit) with a steady solver used in the inner iteration.
8 Nishikawa, JCP2007 Hyperbolic Diffusion System Sushi Burger for Diffusion t U + x F + y G = S, U = u p, F = νp u/t r, G = νq 0, S = 0 p/t r q 0 u/t r q/t r Normal Jacobian: Real eigenvalues: A n = H U = (Fn x + Gn y ) U ν λ 1 =, λ 2 = T r ν T r, λ 3 =0
9 Energy Estimate Integrate over the domain l E =(u, L 2 rp, L 2 rq) x (hyperbolic system): d dt Ω steady EdV = ν Ω (p 2 + q 2 ) dv Ω f E n da E = u 2 + L 2 r(p 2 + q 2 ) f E =( νup, νuq) which reduces to the energy estimate for the Laplace equation: 0= Ω u udv Energy estimate consistent with steady diffusion problem Ω u u n da
10 Diffusion is Hyperbolic Hyperbolic Advection-Diffusion System (JCP2010) Hyperbolic Navier-Stokes System (AIAA ) If you have a good inviscid scheme, you have a very good viscous scheme. - Equal order of accuracy for solution and derivatives. - O(h) time step; O(1/h) condition number.
11 Edge-Based Finite-Volume Method NASA s FUN3D; Software Cradle s SC/Tetra; DLR Tau code, etc. Edge-based finite-volume scheme: k V j du j dt = k {k j } Φ jk A jk + S j V j A jk = n l jk + n r jk n l jk n r jk with the upwind flux at edge midpoint: j Φ jk = 1 2 (H L + H R ) 1 2 A n (U R U L ) Accuracy is determined by the left and right states.
12 First-Order Scheme Home sweet home for diffusion schemes
13 First-Order Scheme for Diffusion Left and right states: R U L = U j, U R = U k j k Discrete Energy Estimate: L j {j} V j de j dt = e b {e b } f E 1 + f E 2 2 ˆn eb A eb j {j} ν(p 2 j + q 2 j )V j ν 2L r e {e} ɛ jk A jk Consistent with continuous energy estimate Dissipation ɛ jk =(u k u j ) 2 + L 2 r [(p k p j,q k q j ) ˆn jk ] 2 0 First-order upwind diffusion scheme is energy-stable for general grids.
14 Second-Order Scheme Even better
15 Second-Order Scheme For triangular/tetrahedral and smooth mixed grids. 1. Compute gradients at nodes (e.g., LSQ). 2. Extrapolate the solution to the midpoint. L R Left and right states: u L = u j (p j,q j ) l jk, u R = u k 1 2 (p k,q k ) l jk j k p L = p j p j l jk, p R = p k 1 2 p k l jk q L = q j q j l jk, q R = q k 1 2 q k l jk l jk =(x k x j,y k y j ) Gradients are not needed for the solution.
16 Taylor Expansion 5 4 steady du j dt = ν( xp + y q) νh 6L r [ ( 2+ 5) x (p x u)+ 2 y (p x u)+ 2 x (q y u)+( 2+ 5) y (q y u) νh2 12 [ xx( x p + y q)+ xy ( x p + y q)+ yy ( x p + y q)] + O(h 3 ), steady dp j dt = 1 T r ( x u p) h2 6T r steady dq j dt = 1 T r ( y u q) h2 6T r Source discretization [ xx (p x u)+ xy (p x u)+ xy (q y u)+ 1 ] 2 ( xxp + yy p + xx q) [ yy (q y u)+ xy (q y u)+ xy (p x u)+ 1 ] 2 ( xxq + yy q + xx p) j 3 ] + O(h 3 ), + O(h 3 ) Second-order accurate for solution and gradients. First-order error makes it stable with forward Euler.
17 Third-Order Scheme A new wave
18 Third-Order Scheme For triangular/tetrahedral grids only. (Katz and Sankaran JCP2011) R L j k 1. 2nd-order gradients at nodes (e.g., LSQ quadratic fit). 2. Extrapolate flux/solution to the midpoint. H L = H j H j l jk, H R = H k 1 2 H k l jk Third-order accuracy on a second-order stencil Source term needs special discretization (NIA CFD Seminar ).
19 Divergence Form of Source Write the source term at each node j as follows: S x F s + y G s F s = 0 (y y j ) q/t r, Gs = 0 (y y j ) p/t r So that (x x j ) q/t r t U + x F + y G = S (x x j ) p/t r t U + x (F F s )+ y (G G s )=0 The system is fully hyperbolic with no source terms. Source term discretization is not needed.
20 Third-Order Scheme (Katz and Sankaran JCP2011) 1. 2nd-order gradients at nodes (e.g., LSQ). 2. Extrapolate flux/solution to the midpoint. 3. Upwind flux for fully hyperbolic system. L R Left and right states: u L = u j (p j,q j ) l jk, u R = u k 1 2 (p k,q k ) l jk j k p L = p j p j l jk, p R = p k 1 2 p k l jk q L = q j q j l jk, q R = q k 1 2 q k l jk Gradients are not needed for the solution.
21 Taylor Expansion (3rd-order) steady du j dt = ν( xp + y q) νh 6L r [ ( 2+ 5) x (p x u)+ 2 y (p x u)+ 2 x (q y u)+( 2+ 5) y (q y u) νh2 12 [ xx( x p + y q)+ xy ( x p + y q)+ yy ( x p + y q)] + O(h 3 ), steady dp j dt = 1 ( x u p) h2 [( xx + xy )(q y u)+ xx (p x u)+ y ( x q y p)] + O(h 3 ), T r 6T r steady dq j dt = 1 ( y u q) h2 [( xy + yy )(p x u)+ yy (q y u) x ( x q y p)] + O(h 3 ) T r 6T r ] Third-order accurate for solution and gradients. First-order error makes it stable with forward Euler.
22 Numerical Experiment Exact solution: u(x, y) = sinh(πx)sin(πy)+sinh(πy)sin(πx) sinh(π) 1 - n x n grids: n = 9, 17, 33, 65, 129, Dirichlet boundary condition. y - 10 neighbors for quadratic fit. (to avoid ill-conditioning of LSQ matrix) - Forward Euler time stepping x - Steady state reached when residual drops below 1.0E-15 - Comparison with the Galerkin scheme
23 Max CFL Number Max CFL numbers determined by Fourier analysis: First-Order Second-Order Third-Order Forward Euler Hyperbolic schemes are stable with O(h) time step: O(1/h) faster than typical schemes with O(h^2) time step.
24 Error Convergence 1 Log 10 (L 1 error of u) u: solution p: x-derivative Galerkin 1st Order Slope 1 Slope 2 Slope Log 10 (h) Log 10 (L 1 error of p) LSQ (Galerkin) 1st Order Slope 1 Slope 2 Slope Log (h) 10 1st-order accurate solution and gradients.
25 Error Convergence 2 Log 10 (L 1 error of u) u: solution p: x-derivative Galerkin 1st Order 2nd Order Slope 1 Slope 2 Slope Log 10 (h) Log 10 (L 1 error of p) LSQ (Galerkin) 1st Order 2nd Order Slope 1 Slope 2 Slope Log 10 (h) 2nd-order accurate solution and gradients.
26 Error Convergence 3 Log 10 (L 1 error of u) u: solution p: x-derivative Galerkin 1st Order 2nd Order 3rd Order Slope 1 Slope 2 Slope Log (h) 10 Log 10 (L 1 error of p) LSQ (Galerkin) 1st Order 2nd Order 3rd Order Slope 1 Slope 2 Slope Log 10 (h) 3rd-order accurate solution and gradients.
27 Cost Comparison Cost per time step (the irregular grid case). Galerkin First-Order Second-Order Third-Order Forward Euler Almost the same. Reality is that hyperbolic schemes are more economical because they converge O(1/h) faster than typical diffusion schemes.
28 Time to Solution O(1/h) acceleration overwhelms the increased cost per time step. CPU Time Galerkin 1st Order 2nd Order 3rd Order Log 10 ( CPU Time ) Galerkin 1st Order 2nd Order 3rd Order Slope 2 Slope Number of nodes (x 10000) Log 10 ( Number of nodes ) Orders of magnitude acceleration in CPU time.
29 Conclusions Diffusion and source recreated hyperbolic. 1. Energy-stable first-order diffusion scheme 2. Equal order of accuracy for sol. and gradients on irregular grids. 3. Third-order diffusion scheme by fully hyperbolic system Third-order scheme is incomparably efficient and accurate.
30 Future Work Uniformly third-order advection-diffusion scheme. Implicit schemes. Time-dependent problems. Third-order scheme for Navier-Stokes (2nd-order scheme in AIAA2011). New system for accurate velocity gradients (vorticity, turb source). Hyperbolic formulation for turbulence models (robust diffusion). Various Other Applications: - High-order residual-distribution schemes, dispersion eq., Incomp. NS at INRIA - 3rd-order active flux scheme at University of Michigan - Entropy-consistent scheme at Universiti Sains Malaysia - Many other potential applications: DG, SV, CESE, SUPG, etc.
31 Declaration of Hyperbolicity We hold these truths to be self-evident, that all PDEs are created equal, that they are endowed by us with certain unalienable Rights, that among these are hyperbolicity, consistent and accurate schemes and the pursuit of robustness.
Divergence 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 informationActive Flux for Advection Diffusion
Active Flux for Advection Diffusion A Miracle in CFD Hiroaki Nishikawa National Institute of Aerospace! NIA CFD Seminar! August 25, 2015 In collaboration with the University of Michigan Supported by NASA
More informationAccuracy-Preserving Source Term Quadrature for Third-Order Edge-Based Discretization
Preprint accepted in Journal of Computational Physics. https://doi.org/0.06/j.jcp.07.04.075 Accuracy-Preserving Source Term Quadrature for Third-Order Edge-Based Discretization Hiroaki Nishikawa and Yi
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 informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
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 informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationAdditive Manufacturing Module 8
Additive Manufacturing Module 8 Spring 2015 Wenchao Zhou zhouw@uark.edu (479) 575-7250 The Department of Mechanical Engineering University of Arkansas, Fayetteville 1 Evaluating design https://www.youtube.com/watch?v=p
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
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 informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationChapter 10 Exercises
Chapter 10 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rl/fdmbook Exercise 10.1 (One-sided and
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
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 informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
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 informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationTAU Solver Improvement [Implicit methods]
TAU Solver Improvement [Implicit methods] Richard Dwight Megadesign 23-24 May 2007 Folie 1 > Vortrag > Autor Outline Motivation (convergence acceleration to steady state, fast unsteady) Implicit methods
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
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 informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Partial Differential Equations (Lecture 1, Week 1) Markus Schmuck Department of Mathematics and Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
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 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 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 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 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 informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
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 informationChapter 2 Finite-Difference Discretization of the Advection-Diffusion Equation
Chapter Finite-Difference Discretization of the Advection-Diffusion Equation. Introduction Finite-difference methods are numerical methods that find solutions to differential equations using approximate
More informationThe Orchestra of Partial Differential Equations. Adam Larios
The Orchestra of Partial Differential Equations Adam Larios 19 January 2017 Landscape Seminar Outline 1 Fourier Series 2 Some Easy Differential Equations 3 Some Not-So-Easy Differential Equations Outline
More informationBorel Summability in PDE initial value problems
Borel Summability in PDE initial value problems Saleh Tanveer (Ohio State University) Collaborator Ovidiu Costin & Guo Luo Research supported in part by Institute for Math Sciences (IC), EPSRC & NSF. Main
More informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationSolving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Solving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics Fikriye Yılmaz 1, Bülent Karasözen
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1
AE/ME 339 Professor of Aerospace Engineering 12/21/01 topic7_ns_equations 1 Continuity equation Governing equation summary Non-conservation form D Dt. V 0.(2.29) Conservation form ( V ) 0...(2.33) t 12/21/01
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 informationFUNDAMENTALS OF FINITE DIFFERENCE METHODS
FUNDAMENTALS OF FINITE DIFFERENCE METHODS By Deep Gupta 3 rd Year undergraduate, Mechanical Engg. Deptt., IIT Bombay Supervised by: Prof. Gautam Biswas, IIT Kanpur Acknowledgements It has been a pleasure
More informationChapter 5 Types of Governing Equations. Chapter 5: Governing Equations
Chapter 5 Types of Governing Equations Types of Governing Equations (1) Physical Classification-1 Equilibrium problems: (1) They are problems in which a solution of a given PDE is desired in a closed domain
More informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Pseudo-Time Integration 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 10 Outline 1
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationPartial Differential Equations II
Partial Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Partial Differential Equations II 1 / 28 Almost Done! Homework
More informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
More informationThird-Order Active-Flux Scheme for Advection Diffusion: Hyperbolic Diffusion, Boundary Condition, and Newton Solver
Preprint accepted for publication in Computers and Fluids, 16 http://d.doi.org/1.116/j.compfluid.15.1. Third-Order Active-Flu Scheme for Advection Diffusion: Hyperbolic Diffusion, Boundary Condition, and
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationComputational Fluid Dynamics I + II. Mathematische Strömungslehre I + II
Computational Fluid Dynamics I + II Mathematische Strömungslehre I + II Prof. Dr.-Ing. D. Hänel Aerodynamisches Institut RWTH Aachen Contents 1 Governing equations of fluid dynamics 1 1.1 General formulation
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationComputation Fluid Dynamics
Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationRunge-Kutta-Chebyshev Projection Method
Runge-Kutta-Chebyshev Projection Method Zheming Zheng Linda Petzold Department of Mechanical Engineering, University of California Santa Barbara, Santa Barbara, CA 93106, USA June 8, 2006 This work was
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
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 Contours of Pressure.8.6.4...4.6.8 Matthew McMullen Antony Jameson Juan J. Alonso Aero Astro Dept. Stanford
More informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
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 informationNumerical Analysis and Computer Science DN2255 Spring 2009 p. 1(8)
Numerical Analysis and Computer Science DN2255 Spring 2009 p. 1(8) Contents Important concepts, definitions, etc...2 Exact solutions of some differential equations...3 Estimates of solutions to differential
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 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 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 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 information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
More informationCORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire
CORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire main ingredients: [LM (M3AS 00, JCP 00)] plane flow: D BGK Model conservative and entropic velocity discretization space discretization: finite
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 informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Discretization of Boundary Conditions Discretization of Boundary Conditions On
More informationFinite Difference Method for PDE. Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36
Finite Difference Method for PDE Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36 1 Classification of the Partial Differential Equations Consider a scalar second order partial
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationURANS Computations of Cavitating Flow around a 2-D Wedge by Compressible Two-Phase Flow Solver
URANS Computations of Cavitating Flow around a 2-D Wedge by Compressible Two-Phase Flow Solver *Yohan Choe 1), Hyeongjun Kim 1), Chongam Kim 2) 1), 2) Department of Aerospace Engineering, Seoul National
More informationChapter 4. Nonlinear Hyperbolic Problems
Chapter 4. Nonlinear Hyperbolic Problems 4.1. Introduction Reading: Durran sections 3.5-3.6. Mesinger and Arakawa (1976) Chapter 3 sections 6-7. Supplementary reading: Tannehill et al sections 4.4 and
More informationOutline. 1 Partial Differential Equations. 2 Numerical Methods for PDEs. 3 Sparse Linear Systems
Outline Partial Differential Equations Numerical Methods for PDEs Sparse Linear Systems 1 Partial Differential Equations 2 Numerical Methods for PDEs 3 Sparse Linear Systems Michael T Heath Scientific
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
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 informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
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 informationMimetic Methods and why they are better
Mimetic Methods and why they are better Blair Perot Theoretical and Computational Fluid Dynamics Laboratory February 25, 2015 Numerical Methods: Are more than Mathematics Math: Accuracy Stability Convergence
More informationSOE3213/4: CFD Lecture 1
What is CFD SOE3213/4: CFD Lecture 1 3d 3d Computational Fluid Dynamics { use of computers to study uid dynamics. Solve the Navier-Stokes Equations (NSE) : r:u = 0 Du Dt = rp + r 2 u + F 4 s for 4 unknowns,
More information2.2. Methods for Obtaining FD Expressions. There are several methods, and we will look at a few:
.. Methods for Obtaining FD Expressions There are several methods, and we will look at a few: ) Taylor series expansion the most common, but purely mathematical. ) Polynomial fitting or interpolation the
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
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 informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More information