Research

# High order of accuracy difference schemes for the inverse elliptic problem with Dirichlet condition

Charyyar Ashyralyyev

Author Affiliations

Department of Mathematical Engineering, Gumushane University, Gumushane, 29100, Turkey

TAU, 2009 Street, 143, Ashgabat, 744000, Turkmenistan

Boundary Value Problems 2014, 2014:5  doi:10.1186/1687-2770-2014-5

 Received: 21 October 2013 Accepted: 11 December 2013 Published: 7 January 2014

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

### Abstract

The overdetermination problem for elliptic differential equation with Dirichlet boundary condition is considered. The third and fourth orders of accuracy stable difference schemes for the solution of this inverse problem are presented. Stability, almost coercive stability, and coercive inequalities for the solutions of difference problems are established. As a result of the application of established abstract theorems, we get well-posedness of high order difference schemes of the inverse problem for a multidimensional elliptic equation. The theoretical statements are supported by a numerical example.

MSC: 35N25, 39A14, 39A30, 65J22.

##### Keywords:
difference scheme; inverse elliptic problem; high order accuracy; well-posedness; stability; almost coercive stability; coercive stability

### 1 Introduction

Many problems in various branches of science lead to inverse problems for partial differential equations [1-3]. Inverse problems for partial differential equations have been investigated extensively by many researchers (see [3-18] and the references therein).

Consider the inverse problem of finding a function u and an element p for the elliptic equation

{ u t t ( t ) + A u ( t ) = f ( t ) + p , 0 < t < T , u ( 0 ) = φ , u ( T ) = ψ , u ( λ ) = ξ , 0 < λ < T (1.1)

in an arbitrary Hilbert space H with a self-adjoint positive definite operator A. Here, λ is a known number, φ, ξ, and ψ are given elements of H.

Existence and uniqueness theorems for problem (1.1) in a Banach space are presented in [5]. The first and second accuracy stable difference schemes for this problem have been constructed in [15]. High order of accuracy stable difference schemes for nonlocal boundary value elliptic problems are presented in [19-21].

Our aim in this work is the construction of the third and fourth order stable accuracy difference schemes for the inverse problem (1.1).

In the present paper, the third and fourth orders of accuracy difference schemes for the approximate solution of problem (1.1) are presented. Stability, almost coercive stability, and coercive stability inequalities for the solution of these difference schemes are established.

In the application, we consider the inverse problem for the multidimensional elliptic equation with Dirichlet condition

{ u t t ( t , x ) r = 1 n ( a r ( x ) u x r ) x r + σ u = f ( t , x ) + p ( x ) , x = ( x 1 , , x n ) Ω , 0 < t < T , u ( 0 , x ) = φ ( x ) , u ( T , x ) = ψ ( x ) , u ( λ , x ) = ξ ( x ) , x Ω ¯ , u ( t , x ) = 0 , x S , 0 t T . (1.2)

Here, Ω = ( 0 , L ) × × ( 0 , L ) is the open cube in the n-dimensional Euclidean space with boundary S, Ω ¯ = Ω S , a r ( x ) ( x Ω ), φ ( x ) , ξ ( x ) , ψ ( x ) ( x Ω ¯ ), f ( t , x ) ( t ( 0 , 1 ) , x Ω ) are given smooth functions, a r ( x ) a > 0 ( x Ω ¯ ), and 0 < λ < T , σ > 0 are given numbers.

The first and second orders of accuracy stable difference schemes for equation (1.2) are presented in [15]. We construct the third and fourth orders of accuracy stable difference schemes for problem (1.2).

The remainder of this paper is organized as follows. In Section 2, we present the third and fourth order difference schemes for problem (1.1) and obtain stability estimates for them. In Section 3, we construct the third and fourth order difference schemes for problem (1.2) and establish their well-posedness. In Section 4, the numerical results are given. Section 5 is our conclusion.

### 2 High order of accuracy difference schemes for (1.1) and stability inequalities

We use, respectively, the third and fourth order accuracy approximate formulas

u ( λ ) = ( 1 2 ( λ τ [ λ τ ] ) + 1 2 ( λ τ [ λ τ ] ) 2 ) u ( ( [ λ τ ] 1 ) τ ) u ( λ ) = + ( 1 ( λ τ [ λ τ ] ) 2 ) u ( [ λ τ ] τ ) u ( λ ) = + ( 1 2 ( λ τ [ λ τ ] ) + 1 2 ( λ τ [ λ τ ] ) 2 ) u ( ( [ λ τ ] + 1 ) τ ) + o ( τ 3 ) , (2.1)

u ( λ ) = ( 1 12 ( λ τ [ λ τ ] ) + 1 12 ( λ τ [ λ τ ] ) 3 ) u ( ( [ λ τ ] 2 ) τ ) u ( λ ) = + ( 8 12 ( λ τ [ λ τ ] ) + 1 2 ( λ τ [ λ τ ] ) 2 u ( λ ) = 1 6 ( λ τ [ λ τ ] ) 3 ) u ( ( [ λ τ ] 1 ) τ ) u ( λ ) = + ( 1 ( λ τ [ λ τ ] ) 2 ) u ( [ λ τ ] τ ) + ( 8 12 ( λ τ [ λ τ ] ) + 1 2 ( λ τ [ λ τ ] ) 2 u ( λ ) = + 1 6 ( λ τ [ λ τ ] ) 3 ) u ( ( [ λ τ ] 1 ) τ ) u ( λ ) = + ( 1 12 ( λ τ [ λ τ ] ) 1 12 ( λ τ [ λ τ ] ) 3 ) u ( ( [ λ τ ] + 2 ) τ ) + o ( τ 4 ) (2.2)

for u ( λ ) . Here, l = [ λ τ ] , [ ] is a notation for the greatest integer function. Applying formulas (2.1) and (2.2) to u ( λ ) = ξ , we get, respectively,

{ τ 2 ( u k + 1 2 u k + u k 1 ) + A u k + τ 2 12 A 2 u k = θ k + p , θ k = f ( t k ) + τ 2 12 ( f ( t k + 1 ) 2 f ( t k ) + f ( t k 1 ) τ 2 + A f ( t k ) ) , t k = k τ , 1 k N 1 , N τ = T , u 0 = φ , u N = ψ , ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) u l 1 + ( 1 1 2 ( λ τ l ) 2 ) u l + ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) u l + 1 = ξ , (2.3)

the third order of accuracy difference problem and

{ τ 2 ( u k + 1 2 u k + u k 1 ) + A u k + τ 2 12 A 2 u k = θ k + p , θ k = f ( t k ) + τ 2 12 ( f ( t k + 1 ) 2 f ( t k ) + f ( t k 1 ) τ 2 + A f ( t k ) ) , t k = k τ , 1 k N 1 , N τ = T , u 0 = φ , u N = ψ , ( 1 12 ( λ τ l ) 1 12 ( λ τ l ) 3 ) u l 2 + ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 + 1 6 ( λ τ l ) 3 ) u l 1 + ( 1 ( λ τ l ) 2 ) u l + ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 1 6 ( λ τ l ) 3 ) u l + 1 + ( 1 12 ( λ τ l ) + 1 12 ( λ τ l ) 3 ) u l + 2 = ξ , (2.4)

the fourth order of accuracy difference problem for inverse problem (1.1).

For solving of problems (2.3) and (2.4), we use the algorithm [14], which includes three stages. For finding a solution { u k } k = 1 N 1 of difference problems (2.3) and (2.4) we apply the substitution

u k = v k + A 1 p . (2.5)

In the first stage, applying approximation (2.5), we get a nonlocal boundary value difference problem for obtaining { v k } k = 0 N . In the second stage, we put k = 0 and find v 0 . Then, using the formula

p = A φ A v 0 , (2.6)

we define an element p. In the third stage, by using approximation (2.5), we can obtain the solution { u k } k = 1 N 1 of difference problems (2.3) and (2.4). In the framework of the above mentioned algorithm for { v k } k = 0 N , we get the following auxiliary nonlocal boundary value difference scheme:

{ τ 2 ( v k + 1 2 v k + v k 1 ) + A v k + τ 2 12 A 2 v k = θ k , θ k = f ( t k ) + τ 2 12 ( f ( t k + 1 ) 2 f ( t k ) + f ( t k 1 ) τ 2 + A f ( t k ) ) , t k = k τ , 1 k N 1 , N τ = T , v 0 ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v l 1 ( 1 1 2 ( λ τ l ) 2 ) v l ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v l + 1 = φ ξ , v N ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v l 1 ( 1 1 2 ( λ τ l ) 2 ) v l ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v l + 1 = ψ ξ (2.7)

for the third order of accuracy difference problem (2.3) and

{ τ 2 ( v k + 1 2 v k + v k 1 ) + A v k + τ 2 12 A 2 v k = θ k , θ k = f ( t k ) + τ 2 12 ( f ( t k + 1 ) 2 f ( t k ) + f ( t k 1 ) τ 2 + A f ( t k ) ) , t k = k τ , 1 k N 1 , N τ = T , v 0 ( 1 12 ( λ τ l ) 1 12 ( λ τ l ) 3 ) v l 2 ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 + 1 6 ( λ τ l ) 3 ) v l 1 ( 1 ( λ τ l ) 2 ) v l ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 1 6 ( λ τ l ) 3 ) v l + 1 ( 1 12 ( λ τ l ) + 1 12 ( λ τ l ) 3 ) v l + 2 = φ ξ , v N ( 1 12 ( λ τ l ) 1 12 ( λ τ l ) 3 ) v l 2 ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 + 1 6 ( λ τ l ) 3 ) v l 1 ( 1 ( λ τ l ) 2 ) v l ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 1 6 ( λ τ l ) 3 ) v l + 1 ( 1 12 ( λ τ l ) + 1 12 ( λ τ l ) 3 ) v l + 2 = ψ ξ (2.8)

for the fourth order of accuracy difference problem (2.4).

For a self-adjoint positive definite operator A, it follows that [22] D = 1 2 ( τ C + 4 C + τ 2 C 2 ) is a self-adjoint positive definite operator, where C = A + τ 2 12 A 2 , R = ( I + τ D ) 1 , I is the identity operator. Moreover, the bounded operator D is defined on the whole space H.

Now we give some lemmas that will be needed below.

Lemma 2.1The following estimates hold[23]:

exp ( k τ A 1 2 ) R k H H M ( δ ) τ τ k , k 1 , ( I R 2 N ) 1 H H M ( δ ) , k τ D R k H H M ( δ ) , R k H H M ( δ ) ( 1 + δ τ ) k , k 1 , δ > 0 , D β ( R k + r R k ) H H M ( δ ) ( r τ ) α ( k τ ) α + β , 1 k < k + r N , 0 α , β 1 .

Lemma 2.2The following estimate holds[23]:

j = 1 N 1 τ D R j H H M ( δ ) Y ( τ , δ ) ,

where

Y ( τ , δ ) = min { ln 1 τ , 1 + | ln B H H | } .

Lemma 2.3For 1 l N 1 , the operator

S 1 = I R 2 N ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) ( R l 1 + R 2 N l + 1 R N l + 1 + R N + l 1 ) ( 1 1 2 ( λ τ l ) 2 ) ( R l + R 2 N l R N l + R N + l ) ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) ( R l + 1 + R 2 N l 1 R N l 1 + R N + l + 1 )

has an inverse such that

G 1 = S 1 1

and the estimate

G 1 H H M ( δ ) (2.9)

is valid.

Proof We have

G 1 G = G 1 G K 1 , (2.10)

where

G = I R 2 N R l + R 2 N l R N l + R N + l , K 1 = ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) ( R l 1 + R 2 N l + 1 R N l + 1 + R N + l 1 ) K 1 = ( 1 1 2 ( λ τ l ) 2 ) ( R l + R 2 N l R N l + R N + l ) ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) K 1 = × ( R l + 1 + R 2 N l 1 R N l 1 + R N + l + 1 ) . (2.11)

Applying estimates of Lemma 2.1, we have

K 1 H H = ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) ( R l 1 + R 2 N l + 1 R N l + 1 + R N + l 1 ) ( 1 1 2 ( λ τ l ) 2 ) ( R l + R 2 N l R N l + R N + l ) ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) ( R l + 1 + R 2 N l 1 R N l 1 + R N + l + 1 ) H H M 1 ( δ ) τ . (2.12)

By using the triangle inequality, formula (2.10), estimates (2.9), (2.12), and Lemma 2.2 of paper [15], we obtain

G 1 H H = G H H + G 1 H H G H H M ( δ ) + G 1 H H M ( δ ) M 1 ( δ ) τ

for any small positive parameter τ. From that follows estimate (2.9). Lemma 2.3 is proved. □

Lemma 2.4For 1 l N 1 , the operator

S 2 = I R 2 N ( 1 12 ( λ τ l ) 1 12 ( λ τ l ) 3 ) × ( R l 2 + R 2 N l + 2 R N l + 2 + R N + l 2 ) ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 + 1 6 ( λ τ l ) 3 ) ( R l 1 + R 2 N l + 1 R N l + 1 + R N + l 1 ) + ( λ τ l ) 2 × ( R l + R 2 N l R N l + R N + l ) ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 1 6 ( λ τ l ) 3 ) × ( R l + 1 + R 2 N l 1 R N l 1 + R N + l + 1 ) ( 1 12 ( λ τ l ) + 1 12 ( λ τ l ) 3 ) × ( R l + 2 + R 2 N l 2 R N l 2 + R N + l + 2 )

has an inverse

G 2 = S 2 1

and the estimate

G 2 H H M ( δ ) (2.13)

is satisfied.

Proof We can get

G 2 G = G 2 G K 2 , (2.14)

where G is defined by formula (2.11) and

K 2 = ( 1 12 ( λ τ l ) 1 12 ( λ τ l ) 3 ) ( R l 2 + R 2 N l + 2 R N l + 2 + R N + l 2 ) ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 + 1 6 ( λ τ l ) 3 ) × ( R l 1 + R 2 N l + 1 R N l + 1 + R N + l 1 ) + ( λ τ l ) 2 ( R l + R 2 N l R N l + R N + l ) ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 1 6 ( λ τ l ) 3 ) × ( R l + 1 + R 2 N l 1 R N l 1 + R N + l + 1 ) ( 1 12 ( λ τ l ) + 1 12 ( λ τ l ) 3 ) ( R l + 2 + R 2 N l 2 R N l 2 + R N + l + 2 ) .

Applying estimates of Lemma 2.1, we have

K 2 H H M 2 ( δ ) τ . (2.15)

Using the triangle inequality, formula (2.14), estimates (2.13), (2.15), and Lemma 2.3 of paper [15], we get

G 2 H H = G H H + G 2 H H G H H M ( δ ) + G 2 H H M ( δ ) M 2 ( δ ) τ

for any small positive parameter τ. From that follows estimate (2.13). Lemma 2.4 is proved. □

Let C τ ( H ) and C τ α , α ( H ) be the spaces of all H-valued grid functions { θ k } k = 1 N 1 in the corresponding norms,

{ θ k } k = 1 N 1 C τ ( H ) = max 1 k N 1 θ k H , { θ k } k = 1 N 1 C τ α , α ( H ) = { θ k } k = 1 N 1 C τ ( H ) { θ k } k = 1 N 1 C τ α , α ( H ) = + sup 1 k < k + n N 1 ( k τ + n τ ) α ( T k τ ) α θ k + n θ k H ( n τ ) α .

Theorem 2.1Assume thatAis a self-adjoint positive definite operator, φ , ψ , ξ D ( A ) and { θ k } k = 1 N 1 C τ α , α ( H ) ( 0 < α < 1 ). Then, the solution ( { u k } k = 1 N 1 , p ) of difference problem (2.3) obeys the following stability estimates:

{ u k } k = 1 N 1 C τ ( H ) M ( δ ) [ φ H + ψ H + ξ H + { θ k } k = 1 N 1 C τ ( H ) ] , (2.16)

A 1 p H M ( δ ) [ φ H + ψ H + ξ H + { θ k } k = 1 N 1 C τ ( H ) ] , (2.17)

p H M ( δ ) [ A φ H + A ψ H + A ξ H + 1 α ( 1 α ) { θ k } k = 1 N 1 C τ α , α ( H ) ] . (2.18)

Proof We will obtain the representation formula for the solution of problem (2.7). Applying the formula [23], we get

v k = ( I R 2 N ) 1 [ ( ( R k R 2 N k ) v 0 + ( R N k R N + k ) v N ) ( R N k R N + k ) ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ ] + ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | k i | R k + i ) θ i τ . (2.19)

By using formula (2.19) and nonlocal boundary conditions

v 0 = ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v l 1 + ( 1 1 2 ( λ τ l ) 2 ) v l v 0 = + ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v l + 1 + φ ξ , v N = v 0 + ψ φ ,

we get the system of equations

v 0 = ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) { ( I R 2 N ) 1 v 0 = × [ ( ( R l 1 R 2 N l + 1 ) v 0 + ( R N l 1 R N + l 1 ) v N ) ( R N l + 1 R N + l 1 ) v 0 = × ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ ] + ( I + τ D ) ( 2 I + τ D ) 1 v 0 = × D 1 i = 1 N 1 ( R | l 1 i | R l 1 + i ) θ i τ } + ( 1 1 2 ( λ τ l ) 2 ) { ( I R 2 N ) 1 v 0 = × [ ( ( R l R 2 N l ) v 0 + ( R N l R N + l ) v N ) ( R N l R N + l ) ( I + τ D ) v 0 = × ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ ] + ( I + τ D ) ( 2 I + τ D ) 1 D 1 v 0 = × i = 1 N 1 ( R | l i | R l + i ) θ i τ } + ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) { ( I R 2 N ) 1 v 0 = ( R N l 1 R N + l + 1 ) ( I + τ D ) ( 2 I + τ D ) 1 D 1 [ ( ( R l + 1 R 2 N l 1 ) v 0 v 0 = + ( R N l 1 R N + l + 1 ) v N ) i = 1 N 1 ( R N i R N + i ) θ i τ ] v 0 = + ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | l + 1 i | R l + 1 + i ) θ i τ } + φ ξ , v N = v 0 + ψ φ . (2.20)

Solving system (2.20), we obtain

v 0 = ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) G 1 ( R N l + 1 R N + l 1 ) ( I + τ D ) v 0 = × ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ + G 1 ( I R 2 N ) ( I + τ D ) v 0 = × ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | l 1 i | R l 1 + i ) θ i τ + ( 1 1 2 ( λ τ l ) 2 ) G 1 v 0 = × ( R N l 1 R N + l + 1 ) ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ v 0 = + G 1 ( I R 2 N ) ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | l + 1 i | R l + 1 + i ) θ i τ v 0 = + ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) G 1 ( R N l 1 R N + l + 1 ) ( I + τ D ) v 0 = × ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ + G 1 ( I R 2 N ) ( I + τ D ) v 0 = × ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | l + 1 i | R l + 1 + i ) θ i τ + G 1 ( I R 2 N ) ( φ ξ ) v 0 = + G 1 ( ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) ( R N l + 1 R N + l 1 ) v 0 = + ( 1 1 2 ( λ τ l ) 2 ) ( R N l 1 R N + l + 1 ) v 0 = + ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) ( R N l 1 R N + l + 1 ) ) ( ψ φ ) , (2.21)

v N = v 0 + ψ φ . (2.22)

Therefore, difference problem (2.7) has a unique solution { v k } k = 0 N which is defined by formulas (2.19), (2.21), and (2.22). Applying formulas (2.19), (2.21), (2.22), and the method of the monograph [23], we get

{ v k } k = 1 N 1 C τ ( H ) M ( δ ) [ φ H + ψ H + ξ H + { θ k } k = 1 N 1 C τ ( H ) ] . (2.23)

The proofs of estimates (2.17), (2.18) are based on formula (2.5) and estimate (2.23). Using formula (2.5) and estimates (2.23), (2.17), we obtain inequality (2.16). Theorem 2.1 is proved. □

Theorem 2.2Suppose thatAis a self-adjoint positive definite operator, φ , ψ , ξ D ( A ) and { θ k } k = 1 N 1 C τ α , α ( H ) ( 0 < α < 1 ). Then, the solution ( { u k } k = 1 N 1 , p ) of difference problem (2.4) obeys the stability estimates (2.16), (2.17), and (2.18).

Proof By using the representation formula (2.19) for the solution of (2.8), formula (2.19), and the nonlocal boundary conditions

v 0 = ( 1 12 ( λ τ l ) 1 12 ( λ τ l ) 3 ) v l 2 v 0 = + ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 + 1 6 ( λ τ l ) 3 ) v l 1 v 0 = + ( 1 ( λ τ l ) 2 ) v l + ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 1 6 ( λ τ l ) 3 ) v l + 1 v 0 = + ( 1 12 ( λ τ l ) + 1 12 ( λ τ l ) 3 ) v l + 2 + φ ξ , v N = v 0 + ψ φ ,

we obtain the system of equations

v 0 = ( 1 12 ( λ τ l ) 1 12 ( λ τ l ) 3 ) { ( I R 2 N ) 1 v 0 = × [ ( ( R l 2 R 2 N l + 2 ) v 0 + ( R N l 2 R N + l 2 ) v N ) ( R N l + 2 R N + l 2 ) ( I + τ D ) v 0 = × ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ ] + ( I + τ D ) ( 2 I + τ D ) 1 D 1 v 0 = × i = 1 N 1 ( R | l 2 i | R l 2 + i ) θ i τ } ( R N l 1 R N + l + 1 ) ( I + τ D ) v 0 = × ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ + ( I + τ D ) ( 2 I + τ D ) 1 D 1 v 0 = × i = 1 N 1 ( R | l + 1 i | R l + 1 + i ) θ i τ + ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 + 1 6 ( λ τ l ) 3 ) v 0 = × { ( I R 2 N ) 1 [ ( ( R l 1 R 2 N l + 1 ) v 0 + ( R N l 1 R N + l 1 ) v N ) v 0 = ( R N l + 1 R N + l 1 ) ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ ] v 0 = + ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | l 2 i | R l 2 + i ) θ i τ } v 0 = + ( 1 ( λ τ l ) 2 ) { ( I R 2 N ) 1 [ ( ( R l R 2 N l ) v 0 + ( R N l R N + l ) v N ) v 0 = ( R N l R N + l ) ( I + τ D ) ( 2 I + τ D ) 1 D 1 v 0 = × i = 1 N 1 ( R N i R N + i ) θ i τ ] + ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | l i | R l + i ) θ i τ } v 0 = + ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 1 6 ( λ τ l ) 3 ) { ( I R 2 N ) 1 v 0 = × [ ( ( R l + 1 R 2 N l 1 ) v 0 + ( R N l 1 R N + l + 1 ) v N ) ( R N l 1 R N + l + 1 ) ( I + τ D ) v 0 = × ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ ] + ( I + τ D ) ( 2 I + τ D ) 1 D 1 v 0 = × i = 1 N 1 ( R | l + 1 i | R l + 1 + i ) θ i τ } + ( 1 12 ( λ τ l ) + 1 12 ( λ τ l ) 3 ) { ( I R 2 N ) 1 v 0 = × [ ( ( R l + 2 R 2 N l 2 ) v 0 + ( R N l 2 R N + l + 2 ) v N ) ( R N l 2 R N + l + 2 ) v 0 = × ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ ] v 0 = + ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | l + 2 i | R l + 2 + i ) θ i τ } + φ ξ , v N = v 0 + ψ φ . (2.24)

Solving system (2.24), we have

v 0 = ( 1 12 ( λ τ l ) 1 12 ( λ τ l ) 3 ) G 2 ( R N l + 2 R N + l 2 ) ( I + τ D ) v 0 = × ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ + G 2 ( I R 2 N ) ( I + τ D ) v 0 = × ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | l 2 i | R l 2 + i ) θ i τ v 0 = + ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 + 1 6 ( λ τ l ) 3 ) G 2 ( R N l + 1 R N + l 1 ) v 0 = × ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ + G 2 ( I R 2 N ) ( I + τ D ) v 0 = × ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | l 1 i | R l 1 + i ) θ i τ + ( 1 ( λ τ l ) 2 ) G 2 v 0 = × ( R N l R N + l ) ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ v 0 = + G 2 ( I R 2 N ) ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | l i | R l + i ) θ i τ v 0 = + ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 1 6 ( λ τ l ) 3 ) G 2 ( R N l 1 R N + l + 1 ) v 0 = × ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ + G 2 ( I R 2 N ) ( I + τ D ) v 0 = × ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | l + 1 i | R l + 1 + i ) θ i τ + ( 1 12 ( λ τ l ) + 1 12 ( λ τ l ) 3 ) v 0 = × G 2 ( R N l 2 R N + l + 2 ) ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R N i R N + i ) θ i τ v 0 = + G 2 ( I R 2 N ) ( I + τ D ) ( 2 I + τ D ) 1 D 1 i = 1 N 1 ( R | l + 2 i | R l + 2 + i ) θ i τ v 0 = + G 2 ( I R 2 N ) ( φ ξ ) + G 2 ( ( 1 12 ( λ τ l ) 1 12 ( λ τ l ) 3 ) ( R N l 2 R N + l + 2 ) v 0 = + ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 + 1 6 ( λ τ l ) 3 ) ( R N l + 1 R N + l 1 ) v 0 = + ( 1 ( λ τ l ) 2 ) ( R N l R N + l ) v 0 = + ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 1 6 ( λ τ l ) 3 ) ( R N l 1 R N + l + 1 ) v 0 = + ( 1 12 ( λ τ l ) + 1 12 ( λ τ l ) 3 ) ( R N l 2 R N + l + 2 ) ) ( ψ φ ) , (2.25)

v N = v 0 + ψ φ . (2.26)

So, the difference problem (2.8) has a unique solution { v k } k = 0 N , which is defined by formulas (2.19), (2.25), and (2.26). By using formulas (2.19), (2.25), (2.26), and the method of the monograph [23], we can get the stability estimate (2.23) for the solution of difference problem (2.8). The proofs of estimates (2.17), (2.18) are based on (2.5) and (2.23). Applying formula (2.5) and estimates (2.23), (2.17), we get estimate (2.16). Theorem 2.2 is proved. □

Theorem 2.3Assume thatAis a self-adjoint positive definite operator, φ , ψ , ξ D ( A ) and { θ k } k = 1 N 1 C τ ( H ) . Then, the solutions ( { u k } k = 1 N 1 , p ) of difference problems (2.3) and (2.4) obey the following almost coercive inequality:

{ τ 2 ( u k + 1 2 u k + u k 1 ) } k = 1 N 1 C τ ( H ) + { ( A + τ 2 12 A 2 ) u k } k = 1 N 1 C τ ( H ) + p H M ( δ ) [ min { ln 1 τ , 1 + | ln B H H | } { θ k } k = 1 N 1 C τ ( H ) + ( A + τ 2 12 A 2 ) φ H + ( A + τ 2 12 A 2 ) ψ H + ( A + τ 2 12 A 2 ) ξ H ] . (2.27)

Theorem 2.4Assume thatAis a self-adjoint positive definite operator, φ , ψ , ξ D ( A ) and { θ k } k = 1 N 1 C τ α , α ( H ) ( 0 < α < 1 ). Then, the solutions ( { u k } k = 1 N 1 , p ) of difference problems (2.3) and (2.4) obey the following coercive inequality:

{ τ 2 ( u k + 1 2 u k + u k 1 ) } k = 1 N 1 C τ α , α ( H ) + { ( A + τ 2 12 A 2 ) u k } k = 1 N 1 C τ α , α ( H ) + p H M ( δ ) [ 1 α ( 1 α ) { θ k } k = 1 N 1 C τ α , α ( H ) + ( A + τ 2 12 A 2 ) φ H + ( A + τ 2 12 A 2 ) ψ H + ( A + τ 2 12 A 2 ) ξ H ] . (2.28)

The proofs of Theorems 2.3 and 2.4 are based on formulas (2.5), (2.19), (2.21), (2.22), (2.25), (2.26), Lemmas 2.1 and 2.2.

### 3 High order of accuracy difference schemes for the problem (1.2) and their well-posedness

Now, we consider problem (1.2). The differential expression [22,23]

A x u ( x ) = r = 1 n ( a r ( x ) u x r ) x r + σ u

defines a self-adjoint strongly positive definite operator A x acting on L 2 ( Ω ¯ ) with the domain

D ( A x ) = { u ( x ) W 2 ( Ω ¯ ) , u ( x ) = 0 , x S } .

The discretization of problem (1.2) is carried out in two steps. In the first step, we define the grid spaces

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

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

A h x u h ( x ) = r = 1 n ( a r ( x ) u x ¯ r h ) x r , j r + σ u h ( x ) (3.1)

acting in the space of grid functions u h ( x ) , satisfying the condition u h ( x ) = 0 for all x S h .

To formulate our results, let L 2 h = L 2 ( Ω ˜ h ) and W 2 h 2 = W 2 2 ( Ω ˜ h ) be spaces of the grid functions ζ h ( x ) = { ζ ( h 1 m 1 , , h n m n ) } defined on Ω ˜ h , equipped with the norms

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

Applying formula (2.5) to A h x , we arrive for v h ( t , x ) functions, at auxiliary nonlocal boundary value problem for a system of ordinary differential equations

{ d 2 v h ( t , x ) d t 2 + A h x v h ( t , x ) = f h ( t , x ) , 0 < t < T , x Ω ˜ h , v h ( 0 , x ) v h ( λ , x ) = φ ( x ) ξ ( x ) , x Ω ˜ h v h ( T , x ) v h ( λ , x ) = ψ ( x ) ξ ( x ) , x Ω ˜ h . (3.2)

We define function p h ( x ) by formula

p h ( x ) = A h x φ h ( x ) A h x v h ( 0 , x ) , x Ω ˜ h . (3.3)

In the second step, auxiliary nonlocal problem (3.2) is replaced by the third order of accuracy difference scheme

v k + 1 h ( x ) 2 v k h ( x ) + v k 1 h ( x ) τ 2 + A h x v k h ( x ) + τ 2 12 ( A h x ) 2 v k h ( x ) = θ k h ( x ) , θ k h ( x ) = f h ( t k , x ) + τ 2 12 ( f h ( t k + 1 , x ) 2 f h ( t k , x ) + f h ( t k 1 , x ) τ 2 + A h x f h ( t k , x ) ) , v 0 h ( x ) ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v l 1 h ( x ) ( 1 1 2 ( λ τ l ) 2 ) v l h ( x ) ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v l + 1 h ( x ) = φ h ( x ) ξ h ( x ) , x Ω ˜ h , v N h ( x ) ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v l 1 h ( x ) ( 1 1 2 ( λ τ l ) 2 ) v l h ( x ) ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v l + 1 h ( x ) = φ h ( x ) ξ h ( x ) , x Ω ˜ h , l = [ λ τ ] , (3.4)

and by the fourth order of accuracy difference scheme

v k + 1 h ( x ) 2 v k h ( x ) + v k 1 h ( x ) τ 2 + A h x v k h ( x ) + τ 2 12 ( A h x ) 2 v k h ( x ) = θ k h ( x ) , θ k h ( x ) = f h ( t k , x ) + τ 2 12 ( f h ( t k + 1 , x ) 2 f h ( t k , x ) + f h ( t k 1 , x ) τ 2 + A h x f h ( t k , x ) ) , t k = k τ , 1 k N 1 , x Ω ˜ h , v 0 h ( x ) ( 1 12 ( λ τ l ) 1 12 ( λ τ l ) 3 ) v l 2 h ( x ) ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 + 1 6 ( λ τ l ) 3 ) v l 1 h ( x ) ( 1 ( λ τ l ) 2 ) v l h ( x ) ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 1 6 ( λ τ l ) 3 ) v l + 1 h ( x ) ( 1 12 ( λ τ l ) + 1 12 ( λ τ l ) 3 ) v l + 2 h ( x ) = φ h ( x ) ξ h ( x ) , x Ω ˜ h , v N h ( x ) ( 1 12 ( λ τ l ) 1 12 ( λ τ l ) 3 ) v l 2 h ( x ) ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 + 1 6 ( λ τ l ) 3 ) v l 1 h ( x ) ( 1 ( λ τ l ) 2 ) v l h ( x ) ( 8 12 ( λ τ l ) + 1 2 ( λ τ l ) 2 1 6 ( λ τ l ) 3 ) v l + 1 h ( x ) ( 1 12 ( λ τ l ) + 1 12 ( λ τ l ) 3 ) v l + 2 h ( x ) = φ h ( x ) ξ h ( x ) , x Ω ˜ h . (3.5)

Let τ and | h | = h 1 2 + + h n 2 be sufficiently small positive numbers.

Theorem 3.1The solutions of difference schemes (3.4) and (3.5) obey the following stability estimates:

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

Theorem 3.2The solutions of difference schemes (3.4) and (3.5) obey the following almost coercive stability estimate:

{ u k + 1 h 2 u k h + u k 1 k τ 2 } 1 N 1 C τ ( L 2 h ) + { ( A + τ 2 12 A 2 ) u k } k = 1 N 1 C τ ( W 2 h 2 ) + p h L 2 h M ( δ ) [ ln ( 1 τ + h ) { f k h } 1 N C τ ( L 2 h ) + ( A + τ 2 12 A 2 ) φ h W 2 h 2 + ( A + τ 2 12 A 2 ) ψ h W 2 h 2 + ( A + τ 2 12 A 2 ) ξ h W 2 h 2 ] .

Theorem 3.3The solutions of difference schemes (3.4) and (3.5) obey the following coercive stability estimate:

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

The proofs of Theorems 3.1-3.3 are based on the abstract Theorems 2.1-2.4, symmetry properties of the operator A h x in L 2 h and the following theorem on the coercivity inequality for the solution of the elliptic difference problem in L 2 h .

Theorem 3.4[24]

For the solution 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:

r = 1 n ( u k h ) x ¯ r x ¯ r , j r L 2 h M ω h L 2 h ,

whereMdoes not depend onhand ω h .

### 4 Numerical results

In this section, by using the third and fourth order of the accuracy approximation, we obtain an approximate solution of the inverse problem

{ 2 u ( t , x ) t 2 x ( u ( t , x ) x ) + u ( t , x ) = f ( t , x ) + p ( x ) , 0 < x < π , 0 < t < T , f ( t , x ) = ( exp ( t ) + 2 t ) sin ( x ) , u ( 0 , x ) = 2 sin ( x ) , 0 x π , u ( T , x ) = ( exp ( T ) + T + 1 ) sin ( x ) , 0 x π , u ( λ , x ) = ( exp ( λ ) + λ + 1 ) sin ( x ) , 0 x π , u ( t , 0 ) = u ( t , π ) = 0 , 0 t T ( T = 2 , λ = 4 7 T ) (4.1)

for the elliptic equation. Note that u ( t , x ) = ( exp ( t ) + t + 1 ) sin ( x ) and p ( x ) = 2 sin ( x ) are the exact solutions of equation (4.1).

For the approximate solution of the nonlocal boundary value problem (3.2), consider the set of grid points

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

which depends on the small parameters τ = T N and h = π M .

Applying approximations (3.4) and (3.5), we get, respectively, the third order of the accuracy difference scheme

{ v n k + 1 2 v n k + v n k 1 τ 2 v n + 1 k 2 v n k + v n 1 k h 2 + v n k τ 2 12 [ 1 h 2 ( v n + 2 k 2 v n + 1 k + v n k h 2 + v n + 1 k ) + 2 h 2 ( v n + 1 k 2 v n k + v n 1 k h 2 + v n k ) 1 h 2 ( v n k 2 v n 1 k + v n 2 k h 2 + v n 1 k ) v n + 1 k 2 v n k + v n 1 k h 2 + v n k ] = ( exp ( t k ) + 2 t k ) sin ( x n ) + τ 2 12 ( exp ( t k ) + 4 t k ) sin ( x n ) , = k = 1 , , N 1 , n = 1 , , M 2 , v 0 k = v M k = 0 , v 1 k = 4 5 v 2 k 1 5 v 3 k , v M 1 k = 4 5 v M 2 k 1 5 v M 3 k , k = 0 , , N , v n 0 ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v n l 1 ( 1 1 2 ( λ τ l ) 2 ) v n l ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v n l + 1 = ( 1 exp ( λ ) λ ) sin ( x n ) , n = 0 , , M , v n N ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v n l 1 ( 1 1 2 ( λ τ l ) 2 ) v n l ( 1 2 ( λ τ l ) + 1 2 ( λ τ l ) 2 ) v n l + 1 = ( exp ( t N ) exp ( λ ) + t N λ ) sin ( x n ) , = n = 0 , , M , (4.2)

and the fourth order of the accuracy difference scheme

{ v n k + 1 2 v n k + v n k 1 τ 2 v n + 1 k 2 v n k + v n 1 k h 2 + v n k τ 2 12 [ 1 h 2 ( v n + 2 k 2 v n + 1 k + v n k h 2 + v n + 1 k ) + 2 h 2 ( v n + 1 k 2 v n k + v n 1 k h 2 + v n k ) 1 h 2 ( v n k 2 v n 1 k + v n 2 k h 2 + v n 1 k ) v n + 1 k 2 v n k + v n 1 k h 2 + v n k ] = ( exp ( t k ) + 2 t k ) sin ( x n ) + τ 2 12 ( exp ( t k ) + 4 t k ) sin ( x n ) , = k = 1 , , N 1 , n = 2 , , M 2 , v 0 k = v M k = 0 , v 1 k = 4 5 v 2 k 1 5 v 3 k , v M 1 k = 4 5 v M 2 k 1 5 v M 3 k , k = 0 , , N , v n 0 ( 1 12 ( λ τ l ) 1 12 ( λ τ l )