arxiv: v1 [math.ap] 3 Jul 2015

Transcription

1 Strong Maximum Principle for Fractional Diffusion Equations and an Application to an Inverse Source Problem Yikan LIU William RUNDELL Masahiro YAMAMOTO arxiv:57.845v [math.ap] 3 Jul 25 Abstract The strong maximum principle is a remarkable characterization of parabolic equations, which is expected to be partly inherited by fractional diffusion equations. Based on the corresponding weak maximum principle, in this paper we establish a strong maximum principle for time-fractional diffusion equations with Caputo derivatives, which is slightly weaker than that for the parabolic case. As a direct application, we give a uniqueness result for a related inverse source problem on the determination of the temporal component of the inhomogeneous term. Keywords Fractional diffusion equation, Caputo derivative, Strong maximum principle, Mittag-Leffler function, Inverse source problem, Fractional Duhamel s principle AMS Subject Classifications 35R, 26A33, 35B5, 35R3 Introduction and main results Let Ω R d (d =,2,3) be an open bounded domain with a smooth boundary (for example, of C class), T > and < α <. Consider the following initial-boundary value problem for a time-fractional diffusion equation t α u(x,t)+au(x,t) = F(x,t) (x Ω, < t T), u(x,) = a(x) (x Ω), (.) u(x,t) = (x Ω, < t T), where t α denotes the Caputo derivative defined by α t f(t) := Γ( α) f (s) (t s) α ds, and Γ( ) denotes the Gamma function. Here A is an elliptic operator defined for f D(A) := H 2 (Ω) H (Ω) as Af(x) = d j (a ij (x) i f(x))+c(x)f(x) (x Ω), (.2) i,j= Manuscript last updated: March 7, 28. Graduate School of Mathematical Sciences, The University of Tokyo, 3-8- Komaba, Meguro-ku, Tokyo , Japan. ykliu@ms.u-tokyo.ac.jp, myama@ms.u-tokyo.ac.jp Department of Mathematics, Texas A&M University, College Station, TX , USA. rundell@math.tamu.edu Y. Liu and M. Yamamoto have been partly supported by the A3 Foresight Program Modeling and Computation of Applied Inverse Problems by Japan Society of the Promotion of Science.

2 2 Y. Liu, W. Rundell and M. Yamamoto where a ij = a ji ( i,j d) and c in Ω. Moreover, it is assumed that a ij C (Ω), c C(Ω) and there exists a constant δ > such that δ d d ξi 2 a ij (x)ξ i ξ j ( x Ω, (ξ,...,ξ d ) R d ). i= i,j= The assumptions on the initial data a and the source term F will be specified later. Fractional diffusion equations, especially the governing equation in (.) with a Caputo derivative in time, have been widely used as model equations for describing the anomalous diffusion phenomena in highly heterogeneous aquifer and complex viscoelastic material (see [, 8,, 2, 22]). Due to its practical applications, (.) has drawn extensive attentions of mathematical researchers during the recent years. In Luchko [7], the generalized solution to (.) with F = was represented by means of the Mittag-Leffler function, and the unique existence of the solution was proved. Sakamoto and Yamamoto [27] investigated the wellposedness and the asymptotic behavior of the solution to(.). Very recently, Gorenflo et al.[] re-defined the Caputo derivative in the fractional Sobolev spaces and investigated(.) from the viewpoint of the operator theory. Regarding numerical treatments, we refer e.g. to [5, 2] for the finite difference method and [2,3] for the finite element method. Meanwhile, (.) has also gained population among the inverse problem school; recent literatures include [4,9,24]. Here we do not intend to enumerate a complete list of related works. It reveals in the existing works that fractional diffusion equations show certain similarities to classical parabolic equations (i.e., α = in (.)), whereas also diverge considerably from their integer prototypes in the senses of the limited smoothing property in space and slow decay in time. Other than the above mentioned aspects, the maximum principle is also one of the remarkable characterizations of parabolic equations, which is not only significant by itself but also applicable in many related problems. However, researches especially on the strong maximum principle for time-fractional diffusion equations with Caputo derivatives are inadequate due to the technical difficulties in treating the fractional derivatives. Luchko [6] established a weak maximum principle for (.) by a key estimate of the Caputo derivative at an extreme point, by which the uniqueness of a classical solution was also proved. On the other hand, both weak and strong maximum principles were recently obtained for time-fractional diffusion equations with Riemann-Liouville derivatives (see Al-Refai and Luchko [3]). In this paper, we are interested in improving the maximum principle for fractional diffusion equations with Caputo derivatives. Based on the weak maximum principle obtained in [6] (see Lemma 2.3), first we establish a strong maximum principle for the initial-boundary value problem (.), which is slightly weaker than that for the parabolic case. Theorem. Let a L 2 (Ω) satisfy a and a, F =, and u be the solution to (.) with d 3. Then for any x Ω, the set E x := {t > ; u(x,t) } is at most a finite set. Remark.2 (a) By the Sobolev embedding and Lemma 2.2(a) in Section 2, we see u C(Ω (, )) in Theorem. and thus the set E x is well-defined. According to the weak maximum principle, it reveals that E x is actually the set of zero points of u(x,t) as a function of t, and Theorem. asserts the strict positivity of u(x,t) for all x Ω and almost all t > except for the finite set E x.

3 Strong Maximum Principle for Fractional Diffusion Equations 3 (b) Note that we have stated Theorem. for spatial dimensions d 3. This can be generalized to arbitrary d provided that our initial data a has sufficient regularity to allow a pointwise definition. For d > 3 this will mean restricting a in a subset of L 2 (Ω). If this is done, then the set E x in Theorem. will again be well-defined and, in addition, in (3.8) we can observe that A 3 can be replaced by any higher power k necessary since the crucial requirement of C (ω) D(Ak ) is satisfied. However, to keep the exposition simpler, we shall make the restriction d 3 throughout the remainder of the paper. (c) It is an immediate consequence of Theorem. that u > a.e. in Ω (, ). To see this, we investigatethe set D := {(x,t) Ω (, ); u(x,t) } and notice D ({x} (, )) = E x. Since the characteristic function χ Ex = a.e. in (, ) by Theorem., it follows from Fubini s theorem that D = χ D (x,t)dxdt = χ Ex (t)dtdx =, Ω (, ) where denotes the Lebesgue measure. (d) If the inhomogeneous term F in Theorem. is allowed to be non-negative in Ω (,T), then it follows immediately from the weak maximum principle that u > a.e. in Ω (,T). So far, we do not know if E x = ( x Ω) although it can be conjectured. Nevertheless, we can prove the following result. Corollary.3 Let a L 2 (Ω) satisfy a > a.e. in Ω, F =, and u be the solution to (.). Then u > in Ω (, ). Theorem. is a weaker result than our expected strong maximum principle, but is sufficient for some application. Next we study an inverse source problem for (.) under the assumption that the inhomogeneous term F takes the form of separation of variables. Problem.4 Let x Ω and T > be arbitrarily given, and u be the solution to (.) with a = and F(x,t) = ρ(t)g(x). Provided that g is known, determine ρ(t) ( t T) by the single point observation data u(x,t) ( t T). The above problem is concerned with the determination of the temporal component ρ in the inhomogeneous term F(x,t) = ρ(t)g(x) in (.). The spatial component g simulates e.g. a source of contaminants which may be dangerous. Although g is usually limited to a small region given by suppg( Ω), its influence may expand wider because ρ(t) is large. We are requested to determine the time-dependent magnitude by the pointwise data u(x,t) ( t T), where x is understood as a monitoring point. For the case of x suppg, we know the both-sided stability estimate as well as the uniqueness for Problem.4 (see Sakamoto and Yamamoto [27]). For the case of x / suppg, there were no published results even on the uniqueness. From the practical viewpoints mentioned above, it is very desirable that x should be spatially far from the location of the source, that is, the case x / suppg should be discussed for the inverse problem. As a direct application of Theorem., we can give an affirmative answer for the uniqueness regarding Problem.4. Theorem.5 Under the same settings in Problem.4, we further assume that ρ C [,T], g D(A ε ) with some ε > (see Section 2 for the definition of D(A ε )), g and g. Then Ω

4 4 Y. Liu, W. Rundell and M. Yamamoto u(x,t) = ( t T) implies ρ(t) = ( t T). In the case of x / suppg, the condition g and g in Ω is essential for the uniqueness. In fact, as a simple counterexample, we consider Ω = (,) and g(x) = sin2πx in Problem.4, where the condition g is not satisfied. In this case, it follows from [27] that u(x,t) = s α E α,α ( 4π 2 s α )ρ(t s)dssin2πx, where E α,α ( ) denotes the Mittag-Leffler function (see (2.)). It is readily seen that u(/2,t) = for t > and any ρ C [,T]. In other words, the data at x = /2 does not imply the uniqueness for Problem.4. For the same kind of inverse source problems for parabolic equations, we refer to Cannon and Esteva [4], Saitoh, Tuan and Yamamoto [25, 26]. The rest of this paper is organized as follows. Section 2 introduces the notations and collects the existing results concerning problem (.). Sections 3 is devoted to the proofs of Theorem. and Corollary.3, and the proof of Theorem.5 is given in Section 4. 2 Preliminaries Tostartwith, wefixsomegeneralsettingsandnotations. LetL 2 (Ω)beausualL 2 -spacewith the inner product (, ) and H (Ω), H2 (Ω) denote the Sobolev spaces (see, e.g., Adams [2]). Let {(λ n,ϕ n )} be the eigensystem of the symmetric uniformly elliptic operator A in (.) such that < λ < λ 2 (the multiplicity is also counted), λ n as n and {ϕ n } H 2 (Ω) H(Ω) forms an orthonormal basis of L 2 (Ω). Then we can define the fractional power A γ for γ as { } D(A γ ) = f L 2 (Ω); λ γ n(f,ϕ n ) 2 <, A γ f := λ γ n(f,ϕ n )ϕ n, and D(A γ ) is a Hilbert space with the norm ( ) /2 f D(Aγ ) = λ γ n (f,ϕ n) 2. For p and a Banach space X, we say that f L p (,T;X) provided ( ) /p T f(,t) p X f L p (,T;X) := dt if p < <. ess sup f(,t) X if p = <t<t Similarly, for t < T, we set f C([t,T];X) := max f(,t) X. t t T In addition, we define C((,T];X) := C([t,T];X), C([, );X) := C([,T];X). <t <T T>

5 Strong Maximum Principle for Fractional Diffusion Equations 5 To represent the explicit solution of (.), we first recall the Mittag-Leffler function (see, e.g., Podlubny [23] and Gorenflo et al. [9]) E α,β (z) := k= which possesses the following properties. z k Γ(αk +β) (z C, α >, β R), (2.) Lemma 2. (a) Let < α < 2 and η >. Then E α, ( η) > and the following expansion holds: ( ) E α, ( η) = Γ( α)η +O η 2 as η. (b) Let < α < 2 and β R be arbitrary. Then there exists a constant C = C(α,β) > such that E α,β ( η) C (η ). +η (c) For any l =,,2,..., there holds E α,+lα (z) = (d) For λ > and α >, we have Γ(+lα) +ze α,+(l+)α(z) (α >, z C). d dt E α,( λt α ) = λt α E α,α ( λt α ) (t > ). We mention that Lemma 2.(a) (b) are well-known results from[23,.2], and (c) (d) follow immediately from direct calculations by definition (2.). Regarding some important existing results of the solution to (.), we state the following two lemmata for later use. Lemma 2.2 Fix T > arbitrarily. Concerning the solution u to (.), we have: (a) Let a L 2 (Ω) and F =. Then there exists a unique solution u C([,T];L 2 (Ω)) C((,T];H 2 (Ω) H (Ω)), which can be represented as u(,t) = E α, ( λ n t α )(a,ϕ n )ϕ n (2.2) in C([,T];L 2 (Ω)) C((,T];H 2 (Ω) H (Ω)), where {(λ n,ϕ n )} Moreover, there exists a constant C = C(Ω,T,α,A) > such that is the eigensystem of A. u(,t) L2 (Ω) C a L2 (Ω), u(,t) H2 (Ω) C a L2 (Ω)t α ( < t T). (2.3) In addition, u : (,T] H 2 (Ω) H (Ω) can be analytically extended to a sector {z C; z, argz < π/2}. (b) Let a = and F L (,T;L 2 (Ω)). Then there exists a unique solution u L 2 (,T; H 2 (Ω) H (Ω)) such that lim t u(,t) L2 (Ω) =. We note that Lemma 2.2 almost coincides with [27, Theorem 2.], but the regularity of the analyticity result stated in Lemma 2.2(a) is stronger. Indeed, one can improve the regularity up to H 2 (Ω) by the same reasoning, but here we omit the details.

6 6 Y. Liu, W. Rundell and M. Yamamoto Lemma 2.3 (Weak maximum principle) Let a L 2 (Ω) and F L (,T;L 2 (Ω)) be nonnegative, and u be the solution to (.). Then there holds u a.e. in Ω (,T). Since we impose a homogeneous Dirichlet boundary condition, Lemma 2.3 is a special case of [6, Theorem 3], but our choices of the initial data, the inhomogeneous term and the elliptic operator are more general. Fortunately, the same argument still works in our settings, which indicates Lemma 2.3 immediately. Again we omit the details here. 3 Proof of Theorem. and Corollary.3 Now we proceed to the proof of the strong maximum principle. Throughout this section, we concentrate on the homogeneous problem, that is, t α v +Av = in Ω (, ), v a v = a in Ω {}, (3.) v = on Ω (, ), where we emphasize the dependency of the solution upon the initial data a by denoting the solution as v a. To begin with, we investigate the Green function of problem (3.). Using the Mittag-Leffler function and the eigensystem {(λ n,ϕ n )}, for N N we set G N (x,y,t) := N E α, ( λ n t α )ϕ n (x)ϕ n (y) (x,y Ω, t > ). According to Lemma 2.2(a), there holds v a (x,t) = lim N Ω G N (x,y,t)a(y)dy in C([, );L 2 (Ω)) C((, );H 2 (Ω) H (Ω)) for any a L2 (Ω) and any t >. Therefore, for any fixed x Ω and t >, we see that v a (x,t) is defined pointwisely and thus G N (x,,t) is weakly convergent to G(x,y,t) := E α, ( λ n t α )ϕ n (x)ϕ n (y) (3.2) as a series with respect to y. In particular, we obtain G(x,,t) L 2 (Ω) for all x Ω and all t >. Moreover, the solution to (3.) can be represented as v a (x,t) = G(x,y,t)a(y)dy (x Ω, t > ). (3.3) Ω Next we show that for arbitrarily fixed x Ω and t >, G(x,,t) a.e. in Ω. Actually, assume that on the contrary there exist x Ω and t > such that the Lebesgue measure of the subdomain ω := {G(x,,t ) < } Ω is positive. Then it is readily seen that v χω (x,t ) = G(x,y,t )χ ω (y)dy <, (3.4) Ω where χ ω is the characteristic function of ω satisfying χ ω L 2 (Ω) and χ ω. On the other hand, Lemma 2.3 and (3.3) imply that v χω (x,t ), which contradicts with (3.4). In summary, we have proved the following lemma.

7 Strong Maximum Principle for Fractional Diffusion Equations 7 Lemma 3. Let G(x,y,t) be the Green function defined in (3.2). Then for arbitrarily fixed x Ω and t >, we have G(x,,t) L 2 (Ω) and G(x,,t) a.e. in Ω. Now we are well prepared to prove the strong maximum principle. Proof of Theorem.. We deal with the homogeneous problem (3.) with the initial data a L 2 (Ω) such that a and a. By Lemma 2.2(a) and the Sobolev embedding H 2 (Ω) C(Ω) for d 3, we have v a C(Ω (, )). Meanwhile, according to the weak maximum principle stated in Lemma 2.3, there holds v a in Ω (, ), indicating E x := {t > ; v a (x,t) } = {t > ; v a (x,t) = }, that is, E x coincides with the zero point set of v a (x,t) as a function of t >. Assume contrarilythat thereexists x Ω suchthat the sete x is notafinite set. Then E x contains at least an accumulation point t [, ]. We treat the cases of t =, t (, ) and t = separately. Case. If t =, then by definition there exists {t i } i= E x such that t i (i ) and u(x,t i ) =. Recall the explicit representation v a (x,t) = E α, ( λ n t α )(a,ϕ n )ϕ n (x ). Then the asymptotic behavior described in Lemma 2.(a) implies v a (x,t) = Γ( α)t α (a,ϕ n ) λ n ( ) (a,ϕ n ) ϕ n (x )+O t 2α λ 2 ϕ n (x ) as t. n Substituting t = t i with sufficiently large i into the above expansion, multiplying both sides by t α i and passing i, we obtain b(x ) =, where b := (a,ϕ n ) λ n ϕ n. Simple calculations reveal that b satisfies the boundary value problem { Ab = a in Ω, b = on Ω. (3.5) Since the coefficient c in the elliptic operator A is non-negative, the weak maximum principle for (3.5) (see Gilbarg and Trudinger [7, Chapter 3]) indicates b in Ω. Moreover, as b attains its minimum at x Ω, the strong maximum principle for (3.5) implies b const. = and thus a = Ab =, which contradicts with the assumption a. Therefore, cannot be an accumulation point of E x. Case 2. Now suppose that the set of zeros E x admits an accumulation point t (, ). By the analyticity of v a : (, ) H 2 (Ω) H (Ω) C(Ω), we see that v a (x,t) is analytic with respect to t >. Therefore, v a (x,t) should vanish identically if its zero points accumulate

8 8 Y. Liu, W. Rundell and M. Yamamoto at some finite and non-zero point t. Then this case reduces to Case and eventually result in a contradiction. Case 3. Since v a (x,t) is not analytic at t =, we shall treat the case of t = separately. Henceforth C > denotes generic constants independent of n N and t, which may change line by line. By definition, there exists {t i } i= E x such that t i (i ) and, in view of the representation (3.3), v a (x,t i ) = G(x,y,t i )a(y)dy = (i =,2,...). Ω SinceG(x,,t i ) bylemma3.anda, wededuceg(x,y,t i )a(y) = foralli =,2,... and almost all y Ω. Since a, it follows that G(x,,t i ) should vanish in the subdomain ω := {a > } whose Lebesgue measure is positive. By the representation (3.2), this indicates E α, ( λ n t α i )ϕ n (x )ϕ n = a.e. in ω (i =,2,...). (3.6) Now we choose ψ C v ψ (x,t) = Ω (ω) arbitrarily as the initial data of (3.) and investigate G(x,y,t)ψ(y)dy = E α, ( λ n t α )(ψ,ϕ n )ϕ n (x ). (3.7) For later convenience, we abbreviate ψ n := (ψ,ϕ n )ϕ n (x ). We shall show that the series in (3.7) is convergent in C[, ). In fact, by the Sobolev embedding H 2 (Ω) C(Ω) for d 3, first we estimate ϕ n (x ) C ϕ n H 2 (Ω) C Aϕ n L 2 (Ω) Cλ n. Next, since ψ C (ω) D(A3 ), we have (ψ,ϕ n ) = (A3 ψ,ϕ n ) λ 3 n A3 ψ L 2 (ω) ϕ n L 2 (Ω) λ 3 n C ψ C 6 (ω) λ 3. (3.8) n Ontheotherhand, weknowλ n n 2/d asn (see, e.g., CourantandHilbert[5]). Therefore, the combination of the above estimates yields E α, ( λ n t α )ψ n C ψ C6 (ω)λ 2 n Cn 4/d as n, where the boundedness of E α, ( λ n t α ) is guaranteed by Lemma 2.(b). Since the restriction d 3 gives 4/d >, we obtain E α, ( λ n t α )ψ n < ( t, ψ C (ω)), which indicates that v ψ (x,t) is well-defined in the sense of C[, ). Meanwhile, by the same reasoning as that for (3.8), for any l =,,2,... we estimate (ψ,ϕ n ) = (Al+3 ψ,ϕ n ) λ l+3 n C ψ C 2(l+3) (ω) λ l+3, n

9 Strong Maximum Principle for Fractional Diffusion Equations 9 implying λ l n ψ n C ψ C 2(l+3) (ω) λ 2 n < ( l =,,..., ψ C (ω)). (3.9) Moreover, since E α,β ( η) is uniformly bounded for all η and all β > by Lemma 2.(b), we further have λ l ne α,β ( λ n t α )ψ n < ( l =,,..., β >, t, ψ C (ω)). (3.) Utilizing Lemma 2.(c) with l =, we treat v ψ (x,t) as v ψ (x,t) = ψ n t α λ n E α,+α ( λ n t α )ψ n, where the boundedness of the involved summations were verified in (3.9) (3.). Taking t = t i and passing i, we obtain ψ n =, implying v ψ (x,t) = ( λ n t α )E α,+α ( λ n t α )ψ n. For t >, we divide the above equality by t α and take l = in Lemma 2.(c) to deduce v ψ (x,t) t α = Γ(+α) λ n ψ n t α λ 2 n E α,+2α( λ n t α )ψ n. Again, we take t = t i and pass i to get λ nψ n = and thus v ψ (x,t) = ( λ n t α ) 2 E α,+2α ( λ n t α )ψ n. Repeating the same process, we can show by induction that v ψ (x,t) = ( λ n t α ) l E α,+lα ( λ n t α )ψ n ( l =,,2,...). (3.) Actually, suppose that (3.) holds for some l N. For t >, we divide (3.) by ( t α ) l and apply Lemma 2.(c) to deduce v ψ (x,t) ( t α ) l = Γ(+lα) λ l n ψ n t α λ l+ n E α,+(l+)α( λ n t α )ψ n, where the boundedness of the involved summations follows from (3.9) (3.). Taking t = t i and passing i, we obtain λl nψ n = and thus v ψ (x,t) = ( λ n t α ) l+ E α,+(l+)α ( λ n t α )ψ n. Now it suffices to prove lim l ( λ n t α ) l E α,+lα ( λ n t α )ψ n = ( t ). (3.2)

10 Y. Liu, W. Rundell and M. Yamamoto In fact, writing η = λ n t α, it turns out that ( λ n t α ) l E α,+lα ( λ n t α ) = ( η) l k= ( η) k Γ(αk ++lα) = ( η) k Γ(αk +) coincides with the summation after the lth term in the series by which E α, ( η) is defined. Noting that the series is uniformly convergent with respect to η, we obtain k=l lim ( λ nt α ) l E α,+lα ( λ n t α ) = ( n =,2,..., t ), l which, together with the boundedness of ψ n, yields (3.2) immediately. Since v ψ (x,t) is independent of l, by (3.) we eventually conclude v ψ (x,t) = E α, ( λ n t α )ψ n = (t, ψ C (ω)). (3.3) Recall that the Laplace transform of E α, ( λ n t α ) reads (see, e.g., Podlubny [23, p.2]) e zt E α, ( λ n t α )dt = zα z α +λ n (Rez > λ /α n ), which is analytically extended to Rez >. Since the series in (3.3) converges in C[, ), we can take the Laplace transform with respect to t in (3.3) to derive that is, z α ψ n z α +λ n = (Rez >, ψ C (ω)), ψ n ζ +λ n = (Reζ >, ψ C (ω)). (3.4) Byasimilarargumentfortheconvergenceof (3.7),weseethattheaboveseriesisalsoconvergent in any compact set in C\{ λ n }, and the analytic continuation in ζ yields that (3.4) holds for ζ C\{ λ n }. Especially, since the first eigenvalue λ is single, we can choose a small circle around λ which does not contain λ n (n 2). Integrating (3.4) on this circle yields ψ = (ϕ,ψ)ϕ (x ) = ( ψ C (ω)). Since ψ C (ω) is arbitrarily chosen, there should hold ϕ (x )ϕ = a.e. in ω. However, this contradicts with the strict positivity of the first eigenfunction ϕ (see, e.g., Evans [6]). Therefore, t = cannot be an accumulation point of E x. In summary, for any x Ω, we have excluded all the possibilities for E x to possess any accumulation point, indicating that E x is at most a finite set. Taking advantage of the Green function introduced in (3.2), it is straightforward to demonstrate Corollary.3. Proof of Corollary.3. Recall that the solution v a allows a pointwise definition if a L 2 (Ω), and v a in Ω (, ) by Lemma 2.3. Assume contrarily that there exists x Ω and t >

11 Strong Maximum Principle for Fractional Diffusion Equations such that v g (x,t ) =. Employing the representation (3.3), we see Ω G(x,y,t )a(y)dy =. Since G(x,,t ) by Lemma 3. and a >, there should be G(x,,t ) =, that is, E α, (λ n t α )ϕ n (x )ϕ n = in Ω. Since {ϕ n } is a complete orthonormal basis in L 2 (Ω), we obtain E α, ( λ n t α )ϕ n(x ) = for all n =,2,..., especially, E α, ( λ t α )ϕ (x ) =. However, it is impossible because E α, ( λ t α ) > according to Lemma 2.(a) and ϕ (x ) >. Therefore, such a pair (x,t ) cannot exist and we complete the proof. 4 Proof of Theorem.5 Now we turn to the proof of Theorem.5, i.e., the uniqueness for Problem.4. Recall that the initial-boundary value problem under consideration is t α u(x,t)+au(x,t) = ρ(t)g(x) (x Ω, < t T), u(x,) = (x Ω), (4.) u(x,t) = (x Ω, < t T), where ρ C [,T], g D(A ε ) with some ε >, g and g. First we establish the following fractional Duhamel s principle for the fractional diffusion equation. Lemma 4. Let u be the solution to (4.), where ρ C [,T] and g D(A ε ) with some ε >. Then u allows the representation u(,t) = µ(t s)v g (,s)ds ( < t T), where v g solves the homogeneous problem (3.) with g as the initial data, and µ(t) := d Γ(α) dt ρ(s) ds ( < t T). (4.2) (t s) α Proof. Henceforth C > denotes generic constants independent of the choice of g. First, since ρg L (,T;L 2 (Ω)), Lemma 2.2(b) indicates that u L 2 (,T;H 2 (Ω) H (Ω)) and lim t u(,t) L 2 (Ω) =. By setting we shall demonstrate ũ(,t) := µ(t s)v(,s)ds, (4.3) u = ũ in L 2 (,T;H 2 (Ω) H (Ω)), lim ũ(,t) L t 2 (Ω) =. Since ρ C [,T], simple calculations for (4.2) yield ( µ(t) = d ρ(s)(t s)α s=t Γ(α) dt α + ) t (t s) α ρ (s)ds s= α = ( d ρ() Γ(α) dt α tα + ) (t s) α ρ (s)ds α

12 2 Y. Liu, W. Rundell and M. Yamamoto = ( ρ() t Γ(α) t α + ρ ) (s) ds, (4.4) (t s) α which further implies µ L (,T), µ(t) Ct α ( < t T). (4.5) Regarding the solution v g to (3.), since g L 2 (Ω), the application of estimate (2.3) in Lemma 2.2(a) yields v g (,t) L 2 (Ω) C g L 2 (Ω), v g (,t) H 2 (Ω) C g L 2 (Ω)t α ( < t T). By the relation (4.3) and the estimate in (4.5), for < t T we estimate ũ(,t) L2 (Ω) ũ(,t) H 2 (Ω) indicating µ(t s) v g (,s) L2 (Ω)ds C g L2 (Ω) µ(t s) v g (,s) H 2 (Ω)ds C g L 2 (Ω) C g L 2 (Ω), s α ds C g L2 (Ω)t α, (t s) α s α ds ũ L (,T;H 2 (Ω) H (Ω)) L 2 (,T;H 2 (Ω) H (Ω)), lim t ũ(,t) L 2 (Ω) =. On the other hand, according to the explicit representation (2.2) and Lemma 2.(d), we obtain t v g (,t) = t α E α,α ( λ n t α )(g,ϕ n )ϕ n. By Lemma 2.(b) and the fact g D(A ε ) with ε >, we estimate for < t T that t v g (,t) L 2 (Ω) t 2(α ) E α,α ( λ n t α )(g,ϕ n ) 2 = t 2(α ) λ ε n E α,α( λ n t α ) 2 λ ε n (g,ϕ n) 2 ( Ct αε ) 2 (λ n t α ) ε 2 +λ n t α λ ε n (g,ϕ n) 2 ( C g D(Aε )t αε ) 2. (4.6) In order to show u = ũ, it suffices to verify that ũ also solves the initial-boundary value problem (4.) which possesses a unique solution (see Lemma 2.2(b)). To calculate t α ũ, first we formally calculate t ũ(,t) = t µ(s)v g (,t s)ds = Then we employ (4.5) and (4.6) to estimate t ũ(,t) L 2 (Ω) µ(t s) s v g (,s) L 2 (Ω)ds+ µ(t) g L 2 (Ω) µ(s) t v g (,t s)ds+µ(t)g. (4.7)

13 Strong Maximum Principle for Fractional Diffusion Equations 3 C g D(Aε ) (t s) α s αε ds+c g L2 (Ω)t α C g D(Aε )t α(+ε) +C g L2 (Ω)t α C g D(Aε )t α ( < t T), implying that the above differentiation makes sense in L 2 (Ω) for < t T. By definition, we have t α ũ(,t) = Γ( α) I := I 2 := Γ( α) Γ( α) (t s) α s ũ(,s) (t s) α ds = I +I 2 g, where s µ(s) (t s) α ds. µ(τ) s v(,s τ)dτds, The governing equation (3.) for v g and formula (4.4) for µ imply respectively I = Γ( α) ( = µ(τ) I 2 = = µ(τ) Γ( α) s v(,s τ) (t s) α dsdτ ) s v(,,s) ((t τ) s) α ds dτ = τ τ µ(τ)av(,t τ)dτ = A Γ( α)γ(α) = Γ( α)γ(α) = ρ()+ ( ρ() (t s) α ( ρ() s α + µ(τ) α t v(,t τ)dτ µ(τ)v(,t τ)dτ = Aũ(,t), s ds (t s) α s α + ρ (τ) (s τ) ρ (τ) dτ α τ ) ds ds (t s) α dτ (s τ) α ρ (τ)dτ = ρ(t). (4.8) Therefore, we conclude α t ũ+aũ = ρg and the proof is completed. At this stage, we can proceed to show Theorem.5 by applying the established strong maximum principle. Completion of Proof of Theorem.5. Let the conditions in the statement of Theorem.5 be valid, namely, it is assumed that ρ C [,T], g D(A ε ) with ε >, g, g, and the solution u to (4.) vanishes in {x } [,T] for some x Ω. According to the fractional Duhamel s principle proved above, we have u(x,t) = µ(t s)v g (x,s)ds = ( t T), where µ was defined in (4.2) and v g solves (3.) with the initial data g. Now the estimate (2.3) in Lemma 2.2(a) and the Sobolev embedding indicate v g (x,t) C v g (,t) H2 (Ω) C g L2 (Ω)t α ( < t T) and thus v g (x, ) L (,T). Meanwhile, (4.5) guarantees µ L (,T). Therefore, the Titchmarshconvolutiontheorem(see[28])impliestheexistenceofT,T 2 satisfyingt +T 2 )

14 4 Y. Liu, W. Rundell and M. Yamamoto T such that µ(t) = for almost all t (,T ) and v g (x,t) = for all t [,T 2 ]. However, since the initial data g of (3.) satisfies g and g, Theorem. asserts that v g (x, ) only admits at most a finite number of zero points, i.e., v g (x, ) > a.e. in (,T). As a result, the only possibility is T 2 = and thus T = T, that is, µ = a.e. in (,T). Finally, it suffices to utilize the following reverse formula t µ(s) ρ(t) = Γ( α) (t s) α ds, which was obtained in (4.8). We apply Young s inequality to conclude t ρ L (,T) = µ(s) Γ( α) (t s) α ds T α Γ(2 α) µ L (,T) =, which finishes the proof. References L (,T) [] Adams E E and Gelhar L W 992 Field study of dispersion in a heterogeneous aquifer 2. Spatial moments analysis Water Resources Res [2] Adams R A 975 Sobolev Spaces (New York: Academic Press) [3] Al-Refai M and Luchko Y 24 Maximum principles for the fractional diffusion equations with the Riemann.Liouville fractional derivative and their applications Fract. Calc. Appl. Anal [4] Cannon J R and Esteva S P 986 An inverse problem for the heat equation Inverse Problems [5] Courant R and Hilbert D 953 Methods of Mathematical Physics (New York: Interscience) [6] Evans L C 2 Partial Differential Equations 2nd ed (Providence, RI: American Mathematical Society) [7] Gilbarg D and Trudinger N S 2 Elliptic Partial Differential Equations of Second Order (Berlin: Springer) [8] Ginoa M, Cerbelli S amd Roman H E 992 Fractional diffusion equation and relaxation in complex viscoelastic materials Phys. A [9] Gorenflo R, Kilbas A A, Mainardi F and Rogosin S V 24 Mittag-Leffler Functions, Related Topics and Applications (Berlin: Springer) [] Gorenflo R, Luchko Y and Yamamoto M 25 Time-fractional diffusion equation in the fractional Sobolev spaces Fract. Calc. Appl. Anal [] Hatano Y and Hatano N 998 Dispersive transport of ions in column experiments: an explanation of long-tailed profiles Water Resources Res [2] Jin B, Lazarov R, Pasciak J and Zhou Z 25 Error analysis of semidiscrete finite element methods for inhomogeneous time-fractional diffusion IMA J. Numer. Anal [3] Jin B, Lazarov R and Zhou Z 23 Error estimates for a semidiscrete finite element method for fractional order parabolic equations SIAM J. Numer. Anal [4] Jin B and Rundell W 25 A tutorial on inverse problems for anomalous diffusion processes Inverse Problems [5] Liu F, Zhuang P, Anh V, Turner I and Burrage K 27 Stability and convergence of the difference methods for the space-time fractional advection-diffusion equation Appl. Math. Comput [6] Luchko Y 29 Maximum principle for the generalized time-fractional diffusion equation J. Math. Anal. Appl [7] Luchko Y 2 Some uniqueness and existence results for the initial-boundary value problems for the generalized time-fractional diffusion equation Comput. Math. Appl [8] Luchko Y and Gorenflo R 999 An operational method for solving fractional differential equations with the Caputo derivatives Acta Math. Vietnam [9] Luchko Y, Rundell W, Yamamoto M and Zuo L 23 Uniqueness and reconstruction of an unknown semilinear term in a time-fractional reaction.diffusion equation Inverse Problems

15 Strong Maximum Principle for Fractional Diffusion Equations 5 [2] Meerschaert M M and Tadjeran C 24 Finite difference approximations for fractional advection-dispersion flow equations J. Comput. Appl. Math [2] Metzler R and Klafter J 2 The random walk s guide to anomalous diffusion: a fractional dynamics approach Phys. Rep [22] Nigmatullin R R 986 The realization of the generalized transfer equation in a medium with fractal geometry Phys. Stat. Sol. B [23] Podlubny I 999 Fractional Differential Equations (San Diego: Academic Press) [24] Rundell W, Xu X and Zuo L 23 The determination of an unknown boundary condition in a fractional diffusion equation Appl. Anal [25] Saitoh S, Tuan V K and Yamamoto M 22 Reverse convolution inequalities and applications to inverse heat source problems JIPAM. J. Inequal. Pure Appl. Math. 3 Article 8 [26] Saitoh S, Tuan V K and Yamamoto M 23 Convolution inequalities and applications JIPAM. J. Inequal. Pure Appl. Math. 4 Article 5 [27] Sakamoto K and Yamamoto M 2 Initial value/boundary value problems for fractional diffusion-wave equations and applications to some inverse problems J. Math. Anal. Appl [28] Titchmarsh E C 926 The zeros of certain integral functions Proc. London Math. Soc