Saul Abarbanel; Half a century of scientific work. Bertil Gustafsson, Uppsala University
|
|
- Elaine Stevens
- 5 years ago
- Views:
Transcription
1 Saul Abarbanel; Half a century of scientific work Bertil Gustafsson, Uppsala University
2 Grew up in Tel Aviv Served in Israeli Army during the War of Independence
3 MIT Ph.D 1959, Theoretical Aerodynamics
4 Post Doc Weizmann Insitute,
5 Tel Aviv University, Professor Head of Appl. Math. Dept., 1964 (As Associate Professor) Dean of Science Vice Rector, Rector Chairman National Research Council Director Sackler Institute of Advanced Studies
6 Visitor ICASE (NASA Langley)
7 Brown University Visitor IBM Distinguished Visiting Research Professor
8 Heat transfer, gas dynamics Most part mathematical analysis, little numerics. Abarbanel: J. Math. and Physics (1960) Time Dependent Temperature Distribution in Radiating Solids. Abarbanel: Israel Journal of Technology (1966) The deflection of confining walls by explosive loads. Abarbanel Zwas: J. Math. Anal. & Appl. (1969) The Motion of Shock Waves and Products of Detonation Confined between a Wall and a Rigid Piston. "...a detailed analytical solution of the piston motion and flow field is carried out..."
9 1969 Construction and analysis of difference methods for PDE Stability of PDE and difference methods Lax Wendroff type methods Compact high-order finite-difference schemes. Method of lines, Runge Kutta methods PML methods
10 Law Wendroff type methods and shocks u t = f(u) x von Neumann Richtmyer (1950): Add viscosity for numerical computation u t = f(u) x + ε 2 u x 2 Difference approximation "may be used for the entire calculation, just as though there were no shocks at all". 1954: Lax defines shocks as viscous limits ε 0 Dissipative difference methods for computation 1960: Lax Wendroff scheme, damping all frequencies 1969: MacCormack scheme, two stage, easier to apply Godunov methods (Riemann solvers), upwind methods, shock fitting
11 Lax-W methods: Possible oscillations near shock 97 il t6t t77 r95
12 Abarbanel Zwas: Math. Comp. (1969): An iterative finite-difference method for hyperbolic systems. Lax Wendroff type methods How to avoid oscillations near shocks? W t + F(W) x = 0 W t + A(W)W x = 0 Lax-W = W n j λ(f n 2 j+1 Fj 1) n + λ2 2 [An (F n j+1/2 j+1 F n j ) A n (F n j 1/2 j W n+1 j F n j 1)]
13 W n+1 = W n + Q W n Modify to W n+1 = W n + Q [θw n+1 + (1 θ)w n ] with iteration W n+1,s+1 = W n +Q [θw n+1,s +(1 θ)w n ], s = 0, 1,..., k 1, W n+1,0 = Analysis for different θ and different k: Courant number λ = t/ x No oscillations for 1 and 2 iterations
14 97 il t6t t77 r95
15 Abarbanel-Goldberg: J. Comp. Phys. (1972) Numerical Solution of Quasi-Conservative Hyperbolic Systems; The Cylindrical Shock Problem. General difference scheme Implicit scheme External: Internal: W t + [F(W)] x = Ψ(x; W) W n+1 = W n + CW n (1) W n+1,s+1 = W n + CW n + θ[cw n+1,s CW n ] W n+1,s+1 = W n + C(1 θ)w n + θcw n+1,s Iterative solver as in Abarbanel Zwas (1969), fixed number of iterations Larger timestep compared to explicit solver.
16 Standard scheme i nt,i iexocl) t1 (opprox.) , ? l.o 7?
17 Internal scheme
18 Use of time-dependent methods for computation of steady state. Abarbanel-Dwoyer-Gottlieb: J. Comp. Phys. (1986) Improving the Convergence Rate to Steady State of Parabolic ADI Methods. u t = u xx + u yy ADI-methods: Peaceman Rachford (1955)... Beam Warming (1976) (1 λδ 2 x)(1 λδ 2 y)(v n+1 v n ) = αλ(δ 2 x + δ 2 y)v n, λ = t/h 2 Improve convergence rate as n by adding extra term (1 λδ 2 x)(1 λδ 2 y)(v n+1 v n ) = αλ(δ 2 x+δ 2 y)v n + γ 4 λ2 δ 2 xδ 2 y(δ 2 x + δ 2 y)v n Fourier analysis. Choose γ to minimize amplification factor. Model equation γ = 0.8 independent of mesh-size.
19 Compact Pade type difference methods Orzag 1971, Kreiss-Oliger 1972: pseudospectral methods high order accuracy. Number of points per wavelength? High order difference methods? Pade (1890): Approximation of functions by rational functions Lele 1992: "Compact Finite Difference Schemes with Spectral-like Resolution" v = u/ x v j+1 + 4v j + v j 1 = 1 h (3u j+1 3u j 1 ) (4th order)
20 Approximation ˆQ(ξ) of ξ in Fourier space 0 ξ π Standard 4th order, standard 6th order, compact 4th order
21 Boundary conditions? Stability? Lele: Numerical computation of eigenvalues of difference operators, fixed x.
22 Carpenter-Gottlieb-Abarbanel, J. Comp. Phys. (1993) The stability of numerical boundary treatments for compact high-order finite-difference schemes. Normal mode stability analysis (GKS). "Weak point: complexity in its application to higher order numerical schemes." Extra consideration: Fixed t: Growing solutions V(t) Ce αt V(0)? Time-stable if α = 0. Analysis and construction of boundary conditions leading to time stability. Extensive thorough analysis, but for scalar case.
23 . SBP-operators (Summation By Parts) Kreiss Scherer (1977) u t = u x, 0 x 1, u(1, t) = g(t), u(x, 0) = f(x) (v, x v) = 1 2 ( v(1) 2 v(0) 2 ) for all v d dt u 2 = u(1, t) 2 u(0, t) 2 SBP: Construct scalar product (u, v) h and a difference operator D such that (v, Dv) h = 1 2 ( v N 2 v 0 2 )
24 Simultaneous Approximation Terms (SAT) Funaro 1988, Funaro Gottlieb 1988: SAT for pseudospectral methods Add penalty term dv dt = Dv τ ( v N g(t) ) w (2) Carpenter-Gottlieb-Abarbanel, J. Comp.Phys. (1994) Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: Methodology and application to high-order compact schemes. Previous article (1993) with stable and time-stable methods are constructed for the scalar case. Use SAT method based on SBP-operators for systems This article: A systematic way of constructing time-stable SAT.
25 Abarbanel Ditkowski, J. Comp. Phys. (1997) Asymptotically Stable Fourth-Order Accurate Schemes for the Diffusion Equation on Complex Shapes 4-th order, nonsymmetric difference operators near boundaries, "SAT-type". Solution bounded by constant independent of t.
26 Method of lines Carpenter-Gottlieb-Abarbanel-Don: SIAM J. Sci. Comput. (1995) The theoretical accuracy of Runge Kutta time discretizations for the initial boundary value problem: A study of the boundary error. u t + u t = 0, 0 x 1, u(0, t) = g(t) Physical boundary condition at each stage of the R-K method (4th order) v 1 0 = g(t + δt Theoretical analysis showing deterioration of accuracy. Use instead derivative boundary conditions derived from original b.c. v ) = g(t) + δt 2 g (t) Full accuracy for the linear case, only 3rd order in nonlinear case.
27 Abarbanel Gottlieb, J. Comp. Phys. (1981): Optimal Time Splitting for Two- and Three-Dimensional Navier-Stokes Equations with Mixed Derivatives (33 pages) Interview by Philip Davis 2003: "Perhaps the most important article" U = [ρ, ρu, ρv, ρw, e] T U t + F x + G y + H z = 0 V = [ρ, u, v, w, p] T V t +AV x +BV y +JV z = CV xx +DV yy +K V zz +E xy V xy +E yz V yz +E zx V xz Similarity transformation such that S 1 MS are symmetric for all matrixes M = A, B,..., E zx
28 U t + (F H + F P + F M ) x + (G H + G P + G M ) y + (H H + H P + H M ) z = 0 U n+2 = [L x ( t x )L y ( t y )L z ( t z )L xyz ( t xyz )L xx ( t xx )L yy ( t yy )L zz ( t zz )] [L zz ( t zz )L yy ( t yy )L xx ( t xx )L xyz ( t xyz )L z ( t z )L y ( t y )L x ( t x )]U n L x..., L xx... MacCormack solvers L xyz MacCormack-like solver
29 Scalar equation: u t = au x + bu y + ju z + cu xx + du yy + ku zz + e xy u xy + e yz u yz + e zx u zx Stability under the standard one-dimensional conditions and t xyz t x. a t x x 1,... c t xx 1,... ( x) 2 2 The same stability result for the Navier-Stokes equations due to symmetric coefficient matrices.
30 Abarbanel-Duth-Gottlieb: Computers & Fluids (1989) Splitting methods for low Mach number Euler and Navier-Stokes equations Stiff system Splitting Symmetrizing Stiffness isolated to linear system ("may be solved implicitly with ease")
31 Abarbanel-Chertock: J. Comp. Phys. (2000) Strict Stability of High-Order Compact Implicit Finite-Difference Schemes: The Role of Boundary Conditions for Hyperbolic PDEs, I,II Derivation of general compact implicit methods.
32 Absorbing boundary conditions Enquist Majda (1977): Wave equation u tt = u xx + u yy, < x, y < Boundary conditions for finite domain x x 0? Fourier transform ω 2 = ξ 2 + η 2 ξ = ±ω 1 η 2 /ω 2, +ω for leftgoing wave Pseudo-differential equation. η/ω small 1 η 2 /ω 2 1 η2 2ω 2 ξω ω η2 = 0 boundary condition at x = x 0 2 u x t = t 2 2 y 0 2
33 Berenger (1994): (Centre d Analyse de Dèfense, France) Perfectly Matched Layers (PML). Outer boundaries of computational domain Absorbing layer y x
34 Maxwell equations 2D W = [E x, E y, H z ] T W t = A W x + B W y + CW Can be symmetrized. PML formulation W b = [E x, E y, H zx, H zy ] T W b t = A b W b x + B W b b y + C bw b
35 Abarbanel-Gottlieb, J. Comp. Phys. (1997) A mathematical analysis of the PML method New system cannot be symmetrized. Shown in the article: Initial value problem weakly well posed: Fourier transform / x iω 1 / y iω 2 Explicit form of transformed system is derived. Ĥ x (t) (αω 1 + βω 2 )t Requires bounded derivatives, but still growth in time.
36 Even worse: Perturbation 0 0 δ δ 0 0 δ δ Compute eigenvalues λ Ill posed! λ 1 ωδ Ŵ(t) e ωδt Similar results for semi-discrete and fully discrete approximations.
37 Abarbanel-Gottlieb, Appl. Numer. Math., 1998 On the construction and analysis of absorbing layers in CEM. New PML type formulation. Introduce new variable polarization current J (Zilkowski 1997) E x t J t = Hz y = σ Hz y J P = J + σe x P t = σp + σ 2 E x Strongly well posed (even when the outer boundary is taken into account). Still another formulation constructed, strongly well posed.
38 Abarbanel-Gottlieb-Hesthaven, J. Comp. Phys., 1999 Well-posed Perfectly Matched Layers for Advective Acoustics Development based on Abarbanel-Gottlieb (1998). "...somewhat lengthy algebraic manipulations..." Strongly well posed Numerical method: 4th order in space, Runge Kutta in time
39 Abarbanel-Gottlieb-Hesthaven, J. Sci. Comp Long Time Behavior of the Perfectly Matched Layer Equations in Computational Electromagnetics PML-method of Abarbanel Gottlieb (1998) shows long time growth (after the initial pulse has left the original domain).
40 0 t 70
41 a l a X "): 0 t 5000
42 Analysis of source of the problem Double eigenvalue, one eigenvector Cure: Split the eigenvalues by introducing small perturbation ε Uncertainty about damping properties in the PML-layer
43 Abarbanel-Quasimov-Tsynkov: J. Sci. Comp. (2009) Long-Time Performance of Unsplit PMLs with Explicit Second Order Schemes. Long-time growth with PML analyzed. Sensitive to choice of numerical method. Perturbation may or may not enter the original domain from PML-layer. "Lacunae based stabilization" by Qasimov-Tsynkov (2008).
44 Last publication: Abarbanel-Ditkowski: Appl. Numer.Math. (2015) Wave propagation in advected acoustics within a non-uniform medium under the effect of gravity. Saul 84 years old.
45
arxiv: v1 [math.na] 21 Nov 2017
High Order Finite Difference Schemes for the Heat Equation Whose Convergence Rates are Higher Than Their Truncation Errors, A. Ditkowski arxiv:7.0796v [math.na] Nov 07 Abstract Typically when a semi-discrete
More informationNUMERICAL SOLUTION OF THE LINEARIZED EULER EQUATIONS USING HIGH ORDER FINITE DIFFERENCE OPERATORS WITH THE SUMMATION BY PARTS PROPERTY
NUMERICAL SOLUTION OF THE LINEARIZED EULER EQUATIONS USING HIGH ORDER FINITE DIFFERENCE OPERATORS WITH THE SUMMATION BY PARTS PROPERTY Stefan Johansson Department of Information Technology Scientific Computing
More informationWell-posedness, stability and conservation for a discontinuous interface problem: an initial investigation.
Well-posedness, stability and conservation for a discontinuous interface problem: an initial investigation. Cristina La Cognata and Jan Nordström Abstract A robust interface treatment for the discontinuous
More informationEfficient wave propagation on complex domains
Center for Turbulence Research Annual Research Briefs 2006 223 Efficient wave propagation on complex domains By K. Mattsson, F. Ham AND G. Iaccarino 1. Motivation and objectives In many applications, such
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 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 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 informationStrict Stability of High-Order Compact Implicit Finite-Difference Schemes: The Role of Boundary Conditions for Hyperbolic PDEs, I
Journal of Computational Physics 160, 42 66 (2000) doi:10.1006/jcph.2000.6420, available online at http://www.idealibrary.com on Strict Stability of High-Order Compact Implicit Finite-Difference Schemes:
More informationStable and high-order accurate finite difference schemes on singular grids
Center for Turbulence Research Annual Research Briefs 006 197 Stable and high-order accurate finite difference schemes on singular grids By M. Svärd AND E. van der Weide 1. Motivation and objectives The
More informationDiagonal-norm upwind SBP operators
Diagonal-norm upwind SBP operators Ken Mattsson June 8, 16 Abstract High-order accurate first derivative finite difference operators are derived that naturally introduce artificial dissipation. The boundary
More informationSpectral analysis of the incompressible Navier-Stokes equations with different boundary conditions
Spectral analysis of the incompressible Navier-Stokes equations with different boundary conditions Cristina La Cognata, Jan Nordström Department of Mathematics, Computational Mathematics, Linköping University,
More informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
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 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 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 informationFDM 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 Solution of Initial Boundary Value Problems. Jan Nordström Division of Computational Mathematics Department of Mathematics
Numerical Solution of Initial Boundary Value Problems Jan Nordström Division of Computational Mathematics Department of Mathematics Overview Material: Notes + GUS + GKO + HANDOUTS Schedule: 5 lectures
More informationOptimal diagonal-norm SBP operators
Optimal diagonal-norm SBP operators Ken Mattsson 1, Martin Almquist 1 and Mark H. Carpenter 2 1 Department of Information Technology, Uppsala University 2 Computational Aerosciences Branch, NASA Langley
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 informationSome notes about PDEs. -Bill Green Nov. 2015
Some notes about PDEs -Bill Green Nov. 2015 Partial differential equations (PDEs) are all BVPs, with the same issues about specifying boundary conditions etc. Because they are multi-dimensional, they can
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 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 information1 PART1: Bratu problem
D9: Advanced Numerical Analysis: 3 Computer Assignment 3: Recursive Projection Method(RPM) Joakim Möller PART: Bratu problem. Background The Bratu problem in D is given by xx u + yy u + λe u =, u Γ = ()
More information1. Numerical Analysis of Spectral Methods/Theory and Applications, with S. Orszag, CBMS-SIAM No. 26, 1977, 170 pages.
DAVID GOTTLIEB Ford Foundation Professor Professor of Applied Mathematics Brown University Publications Updated February 2007 Books 1. Numerical Analysis of Spectral Methods/Theory and Applications, with
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
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 informationHigh Order Difference Approximations for the Linearized Euler Equations
High Order Difference Approximations for the Linearized Euler Equations Stefan Johansson Abstract The computers of today make it possible to do direct simulation of aeroacoustics, which is very computationally
More informationAE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationOn Error Bounds of Finite Difference Approximations to Partial Differential EquationsTemporal Behavior and Rate of Convergence
Journal of Scientific Computing, Vol. 15, No. 1, 2000 On Error Bounds of Finite Difference Approximations to Partial Differential EquationsTemporal Behavior and Rate of Convergence Saul Abarbanel, 1 Adi
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 informationCONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION
Journal of Computational Acoustics, Vol. 8, No. 1 (2) 139 156 c IMACS CONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION MURTHY N. GUDDATI Department of Civil Engineering, North Carolina
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 informationWell-posedness, Stability and Conservation for a Discontinuous Interface Problem
Department of Mathematics Well-posedness, Stability and Conservation for a Discontinuous Interface Problem Cristina La Cognata and Jan Nordström LiTH-MAT-R--214/16--SE Department of Mathematics Linköping
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 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 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 informationComputational Domain
Long Time Behavior of the Perfectly Matched Layer Equations in Computational Electromagnetics S. Abarbanel y, D. Gottlieb z, and J. S. Hesthaven z 1 y Department of Applied Mathematics, Tel Aviv University,
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 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 in a 2D Box Region Let us consider a model problem of parabolic
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 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 informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
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 informationA Stable and Accurate Davies-like Relaxation Procedure using Multiple Penalty Terms for Lateral Boundary Conditions
A Stable and Accurate Davies-like Relaxation Procedure using Multiple Penalty Terms for Lateral Boundary Conditions Hannes Frenander Division of Computational Mathematics, Department of Mathematics, Linköping
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationApplication of Computational Fluid Dynamics (CFD) Based Technology to Computational Electromagnetics Ramesh K. Agarwal
Application of Computational Fluid Dynamics (CFD) Based Technology to Computational Electromagnetics Ramesh K. Agarwal IEEE Distinguished Lecturer The William Palm Professor of Engineering Washington University
More informationCharacterizing the Accuracy of Summation-by-Parts Operators for Second-Derivatives with Variable-Coefficients
1 Characterizing the Accuracy of Summation-by-Parts Operators for Second-Derivatives with Variable-Coefficients by Guang Wei Yu Supervisor: D. W. Zingg April 13, 2013 Abstract This paper presents the
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 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 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 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 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 informationA Study of Transonic Flow and Airfoils. Presented by: Huiliang Lui 30 th April 2007
A Study of Transonic Flow and Airfoils Presented by: Huiliang Lui 3 th April 7 Contents Background Aims Theory Conservation Laws Irrotational Flow Self-Similarity Characteristics Numerical Modeling Conclusion
More informationThe PML Method: Continuous and Semidiscrete Waves.
Intro Continuous Model. Finite difference. Remedies. The PML Method: Continuous and Semidiscrete Waves. 1 Enrique Zuazua 2 1 Laboratoire de Mathématiques de Versailles. 2 Universidad Autónoma, Madrid.
More informationBoundary Conditions for a Divergence Free Velocity-Pressure Formulation of the Incompressible Navier-Stokes Equations
Boundary Conditions for a Divergence Free Velocity-Pressure Formulation of the Incompressible Navier-Stokes Equations Jan Nordström, Ken Mattsson and Charles Swanson May 5, 6 Abstract New sets of boundary
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 informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
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 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 informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
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 informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationConjugate Heat Transfer for the Unsteady Compressible Navier-Stokes Equations Using a Multi-block Coupling
Conjugate Heat Transfer for the Unsteady Compressible Navier-Stokes Equations Using a Multi-block Coupling Jan Nordström Department of Mathematics, Linköping University, SE-58 83 Linköping, Sweden Jens
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
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 informationKrylov single-step implicit integration factor WENO methods for advection-diffusion-reaction equations
Accepted Manuscript Krylov single-step implicit integration factor WENO methods for advection diffusion reaction equations Tian Jiang, Yong-Tao Zhang PII: S0021-9991(16)00029-2 DOI: http://dx.doi.org/10.1016/j.jcp.2016.01.021
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 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 informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
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 informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
More informationTwo-Dimensional Riemann Solver for Euler Equations of Gas Dynamics
Journal of Computational Physics 167, 177 195 (2001) doi:10.1006/jcph.2000.6666, available online at http://www.idealibrary.com on Two-Dimensional Riemann Solver for Euler Equations of Gas Dynamics M.
More informationApplication of Nodal Discontinuous Glaerkin Methods in Acoustic Wave Modeling
1 Application of Nodal Discontinuous Glaerkin Methods in Acoustic Wave Modeling Xin Wang ABSTRACT This work will explore the discontinuous Galerkin finite element method (DG-FEM) for solving acoustic wave
More informationX i t react. ~min i max i. R ij smallest. X j. Physical processes by characteristic timescale. largest. t diff ~ L2 D. t sound. ~ L a. t flow.
Physical processes by characteristic timescale Diffusive timescale t diff ~ L2 D largest Sound crossing timescale t sound ~ L a Flow timescale t flow ~ L u Free fall timescale Cooling timescale Reaction
More informationICASE ON THE REMOVAL OF BOUNDARY ERRORS CAUSED BY RUNGE-KUTTA INTEGRATION OF NON-LINEAR PARTIAL DIFFERENTIAL EQUATIONS
DTIC pi FCTF NASA Contractor Report 194989 %-- Q^Q g7 1994 1 ICASE Report No. 94-79 ICASE ON THE REMOVAL OF BOUNDARY ERRORS CAUSED BY RUNGE-KUTTA INTEGRATION OF NON-LINEAR PARTIAL DIFFERENTIAL EQUATIONS
More informationRunge Kutta Chebyshev methods for parabolic problems
Runge Kutta Chebyshev methods for parabolic problems Xueyu Zhu Division of Appied Mathematics, Brown University December 2, 2009 Xueyu Zhu 1/18 Outline Introdution Consistency condition Stability Properties
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
More informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationThird and Fourth Order Accurate Schemes for Hyperbolic Equations of Conservation Law Form*
MATHEMATICS OF COMPUTATION, VOLUME 25, NUMBER 114, APRIL, 1971 Third and Fourth Order Accurate Schemes for Hyperbolic Equations of Conservation Law Form* By Gideon Zwas and Saul Abarbanel Abstract. It
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationModule 2: Governing Equations and Hypersonic Relations
Module 2: Governing Equations and Hypersonic Relations Lecture -2: Mass Conservation Equation 2.1 The Differential Equation for mass conservation: Let consider an infinitely small elemental control volume
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 informationLASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE
LASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE H. M. Al-Qahtani and S. K. Datta University of Colorado Boulder CO 839-7 ABSTRACT. An analysis of the propagation of thermoelastic waves
More information1 Upwind scheme for advection equation with variable. 2 Modified equations: numerical dissipation and dispersion
1 Upwind sceme for advection equation wit variable coefficient Consider te equation u t + a(x)u x Applying te upwind sceme, we ave u n 1 = a (un u n 1), a 0 u n 1 = a (un +1 u n ) a < 0. CFL condition
More informationStable boundary treatment for the wave equation on second-order form
Stable boundary treatment for the wave equation on second-order form Ken Mattsson Frank Ham Gianluca Iaccarino June 24, 2008 Abstract A stable and accurate boundary treatment is derived for the secondorder
More informationHIGH ORDER LAX-WENDROFF-TYPE SCHEMES FOR LINEAR WAVE PROPAGATION
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 HIGH ORDER LAX-WENDROFF-TYPE SCHEMES FOR LINEAR WAVE PROPAGATION
More informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
More informationNew Diagonal-Norm Summation-by-Parts Operators for the First Derivative with Increased Order of Accuracy
AIAA Aviation -6 June 5, Dallas, TX nd AIAA Computational Fluid Dynamics Conference AIAA 5-94 New Diagonal-Norm Summation-by-Parts Operators for the First Derivative with Increased Order of Accuracy David
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 informationarxiv: v1 [physics.comp-ph] 30 Sep 2015
On the quasi-unconditional stability of BDF-ADI solvers for the compressible Navier-Stokes equations arxiv:1509.09213v1 [physics.comp-ph] 30 Sep 2015 Oscar P. Bruno and Max Cubillos Abstract The companion
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
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 informationLecture 17: Initial value problems
Lecture 17: Initial value problems Let s start with initial value problems, and consider numerical solution to the simplest PDE we can think of u/ t + c u/ x = 0 (with u a scalar) for which the solution
More informationImplicit kinetic relaxation schemes. Application to the plasma physic
Implicit kinetic relaxation schemes. Application to the plasma physic D. Coulette 5, E. Franck 12, P. Helluy 12, C. Courtes 2, L. Navoret 2, L. Mendoza 2, F. Drui 2 ABPDE II, Lille, August 2018 1 Inria
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 information