This book covers high order finite difference methods for time dependent PDE. It gives an overview of the basic theory and construction principles by using model examples. The book also contains a general presentation of the techniques and results for well-posedness and stability, with inclusion of the three fundamental methods of analysis both for PDE in its original and discretized form: the Fourier transform, the eneregy method and the Laplace transform.
This book introduces finite difference methods for both ordinary differential equations (ODEs) and partial differential equations (PDEs) and discusses the similarities and differences between algorithm design and stability analysis for different types of equations. A unified view of stability theory for ODEs and PDEs is presented, and the interplay between ODE and PDE analysis is stressed. The text emphasizes standard classical methods, but several newer approaches also are introduced and are described in the context of simple motivating examples.
Praise for the First Edition ". . . fills a considerable gap in the numerical analysis literature by providing a self-contained treatment . . . this is an important work written in a clear style . . . warmly recommended to any graduate student or researcher in the field of the numerical solution of partial differential equations." —SIAM Review Time-Dependent Problems and Difference Methods, Second Edition continues to provide guidance for the analysis of difference methods for computing approximate solutions to partial differential equations for time-dependent problems. The book treats differential equations and difference methods with a parallel development, thus achieving a more useful analysis of numerical methods. The Second Edition presents hyperbolic equations in great detail as well as new coverage on second-order systems of wave equations including acoustic waves, elastic waves, and Einstein equations. Compared to first-order hyperbolic systems, initial-boundary value problems for such systems contain new properties that must be taken into account when analyzing stability. Featuring the latest material in partial differential equations with new theorems, examples, and illustrations,Time-Dependent Problems and Difference Methods, Second Edition also includes: High order methods on staggered grids Extended treatment of Summation By Parts operators and their application to second-order derivatives Simplified presentation of certain parts and proofs Time-Dependent Problems and Difference Methods, Second Edition is an ideal reference for physical scientists, engineers, numerical analysts, and mathematical modelers who use numerical experiments to test designs and to predict and investigate physical phenomena. The book is also excellent for graduate-level courses in applied mathematics and scientific computations.
The book contains a selection of high quality papers, chosen among the best presentations during the International Conference on Spectral and High-Order Methods (2014), and provides an overview of the depth and breadth of the activities within this important research area. The carefully reviewed selection of papers will provide the reader with a snapshot of the state-of-the-art and help initiate new research directions through the extensive biography.
The book contains a selection of high quality papers, chosen among the best presentations during the International Conference on Spectral and High-Order Methods (2009), and provides an overview of the depth and breadth of the activities within this important research area. The carefully reviewed selection of the papers will provide the reader with a snapshot of state-of-the-art and help initiate new research directions through the extensive bibliography.
The book introduces modern high-order methods for computational fluid dynamics. As compared to low order finite volumes predominant in today's production codes, higher order discretizations significantly reduce dispersion errors, the main source of error in long-time simulations of flow at higher Reynolds numbers. A major goal of this book is to teach the basics of the discontinuous Galerkin (DG) method in terms of its finite volume and finite element ingredients. It also discusses the computational efficiency of high-order methods versus state-of-the-art low order methods in the finite difference context, given that accuracy requirements in engineering are often not overly strict. The book mainly addresses researchers and doctoral students in engineering, applied mathematics, physics and high-performance computing with a strong interest in the interdisciplinary aspects of computational fluid dynamics. It is also well-suited for practicing computational engineers who would like to gain an overview of discontinuous Galerkin methods, modern algorithmic realizations, and high-performance implementations.
This book features a selection of high-quality papers chosen from the best presentations at the International Conference on Spectral and High-Order Methods (2016), offering an overview of the depth and breadth of the activities within this important research area. The carefully reviewed papers provide a snapshot of the state of the art, while the extensive bibliography helps initiate new research directions.
In this thesis, we use finite difference operators with the Summation-By-Partsproperty (SBP) and a weak boundary treatment, known as SimultaneousApproximation Terms (SAT), to construct high-order accurate numerical schemes.The SBP property and the SAT’s makes the schemes provably stable. The numerical procedure is general, and can be applied to most problems, but we focus on hyperbolic problems such as the shallow water, Euler and wave equations. For a well-posed problem and a stable numerical scheme, data must be available at the boundaries of the domain. However, there are many scenarios where additional information is available inside the computational domain. In termsof well-posedness and stability, the additional information is redundant, but it can still be used to improve the performance of the numerical scheme. As a first contribution, we introduce a procedure for implementing additional data using SAT’s; we call the procedure the Multiple Penalty Technique (MPT). A stable and accurate scheme augmented with the MPT remains stable and accurate. Moreover, the MPT introduces free parameters that can be used to increase the accuracy, construct absorbing boundary layers, increase the rate of convergence and control the error growth in time. To model infinite physical domains, one need transparent artificial boundary conditions, often referred to as Non-Reflecting Boundary Conditions (NRBC). In general, constructing and implementing such boundary conditions is a difficult task that often requires various approximations of the frequency and range of incident angles of the incoming waves. In the second contribution of this thesis,we show how to construct NRBC’s by using SBP operators in time. In the final contribution of this thesis, we investigate long time error bounds for the wave equation on second order form. Upper bounds for the spatial and temporal derivatives of the error can be obtained, but not for the actual error. The theoretical results indicate that the error grows linearly in time. However, the numerical experiments show that the error is in fact bounded, and consequently that the derived error bounds are probably suboptimal.
