Compact Difference Schemes on a Three-Point Stencil for Second-Order Hyperbolic Equations

We consider compact difference schemes of approximation order \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4+2 $$\end{document} on a three-point spatial stencil for the Klein–Gordon equations with constant and variable coefficients. New compact schemes are proposed for one type of second-order quasilinear hyperbolic equations. In the case of constant coefficients, we prove the strong stability of the difference solution under small perturbations of the initial conditions, the right-hand side, and the coefficients of the equation. A priori estimates are obtained for the stability and convergence of the difference solution in strong mesh norms.


INTRODUCTION
Improving the accuracy of numerical methods for solving problems of mathematical physics on minimal stencils has always been a topical problem in numerical analysis (see, e.g., [1][2][3][4]). A special place among the methods for constructing difference schemes of higher approximation order is occupied by so-called compact schemes, which are written on a stencil that differs insignificantly from the ones traditional for the equation in question [5]. The fundamental paper [6] on this topic for classical equations of mathematical physics with a self-adjoint elliptic operator was published by Samarskii more than 50 years ago. Compact difference schemes for other classes of equations, including convection-diffusion type equations, nonlinear equations without mixed derivatives, and aerohydrodynamic problems were constructed in [7][8][9].
In the present paper, we study compact difference schemes of approximation order 4 + 2 on the usual three-point stencil for various types of the Klein-Gordon equation. This equation plays an important role in mathematical physics; in particular, it is used in the study of solitons and in condensed matter physics [10]. Compact difference schemes for this equation were constructed and studied, e.g., in [9,11]. Some results on the topic were also announced in [12,13]. For this equation, although the differential and difference problems with variable coefficients are linear, one fails to obtain relevant a priori estimates by applying well-known results from Samarskii's stability theory of three-level operator-difference schemes [1,Ch. VI,Sec. 3]. In the present paper, we apply the method of energy inequalities to compact difference schemes approximating the Klein-Gordon equations with variable coefficients to obtain a priori estimates for the stability and convergence of the difference solution in the mesh norms of L 2 (ω h ), W 1 2 (ω h ), C(ω h ), or L ∞ (ω h ). We use a numerical experiment for a quasilinear equation that is a corollary of the system of gasdynamic equations as an example shown how to use the Runge rule for determining various orders of the convergence rate of the solution of a difference scheme in the case of two independent variables.

NECESSARY AND SUFFICIENT CONDITIONS FOR THE STABILITY OF TWO-AND THREE-LEVEL OPERATOR-DIFFERENCE SCHEMES
When studying compact difference schemes approximating linear Klein-Gordon equations, it is natural to use Samarskii's general theory [1, Ch. VI, Sec. 3] of operator-difference schemes. Below, we suggest to use some other canonical forms of operator-difference schemes, for which the stability conditions are much simpler and do not contain conditions on the relationship between the operators. Consider a real finite-dimensional Euclidean space H and a time mesh ω τ = {t n = nτ , n = 0, . . . , N 0 , τ N 0 = T } = ω τ ∪ {0}. Let the inner product on H be denoted by (·, ·), and let A, B, D : H → H be linear operators independent of τ and t n .
Consider the Cauchy problem for the two-level operator-difference scheme where y n = y(t n ) ∈ H is the desired vector function, ϕ n = ϕ(t n ) and u 0 are given, and y (σ) = σy n+1 + (1 − σ)y n , y n ∈ H. Then the Samarskii criterion (see [1, p. 333]) can be stated for ϕ(t) ≡ 0 as follows.
are necessary and sufficient for the stability of the solution of the difference problem (1), (2) in H A with respect to the initial data, i.e., for the estimate We will use the following canonical form for three-level operator-difference schemes: where In the present paper, we use the notation in [1,2]. In what follows, we assume that D and A 1 are positive, self-adjoint, and constant operators, Then the following assertion holds.
Theorem 2. The difference scheme (3), (4) is stable in H A1 with respect to the initial data and the right-hand side, and one has the a priori estimate in the sense of the inner product on H, we arrive at the energy relation This implies the desired estimate.

Statement of the Problem and the Difference Scheme
In the domain Q T = {(x, t) : 0 ≤ x ≤ l, 0 ≤ t ≤ T }, we consider the initial-boundary value problem for the Klein-Gordon equation with constant coefficients u(0, t) = µ 1 (t), u(l, t) = µ 2 (t). (8) Note that Eq. (6) is a generalization of the wave equation and is used to describe rapidly moving particles with nonzero rest mass. From now on, we assume that there exists a unique solution of the differential problem and that all of its derivatives needed in the exposition are continuous in Q T .
On the uniform mesh of grid points , we replace the differential problem by the difference problem whereŷ = y n+1 ,y = y n−1 , y n i = y(x i , t n ), Just as in the monograph [1, p. 309], one can readily show that the discrepancy Λu + ϕ and the approximation error in the second initial condition satisfy the a priori estimates i.e., the difference scheme (9)-(11) approximates the original differential problem with the fourth order in space and the second order in time.

Stability with Respect to the Initial Data and the Right-Hand Side
To study these issues in the linear case, one usually applies the theory of three-level operatordifference schemes developed by Samarskii [1, Ch. VI, Sec. 3]. The requirement for the main spatial operator to be self-adjoint, results in stringent restrictions on the homogeneity of the boundary conditions. To avoid this, consider the perturbed solutionỹ produced by the difference scheme (9)-(11) with perturbed righthand sidef and perturbed initial conditionsũ 0 ,ũ 1 . Then the problem for the perturbation y =ỹ −y can be written in the operator form Here The operator A defined by relation (14) is constant, positive, and self-adjoint, 0 < A * = A : H → H; here H is the space of mesh functions defined on ω h and vanishing for x = 0 and x = l. Therefore, the operators D and A 1 are positive, constant, and self-adjoint as well, D = D * > E and A 1 = A * 1 > A. Note the following well-known properties of the operator A [1, Ch. II, Sec. 3]: We need the following assertion.
Assume that the operators D and A 1 in the canonical form (15), (16) are constant, positive, and self-adjoint in H and satisfy the inequality Then the solution of the scheme (15), (16) satisfies the a priori estimate Let us apply the estimate (21) to the stability analysis of the compact difference scheme (9)- (11). In what follows, we assume that Then condition (19) is satisfied. Indeed, Thus, the following assertion holds.
Theorem 3. Let condition (22) be satisfied. Then the difference scheme with the inhomogeneous boundary conditions (9)-(11) is conditionally stable with respect to the initial data and the right-hand side, and its solution satisfies the a priori estimate Remark. In the linear case, the change of variables reduces the original differential problem to a problem with homogeneous boundary conditions for v. Further, the stability of the solution of difference schemes for these problems with respect to the boundary conditions and the right-hand side takes place simultaneously for both problems. However, this change of variables may fail to produce the desired result for the nonlinear Klein-Gordon equation. Thus, we have proved the stability of solution of the difference scheme in the weak L 2 (ω h ) energy norm under rather weak constraints on the function weight σ. Such estimates are useful when studying the convergence of difference schemes with generalized solutions [14]. To produce a priori estimates in the stronger W 1 2 (ω h ) and C(ω h ) norms, we will use Theorem 2 and the embedding (see [1, p. 107 Let us write the difference scheme for the perturbation y =ỹ − y in the canonical form where Then the following assertion holds. is stable with respect to the initial data and the right-hand side, and for all n = 2, . . . , N 0 one has the a priori estimates (5) is satisfied; therefore, by Theorem 2, Since A 1 > A and D −1 ≤ 2E, we have the estimate Consequently, it follows from the lemma and the embedding (23) that The proof of the theorem is complete.

Strong Stability
When studying the well-posedness of difference schemes, attention is mainly paid to the stability of the solution with respect to the initial data and the right-hand side [1,2]. However, when solving a differential problem numerically, it may turn out that the coefficients of the equation are given approximately rather than exactly. This shows how important it is to study schemes with perturbed coefficients. By strong stability we mean the stability of the solution of the difference problem with respect to small perturbations in the initial conditions, right-hand side, and coefficients of the equation [15].
Along with the difference scheme (9)-(11), consider the corresponding perturbed problem We subtract Eqs. (9)-(11) from the respective equations (27)-(29) and obtain a problem for the perturbation y =ỹ − y, which we write in the operator form (15), (16), Here D = D * > E andÃ 1 =Ã * 1 > A ≥ δE are positive constant self-adjoint operators. Note that condition (19) is satisfied for h/ √ 6 ≤ τ ≤ 1/m; therefore, we obtain the inequality for y based on the a priori estimate (21). By virtue of relations (17) and (18), we have A < 4/h 2 , and consequently, the estimate (21) implies the estimate Thus, we are now in a position to state the strong stability theorem.
Theorem 5. Assume that the grid increments satisfy the condition Then the solution of the difference scheme (9)-(11) is strongly stable, and the perturbation of the solution satisfies the a priori estimate K s |m − m| , n = 1, . . . , N 0 − 1.
Since problems (9)-(11) and (30)-(32) are identical, we can use Theorem 4 to estimate the method error. Proof. Indeed, the a priori estimates (12), (13), and (26) imply the inequality Consequently, the difference solution converges to the exact solution with the fourth order in space and the second order in time.

Stability
To avoid cumbersome calculations, we restrict ourselves to the case in which the coefficient k = k(x) depends on the spatial variable alone. Consider the perturbed solutionỹ produced by the difference scheme (36)-(38) with perturbed right-hand sidef and perturbed initial conditionsũ 0 andũ 1 . Then the problem for the perturbation y =ỹ − y acquires the form Here u 0 =ũ 0 − u 0 , u 1 =ũ 1 − u 1 , and ϕ =φ − ϕ.
Unfortunately, although the difference problem is linear, the Samarskii theory of three-level operator-difference schemes [1] does not apply. We will use the method of energy inequalities; in addition to some well-known facts from the theory of difference schemes such as the first difference Green's formula, the Cauchy-Schwarz inequality with ε [1, Ch. II, Sec. 3], and the difference analog of Gronwall's lemma [16, Ch. III, Sec. 1], we also need the following assertion.

Lemma 2. The expression
Proof. It suffices to show that the expression is nonnegative. In view of the obvious inequalities Under the first condition in (42), the inequality holds. Consequently, we have under the second condition in (42). The proof of the lemma is complete.
To obtain an a priori estimate for y, we multiply the difference equation (39) by 2τ y• t in the sense of the inner product and apply the first Green's difference formula to obtain the energy relation Consider the terms other than Q n and Q n+1 in (44). We apply the Cauchy-Schwarz inequality with ε and readily obtain the estimates where c > 0 is a case-specific constant depending on m, ε, and max We take into account inequalities (45)-(47) in (44) and, under conditions (42), arrive at the recursion relation Q n+1 ≤ (1 + τ c)Q n + τ c ϕ 2 ≤ e cτ Q n + τ c ϕ 2 .
Thus, the following assertion holds.
Theorem 7. Let the condition τ ≥ max 1, 2 3k 1 h be satisfied. Then one has the estimate which implies the ρ-stability of the solution of the difference scheme (36)-(38) with respect to the initial data and the right-hand side in the L 2 (ω h ), W 1 2 (ω h ), and C(ω h ) mesh norms. The proof of the theorem follows from inequality (48), the Gronwall lemma, and the embed-

Convergence of the Difference Scheme in the C(ω h ) Mesh Norm
We substitute z + u for y into the difference equations (36)-(38), where u is the solution of problem (33)-(35), and obtain the problem for the error z. Problems (50)-(52) and (39)-(41) are identical. Therefore, we can apply Theorem 7 to estimate the method error. Then, in accordance with inequality (49), we obtain the estimate where M 1 = const > 0. Thus, we are now in a position to formulate the convergence theorem. where M 2 = const > 0.

QUASILINEAR KLEIN-GORDON EQUATIONS
For equations with variable coefficients, the compact difference schemes proposed by Samarskii have the disadvantage that it is impossible to generalize them to the case of quasilinear equations, because the corresponding stencil functional must be calculated at a half-integer point, which does not exist in the quasilinear case. Nevertheless, compact schemes of order 4+2, similar to schemes for the case of constant coefficients, can also be constructed for the quasilinear Klein-Gordon equations with the condition φ u = k(u) ≥ k 1 > 0.
The difference scheme of approximation order 4 + 2 on the standard stencil has the form To implement this scheme, one needs to use the Newton iteration method.

COMPUTATIONAL EXPERIMENT
Here we present the results of numerical calculations obtained when solving an initial-boundary value problem for the equation with parameter values γ = 5, l = 2, and T = 2. The initial and boundary conditions are determined from the exact solution .
Equation (53) is a corollary of the system of equations of gas dynamics in Lagrangian variables with the equation of state for a polytropic gas, where u = η is the specific volume.
On each layer, this process stops when the following condition is satisfied for some M : The order of convergence in the time and space variables in the L ∞ = C norm is determined by the formulas Since the difference solution converges to the exact solution with the fourth order in space and the second order in time, we select steps h and τ satisfying the inequality h 4 ≤ τ 2 to verify the rate   of convergence in the time variable. Then we obtain an O(τ 2 ) scheme and deal with the second Runge rule (55). In a similar way, when considering the order with respect to h, we take care that the inequality h 4 ≥ τ 2 holds in the calculations. Then we can apply the first Runge rule (55).
Tables 1 and 2 list the rates of convergence of the approximate solution to the exact one.
Hence we see that the constructed difference scheme has the fourth order of accuracy in the spatial variable and the second order of accuracy in the time variable.
In addition, in Figs. 1 and 2 one can clearly see by colors and figures that the approximate solution converges to the exact solution as the grid increments h and τ are refined. For convenience of visual observation, the above results were produced in the domain 0 ≤ x ≤ 10, 0 ≤ t ≤ 10.
The computational experiment confirms our theoretical conclusions.

OPEN ACCESS
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.