Even Higher Order Fractional Initial Boundary Value Problem with Nonlocal Constraints of Purely Integral Type

Abstract: In this paper, the a priori estimate method, the so-called energy inequalities method based on some functional analysis tools is developed for a Caputo time fractional 2mth order diffusion wave equation with purely nonlocal conditions of integral type. Existence and uniqueness of the solution are proved. The proofs of the results are based on some a priori estimates and on some density arguments.


Introduction
Classical initial boundary value problems for partial differential equations with integer and noninteger order have been widely studied during the last three decades by using different methods.One of the most important methods used and applied to linear and nonlinear partial differential equations with integer order supplemented with classical conditions is the functional analysis method.However, for equations with Caputo time fractional order and nonlocal conditions, there are only a few results obtained by using the mentioned method.The Caputo fractional derivative is a nonlocal operator since it is an integral which is a nonlocal operator.Caputo time fractional derivative can be used to model systems with memory, since it requires all the past history.Time fractional order partial differential equations play a great role in reducing the errors coming from the neglected parameters while modeling real life phenomena.
One of the most important classes of the above equations are the fractional diffusion-wave equations that have been studied and used in different branches of Science.Problem (1) constitutes a large class of time fractional diffusion wave equations of even order such as second and fourth order time fractional wave equations that have numerous applications in physics and engineering as mentioned below.In our problem, local conditions at 0 and 1 are replaced by other conditions on the moments of order 1, 2, ..., 2m − 1 which are non-local integral conditions.Although mathematical models in two and three-dimensions are of big significance for applications, most of the recent research articles are devoted to the fractional order diffusion wave equations in one-dimensional settings.These equations model, for example, propagation of mechanical waves in viscoelastic media [1][2][3][4], a non-Markovian diffusion process with memory [5], and a model governing the propagation of mechanical diffusive waves in viscoelastic media that exhibit a power-law creep [2][3][4].
In the literature, many researchers used the functional analysis method to investigate the well posedness of initial boundary value problems for partial differential equations with time and space integer order having nonlocal conditions-we cite, for example, the references [14][15][16][17].For the fractional diffusion wave equations case with higher order derivatives and classical boundary conditions, there are only few papers dealing with the existence and uniqueness of solution such as [18][19][20].In this paper, an initial boundary value problem with purely nonlocal constraints of integral type for a Caputo time fractional 2mth order diffusion wave equation is studied by applying the functional analysis method, the so-called energy inequality method based mainly on some a priori estimates and on the density of the range of the operator generated by the studied problem.This work can be considered as a contribution to the development of the functional analysis method used to prove the well posedness of problems with fractional order.The obtained results show the efficiency of this method to study the existence and uniqueness of solution for the time fractional order differential equations with nonlocal conditions.This paper is organized as follows: in Section 2, we set our fractional initial boundary value problem.In Section 3, we give some preliminaries concerning the used function spaces, some useful tools and write down the given problem in its operator form.In Section 4, we establish an a priori estimate for the solution and deduce some consequences about the uniqueness of the solution and its dependence on the free term and the given data.Section 5 provides proofs of the main result concerning the solvability of the posed problem.We end our problem with conclusions.

Problem Setting
In the domain Q = (0, 1) × (0, T) where 0 ≤ T < ∞, we consider the time fractional initial boundary problem of higher order with purely integral conditions where θ(t), f (x, t), g(x) and h(x) are given functions that satisfy certain conditions which will be specified later on, and the operator ∂ α+1 t denotes the Caputo left fractional derivative of order 1 + α with 0 < α < 1 defined by (see [21]) where Γ(1 − α) is the Gamma function.
The Riemann-Liouville integral of order 0 < α < 1 is defined by (see [21]) Different properties of the Caputo fractional derivative and Riemann fractional-Liouville integral can be found in [21][22][23] and the references therein.

Preliminaries
In this section, we introduce some important lemmas and inequalities needed throughout the sequel, and write the posed problem in its operator form.
Lemma 1 (Poincare type inequality).For m ∈ N, we have where Lemma 2 ([24]).For any absolutely continuous function J(s) on the interval [0, T], the following inequality hold Lemma 3 ([25]).Let ϕ(t) be nonnegative and absolutely continuous on [0, T], and, for almost all t ∈ [0, T], satisfies the inequality where the functions C(t) and B(t) are summable and nonnegative on [0, T].Then, Lemma 3 can be generalized as Lemma 4 ([24]).Let a nonnegative absolutely continuous function Z (t) satisfy the inequality for almost all t ∈ [0, T], where c 1 is a positive constant and c 2 (t) is an integrable nonnegative function on where are Mittag-Leffler functions.
A special case of ( 11) is the Cauchy inequality with ε: The solution of the problem (1) can be regarded as the solution of operator equation where M = (L, l 1 , l 2 ), and M : B −→ Y is an unbounded operator with domain of definition such that v satisfies the initial conditions and where B is a Banach space of functions v endowed with the finite norm and Y is Hilbert space constituting of the elements W = ( f , g, h) equipped with the norm Here, L denotes the time fractional differential operator

A Priori Estimate for the Solution and Uniqueness
To prove the uniqueness of solution of problem (1), we establish an energy inequality for the solution from which we deduce the uniqueness of solution of the posed problem.

Theorem 1. Assume that the function θ(t) satisfies the conditions
where c 1 , c 2 , c 3 , and c 4 are positive constants.Then, for any v ∈ D(M), there exists a positive constant K such that the following a priori estimate is satisfied: where K = K(η, δ, ρ) is given by with η, δ and ρ are respectively given by (30), (34) and (37).
Proof.For v ∈ D(L), we consider the scalar product in L 2 (0, 1) of the differential equation in problem (1) and the integrodifferential operator Nv = 2(−1) m I 2m x v t , we have We separately consider the inner products on the left-hand side of Equation ( 20) and we integrate by parts and taking into account boundary and initial conditions in Problem (1), we obtain Substitution of ( 21) and ( 22) into (20) yields By Lemmas 1 and 2 and inequality (12), identity (23) reduces to Replacing t by τ , integrating with respect to τ from zero to t and using given conditions, we obtain The second term on the left-hand side can be evaluated as Hence, inequality (25) becomes evoking conditions (17) and using (28), we infer from (26) that If, in Lemma 3, we set , and ϕ(0) = 0, (31) then it yields Consequently, (30) transforms to where Now, by dropping the second term on the left-hand side of (33) then setting Z (t) = in Lemma 4, we obtain Combination of ( 33) and (36) leads to where It is obvious that Then, it follows from (37) and (38) that where Observe that the right-hand side of (39) is independent of the variable t, so we are allowed to take the least upper bound of the left-hand side with respect to t over [0, T], and the a priori estimate (18) then follows and from which we deduce the uniqueness and continuous dependence of the solution on the input data of problem (1).

Existence of Solution
In this section, we prove the main result concerning the existence of solution of problem (1) The a priori estimate (18) shows that the unbounded operator M has an inverse M −1 : R(M) →B.Since R(M) is a subset of Y, we therefore can construct its closure M in a manner that the estimate (18) holds for this extension and R(M) coincides with the whole space B. Hence, the following: Corollary 1.The operator M : B→ Y admits a closure (the proof is similar to that in [14]).
Estimate (18) can be then extended to for all v ∈ D(M).
It follows from (40) that the strong solution of problem ( 1) is unique, that is, Mv = Y.We also deduce from estimate (40) the following: We are now ready to give the result of existence of the solution of problem (1).
Theorem 2. Suppose that conditions of Theorem 4.1 are satisfied.Then, for all W = ( f , g, h) ∈ Y, there exists a unique strong solution v = M −1 W = M −1 W of problem (1).
Proof.Estimate (40) asserts that, if a strong solution of (1) exists, it is unique and depends continuously on the data.Corollary 2 says that, in order to prove that problem (1) admits a strong solution for any W = ( f , g, h) ∈ Y, it suffices to show that the closure of the range of the operator M is dense in Y.
To establish the existence of the strong solution of problem (1), we use a density argument.That is, we show that the range R(M) of the operator M is dense in the space Y for every element v in the Banach space B. For this, we consider the following special case of density.
Theorem 3. Suppose that conditions of Theorem 1 are satisfied.If for all functions v ∈ D(M) such that then ψ vanishes a.e in Q.

Conclusions
The existence and uniqueness of a generalized solution for a higher order fractional diffusion wave equation in Caputo sense subject to initial and weighted integral boundary conditions are established.It is found that the method of energy inequalities is successfully applied to obtaining a priori bounds for the solution of fractional initial boundary value problems of higher order with nonlocal constraints as in the classical case.The obtained results will contribute to the development of the functional analysis method and enrich the existing non-extensive literature on the non local fractional mixed problems in the Caputo sense.