Skip to main content

On a difference scheme of second order of accuracy for the Bitsadze-Samarskii type nonlocal boundary-value problem

Abstract

In this study, the Bitsadze-Samarskii type nonlocal boundary-value problem with integral condition for an elliptic differential equation in a Hilbert space H with self-adjoint positive definite operator A is considered. The second order of the accuracy difference scheme for the approximate solutions of this nonlocal boundary-value problem is presented. The well-posedness of this difference scheme in Hölder spaces with a weight is proved. The theoretical statements for the solution of this difference scheme are supported by the results of numerical example.

1 Introduction

In 1969 Bitsadze and Samarskii [1] stated and studied a new problem in which a nonlocal condition is related to the values of the solution on parts of the boundary and on an interior curve for a uniformly elliptic equation. Furthermore, in [216] the Bitsadze-Samarskii type nonlocal boundary-value problems were investigated for the various differential and difference equations of elliptic type. The role played by coercive inequalities in the study of local boundary-value problems for elliptic differential equations is well known [17]. Methods of solutions of elliptic differential and difference equations have been studied extensively by many researchers (see [1827] and the references therein). In the present paper we consider the Bitsadze-Samarskii type nonlocal boundary-value problem with integral condition,

{ d 2 u ( t ) d t 2 + A u ( t ) = f ( t ) , 0 < t < 1 , u ( 0 ) = φ , u ( 1 ) = 0 1 ρ ( λ ) u ( λ ) d λ + ψ
(1)

for the differential equation of elliptic type in a Hilbert space H with the self-adjoint positive definite operator A with a closed domain D(A)H. Here, let f(t) be a given abstract continuous function defined on [0,1] with values in H, φ, and ψ are elements of D(A) and ρ(t) is a scalar function. A function u(t) is called a solution of problem (1) if the following conditions are satisfied:

  1. i.

    u(t) is a twice continuously differentiable on the segment [0,1].

  2. ii.

    The element u(t) belongs to D(A) for all t[0,1], and the function Au(t) is continuous on the segment [0,1].

  3. iii.

    u(t) satisfies the equation and nonlocal boundary conditions (1).

The paper is organized as follows. In Section 2 the second order of the accuracy difference scheme for the approximate solution (1) is presented. The stability, the almost coercive stability, and the coercive stability estimates for the solution of the difference scheme for an approximate solution of the nonlocal boundary-value problem with integral condition for elliptic equations are obtained. Section 3 contains the applications of Section 2. The final section is devoted to the numerical result. Theoretical statements for the solution of the second order of the accuracy difference scheme is supported by a numerical experiment.

2 The second order of the accuracy difference scheme

Let us associate the nonlocal boundary-value problem (1) with the corresponding difference problem,

{ 1 τ 2 [ u k + 1 2 u k + u k 1 ] + A u k = φ k , φ k = f ( t k ) , t k = k τ , 1 k N 1 , N τ = 1 , u 0 = φ , u N = j = 1 N ρ ( t j τ 2 ) ( u j + u j 1 2 ) τ + ψ .
(2)

We will study the problem (2) under the following assumption:

j = 1 N |ρ ( t j τ 2 ) |τ<1.
(3)

It is well known [28] that for a self-adjoint positive definite operator A it follows that B= 1 2 (τA+ 4 A + τ 2 A 2 ) is self-adjoint positive definite and R= ( I + τ B ) 1 , which is defined on the whole space H is a bounded operator. Here, I is the unit operator. Furthermore, we have

{ ( I R 2 N ) 1 H H M ( δ ) , R k H H M ( δ ) ( 1 + δ τ ) k , k τ B R k H H M ( δ ) , k 1 , δ > 0 , B β ( R k + r R k ) H H M ( δ ) ( r τ ) α ( k τ ) α + β , 1 k < k + r N , 0 α , β 1 .
(4)

In this paper, positive constants, which can differ in time (hence they are not a subject of precision considerations) will be indicated with M. On the other hand M(α,β,) is used to focus on the fact that the constant depends only on α,β, .

Lemma 1 The operator

I j = 1 N ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( R N j R N + j + R N j + 1 R N + j 1 )

has an inverse

S τ = ( I j = 1 N ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( R N j R N + j + R N j + 1 R N + j 1 ) ) 1

and the following estimate is satisfied:

S τ H H M(δ)τ,
(5)

where M does not depend on τ.

The proof of the estimate (5) is based on the estimate

( I L ) u , u ( 1 j = 1 N ρ ( t j τ 2 ) τ ) u,u.
(6)

Here

L= j = 1 N ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( R N j R N + j + R N j + 1 R N + j 1 ) .

The estimate (6) follows from the spectral representation of A and the Cauchy inequality.

Theorem 2 For any φ k , 1kN1, the solution of the problem (2) exists and the following formula holds:

u k = ( I R 2 N ) 1 { ( R k R 2 N k ) φ + ( R N k R N + k ) u N ( R N k R N + k ) ( I + τ B ) ( 2 I + τ B ) 1 B 1 i = 1 N 1 ( R N 1 i R N 1 + i ) φ i τ } + ( I + τ B ) ( 2 I + τ B ) 1 B 1 i = 1 N 1 ( R | k i | 1 R k + i 1 ) φ i τ
(7)

for k=0,,N1,

u N = S τ ( j = 1 N ρ ( t j τ 2 ) τ 2 [ { ( I R 2 N ) 1 ( R j R 2 N j + R j 1 R 2 N j + 1 ) φ ( R N j R N + j + R N j + 1 R N + j 1 ) ( I + τ B ) ( 2 I + τ B ) 1 × i = 1 N 1 B 1 ( R N 1 i R N 1 + i ) φ i τ } + ( I + τ B ) ( 2 I + τ B ) 1 B 1 × i = 1 N 1 ( R | j i | 1 R j + i 1 + R | j 1 i | 1 R j + i 2 ) φ i τ ] + ψ )

for k=N.

Proof

{ 1 τ 2 [ u k + 1 2 u k + u k 1 ] + A u k = φ k , 1 k N 1 , N τ = 1 , u 0 = φ , u N  are given
(8)

has a solution and the following formula holds [29]:

u k = ( I R 2 N ) 1 { ( R k R 2 N k ) φ + ( R N k R N + k ) u N ( R N k R N + k ) ( I + τ B ) ( 2 I + τ B ) 1 B 1 i = 1 N 1 ( R N 1 i R N 1 + i ) φ i τ } + ( I + τ B ) ( 2 I + τ B ) 1 B 1 i = 1 N 1 ( R | k i | 1 R k + i 1 ) φ i τ .
(9)

Applying formula (9) and the nonlocal boundary condition

u N = j = 1 N ρ ( t j τ 2 ) ( u j + u j 1 2 ) τ+ψ,

we obtain

u N = ( I j = 1 N ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( R N j R N + j + R N j + 1 R N + j 1 ) ) 1 × ( j = 1 N ρ ( t j τ 2 ) τ 2 [ ( I R 2 N ) 1 { ( R j R 2 N j + R j 1 R 2 N j + 1 ) φ ( R N j R N + j + R N j + 1 R N + j 1 ) × ( I + τ B ) ( 2 I + τ B ) 1 B 1 i = 1 N 1 ( R N 1 i R N 1 + i ) φ i τ } + ( I + τ B ) ( 2 I + τ B ) 1 B 1 × i = 1 N 1 ( R | j i | 1 R j + i 1 + R | j 1 i | 1 R j + i 2 ) φ i τ ] + ψ ) .

Since the operator

I j = 1 N ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( R N j R N + j + R N j + 1 R N + j 1 )

has an inverse S τ , it follows that

u N = S τ ( j = 1 N ρ ( t j τ 2 ) τ 2 [ ( I R 2 N ) 1 { ( R j R 2 N j + R j 1 R 2 N j + 1 ) φ ( R N j R N + j + R N j + 1 R N + j 1 ) × ( I + τ B ) ( 2 I + τ B ) 1 B 1 i = 1 N 1 ( R N 1 i R N 1 + i ) φ i τ } + ( I + τ B ) ( 2 I + τ B ) 1 B 1 × i = 1 N 1 ( R | j i | 1 R j + i 1 + R | j 1 i | 1 R j + i 2 ) φ i τ ] + ψ ) .

Theorem 2 is proved. □

Let F( [ 0 , 1 ] τ ,H) be the linear space of the mesh functions φ τ = { φ k } 1 N 1 with values in the Hilbert space H. We denote by C( [ 0 , 1 ] τ ,H) and C 01 α ( [ 0 , 1 ] τ ,H), 0<α<1, Banach spaces with the norms

φ τ C ( [ 0 , 1 ] τ , H ) = max 1 k N 1 φ k H , φ τ C 01 α ( [ 0 , 1 ] τ , H ) = φ τ C ( [ 0 , 1 ] τ , H ) φ τ C 01 α ( [ 0 , 1 ] τ , H ) = + sup 1 k k + r N 1 ( ( N k ) τ ) α ( ( k + r ) τ ) α ( r τ ) α φ k + r φ k H .

The nonlocal boundary-value problem (2) is said to be stable in F( [ 0 , 1 ] τ ,H) if we have the inequality

u τ F ( [ 0 , 1 ] τ , H ) M(δ) [ φ τ F ( [ 0 , 1 ] τ , H ) + φ H + ψ H ] .

Theorem 3 The solutions of the difference scheme (2) under the assumption (3) satisfy the stability estimate

u τ C ( [ 0 , 1 ] τ , H ) M(δ) [ φ τ C ( [ 0 , 1 ] τ , H ) + φ H + ψ H ] .
(10)

Proof By [29],

u τ C ( [ 0 , 1 ] τ , H ) M(δ) [ φ τ C ( [ 0 , 1 ] τ , H ) + φ H + u N H ]
(11)

is proved for the solution of difference scheme (8). Then the proof of (10) is based on (11) and on the estimate

u N H M(δ) [ φ τ C ( [ 0 , 1 ] τ , H ) + φ H + ψ H ] .

Using the formula (7) and the estimates (4), (5), we get

u N H S τ H H ( j = 1 N | ρ ( t j τ 2 ) | τ 2 [ ( I R 2 N ) 1 H H { ( R j H H + R 2 N j H H + R j 1 H H + R 2 N j + 1 H H ) φ H + ( R N j H H + R N + j H H + R N j + 1 H H + R N + j 1 H H ) ( I + τ B ) ( 2 I + τ B ) 1 H H × B 1 H H i = 1 N 1 τ ( R N i 1 H H + R N + i 1 H H ) φ i H } + ( I + τ B ) ( 2 I + τ B ) 1 H H B 1 H H × ( i = 1 j 1 τ R j i 1 H H φ i H + i = 1 j 1 τ R j i 2 H H φ i H + i = j N 1 τ R i j 2 H H φ i H + i = j N 1 τ R i j 1 H H φ i H + i = 1 N 1 τ R j + i 1 H H φ i H + i = 1 N 1 τ R j + i 2 H H φ i H ) ] + ψ H ) M ( δ ) [ φ τ C ( [ 0 , 1 ] τ , H ) + φ H + ψ H ] .

Theorem 3 is proved. □

Theorem 4 The solutions of the difference problem (2) in C( [ 0 , 1 ] τ ,H) under the assumption (3) obey the almost coercive inequality

{ τ 2 ( u k + 1 2 u k + u k 1 ) } 1 N 1 C ( [ 0 , 1 ] τ , H ) + { A u k } 1 N C ( [ 0 , 1 ] τ , H ) M ( δ ) [ min { ln 1 τ , 1 + | ln B H H | } φ τ C ( [ 0 , 1 ] τ , H ) + A φ H + A ψ H ] .

Proof By [29],

{ τ 2 ( u k + 1 2 u k + u k 1 ) } 1 N 1 C ( [ 0 , 1 ] τ , H ) + { A u k } 1 N C ( [ 0 , 1 ] τ , H ) M ( δ ) [ min { ln 1 τ , 1 + | ln B H H | } φ τ C ( [ 0 , 1 ] τ , H ) + A φ H + A u N H ]

is proved for the solution of the boundary-value problem (8). Using the estimates (4), (5) and the formula (7), we obtain

A u N H M ( δ ) ( min { ln 1 τ , 1 + | ln B H H | } φ τ C ( [ 0 , 1 ] τ , H ) + A φ H + A ψ H )
(12)

for the solution of difference scheme (2). Applying formula (7) and A= B 2 R, we get

A u N = J 1 + J 2 ,

where

J 1 = S τ ( j = 1 N ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( R j R 2 N j + R j 1 R 2 N j + 1 ) A φ + A ψ ) ,
(13)
J 2 = S τ ( j = 1 N ρ ( t j τ 2 ) τ 2 [ ( I R 2 N ) 1 { ( R N j R N + j + R N j + 1 R N + j 1 ) ( I + τ B ) J 2 = × ( 2 I + τ B ) 1 i = 1 N 1 B ( R N i R N + i ) φ i τ } + ( I + τ B ) ( 2 I + τ B ) 1 J 2 = × i = 1 j 1 B ( R j i R j + i + R j i 1 R j + i 1 ) φ i τ + ( I + τ B ) ( 2 I + τ B ) 1 J 2 = × ( i = j N 1 B ( R i j R j + i + R i j 1 R j + i 1 ) φ i τ ) ] ) .
(14)

To this end it suffices to show that

J 1 H M(δ) [ A φ H + A ψ H ]
(15)

and

J 2 H M(δ)min { ln 1 τ , 1 + | ln B H H | } φ τ C ( [ 0 , 1 ] τ , H ) .
(16)

The estimate (15) follows from formula (13) and the estimates (4), (5). Using formula (14) and the estimates (4), (5), we obtain

J 2 H S τ H H ( j = 1 N | ρ ( t j τ 2 ) | τ 2 × [ ( I R 2 N ) 1 H H { ( R N j H H + R N + j H H + R N j + 1 H H + R N + j 1 H H ) ( I + τ B ) ( 2 I + τ B ) 1 H H × i = 1 N 1 ( ( I R ) R N i 1 H H + ( I R ) R N + i 1 H H ) φ i H } + ( I + τ B ) ( 2 I + τ B ) 1 H H × ( i = 1 j ( I R ) R j i 1 H H φ i H + i = 1 j ( I R ) R i + j 1 H H φ i H + i = 1 j ( I R ) R j i 2 H H φ i H + i = 1 j ( I R ) R i + j 2 H H φ i H ) + ( I + τ B ) ( 2 I + τ B ) 1 H H × ( i = j + 1 N 1 ( I R ) R i j 1 H H φ i H + i = j + 1 N 1 ( I R ) R i + j 1 H H φ i H + i = j + 1 N 1 ( I R ) R i j 2 H H φ i H + i = j + 1 N 1 ( I R ) R i + j 2 H H φ i H ) ] ) M ( δ ) min { ln 1 τ , 1 + | ln B H H | } φ τ C ( [ 0 , 1 ] τ , H ) .

From the last estimate and the estimate (15) follows the estimate (12). Theorem 4 is proved. □

Theorem 5 The difference problem (2) is well posed in the Hölder spaces C 01 α ( [ 0 , 1 ] τ ,H) under the assumption (3) and the following coercivity inequality holds:

{ τ 2 ( u k + 1 2 u k + u k 1 ) } 1 N 1 C 01 α ( [ 0 , 1 ] τ , H ) + { A u k } 1 N C 01 α ( [ 0 , 1 ] τ , H ) M ( δ ) [ 1 α ( 1 α ) φ τ C 01 α ( [ 0 , 1 ] τ , H ) + A φ H + A ψ H ] .
(17)

Proof By [29],

{ τ 2 ( u k + 1 2 u k + u k 1 ) } 1 N 1 C 01 α ( [ 0 , 1 ] τ , H ) + { A u k } 1 N 1 C 01 α ( [ 0 , 1 ] τ , H ) M ( δ ) 1 α ( 1 α ) φ τ C 01 α ( [ 0 , 1 ] τ , H ) + M ( δ ) [ A φ H + A u N H ]
(18)

is proved for the solution of difference scheme (8). Then the proof of (17) is based on (18) and on the estimate

A u N H M(δ) 1 α ( 1 α ) φ τ C 01 α ( [ 0 , 1 ] τ , H ) +M(δ) [ A φ H + A ψ H ] .

Applying the triangle inequality, formula (7) and the estimate (15), we get

A u N H J 1 H + J 2 H J 2 H +M(δ) [ A φ H + A ψ H ] .

To this end it suffices to show that

J 2 H M(δ) 1 α ( 1 α ) φ τ C 01 α ( [ 0 , 1 ] τ , H ) .
(19)

Applying formula (14), we get

J 2 = S τ j = 1 N ρ ( t j τ 2 ) τ 2 × ( I R 2 N ) 1 { ( R N j R N + j + R N j + 1 R N + j 1 ) τ 2 ( I R ) 2 × i = 1 j 1 τ 2 ( R N i R N + i ) ( I R 2 ) 1 ( φ i φ j ) + ( ( R N j R N + j + R N j + 1 R N + j 1 ) ) × τ 2 ( I R ) 2 i = j + 1 N 1 τ 2 ( R N i R N + i ) ( I R 2 ) 1 ( φ i φ j ) + ( I R 2 N ) τ 2 ( I R ) 2 i = 1 j 1 τ 2 ( R j i R j + i + R j i 1 R j + i 1 ) × ( I R 2 ) 1 ( φ i φ j ) + ( I R 2 N ) τ 2 ( I R ) 2 × i = j + 1 N 1 τ 2 ( R i j R j + i + R j i 1 R j + i 1 ) ( I R 2 ) 1 ( φ i φ j ) ( R N j R N + j + R N j + 1 R N + j 1 ) τ 2 ( I R ) 2 × i = 1 j 1 τ 2 ( R N i R N + i ) ( I R 2 ) 1 φ j τ 2 ( I R ) 2 × ( R N j R N + j + R N j + 1 R N + j 1 ) i = j + 1 N 1 τ 2 ( R N i R N + i ) ( I R 2 ) 1 φ j + ( I R 2 N ) τ 2 ( I R ) 2 i = 1 j 1 τ 2 ( R j i R j + i + R j i 1 R j + i 1 ) ( I R 2 ) 1 φ j + ( I R 2 N ) τ 2 ( I R ) 2 i = j + 1 N 1 τ 2 ( R i j R j + i + R j i 1 R j + i 1 ) ( I R 2 ) 1 φ j } = z = 2 4 J 2 z ,
(20)

where

J 2 2 = S τ j = 1 N ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( I R ) ( I + R ) 1 J 2 2 = × ( ( I + R ) ( R 2 N j 1 + R 2 N + j 1 R 2 N + j R 2 N j + R j R j 2 ) 2 R 2 N + j 2 ) φ j , J 2 3 = S τ j = 1 N ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( I R ) ( I + R ) R 1 ( I R 2 N 2 j + 1 ) J 2 3 = × i = 1 j 1 R j i ( I R 2 i ) ( I + R ) 1 ( φ i φ j ) = J 2 3 , 1 + J 2 3 , 2 , J 2 3 , 1 = S τ j = 1 [ N 2 ] ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( I R ) ( I + R ) R 1 ( I R 2 N 2 j + 1 ) J 2 3 , 1 = × i = 1 j 1 R j i ( I R 2 i ) ( I + R ) 1 ( φ i φ j ) , J 2 3 , 2 = S τ j = [ N 2 ] + 1 N ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( I R ) ( I + R ) R 1 ( I R 2 N 2 j + 1 ) J 2 3 , 2 = × i = 1 j 1 R j i ( I R 2 i ) ( I + R ) 1 ( φ i φ j ) , J 2 4 = S τ j = 1 N ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( I R ) ( I + R ) ( I R 2 j 1 ) J 2 4 = × i = j + 1 N 1 R i j ( I R 2 N 2 i ) ( I + R ) 1 ( φ i φ j ) = J 2 4 , 1 + J 2 4 , 2 , J 2 4 , 1 = S τ j = 1 [ N 2 ] ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( I R ) ( I + R ) ( I R 2 j 1 ) J 2 4 , 1 = × i = j + 1 N 1 R i j ( I R 2 N 2 i ) ( I + R ) 1 ( φ i φ j ) , J 2 4 , 2 = S τ j = [ N 2 ] + 1 N ρ ( t j τ 2 ) τ 2 ( I R 2 N ) 1 ( I R ) ( I + R ) ( I R 2 j 1 ) J 2 4 , 2 = × i = j + 1 N 1 R i j ( I R 2 N 2 i ) ( I + R ) 1 ( φ i φ j ) .

Second, let us estimate J 2 m for any m=2,,4 separately. We start with J 2 2 . Using estimates (4), (5) and the definition of the norm of the space C 01 α ( [ 0 , 1 ] τ ,H), we get

J 2 2 H S τ H H j = 1 N | ρ ( t j τ 2 ) | τ 2 ( I R 2 N ) 1 H H × I R H H ( I + R ) 1 H H ( I + R H H ( R 2 N j 1 H H + R 2 N + j 1 H H + R 2 N + j H H + R 2 N j H H + R j H H + R j 2 H H ) + R 2 N + j 2 H H ) φ j H M 1 ( δ ) j = 1 N | ρ ( t j τ 2 ) | τ 2 max 1 j N φ j H M 1 ( δ ) j = 1 N | ρ ( t j τ 2 ) | τ 2 φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

From (3) it follows that

J 2 2 H M 2 (δ) φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

Now, let us estimate J 2 3 , 1 . Using estimates (4), (5) and the definition of the norm of the space C 01 α ( [ 0 , 1 ] τ ,H), we obtain

J 2 3 , 1 H S τ H H j = 1 [ N 2 ] | ρ ( t j τ 2 ) | τ 2 ( I R 2 N ) 1 H H × i = 1 j 1 R j i ( I R 2 N 2 j + 1 ) ( I R ) H H I R 2 i H H R 1 H H φ i φ j H M ( δ ) j = 1 [ N 2 ] | ρ ( t j τ 2 ) | τ 2 i = 1 j 1 τ ( ( j i ) τ ) α ( j i ) τ ( ( N i ) τ ) α ( j τ ) α φ τ C 01 α ( [ 0 , 1 ] τ , H ) M ( δ ) j = 1 [ N 2 ] | ρ ( t j τ 2 ) | τ 2 ( j τ ) α ( ( N j ) τ ) α i = 1 j 1 τ ( ( j i ) τ ) 1 α φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

The sum

i = 1 j 1 τ ( ( j i ) τ ) 1 α

is the lower Darboux integral sum for the integral

0 j τ d s ( j τ s ) 1 α .

It follows that

J 2 3 , 1 H M(δ) j = 1 [ N 2 ] | ρ ( t j τ 2 ) | τ α ( ( N j ) τ ) α φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

By the lower Darboux integral sum for the integral it concludes that

J 2 3 , 1 H M(δ) 2 α 2 α ( 1 α ) j = 1 [ N 2 ] |ρ ( t j τ 2 ) |τ φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

For J 2 3 , 2 , applying (4), (5) and the definition of the norm of the space C 01 α ( [ 0 , 1 ] τ ,H), we get

J 2 3 , 2 H S τ H H j = [ N 2 ] + 1 N | ρ ( t j τ 2 ) | τ 2 ( I R 2 N ) 1 H H × i = 1 j 1 R j i ( I R 2 N 2 j + 1 ) ( I R ) H H I R 2 i H H R 1 H H φ i φ j H M ( δ ) j = [ N 2 ] + 1 N | ρ ( t j τ 2 ) | τ 2 × i = 1 j 1 ( 2 τ ( N j + 1 ) ) α ( ( j i ) τ ) 1 α ( ( 2 N j i + 1 ) τ ) α ( j τ ) α ( ( N i ) τ ) α φ τ C 01 α ( [ 0 , 1 ] τ , H ) M ( δ ) j = [ N 2 ] + 1 N | ρ ( t j τ 2 ) | τ 2 α ( ( N j + 1 ) τ ) α ( ( N j ) τ ) α ( j τ ) α × i = 1 j 1 τ ( ( j i ) τ ) 1 α ( ( N j i + N + 1 ) τ ) α φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

The sum

i = 1 j 1 τ ( ( j i ) τ ) 1 α ( ( N j i + N + 1 ) τ ) α

is the lower Darboux integral sum for the integral

0 j τ d s ( j τ s ) 1 α ( N τ j τ s + τ + N τ ) α .

Since

0 j τ d s ( j τ s ) 1 α ( N τ j τ s + N τ + τ ) α 1 ( N τ j τ + τ ) α 0 j τ d s ( j τ s ) 1 α M α ( j τ ) α ,

it follows that

J 2 3 , 2 H M(δ) j = [ N 2 ] + 1 N |ρ ( t j τ 2 ) |τ 2 α ( j τ ) α ( N τ j τ + τ ) α α ( j τ ) α φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

By the lower Darboux integral sum for the integral it follows that

J 2 3 , 2 H M ( δ ) 2 2 α 2 α ( 1 α ) j = [ N 2 ] + 1 N |ρ ( t j τ 2 ) |τ φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

Combining J 2 3 , 1 and J 2 3 , 2 , we get

J 2 3 H M 3 ( δ ) α ( 1 α ) j = 1 N |ρ ( t j τ 2 ) |τ φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

From (3) it follows that

J 2 3 H M 4 ( δ ) α ( 1 α ) φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

Next, let us estimate J 2 4 , 1 . Using the estimates (4), (5), and the definition of the norm space C 01 α ( [ 0 , 1 ] τ ,H), we obtain

J 2 4 , 1 H S τ H H j = 1 [ N 2 ] | ρ ( t j τ 2 ) | τ 2 ( I R 2 N ) 1 H H I R 2 j 1 H H × i = j + 1 N 1 R i j ( I R 2 N 2 i ) ( I R ) H H φ i φ j H M ( δ ) j = 1 [ N 2 ] | ρ ( t j τ 2 ) | τ 2 × i = j + 1 N 1 2 α ( ( i j ) τ ) α ( N τ i τ ) α ( ( i j ) τ ) ( ( N j ) τ ) α ( i τ ) α ( ( 2 N j i ) τ ) α φ τ C 01 α ( [ 0 , 1 ] τ , H ) M ( δ ) j = 1 [ N 2 ] | ρ ( t j τ 2 ) | τ ( N τ j τ ) α 2 ( N τ j τ ) α ( j τ ) α × i = j + 1 N 1 τ ( ( i j ) τ ) 1 α ( ( N j i + N ) τ ) α φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

The sum

i = j + 1 N 1 τ ( ( i j ) τ ) 1 α

is the lower Darboux integral sum for the integral

j τ 1 d s ( s j τ ) 1 α .

Since

j τ 1 d s ( s j τ ) 1 α ( 2 N τ j τ s ) α 1 ( N τ j τ ) α j τ 1 d s ( s j τ ) 1 α ( N τ j τ ) α α ( N τ j τ ) α ,

it follows that

J 2 4 , 1 H M(δ) j = 1 [ N 2 ] | ρ ( t j τ 2 ) | τ ( j τ ) α α φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

By the lower Darboux integral sum for the integral it follows that

J 2 4 , 1 H M ( δ ) 2 2 α 1 α ( 1 α ) j = 1 [ N 2 ] |ρ ( t j τ 2 ) |τ φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

Finally, let us estimate J 2 4 , 2 . Using the estimates (4), (5), and the definition of the norm space C 01 α ( [ 0 , 1 ] τ ,H), we get

J 2 4 , 2 H S τ H H j = [ N 2 ] + 1 N | ρ ( t j τ 2 ) | τ 2 ( I R 2 N ) 1 H H I R 2 j 1 H H × i = j + 1 N 1 R i j ( I R 2 N 2 i ) ( I R ) H H φ i φ j H M ( δ ) j = [ N 2 ] + 1 N | ρ ( t j τ 2 ) | τ 2 × i = j + 1 N 1 τ ( ( i j ) τ ) α ( ( i j ) τ ) ( ( N i ) τ ) α ( i τ ) α φ τ C 01 α ( [ 0 , 1 ] τ , H ) M ( δ ) j = [ N 2 ] + 1 N | ρ ( t j τ 2 ) | τ 2 ( ( N j ) τ ) α ( j τ ) α i = j + 1 N 1 τ ( ( i j ) τ ) 1 α φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

The sum

i = j + 1 N 1 τ ( ( i j ) τ ) 1 α

is the lower Darboux integral sum for the integral

j τ 1 d s ( s j τ ) 1 α .

Thus, we show that

J 2 4 , 2 H M(δ) j = [ N 2 ] + 1 N | ρ ( t j τ 2 ) | τ ( j τ ) α α φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

By the lower Darboux integral sum for the integral it follows that

J 2 4 , 2 H M(δ) 2 α 2 α ( 1 α ) j = [ N 2 ] + 1 N |ρ ( t j τ 2 ) |τ φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

Combining J 2 4 , 1 and J 2 4 , 2 , we get

J 2 4 H M(δ) ( 2 α 2 α ( 1 α ) + 2 2 α 1 α ( 1 α ) ) j = 1 N |ρ ( t j τ 2 ) |τ φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

From (3) it follows that

J 2 4 H M 5 ( δ ) α ( 1 α ) φ τ C 01 α ( [ 0 , 1 ] τ , H ) .

Combining estimates for J 2 m , m=2,,4 we get the estimate (19). Theorem 5 is proved. □

3 Applications

Now, the application of Theorems 3-5 will be given. First, we consider the mixed boundary-value problem for elliptic equation

{ u t t ( a ( x ) u x ) x + δ u = f ( t , x ) , 0 < t < 1 , 0 < x < 1 , u ( t , 0 ) = u ( t , 1 ) , u x ( t , 1 ) = u x ( t , 0 ) , 0 t 1 , u ( 0 , x ) = φ ( x ) , u ( 1 , x ) = 0 1 ρ ( λ ) u ( λ , x ) d λ + ψ ( x ) , 0 x 1 ,
(21)

where a(x), φ(x), ψ(x) and f(t,x) are given sufficiently smooth functions and a(x)a>0, a(1)=a(0), δ=const>0. The discretization of problem (21) is carried out in two steps. In the first step, let us define the grid space

[ 0 , 1 ] h ={x: x n =nh,0nM,Mh=1}.

We introduce the Hilbert space L 2 h = L 2 ( [ 0 , 1 ] h ) of the grid functions φ h (x)= { φ n } n = 1 M 1 defined on [ 0 , 1 ] h , equipped with the norms

φ h L 2 h = ( x [ 0 , 1 ] h | φ h ( x ) | 2 h ) 1 2 , φ h W 2 h 2 = φ h L 2 h + ( x [ 0 , 1 ] h | ( φ h ( x ) ) x | 2 h ) 1 / 2 + ( x [ 0 , 1 ] h | ( φ h ( x ) ) x ¯ , x | 2 h ) 1 / 2 .

To the differential operator A generated by the problem (21) we assign the difference operator A h x by the formula

A h x φ h (x)= { ( a ( x ) φ x ¯ ) x , n + δ φ n } 1 M 1 ,
(22)

acting in the space of the grid functions φ h (x)= { φ n } 1 M 1 satisfying the conditions φ 0 = φ M , φ 1 φ 0 = φ M φ M 1 . It is know that A h x is a self-adjoint positive definite operator in L 2 h . With the help of A h x , we arrive at the nonlocal boundary-value problem

{ d 2 u h ( t , x ) d t 2 + A h x u h ( t , x ) = f h ( t , x ) , 0 < t < 1 , x [ 0 , 1 ] h , u h ( 0 , x ) = φ h ( x ) ; u h ( 1 , x ) = 0 1 ρ ( t ) u h ( t , x ) d t + ψ h ( x ) , x [ 0 , 1 ] h
(23)

for an infinite system of ordinary differential equations. Therefore, in the second step, equation (23) is replaced by the difference scheme (2), and we get the following difference scheme:

{ u k + 1 h ( x ) 2 u k h ( x ) + u k 1 h ( x ) τ 2 + A h x u k h ( x ) = φ k h ( x ) , φ k h ( x ) = f h ( t k , x ) , t k = k τ , 1 k N 1 , N τ = 1 , x [ 0 , 1 ] h , u 0 h ( x ) = φ h ( x ) ; u N h ( x ) = j = 1 N ρ ( t j τ 2 ) τ ( u j h ( x ) + u j 1 h ( x ) 2 ) + ψ h ( x ) , x [ 0 , 1 ] h
(24)

for the numerical solution of (21).

Theorem 6 Let τ and h be sufficiently small positive numbers. Then under the assumption (3), the solution of the difference scheme (24) satisfies the following stability and almost coercivity estimates:

max 1 k N 1 u k h L 2 h M ( δ ) [ max 1 k N 1 φ k h L 2 h + ψ h L 2 h + φ h L 2 h ] , max 1 k N 1 τ 2 ( u k + 1 h 2 u k h + u k 1 h ) L 2 h + max 1 k N 1 ( u k h ) W 2 h 2 M ( δ ) [ ln 1 τ + | h | max 1 k N 1 φ k h L 2 h + φ h W 2 h 2 + ψ h W 2 h 2 ] .

The proof of Theorem 6 is based on Theorems 3 and 4, on the estimate

min { ln 1 τ , 1 + | ln B h x L 2 h L 2 h | } Mln 1 τ + | h | ,
(25)

and on the symmetry properties of the difference operator A h x defined by the formula (22) in L 2 h .

Theorem 7 Let τ and |h| be sufficiently small positive numbers. Then under the assumption (3), the solution of the difference scheme (24) satisfies the following coercivity estimate:

{ τ 2 ( u k + 1 h 2 u k h + u k 1 h ) } 1 N 1 C 01 α ( [ 0 , 1 ] τ , L 2 h ) + { u k h } 1 N 1 C 01 α ( [ 0 , 1 ] τ , W 2 h 2 ) M ( δ ) [ φ h W 2 h 2 + ψ h W 2 h 2 + 1 α ( 1 α ) { φ k h } 1 N 1 C 01 α ( [ 0 , 1 ] τ , L 2 h ) ] .

The proof of Theorem 7 is based on Theorem 5 and the symmetry properties of the difference operator A h x defined by formula (22).

Second, let Ω be the unit open cube in R n (x=( x 1 ,, x n ):0< x k <1,1kn) with boundary S, Ω ¯ =ΩS. In [0,1]×Ω, the Dirichlet-Bitsadze-Samarskii type mixed boundary-value problem for the multidimensional elliptic equation

{ u t t r = 1 n ( a r ( x ) u x r ) x r = f ( t , x ) , 0 < t < 1 , x = ( x 1 , , x n ) Ω , u ( 0 , x ) = φ ( x ) , u ( 1 , x ) = 0 1 ρ ( λ ) u ( λ , x ) d λ + ψ ( x ) , x Ω ¯ , u ( t , x ) | x S = 0 , x Ω ¯
(26)

is considered. We will study the problem (26) under the assumption (3). Here, a r (x) (xΩ), ψ(x), φ(x) (x Ω ¯ ) and f(t,x) (t(0,1), xΩ) are smooth functions and a r (x)a>0. The discretization of problem (26) is carried out in two steps. In the first step let us define the grid sets

Ω ¯ h = { x = x m = ( h 1 m 1 , , h m m m ) , m = ( m 1 , , m m ) , 0 m r N r , h r N r = 1 , r = 1 , , m } , Ω h = Ω ˜ h Ω , S h = Ω ˜ h S .

We introduce the Hilbert space L 2 h = L 2 ( Ω ˜ h ) of the grid functions φ h (x)={φ( h 1 m 1 ,, h m m m )} defined on Ω ˜ h , equipped with the norms

φ h L 2 h = ( x Ω ¯ h | φ h ( x ) | 2 h 1 h m ) 1 / 2 , φ h W 2 h 2 = φ h L 2 h + ( x Ω ˜ h r = 1 m | ( φ h ( x ) ) x r | 2 h 1 h m ) 1 / 2 φ h W 2 h 2 = + ( x Ω ˜ h r = 1 m | ( φ h ( x ) ) x r x ¯ r , m r | 2 h 1 h m ) 1 / 2 .

To the differential operator A generated by the problem (26), we assign the difference operator A h x by the formula

A h x u h (x)= r = 1 m ( a r ( x ) u x ¯ r h ) x r , j r ,
(27)

acting in the space of the grid functions u h (x), satisfying the conditions u h (x)=0 for all x S h . It is known that A h x is a self-adjoint positive definite operator in L 2 h . With the help of A h x , we arrive at the nonlocal boundary-value problem for an infinite system of ordinary differential equations

{ d 2 u h ( t , x ) d t 2 + A h x u h ( t , x ) = f h ( t , x ) , 0 < t < 1 , x Ω ˜ h , u h ( 0 , x ) = φ h ( x ) ; u h ( 1 , x ) = 0 1 ρ ( t ) u h ( t , x ) d t + ψ h ( x ) , x Ω ˜ h .
(28)

In the second step, (28) is replaced by the difference scheme (2), and we get the following difference scheme:

{ u k + 1 h ( x ) 2 u k h ( x ) + u k 1 h ( x ) τ 2 + A h x u k h ( x ) = φ k h ( x ) , φ k h ( x ) = f h ( t k , x ) , x Ω h , t k = k τ , 1 k N 1 , N τ = 1 , u 0 h ( x ) = φ h ( x ) , x Ω ˜ h , u N h ( x ) = j = 1 N ρ ( t j τ 2 ) τ ( u j h ( x ) + u j 1 h ( x ) 2 ) + ψ h ( x ) , x Ω ˜ h
(29)

for the numerical solution of (26).

Theorem 8 Let τ and |h|= h 1 2 + + h n 2 be sufficiently small positive numbers. Then under the assumption (3) the solution of the difference scheme (29) satisfies the following stability estimates:

max 1 k N 1 u k h L 2 h M(δ) [ max 1 k N 1 φ k h L 2 h + ψ h L 2 h + φ h L 2 h ] .

The proof of Theorem 8 is based on Theorem 3 and the symmetry properties of the difference operator A h x defined by (27) in L 2 h .

Theorem 9 Let τ and |h| be sufficiently small positive numbers. Then under the assumption (3) the solution of the difference scheme (29) satisfies the following almost coercivity estimates:

max 1 k N 1 τ 2 ( u k + 1 h 2 u k h + u k 1 h ) L 2 h + max 1 k N 1 u k h W 2 h 2 M ( δ ) [ ln 1 τ + | h | max 1 k N 1 φ k h L 2 h + φ h W 2 h 2 + ψ h W 2 h 2 ] .

The proof of Theorem 9 is based on Theorem 4, on the estimate (25), on the symmetry properties of the difference operator A h x defined by (27) in L 2 h , and on the following theorem on the coercivity inequality for the solution of the elliptic difference problem in  L 2 h .

Theorem 10 For the solutions of the elliptic difference problem

A h x u h (x)= ω h (x),x Ω h , u h (x)=0,x S h

the following coercivity inequality holds [30]:

u h W 2 h 2 M(δ) ω h L 2 h .

Theorem 11 Let τ and |h| be sufficiently small positive numbers. Then under the assumption (3) the solution of the difference scheme (29) satisfies the following coercivity stability estimate:

{ τ 2 ( u k + 1 h 2 u k h + u k 1 h ) } 1 N 1 C 01 α ( [ 0 , 1 ] τ , L 2 h ) + { u k h } 1 N 1 C 01 α ( [ 0 , 1 ] τ , W 2 h 2 ) M ( δ ) [ φ h W 2 h 2 + ψ h W 2 h 2 + 1 α ( 1 α ) { φ k h } 1 N 1 C 01 α ( [ 0 , 1 ] τ L 2 h ) ] .

The proof of Theorem 11 is based on Theorem 5, on the symmetry properties of the difference operator A h x defined by the formula (27), and on Theorem 10 on the coercivity inequality for the solution of the elliptic difference equation in L 2 h .

4 Numerical results

We consider the Bitsadze -Samarskii type nonlocal boundary problem for the elliptic equation

{ 2 u ( t , x ) t 2 2 u ( t , x ) x 2 + u = π 2 exp ( t ) sin ( π x ) , 0 < t < 1 , 0 < x < 1 , u ( 0 , x ) = sin ( π x ) , u ( 1 , x ) = 0 1 e λ u ( λ , x ) d λ + ( exp ( 1 ) + 1 2 exp ( 2 ) 1 2 ) sin ( π x ) , 0 < x < 1 , u ( t , 0 ) = u ( t , 1 ) = 0 , 0 < t < 1 .
(30)

The exact solution of this problem is

u(t,x)=exp(t)sin(πx).

In the present part for the approximate solutions of the Bitsadze-Samarskii type nonlocal boundary-value problem (30), we will use the first and second orders of the accuracy difference schemes with grid intervals τ= 1 N , h= 1 M for t and x, respectively. For the approximate solution of the nonlocal boundary Bitsadze-Samarskii type problem (30), we consider the set [ 0 , 1 ] τ × [ 0 , 1 ] h of a family of grid points depending on the small parameters τ and h,

[ 0 , 1 ] τ × [ 0 , 1 ] h = { ( t k , x n ) : t k = k τ , 1 k N 1 , N τ = 1 , x n = n h , 1 n M 1 , M h = 1 } .

Applying the first order of the accuracy difference scheme from [31] and the second order of the accuracy difference scheme (2) for the approximate solutions of the problem, we have the second-order difference equations with respect to n with matrix coefficients. To solve these difference equations, we have applied the procedure of a modified Gauss elimination method for the difference equations with respect to n with matrix coefficients. To obtain the solution of (2), we use MATLAB programming. The errors are computed by

E M N = max 1 k N 1 ( n = 1 M 1 | u ( t k , x n ) u n k | 2 h ) 1 2

of numerical solutions for different values of M and N, where u( t k , x n ) represents the exact solution and u n k represents the numerical solution at ( t k , x n ). The results are shown in Table 1, respectively.

Table 1 The errors for first- and second-order difference scheme

5 Conclusion

In this paper, the second order of the accuracy difference scheme for the approximate solution of the Bitsadze-Samarskii type nonlocal boundary-value problem with the integral condition for elliptic equations is presented. Theorems on the stability estimates, almost coercive stability estimates, and coercive stability estimates for the solution of difference scheme for elliptic equations are proved. The theoretical statements for the solution of this difference scheme are supported by the result of a numerical example. As can be seen from Table 1, the second order of the accuracy difference scheme is more accurate than the first order of the accuracy difference scheme.

References

  1. Bitsadze AV, Samarskii AA: On some simplest generalizations of linear elliptic problems. Dokl. Akad. Nauk SSSR 1969, 185(4):739-740.

    MathSciNet  Google Scholar 

  2. Samarskii AA: Some problems in differential equation theory. Differ. Uravn. 1980, 16(11):1925-1935.

    MathSciNet  Google Scholar 

  3. Sapagovas MP: Difference method of increased order of accuracy for the Poisson equation with nonlocal conditions. Differ. Equ. 2008, 44(7):1018-1028. 10.1134/S0012266108070148

    Article  MATH  MathSciNet  Google Scholar 

  4. Ashyralyev A: A note on the Bitsadze-Samarskii type nonlocal boundary value problem in a Banach space. J. Math. Anal. Appl. 2008, 344(1):557-573. 10.1016/j.jmaa.2008.03.008

    Article  MATH  MathSciNet  Google Scholar 

  5. Ashyralyev A, Ozturk E: The numerical solution of Bitsadze-Samarskii nonlocal boundary value problems with the Dirichlet-Neumann condition. Abstr. Appl. Anal. 2012., 2012: Article ID 730804

    Google Scholar 

  6. Ashyralyev A, Ozturk E: On a difference scheme of fourth-order of accuracy for the Bitsadze-Samarskii type nonlocal boundary value problem. Math. Methods Appl. Sci. 2013, 36(8):936-955. 10.1002/mma.2650

    Article  MATH  MathSciNet  Google Scholar 

  7. Ashyralyev A, Ozturk E: On Bitsadze-Samarskii type nonlocal boundary value problems for elliptic differential and difference equations: well-posedness. Appl. Math. Comput. 2012, 219(3):1093-1107. 10.1016/j.amc.2012.07.016

    Article  MATH  MathSciNet  Google Scholar 

  8. Volkova EA, Dosiyev AA, Buranay SC: On the solution of a nonlocal problem. Comput. Math. Appl. 2013, 66(3):330-338. 10.1016/j.camwa.2013.05.010

    Article  MathSciNet  Google Scholar 

  9. Ashyralyev A, Tetikoğlu FS: A note on Bitsadze-Samarskii type nonlocal boundary value problems: well-posedness. Numer. Funct. Anal. Optim. 2013, 34(9):939-975. 10.1080/01630563.2012.738458

    Article  MATH  MathSciNet  Google Scholar 

  10. Berikelashvili G: On a nonlocal boundary value problem for a two-dimensional elliptic equation. Comput. Methods Appl. Math. 2003, 3(1):35-44.

    MATH  MathSciNet  Google Scholar 

  11. Gordeziani DG: On a method of resolution of Bitsadze-Samarskii boundary value problem. Abst. Rep. Inst. Appl. Math. Tbilisi State Univ. 1970, 2: 38-40.

    Google Scholar 

  12. Kapanadze DV: On the Bitsadze-Samarskii nonlocal boundary value problem. Differ. Equ. 1987, 23(3):543-545.

    MATH  MathSciNet  Google Scholar 

  13. Il’in VA, Moiseev EI: Two-dimensional nonlocal boundary value problems for Poisson’s operator in differential and difference variants. Mat. Model. 1990, 1(2):139-159.

    MathSciNet  Google Scholar 

  14. Berikelashvili GK: On the convergence rate of the finite-difference solution of a nonlocal boundary value problem for a second-order elliptic equation. Differ. Equ. 2003, 39(7):945-953.

    Article  MATH  MathSciNet  Google Scholar 

  15. Skubaczewski AL: Solvability of elliptic problems with Bitsadze-Samarskii boundary conditions. Differ. Uravn. 1985, 21(4):701-706.

    Google Scholar 

  16. Ashyralyev A, Tetikoğlu FS: FDM for elliptic equations with Bitsadze-Samarskii-Dirichlet conditions. Abstr. Appl. Anal. 2012., 2012: Article ID 454831

    Google Scholar 

  17. Ladyzhenskaya OA, Ural’tseva NN: Linear and Quasilinear Equations of Elliptic Type. Nauka, Moscow; 1973. (Russian)

    MATH  Google Scholar 

  18. Gorbachuk VL, Gorbachuk ML: Boundary Value Problems for Differential-Operator Equations. Naukova Dumka, Kiev; 1984. (Russian)

    Google Scholar 

  19. Grisvard P: Elliptic Problems in Nonsmooth Domains. Pitman, London; 1985.

    MATH  Google Scholar 

  20. Agmon S: Lectures on Elliptic Boundary Value Problems. Van Nostrand, Princeton; 1965.

    MATH  Google Scholar 

  21. Krein SG: Linear Differential Equations in Banach Space. Nauka, Moscow; 1966. (Russian)

    Google Scholar 

  22. Skubachevskii AL Operator Theory - Advances and Applications. In Elliptic Functional Differential Equations and Applications. Birkhäuser, Basel; 1997.

    Google Scholar 

  23. Agarwal R, Bohner M, Shakhmurov VB: Maximal regular boundary value problems in Banach-valued weighted spaces. Bound. Value Probl. 2005, 18: 9-42.

    MathSciNet  Google Scholar 

  24. Sobolevskii PE: Well-posedness of difference elliptic equations. Discrete Dyn. Nat. Soc. 1997, 1(4):219-231.

    Article  MATH  Google Scholar 

  25. Ashyralyev A: Well-posed solvability of the boundary value problem for difference equations of elliptic type. Nonlinear Anal., Theory Methods Appl. 1995, 24(2):251-256. 10.1016/0362-546X(94)E0003-Y

    Article  MATH  MathSciNet  Google Scholar 

  26. Ashyralyev A, Altay N: A note on the well-posedness of the nonlocal boundary value problem for elliptic difference equations. Appl. Math. Comput. 2006, 175(1):49-60. 10.1016/j.amc.2005.07.013

    Article  MATH  MathSciNet  Google Scholar 

  27. Agmon S, Douglis SA, Nirenberg L: Estimates near the boundary for solutions of elliptic partial differential equations satisfying general boundary conditions II. Commun. Pure Appl. Math. 1964, 17: 35-92. 10.1002/cpa.3160170104

    Article  MATH  MathSciNet  Google Scholar 

  28. Sobolevskii PE: On elliptic equations in a Banach space. Differ. Uravn. 1969, 4(7):1346-1348. (Russian)

    MathSciNet  Google Scholar 

  29. Ashyralyev A, Sobolevskii PE: New Difference Schemes for Partial Differential Equations. Birkhäuser, Basel; 2004.

    Book  MATH  Google Scholar 

  30. Sobolevskii PE: On Difference Method for Approximate Solution of Differential Equations. Izdat. Voronezh. Gosud. Univ., Voronezh; 1975.

    Google Scholar 

  31. Ozturk, E: Nonlocal boundary value problems for elliptic differential and difference equations. PhD thesis, Mathematics Department, Uludag University (2013) (Turkish)

Download references

Acknowledgements

The authors would like to thank Prof. Dr. PE Sobolevskii for his helpful suggestions on the improvement of this paper.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Elif Ozturk.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

EO carried out the studies, participated in the sequence alignment and drafted the manuscript and AA carried out the studies, participated in the sequence alignment. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Ashyralyev, A., Ozturk, E. On a difference scheme of second order of accuracy for the Bitsadze-Samarskii type nonlocal boundary-value problem. Bound Value Probl 2014, 14 (2014). https://doi.org/10.1186/1687-2770-2014-14

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2014-14

Keywords