FDM for wave equations
|
|
- Rhoda O’Neal’
- 5 years ago
- Views:
Transcription
1 FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D
2 Finite difference discretization
3 Finite difference discretization Discretization Explicit method
4 Finite difference discretization Local truncation error Consistency: yes!! Stability condition see details in class or as an exercise Convergence rate based on Lax equivalence theorem
5 FDM for wave equation u 2u + u Implicit scheme exercise n+ 1 n n 1 2 j j j ν θ n+ n+ n+ n n n n n n uj 1 uj uj 1 θ uj 1 uj uj 1 θ = uj+ 1 uj + uj 1 k h [ ( 2 ) (1 2 )( 2 ) ( 2 )] Efficient computation & error estimates Extension to 2D & 3D cases With variable coefficients & with other BCs
6 For nonlinear wave equation Nonlinear wave equation d u( xt, ) uxt (, ) + f( u) = 0, x, t> 0 tt ux (,0) = φ( x), ux (,0) = ϕ( x), x Sine-Gordon equation: Klein-Gordon equation: Soliton solution & energy conservation t f( u) = sin( u) f( u) = u+ gu ( ) Et = uxt + uxt + Fuxt dx E 2 2 () t (,) (,) ((,)) (0) d u 2 2 : = ϕ( ) + φ( ) + ( φ( )), ( ) = ( ) d x x F x dx F u f s ds 0 d
7 For nonlinear wave equation Numerical methods Semi-implicit scheme n+ 1 n n 1 uj 2uj + uj 1 n+ 1 n+ 1 n+ 1 n 1 n 1 n 1 n u 2 2 j+ 1 2uj + uj 1 + uj+ 1 2 uj + uj 1 + f( uj) = 0 k 2h Energy conservative scheme u 2 u + u Fu ( ) Fu ( ) = 0 n+ 1 n n 1 n+ 1 n 1 j j j 1 n+ 1 n+ 1 n+ 1 n 1 n 1 n 1 j j u 2 2 j 1 uj uj 1 uj 1 uj uj 1 n 1 n 1 k 2h uj uj Nonrelativistic limit regime -- current research topics!!! 2 1 ε tt ( xt, ) uxt (, ) + u+ f( u) = 0, 0< ε 1 2 ε
8 For nonlinear Schrodinger equation The nonlinear Schr o dinger equation (NLS) 1 i t (,) xt = + V( x) ψ ψ ψ β ψ ψ t : time & : spatial coordinate ψ (,) xt V( x ) : complex-valued wave function : real-valued external potential β : coupling constant =0: linear d x ( R ) >0: repulsive interaction <0: attractive interaction
9 Modelling In quantum physics & nonlinear optics: Interaction between particles with quantum effect Bose-Einstein condensation (BEC): bosons at low temperature Superfluidity Propagation of laser beams, In plasma physics: wave interaction between electrons and ions Zakharov system,.. In quantum chemistry: chemical interaction based on the first principle Schr o dinger-poisson system
10 Modelling In particle physics: interaction between fundamental particles Klein-Gordon-Schr o dinger (KGS) system Maxwell-Dirac (MD) system In biology: protein folding, Schr o dinger-poisson system,.. In materials science: First principle computation Semiconductor industry..
11 Numerical difficulties for dynamics Dispersive & nonlinear Solution and/or potential are smooth but may oscillate wildly Keep the properties of NLS on the discretized level Time reversible & time transverse invariant Mass & energy conservation Dispersion relation In high dimensions: many-body problems Design efficient & accurate numerical algorithms Explicit vs implicit (or computation cost) Spatial/temporal accuracy, Stability Resolution in strong interaction regime: β 1
12 Algorithms for computing dynamics Different methods Crank-Nicolson finite difference method (CNFD) Time-splitting spectral method (TSSP) Leap-frog (or RK4) + FD (or spectral) methods.. Time-splitting spectral method (TSSP) i tψ( xt, ) = ( A+ B) ψ with A=, B= V( x) + β ψ 2 ia t ib t 2 e e ψ ( xt, n) + O(( t) ) i( A+ B) t ia t/2 ib t ia t/2 3 ψ( xt, n+ 1) = e ψ( xt, n) e e e ψ( xt, n) + O(( t) ) 5 + O(( t) )
13 Algorithms for computing dynamics For, apply time-splitting technique [ tn, t n + 1] Step 1: Discretize by spectral method & integrate in phase space exactly 1 2 i tψ(,) xt = ψ 2 Step 2: solve the nonlinear ODE analytically 2 i tψ(,) xt = V( x) ψ(,) xt + β ψ(,) xt ψ(,) xt = = i ψ xt = V xψ xt + β ψ xt ψ xt t 2 t( ψ( xt, ) ) 0 ψ( xt, ) ψ( xt, n) 2 (, ) ( ) (, ) (, n) (, ) 2 i( t tn)[ V( x) + βψ ( xt, n) ] ψ( xt, ) = e ψ( xt, n) Use 2 nd order Strang splitting (or 4 th order time-splitting)
14 An algorithm for the dynamics of NLS ψ ψ ψ ( x) = e ψ ( x) ( x) = F e F( ψ ) ( x) = e ψ ( x) n 2 (1) i tv [ ( x) + βψ ( x) ]/2 n (2) 1 ia t (1) (2) 2 n+ 1 i tv [ ( x) + βψ ( x) ]/2 (2) Different spectral transforms for F: Fourier transform, Laguerre transform, Hermite transform, Gaussian transform, Different ways of transform (1) (2) ( d) For small d: ψ( x1, x2,..., xd) = ψ j1, j2,..., jφ ( 1 1) ( 2 2)... ( ) d j x φj x φj x d d N 1 j1, j2,..., jd N For large d: d ψ( x, x,..., x ) = ψ φ ( x ) φ ( x )... φ ( x ) N(log( N)) (1) (2) ( d) 1 2 d j1, j2,..., jd j1 1 j2 2 jd d j j... j N 1 2 d d
15 Properties of the method Explicit & computational cost per time step: Time reversible: yes Time transverse invariant: yes Mass conservation: yes Stability: yes Dispersion relation without potential: yes Accuracy Spatial: spectral order; Temporal: 2 nd or 4 th order OM ( ln M) Best resolution in strong interaction regime: β 1
16 For hyperbolic conservation laws Consider 1D conservation law Where f is a convex function Three types of problems Cauchy problem -- defined in the whole space Riemann problem initial data has discontinuity Initial-boundary problem
17 Special cases Linear case, i.e. f(u) = a u Exact solution Burgers equation It can generate shock solution for smooth initial data
18 Special cases Riemann problem Shock solution Shock speed satisfies the Rankine-Hugoniot condition Rarefaction solution
19 Properties Conservation Characteristics Line for linear case Curve for nonlinear problem Dependence region Wind direction Shock formation
20 Finite difference discretization Numerical difficulties Discontinuous solution Keep physical laws in the discretized level Boundary condition setups must according to characteristics.. Numerical methods Shock tracking random selection method Shock capturing Upwind method, Lax-Friedrichs method,. Total variation diminishing (TVD) schemes, Godunov scheme, ENO or WENO scheme Gas-kinetic scheme, relaxation scheme,
21 Finite difference discretization
22 A fake method Forward Euler + 2 nd central finite difference Stability von Neumann method Assume Plugging into the difference equation Amplification factor It is unconditionally unstable!!!
23 Numerical method Lax-Friedrichs method For linear problem, e.g. f(u) = 3u For Burgers equation Convergence Local truncation error -- consistency Conservation: yes!!
24 Numerical method CFL (Courant-Friedrichs-Lax, 1928) condition--- stability For linear problem For nonlinear case TVD property Total variation Lax-Friedrichs is a TVD scheme under the stability condition Convergence Lax equivalent theorem
25 High resolution scheme Godunov (Godunov, 1959) scheme Use the conservation of the equation & (approximate) Riemann solver
26 High resolution scheme Integrate exactly Conservative scheme
27 Second-order method
28 Second-order method Choice of slope limiter Roe s superbee limiter Van Leer limiter..
29 Examples
30 Numerical methods Properties of Godunov schemes Local truncation error: Consistency: yes!! Explicit, conservative & TVD method CFL (stability) condition: Convergent under the stability condition More topics High dimensions dimension-by-dimension System of conservation laws Euler system, shallow water,...
31 Finite element method FEM for Poisson equation Variational formulation Sobolev spaces: { x } j 1 1 { φ( ) φ ( ), φ Γ 0 }, g { φ( ) φ ( ), φ Γ } Variational problem ux ( ) = f( x) in Ω ux ( ) = gx ( ) on Γ = Ω (P1) L ( Ω ) = φ( x) φ( x) dx <, ( φϕ, ) = φ( x) ϕ( x) dx, φ = ( φφ, ) 2 2 1/2 2 L Ω Ω H x L L j d ( Ω ) = φ( ) φ ( Ω), φ ( Ω ), = 1,...,, φ = φ + φ, φ : = φ V = x H Ω = V = x H Ω = g H L L H L Find u V such that auv (, ) = f( v) v V g a( u, v) = u( x) v( x) dx, f ( v) = f ( x) v( x) dx (P2) Ω Ω 0
32 Finite element method Minimization formulation Theorem: Find u V such that (P3) 1 Eu ( ) = min Ev ( ) Ev ( ) = avv (, ) f( v) vv g 2 Partition of the domain In 1D: intervals g In 2D: triangles or rectangles or quadrilaterals, etc. In 3D: tetrahedrons, etc. If u C ( P1) ( P2) ( P3) ( P1) T : a = x < x < < x = b, T = [ x, x ], h = x x, h = max h h 0 1 N i i i+ 1 i i+ 1 i i 0 i N 1 2
33 Finite element method Finite element approximation Finite element subspaces { φ( ) ( ) φ ( ), 0,1,..., 1; φ( ) αφ, ( ) β} i { φ ( ) ( ) φ T ( i), 0,1,..., 1; φ ( ) φ ( ) 0} V = x C Ω P T i = M a = b = h h h k h h g T i V = x C Ω P T i = M a = b = h h h k h h 0 Finite element approximation Find u h Error Estimates (Cea lemma) h 0 i V h g such that au v = f v v V h h h h h h (, ) ( ) 0 (P ) u u C u v C u u Ch v = a v v V V h h h k 1/2 h h inf π, (, ), π : V h V V V g g v V
34 Finite element method How to form the linear system (for k=1 as example) Total stiff matrix & load vector h h h h h h h h T a( u, v ) = u ( x) v ( x) dx = u ( x) v ( x) dx : = u ( x) v ( x) dx = V AU Ω h h h h T f( v) = f( xv ) ( xdx ) = f( xv ) ( xdx ) : = f( xv ) ( xdx ) = V F Ω i h i Element stiff matrix & load vector i h T T T T T T T T T T i u a u v u x v x dx v v N x N x dx V AU xi+ 1 xi+ 1 h h h h T i T T(, ): = ( ) ( ) ( 1) ( ) ( ) : i x x = i i+ x i x i i i i u = x x i+ 1 i xi+ 1 xi+ 1 h h T xi+ 1 x x xi T( ): ( ) ( ) ( 1) ( ) ( ) ( ) i i i+ i i i i x x x i+ 1 xi xi+ 1 xi f v = f xv xdx= v v f xn xdx= V F N x = i i i i h x i+ 1 x i i h x i+ 1 x i T
35 Finite element method First compute element stiff matrix & load vector Assembling the total stiff matrix & load vector Plugging into the boundary condition How to solve the linear system: spare & large PCG or GMRES Multigrid,.. FEM software packages: ANSTRAN: ADIAN: More:
FDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationNumerical aspects of the nonlinear Schro dinger equation. Christophe Besse
Numerical aspects of the nonlinear Schro dinger equation by Christophe Besse Laboratoire Paul Painleve, Universite Lille 1, CNRS, Inria Simpaf Team Lille Outline of the talk 1 Motivation 2 Numerical methods
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
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 informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
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 informationLast time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler:
Lecture 7 18.086 Last time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler: U j,n+1 t U j,n = U j+1,n 2U j,n + U j 1,n x 2 Expected accuracy: O(Δt) in time,
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationQuantized Vortex Stability and Dynamics in Superfluidity and Superconductivity
Quantized Vortex Stability and Dynamics in Superfluidity and Superconductivity Weizhu Bao Department of Mathematics National University of Singapore Email: matbaowz@nus.edu.sg URL: http://www.math.nus.edu.sg/~bao
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
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 informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
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 informationNUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS JOHN A. TRANGENSTEIN Department of Mathematics, Duke University Durham, NC 27708-0320 Ш CAMBRIDGE ЩР UNIVERSITY PRESS Contents 1 Introduction
More informationComputer Physics Communications. Computational methods for the dynamics of the nonlinear Schrödinger/Gross-Pitaevskii equations
Computer Physics Communications (13) 1 3 Computer Physics Communications Computational methods for the dynamics of the nonlinear Schrödinger/Gross-Pitaevskii equations Xavier Antoine a,b, Weizhu Bao 1c,d,
More informationStrong Stability-Preserving (SSP) High-Order Time Discretization Methods
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods
More informationTyn Myint-U Lokenath Debnath. Linear Partial Differential Equations for Scientists and Engineers. Fourth Edition. Birkhauser Boston Basel Berlin
Tyn Myint-U Lokenath Debnath Linear Partial Differential Equations for Scientists and Engineers Fourth Edition Birkhauser Boston Basel Berlin Preface to the Fourth Edition Preface to the Third Edition
More informationCLASSROOM NOTES PART II: SPECIAL TOPICS. APM526, Spring 2018 Last update: Apr 11
CLASSROOM NOTES PART II: SPECIAL TOPICS APM526, Spring 2018 Last update: Apr 11 1 Function Space Methods General Setting: Projection into finite dimensional subspaces t u = F (u), u(t = 0) = u I, F : B
More informationNumerical Methods for Partial Differential Equations: an Overview.
Numerical Methods for Partial Differential Equations: an Overview math652_spring2009@colorstate PDEs are mathematical models of physical phenomena Heat conduction Wave motion PDEs are mathematical models
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 informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
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 informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationNumerical Methods for the Dynamics of the Nonlinear Schrodinger / Gross-Pitaevskii Equations. Weizhu Bao.
Numerical Methods for the Dynamics of the Nonlinear Schrodinger / Gross-Pitaevskii Equations Weizhu Bao Department of Mathematics National University of Singapore Email: matbaowz@nus.edu.sg URL: http://www.math.nus.edu.sg/~bao
More informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8,
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 informationMath Partial Differential Equations 1
Math 9 - Partial Differential Equations Homework 5 and Answers. The one-dimensional shallow water equations are h t + (hv) x, v t + ( v + h) x, or equivalently for classical solutions, h t + (hv) x, (hv)
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationSurprising Computations
.... Surprising Computations Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca www.cs.ubc.ca/ ascher/ Uri Ascher (UBC) Surprising Computations Fall 2012 1 / 67 Motivation.
More informationAttribution-NonCommercial-NoDerivs 2.0 KOREA. Share copy and redistribute the material in any medium or format
Attribution-NonCommercial-NoDerivs. KOREA You are free to : Share copy and redistribute the material in any medium or format Under the follwing terms : Attribution You must give appropriate credit, provide
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 informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationNumerical Methods for Hyperbolic Conservation Laws Lecture 4
Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Wen Shen Department of Mathematics, Penn State University Email: wxs7@psu.edu Oxford, Spring, 018 Lecture Notes online: http://personal.psu.edu/wxs7/notesnumcons/
More informationStability and instability of solitons in inhomogeneous media
Stability and instability of solitons in inhomogeneous media Yonatan Sivan, Tel Aviv University, Israel now at Purdue University, USA G. Fibich, Tel Aviv University, Israel M. Weinstein, Columbia University,
More informationA parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows.
A parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows Tao Xiong Jing-ei Qiu Zhengfu Xu 3 Abstract In Xu [] a class of
More informationNumerical Methods for Partial Differential Equations CAAM 452. Spring 2005
Numerical Methods for Partial Differential Equations Instructor: Tim Warburton Class Location: Duncan Hall 1046 Class Time: 9:5am to 10:40am Office Hours: 10:45am to noon in DH 301 CAAM 45 Spring 005 Homeworks
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 informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationPartial Differential Equations and the Finite Element Method
Partial Differential Equations and the Finite Element Method Pavel Solin The University of Texas at El Paso Academy of Sciences ofthe Czech Republic iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC, PUBLICATION
More informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
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 informationModeling, Analysis and Simulation for Degenerate Dipolar Quantum Gas
Modeling, Analysis and Simulation for Degenerate Dipolar Quantum Gas Weizhu Bao Department of Mathematics National University of Singapore Email: matbaowz@nus.edu.sg URL: http://www.math.nus.edu.sg/~bao
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 informationComputational Astrophysics 7 Hydrodynamics with source terms
Computational Astrophysics 7 Hydrodynamics with source terms Oscar Agertz Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term
More informationGas Dynamics Equations: Computation
Title: Name: Affil./Addr.: Gas Dynamics Equations: Computation Gui-Qiang G. Chen Mathematical Institute, University of Oxford 24 29 St Giles, Oxford, OX1 3LB, United Kingdom Homepage: http://people.maths.ox.ac.uk/chengq/
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationSolving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation
Solving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation W.L. Jin and H.M. Zhang August 3 Abstract: In this paper we study the Payne-Whitham (PW) model as
More informationA Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws
A Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws Mehdi Dehghan, Rooholah Jazlanian Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University
More informationYINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
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 informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationThe RAMSES code and related techniques 4. Source terms
The RAMSES code and related techniques 4. Source terms Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term - Relaxation towards
More informationFundamentals Physics
Fundamentals Physics And Differential Equations 1 Dynamics Dynamics of a material point Ideal case, but often sufficient Dynamics of a solid Including rotation, torques 2 Position, Velocity, Acceleration
More informationNumerical Solution of partial differential equations
G. D. SMITH Brunei University Numerical Solution of partial differential equations FINITE DIFFERENCE METHODS THIRD EDITION CLARENDON PRESS OXFORD Contents NOTATION 1. INTRODUCTION AND FINITE-DIFFERENCE
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 informationNon-linear Methods for Scalar Equations
Non-linear Methods for Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 3, 04 / 56 Abstract
More informationAdvanced numerical methods for transport and reaction in porous media. Peter Frolkovič University of Heidelberg
Advanced numerical methods for transport and reaction in porous media Peter Frolkovič University of Heidelberg Content R 3 T a software package for numerical simulation of radioactive contaminant transport
More informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
More informationMethod of Lines. Received April 20, 2009; accepted July 9, 2009
Method of Lines Samir Hamdi, William E. Schiesser and Graham W. Griffiths * Ecole Polytechnique, France; Lehigh University, USA; City University,UK. Received April 20, 2009; accepted July 9, 2009 The method
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
More informationSolutions to ODEs and PDEs: Numerical analysis using R (Cambridge University Press, ISBN: )
Solutions to ODEs and PDEs: Numerical analysis using R (Cambridge University Press, ISBN: 978-1-107-11561-3) Graham W. Griffiths City University, UK. Chapter Abstracts Chapter 1 ODE integration methods
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationA semi-lagrangian finite difference WENO scheme for scalar nonlinear conservation laws
A semi-lagrangian finite difference WENO scheme for scalar nonlinear conservation laws Chieh-Sen Huang a,, Todd Arbogast b,, Chen-Hui Hung c a Department of Applied Mathematics, National Sun Yat-sen University,
More informationLinear Partial Differential Equations for Scientists and Engineers
Tyn Myint-U Lokenath Debnath Linear Partial Differential Equations for Scientists and Engineers Fourth Edition Birkhäuser Boston Basel Berlin Tyn Myint-U 5 Sue Terrace Westport, CT 06880 USA Lokenath Debnath
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationFinite Differences: Consistency, Stability and Convergence
Finite Differences: Consistency, Stability and Convergence Varun Shankar March, 06 Introduction Now that we have tackled our first space-time PDE, we will take a quick detour from presenting new FD methods,
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationResearch Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation with Inconsistent Initial and Boundary Conditions
International Mathematics and Mathematical Sciences Volume 212, Article ID 82197, 13 pages doi:1.1155/212/82197 Research Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation
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 informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationThe Nonlinear Schrodinger Equation
Catherine Sulem Pierre-Louis Sulem The Nonlinear Schrodinger Equation Self-Focusing and Wave Collapse Springer Preface v I Basic Framework 1 1 The Physical Context 3 1.1 Weakly Nonlinear Dispersive Waves
More informationApplication of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms
Future Generation Computer Systems () 65 7 Application of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms R. Naidoo a,b, S. Baboolal b, a Department
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 informationVISCOUS FLUX LIMITERS
VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,
More informationarxiv: v1 [physics.comp-ph] 22 Feb 2013
Numerical Methods and Causality in Physics Muhammad Adeel Ajaib 1 University of Delaware, Newark, DE 19716, USA arxiv:1302.5601v1 [physics.comp-ph] 22 Feb 2013 Abstract We discuss physical implications
More informationNo-hair and uniqueness results for analogue black holes
No-hair and uniqueness results for analogue black holes LPT Orsay, France April 25, 2016 [FM, Renaud Parentani, and Robin Zegers, PRD93 065039] Outline Introduction 1 Introduction 2 3 Introduction Hawking
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationNumerical Schemes Applied to the Burgers and Buckley-Leverett Equations
University of Reading Numerical Schemes Applied to the Burgers and Buckley-Leverett Equations by September 4 Department of Mathematics Submitted to the Department of Mathematics, University of Reading,
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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationFinite volume method for conservation laws V
Finite volume method for conservation laws V Schemes satisfying entropy condition Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
More informationApproximations of diffusions. Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine
Lecture 3b Approximations of diffusions Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.1 04/10/2018 1 Learning objectives Become aware of the existence of stability conditions for the
More informationThe Discontinuous Galerkin Method for Hyperbolic Problems
Chapter 2 The Discontinuous Galerkin Method for Hyperbolic Problems In this chapter we shall specify the types of problems we consider, introduce most of our notation, and recall some theory on the DG
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem and Its Difference Approximations 1-D Initial Boundary Value
More informationNumerical methods Revised March 2001
Revised March 00 By R. W. Riddaway (revised by M. Hortal) Table of contents. Some introductory ideas. Introduction. Classification of PDE's.3 Existence and uniqueness.4 Discretization.5 Convergence, consistency
More informationPart 1. The diffusion equation
Differential Equations FMNN10 Graded Project #3 c G Söderlind 2016 2017 Published 2017-11-27. Instruction in computer lab 2017-11-30/2017-12-06/07. Project due date: Monday 2017-12-11 at 12:00:00. Goals.
More informationLecture I: numerical computation of the dynamics of Schrödinger equations. Extension to the GPE for BECs
/67 Lecture I: numerical computation of the dynamics of Schrödinger equations. Extension to the GPE for BECs Xavier Antoine Institut Elie Cartan de Lorraine (IECL) Université de Lorraine INRIA-CORIDA Team
More informationAMath 574 February 11, 2011
AMath 574 February 11, 2011 Today: Entropy conditions and functions Lax-Wendroff theorem Wednesday February 23: Nonlinear systems Reading: Chapter 13 R.J. LeVeque, University of Washington AMath 574, February
More information