Numerical Methods for Wave Equations: Part I: Smooth Solutions

Total Page:16

File Type:pdf, Size:1020Kb

Numerical Methods for Wave Equations: Part I: Smooth Solutions WPPII Computational Fluid Dynamics I Numerical Methods for Wave Equations: Part I: Smooth Solutions Instructor: Hong G. Im University of Michigan Fall 2001 Computational Fluid Dynamics I WPPII Outline Solution Methods for Wave Equation Part I • Method of Characteristics • Finite Volume Approach and Conservative Forms • Methods for Continuous Solutions - Central and Upwind Difference - Stability, CFL Condition - Various Stable Methods Part II • Methods for Discontinuous Solutions - Burgers Equation and Shock Formation - Entropy Condition - Various Numerical Schemes WPPII Computational Fluid Dynamics I Method of Characteristics WPPII Computational Fluid Dynamics I 1st Order Wave Equation ∂f ∂f + c = 0 ∂t ∂x The characteristics for this equation are: dx df = c; = 0; f dt dt t f x WPPII Computational Fluid Dynamics I 1-D Wave Equation (2nd Order Hyperbolic PDE) ∂ 2 f ∂ 2 f − c2 = 0 ∂t 2 ∂x2 Define ∂f ∂f v = ; w = ; ∂t ∂x which leads to ∂v ∂w − c2 = 0 ∂t ∂x ∂w ∂v − = 0 ∂t ∂x WPPII Computational Fluid Dynamics I In matrix form v 0 − c2 v t + x = 0 or u + Au = 0 − t x wt 1 0 wx Can it be transformed into the form v λ 0v t + x = 0 or u + λu = 0 ? λ t x wt 0 wx Find the eigenvalue, eigenvector ()AT − λI q = 0 − λ −1 l 1 = 0 − 2 − λ c l2 WPPII Computational Fluid Dynamics I Eigenvalue dx AT − λI = 0; λ2 − c2 = 0; λ = = ±c dt Eigenvector 1 For λ = +c, − cl − l = 0; q = 1 2 1 − c 1 For λ = −c, cl − l = 0; q = 1 2 2 c The solution (v,w) is governed by ODE’s along the characteristic lines dx / dt = ±c WPPII Computational Fluid Dynamics I On dx / dt = +c, T ⋅()+ = ()− 2 − ()− q1 ut Au x vt c wx c wt vx dv dw = ()()v + cv − c w + cw = − c = 0 t x t x dt dt ≡ − If r1 v cw (Riemann Invariant), we have d(v − cw) = 0 on dx = +cdt (If c = const) dt WPPII Computational Fluid Dynamics I On dx / dt = −c, T ⋅()+ = ()− 2 + ()− q2 ut Au x vt c wx c wt vx dv dw = ()()v − cv + c w − cw = + c = 0 t x t x dt dt ≡ + If r1 v cw (Riemann Invariant), we have d(v + cw) = 0 on dx = −cdt (If c = const) dt WPPII Computational Fluid Dynamics I 3 dx / dt = −c P 2 dx / dt = +c dv dw − c = 0 on dx = +cdt dt dt 1 dv + dw = = − ∂f ∂f c 0 on dx cdt f , , known dt dt ∂t ∂x = + df ft dt f x dx WPPII Computational Fluid Dynamics I Finite difference approximation to determine f (P) dx / dt = +c []= = v ft , w f x 2 P dx / dt = −c 1 − = [][]+ − x(P) x(1) c(P) c(1) t(P) t(1) dx 2 on = +c 1 1 dt v(P) − v(1) = [][]c(P) + c(1) w(P) − w(1) 2 1 − = − [][]+ − x(P) x(2) c(P) c(2) t(P) t(2) dx 2 on = −c 1 dt v(P) − v(2) = − [][]c(P) + c(2) w(P) − w(2) 2 1 1 f (P) − f (1) = [][]v(P) + v(1) t(P) − t(1) + [][]w(P) + w(1) x(P) − x(1) 2 2 Nonlinear system – iterative procedure WPPII Computational Fluid Dynamics I Numerical Methods: Finite Volume Approach WPPII Computational Fluid Dynamics I When using finite volume approximations, we work directly with the integral form of the conservation principles. The average values of f over a small volume are stored f f j−1 f j f j+1 x xj−1/2 xj+1/2 WPPII Computational Fluid Dynamics I In finite volume method, equations in conservative forms are needed in order to satisfy conservation properties. As an example, consider a 1-D equation ∂f (x,t) ∂F[ f (x,t)] + = 0 ∂t ∂x where F denotes a general advection/diffusion term, e.g. 1 ∂f F = f 2 ; F = µ(x) 2 ∂x WPPII Computational Fluid Dynamics I Integrating over the domain L, ∂f ∂F dx + dx = 0 ∫ ∂ ∫ ∂ L t L x F(L) − F(0) = 0 d ⇒ ∫ fdx = 0 dt L If F = 0 at the end points of the domain, f is conserved. WPPII Computational Fluid Dynamics I f j−1 In discretized form: f j f j+1 − ∂F F + F − ∫ dx = ∑ j 1/ 2 j 1/ 2 ∂x ∆x L xj−1/2 xj+1/2 x 1 = []+ F − − F − + F + − F − + F + − F + + ∆x ! j 1/ 2 j 3/ 2 j 1/ 2 j 1/ 2 j 3/ 2 j 1/ 2 ! 1 = []F − F ∆x L 0 WPPII Computational Fluid Dynamics I Examples of Conservative Form ∂f ∂ 1 ∂f ∂ ∂f + f 2 = 0, + µ(x) = 0 ∂t ∂x 2 ∂t ∂x ∂x ∂ ∂f Discretize µ(x) ∂x ∂x 1 ∂f ∂f µ − µ = ∂ ∂ h x j+1/ 2 x j−1/ 2 1 []()()()()µ + µ f − f − µ + µ f − f 2h2 j+1 j j+1 j j j−1 j j−1 Conservative WPPII Computational Fluid Dynamics I Examples of Non-conservative Form ∂f ∂f ∂f ∂ 2 f ∂µ ∂f + f = 0, + µ + = 0 ∂t ∂x ∂t ∂x2 ∂x ∂x ∂ 2 ∂µ ∂ Discretize µ f + f ∂x2 ∂x ∂x 1 µ ()+ − + 1 ()()µ − µ − = j f j+1 f j−1 2 f j j+1 j−1 f j+1 f j−1 h2 4 1 ()µ + µ − µ + 1 ()µ − µ − µ + µ j f j+1 j f j−1 2 j f j j+1 f j+1 j+1 f j−1 j−1 f j+1 j−1 f j−1 h2 4 Non-conservative WPPII Computational Fluid Dynamics I Finite Volume Method for Conservative Equations ∂f ∂F + = 0 ∂t ∂x F = Uf Advection ∂f F = D Diffusion ∂x ∂f F = Uf − D Advection/Diffusion ∂x WPPII Computational Fluid Dynamics I Finite Volume Formulation f j−1 f j ∂f ∂F f + + = 0 j 1 ∂t ∂x ∂f ∂F j−1 j j+1 x ∫ dx + ∫ dx = 0 ∂ ∂ Fj−1/2 Fj+1/2 ∆x t ∆x x d + − = ∫ f dx Fj+1/ 2 Fj−1/ 2 0 dt ∆x d (hf ) + F − F = 0 dt j j+1/ 2 j−1/ 2 WPPII Computational Fluid Dynamics I 1-D Advection-Diffusion Equation ∂f ∂ ∂f + Uf − D = 0 ∂t ∂x ∂x FVM Equation: d ∂f ∂f h f = −U ( f − f ) + D − j j+1/ 2 j−1/ 2 ∂ ∂ dt x j+1/ 2 x j−1/ 2 Approximating f ≈ 1 ( f n + f n ) j+1/ 2 2 j+1 j ∂ f n − f n f ≈ j+1 j ∂ x j+1/ 2 h d 1 + f ≈ ( f n 1 − f n ) dt j ∆t j j WPPII Computational Fluid Dynamics I Substituting: h + ( f n 1 − f n ) = ∆t j j f n − f n f n − f n − []1 n + n − 1 n + n + j+1 j − j j−1 U ( f j+1 f j ) ( f j f j−1 ) D 2 2 h h Rearranging the terms: f n+1 − f n j j U n n D ()n n n + ( f + − f − ) = f + − 2 f + f − ∆t 2h j 1 j 1 h2 j 1 j j 1 Which is exactly the same as the finite difference equation if we take the average value to be the same as the value in the center of the cell WPPII Computational Fluid Dynamics I Numerical Methods for 1-D Advection Equation: Stability Consideration (Finite Difference Approach) WPPII Computational Fluid Dynamics I We will start by examining the linear advection equation: ∂f ∂f +U = 0 ∂t ∂x The characteristic for this equation are: dx df = U; = 0; dt dt Showing that the initial conditions t f are simply advected by a constant velocity U f x WPPII Computational Fluid Dynamics I Finite difference equation ∂f ∂f +U = 0 ∂t ∂x A forward in time, centered in space (FTCS) discretization yields ∆ n+1 n t n n f = f − U ( f + − f − ) j j 2h j 1 j 1 n+1 n j−1 j j+1 WPPII Computational Fluid Dynamics I This scheme is O(∆t, ∆x2) accurate, but a stability analysis shows that the error grows as + ε n 1 U∆t = 1− i sin kh ε n 2h n+1 Since the amplification ε factor has the form 1+i() ε n the absolute value of this complex number is always U∆t i sin kh larger than unity and the 2h method is unconditionally unstable for this case. 1 WPPII Computational Fluid Dynamics I Alternative Scheme: Upwind Difference f j−1 f j f j+1 ∂f ∂f +U = 0 ∂t ∂x U A forward in time but “upwind” (windward) in space discretization yields + ∆t f n 1 = f n − U ( f n − f n ) j j h j j−1 n+1 This scheme is O(∆t, ∆x) n accurate. j−1 j WPPII Computational Fluid Dynamics I To examine the stability we use the von Neuman’s method: n = n + ε n f j f j j Substituting into the modified equation, + ε n 1 −ε n U j j + (ε n −ε n ) = 0 ∆t h j j−1 ε n = ε n ikx j j e ε n+1 −ε n ε n +U (1− e−ikh ) = 0 ∆t h WPPII Computational Fluid Dynamics I ε n+1 −ε n ε n +U (1− e−ikh ) = 0 ∆t h Amplification factor n+1 ε U∆t − U∆t U∆t =1− (1− e ikh ) =1− (1− cos kh) − isin kh ε n h h h or + ε n 1 2U∆t kh U∆t =1− sin 2 − isin kh ε n h 2 h WPPII Computational Fluid Dynamics I Amplification Factor − U∆t G =1− λ(1− e ikh ), λ = h Im(G) Stability Condition: λ <1 1 U∆t Stable ≤ 1 G h kh Re(G) 1 CFL Condition (Courant-Friedrichs-Lewy; 1932) 1−λ λ Computational Fluid Dynamics I WPPII Stability Implication of the CFL Condition (Hirsch, vol.1, p.
Recommended publications
  • Numerical Convection Algorithms and Their Role in Eulerian CFD Reactor Simulations
    International Journal of Chemical Reactor Engineering Volume 1 2003 Article A1 Numerical Convection Algorithms and Their Role in Eulerian CFD Reactor Simulations Hugo A. Jakobsen∗ ∗Norwegian University of Science and Technology, [email protected] Copyright c 2003 by the authors. All rights reserved. No part of this publication may be reproduced, stored in a retrieval system, or transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise, without the prior written permission of the publisher, bepress. Numerical Convection Algorithms and Their Role in Eulerian CFD Reactor Simulations Hugo A. Jakobsen Abstract In this paper a comparative convection algorithm study is presented. The performance of a large number of schemes is compared evaluating the predicted solutions for a standard benchmarking test problem. The nature of the errors caused by the numerical approximations to the convection term is highlighted. Although there is no algorithm that performs the best in general, several conclu- sions can be made. The tests performed show that the 1st order upwind scheme and several variations of this scheme are very diffusive and should be avoided. Most stable 2nd order schemes seem to be much more accurate, whereas the ac- curacy gained by higher order schemes (3rd order and 4th order) may be a little more costly. Implicit time integration schemes are usually not as efficient as the corresponding explicit schemes due to the computational time required on the it- erative process. With larger time steps the accuracy of implicit schemes decrease rapidly. The choice of proper higher order schemes (2nd order schemes) is then seemingly determined by the trade off between accuracy and computational time.
    [Show full text]
  • Upwind Schemes for the Wave Equation in Second-Order Form
    Upwind schemes for the wave equation in second-order form Jeffrey W. Banksa,1,∗, William D. Henshawa,1 aCenter for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore, CA 94551, USA Abstract We develop new high-order accurate upwind schemes for the wave equation in second-order form. These schemes are developed directly for the equations in second-order form, as opposed to transforming the equations to a first-order hyperbolic system. The schemes are based on the solution to a local Riemann-type problem that uses d’Alembert’s exact solution. We construct conservative finite difference approximations, although finite volume approximations are also possible. High-order accuracy is obtained using a space-time procedure which requires only two discrete time levels. The advantages of our approach include efficiency in both memory and speed together with accuracy and robustness. The stability and accuracy of the approximations in one and two space dimensions are studied through normal-mode analysis. The form of the dissipation and dispersion introduced by the schemes is elucidated from the modified equations. Upwind schemes are implemented and verified in one dimension for approximations up to sixth-order accuracy, and in two dimensions for approximations up to fourth-order accuracy. Numerical computations demonstrate the attractive properties of the approach for solutions with varying degrees of smoothness. Keywords: second-order wave equations, upwind discretization, Godunov methods, high-order accurate, finite-difference, finite-volume Contents 1 Introduction 2 2 Construction of an upwind scheme for the second-order wave equation in one dimension 3 3 A general construction for high-order upwind schemes – the one-dimensional case 8 3.1 The upwind flux for the second-order system .
    [Show full text]
  • Mafelap 2013
    MAFELAP 2013 Conference on the Mathematics of Finite Elements and Applications 10{14 June 2013 Abstracts MAFELAP 2013 The organisers of MAFELAP 2013 are pleased to acknowledge the nancial support given to the conference by the Institute of Mathematics and its Applications (IMA) in the form of IMA Stu- dentships. Contents of the MAFELAP 2013 Abstracts Invited, parallel and mini-symposium order Invited talks Finite Element Methods in Coastal Ocean Modeling: Successes and Challenges Clint N Dawson ZIENKIEWICZ LECTURE . 53 Forty years of the Crouzeix-Raviart element Susanne C. Brenner BABUSKAˇ LECTURE . 27 Advances in reducing the mesh burden in computational mechanics applications to fracture and surgical simulation St´ephaneP.A. Bordas ..........................................................................25 Fields, control fields, and commuting diagrams in isogeometric analysis Annalisa Buffa, Giancarlo Sangalli and Rafael V´azquez ..........................................................................30 The inf-sup constant of the divergence Martin Costabel ..........................................................................49 Finite element methods for surface PDEs Charles M. Elliott ..........................................................................64 A stochastic collocation approach to PDE-constrained optimization for random data identification problems Max Gunzburger ..........................................................................98 Time Domain Integral Equations for Computational Electromagnetism Peter
    [Show full text]
  • 65 Numerical Analysis
    e Q (e t o 5 M SectionsSet 1Q (Section 65)MR September 2012 65 NUMERICAL ANALYSIS MR2918625 65-06 FRecent advances in scientific computing and matrix analysis. Proceedings of the International Workshop held at the University of Macau, Macau, December 28{30, 2009. Edited by Xiao-Qing Jin, Hai-Wei Sun and Seak-Weng Vong. International Press, Somerville, MA; Higher Education Press, Beijing, 2011. xii+126 pp. ISBN 978-1-57146-202-2 Contents: Zheng-jian Bai and Xiao-qing Jin [Xiao Qing Jin1], A note on the Ulm-like method for inverse eigenvalue problems (1{7) MR2908437; Che-man Cheng [Che-Man Cheng], Kin-sio Fong [Kin-Sio Fong] and Io-kei Lok [Io-Kei Lok], Another proof for commutators with maximal Frobenius norm (9{14) MR2908438; Wai-ki Ching [Wai-Ki Ching] and Dong-mei Zhu [Dong Mei Zhu1], On high-dimensional Markov chain models for categorical data sequences with applications (15{34) MR2908439; Yan-nei Law [Yan Nei Law], Hwee-kuan Lee [Hwee Kuan Lee], Chao-qiang Liu [Chaoqiang Liu] and Andy M. Yip, An additive variational model for image segmentation (35{48) MR2908440; Hai-yong Liao [Haiyong Liao] and Michael K. Ng, Total variation image restoration with automatic selection of regularization parameters (49{59) MR2908441; Franklin T. Luk and San-zheng Qiao [San Zheng Qiao], Matrices and the LLL algorithm (61{69) MR2908442; Mila Nikolova, Michael K. Ng and Chi-pan Tam [Chi-Pan Tam], A fast nonconvex nonsmooth minimization method for image restoration and reconstruction (71{83) MR2908443; Gang Wu [Gang Wu1], Eigenvalues of certain augmented complex stochastic matrices with applications to PageRank (85{92) MR2908444; Yan Xuan and Fu-rong Lin, Clenshaw-Curtis-rational quadrature rule for Wiener-Hopf equations of the second kind (93{110) MR2908445; Man-chung Yeung [Man-Chung Yeung], On the solution of singular systems by Krylov subspace methods (111{116) MR2908446; Qi- fang Yu, San-zheng Qiao [San Zheng Qiao] and Yi-min Wei, A comparative study of the LLL algorithm (117{126) MR2908447.
    [Show full text]
  • Numerical Methods for the Problem of Spurious Oscillations Arising from the Black-Scholes Equation
    NUMERICAL METHODS FOR THE PROBLEM OF SPURIOUS OSCILLATIONS ARISING FROM THE BLACK-SCHOLES EQUATION Zewen Shen July, 2019 ABSTRACT The Black-Scholes equation is formally parabolic. However, the numerical approximation behaves as if it was hyperbolic when it is convection dominated. Spurious oscillations appear when a conventional finite difference discretization is applied, leading to approximations that violate the financial rule. While the upwind scheme can eliminate the oscillations, the numerical diffusion it introduced impaired the accuracy of the solution. Thus, we illustrated that the total variation diminishing scheme with flux limiters is one of the best solutions to this problem. In the end, we introduced the discontinuous Galerkin method for further study. Acknowledgements I would like to express my deep appreciation to my research project supervisor, Professor Christina Christara, for introducing me to the field of mathematical finance, and also for her patience and guidance during the research. Under her supervision, this research project becomes a valuable learning experience during my undergraduate studies. Finally, I would like to thank my family for their endless support and love. 2 Contents 1 Introduction 3 1.1 Black-Scholes equation . .3 1.2 Motivations . .4 1.3 Spurious Oscillations in the Numerical Solution of the Convection Dominated B-S Equation . .5 2 Upwind Scheme 7 2.1 Discretization . .7 2.2 Numerical Experiments . .7 2.3 The Flaw Behind the Upwind Scheme: Numerical Diffusion . .8 3 Total Variation Diminishing Scheme with van Leer Flux Limiter 9 3.1 Finite Volume Discretization of the Black-Scholes Equation . .9 3.2 Total Variation Diminishing Scheme .
    [Show full text]
  • Mafelap 2013
    MAFELAP 2013 Conference on the Mathematics of Finite Elements and Applications 10{14 June 2013 Abstracts MAFELAP 2013 The organisers of MAFELAP 2013 are pleased to acknowledge the nancial support given to the conference by the Institute of Mathematics and its Applications (IMA) in the form of IMA Stu- dentships. Contents of the MAFELAP 2013 Abstracts Invited, parallel and mini-symposium order Invited talks Finite Element Methods in Coastal Ocean Modeling: Successes and Challenges Clint N Dawson ZIENKIEWICZ LECTURE . 1-1 Forty years of the Crouzeix-Raviart element Susanne C. Brenner BABUSKAˇ LECTURE . 1-1 Advances in reducing the mesh burden in computational mechanics applications to fracture and surgical simulation St´ephaneP.A. Bordas .........................................................................1-2 Fields, control fields, and commuting diagrams in isogeometric analysis Annalisa Buffa, Giancarlo Sangalli and Rafael V´azquez .........................................................................1-4 The inf-sup constant of the divergence Martin Costabel .........................................................................1-5 Finite element methods for surface PDEs Charles M. Elliott .........................................................................1-6 A stochastic collocation approach to PDE-constrained optimization for random data identification problems Max Gunzburger .........................................................................1-7 Time Domain Integral Equations for Computational Electromagnetism Peter
    [Show full text]
  • Numerical Treatment of a Modified Maccormack Scheme in a Nondimensional Form of the Water Quality Models in a Nonuniform Flow Stream
    Hindawi Publishing Corporation Journal of Applied Mathematics Volume 2014, Article ID 274263, 8 pages http://dx.doi.org/10.1155/2014/274263 Research Article Numerical Treatment of a Modified MacCormack Scheme in a Nondimensional Form of the Water Quality Models in a Nonuniform Flow Stream Nopparat Pochai1,2 1 Department of Mathematics, Faculty of Science, King Mongkut’s Institute of Technology Ladkrabang, Bangkok 10520, Thailand 2 Centre of Excellence in Mathematics, Commission on Higher Education (CHE), Si Ayutthaya Road, Bangkok 10400, Thailand Correspondence should be addressed to Nopparat Pochai; nop [email protected] Received 1 September 2013; Revised 17 December 2013; Accepted 17 December 2013; Published 23 February 2014 Academic Editor: Lu´ıs Godinho Copyright © 2014 Nopparat Pochai. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Two mathematical models are used to simulate water quality in a nonuniform flow stream. The first model is the hydrodynamic model that provides the velocity field and the elevation of water. The second model is the dispersion model that provides the pollutant concentration field. Both models are formulated in one-dimensional equations. The traditional Crank-Nicolson method is also used in the hydrodynamic model. At each step, the flow velocity fields calculated from the first model are the input into the second model as the field data. A modified MacCormack method is subsequently employed in the second model. This paper proposes a simply remarkable alteration to the MacCormack method so as to make it more accurate without any significant loss of computational efficiency.
    [Show full text]
  • Large Eddy Simulation of Turbulent Flows Using Finite Volume Methods with Structured, Unstructured, and Zonal Embedded Grids Todd Allen Simons Iowa State University
    Iowa State University Capstones, Theses and Retrospective Theses and Dissertations Dissertations 1998 Large eddy simulation of turbulent flows using finite volume methods with structured, unstructured, and zonal embedded grids Todd Allen Simons Iowa State University Follow this and additional works at: https://lib.dr.iastate.edu/rtd Part of the Mechanical Engineering Commons Recommended Citation Simons, Todd Allen, "Large eddy simulation of turbulent flows using finite volume methods with structured, unstructured, and zonal embedded grids " (1998). Retrospective Theses and Dissertations. 11892. https://lib.dr.iastate.edu/rtd/11892 This Dissertation is brought to you for free and open access by the Iowa State University Capstones, Theses and Dissertations at Iowa State University Digital Repository. It has been accepted for inclusion in Retrospective Theses and Dissertations by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. INFORMATION TO USERS This manuscript has been reproduced from the microfilm master. UMI films the text directly from the original or copy submitted. Thus, some thesis and dissertation copies are in typewriter face, while others may be from any type of computer printer. The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleedthrough, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send UMI a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion.
    [Show full text]
  • Computational Modelling of Flow and Transport
    Course CIE4340 Computational Modelling of Flow and Transport January 2015 Faculty of Civil Engineering and Faculty Geosciences M. Zijlema Delft University of Technology Artikelnummer 06917300083 CIE4340 Computational modelling of flow and transport by : M. Zijlema mail address : Delft University of Technology Faculty of Civil Engineering and Geosciences Environmental Fluid Mechanics Section P.O. Box 5048 2600 GA Delft The Netherlands website : http://fluidmechanics.tudelft.nl Copyright (c) 2011-2015 Delft University of Technology iv Contents 1 Introduction 1 1.1Amathematicalmodel............................. 2 1.2Acomputationalmodel............................. 4 1.3Scopeofthislecturenotes........................... 7 2 Numerical methods for initial value problems 9 2.1Firstorderordinarydifferentialequation................... 9 2.1.1 Timeintegration............................ 10 2.1.2 Convergence, consistency and stability . .............. 14 2.1.3 Two-stepschemes............................ 23 2.2Systemoffirstorderdifferentialequations.................. 26 2.2.1 ThesolutionofsystemoffirstorderODEs.............. 27 2.2.2 Consistencyandtruncationerror................... 33 2.2.3 Stability and stability regions . .................. 36 2.2.4 StiffODEs................................ 39 2.3Concludingremarks............................... 42 3 Numerical methods for boundary value problems 43 3.1IntroductiontoPDEs.............................. 43 3.1.1 ExamplesofPDEs........................... 43 3.1.2 RelevantnotionsonPDEs......................
    [Show full text]
  • A Complete Bibliography of Publications in the Journal of Scientific Computing
    A Complete Bibliography of Publications in The Journal of Scientific Computing Nelson H. F. Beebe University of Utah Department of Mathematics, 110 LCB 155 S 1400 E RM 233 Salt Lake City, UT 84112-0090 USA Tel: +1 801 581 5254 FAX: +1 801 581 4148 E-mail: [email protected], [email protected], [email protected] (Internet) WWW URL: http://www.math.utah.edu/~beebe/ 13 May 2021 Version 1.22 Title word cross-reference −1 [1754]. 0 [1255]. 1 [812, 955, 1192, 1359, 1558, 1641, 1795, 1858, 2652, 2707]. 2 [177, 316, 384, 486, 593, 600, 658, 742, 769, 823, 865, 1174, 1260, 1276, 1306, 1310, 1658, 1748, 1779, 2002, 2252, 2581, 2583, 2597, 2644]. 3 [248, 263, 336, 384, 403, 530, 600, 646, 869, 1002, 1006, 1021, 1061, 1078, 1236, 1276, 1310, 1339, 1579, 1737, 1779, 1780, 1838, 1888, 2107,2529, 2581, 2746, 2755]. −1 4 [2002]. 5 [1629]. 7 [982]. p [1628]. A [276, 480, 629, 2767]. A [629]. α [2740]. Ax = b [276]. Ax = b [228]. B [1612, 1776]. C [1612, 1776]. C0 [440, 1161,1273, 1639, 1729, 1974, 2013]. C1 [578]. D [737, 1400, 2024, 2630]. 2 q ∆ · B [582]. Div(Y )=f [571]. ` [2316]. ` [2316]. `0 [2353]. `1 rot [1190, 1318, 1583, 2039, 2713]. `2 [1190, 1628]. [104, 247]. EQ1 [1195, 1739]. 2 1 f [948]. Fj (x) [278, 283]. f 2 L [571]. fp [790]. G [1373]. H [473, 1194,1213, 1491, 1733, 1951, 1999, 2064, 2253, 2535]. h=p [2205]. H1 [456, 1573,2267]. H1=2 [631]. H2 [742, 1355].
    [Show full text]
  • Numerical Analysis of Convection–Diffusion Using a Modified Upwind
    mathematics Article Numerical Analysis of Convection–Diffusion Using a Modified Upwind Approach in the Finite Volume Method Arafat Hussain 1, Zhoushun Zheng 1,* and Eyaya Fekadie Anley 1,2 1 School of Mathematics and Statistics, Central South University, Changsha 410083, China; [email protected] (A.H.); [email protected] (E.F.A.) 2 Department of Mathematics, College of Natural and Computational Science, Arba-Minch University, Arba-Minch 21, Ethiopia * Correspondence: [email protected] Received: 14 September 2020; Accepted: 5 October 2020; Published: 28 October 2020 Abstract: The main focus of this study was to develop a numerical scheme with new expressions for interface flux approximations based on the upwind approach in the finite volume method. Our new proposed numerical scheme is unconditionally stable with second-order accuracy in both space and time. The method is based on the second-order formulation for the temporal approximation, and an upwind approach of the finite volume method is used for spatial interface approximation. Some numerical experiments have been conducted to illustrate the performance of the new numerical scheme for a convection–diffusion problem. For the phenomena of convection dominance and diffusion dominance, we developed a comparative study of this new upwind finite volume method with an existing upwind form and central difference scheme of the finite volume method. The modified numerical scheme shows highly accurate results as compared to both numerical schemes. Keywords: finite volume method; second-order; convection–diffusion problem; stability; convergence order MSC: 26A33; 35R11; 65M12 1. Introduction Partial differential equations (PDEs) have a vital role in describing different phenomena in real life.
    [Show full text]
  • Numerical Methods for Inviscid Flow Equations
    CHAPTER SIX NUMERICAL METHODS FOR INVISCID FLOW EQUATIONS 6.1 INTRODUCTION The Navier-Stokes equations govern the flows commonly encountered in both internal and external applications. Computing a solution of the Navier-Stokes equations is often difficult or at least impractical and, in many of these applications, unnecessary. Results obtained from a solution of the Euler equations are particularly useful in preliminary design work, where information on pressure alone is desired. In problems where heat transfer and skin friction are required, a solution of the boundary-layer equations usually provides an adequate approximation. However, the outer-edge conditions, including the pressure, must be specified from the inviscid solution as the first step in such an analysis. The Euler equations are also of interest because many of the major elements of fluid dynamics are incorporated in them. For example, fluid flows frequently have internal discontinuities such as shock waves or contact surfaces, Solutions relating the end states across a shock are given by the Rankine- Hugoniot relations; these relations are contained in solutions of the Euler equations. The Euler equations govern the motion of an inviscid nonheat-conducting gas and have a different character in different flow regimes. If the time- dependent terms are retained, the resulting unsteady equations are hyperbolic for all Mach numbers, and solutions can be obtained using time-marching procedures. The situation is very different when a steady flow is assumed. In this case, the Euler equations are elliptic when the flow is subsonic, and hyperbolic 351 352 APPLICATION OF NUMERICAL METHODS when the flow is supersonic.
    [Show full text]