Abstract
The second order of accuracy difference scheme generated by CrankNicholson difference scheme for approximately solving multipoint nonlocal boundary value problem is considered. Wellposedness of this difference scheme in Hölder spaces is established. Furthermore, as applications, coercivity estimates in Hölder norms for approximate solutions of the multipoint nonlocal boundary value problems for mixed type equations are obtained. Moreover, the method is illustrated by numerical examples.
Keywords:
difference scheme; ellipticparabolic equation; wellposedness1 Introduction
In recent years, more and more mathematicians have been studying nonlocal problems for ordinary differential equations and partial differential equations because of their existence in many applied problems included in applied sciences. Theory and numerical methods of solutions of the nonlocal boundary value problems for these partial differential equations were investigated by many researchers (see, e.g., [113] and the references therein). Several types of problems in fluid mechanics, other areas of physics, and mathematical biology led to partial differential equations of ellipticparabolic type (see, [1418]). The purpose of this paper is to study the second order of accuracy difference schemes of ellipticparabolic problem with nonlocal boundary value problems.
In [19], we established the wellposedness of multipoint nonlocal boundary value problem
in a Hilbert space H with the selfadjoint positive definite operator A under assumption
The wellposedness of multipoint nonlocal boundary value problem (1) in Hölder spaces with a weight was established. In applications, coercivity inequalities for the solutions of nonlocal boundary value problems for ellipticparabolic equations were obtained.
In [20], we studied the wellposedness of the first order of accuracy difference scheme for the approximate solution of boundary value problem (1) under assumption (2).
In the present paper, we consider the second order of accuracy difference scheme generated by CrankNicholson difference scheme
for the approximate solution of boundary value problem (1) under assumption (2).
The wellposedness of difference scheme (3) in Hölder spaces is established. In applications, the stability, almost coercivity stability, coercivity stability estimates for solutions of the second order of accuracy difference scheme for ellipticparabolic equations are obtained. Furthermore, the theoretical statements for the solution of the first and second order of accuracy schemes for onedimensional ellipticparabolic differential equation are supported by the results of a numerical example.
2 Main theorems
Let us give some auxiliary lemmas we need below. Throughout the paper, H is a Hilbert space and we denote , where A is a selfadjoint positive definite operator. Then, it is clear that B is a selfadjoint positive definite operator and , where , and which is defined on the whole space H is a bounded operator. The following operators
and
exist and are bounded for a selfadjoint positive operator A. Here,
and I is the identity operator.
Lemma 2.1For any, and, , the solution of problem (3) exists and the following formulas hold:
Proof Clearly, the solution formula of the problem
is [22]:
for any and γ. Equation (9) and the fact that yield Equation (6).
The solution of the problem
satisfies the following formula [21]:
Equation (5) follows from Equations (9) and (11), initial condition , and
Finally, let us obtain formula (7). Combining (5), (6), and the condition
we get
From Equation (4), it follows that
This finishes the proof of Lemma 2.1. □
Here, we study wellposedness of problem (3). First, we give some necessary estimates for , and .
Lemma 2.2For a selfadjoint positive operatorAthe following estimates are satisfied[21,22,24]:
From these estimates, it follows that
Now, we study wellposedness of problem (3). Let be the linear space of mesh functions defined on with values in the Hilbert space H. Next, on we denote , , , , , and , Banach spaces with the norms
Theorem 2.1Nonlocal boundary value problem (3) is stable innorm.
Proof By [21], we have
for the solution of boundary value problem (10).
By [22], we have
for the solution of inverse Cauchy difference problem (8).
Then, the proof of Theorem 2.1 is based on stability inequalities (14), (15), and on estimates
for the solution of boundary value problem (3). Estimates (16) and (17) follow from formula (7) and estimates (12), (13). Theorem 2.1 is proved. □
Theorem 2.2Assume thatand. Then, for the solution of difference problem (3), we have the following almost coercivity inequality:
whereis independent not only of, , φbut also ofτ.
Proof By [24], we have
for the solution of inverse Cauchy difference problem (8).
By [21], we have
for the solution of boundary value problem (10). Then, the proof of Theorem 2.2 is based on almost coercivity inequalities (18), (19), and on the estimates
for the solution of boundary value problem (3). The proof of these estimates follows the scheme of papers [21,24] and relies on both formula (7) and estimates (12), (13). This finalizes the proof of Theorem 2.2. □
Theorem 2.3Let the assumptions of Theorem 2.2 be satisfied. Then, boundary value problem (3) is wellposed in Hölder spaces, and the following coercivity inequalities hold:
whereMis independent of, , φ, τandα.
Proof By [24],
for the solution of inverse Cauchy difference problem (8).
By [21], we have
for the solution of boundary value problem (10). Then, the proof of Theorem 2.3 is based on coercivity inequalities (20)(22) and estimates
for the solution of difference scheme (3). The proof of these estimates follows the scheme of the papers [21,24] and relies on both formula (7) and estimates (12), (13). This is the end of the proof of Theorem 2.3. □
3 Application
Now, the application of the abstract result is considered. In , let us consider the boundary value problem for multidimensional ellipticparabolic equation
where (), (, ), (, ), and (, ) are given smooth functions. Here, Ω is the unit open cube in the ndimensional Euclidean space (, ) with boundary S, , and .
The discretization of problem (23) is carried out in two steps. In the first step, the grid sets
are defined. To the differential operator A generated by problem (23), we assign the difference operator by formula
acting in the space of grid functions , satisfying the conditions for all . With the help of , we arrive at the nonlocal boundary value problem
for an infinite system of ordinary differential equations (see [21]).
Secondly, problem (25) is replaced by difference scheme (3), so that the following second order of accuracy difference scheme
To formulate the results, we introduce the spaces , , and of the grid functions defined on , equipped with the norms
and
Theorem 3.1Letτandbe sufficiently small positive numbers. Then, solutions of difference scheme (26) satisfy the following stability and almost coercivity estimates:
whereMis independent not only ofτ, h, but also of, and, .
The proof of Theorem 3.1 is based on Theorem 2.1, Theorem 2.2, the symmetry properties of the difference operator defined by formula (24) in , the estimate
and the following theorem in :
Theorem 3.2For the solution of the elliptic difference problem
the following coercivity inequality holds[23]:
Theorem 3.3Letτandbe sufficiently small positive numbers. Then, the solutions of difference scheme (26) satisfy the following coercivity stability estimates:
Here, Mis independent not only ofτ, h, but also of, and, .
The proof of Theorem 3.3 is based on the abstract Theorem 2.3, Theorem 3.2, and the symmetry properties of the difference operator defined by formula (24).
4 Numerical Analysis
The theoretical statements for the solution of these difference schemes are supported by the results of numerical experiments of the nonlocal boundary value problem
for the ellipticparabolic equation. The exact solution of this problem is
For the comparison, the errors computed by the following formula
are recorded for different values of N and M, where represents the exact solution and represents the numerical solution at . The results are shown in Table 1 for respectively.
Table 1. Error analysis for the solution
Therefore, the results indicate that the second order of accuracy difference scheme is more accurate than the first order of accuracy difference scheme.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
The manuscript was drafted by OG and it is based on his PhD thesis. AA is the supervisor of the thesis and gave detailed comments on the manuscript. All authors read and approved the final manuscript.
Acknowledgement
The authors are very grateful to Prof. P. E. Sobolevskii (Jerusalem, Israel) for valuable comments to the improvement of this article.
References

Salakhitdinov, MS: Equations of MixedComposite Type, FAN, Tashkent (1974)

Chipot, M, Lovat, B: Existence and uniqueness results for a class of nonlocal elliptic and parabolic problems. Dyn. Contin. Discrete Impuls. Syst.. 8(1), 35–51 (2001)

Gulin, AV, Ionkin, NI, Morozova, VA: Stability of a nonlocal twodimensional finitedifference problem. Differ. Equ.. 37(7), 970–978 (2001). Publisher Full Text

Agarwal, RP, Bohner, M, Shakhmurov, VB: Maximal regular boundary value problems in Banachvalued weighted space. Bound. Value Probl.. 1, 9–42 (2005)

Gordeziani, DG, Avalishvili, GA: Timenonlocal problems for Schrödingertype equations: I. Problems in abstract spaces. Differ. Equ.. 41(5), 703–711 (2005). Publisher Full Text

Gordeziani, N, Natalini, P, Ricci, PE: Finitedifference methods for solution of nonlocal boundary value problems. Comput. Math. Appl.. 50, 1333–1344 (2005). Publisher Full Text

Berdyshev, AS, Karimov, ET: Some nonlocal problems for the parabolichyperbolic type equation with noncharacteristic line of changing type. Cent. Eur. J. Math.. 4(2), 183–193 (2006). Publisher Full Text

MartínVaquero, J, VigoAguiar, J: On the numerical solution of the heat conduction equations subject to nonlocal conditions. Appl. Numer. Math.. 59, 2507–2514 (2009). Publisher Full Text

Ashyralyev, A, Gercek, O: Finite difference method for multipoint nonlocal ellipticparabolic problems. Comput. Math. Appl.. 2010(7), 2043–2052 (2010)

Gordeziani, DG, Berikelashvili, GK: On a nonlocal generalization of the biharmonic Dirichlet problem. Differ. Equ.. 46(3), 321–328 (2010). Publisher Full Text

Gavrilyuk, IP, Makarov, VL, Sytnyk, DO, Vasylyk, VB: Exponentially convergent method for the mpoint nonlocal problem for a first order differential equation in Banach space. Numer. Funct. Anal. Optim.. 31(1), 1–21 (2011)

Jangveladze, TA, Lobjanidze, GB: On a nonlocal boundary value problem for a fourthorder ordinary differential equation. Differ. Equ.. 47(2), 335–343 (2011)

MartínVaquero, J, QueirugaDios, A, Encinas, AH: Numerical algorithms for diffusionreaction problems with nonclassical conditions. Appl. Math. Comput.. 218(9), 5487–5495 (2012). Publisher Full Text

Vragov, VN: Boundary Value Problems for Nonclassical Equations of Mathematical Physics. Textbook for Universities, NSU Publishing house, Novosibirsk (1983)

Kröner, D, Rodrigues, JF: Global behaviour for bounded solutions of a porous media equation of ellipticparabolic type. J. Math. Pures Appl.. 64(2), 105–120 (1985)

Nakhushev, AM: Equations of Mathematical Biology. Textbook for Universities, Vysshaya Shkola, Moskow (1995)

Ashyralyev, A, Soltanov, H: On one difference schemes for an abstract nonlocal problem generated by the investigation of the motion of gas on the homogeneous space. In: Modeling Processes of Exploitation of Gas Places and Applied Problems of Theoretical Gasohydrodynamics, pp. 147154. Ashgabat (1998)

Diaz, JI, Lerena, MB, Padial, JF, Rakotoson, JM: An ellipticparabolic equation with a nonlocal term for the transient regime of a plasma in a stellarator. J. Differ. Equ.. 198(2), 321–355 (2004). Publisher Full Text

Gercek, O: Difference schemes of nonlocal boundary value problems for ellipticparabolic differential equations. PhD thesis, Yildiz Teknik University (2010)

Ashyralyev, A, Gercek, O: On multipoint nonlocal ellipticparabolic difference problems. Vestn. Odessa Nat. Univ. Math. Mech.. 15(19), 135–156 (2010)

Sobolevskii, PE: The coercive solvability of difference equations. Dokl. Akad. Nauk SSSR. 201(5), 1063–1066 (1971)

Sobolevskii, PE: On the stability and convergence of the CrankNicolson scheme. VariationalDifference Methods in Mathematical Physics, pp. 146–151. Vychisl. Tsentr Sibirsk. Otdel. Akad. Nauk SSSR, Novosibirsk (1974)

Sobolevskii, PE: Difference Methods for the Approximate Solution of Differential Equations, Voronezh State University Press, Voronezh (1975)

Ashyralyev, A, Sobolevskii, PE: Coercive stability of a CrankNicolson difference scheme in $\widetilde{C}_{0}^{\alpha}$C˜0α spaces. In: Approximate Methods for Investigations of Differential Equations and Their Applications, pp. 1624. Kuybishev (1982)