Skip to main content

Well-posedness of fractional parabolic equations

Abstract

In the present paper, we consider the abstract Cauchy problem for the fractional differential equation

d u ( t ) d t + D t 1 2 u(t)+A(t)u(t)=f(t),0<t<1,u(0)=0
(1)

in an arbitrary Banach space E with the strongly positive operators A(t). The well-posedness of this problem in spaces of smooth functions is established. The coercive stability estimates for the solution of problems for 2m th order multidimensional fractional parabolic equations and one-dimensional fractional parabolic equations with nonlocal boundary conditions in a space variable are obtained. The stable difference scheme for the approximate solution of this problem is presented. The well-posedness of the difference scheme in difference analogues of spaces of smooth functions is established. In practice, the coercive stability estimates for the solution of difference schemes for the fractional parabolic equation with nonlocal boundary conditions in a space variable and the 2m th order multidimensional fractional parabolic equation are obtained.

MSC:65M12, 65N12.

1 Introduction

It is known that differential equations involving derivatives of noninteger order have shown to be adequate models for various physical phenomena in areas like rheology, damping laws, diffusion processes, etc. Methods of solutions of problems for fractional differential equations have been studied extensively by many researchers (see, e.g., [1–43] and the references given therein).

The role played by coercive stability inequalities (well-posedness) in the study of boundary value problems for parabolic partial differential equations is well known (see, e.g., [44–51]). In the present paper, the initial value problem

d u ( t ) d t + D t 1 2 u(t)+A(t)u(t)=f(t),0<t<1,u(0)=0
(2)

for the fractional differential equation in an arbitrary Banach space E with the linear (unbounded) operators A(t) is considered. Here u(t) and f(t) are the unknown and the given functions, respectively, defined on [0,T] with values in E. The derivative u ′ (t) is understood as the limit in the norm of E of the corresponding ratio of differences. A(t) is a given closed linear operator in E with the domain D(A(t))=D, independent of t and dense in E. Finally, u(0)=0.

Here D t 1 2 = D 0 + 1 2 is the standard Riemann-Liouville derivative of order 1 2 . This fractional differential equation corresponds to the Basset problem [9]. It represents a classical problem in fluid dynamics where the unsteady motion of a particle accelerates in a viscous fluid due to the gravity of force. Recently, fractional Basset equations with independent in t operator coefficients A(t)=A have been studied extensively (see, e.g., [52–56] and the references given therein).

In the present paper, the well-posedness of problem (2) with dependent in t operator coefficients A(t) in spaces of smooth functions is established. In practice, the coercive stability estimates for the solution of problems for 2m th order multidimensional fractional parabolic equations and one-dimensional fractional parabolic equations with nonlocal boundary conditions in a space variable are obtained. The stable difference scheme for the approximate solution of initial value problem (2)

{ τ − 1 ( u k − u k − 1 ) + A k u k + 1 π ∑ m = 1 k Γ ( k − m + 1 2 ) ( k − m ) ! u m − u m − 1 τ 1 2 = f k , f k = f ( t k ) , A k = A ( t k ) , t k = k τ , 1 ≤ k ≤ N , N τ = 1 , u 0 = 0
(3)

is presented. Here Γ(k−m+ 1 2 )= ∫ 0 ∞ t k − m − 1 2 e − t dt.

The paper is organized as follows. The well-posedness of problem (2) in spaces of smooth functions is established in Section 2. In Section 3 the coercive stability estimates for the solution of problems for 2m th order multidimensional fractional parabolic equations and one-dimensional fractional parabolic equations with nonlocal boundary conditions are obtained. The well-posedness of (3) in difference analogues of spaces of smooth functions is established and the coercive stability estimates for the solution of difference schemes for the fractional parabolic equation with nonlocal boundary conditions in a space variable and the 2m th order multidimensional fractional parabolic equation are obtained in Section 4.

2 The well-posedness of problem (2)

A function u(t) is called a solution of problem (2) if the following conditions are satisfied:

  1. (i)

    u(t) is continuously differentiable on the segment [0,1]. The derivatives at the endpoints of the segment are understood as the appropriate unilateral derivatives.

  2. (ii)

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

  3. (iii)

    u(t) satisfies the equation and the initial condition (2).

A solution of problem (2) defined in this manner will from now on be referred to as a solution of problem (2) in the space C(E)=C([0,1],E) of all continuous functions φ(t) defined on [0,1] with values in E equipped with the norm

∥ φ ∥ C ( E ) = max 0 ≤ t ≤ 1 ∥ φ ( t ) ∥ E .
(4)

In this paper, positive constants, which can differ in time, are indicated with an M. On the other hand, M(α,β,…) is used to focus on the fact that the constant depends only on α,β,… .

The well-posedness in C(E) of boundary value problem (2) means that the coercive inequality

∥ u ′ ∥ C ( E ) + ∥ A ( ⋅ ) u ∥ C ( E ) ≤M ∥ f ∥ C ( E )
(5)

is true for its solution u(t)∈C(E).

Suppose that for each t∈[0,1] the operator −A(t) generates an analytic semigroup exp{−sA(t)} (s≥0) with an exponentially decreasing norm, when s→+∞, i.e., the following estimates

∥ exp ( − s A ( t ) ) ∥ E → E , ∥ s A ( t ) exp ( − s A ( t ) ) ∥ E → E ≤M e − δ s (s>0)
(6)

hold for some M∈[1,+∞), δ∈(0,+∞). From this inequality it follows the operator A − 1 (t) exists and is bounded, and hence A(t) is closed in C(E).

Suppose that the operator A(t) A − 1 (s) is Hölder continuous in t in the uniform operator topology for each fixed s, that is,

∥ [ A ( t ) − A ( τ ) ] A − 1 ( s ) ∥ E → E ≤M | t − τ | ε ,0<ε≤1,0≤t,s,τ≤1.
(7)

An operator-valued function v(t,s), defined and strongly continuous jointly in t, s for 0≤s<t≤1, is called a fundamental solution of (2) if

  1. (1)

    the operator v(t,s) is strongly continuous in t and s for 0≤s<t≤T,

  2. (2)

    the following identity holds:

v(t,s)=v(t,Ï„)v(Ï„,s),v(t,t)=Ifor 0≤s≤τ≤t≤1,
  1. (3)

    the operator v(t,s) maps the region D into itself. The operator u(t,s)=A(t)v(t,s) A − 1 (s) is bounded and strongly continuous in t and s for 0≤s<t≤1,

  2. (4)

    on the region D the operator v(t,s) is strongly differentiable relative to t and s, while

    ∂ v ( t , s ) ∂ t =−A(t)v(t,s)
    (8)

and

∂ v ( t , s ) ∂ s =v(t,s)A(s).
(9)

Now, let us obtain the representation for the solution of problem (2). The initial value problem

d u d t +A(t)u(t)=F(t),0<t<1,u(0)= u 0
(10)

has a unique solution [54] and the following formula holds:

u(t)=v(t,0) u 0 + ∫ 0 t v(t,s)F(s)ds.
(11)

Using u(0)=0 and the formula F(s)=f(s)− D s 1 2 u(s), we get

u(t)=− ∫ 0 t v(t,s) D s 1 2 u(s)ds+ ∫ 0 t v(t,s)f(s)ds.
(12)

Now, we will give a series of interesting lemmas and estimates concerning the fundamental solution v(t,s) of (2) which will be useful in the sequel.

Lemma 2.1 For any 0≤s<t≤1 and u∈D, the following identities hold:

(13)
(14)
(15)
(16)

Lemma 2.2 For any 0≤s<t≤t+r≤1, 0≤α≤1 and 0≤ε≤1, the following estimates hold:

(17)
(18)
(19)
(20)
(21)

Theorem 2.1 Let A(t) be a strongly positive operator in a Banach space E and f(t)∈C(E). Then for the solution u(t) in C(E) of initial value problem (2), the following stability inequality holds:

∥ D t 1 2 u ∥ C ( E ) + ∥ u ′ + A ( ⋅ ) u ∥ C ( E ) ≤M ∥ f ∥ C ( E ) .
(22)

Proof Using formula (12), we get

u ′ (t)=− D t 1 2 u(t)+f(t)+ ∫ 0 t A(t)v(t,s) D s 1 2 u(s)ds− ∫ 0 t A(t)v(t,s)f(s)ds.
(23)

Applying formula (23) and the formula

D t 1 2 u(t)= ∫ 0 t u ′ ( p ) d p π ( t − p ) 1 2 ,
(24)

we obtain

D t 1 2 u ( t ) = ∫ 0 t 1 π ( t − s ) 1 2 ( − D s 1 2 u ( s ) + f ( s ) ) d s + ∫ 0 t ∫ s t 1 π ( t − p ) 1 2 A ( p ) v ( p , s ) d p D s 1 2 u ( s ) d s − ∫ 0 t ∫ s t 1 π ( t − p ) 1 2 A ( p ) v ( p , s ) d p f ( s ) d s .
(25)

Let us first obtain the estimate

∥ ∫ s t 1 π ( t − p ) 1 2 A ( p ) v ( p , s ) d p ∥ E → E ≤ M t − s
(26)

for any 0≤s<t≤1. We have that

∫ s t 1 π ( t − p ) 1 2 A(p)v(p,s)dp= ∫ t + s 2 t 1 π ( t − p ) 1 2 A(p)v(p,s)dp
(27)
+ ∫ s t + s 2 1 π ( t − p ) 1 2 A(p)v(p,s)dp= J 1 + J 2 .
(28)

Applying estimate (21), we get

∥ J 1 ∥ E → E ≤M ∫ t + s 2 t 1 π ( t − p ) 1 2 1 p − s dp≤ 2 M t − s ∫ t + s 2 t 1 π ( t − p ) 1 2 dp= M 1 t − s .
(29)

Now, we will estimate J 2 . We have that

J 2 = 1 π t − s I−v ( t + s 2 , s ) 2 π t − s + ∫ s t + s 2 1 2 π ( t − p ) 3 2 v(p,s)dp.
(30)

Applying estimate (17), we get

∥ J 2 ∥ E → E ≤ 1 t − s +M 2 t − s +M ∫ s t + s 2 1 2 π ( t − p ) 3 2 dp≤ M 2 t − s .
(31)

Estimate (26) follows from estimates (29) and (31).

Now, let us first estimate z(t)= ∥ D t 1 2 u ( t ) ∥ E . Applying the triangle inequality and estimate (26), we get

z(t)≤M ∫ 0 t 1 t − s z(s)ds+M ∫ 0 t 1 t − s ∥ f ( s ) ∥ E ds
(32)

for any t∈[0,1]. Applying the above inequality and the integral inequality, we obtain

max 0 ≤ t ≤ 1 z(t)≤M max 0 ≤ t ≤ 1 ∥ f ( t ) ∥ E .
(33)

Using the triangle inequality and equation (2), we get

max 0 ≤ t ≤ 1 ∥ u t + A ( t ) u ( t ) ∥ E ≤ [ max 0 ≤ t ≤ 1 ∥ f ( t ) ∥ E + max 0 ≤ t ≤ 1 ∥ D t 1 2 u ( t ) ∥ E ]
(34)
≤ M 1 max 0 ≤ t ≤ 1 ∥ f ( t ) ∥ E .
(35)

Estimate (22) follows from estimates (33) and (35). Theorem 2.1 is proved. □

With the help of A(t), we introduce the fractional spaces E α (E,A(t)), 0<α<1, consisting of all v∈E for which the following norms are finite:

∥ v ∥ E α = sup z > 0 z 1 − α ∥ A ( t ) exp { − z A ( t ) } v ∥ E .
(36)

From (6) and (7) it follows that

Theorem 2.2 E α (E,A(t))= E α (E,A(0)) for all 0<α<1 and 0≤t≤1.

Problem (2) is not well posed in C(E) for arbitrary E. It turns out that a Banach space E can be restricted to a Banach space E ′ in such a manner that the restricted problem (2) in E ′ will be well posed in C( E ′ ). The role of E ′ will be played here by the fractional spaces E α = E α (A(t),E) (0<α<1).

Theorem 2.3 Suppose f(t)∈C( E α ) (0<α<1). Suppose that assumptions (6) and (7) hold and 0<α≤ε<1. Then for the solution u(t) in C( E α ) of problem (2), the coercive inequality

∥ u ′ ∥ C ( E α ) + ∥ A ( ⋅ ) u ∥ C ( E α ) ≤M α − 1 ( 1 − α ) − 1 ∥ f ∥ C ( E α )
(37)

holds.

Proof

By Theorem 2.1,

∥ D t 1 2 u ∥ C ( E α ) ≤M ∥ f ∥ C ( E α )
(38)

for the solution of initial value problem (2). The proof of the estimate

∥ A ( ⋅ ) u ∥ C ( E α ) ≤M α − 1 ( 1 − α ) − 1 ∥ f ∥ C ( E α )
(39)

for the solution of initial value problem (2) is based on formula (12), estimate (38) and the following estimates [54]:

(40)
(41)

Using equation (2) and the triangle inequality, we get

max 0 ≤ t ≤ 1 ∥ u ′ ( t ) ∥ E α ≤ [ max 0 ≤ t ≤ 1 ∥ f ( t ) ∥ E α + max 0 ≤ t ≤ 1 ∥ A ( t ) u ( t ) ∥ E α + max 0 ≤ t ≤ 1 ∥ D t 1 2 u ( t ) ∥ E α ] ≤ M 1 α − 1 ( 1 − α ) − 1 max 0 ≤ t ≤ 1 ∥ f ( t ) ∥ E α .
(42)

Estimate (37) follows from estimates (39) and (42). Theorem 2.3 is proved. □

Let us give, without proof, the following result.

Theorem 2.4 Suppose that assumption (6) holds. Suppose that the operator A(t) A − 1 (s) is Hölder continuous in t in the uniform operator topology for each fixed s, that is,

∥ [ A ( t ) − A ( τ ) ] A − 1 ( s ) ∥ E α → E α ≤M | t − τ | ε ,0<ε≤1,
(43)

where M and ε are positive constants independent of t, s and τ for 0≤t,s,τ≤T. Suppose f(t)∈C( E α ) (0<α<1). Then for the solution u(t) in C( E α ) of problem (2), the coercive inequality

∥ u ′ ∥ C ( E α ) + ∥ A ( ⋅ ) u ∥ C ( E α ) ≤M α − 1 ( 1 − α ) − 1 ∥ f ∥ C ( E α )
(44)

holds.

3 Applications

Now, we consider the applications of Theorems 2.1, 2.3 and 2.4.

First, the Cauchy problem on the range {0≤t≤1,x∈ R n } for the 2m-order multidimensional fractional parabolic equation is considered:

{ ∂ v ( t , x ) ∂ t + D t 1 2 v ( t , x ) + ∑ | r | = 2 m a r ( t , x ) ∂ | r | v ( t , x ) ∂ x 1 r 1 ⋯ ∂ x n r n + σ v ( t , x ) = f ( t , x ) , 0 < t < 1 , v ( 0 , x ) = 0 , x ∈ R n , | r | = r 1 + ⋯ + r n ,
(45)

where a r (t,x) and f(t,x) are given as sufficiently smooth functions. Here, σ is a sufficiently large positive constant.

Let us consider a differential operator with constant coefficients of the form

B= ∑ | r | = 2 m b r ∂ r 1 + ⋯ + r n ∂ x 1 r 1 ⋯ ∂ x n r n ,
(46)

acting on functions defined on the entire space R n . Here r∈ R n is a vector with nonnegative integer components, |r|= r 1 +⋯+ r n . If φ(y) (y=( y 1 ,…, y n )∈ R n ) is an infinitely differentiable function that decays at infinity together with all its derivatives, then by means of the Fourier transformation, one establishes the equality

F( B φ )(ξ)=B(ξ)F(φ)(ξ).
(47)

Here the Fourier transform operator is defined by the following rule:

(48)
(49)

The function B(ξ) is called the symbol of the operator B and is given by

B(ξ)= ∑ | r | = 2 m b r ( i ξ 1 ) r 1 ⋯ ( i ξ n ) r n .
(50)

We will assume that the symbol

B t , x (ξ)= ∑ | r | = 2 m a r (t,x) ( i ξ 1 ) r 1 ⋯ ( i ξ n ) r n ,ξ=( ξ 1 ,…, ξ n )∈ R n
(51)

of the differential operator of the form

B t , x = ∑ | r | = 2 m a r (t,x) ∂ | r | ∂ x 1 r 1 ⋯ ∂ x n r n
(52)

acting on functions defined on the space R n , satisfies the inequalities

0< M 1 | ξ | 2 m ≤ ( − 1 ) m B t , x (ξ)≤ M 2 | ξ | 2 m <∞
(53)

for ξ≠0. Problem (45) has a unique smooth solution. This allows us to reduce problem (45) to the abstract Cauchy problem (2) in a Banach space E= C μ ( R n ) of all continuous bounded functions defined on R n satisfying the Hölder condition with the indicator μ∈(0,1) with a strongly positive operator A t , x = B t , x +δI defined by (52) (see [57, 58]).

Theorem 3.1 For the solution of boundary problem (45), the following estimates are satisfied:

(54)
(55)

The proof of Theorem 3.1 is based on the abstract Theorems 2.1, 2.3, 2.4 and the coercivity inequality for an elliptic operator A t , x in C μ ( R n ) and on the following theorem on the structure of the fractional spaces E α ( C μ ( R n ), A t , x ).

Theorem 3.2 E α ( C μ ( R n ), A t , x )= C 2 m α + μ ( R n ) for all 0<2mα+μ<1 and 0≤t≤1 [59].

Second, we consider the mixed boundary value problem for the fractional parabolic equation

{ ∂ v ( t , x ) ∂ t + D t 1 2 v ( t , x ) − a ( t , x ) ∂ 2 v ( t , x ) ∂ x 2 + σ v ( t , x ) = f ( t , x ) , 0 < t < 1 , 0 < x < 1 , v ( 0 , x ) = 0 , 0 ≤ x ≤ 1 , u ( t , 0 ) = u ( t , 1 ) , u x ( t , 0 ) = u x ( t , 1 ) , 0 ≤ t ≤ 1 ,
(56)

where a(t,x) and f(t,x) are given sufficiently smooth functions and a(t,x)≥a>0. Here, σ is a sufficiently large positive constant.

We introduce the Banach spaces C β [0,1] (0<β<1) of all continuous functions φ(x) satisfying the Hölder condition for which the following norms are finite:

∥ φ ∥ C β [ 0 , 1 ] = ∥ φ ∥ C [ 0 , 1 ] + sup 0 ≤ x < x + τ ≤ 1 | φ ( x + τ ) − φ ( x ) | τ β ,
(57)

where C[0,1] is the space of all continuous functions φ(x) defined on [0,1] with the usual norm

∥ φ ∥ C [ 0 , 1 ] = max 0 ≤ x ≤ 1 | φ ( x ) | .
(58)

It is known that the differential expression [60]

A t , x v=−a(t,x) v ′ ′ (t,x)+σv(t,x)
(59)

defines a positive operator A t , x acting in C β [0,1] with the domain C β + 2 [0,1] and satisfying the conditions v(t,0)=v(t,1), v x (t,0)= v x (t,1). Therefore, we can replace the mixed problem (56) by the abstract boundary value problem (2). Using the results of Theorems 2.1, 2.3, 2.4, we can obtain the following theorem.

Theorem 3.3 For the solution of mixed problem (56), the following estimates are valid:

(60)
(61)

The proof of Theorem 3.3 is based on abstract Theorems 2.1, 2.3, 2.4 and on the following theorem on the structure of the fractional spaces E α (C[0,1], A t , x ).

Theorem 3.4 E α (C[0,1], A t , x )= C 2 α [0,1] for all 0<α< 1 2 , 0≤t≤1 [60].

4 The well-posedness of problem (3)

Let us first obtain the representation for the solution of problem (3). It is clear that the first order of accuracy difference scheme

τ − 1 ( u k − u k − 1 )+ A k u k = F k ,1≤k≤N,Nτ=1, u 0 =0
(62)

has a solution and the following formula holds:

u k = ∑ s = 1 k u τ (k,s) F s τ,1≤k≤N,
(63)

where

u τ (k,j)={ R k ⋯ R j + 1 , k > j , I , k = j .
(64)

Here R k = ( I + τ A k ) − 1 . Denote that

D τ 1 2 u k = 1 π ∑ m = 1 k Γ ( k − m + 1 2 ) ( k − m ) ! u m − u m − 1 τ 1 2 .
(65)

Applying the formula F k = f k − D τ 1 2 u k , we get

u k =− ∑ s = 1 k u τ (k,s) D τ 1 2 u s τ+ ∑ s = 1 k u τ (k,s) f s τ,1≤k≤N.
(66)

So, formula (66) gives the representation for the solution of problem (3).

Let F τ (E) be the linear space of mesh functions φ τ = { φ k } 1 N with values in the Banach space E. Next on F τ (E) we introduce the Banach space C τ (E)=C( [ 0 , 1 ] τ ,E) with the norm

∥ φ τ ∥ C τ ( E ) = max 1 ≤ k ≤ N ∥ φ k ∥ E .
(67)

Theorem 4.1 Let A(t) be a strongly positive operator in a Banach space E. Then for the solution u Ï„ = { u k } 1 N in C Ï„ (E) of initial value problem (3), the stability inequality

∥ { D τ 1 2 u k } 1 N ∥ C τ ( E ) + ∥ { τ − 1 ( u k − u k − 1 ) + A k u k } 1 N ∥ C τ ( E ) ≤M ∥ f τ ∥ C τ ( E )
(68)

holds.

Proof Using formula (66), we get

τ − 1 ( u k − u k − 1 )=− D τ 1 2 u k + ∑ s = 1 k A k u τ (k,s) D τ 1 2 u s τ+ f k − ∑ s = 1 k A k u τ (k,s) f s τ.
(69)

Applying formulas (69) and (65), we obtain

D τ 1 2 u k = 1 π ∑ m = 1 k Γ ( k − m + 1 2 ) ( k − m ) ! τ 1 2 [ − D τ 1 2 u m + f m ] + 1 π ∑ m = 1 k Γ ( k − m + 1 2 ) ( k − m ) ! [ ∑ s = 1 m A k u τ ( k , s ) D τ 1 2 u s τ 3 2 − ∑ s = 1 m A k u τ ( k , s ) f s τ 3 2 ] = 1 π ∑ m = 1 k Γ ( k − m + 1 2 ) ( k − m ) ! τ 1 2 [ − D τ 1 2 u m + f m ] + 1 π ∑ s = 1 k ∑ m = s k Γ ( k − m + 1 2 ) ( k − m ) ! A m u τ ( m , s ) D τ 1 2 u s τ 3 2 − 1 π ∑ s = 1 k ∑ m = s k Γ ( k − m + 1 2 ) ( k − m ) ! A m u τ ( m , s ) f s τ 3 2 .
(70)

Let us first obtain the estimate

∥ 1 π ∑ m = s k Γ ( k − m + 1 2 ) ( k − m ) ! A m u τ ( m , s ) τ 1 2 ∥ E → E ≤ M ( k − s ) τ
(71)

for any 1≤s<k≤N. We have that

(72)

Using estimates

∥ A k u τ ( k , s ) ∥ E → E ≤ M ( k − s + 1 ) τ , ∥ u τ ( k , s ) ∥ E → E ≤M,1≤k≤N
(73)

and the following elementary inequality:

Γ ( k − m + 1 2 ) ( k − m ) ! ≤ 1 k − m ,0≤m<k,
(74)

we obtain

∥ J 1 ∥ E → E ≤ 1 π ∑ m = [ s + k 2 ] k Γ ( k − m + 1 2 ) ( k − m ) ! ∥ A m u τ ( m , s ) ∥ E → E τ 1 2
(75)
≤M 1 π ∑ m = [ s + k 2 ] k Γ ( k − m + 1 2 ) ( k − m ) ! 1 ( m − s + 1 ) τ τ 1 2
(76)
≤ 2 M ( k − s ) τ 1 π ∑ m = [ s + k 2 ] k τ ( k − m ) τ ≤ M 1 ( k − s ) τ .
(77)

Now, we will estimate J 2 . We have that

J 2 = 1 π Γ ( k − s + 1 2 ) ( k − s ) ! τ − 1 2 − 1 π Γ ( k − [ s + k 2 ] + 3 2 ) ( k − [ s + k 2 ] + 1 ) ! u τ ( [ s + k 2 ] , s ) τ − 1 2
(78)
+ 1 π ∑ m = s + 1 [ s + k 2 ] − 1 [ Γ ( k − m + 1 2 ) ( k − m ) ! − Γ ( k − m + 3 2 ) ( k − m + 1 ) ! ] u τ (m−1,s) τ − 1 2 .
(79)

Applying estimates (73) and (74), we get

∥ J 2 ∥ E → E ≤ 1 π 1 ( k − s ) τ + 1 π ∥ u τ ( [ s + k 2 ] , s ) ∥ E → E 1 ( k − [ s + k 2 ] + 1 ) τ
(80)
+ 1 π ∑ m = s + 1 [ s + k 2 ] − 1 | Γ ( k − m + 1 2 ) ( k − m ) ! − Γ ( k − m + 3 2 ) ( k − m + 1 ) ! | ∥ u τ ( m − 1 , s ) ∥ E → E τ − 1 2
(81)
≤ 1 π 1 ( k − s ) τ + 1 π M 2 ( k − s ) τ
(82)
+M 1 2 π ∑ m = s + 1 [ s + k 2 ] − 1 τ ( k − m + 1 ) τ ( k − m ) τ ≤ M 2 ( k − s ) τ .
(83)

Estimate (71) follows from estimates (75) and (80).

Now, let us first estimate z k = ∥ D τ 1 2 u k ∥ E . Applying the triangle inequality and estimate (71), we get

z k ≤ 1 π ∑ m = 1 k Γ ( k − m + 1 2 ) ( k − m ) ! τ 1 2 [ z m + ∥ f m ∥ E ]
(84)
+ 1 π ∑ s = 1 k ∥ ∑ m = s k Γ ( k − m + 1 2 ) ( k − m ) ! A m u τ ( m , s ) ∥ E → E z s τ 3 2
(85)
+ 1 π ∑ s = 1 k ∥ ∑ m = s k Γ ( k − m + 1 2 ) ( k − m ) ! A m u τ ( m , s ) ∥ E → E ∥ f s ∥ E τ 3 2
(86)
≤ M 3 ∑ s = 1 k − 1 1 ( k − s ) τ τ [ z s + ∥ f s ∥ E ] + M 4 [ z k + ∥ f k ∥ E ] τ 1 2
(87)

for any k=1,…,N. Applying the above inequality and the difference analogue of the integral inequality, we obtain

∥ { D τ 1 2 u k } 1 N ∥ C τ ( E ) ≤M ∥ { f k } 1 N ∥ C τ ( E ) .
(88)

Using the triangle inequality and equation (3), we get

(89)

Estimate (68) follows from estimates (88) and (89). Theorem 4.1 is proved. □

With the help of A(t), we introduce the fractional spaces E α ′ = E α ′ (E,A(t)), 0<α<1, consisting of all v∈E for which the following norms are finite:

∥ v ∥ E α ′ = sup λ > 0 λ α ∥ A ( t ) ( λ + A ( t ) ) − 1 v ∥ E .
(90)

From (73) it follows that

Theorem 4.2 E α ′ (E,A(t))= E α ′ (E,A(0)) for all 0<α<1 and 0≤t≤1.

Problem (3) is not well posed in C τ (E) for arbitrary E. It turns out that a Banach space E can be restricted to a Banach space E ′ in such a manner that the restricted problem (3) in E ′ will be well posed in C( E ′ ). The role of E ′ will be played here by the fractional spaces E α = E α (A(t),E) (0<α<1).

Theorem 4.3 Suppose that assumptions (6) and (7) hold and 0<α≤ε<1. Then for the solution u τ = { u k } 1 N in C τ ( E α ′ ) of initial value problem (3), the coercive stability inequality

(91)

holds.

Proof

By Theorem 4.1,

∥ { D τ 1 2 u k } 1 N ∥ C τ ( E α ′ ) ≤M ∥ f τ ∥ C τ ( E α ′ )
(92)

for the solution of initial value problem (3). The proof of the estimate

∥ { A k u k } 1 N ∥ C τ ( E α ′ ) ≤M α − 1 ( 1 − α ) − 1 ∥ f τ ∥ C τ ( E α ′ )
(93)

for the solution of initial value problem (3) is based on estimate (92) and the following estimates [51]:

(94)
(95)

Using the triangle inequality and equation (3), we get

(96)

Estimate (91) follows from estimates (93) and (96). Theorem 4.3 is proved. □

Let us give, without proof, the following result.

Theorem 4.4 Suppose that assumptions (6) and (43) hold. Then for the solution u τ = { u k } 1 N in C τ ( E α ′ ) of initial value problem (3), the coercive stability inequality

(97)

holds.

Note that by passing to the limit for τ→0, one can recover Theorems 2.1-2.3 and 2.4.

5 Applications

Now, we consider the applications of Theorems 4.1, 4.3 and 4.4.

First, initial value problem (45) is considered. The discretization of problem (45) is carried out in two steps. In the first step, the grid space R h n (0<h≤ h 0 ) is defined as the set of all points of the Euclidean space R n whose coordinates are given by

x k = s k h, s k =0,±1,±2,…,k=1,…,n.
(98)

The difference operator A h t , x = B h t , x +σ I h is assigned to the differential operator A x = B x +σI, defined by (52). The operator

B h t , x = h − 2 m ∑ 2 m ≤ | s | ≤ S b s t , x Δ 1 − s 1 Δ 1 + s 2 ⋯ Δ n − s 2 n − 1 Δ n + s 2 n
(99)

acts on functions defined on the entire space R h n . Here s∈ R 2 n is a vector with nonnegative integer coordinates,

Δ k ± f h (x)=± ( f h ( x ± e k h ) − f h ( x ) ) ,
(100)

where e k is the unit vector of the axis x k .

An infinitely differentiable function φ(x) of the continuous argument x∈ R n that is continuous and bounded together with all its derivatives is said to be smooth. We say that the difference operator A h t , x is a λ th order (λ>0) approximation of the differential operator A t , x if the inequality

sup x ∈ R h n | A h t , x φ(x)− A t , x φ(x)|≤M(φ) h λ
(101)

holds for any smooth function φ(x). The coefficients b s t , x are chosen in such a way that the operator A h t , x approximates in a specified way the operator A t , x . It will be assumed that the operator A h t , x approximates the differential operator A t , x with any prescribed order [57, 58].

The function A t , x (ξh,h) is obtained by replacing the operator Δ k ± in the right-hand side of equality (99) with the expression ±(exp{±i ξ k h}−1), respectively, and is called the symbol of the difference operator B h t , x .

It will be assumed that for | ξ k h|≤π and fixed x, the symbol A t , x (ξh,h) of the operator B h t , x = A h t , x −σ I h satisfies the inequalities

( − 1 ) m A t , x (ξh,h)≥M | ξ | 2 m , | arg A t , x ( ξ h , h ) | ≤ϕ< ϕ 0 ≤ π 2 .
(102)

Suppose that the coefficient b s x of the operator B h t , x = A h t , x −σ I h is bounded and satisfies the inequalities

| b s t , x + e k h − b s t , x | ≤M h ϵ ,x∈ R h n ,ϵ∈(0,1].
(103)

With the help of A h t , x , we arrive at the nonlocal boundary value problem

{ d v h ( t , x ) d t + D t 1 2 v h ( t , x ) + A h t , x v h ( t , x ) = f h ( t , x ) , 0 < t < 1 , x ∈ R h n , v h ( 0 , x ) = 0 , x ∈ R h n
(104)

for an infinite system of ordinary differential equations.

In the second step, problem (104) is replaced by the difference scheme

{ u k h ( x ) − u k − 1 h ( x ) τ + 1 π ∑ m = 1 k Γ ( k − m + 1 2 ) ( k − m ) ! u m h − u m − 1 h τ 1 2 + A h k , x u k h = f k h ( x ) , f k h ( x ) = f h ( t k , x ) , t k = k τ , 1 ≤ k ≤ N − 1 , N τ = 1 , x ∈ R h n , u 0 h ( x ) = 0 , x ∈ R h n .
(105)

Based on the number of corollaries of the abstract theorems given in the above, to formulate the result, one needs to introduce the spaces C h =C( R h n ) and C h β = C β ( R h n ) of all bounded grid functions u h (x) defined on R h n , equipped with the norms

(106)
(107)

Theorem 5.1 Suppose that assumptions (102) and (103) for the operator A h k , x hold. Then, the solutions of difference scheme (105) satisfy the following stability estimates:

max 1 ≤ k ≤ N ∥ D τ 1 2 u k h ∥ C h μ ≤ M 1 ( μ ) max 1 ≤ k ≤ N ∥ f k h ∥ C h μ , 0 ≤ μ ≤ 1 , ∥ { τ − 1 ( u k h − u k − 1 h ) } 1 N ∥ C τ ( C h μ + 2 m α ) + ∥ { A k u k } 1 N ∥ C τ ( C h μ + 2 m α ) ≤ M ( α , μ ) max 1 ≤ k ≤ N ∥ f k h ∥ C h μ + 2 m α , 0 < 2 m α + μ < 1 .
(108)

The proof of Theorem 5.1 is based on the abstract Theorems 4.1, 4.3, 4.4 and the strong positivity of the operator A h x defined by (114) in C h μ and on the following two theorems on the coercivity inequality for the solution of the elliptic difference equation in C h β and on the structure of the fractional space E α ′ ( C h , A h x ).

Theorem 5.2 Suppose that assumptions (102) and (103) for the operator A h k , x hold. Then for the solutions of the elliptic difference equation

A h k , x u h (x)= ω h (x),x∈ R h n ,
(109)

the estimates [54]

∑ 2 m ≤ | s | ≤ S h − 2 m ∥ Δ 1 − s 1 Δ 1 + s 2 ⋯ Δ n − s 2 n − 1 Δ n + s 2 n u h ∥ C h β ≤M(σ,β) ∥ ω h ∥ C h β
(110)

are valid.

Theorem 5.3 Suppose that assumptions (102) and (103) for the operator A h k , x hold. Then for any 0<α< 1 2 m , the norms in the spaces E α ′ ( C h , A h x ) and C h 2 m α are equivalent uniformly in h [51].

Second, we consider mixed boundary value problem (56). The discretization of problem (56) is carried out in two steps. In the first step, let us define the grid space

[ 0 , 1 ] h ={x: x r =rh,0≤r≤K,Kh=1}.
(111)

We introduce the Banach space C h β = C β ( [ 0 , 1 ] h ) (0<β<1) of the grid functions φ h (x)= { φ r } 1 K − 1 defined on [ 0 , 1 ] h , equipped with the norm

∥ φ h ∥ C h β = ∥ φ h ∥ C h + sup 1 ≤ k < k + τ ≤ K − 1 | φ k + r − φ k | τ β ,
(112)

where C h =C( [ 0 , 1 ] h ) is the space of the grid functions φ h (x)= { φ r } 1 K − 1 defined on [ 0 , 1 ] h , equipped with the norm

∥ φ h ∥ C h = max 1 ≤ k ≤ K − 1 | φ k |.
(113)

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

A h t , x φ h (x)= { − ( a ( t , x ) φ x − ) x , r + δ φ r } 1 K − 1 ,
(114)

acting in the space of grid functions φ h (x)= { φ r } 0 K satisfying the conditions φ 0 = φ K , φ 1 − φ 0 = φ K − φ K − 1 . With the help of A h x , we arrive at the initial boundary value problem

{ d v h ( t , x ) d t + D t 1 2 v h ( t , x ) + A h t , x v h ( t , x ) = f h ( t , x ) , 0 < t < 1 , x ∈ [ 0 , 1 ] h , v h ( 0 , x ) = 0 , x ∈ [ 0 , 1 ] h
(115)

for an infinite system of ordinary fractional differential equations. In the second step, we replace problem (115) by difference scheme (3)

(116)

Theorem 5.4 Let Ï„ and h be sufficiently small numbers. Then, the solutions of difference scheme (116) satisfy the following stability estimates:

max 1 ≤ k ≤ N ∥ D τ 1 2 u k h ∥ C h μ ≤ M 1 ( μ ) max 1 ≤ k ≤ N ∥ f k h ∥ C h μ , 0 ≤ μ ≤ 1 , ∥ { τ − 1 ( u k h − u k − 1 h ) } 1 N ∥ C τ ( C h μ + 2 α ) + ∥ { A k u k } 1 N ∥ C τ ( C h μ + 2 α ) ≤ M ( α , μ ) max 1 ≤ k ≤ N ∥ f k h ∥ C h μ + 2 α , 0 < 2 α + μ < 1 .
(117)

The proof of Theorem 5.4 is based on the abstract Theorems 4.1, 4.3, 4.4 and the strong positivity of the operator A h t , x defined by (114) in C h μ and on the following theorem on the structure of the fractional space E α ′ ( C h , A h t , x ).

Theorem 5.5 For any 0<α< 1 2 , the norms in the spaces E α ′ ( C h , A h t , x ) and C h 2 α are equivalent uniformly in h and t∈[0,1] [60].

References

  1. Podlubny I Mathematics in Science and Engineering 198. In Fractional Differential Equations. Academic Press, San Diego; 1999.

    Google Scholar 

  2. Samko SG, Kilbas AA, Marichev OI: Fractional Integrals and Derivatives. Gordon & Breach, Yverdon; 1993.

    MATH  Google Scholar 

  3. Kilbas AA, Sristava HM, Trujillo JJ 204. In Theory and Applications of Fractional Differential Equations. Elsevier, Amsterdam; 2006.

    Google Scholar 

  4. El-Mesiry AEM, El-Sayed AMA, El-Saka HAA: Numerical methods for multi-term fractional (arbitrary) orders differential equations. Appl. Math. Comput. 2005, 160(3):683-699. 10.1016/j.amc.2003.11.026

    Article  MATH  MathSciNet  Google Scholar 

  5. Diethelm K: The Analysis of Fractional Differential Equations. Springer, Berlin; 2010.

    Book  MATH  Google Scholar 

  6. Pedas A, Tamme E: On the convergence of spline collocation methods for solving fractional differential equations. J. Comput. Appl. Math. 2011, 235(12):3502-3514. 10.1016/j.cam.2010.10.054

    Article  MATH  MathSciNet  Google Scholar 

  7. Diethelm K, Ford NJ: Multi-order fractional differential equations and their numerical solution. Appl. Math. Comput. 2004, 154(3):621-640. 10.1016/S0096-3003(03)00739-2

    Article  MATH  MathSciNet  Google Scholar 

  8. El-Sayed AMA, Gaafar FM: Fractional order differential equations with memory and fractional-order relaxation-oscillation model. Pure Math. Appl. 2001, 12(3):296-310.

    MathSciNet  Google Scholar 

  9. Basset AB: On the descent of a sphere in a viscous liquid. Q. J. Math. 1910, 42: 369-381.

    Google Scholar 

  10. Mainardi F: Fractional calculus: some basic problems in continuum and statistical mechanics. In Fractals and Fractional Calculus in Continuum Mechanics. Edited by: Carpinteri A, Mainardi F. Springer, New York; 1997:291-348.

    Chapter  Google Scholar 

  11. Momani S, Al-Khaled K: Numerical solutions for systems of fractional differential equations by the decomposition method. Appl. Math. Comput. 2005, 162(3):1351-1365. 10.1016/j.amc.2004.03.014

    Article  MATH  MathSciNet  Google Scholar 

  12. Bagley RL, Torvik PJ: On the appearance of the fractional derivative in the behavior of real materials. J. Appl. Mech. 1984, 51: 294-298. 10.1115/1.3167615

    Article  MATH  Google Scholar 

  13. Matignon D: Stability results for fractional differential equations with applications to control processing. 2. Computational Engineering in System Application, Lille, France, 1996.

    Google Scholar 

  14. Lavoie JL, Osler TJ, Tremblay R: Fractional derivatives and special functions. SIAM Rev. 1976, 18(2):240-268. 10.1137/1018042

    Article  MATH  MathSciNet  Google Scholar 

  15. Tarasov VE: Fractional derivative as fractional power of derivative. Int. J. Math. 2007, 18: 281-299. 10.1142/S0129167X07004102

    Article  MATH  MathSciNet  Google Scholar 

  16. Ashyralyev A, Dal F, Pinar Z: On the numerical solution of fractional hyperbolic partial differential equations. Math. Probl. Eng. 2009., 2009: Article ID 730465

    Google Scholar 

  17. Ashyralyev A: A note on fractional derivatives and fractional powers of operators. J. Math. Anal. Appl. 2009, 357(1):232-236. 10.1016/j.jmaa.2009.04.012

    Article  MATH  MathSciNet  Google Scholar 

  18. Podlubny, I, El-Sayed, AMA: On two definitions of fractional calculus. Slovak Academy of Science, Institute of Experimental Physics (1996)

    Google Scholar 

  19. Yakar A, Koksal ME: Existence results for solutions of nonlinear fractional differential equations. Abstr. Appl. Anal. 2012., 2012: Article ID 267108

    Google Scholar 

  20. De la Sen M: Positivity and stability of the solutions of Caputo fractional linear time-invariant systems of any order with internal point delays. Abstr. Appl. Anal. 2011., 2011: Article ID 161246

    Google Scholar 

  21. Yuan C:Two positive solutions for (n−1,1)-type semipositone integral boundary value problems for coupled systems of nonlinear fractional differential equations. Commun. Nonlinear Sci. Numer. Simul. 2012, 17(2):930-942. 10.1016/j.cnsns.2011.06.008

    Article  MATH  MathSciNet  Google Scholar 

  22. De la Sen M, Agarwal RP, Ibeas A, et al.: On the existence of equilibrium points, boundedness, oscillating behavior and positivity of a SVEIRS epidemic model under constant and impulsive vaccination. Adv. Differ. Equ. 2011., 2011: Article ID 748608

    Google Scholar 

  23. De la Sen M: About robust stability of Caputo linear fractional dynamic systems with time delays through fixed point theory. Fixed Point Theory Appl. 2011., 2011: Article ID 867932

    Google Scholar 

  24. Yuan C:Multiple positive solutions for semipositone (n,p)-type boundary value problems of nonlinear fractional differential equations. Anal. Appl. 2011, 9(1):97-112. 10.1142/S0219530511001753

    Article  MATH  MathSciNet  Google Scholar 

  25. Yuan C:Multiple positive solutions for (n−1,1)-type semipositone conjugate boundary value problems for coupled systems of nonlinear fractional differential equations. Electron. J. Qual. Theory Differ. Equ. 2011., 2011: Article ID 13

    Google Scholar 

  26. Agarwal RP, Belmekki M, Benchohra M: A survey on semilinear differential equations and inclusions involving Riemann-Liouville fractional derivative. Adv. Differ. Equ. 2009., 2009: Article ID 981728

    Google Scholar 

  27. Agarwal RP, de Andrade B, Cuevas C: On type of periodicity and ergodicity to a class of fractional order differential equations. Adv. Differ. Equ. 2010., 2010: Article ID 179750

    Google Scholar 

  28. Agarwal RP, de Andrade B, Cuevas C: Weighted pseudo-almost periodic solutions of a class of semilinear fractional differential equations. Nonlinear Anal., Real World Appl. 2010, 11: 3532-3554. 10.1016/j.nonrwa.2010.01.002

    Article  MATH  MathSciNet  Google Scholar 

  29. Berdyshev AS, Cabada A, Karimov ET: On a non-local boundary problem for a parabolic-hyperbolic equation involving a Riemann-Liouville fractional differential operator. Nonlinear Anal. 2011, 75(6):3268-3273.

    Article  MathSciNet  Google Scholar 

  30. Ashyralyev A, Amanov D: Initial-boundary value problem for fractional partial differential equations of higher order. Abstr. Appl. Anal. 2012., 2012: Article ID 973102

    Google Scholar 

  31. 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., 2012: Article ID 594802

    Google Scholar 

  32. Ashyralyev A, Hicdurmaz B: A note on the fractional Schrödinger differential equations. Kybernetes 2011, 40(5-6):736-750.

    MathSciNet  Google Scholar 

  33. Ashyralyev A, Dal F, Pinar Z: A note on the fractional hyperbolic differential and difference equations. Appl. Math. Comput. 2011, 217(9):4654-4664. 10.1016/j.amc.2010.11.017

    Article  MATH  MathSciNet  Google Scholar 

  34. Kirane M, Laskri Y: Nonexistence of global solutions to a hyperbolic equation with a space-time fractional damping. Appl. Math. Comput. 2005, 167(2):1304-1310. 10.1016/j.amc.2004.08.038

    Article  MATH  MathSciNet  Google Scholar 

  35. Kirane M, Laskri Y, Tatar N-e: Critical exponents of Fujita type for certain evolution equations and systems with spatio-temporal fractional derivatives. J. Math. Anal. Appl. 2005, 312(2):488-501. 10.1016/j.jmaa.2005.03.054

    Article  MATH  MathSciNet  Google Scholar 

  36. Kirane M, Malik SA: The profile of blowing-up solutions to a nonlinear system of fractional differential equations. Nonlinear Anal., Theory Methods Appl. 2010, 73(12):3723-3736. 10.1016/j.na.2010.06.088

    Article  MATH  MathSciNet  Google Scholar 

  37. Araya D, Lizama C: Almost automorphic mild solutions to fractional differential equations. Nonlinear Anal. 2008, 69(11):3692-3705. 10.1016/j.na.2007.10.004

    Article  MATH  MathSciNet  Google Scholar 

  38. N’Guerekata GM: A Cauchy problem for some fractional abstract differential equation with nonlocal conditions. Nonlinear Anal. 2009, 70: 1873-1876. 10.1016/j.na.2008.02.087

    Article  MATH  MathSciNet  Google Scholar 

  39. Mophou GM, N’Guerekata GM: Mild solutions for semilinear fractional differential equations. Electron. J. Differ. Equ. 2009., 2009: Article ID 21

    Google Scholar 

  40. Mophou GM, N’Guerekata GM: Existence of mild solution for some fractional differential equations with nonlocal conditions. Semigroup Forum 2009, 79(2):315-322. 10.1007/s00233-008-9117-x

    Article  MATH  MathSciNet  Google Scholar 

  41. Lakshmikantham V: Theory of fractional differential equations. Nonlinear Anal. 2008, 69(10):3337-3343. 10.1016/j.na.2007.09.025

    Article  MATH  MathSciNet  Google Scholar 

  42. Lakshmikantham V, Devi JV: Theory of fractional differential equations in Banach spaces. Eur. J. Pure Appl. Math. 2008, 1: 38-45.

    MATH  MathSciNet  Google Scholar 

  43. Lakshmikantham V, Vatsala A: Theory of fractional differential inequalities and applications. Commun. Appl. Anal. 2007, 11(3-4):395-402.

    MATH  MathSciNet  Google Scholar 

  44. Clement P, Guerre-Delabrire S: On the regularity of abstract Cauchy problems and boundary value problems. Atti Accad. Naz. Lincei, Rend. Lincei, Mat. Appl. 1998, 9(4):245-266.

    MATH  Google Scholar 

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

    MathSciNet  Google Scholar 

  46. Shakhmurov VB: Coercive boundary value problems for regular degenerate differential-operator equations. J. Math. Anal. Appl. 2004, 292(2):605-620. 10.1016/j.jmaa.2003.12.032

    Article  MATH  MathSciNet  Google Scholar 

  47. Lunardi A: Analytic semigroups and optimal regularity in parabolic problems. In Operator Theory Advances and Applications. Birkhäuser, Basel; 1995.

    Google Scholar 

  48. Ashyralyev A, Hanalyev A, Sobolevskii PE: Coercive solvability of nonlocal boundary value problem for parabolic equations. Abstr. Appl. Anal. 2001, 6(1):53-61. 10.1155/S1085337501000495

    Article  MATH  MathSciNet  Google Scholar 

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

    MathSciNet  Google Scholar 

  50. Sobolevskii PE: Some properties of the solutions of differential equations in fractional spaces. 14. Trudy Nauchn.-Issled. Inst. Mat. Voronezh. Gos. Univ. No. 1975, 68-74. (Russian)

    Google Scholar 

  51. Ashyralyev A, Sobolevskii PE: New difference schemes for partial differential equations. 148. In Operator Theory: Advances and Applications. Birkhäuser, Basel; 2004.

    Google Scholar 

  52. Ashyralyev A: Well-posedness of the Basset problem in spaces of smooth functions. Appl. Math. Lett. 2011, 24: 1176-1180. 10.1016/j.aml.2011.02.002

    Article  MATH  MathSciNet  Google Scholar 

  53. Ashyralyev A, Cakir Z: On the numerical solution of fractional parabolic partial differential equations. AIP Conf. Proc. 2011, 1389: 617-620.

    Article  Google Scholar 

  54. Ashyralyev A: Well-posedness of parabolic differential and difference equations with the fractional differential operator. Malays. J. Math. Sci. 2012, 6(S):73-89.

    MathSciNet  Google Scholar 

  55. Cakir Z: Stability of difference schemes for fractional parabolic PDE with the Dirichlet-Neumann conditions. Abstr. Appl. Anal. 2012., 2012: Article ID 463746

    Google Scholar 

  56. Ashyralyev A, Cakir Z: On the numerical solution of fractional parabolic partial differential equations with the Dirichlet condition. Discrete Dyn. Nat. Soc. 2012., 2012: Article ID 696179

    Google Scholar 

  57. Smirnitskii YA, Sobolevskii PE: Positivity of multidimensional difference operators in the C -norm. Usp. Mat. Nauk 1981, 36(4):202-203. (Russian)

    Google Scholar 

  58. Smirnitskii, YA: Fractional powers of elliptic difference operators. PhD thesis, Voronezh State University, Voronezh (1983) (Russian)

    Google Scholar 

  59. Ashyralyev A, Sobolevskii PE: Well-posedness of parabolic difference equations. In Operator Theory Advances and Applications. Birkhäuser, Basel; 1994.

    Google Scholar 

  60. Ashyralyev A: Fractional spaces generated by the positive differential and difference operator in a Banach space. In Proceedings of the Conference ‘Mathematical Methods and Engineering’. Edited by: Taş K, Tenreiro Machado JA, Baleanu D. Springer, Berlin; 2007:13-22.

    Google Scholar 

Download references

Acknowledgements

The author would like to thank Prof. P. E. Sobolevskii for his helpful suggestions to the improvement of this paper.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Allaberen Ashyralyev.

Additional information

Competing interests

The author declares that they have no competing interests.

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. Well-posedness of fractional parabolic equations. Bound Value Probl 2013, 31 (2013). https://doi.org/10.1186/1687-2770-2013-31

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-31

Keywords