arxiv: v1 [math.ap] 8 Jan 2017

Transcription

1 NON-LOCAL INITIAL PROBLE FOR SECOND ORDER TIE-FRACTIONAL AND SPACE-SINGULAR EQUATION ERKINJON KARIOV, URAT ACHUEV, AND ICHAEL RUZHANSKY arxiv:7.94v [math.ap] 8 Jan 27 Abstract. In this wor, we consider an initial problem for second order partial differential equations with Caputo fractional derivatives in the time-variable and Bessel operator in the space-variable. For non-local boundary conditions, we present a solution of this problem in an explicit form representing it by the Fourier-Bessel series. The obtained solution is written in terms of multinomial ittag-leffler functions and first ind Bessel functions.. Introduction and formulation of a problem It is well-nown that partial differential equations are playing a ey role in constructing mathematical models for many real-life processes. Especially, during the last decades, many applications of various inds of fractional differential equations became target of intensive research due to both theoretical and practical reasons, see e.g. [6] for an overview. any inds of boundary problems, including direct [] and inverse problems [4], were formulated for different type of PDEs of integer order including several differential operators of fractional order. We note wors [3]-[9] devoted to studying partial differential equations with multiple Caputo derivatives. Precisely, in [9] the authors studied fractional differential equations with Caputo fractional derivatives and using the operational method, solutions of initial boundary problem for those equations were obtained in an explicit form involving a multinomial ittag-leffler function. Certain properties of this function were obtained by Li, Liu and Yamamoto [7] and applied to studying initial-boundary problems for time-fractional diffusion equations with positive constant coefficients. Later, Liu [8] established strong maximum principle for fractional diffusion equations with multiple Caputo derivatives and investigated a related inverse problem. Daftardar-Gejji and Bhaliar [3], using the method of separation of variables, solved some boundary-value problems for multi-term fractional diffusion-wave equation. We also note wors related to the Bessel operator. In [], the initial inverse problem for the heat equation with Bessel operator was investigated. Inverse initial and inverse source problems for time-fractional diffusion equation with zero order Bessel operator were recently studied in [2]. Direct and inverse problems for PDEs containing two-term time fractional Caputo derivatives of orders up to, and Bessel operator of order ν were investigated in []. Date: February 5, athematics Subject Classification. 35R, 33E2. Key words and phrases. fractional derivatives, Cauchy problem, Bessel operator. The last author was supported in parts by the EPSRC grant EP/K3947/ and by the Leverhulme Grant RPG No new data was collected or generated during the course of research.

2 2 ERKINJON KARIOV, URAT ACHUEV, AND ICHAEL RUZHANSKY The consideration of non-local initial conditions is often justified by practical usage in certain real-life processes. For instance, when the initial temperature for the heat equation is not given instantly, but there is an information related with the temperature on a certain time interval that can be described by a non-local initial condition in a simple form. Boundary-value problems with non-local initial conditions were considered in wors [2] for reaction-diffusion equations, in [6] for heat equation, in [3]-[4] for degenerate parabolic equations, and in [5] for a mixed parabolic equation. For integer orders much more is nown, and for a review of different questions of time decay of solutions for hyperbolic equations with integer order derivatives we can refer to [5]. In the present wor we deal with the non-local initial boundary problem for multiterm time fractional PDE with Bessel operator of order ν. We use Fourier-Bessel series expansion in order to find the explicit solution for the considered problem, yielding also its existence. Because of the singularities in the Bessel operator such conditions appear naturally in space variables. We note that most of the current literature deals with diffusion type equations considering time-derivatives of orders up to, see e.g. [7]. One of the novelties of the present paper is that we consider wave type equations allowing fractional derivatives of order up to 2, with additional fractional dissipation type terms. If there are multiple fractional time-derivative terms, a multinomial ittag-leffler function appears in the representation of solutions. Let us now describe the problem in more detail. We consider the equation L(u) B ν (u) = f(t,x) (.) in a rectangular domain D = {(x,t) : < x <, < t < T}, T >, where f(t,x) is a given function, n L(u) = tu(t,x) α λ i α i t u(t, x) (.2) i= is the time component of the equation, with orders and < α i, α i α 2, n N, λ i R, B ν (u) = u xx (t,x)+ x u x(t,x) ν2 x2u(t,x) (.3) is the Bessel part of the equation with ν >. Here t g () (z) Γ( α) (t z) α +dz, α / N, β t g(t) = d g(t), α N, dt g(t), α =, is a fractional differential operator of Caputo type, where = [α]+, and [α] is the integer part of α. We can refer to [6] for further details on the Caputo fractional derivative operators.

3 INITIAL PROBLE FOR TIE-FRACTIONAL AND SPACE-SINGULAR EQUATION 3 The non-local initial boundary problem for equation (.)-(.3) is formulated as follows: Problem: Let R. To find a solution u(t,x) of equation (.)-(.3) in D, which satisfies (i) regularity conditions u W with W = u(t,x) : u C(D), u xx, t α u C(D), x u(t,x) dx < + ; (.4) (ii) boundary and non-local initial conditions lim x(t,x) =, u(t,) =, x (.5) u(,x)+u(t,x) =, x, [α] u t (,t) =, < x <. (.6) 2. ain result The main result of this note is the following well-posedness theorem for the initial problem (.)-(.6). The interesting part are the conditions on f allowing one to handle the singularities in the coefficients of the Bessel operator, and the nonresonance conditions (2.) relating the parameter with the length T of the time interval, coefficients and fractional orders of time-derivatives, through the multinomial ittag-leffler function. Theorem 2.. Assume that f(x,t) is differentiable four times with respect to x; f(,t) = f (,t) = f (,t) = f (,t) =, f(,t) = f (,t) = f (,t) = ; 4 f(x,t) is bounded; x f(x, t) is continuous and continuously differentiable with respect to t, and non-resonance conditions (2.) E (α α,α α 2,...,α),(λ T α α,...,λn T α αn, γ 2Tα ) hold for all =,2,... Then there exists a unique solution of the problem (.)- (.6), and it can be written in the following form: u(t,x) = = t,...,λ n z α αn, γ 2 zα )f (t z)dz +E (α α,α α 2,...,α),(λ T α α,...,λn T α αn, γ 2 Tα ) T,...,λ n z α αn, γ 2 zα )f (T z)dz+ +E (α α,α α 2,...,α),(λ t α α,...,λ n t α αn, γ 2 t α ) ] J ν (γ x). (2.2)

4 4 ERKINJON KARIOV, URAT ACHUEV, AND ICHAEL RUZHANSKY The numbers γ and the functions appearing in the formula (2.2) are explained in Section 2.. Briefly, here E (α α,α α 2,...,α α n,α),α( ) is the multinomial ittag-leffler function, J ν are the first ind Bessel functions, γ are their positive zeros, λ are coefficients in the operator (.2), and f are Bessel expansions of f. For the proof of Theorem 2. we start finding a formal solution in a series form and the convergence of the appearing series will be shown in Section Representation of a solution. We start by a formal discussion of the representation of solutions in (2.2). Let ( ) i (z/2) 2i+ν J ν (z) = i!(i+ν)! i= be the Bessel function of the first ind (see e.g. [8]). It is nown that for ν >, the Bessel function J ν (z) has countably many zeros, moreover, they are real and have pairwise opposite signs. Denote the th positive root of the equation J ν (z) = by γ, =,2,... For large, we have (see [7]) γ π + νπ 2 π 4. We now expand functions u(t,x) and f(t,x) in the Fourier-Bessel series (see e.g. [7]), writing them in the form where U (t) = u(t,x) = f(t,x) = U (t)j ν (γ x), (2.3) = f (t)j ν (γ x), (2.4) = 2 u(x,t)xj Jν+(γ 2 ν (γ x)dx, (2.5) ) 2 f (t) = Jν+ 2 (γ f(x,t)xj ν (γ x)dx. (2.6) ) It is nown that if a function g = g(x) is piecewise continuous on [,l] and satisfies l x g(x) dx < +, then for ν > /2, the Fourier-Bessel series converges at every point x (,l), see e.g. [7]. Since we are looing for a function u(t,x) W, it satisfies these required conditions in order to be represented by a Fourier-Bessel series. We substitute (2.3)-(2.4) into equation (.) and obtain the eigenvalue problem L(U )+γ 2 U (t) = f (t). (2.7)

5 INITIAL PROBLE FOR TIE-FRACTIONAL AND SPACE-SINGULAR EQUATION 5 According to [9], the solution for equation (2.7) satisfying initial conditions U () = A, [α]u () =, (2.8) can be represented in the form U (t) = t,...,λ n z α αn, γ 2 z α ) f (t z)dz +AU (t), (2.9) where U (t) = + n λ i t α α i E (α α,α α 2,...,α α n,α),+α α i (λ t α α,...,λ n t α αn, γ 2 tα ) i= γ 2 t α E (α α,α α 2,...,α α n,α),+α(λ t α α,...,λ n t α αn, γ 2 t α ), (2.) and E (a,a 2,...,a n),b(z,z 2,...,z n ) = = l +l l n = l,...,l n n z l i i! i= l!...l n! Γ(b+ n a i l i ) i= (2.) is the multinomial ittag-leffler function ([9]). From (2.9) we find that U (T) = T,...,λ n z α αn, γ 2 zα ) f (T z)dz +AU (T). (2.2) Considering U () = A, from the first condition in (.5), we get the relation Using (2.2) we find A to be A+U (T) =. A = +U (T) T,...,λ n z α αn, γ 2 zα ) f (T z)dz.

6 6 ERKINJON KARIOV, URAT ACHUEV, AND ICHAEL RUZHANSKY Substituting this value of A into (2.9), we rewrite it as U (t) = t U (t) +U (T),...,λ n z α αn, γ 2 z α )f (t z)dz T If we use the formula (see [8]),...,λ n z α αn, γ 2 zα ) f (T z)dz. (2.3) n+ + λ j t α α j E (α α,α α 2,...,α α n+ ),+α α j (λ t α α,...,λ n+ t α α n+ ) j= = E (α α,α α 2,...,α α n+ ),(λ t α α,...,λ n+ t α α n+ ), representation (2.) can be rewritten as Denoting F (t) = t U (t) = E (α α,α α 2,...,α α n+ ),(λ t α α,...,λ n t α αn, γ 2 tα ).,...,λ n z α αn, γ 2 zα )f (t z)dz, U (t) = E (α α,α α 2,...,α),(λ t α α,...,λ n t α αn, γ 2 t α ), (2.4) we rewrite the function (2.3) as U (t) = F (t) +U (T) F (T)U (t). (2.5) We note that the above expression is well-defined in view of the non-resonance conditions (2.), that is U (T) = E (α α,α α 2,...,α),(λ T α α,...,λn T α αn, γ 2 Tα ) holds for all. Finally, based on (2.5), we rewrite our formal solution as u(t,x) = = [ ] F (t) +U (T) F (T)U (t) J ν (γ x). (2.6) 2.2. Justification of formal solution. In this section we prove convergence of the obtained infinite series corresponding to functions u(t,x), u xx (t,x) and α tu(t,x). In order to prove the convergence of these series, we use the estimate for the ittag-leffler function (2.), obtained in [7, Lemma 3.2], of the form E (α α,α α 2,...,α α n),ρ(z,z 2,...,z n ) C + z. (2.7)

7 INITIAL PROBLE FOR TIE-FRACTIONAL AND SPACE-SINGULAR EQUATION 7 Let us first prove the convergence of series (2.6). For this, we collect several other nown estimates. First, we use the following theorem on the estimate of the Fourier-Bessel coefficient: Theorem 2.2 ([7, p. 23]). Let f(x) be a function defined on the interval [,] such that f(x) is differentiable 2s times (s N) and such that () f() = f () =... = f (2s ) () = ; (2) f (2s) (x) is bounded (this derivative may not exist at certain points); (3) f() = f () =... = f (2s 2) () =. Then the following inequality is satisfied by the Fourier-Bessel coefficients of f(x): c f γ 2s /2 where γ is the th positive zero of of the function J ν (x). In our case we have f (t) = F (t) (c = const), (2.8) +U (T) F (T)U (t). According to(2.4), using estimate(2.7) and imposing conditions()-(3) of Theorem (2.2) in the case s = on the function f(t,x), we get F (t) +U (T) F (T)U (t) C. (2.9) γ 3/2 The series (2.5) then converges absolutely and uniformly on [,] in view of the following theorem: Theorem 2.3 ([7, p. 225]). If ν, C >, and if the constants c satisfy for some ε >, then the series c C γ +ε c J ν (γ x)+c 2 J ν (γ 2 x)+c 3 J ν (γ 3 x)+... converges absolutely and uniformly on [, ]. According to [7, Theorem 2, p. 236], sufficient condition for differentiating the series (2.3) twice term by term, i.e. for the validity of the relation u xx (t,x) = U (t)γ 2 J ν (γ x) (2.2) will be = U (t) C 2, γ 5/2+ǫ, (2.2) where U are as in (2.5), and C 2 and ǫ are positive constants. Based on this estimation and considering also (see [7, p. 233]) J ν (γ x) C 3 γ x, (C 3 = const),

8 8 ERKINJON KARIOV, URAT ACHUEV, AND ICHAEL RUZHANSKY we deduce (see e.g. [7, p. 236]) U (t)γj 2 ν(γ x) C 2 γ 2+ǫ x + C3( ν +ν2 ) 2 γ 3+ǫ x 2 x + C 4 γ +ǫ, x which provides the convergence of the series (2.2). We note that if we impose conditions on the given function f(x,t) of the form f(x,t) is differentiable four times with respect to x; f(,t) = f (,t) = f (,t) = f (,t) =, f(,t) = f (,t) = f (,t) = ; is bounded, 4 f(x,t) x one can see that Theorem 2.3 implies estimation (2.2). The convergence of series corresponding to α t u(x,t), u x(x,t) can be shown in similar ways, completing the proof of Theorem 2.. References [] P. Agarwal, E. Karimov,. amchuev and. Ruzhansy, On boundary-value problems for a partial differential equation with Caputo and Bessel operators, arxiv:6.624 [2] F. Al-usalhi, N. Al-Salti and S. Kerbal, Inverse problems of a fractional differential equation with Bessel operator, arxiv: v [3] V. Daftardar-Gejji, S. Bhalear, Boundary value problems for multi-term fractional differential equations, J. ath. Anal. Appl., 345 (28), [4] V. Isaov, Inverse problems for partial differential equations (Second edition), Springer, New Yor, 26. [5] A. A. Kerefov. The Gevrey problem for a certain mixed-parabolic equation. (Russian) Diff. Uravn., 3() (977), [6] A. A. Kilbas, H.. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Elsevier, Amsterdam, 26. [7] Z. Li, Y. Liu and. Yamamoto, Initial-boundary value problems for multi-term time-fractional diffusion equations with positive constant coefficients, Appl. ath. Comput. 257 (25), [8] Y. Liu, Strong maximum principle for multi-term time-fractional diffusion equations and its application to an inverse source problem. arxiv:5.6878v [9] Y. Lucho, R. Gorenflo, An operational method for solving fractional differential equations with the Caputo derivatives, Acta ath. Vietnam., 24 (999), [] K. asood, S. essaoudi and F. D. Zaman, Initial inverse problem in heat equation with Bessel operator, International Journal of Heat and ass Transfer, 45 (22), [] A.. Nahushev, Drobniye ischisleniya i yego primenenie, in: Fractional Calculus and its Applications. Fizmatlit, oscow, 23 (in Russian). [2] C. V. Pao, Reaction diffusion equations with nonlocal boundary and nonlocal initial conditions, J. ath. Anal. Appl., 95 (995), [3] J.. Rassias, E. T. Karimov, Boundary-value problems with non-local initial condition for degenerate parabolic equations, Contemp. Anal. Appl. ath., () (23), [4] J.. Rassias, E. T. Karimov, Boundary-value problems with non-local initial condition for parabolic equations with parameter, European J. Pure Appl. ath., 3(6) (2), [5]. Ruzhansy, J. Smith, Dispersive and Strichartz estimates for hyperbolic equations with constant coefficients, SJ emoirs, 22, athematical Society of Japan, Toyo, 2. [6] N. N. Shopolov, ixed problem with non-local initial condition for a heat conduction equation, Reports of Bulgarian Academy of Sciences, 3(7) (98), [7] G. P. Tolstov, Fourier series, Prentice Hall, Inc., Englewood Cliffs, N.J., 962. [8] G. N. Watson, A treatise on the theory of Bessel functions, Cambridge University Press, Second Edition, 966.

9 INITIAL PROBLE FOR TIE-FRACTIONAL AND SPACE-SINGULAR EQUATION 9 Erinjon Karimov: Institute of athematics National University of Uzbeistan Tashent, 25 Uzbeistan address erinjon@gmail.com urat amchuev: Department of Theoretical and athematical Physics Institute of Applied athematics and Automation Shortanova str. 89-A, Nalchi, 36 Kabardino-Balar Republic Russia address mamchuev@rambler.ru ichael Ruzhansy: Department of athematics Imperial College London 8 Queen s Gate, London, SW7 2AZ United Kingdom address m.ruzhansy@imperial.ac.u