Research

# On solvability of a nonlocal problem for the Laplace equation with the fractional-order boundary operator

Moldir A Muratbekova*, Kanat M Shinaliyev and Batirkhan K Turmetov

Author Affiliations

Department of Mathematics, Akhmet Yasawi International Kazakh-Turkish University, B. Sattarkhanov Street 29, Turkistan, 161200, Kazakhstan

For all author emails, please log on.

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

The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2014/1/29

 Received: 15 August 2013 Accepted: 9 January 2014 Published: 30 January 2014

© 2014 Muratbekova et al.; licensee Springer.

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

In the present work, we study properties of some integro-differential operators of the Hadamard-Marchaud type in the class of harmonic functions. As an application of these properties, we consider the question of the solvability of a nonlocal boundary value problem for the Laplace equation in the unit ball.

MSC: 35J05, 35J25, 26A33.

##### Keywords:
Hadamard-Marchaud operator; fractional derivative; nonlocal problem

### 1 Introduction

Let Ω = { x R n : | x | < 1 } be the unit ball, n 2 . A paradigm in the theory of elliptic partial differential equations and harmonic functions is the Laplace equation

Δ u ( x ) = 0 , x Ω . (1)

If we prescribe the values of the solution at the boundary Ω of Ω, then we can solve equation (1) uniquely. Of course, one can consider many other boundary conditions such as Neumann’s boundary conditions.

In some applied problems of hydrodynamics [1], it is necessary to prescribe the value of a fractional derivative of the solution on the boundary. Fractional differential equations and boundary value problems involving fractional derivatives appear in many applied problems ranging from the spring-pot model [2] to geology [3] or from nonlinear circuits [4] to alternative models to differential equations [5].

Hence, in this paper we study the Laplace equation concentrating on some conditions on the boundary involving derivatives of fractional order.

Note that numerous works of authors [6-13] were dedicated to the research questions of the solvability of boundary value problems for partial differential equations with boundary operators of high (whole and fractional) order. In the paper of A.N. Tikhonov [14] boundary value problems with boundary conditions containing derivatives of higher order have been investigated for the heat equation. Research questions as regards the solvability of similar problems for higher-order equations with boundary operators of whole and fractional order were carried out in [15,16]. Later in [17], these results were generalized for partial differential equations of fractional order. In [18-20] questions about the solvability of boundary value problems with boundary operators of high order were studied for the Laplace equation. In the studies of these authors the exact conditions for the solvability have been established and the integral representations of solutions of the studied problems have been found. The cycle of studies by the authors [21-26] is devoted to the study of the existence and smoothness of solutions of boundary value problems for the second-order elliptic equations with boundary operators of fractional order. In the paper mentioned above local boundary value problems with boundary operators of fractional order in the Riemann-Liouville or Caputo sense are studied. In this paper we study nonlocal problems with boundary operators of fractional-order derivatives of Hadamard type. Definitions of Hadamard operators, a statement of the main problems, and the history of the questions on this topic are in Section 3.

The organization of this paper is as follows. In Section 2, we present the operators of integration and differentiation in the Hadamard sense and some modifications. In the third section we provide a formulation of the basic problem of this paper and some historical information as regards nonlocal boundary value problems. In the fourth section we study the properties of integral and differential Hadamard-Marchaud operators in the class of harmonic functions in the ball. In Section 5 we provide some auxiliary propositions. Finally, Section 6 is devoted to the study of the fundamental problem, where we formulate and prove the main statement of the paper.

### 2 Definition of Hadamard operators of integration and differentiation and some modifications

In this section, we give a statement on the operators of fractional differentiation in the sense of Hadamard, Hadamard-Marchaud, and their modifications.

For any positive α, fractional integrals and derivatives of the order α in the sense of Hadamard are defined by the following formulas [27]:

J α φ ( t ) = 1 Γ ( α ) 0 t ( ln t s ) ( α + 1 ) φ ( s ) s d s , t > 0 , (2)

D α φ ( t ) = δ m + 1 J m + 1 α φ ( t ) , (3)

where δ = t d d t is the Dirac operator, m = [ α ] is the integral part of α.

If 0 < α < 1 , then, in the class of sufficiently ‘good’ functions, operator (3) can be reduced to the following form [27]:

D α φ ( t ) = δ J 1 α φ ( t ) = α Γ ( 1 α ) 0 1 φ ( t ) φ ( s t ) s | ln s | α + 1 d s . (4)

This operator is said to be the differentiation operator of order α in the sense of Hadamard-Marchaud.

In [28], the following modification of the Hadamard-Marchaud operator was considered:

D μ α φ ( t ) = α Γ ( 1 α ) 0 1 φ ( t ) φ ( s t ) s 1 μ | ln s | α + 1 d s + μ α φ ( t ) , μ 0 . (5)

In [18], in the class of harmonic functions in a ball, the properties and applications of the operators in the form of

δ μ = r r + μ , δ μ m = ( r r + μ ) m (6)

are considered. Here μ > 0 , r = | x | , x = ( x 1 , , x n ) , and r r is a differential operator in the form of r r = i = 1 n x i x i .

Let u ( x ) be a harmonic function in the domain Ω, and let 0 α , 0 < μ be arbitrary real numbers. Let us consider a modification of the Bavrin operator (6).

Introduce the operators

J μ α [ u ] ( x ) = { u ( x ) , α = 0 , 1 Γ ( α ) 0 1 | ln s | α 1 s μ 1 u ( s x ) d s , α > 0 , D μ α [ u ] ( x ) = { u ( x ) , α = 0 , α Γ ( 1 α ) 0 1 u ( x ) u ( s x ) s 1 μ | ln s | α + 1 d s + μ α u ( x ) , 0 < α < 1 , ( r r + μ ) m D μ α m [ u ] ( x ) , m α < m + 1 , m = 1 , 2 , .

If α = m , μ > 0 , then we obtain the Bavrin operator D μ m [ u ] ( x ) = ( r r + μ ) m .

### 3 Statement of the problem

Let t k : Ω Γ k Ω ¯ 1 Ω , k = 1 , 2 ,  , Γ k , be continuous mappings, and let a k ( x ) be continuous functions satisfying the condition

k = 1 | a k ( x ) | < , x Ω . (7)

We assume that the series (7) converges uniformly on Ω.

Further, let μ > 0 , 0 β α , and α + β 0 , i.e.α and β are not equal to zero simultaneously.

Consider the following boundary value problem:

Δ u ( x ) = 0 , x Ω , (8)

D μ α [ u ] ( x ) k = 1 a k ( x ) D μ β [ u ] ( t k ( x ) ) = f ( x ) , x Ω , (9)

where f ( x ) C ( Ω ) .

A harmonic function u ( x ) from the class C 2 ( Ω ) C ( Ω ¯ ) , such that D μ α u ( x ) C ( Ω ¯ ) and condition (9) is realized in the classical sense, will be called a solution of problem (8)-(9).

The above-mentioned problem is a simple generalization of Bitsadze-Samarskii’s nonlocal problem [29]. For convenience of the reader, we formulate the Bitsadze-Samarskii problem.

Let D be a finite simply-connected domain of the plane of complex variables z = x + i y with the smooth boundary S = D , and let S 1 be a closed simple smooth curve lying in D.

We denote by τ = ν ( t ) , t S , τ S 1 a diffeomorphism between S and S 1 .

Formulation of the problem: We are to find a harmonic function u ( z ) = u ( x , y ) in D, which is continuous in D ¯ = D S and satisfies the boundary condition

u ( t ) u ( ν ( t ) ) = f ( t ) , t S ,

where f ( t ) is a given function.

Similar problems with operators of integer order were considered in [30-32], and for operators of fractional order with fractional-order derivatives in the sense of Riemann-Liouville and Caputo in [33-41]. It should also be noted that some questions of solvability of nonlocal problems for fractional-order equations in the one-dimensional case were studied in [42-44].

### 4 Properties of operators J μ α and D μ α

In this section, we study some properties of the operators J μ α and D μ α in the class of harmonic functions. Further, for convenience, we shall take everywhere γ = α m .

Lemma 1Let α > 0 , μ > 0 , and H k ( x ) be a homogeneous harmonic polynomial of the power k N 0 = { 0 , 1 , } . Then the following equalities are correct:

J μ α [ H k ] ( x ) = ( k + μ ) α H k ( x ) , μ > 0 , k N 0 , (10)

D μ α [ H k ] ( x ) = ( k + μ ) α H k ( x ) , k N 0 , μ > 0 . (11)

Proof Let μ > 0 . Then, using homogeneity of the polynomial H k ( x ) , we obtain

J μ α [ H k ] ( x ) = 1 Γ ( α ) 0 1 | ln s | α 1 s μ 1 H k ( s x ) d s = H k ( x ) Γ ( α ) 0 1 | ln s | α 1 s k + μ 1 d s .

The value of the last integral can easily be calculated with the help of the change of variables z = ln s . In fact,

1 Γ ( α ) 0 1 | ln s | α 1 s k + μ 1 d s = 1 Γ ( α ) 0 z α 1 e ( k + μ ) z d z = ( k + μ ) α Γ ( α ) 0 t α 1 e t d t = ( k + μ ) α .

The equality (10) is proved.

Further, note that the relation

δ μ m H k ( x ) ( r r + μ ) m H k ( x ) = ( k + μ ) m H k ( x ) , k N 0 , μ 0 , (12)

holds for the operator δ μ m .

Now, let us study actions of the operator D μ γ to the functions H k ( x ) . Using the definition of D μ γ and the homogeneity of H k ( x ) , we have

D μ γ [ H k ] ( x ) = γ Γ ( 1 γ ) 0 1 H k ( x ) H k ( s x ) s 1 μ | ln s | γ + 1 d s + μ γ H k ( x ) = γ H k ( x ) Γ ( 1 γ ) 0 1 ( s μ 1 s k + μ 1 ) | ln s | ( γ + 1 ) d s + μ γ H k ( x ) .

Denoting I = 0 1 ( s μ s k + μ ) ( ln 1 s ) ( γ + 1 ) d s s and integrating by parts, we get

I = 1 γ 0 1 [ μ s μ 1 ( k + μ ) s k + μ 1 ] ( ln 1 s ) γ d s .

After the change of variables z = ln s , as in the proof of equality (10), we easily obtain

I = Γ ( 1 γ ) γ μ γ + Γ ( 1 γ ) γ ( k + μ ) γ ,

which implies D μ γ [ H k ] ( x ) = ( k + μ ) γ H k ( x ) .

Further, taking into account fulfilling of equality (12), we obtain in the general case for m α < m + 1 :

D μ α [ H k ] ( x ) = δ μ m D μ γ [ H k ] ( x ) = ( k + μ ) α H k ( x ) .

The lemma is proved. □

Lemma 2Let α > 0 , μ > 0 , and u ( x ) be a harmonic function in the ball Ω. Then the functions D μ α [ u ] ( x ) and J μ α [ u ] ( x ) are also harmonic in Ω.

Proof Let u ( x ) be a harmonic function in the ball Ω. Then it is known [45] that the function u ( x ) is represented in the form of the series

u ( x ) = k = 0 i = 1 H k u k ( i ) H k ( i ) ( x ) , (13)

where { H k ( i ) ( x ) : i = 1 , , h k } is a complete system of homogeneous harmonic polynomials of power k, and u k ( i ) are coefficients of the expansion (13). Applying formally the operator D μ α to the series (13) and taking into account equality (11), we obtain

D μ α [ u ] ( x ) = k = 0 i = 1 H k ( k + μ ) α u k ( i ) H k ( i ) ( x ) . (14)

Now let us check convergence of the series (13) and (14). The following asymptotical estimate is valid for h k :

h k = ( 1 + 2 k ( n 2 ) ) C k + n 3 n 3 2 k n 2 ( n 2 ) ! , k ( n 3 ) .

Moreover, the series (13) converges absolutely and uniformly by x at | x | ρ < 1 , hence, for any | x | ρ and any i , k N 0 , the equalities | u k ( i ) H k ( i ) ( x ) | c ρ hold. Since lim k ( k + μ ) α k = 1 , we have for | x | r ρ and r < 1

k = 0 i = 1 h k ( k + μ ) α | u k ( i ) H k ( i ) ( x ) | c ρ k = 0 [ ( k + μ ) α h k ] r k < .

Therefore, the series (14) converges absolutely and uniformly by x at | x | r ρ , where r < 1 , ρ < 1 , and its sum is a harmonic function. By virtue of the arbitrariness of r < 1 and ρ < 1 , the function D μ α [ u ] ( x ) is defined in the whole ball Ω.

Let us study the function J μ α . Applying formally the operator J μ α to the series (13), taking into account equality (10), we obtain

J μ α [ u ] ( x ) = k = 1 i = 1 h k ( k + μ ) α u k ( i ) H k ( i ) ( x ) .

Convergence of this series can be checked as in the case of series (14), and that is why J μ α [ u ] ( x ) is a harmonic function in the ball Ω. The lemma is proved. □

Now we show that the function u ( x ) can be represented in terms of the function D μ α [ u ] ( x ) .

Lemma 3Let α > 0 , μ > 0 , and u ( x ) be a harmonic function in the domain Ω. Then for any x Ω the equality

u ( x ) = 1 Γ ( α ) 0 1 | ln s | α 1 s μ 1 D μ α [ u ] ( s x ) d s

is valid.

Proof Let μ > 0 . Represent a harmonic function u ( x ) in the form of the series (13) and transform it to the form of

u ( x ) = k = 0 i = 1 h k ( k + μ ) α u k ( i ) ( k + μ ) α H k ( i ) ( x ) . (15)

Further, taking into account equalities (10)-(11), and the absolute and uniform convergence of the series (15) by x at | x | ρ < 1 , it can be reduced to the form of

u ( x ) = k = 0 i = 1 h k ( k + μ ) α u k Γ ( α ) 0 1 | ln s | α 1 s μ 1 H k ( i ) ( s x ) d s = 1 Γ ( α ) 0 1 | ln s | α 1 s μ 1 D μ α [ k = 0 i = 1 h k u k ( i ) H k ( i ) ( s x ) ] d s = 1 Γ ( α ) 0 1 | ln s | α 1 s μ 1 D μ α [ u ] ( s x ) d s .

The lemma is proved. □

One can similarly prove the following lemma.

Lemma 4Let 0 < β < α , μ > 0 , and u ( x ) be a harmonic function in the domain Ω. Then for any x Ω the equality

D μ β [ u ] ( x ) = 1 Γ ( α β ) 0 1 | ln s | α β 1 s μ 1 D μ α [ u ] ( s x ) d s (16)

is valid.

Lemma 5Let α > 0 , μ > 0 , and u ( x ) be a harmonic function in the domain Ω. Then the following equalities hold:

J μ α [ D μ α [ u ] ] ( x ) = D μ α [ J μ α [ u ] ] ( x ) = u ( x ) .

Proof Let μ > 0 . Applying the operator J μ α to the function D μ α [ u ] ( x ) , we obtain

J μ α [ D μ α [ u ] ] ( x ) = 1 Γ ( α ) 0 1 | ln s | α 1 s μ 1 D μ α [ u ] ( s x ) d s .

By virtue of Lemma 3, the value of the last integral is equal to u ( x ) , i.e. J μ α [ D μ α [ u ] ] ( x ) = u ( x ) .

To prove the second equality, apply the operator D μ γ to the function J μ α [ u ] ( x ) . We get

D μ γ [ J μ α [ u ] ] ( x ) = γ Γ ( 1 γ ) 0 1 J μ α [ u ] ( x ) J μ α [ u ] ( s x ) s 1 μ ln | s | γ + 1 d s + μ γ J μ α [ u ] ( x ) = γ Γ ( 1 γ ) 0 1 1 s 1 μ ln | s | γ + 1 1 Γ ( α ) 0 1 u ( τ x ) u ( τ s x ) τ 1 μ ln | τ | 1 α d τ d s + μ γ J μ α [ u ] ( x ) = 1 Γ ( α ) 0 1 1 | ln τ | 1 α τ 1 μ γ Γ ( 1 γ ) 0 1 u ( τ x ) u ( τ s x ) | ln s | ( γ + 1 ) s 1 μ d s d τ + μ γ J μ α [ u ] ( x ) = 1 Γ ( α ) 0 1 D μ γ [ u ] ( τ x ) μ γ u ( τ x ) | ln τ | 1 α τ 1 μ d τ + μ γ J μ α [ u ] ( x ) = 1 Γ ( α ) 0 1 D μ γ [ u ] ( τ x ) d τ | ln τ | 1 α τ 1 μ μ γ 1 Γ ( α ) 0 1 u ( τ x ) d τ | ln τ | 1 α τ 1 μ + μ γ J μ α [ u ] ( x ) = 1 Γ ( α ) 0 1 D μ γ [ u ] ( τ x ) | ln τ | 1 α τ 1 μ d τ μ γ J μ α [ u ] ( x ) + μ γ J μ α [ u ] ( x ) = 1 Γ ( α ) 0 1 D μ γ [ u ] ( τ x ) | ln τ | 1 α τ 1 μ d τ .

Then, in the general case,

D μ α [ J μ α [ u ] ] ( x ) = δ μ m D μ γ [ J μ α [ u ] ] ( x ) = 1 Γ ( α ) 0 1 δ μ m D μ γ [ u ] ( τ x ) | ln τ | 1 α τ 1 μ d τ = 1 Γ ( α ) 0 1 | ln τ | α 1 τ μ 1 D μ α [ u ] ( τ x ) d τ = u ( x ) .

The lemma is proved. □

### 5 Some auxiliary propositions

Let t k and a k ( x ) satisfy the conditions from Section 2.

Consider the following problem in the domain Ω:

Δ u ( x ) = 0 , x Ω , (17)

v ( x ) k = 1 a k ( x ) J μ α β [ u ] ( t k ( x ) ) = f ( x ) , x Ω , (18)

where μ > 0 , 0 β α , and α + β 0 , i.e.α and β are not equal to zero simultaneously, f ( x ) C ( Ω ) .

A harmonic function v ( x ) from the class C 2 ( Ω ) C ( Ω ¯ ) , satisfying condition (18) in the classical case, will be called a solution of problem (17)-(18).

It should be noted that problem (17)-(18) was investigated for the case of α = β in [30].

Let us investigate uniqueness for the solution of problem (17)-(18). The following statement holds.

Lemma 6Let Γ k Ω ¯ 1 Ω , k = 1 , 2 ,  , a k ( x ) , k = 0 , 1 , 2 ,  , be continuous functions satisfying the condition

k = 1 | a k ( x ) | μ α β , x Ω , (19)

and let a solution of problem (17)-(18) exist.

Then:

(1) If

k = 1 a k ( x ) μ α β , x Ω , (20)

then the solution of problem (17)-(18) is unique.

(2) If

k = 1 a k ( x ) μ α β , x Ω , (21)

then the solution of problem (17)-(18) is unique up to a constant summand.

Proof Let v ( x ) be the solution of problem (17)-(18) at f ( x ) = 0 .

Denote M = | v ( x 0 ) | = max Ω | v ( x ) | , x 0 Ω .

Then if v ( x ) const , then, by virtue of the maximum principle for harmonic functions [46], the inequality | v ( x ) | < M holds for any x Ω .

The boundary condition (18) at f ( x ) = 0 implies

M = | k = 1 a k ( x 0 ) J μ α β [ v ] [ t k ( x 0 ) ] | k = 1 | a k ( x 0 ) | | J μ α β [ v ] ( t k ( x 0 ) ) | k = 1 | a k ( x 0 ) | 1 Γ ( α β ) 0 1 s μ 1 | ln s | ( α β 1 ) | v ( s t k ( x 0 ) ) | d s .

Further, since t k : Ω Γ k Ω ¯ 1 Ω , k = 1 , 2 ,  , Γ k , then t k ( x 0 ) Ω , and for any s [ 0 , 1 ] , s t k ( x 0 ) Ω . Therefore | v ( s t k ( x 0 ) ) | < M .

Hence,

M < M k = 1 | a k ( x 0 ) | 1 Γ ( α β ) 0 1 s μ 1 | ln s | α β 1 d s = M μ β α k = 1 | a k ( x 0 ) | .

If now condition (19) is realized , then μ β α k = 1 | a k ( x 0 ) | 1 , and we obtain from this the contradiction M < M .

Hence, if condition (19) holds, it is necessary that v ( x ) = C const . Since J μ α β [ C ] = μ β α C , substituting the function v ( x ) = C into the boundary condition (18), for f ( x ) = 0 we have

C μ β α C k = 1 a k ( x ) = 0 , x Ω .

The last equality is equivalent to the equality

C [ 1 μ β α k = 1 a k ( x ) ] = 0 .

We obtain from this the result that either C = 0 or k = 1 a k ( x ) = μ α β .

Thus, if conditions (19) and (20) are fulfilled, we obtain C = 0 , i.e. v ( x ) 0 .

If the conditions (21) are fulfilled, then any constant is a solution of the homogeneous problem (17)-(18). In fact, substituting v ( x ) C into equation (18), we obtain

C k = 1 a k ( x ) μ β α C = C [ 1 μ β α k = 1 a k ( x ) ] = C [ 1 μ β α μ α β ] = 0 .

The lemma is proved. □

Now investigate existence of a solution of problem (17)-(18). Let 0 β α and let P ( x , y ) = 1 ω n 1 | x | 2 | x y | n be the Poisson kernel of the Dirichlet problem, and ω n the area of the unit sphere.

Introduce the function

P α β , μ ( x , y ) = { P ( x , y ) , if  α = β , 1 Γ ( α β ) 0 1 s μ 1 | ln s | α β 1 P ( s x , y ) d s , if  α > β , (22)

and consider the equation

ψ ( x ) Ω [ k = 1 a k ( y ) P α β , μ ( t k ( y ) , x ) ] ψ ( y ) d S y = 0 . (23)

The following statement holds.

Lemma 7Let Γ k Ω ¯ 1 Ω , a k ( x ) , k = 1 , 2 ,  , be continuous functions satisfying the condition (19). Then:

(1) If the condition (20) is realized, then problem (17)-(18) is uniquely solvable at any f ( x ) C ( Ω ) .

(2) If the condition (21) is realized, then problem (17)-(18) is solvable if the following condition is realized:

Ω f ( x ) ψ 0 ( x ) d s x = 0 , (24)

where the function ψ 0 ( x ) is a solution of equation (23), moreover the number of independent solutions of this equation under these conditions is equal to 1.

Proof Since v ( x ) is a harmonic function, a solution of problem (17)-(18) can be found in the form of the Poisson integral v ( x ) = Ω P ( x , y ) ν ( x ) d s x where ν ( x ) is an unknown function. Substituting this function into the boundary condition (18), we obtain the integral equation with respect to the unknown function ν ( x ) ,

ν ( x ) Ω [ k = 1 a k ( x ) P α β , μ ( t k ( x ) , y ) ] ν ( y ) d S y = f ( x ) , x Ω . (25)

Designate

K ( x , y ) = k = 1 a k ( x ) P α β , μ ( t k ( x ) , y ) .

Then equation (25) can be rewritten in the form of

ν ( x ) + Ω K ( x , y ) ν ( y ) d S y = f ( x ) , x Ω . (26)

To investigate the solvability of the integral equation (26), we study the properties of the kernel K ( x , y ) . We show that K ( x , y ) is a continuous function on Ω × Ω .

In fact, since t k ( x ) Ω ¯ 1 Ω , we obtain | s t k ( x ) y | > 0 for all x , y Ω , and therefore the function P ( s t k ( x ) , y ) is continuous on Ω × Ω . Further, the function s μ 1 | ln s | α β 1 has an integrable singularity, and that is why the function P α β , μ ( t k ( x ) , y ) is continuous on Ω × Ω . Then by virtue of the uniform convergence of the series k = 1 a k ( x ) , the kernel K ( x , y ) is also a continuous function on Ω × Ω .

Hence, one can apply Fredholm theory to equation (26). Since in the case of f ( x ) = 0 and fulfillment of the condition (20), the solution of problem (17)-(18) can only be v ( x ) 0 , for f ( x ) = 0 the integral equation (26) has only a trivial solution.

Hence, for any f ( x ) C ( Ω ) the solution of equation (26) exists, is unique, and belongs to the class C ( Ω ) . Using this solution, we construct the function v ( x ) which will satisfy all the conditions of problem (17)-(18).

If the condition (21) is valid, then v ( x ) = C satisfies the condition (18) at f ( x ) 0 , i.e. the corresponding homogeneous equation (26) has the nonzero solution v ( x ) C . Then the adjoint homogeneous equation has also a nonzero solution, and that is why in this case fulfillment of the condition (24) is necessary and sufficient for solvability of problem (17)-(18). The lemma is proved. □

### 6 Study of the basic problem

We now formulate the basic statement.

Theorem 1Let μ > 0 , 0 β α , α + β 0 , Γ k Ω ¯ 1 Ω , a k ( x ) , k = 1 , 2 ,  , be continuous functions satisfying the condition (19). Then:

(1) If the condition (20) is fulfilled, then problem (8)-(9) is uniquely solvable at any f ( x ) C ( Ω ) .

(2) If the condition (21) is fulfilled, then the condition (24) is necessary and sufficient for solvability of problem (8)-(9) where the function ψ 0 ( x ) is a solution of equation (23). If a solution of the problem exists, then it is unique up to the constant summand.

(3) If a solution of problem (8)-(9) exists, then it is represented in the form of u ( x ) = J μ α [ v ] ( x ) , where v ( x ) is a solution of problem (17)-(18).

Proof (1) Let a solution u ( x ) of problem (8)-(9) exist. Apply to this function the operator  D μ α and denote v ( x ) = D μ α [ u ] ( x ) . Take the problem which the function v ( x ) satisfies. Since by Lemma 2, in the case of harmonicity of the function u ( x ) , the function D μ α [ u ] ( x ) is also harmonic in Ω, and the function v ( x ) is harmonic.

Further, since according to Lemma 5 the equality D μ β [ u ] ( x ) = J μ α β [ D μ α [ u ] ] ( x ) holds, the boundary condition of problem (8)-(9),

D μ α [ u ] ( x ) k = 1 a k ( x ) D μ β [ u ] ( t k ( x ) ) = f ( x ) ,

with respect to the function v ( x ) will be rewritten in the form of

v ( x ) k = 1 a k ( x ) J μ α β [ v ] ( t k ( x ) ) = f ( x ) , x Ω .

In addition, since D μ α [ u ] ( x ) C ( Ω ¯ ) , we have v ( x ) C ( Ω ¯ ) . Thus, if u ( x ) is a solution of problem (8)-(9), then the function v ( x ) = D μ α [ u ] ( x ) will be a solution of problem (17)-(18).

Now, let the conditions (19) and (20) be realized. Then by Lemmas 6 and 7, for any f ( x ) C ( Ω ) the solution of problem (17)-(18) exists, is unique, and designate u ( x ) = J μ α [ v ] ( x ) , x Ω ¯ . Then we have by Lemma 5 D μ α [ u ] ( x ) = D μ α [ J μ α [ v ] ] ( x ) = v ( x ) in Ω, and therefore we get D μ α [ u ] ( x ) C ( Ω ¯ ) . Harmonicity of the function u ( x ) follows from Lemma 2, and fulfillment of the conditions (9) can be checked immediately:

D μ α [ u ] ( x ) k = 1 a k ( x ) D μ β [ u ] ( t k ( x ) ) = v ( x ) k = 1 a k ( x ) J μ α β [ D μ α [ u ] ] ( t k ( x ) ) = f ( x ) , x Ω .

The first statement of the theorem is proved.

(2) Let now the condition (21) be fulfilled, and let the solution u ( x ) of problem (8)-(9) exist. Consider the function v ( x ) = D μ α [ u ] ( x ) . As in the first case, we show that the function v ( x ) satisfies the conditions of problem (17)-(18). Then according to Lemma 7, fulfillment of the condition (24) is necessary. Thus, we prove the necessity of the condition (24) at fulfillment of the equality (21).

We show that if the equality (21) is fulfilled, then the condition (24) is also sufficient for the existence of the solution of problem (8)-(9).

In fact, if the conditions (21) and (24) are realized, a solution of problem (17)-(18) exists, is unique up to constant summand, and v ( x ) C ( Ω ¯ ) . Then, similarly to the proof of the first statement of the theorem, the function u ( x ) = J μ α [ v ] ( x ) satisfies all the conditions of problem (8)-(9). The theorem is proved. □

Remark 1 One can show that in the case of k = 1 | a k ( x ) | > μ α β , the corresponding homogeneous problem (8)-(9) has nontrivial solutions.

Example 1 Let a 1 0 , a k = 0 , k = 2 , 3 ,  , and t 1 ( x ) = δ x , 0 < δ < 1 . Further, let H k ( x ) be a homogeneous harmonic polynomial of the power k. By virtue of the equality (11), we have D μ α [ H k ] ( x ) = ( k + μ ) α H k ( x ) .

Then

H k ( t 1 ( x ) ) = H k ( δ x ) = δ k H k ( x ) , D μ β H k ( δ x ) = δ k ( k + μ ) β H k ( x )

and

D μ α [ H k ] ( x ) a 1 D μ β H k ( t 1 ( x ) ) = ( k + μ ) α H k ( x ) a 1 δ k ( k + μ ) β H k ( x ) = [ ( k + μ ) α a 1 δ k ( k + μ ) β ] H k ( x ) .

Hence, for a 1 = δ k ( k + μ ) α β the harmonic polynomial H k ( x ) will be the solution of the homogeneous problem (8)-(9). If δ is a number close to zero, then we have a 1 > 1 .

If the dimension of the space n = 3 , then the number of these polynomials is equal to 2 k + 1 [47].

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

All authors completed the paper together. All authors read and approved the final manuscript.

### Acknowledgements

This work has been supported by the MON Republic of Kazakhstan under Research Grant No. 0713/GF.

### References

1. Serbina, LI: A model of substance transfer in fractal media. Math. Model.. 15, 17–28 (2003)

2. Koeller, RC: Application of fractional calculus to the theory of viscoelasticity. J. Appl. Mech.. 51, 299–307 (1984). Publisher Full Text

3. Heymans, N, Podlubny, I: Physical interpretation of initial condition for fractional differential equations with Riemann-Liouville fractional derivatives. Rheol. Acta. 45, 765–771 (2006). Publisher Full Text

4. Arena, P, Caponetto, R, Fortuna, L, Porto, D: Nonlinear Noninteger Order Circuits and Systems: an Introduction, World Scientific, River Edge (2000)

5. Bonilla, B, Rivero, M, Rodrigues-Germa, L, Trujillo, JJ: Fractional differential equations as alternative models to nonlinear differential equations. Appl. Math. Comput.. 187, 79–88 (2007). Publisher Full Text

6. Amanov, D, Ashyralyev, A: Initial-boundary value problem for fractional partial differential equations of higher order. Abstr. Appl. Anal. (2012). Publisher Full Text

7. Ashyralyev, A: A note on fractional derivatives and fractional powers of operators. J. Math. Anal. Appl.. 357(1), 232–236 (2009). Publisher Full Text

8. Ashyralyev, A: Well-posedness of fractional parabolic equations. Bound. Value Probl. (2013). BioMed Central Full Text

9. Ashyralyev, A, Dal, F: Finite difference and iteration methods for fractional hyperbolic partial differential equations with the Neumann condition. Discrete Dyn. Nat. Soc. (2012). Publisher Full Text

10. Ashyralyev, A, Dal, F, Pinar, Z: A note on the fractional hyperbolic differential and difference equations. Appl. Math. Comput.. 217(9), 4654–4664 (2011). Publisher Full Text

11. Ashyralyev, A, Hicdurmaz, B: A note on the fractional Schrodinger differential equations. Kybernetes. 40(5-6), 736–750 (2011)

12. Ashyralyev, A, Cakir, Z: FDM for fractional parabolic equations with the Neumann condition. Adv. Differ. Equ. (2013). BioMed Central Full Text

13. Yakar, A, Koksal, ME: Existence results for solutions of nonlinear fractional differential equations source. Abstr. Appl. Anal. (2012). Publisher Full Text

14. Tikhonov, AN: On boundary conditions containing derivatives of order higher than the order of the equation. Mat. Sb.. 26(1), 35–56 (1950) (Russian)

15. Amanov, D, Yuldasheva, AV: Solvability and spectral properties of boundary value problems for equations of even order. Malays. J. Math. Sci.. 3(2), 227–248 (2009)

16. Amanov, D, Murzambetova, MB: A boundary value problem for a fourth order partial differential equation with the lowest term. Vestn. Udmurtsk. Univ. Mat. Mekh. Komp. Nauki. 2013(1), 3–10 (2013) (Russian)

17. Amanov, D: Boundary value problem for the higher order equation with fractional derivative. Malays. J. Math. Sci.. 7(2), 203–217 (2013)

18. Bavrin, II: Operators for harmonic functions and their applications. Differ. Uravn.. 21, 9–15 (1985) (Russian)

19. Sokolovskii, VB: On a generalization of the Neumann problem. Differ. Uravn.. 24, 714–716 (1988) (Russian)

20. Bitsadze, AV: On the Neumann problem for harmonic functions. Dokl. Akad. Nauk SSSR. 311, 11–13 (1990) (Russian)

21. Berdyshev, AS, Turmetov, BK, Kadirkulov, BJ: Some properties and applications of the integrodifferential operators of Hadamard-Marchaud type in the class of harmonic functions. Sib. Math. J.. 53(4), 600–610 (2012) (Russian)

(Russian)

Publisher Full Text

22. Kirane, M, Tatar, NE: Nonexistence for the Laplace equation with a dynamical boundary condition of fractional type. Sib. Math. J.. 48(5), 1056–1064 (2007) (Russian)

23. Torebek, BT, Turmetov, BK: On solvability of a boundary value problem for the Poisson equation with the boundary operator of a fractional order. Bound. Value Probl.. 2013, (2013) Article ID 93

Article ID 93

BioMed Central Full Text

24. Turmetov, BK: On a boundary value problem for a harmonic equation. Difer. Uravn.. 32(8), 1089–1092 (1996) (Russian)

25. Turmetov, BK: On smoothness of a solution to a boundary value problem with fractional order boundary operator. Mat. Tr.. 7(1), 189–199 (2004) (Russian)

26. Umarov, SR, Luchko, YF, Gorenflo, R: On boundary value problems for elliptic equations with boundary operators of fractional order. Fract. Calc. Appl. Anal.. 3(4), 454–468 (2000)

27. Samko, SG, Kilbas, AA, Marichev, OI: Fractional Integrals and Derivatives, Theory and Applications, Gordon & Breach, Amsterdam (1993)

28. Kilbas, AA, Tityura, AA: Fractional derivative of the Marsho-Hadamard type and inversion of fractional integrals. Dokl. Nat. Akad. Nauk Belarusi. 50(4), 5–10 (2006) (Russian)

29. Bitsadze, AV, Samarskii, AA: Some elementary generalizations of linear elliptic boundary value problems. Dokl. Akad. Nauk SSSR. 185, 739–740 (1969) (Russian)

30. Pulatov, AK: A Bitsadze-Samarskii problem. Differ. Uravn.. 25, 537–540 (1989) (Russian)

31. Skubachevskii, AL: Nonclassical boundary value problems. I. Sovrem. Mat. Fundam. Napravl.. 26, 3–132 (2007) (Russian)

32. Skubachevskii, AL: Nonclassical boundary value problems. II. Sovrem. Mat. Fundam. Napravl.. 33, 3–179 (2009) (Russian)

33. Ashyralyev, A: On well-posedness of the nonlocal boundary value problem for elliptic equations. Numer. Funct. Anal. Optim.. 24(1-2), 1–15 (2003). Publisher Full Text

34. Ashyralyev, A: A note on the Bitsadze-Samarskii type nonlocal boundary value problem in a Banach space. J. Math. Anal. Appl.. 344(1), 557–573 (2008). Publisher Full Text

35. Ashyralyev, A, Sharifov, YA: Existence and uniqueness of solutions of the system of nonlinear fractional differential equations with nonlocal and integral boundary conditions. Abstr. Appl. Anal. (2012). Publisher Full Text

36. Ashyralyev, A, Ozger, F: The hyperbolic-elliptic equation with the nonlocal condition. Math. Methods Appl. Sci. (2013). Publisher Full Text

37. Ashyralyev, A, Ozturk, E: The numerical solution of the Bitsadze-Samarskii nonlocal boundary value problems with the Dirichlet-Neumann condition. Abstr. Appl. Anal. (2012). Publisher Full Text

38. Ashyralyev, A, Ozturk, E: On Bitsadze-Samarskii type nonlocal boundary value problems for elliptic differential and difference equations: well-posedness. Appl. Math. Comput.. 219(3), 1093–1107 (2012). Publisher Full Text

39. Ashyralyev, A, Tetikoglu, FS: A note on Bitsadze-Samarskii type nonlocal boundary value problems: well-posedness. Numer. Funct. Anal. Optim.. 34(9), 939–975 (2013). Publisher Full Text

40. Berdyshev, AS, Kadirkulov, BJ, Nieto, JJ: Solvability of an elliptic partial differential equation with boundary condition involving fractional derivatives. Complex Var. Elliptic Equ. (2013). Publisher Full Text

41. Karachik, VV, Turmetov, BK, Torebek, BT: On some integro-differential operators in the class of harmonic functions and their applications. Sib. Adv. Math.. 22(2), 115–134 (2012). Publisher Full Text

42. Mahmudov, NI, Zorlu, S: Approximate controllability of fractional integro-differential equations involving nonlocal initial conditions. Bound. Value Probl. (2013). BioMed Central Full Text

43. Qu, H, Liu, X: Existence of nonnegative solutions for a fractional m-point boundary value problem at resonance. Bound. Value Probl.. 2013, (2013) Article ID 127

Article ID 127

BioMed Central Full Text

44. Hilmi, E, Kilicman, A: Nonlocal boundary value problems for impulsive fractional integro-differential equations in Banach spaces. Bound. Value Probl.. 2012, (2012) Article ID 145

Article ID 145

BioMed Central Full Text

45. Stein, EM, Weiss, G: Introduction to Fourier Analysis on Euclidean Spaces, Princeton University Press, Princeton (1971)

46. Vladimirov, VS: Uravneniya mathematicheskoi fiziki, Nauka, Moscow (1981)

47. Bitsadze, AV: Uravneniya mathematicheskoi fiziki, Nauka, Moscow (1982)