Open AccessON COMPACT 4TH ORDER FINITE-DIFFERENCE SCHEMES FOR THE WAVE EQUATION
Author(s) -
А. А. Злотник,
Olga Kireeva
Publication year - 2021
Publication title -
mathematical modelling and analysis/mathematical modeling and analysis
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.491
H-Index - 25
eISSN - 1648-3510
pISSN - 1392-6292
DOI - 10.3846/mma.2021.13770
Subject(s) - mathematics , smoothness , corollary , mathematical analysis , stability (learning theory) , term (time) , operator (biology) , polygon mesh , boundary value problem , order (exchange) , wave equation , finite difference , initial value problem , pure mathematics , geometry , computer science , biochemistry , chemistry , physics , finance , repressor , quantum mechanics , machine learning , transcription factor , economics , gene
We consider compact finite-difference schemes of the 4th approximation order for an initial-boundary value problem (IBVP) for the n-dimensional nonhomogeneous wave equation, n≥ 1. Their construction is accomplished by both the classical Numerov approach and alternative technique based on averaging of the equation, together with further necessary improvements of the arising scheme for n≥ 2. The alternative technique is applicable to other types of PDEs including parabolic and time-dependent Schro¨dinger ones. The schemes are implicit and three-point in each spatial direction and time and include a scheme with a splitting operator for n≥ 2. For n = 1 and the mesh on characteristics, the 4th order scheme becomes explicit and close to an exact four-point scheme. We present a conditional stability theorem covering the cases of stability in strong and weak energy norms with respect to both initial functions and free term in the equation. Its corollary ensures the 4th order error bound in the case of smooth solutions to the IBVP. The main schemes are generalized for non-uniform rectangular meshes. We also give results of numerical experiments showing the sensitive dependence of the error orders in three norms on the weak smoothness order of the initial functions and free term and essential advantages over the 2nd approximation order schemes in the non-smooth case as well.