COMPLETELY MONOTONE FUNCTIONS IN THE STUDY OF A CLASS OF FRACTIONAL - PowerPoint PPT Presentation

COMPLETELY MONOTONE FUNCTIONS IN THE STUDY OF A CLASS OF FRACTIONAL EVOLUTION EQUATIONS Emilia Bazhlekova Joint Seminar of Analysis, Geometry and Topology Department Institute of Mathematics and Informatics Bulgarian Academy of Sciences

Completely monotone functions ( CMF ) Bulgarian contributions: N. Obreshkov, Y. Tagamlitski, Bl. Sendov, H. Sendov A function f : (0 , ∞ ) → R is called completely monotone if it is of class C ∞ and ( − 1) n f ( n ) ( t ) ≥ 0 , for all t > 0 , n = 0 , 1 , ... Elementary examples: � b + µt − 1 � t − 1 ; e f ( t ) , f ∈ CMF ; e − λt ; ( λ + µt ) − ν ; ln ; where λ, µ, ν > 0 , b ≥ 1 . Bernstein’s theorem: f ( t ) ∈ CMF iff � ∞ e − tx dg ( x ) , f ( t ) = 0 where g ( x ) is nondecreasing and the integral converges for 0 < t < ∞ . Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 2/26

Bernstein functions ( BF ) and some useful properties A C ∞ function f : (0 , ∞ ) → R is called a Bernstein function if f ( t ) ≥ 0 and f ′ ( t ) ∈ CMF . Proposition: (a) The class CMF is closed under pointwise addition and multiplication; The class BF is closed under pointwise addition, but, in general not under multiplication; (b) If f ∈ CMF and ϕ ∈ BF , then the composite function f ( ϕ ) ∈ CMF ; (c) If f ∈ BF , then f ( t ) /t ∈ CMF ; (d) Let f ∈ L 1 loc ( R + ) be a nonnegative and nonincreasing function, such that lim t → + ∞ f ( t ) = 0 . Then ϕ ( s ) = s � f ( s ) ∈ BF ; loc ( R + ) and f ∈ CMF , then � (e) If f ∈ L 1 f ( s ) admits analytic extension to the sector | arg s | < π and | arg � f ( s ) | ≤ | arg s | , | arg s | < π. Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 3/26

The operators of fractional integration and differentiation J α t - the Riemann-Liouville fractional integral of order α > 0 : � t 1 J α ( t − τ ) α − 1 f ( τ ) dτ, t f ( t ) := α > 0 , Γ( α ) 0 where Γ( · ) is the Gamma function. D α t - the Riemann-Liouville fractional derivative C D α t - the Caputo fractional derivative D 1 t = C D 1 t = J 1 − α D 1 t = D 1 t J 1 − α C D α D α t = d/dt ; t , , α ∈ (0 , 1) . t t Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 4/26

Mittag-Leffler function Fractional relaxation equation ( λ > 0 , 0 < α ≤ 1 ): C D α t u ( t ) + λu ( t ) = f ( t ) , t > 0 , u (0) = c 0 . The solution is given by: � t τ α − 1 E α,α ( − λτ α ) f ( t − τ ) dτ. u ( t ) = c 0 E α ( − λt α ) + 0 Mittag-Leffler function ( α, β ∈ R , α > 0 ): � ∞ ( − t ) k E α,β ( − t ) = Γ( αk + β ) , E α ( − t ) = E α, 1 ( − t ) . k =0 E 1 ( − t ) = e − t ∈ CMF E α ( − t ) ∈ CMF , iff 0 < α < 1 (Pollard, 1948) E α,β ( − t ) ∈ CMF , iff 0 ≤ α ≤ 1 , α ≤ β (Schneider, 1996; Miller, 1999) Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 5/26

Plots of E α ( − t α ) for different values of α Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 6/26

Plots of t α − 1 E α,α ( − t α ) for different values of α Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 7/26

Fractional evolution equation of distributed order Two alternative forms: � 1 µ ( β ) C D β t u ( t ) dβ = Au ( t ) , t > 0 , (1) 0 and � 1 µ ( β ) D β u ′ ( t ) = t Au ( t ) dβ, t > 0 , (2) 0 A - closed linear unbounded operator densely defined in a Banach space X Initial condition: u (0) = a ∈ X Reference: E. Bazhlekova, Completely monotone functions and some classes of fractional evolution equations, preprint, 2015, arXiv:1502.04647 Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 8/26

Two cases for the weight function µ : • discrete distribution m � µ ( β ) = δ ( β − α ) + b j δ ( β − α j ) , (3) j =1 where 1 > α > α 1 ... > α m > 0 , b j > 0 , j = 1 , ..., m, m ≥ 0 , and δ is the Dirac delta function; • continuous distribution µ ∈ C [0 , 1] , µ ( β ) ≥ 0 , β ∈ [0 , 1] , (4) and µ ( β ) � = 0 on a set of a positive measure. Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 9/26

Discrete distribution: Multi-term time-fractional equations in the Caputo sense m � α j C D α C D t u ( t ) + b j t u ( t ) = Au ( t ) , t > 0 , (5) j =1 and in the Riemann-Liouville sense m � α j u ′ ( t ) = D α t Au ( t ) + b j D t Au ( t ) , t > 0 (6) j =1 If m = 0 (single-term equations): problem (5) is equivalent to (6) with α replaced by 1 − α . All problems are generalizations of the classical abstract Cauchy problem u ′ ( t ) = Au ( t ) , t > 0; u (0) = a ∈ X. (7) Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 10/26

Solution u ( t ) of (5) with A = − 1 for: m = 1 , α = 0 . 75 , α 1 = 0 . 25 , m = 0 , α = 0 . 25 m = 0 , α = 0 . 75 . Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 11/26

Solution u ( t ) of (5) with A = − 1 for: m = 2 , α = 0 . 75 , α 1 = 0 . 5 , α 2 = 0 . 25 m = 1 , α = 0 . 75 , α 1 = 0 . 25 , m = 0 , α = 0 . 25 m = 0 , α = 0 . 75 . Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 12/26

Unified approach to the four problems Rewrite problems (1) and (2) as an abstract Volterra integral equation � t u ( t ) = a + k ( t − τ ) Au ( τ ) dτ, t ≥ 0; a ∈ X, 0 where k 1 ( s ) = ( h ( s )) − 1 , � � k 2 ( s ) = h ( s ) /s, In the continuous distribution case: � 1 µ ( β ) s β dβ. h ( s ) = 0 In the discrete distribution case: m � h ( s ) = s α + b j s α j . j =1 Define g i ( s ) = 1 / � k i ( s ) , i = 1 , 2 . Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 13/26

Particular cases In the single-term case: k 1 ( t ) = t α − 1 t − α g 1 ( s ) = s α , g 2 ( s ) = s 1 − α , Γ( α ) , k 2 ( t ) = Γ(1 − α ) , In the double-term case: t − α t − α 1 k 1 ( t ) = t α − 1 E α − α 1 ,α ( − b 1 t α − α 1 ) , k 2 ( t ) = Γ(1 − α ) + b 1 Γ(1 − α 1 ) , s g 1 ( s ) = s α + b 1 s α 1 , g 2 ( s ) = s α + b 1 s α 1 = s � k 1 ( s )!!! In the case of continuous distribution in its simplest form: µ ( β ) ≡ 1 . g 1 ( s ) = s − 1 g 2 ( s ) = s log s log s , s − 1 . Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 14/26

Properties of the kernels Theorem. Let µ ( β ) be either of the form (3) or of the form (4) with the additional assumptions µ ∈ C 3 [0 , 1] , µ (1) � = 0 , and µ (0) � = 0 or µ ( β ) = aβ ν as β → 0 , where a, ν > 0 . Then for i = 1 , 2 , : (a) k i ∈ L 1 loc ( R + ) and lim t → + ∞ k i ( t ) = 0 ; (b) k i ( t ) ∈ CMF for t > 0 ; (c) k 1 ∗ k 2 ≡ 1 ; (d) g i ( s ) ∈ BF for s > 0 ; (e) g i ( s ) /s ∈ CMF for s > 0 ; (f) g i ( s ) admits analytic extension to the sector | arg s | < π and | arg g i ( s ) | ≤ | arg s | , | arg s | < π. In the discrete distribution case a stronger inequality holds: | arg g i ( s ) | ≤ α | arg s | , | arg s | < π. Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 15/26

The classical abstract Cauchy problem: u ′ ( t ) = Au ( t ) , t > 0; u (0) = a ∈ X. Main result: Assume that the classical Cauchy problem is well-posed with solution u ( t ) satisfying � u ( t ) � ≤ M � a � , t ≥ 0 . Then any of the problems � 1 µ ( β ) C D β t u ( t ) dβ = Au ( t ) , t > 0 , u (0) = a ∈ X, 0 � 1 µ ( β ) D β u ′ ( t ) = t Au ( t ) dβ, t > 0 , u (0) = a ∈ X 0 is well-posed with solution satisfying the same estimate. Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 16/26

The classical abstract Cauchy problem: u ′ ( t ) = Au ( t ) , t > 0; u (0) = a ∈ X. T ( t ) - solution operator (defined by T ( t ) a = u ( t ) , t ≥ 0 ); R ( s, A ) - resolvent operator of A : � ∞ R ( s, A ) = ( s − A ) − 1 = e − st T ( t ) dt, s > 0 , 0 The Hille-Yosida theorem states that the classical Cauchy problem is well-posed with solution operator T ( t ) such that � T ( t ) � ≤ M, t ≥ 0 iff R ( s, A ) is well defined for s ∈ (0 , ∞ ) and � R ( s, A ) n � ≤ M/s n , s > 0 , n ∈ N . Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 17/26

Abstract Volterra integral equation � t u ( t ) = a + k ( t − τ ) Au ( τ ) dτ, t ≥ 0; a ∈ X, 0 The Laplace transform of the solution operator S ( t ) � ∞ e − st S ( t ) dt, H ( s ) = s > 0 0 is given by H ( s ) = g ( s ) g ( s ) = 1 / � s R ( g ( s ) , A ) , k ( s ) . The Generation Theorem (Pruss, 1993) states that the integral equation is well- posed with solution operator S ( t ) satisfying � S ( t ) � ≤ M, t ≥ 0 , iff � H ( n ) ( s ) � ≤ M n ! s n +1 , for all s > 0 , n ∈ N 0 . Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 18/26

Main result Theorem. Suppose that the classical Cauchy problem is well-posed with solution u ( t ) satisfying � u ( t ) � ≤ M � a � , t ≥ 0 . Then problems (1) and (2) are well-posed and their solutions satisfy the same estimate. Proof: We know � R ( s, A ) n � ≤ M/s n , s > 0 , n ∈ N . We have to prove � H ( n ) ( s ) � ≤ M n ! s n +1 , for all s > 0 , n ∈ N 0 , where H ( s ) = g ( s ) s R ( g ( s ) , A ) , and g ( s ) = 1 / � k ( s ) , R ( s, A ) = ( s − A ) − 1 . Joint Seminar of Analysis, Geometry and Topology Dept. 24.03.2015, p. 19/26