Next Article in Journal
Left Riemann–Liouville Fractional Sobolev Space on Time Scales and Its Application to a Fractional Boundary Value Problem on Time Scales
Previous Article in Journal
Existence and U-H Stability Results for Nonlinear Coupled Fractional Differential Equations with Boundary Conditions Involving Riemann–Liouville and Erdélyi–Kober Integrals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hadamard-Type Fractional Integro-Differential Problem: A Note on Some Asymptotic Behavior of Solutions

by
Ahmad Mugbil
1,* and
Nasser-Eddine Tatar
2
1
Prep Year Math Program, College of General Studies, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia
2
Department of Mathematics, IRC for Intelligent Manufacturing and Robotics, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(5), 267; https://doi.org/10.3390/fractalfract6050267
Submission received: 30 March 2022 / Revised: 10 May 2022 / Accepted: 13 May 2022 / Published: 15 May 2022
(This article belongs to the Topic Fractional Calculus: Theory and Applications)

Abstract

:
As a follow-up to the inherent nature of Hadamard-Type Fractional Integro-differential problem, little is known about some asymptotic behaviors of solutions. In this paper, an integro-differential problem involving Hadamard fractional derivatives is investigated. The leading derivative is of an order between one and two whereas the nonlinearities may contain fractional derivatives of an order between zero and one as well as some non-local terms. Under some reasonable conditions, we prove that solutions are asymptotic to logarithmic functions. Our approach is based on a generalized version of Bihari–LaSalle inequality, which we prove. In addition, several manipulations and crucial estimates have been used. An example supporting our findings is provided.

1. Introduction

Of concern is the following general class of initial value problems modelled by:
H D t 0 α u ( t ) = f t , H D t 0 α 1 u ( t ) , t 0 t h t , s , H D t 0 α 2 u ( s ) d s , t > t 0 > 0 , H I t 0 1 α u ( t 0 + ) = u 1 , H D t 0 α u ( t 0 + ) = u 2 , u 1 , u 2 R ,
where H D t 0 α , H D t 0 α 1 and H D t 0 α 2 are the Hadamard fractional derivatives of orders α , α 1 and α 2 , respectively, 0 α 1 < α < 1 and 0 α 2 < α < 1 . The operator H I t 0 ρ is the Hadamard fractional integral of order ρ 0 . The definitions of these operators are given in Section 2.
We shall investigate the asymptotic behavior of solutions for Problem (1). Sufficient conditions on the nonlinear source term guaranteeing the convergence of solutions to logarithmic functions, for large values of time, are established. The importance of using analytical techniques to study the asymptotic behavior of solutions for Problem (1) arises from the lack of explicit solutions.
It is known that solutions for many kinds of (integer-order) ordinary differential equations may approach a certain function as time goes to infinity; in particular they may decay to zero, oscillate, or blow up in finite time. Many results in this regard exist in the literature. For example, we refer the reader to the papers [1,2,3,4,5,6], in which various classes of linear and nonlinear ordinary differential equations have been studied. Generalizing the existing results from integer orders to non-integer fractional orders is of great importance due to their numerous applications; see for instance, [7,8,9,10]. Unfortunately, imitating the techniques verbatim is not straightforward. Many difficulties arise when trying to do so. Some of these difficulties are due to the nature of the fractional derivatives themselves as they involve by definition all the past memory of solutions as well as nonregular kernels. In addition, many fundamental properties of integer-order derivatives are not valid for fractional-order derivatives. The chain rule is an example of such invalid properties. We will go around these difficulties by utilizing some adequate estimations, like desingularization methods, to deal with singular terms and by modifying and/or generalizing some versions of Bihari–LaSalle inequality.
The study of the asymptotic behavior of solutions for fractional differential equations, with Riemann–Liouville or Caputo fractional derivatives, has been investigated by many researchers, see e.g., [11,12,13,14,15,16,17,18,19,20]. The authors of [15] considered the fractional differential equation:
D 0 + α u t + f ( t , u ) = 0 , t > 0 ,
under the condition
f ( t , u ) ϕ t , u ( 1 + t ) α , t 0 , u R ,
where D 0 + α is the Riemann–Liouville fractional derivative of order α , 0 < α < 1 . The function f : [ 0 , ) × R R is continuous and the function ϕ : [ 0 , ) × [ 0 , ) [ 0 , ) is continuous in each argument and nondecreasing in the second one. They proved that the solutions can be represented asymptotically as a 1 + a 2 t α + O ( t α 1 ) , a 1 , a 2 R .
In [20], the case when the source function f depends on the solution and its sub-first-order fractional derivative has been considered, namely,
C D t 0 + α u t = f t , u t , C D t 0 + α 1 u t , t > t 0 > 0 , u t 0 = u 0 ,
where C D t 0 + α is the Caputo fractional derivative of order α , 0 < α 1 < α < 1 . The authors showed that any global solution of Problem (2) is asymptotic to c t α 1 for some real number c.
The present authors investigated the boundedness, power-type decay and asymptotic behavior of solutions for the initial value problems:
C D 0 + α u ( t ) = f t , u ( t ) , 0 t h t , s , u ( s ) d s , t > 0 , 0 < α 1 u ( 0 ) = u 0 , u 0 R ,
and
D 0 + α + 1 u ( t ) = f t , ( D 0 + α 1 u ) ( t ) , 0 t h t , s , D 0 + α 2 u ( s ) d s , t > 0 , 0 α 1 , α 2 α < 1 I 0 + 1 α u ( 0 + ) = u 0 , D 0 + α u ( 0 + ) = u 1 , u 0 , u 1 R .
Several different classes of source functions f such as
f ( t , u , υ ) k ( t ) P ( u ) + l ( t ) Q ( υ ) , f ( t , u , υ ) k ( t ) P ( t α 2 u ) Q ( υ ) , or f ( t , u , υ ) k ( t ) P t 1 α + α 1 u + l ( t ) Q ( υ ) ,
and on the kernel h such as:
h ( t , s , u ) w ( s ) K ( u ) , h ( t , s , u ) w ( s ) K ( s α 2 u ) , h ( t , s , u ) w ( t , s ) K ( u ) , or h ( t , s , u ) w ( s ) K t 1 α + α 2 u ,
for some functions k , l , P , Q , w and K have been treated, see [11,12,21,22,23].
For fractional differential equations with Hadamard-type fractional derivatives, we found relatively few results in the literature tackling the long-time behavior of solutions of fractional initial value problems, see, [24,25,26,27,28]. The authors of [28] studied the stability and decay rate of the zero solution of the fractional differential problem:
H D t 0 α u ( t ) = f u t , t > t 0 > 0 , 0 < α < 1 H I t 0 1 α u ( t 0 + ) = u 0 , u 0 R .
They considered first the linear case, f u = c u , c R , and established a criteria for the decay rate of solutions and the Lyapunov stability of the zero equilibrium. For the nonlinear case, they obtained the stability and decay rate of the hyperbolic zero equilibrium. They used a modified Laplace transform to express solutions by Mittag–Leffler functions and then used the asymptotic expansions of these functions to discuss the stability and logarithmic decay of the solutions. In [25], the authors discussed the stability of logarithmic type for the initial value problem:
H D t 0 α u ( t ) = f t , u t , H D t 0 α 1 u t , t > t 0 > 0 , 0 < α 1 < α < 1 , H D t 0 α 1 u ( t 0 + ) = u 0 , u 0 R .
Under some sufficient growth conditions of f, it has been shown that the solutions decay to zero as the logarithmic function ln t t 0 α 1 . Recently, the same authors considered Problem (4) in [26] with d d t H D t 0 α u ( t ) on the left hand side and the additional initial condition H D t 0 α u ( t 0 + ) = u 1 R . They showed that solutions approach a logarithmic function as time goes to infinity.
To the best of our knowledge, the long-time behavior of solutions for the class of fractional integrodifferential equations with Hadamard fractional derivatives (1) has not been investigated so far. In this paper, we prove that solutions of (1) are asymptotic to the logarithmic function ln t t 0 α where α is the order of the involved fractional derivative. Under sufficient growth conditions, we show that there exist a real number r such that any solution for (1) in the space u C 1 α , ln α + 1 t 0 , , see (16), has the following property lim t u t ln t t 0 α = r .
The rest of this paper is organized as follows. In the next section, Section 2, we give some notations from fractional calculus and present some preliminary results. In Section 3, we introduce and prove our main results. Section 4 is devoted to an example that supports our results. A brief conclusion is presented at the end of the study in Section 5.

2. Preliminaries

This section is devoted to briefly introduce some basic definitions, notions, and properties from fractional calculus and fractional differential equations theory which will be used in further considerations.
Definition 1 
([7]). We denote by C γ , ln t 0 , T , 0 γ < 1 , the following weighted space of continuous functions:
C γ , ln t 0 , T = ξ : t 0 , T R : ln t t 0 γ ξ ( t ) C t 0 , T ,
with the norm
ξ C γ , ln = ln t t 0 γ ξ ( t ) C 0 , ln ,
where C t 0 , T = C 0 , ln t 0 , T is the space of continuous functions on t 0 , T .
Definition 2 
([7]). Let δ = t d d t be the δ-derivative. For n N and 0 γ < 1 , the weighted space of continuously δ-differentiable functions up to order n 1 with nth δ-derivative in C γ , ln t 0 , T , is denoted by C δ , γ n t 0 , T and defined by:
C δ , γ n t 0 , T = ξ : t 0 , T R | δ i ξ C t 0 , T , i = 0 , 1 , 2 , , n 1 , δ n ξ C γ , ln t 0 , T ,
with the norm
ξ C δ , γ n = i = 0 n 1 δ i ξ C + δ n ξ C γ , ln .
In particular, C δ , γ 0 t 0 , T = C γ , ln t 0 , T .
A characterization of the space C δ , γ n t 0 , T is given as follows [7]: The functions ξ in the space C δ , γ n t 0 , T , n N and 0 γ < 1 , can be represented as:
ξ ( t ) = 1 ( n 1 ) ! t 0 t ln t s n 1 h ( s ) s d s + i = 0 n 1 b i ln t t 0 i ,
where h C γ , ln t 0 , T and b i , i = 0 , 1 , 2 , , n 1 , are arbitrary constants. In fact, h ( t ) = δ n ξ t and b i = δ i ξ t 0 i ! .
Definition 3 
([7]). The Hadamard left-sided fractional integral of order α > 0 is defined by:
H I t 0 α w ( t ) = 1 Γ ( α ) t 0 t ln t s α 1 w ( s ) s d s , t > t 0 ,
provided the right-hand side exists. We define H I t 0 0 w = w . The function Γ is the Euler gamma function defined by Γ α = 0 t α 1 e t d t , α > 0 , where t α 1 = e α 1 ln t .
Definition 4 
([7]). The Hadamard left-sided fractional derivative of order α > 0 , is defined by:
H D t 0 α w ( t ) = 1 Γ ( n α ) δ n t 0 t ln t s n α 1 w ( s ) s d s , t > t 0 ,
that is,
H D t 0 α w t = δ n H I t 0 n α w ( t ) , t > t 0 ,
where δ n = t d d t n , n = α . In particular, when α = m N 0 , we have H D t 0 m w = δ m w .
The next lemma shows that the Hadamard fractional derivative (or integral) of a logarithmic function results in a multiple of the same logarithmic function with the order of the fractional derivative (or integral) subtracted from (or added to) its power.
Lemma 1 
([7]). If α > 0 , β > 0 , then:
H I t 0 α ln s t 0 β 1 t = Γ β Γ ( β + α ) ln t t 0 β + α 1 , t > t 0 ,
H D t 0 α ln s t 0 β 1 t = Γ β Γ ( β α ) ln t t 0 β α 1 , t > t 0 .
In particular, when 0 < α < 1 , β = 1 , then:
H I t 0 α 1 t = 1 Γ ( 1 + α ) ln t t 0 α , t > t 0 ,
H D t 0 α 1 t = 1 Γ ( 1 α ) ln t t 0 α , t > t 0 .
The last property shows that the Hadamard-type derivative of a constant is not zero.
The composite of the Hadamard operators of fractional differentiation and integration with different orders is given next.
Lemma 2 
([7]). Let 0 < β < α and 0 γ < 1 , then:
H D t 0 H β I t 0 α w = H I t 0 α β w
at every point in t 0 , T if w C γ , ln t 0 , T and at every point in t 0 , T if w C t 0 , T . In particular, when α > β = m N , we obtain:
H D t 0 H m I t 0 α w = H I t 0 α m w .
The Hadamard differentiation operator is the left inverse to the associated Hadamard integration operator [7]. That is, H D t 0 H α I t 0 α w = w at every point in t 0 , T if w C γ , ln t 0 , T . This property is not valid when the Hadamard fractional derivative and the Hadamard fractional integral are inverted as shown in the lemma below.
Lemma 3 
([7]). Let α 0 and n = α . If w C γ , ln t 0 , T , 0 γ < 1 and H I t 0 n α w C δ , γ n t 0 , T , then:
H I t 0 H α D t 0 α w t = w t i = 1 n δ n i H I t 0 n α w ( t 0 + ) Γ ( α i + 1 ) ln t t 0 α i ,
for all t t 0 , T . In particular, for 0 α < 1 , we have:
H I t 0 H α D t 0 α w t = w t H I t 0 1 α w t 0 + Γ ( α ) ln t t 0 α 1
at every point in t 0 , T if w C γ , ln t 0 , T and at every point in t 0 , T if w C t 0 , T and H I t 0 1 α w C δ , γ 1 t 0 , T .
For more about Hadamard fractional integral and derivative, we refer to the books [7,29,30].
The limit of the ratio of the Hadamard fractional integral H I t 0 α + 1 s f s , u ( s ) , v s t and the power function 1 Γ ( α + 1 ) ln t t 0 α as t is treated in the lemma below.
Lemma 4. 
Let f L 1 t 0 , , t 0 > 0 . Suppose that u and v are real-valued functions defined on t 0 , , then:
lim t H I t 0 α + 1 s f s , u ( s ) , v s t 1 Γ ( α + 1 ) ln t t 0 α = lim t 1 ln t t 0 α t 0 t ln t s α f s , u ( s ) , v s d s = t 0 f s , u ( s ) , v s d s , α > 0 .
Proof. 
It is easy to see that:
H I t 0 α + 1 s f s , u ( s ) , v s t 1 Γ ( α + 1 ) ln t t 0 α t 0 f s , u ( s ) , v s d s = 1 ln t t 0 α t 0 t ln t s α f s , u ( s ) , v s d s t 0 f s , u ( s ) , v s d s = t 0 t ln t t 0 ln s t 0 ln t t 0 α f s , u ( s ) , v s d s t 0 f s , u ( s ) , v s d s = t 0 1 ln s t 0 ln t t 0 α χ t 0 , t s 1 f s , u ( s ) , v s d s t 0 1 ln s t 0 ln t t 0 α χ t 0 , t s 1 f s , u ( s ) , v s d s ,
where
χ t 0 , t s = 1 if t 0 s t , 0 , otherwise . .
As f L 1 t 0 , and
lim t 1 ln s t 0 ln t t 0 α χ t 0 , t s = 1 , for s < t ,
we get from the Dominated Convergence Theorem [31],
lim t H I t 0 α + 1 s f s , u ( s ) , v s t 1 Γ ( α + 1 ) ln t t 0 α t 0 f s , u ( s ) , v s d s lim t t 0 1 ln s t 0 ln t t 0 α χ t 0 , t s 1 f s , u ( s ) , v s d s = t 0 lim t 1 ln s t 0 ln t t 0 α χ t 0 , t s 1 f s , u ( s ) , v s d s = 0 .
This completes the proof. □
The next lemma can be considered as a Hadamard fractional version of L’Hôpital’s rule when applied to the solution of problem (1).
Lemma 5. 
Let u be a solution of problem (1) with f L 1 0 , . Then,
lim t u ( t ) ln t t 0 α = lim t H D t 0 α u ( t ) Γ ( α + 1 ) = 1 Γ ( α + 1 ) u 2 + t 0 f s , H D t 0 α 1 u ( s ) , t 0 s h ( s , τ , H D t 0 α 2 u ( τ ) ) d τ d s .
Proof. 
Integrating both sides of the equation in (1) over the interval t 0 , t , we obtain:
H D t 0 α u ( t ) = u 2 + t 0 t f s , H D t 0 α 1 u ( s ) , t 0 s h ( s , τ , H D t 0 α 2 u ( τ ) ) d τ d s .
Apply H I t 0 α to (11) and use Lemmas 1 and 3, to find that:
u ( t ) u 1 ln t t 0 α 1 Γ ( α ) = u 2 ln t t 0 α Γ ( α + 1 ) + 1 Γ ( α ) t 0 t ln t s α 1 × t 0 s f τ , H D t 0 α 1 u ( τ ) , t 0 τ h ( τ , σ , H D t 0 α 2 u ( σ ) ) d σ d τ d s s ,
for all t > t 0 > 0 . Reorder the double integral on the right-hand side and integrate by substitution to have:
u ( t ) = u 1 ln t t 0 α 1 Γ ( α ) + u 2 ln t t 0 α Γ ( α + 1 ) + 1 Γ ( α ) t 0 t τ t ln t s α 1 × f τ , H D t 0 α 1 u ( τ ) , t 0 τ h ( τ , σ , H D t 0 α 2 u ( σ ) ) d σ d s s d τ = u 1 ln t t 0 α 1 Γ ( α ) + u 2 ln t t 0 α Γ ( α + 1 ) + 1 Γ ( α + 1 ) t 0 t ln t τ α × f τ , H D t 0 α 1 u ( τ ) , t 0 τ h ( τ , σ , H D t 0 α 2 u ( σ ) ) d σ d τ .
Dividing both sides of (13) by ln t t 0 α gives:
u ( t ) ln t t 0 α = u 1 Γ ( α ) ln t t 0 + u 2 Γ ( α + 1 ) + 1 Γ ( α + 1 ) ln t t 0 α t 0 t ln t τ α × f τ , H D t 0 α 1 u ( τ ) , t 0 τ h ( τ , σ , H D t 0 α 2 u ( σ ) ) d σ d τ ,
for all t > t 0 > 0 . Taking the limit of the resulting ratio at infinity and applying Lemma 4 leads to the desired result. □
In the next lemma, we recall Bihari–LaSalle inequality which is a nonlinear generalization of the well-known Grönwall–Bellman inequality.
Lemma 6 
([32,33]). Suppose that w and g are nonnegative continuous functions on t 0 , and φ is a positive function on 0 , , continuous and nondecreasing on 0 , . If
w ( t ) c + t 0 t g ( τ ) φ ( w ( τ ) ) d τ , t t 0 , ,
where c > 0 , then
w t Φ 1 Φ ( c ) + t 0 t g ( τ ) d τ , t t 0 , T 1 ,
where Φ 1 is the inverse function of Φ ,
Φ ( x ) = x 0 x d s φ ( s ) , x > 0 , x 0 > 0 ,
and T 1 is chosen so that Φ ( c ) + t 0 t g ( τ ) d τ is in the domain of Φ 1 for all t t 0 , T 1 .
The next nonlinear integral inequality can be considered a generalization of Bihari–LaSalle inequality that has been recalled in Lemma 6.
Lemma 7 
([34]). Assume that w and η j , j = 1 , , n are nonnegative continuous functions on [ t 0 , T ] and φ j , j = 1 , , n are nonnegative, continuous and nondecreasing on 0 , such that φ 1 φ 2 φ n (that is, φ n / φ n 1 , , φ 2 / φ 1 are nondecreasing functions). Assume further that c is a positive constant and
w ( t ) c + j = 1 n t 0 t η j ( s ) φ j ( w ( s ) ) d s , t [ t 0 , T ] ,
then, for all t t 0 , T 1 ,
w t Φ n 1 Φ n ( c n 1 ) + t 0 t η n ( s ) d s ,
where
1
Φ j 1 is the inverse function of Φ j and Φ j ( x ) = x j x d s φ j ( s ) , x > 0 ,   x j > 0 , j = 1 , , n .
2
The constants c j are given by c 0 = c and c j = Φ j 1 Φ j ( c j 1 ) + t 0 T 1 η j ( τ ) d τ , j = 1 , , n 1 .
3
The number T 1 t 0 , T is the largest number such that:
t 0 T 1 η j ( τ ) d τ c j 1 d τ φ j ( τ ) , j = 1 , , n .
Lemma 8 
([34]). Suppose that w and η j , j = 1 , 2 , 3 are nonnegative continuous functions on [ t 0 , T ] and φ j , j = 1 , 2 , 3 are nonnegative, continuous and nondecreasing on 0 , such that φ 1 φ 2 φ 3 , (that is, φ 3 / φ 2 and φ 2 / φ 1 are nondecreasing functions). Assume further that c is a positive constant. If
w ( t ) c + t 0 t η 1 ( τ ) φ 1 ( w ( τ ) ) d τ + t 0 t η 2 ( τ ) φ 2 t 0 τ η 3 ( s ) φ 3 ( w ( s ) ) d s d τ ,
then, for all t t 0 , T 1 ,
w ( t ) Φ 3 1 Φ 3 ( c 2 ) + t 0 t η 3 ( τ ) d τ ,
where the functions Φ j 1 , Φ j , j = 1 , 2 , 3 and the constants c 0 , c 1 , c 2 are the same as those given in Lemma 7.
Remark 1. 
By considering the following functions,
ϖ 1 ( t ) : = max τ 0 , t φ 1 τ , ϖ j ( t ) : = max τ 0 , t φ j τ ϖ j 1 ( τ ) ϖ j 1 ( t ) , j = 2 , 3 ,
we can drop the ordering and monotonicity requirements in Lemma 8. It is clear that φ j t ϖ j t , ϖ j are nonnegative and nondecreasing functions on 0 , for all t 0 , , j = 1 , 2 , 3 and ϖ 1 ϖ 2 ϖ 3 .

3. Main Results

In this section, we study the asymptotic behavior of continuable solutions for the problem (1) in the space C 1 α , ln α + 1 t 0 , T , 0 < t 0 < T , defined by:
C 1 α , ln α + 1 t 0 , T = u : 0 , T R | u C 1 α , ln t 0 , T , d H D t 0 α u t d t C 1 α , ln t 0 , T ,
where the space C 1 α , ln t 0 , T is defined in (5).
The following two types of functions are used repeatedly in the rest of this paper.
Definition 5. 
A function g is said to be Ω ϰ -type function if it is continuous, nonnegative on t 0 , and
t 0 ln t t 0 ϰ g t d t < , ϰ = 0 or 1 .
Definition 6. 
A function ψ is said to be Ψ-type function if it is continuous, positive, nondecreasing on 0 , ,
ψ a b ψ a b for all 1 , b 1 ,
and
ρ t d s ψ ( s ) as t for any ρ > 0 .
The two classes of functions introduced above are not empty, see the example provided in Section 4.
Consider the following nonlinear inequality:
w t a 1 + ln t t 0 a 2 + a 3 t 0 t g 1 ( s ) φ 1 w s + g 2 ( s ) φ 2 t 0 s g 3 ( τ ) φ 3 w τ d τ d s ,
for all t > t 0 , where a 1 , a 2 and a 3 are positive constants, g 1 , g 3 are Ω 1 -type functions, g 2 is Ω 0 -type function and φ j ,   j = 1 , 2 , 3 are Ψ -type functions with φ 1 φ 2 φ 3 .
Let
K : = Φ 3 b 3 + t 0 t 0 e g 3 ( s ) d s ,
and
K 1 : = Φ 3 ( c 2 ) + t 0 e ln t t 0 g 3 ( t ) d t ,
where
b 3 = Φ 2 1 Φ 2 ( b 2 ) + a 3 t 0 t 0 e g 2 ( s ) d s , b 2 = Φ 1 1 Φ 1 a 1 + a 2 + a 3 t 0 t 0 e g 1 ( s ) d s , c 2 = Φ 2 1 Φ 2 ( c 1 ) + a 3 t 0 e g 2 ( t ) d t , c 1 = Φ 1 1 Φ 1 ( K 2 ) + a 3 t 0 e ln t t 0 g 1 ( t ) d t , K 2 = a 1 + a 2 + a 3 φ 1 Φ 3 1 K t 0 t 0 e g 1 ( s ) d s + a 3 φ 2 φ 3 Φ 3 1 K t 0 t 0 e g 3 ( τ ) d τ t 0 t 0 e g 2 ( s ) d s ,
and Φ j 1 is the inverse functions of Φ j x = ρ x d s φ j ( s ) , j = 1 , 2 , 3 , x > ρ > 0 .
A generalized version of Lemma 8 is introduced below. We shall provide an estimate for an integral term that arises later in our present problem. Although this estimate is not the best possible, it ensures that such an integral is bounded, which is exactly what we will need to prove our results below.
Theorem 1. 
Suppose that w ( t ) is a continuous nonnegative function on t 0 , satisfying Inequality (17) for all t > t 0 , a 1 , a 2 and a 3 are positive constants, g 1 , g 3 are Ω 1 -type functions, g 2 is Ω 0 -type function and φ j ,   j = 1 , 2 , 3 , are Ψ-type functions with φ 1 φ 2 φ 3 . Then,
w ( t ) Φ 3 1 K for all t 0 < t < t 0 e ,
and
w ( t ) Φ 3 1 K 1 ln t t 0 for all t t 0 e ,
where K and K 1 are given in (18) and (19), respectively.
Proof. 
For t 0 < t < t 0 e , the inequality (17) becomes:
w t a 1 + a 2 + a 3 t 0 t g 1 ( s ) φ 1 w s + g 2 ( s ) φ 2 t 0 s g 3 ( τ ) φ 3 w τ d τ d s , t 0 < t < t 0 e ,
and Lemma 8 is applicable. The relation (21) follows immediately.
For the case t t 0 e , we have ln t t 0 1 and (17) may be rewritten as:
w ( t ) ln t t 0 a 1 + a 2 + a 3 t 0 t g 1 ( s ) φ 1 w s + g 2 ( s ) φ 2 t 0 s g 3 ( τ ) φ 3 w τ d τ d s a 1 + a 2 + a 3 t 0 t 0 e g 1 ( s ) φ 1 w s + g 2 ( s ) φ 2 t 0 s g 3 ( τ ) φ 3 w τ d τ d s + a 3 t 0 e t g 1 ( s ) φ 1 w s + g 2 ( s ) φ 2 t 0 s g 3 ( τ ) φ 3 w τ d τ d s .
In virtue of the estimate (21) together with the continuity and monotonicity of the functions φ i , i = 1 , 2 , 3 , we get:
w ( t ) ln t t 0 K 2 + a 3 t 0 e t g 1 ( s ) φ 1 w s + g 2 ( s ) φ 2 t 0 s g 3 ( τ ) φ 3 w τ d τ d s ,
where K 2 is the constant given in (20).
Let
z : = z 1 + z 2 + z 3 ,
where
z 1 t : = K 2 + a 3 t 0 e t g 1 ( s ) φ 1 w s d s , t t 0 e , z 2 t : = a 3 t 0 e t g 2 ( s ) φ 2 z 3 s d s , t t 0 e , z 3 t : = t 0 t g 3 ( s ) φ 3 w s d s , t > t 0 .
It is clear that
w ( t ) ln t t 0 < z t for all t t 0 e .
In light of the types of the functions g i , and φ i , i = 1 , 2 , 3 , see Definitions 5 and 6, differentiating z yields:
z t = a 3 g 1 ( t ) φ 1 w t + a 3 g 2 ( t ) φ 2 z 3 t + g 3 ( t ) φ 3 w t a 3 g 1 ( t ) ln t t 0 φ 1 w ( t ) ln t t 0 + a 3 g 2 ( t ) φ 2 z t + g 3 ( t ) ln t t 0 φ 3 w ( t ) ln t t 0 a 3 ln t t 0 g 1 ( t ) φ 1 z t + a 3 g 2 ( t ) φ 2 z t + ln t t 0 g 3 ( t ) φ 3 z t , t t 0 e .
Now, we integrate both sides of (29) over the interval t 0 e , t to find:
z t z t 0 e + a 3 t 0 e t ln s t 0 g 1 ( s ) φ 1 ( z s ) d s + a 3 t 0 e t g 2 ( s ) φ 2 ( z s ) d s + t 0 e t ln s t 0 g 3 ( s ) φ 3 ( z s ) d s .
Applying Lemma 7 with η 1 ( t ) = a 3 ln t t 0 g 1 ( t ) , η 2 ( t ) = a 3 g 2 ( t ) , η 3 ( t ) = ln t t 0 g 3 ( t ) and T 1 = (by assumption ρ d τ φ i ( τ ) = , i = 1 , 2 , 3 for any ρ > 0 ), the inequality (30) leads to:
z ( t ) Φ 3 1 K 1 , for all t t 0 e ,
as desired. □
From now on, we assume that the functions f and h on the right hand-side of (1) , satisfy the conditions below:
(A1)
f ( t , u , υ ) : t 0 , × R 2 R is a C 1 α , ln t 0 , function in E 1 = { ( t , u , υ ) : t > t 0 > 0 , u , υ C 1 α , ln t 0 , } .
(A2)
h ( t , s , u ) is a continuous function in E 2 = { ( t , s , u ) : t 0 s < t < , u C 1 α , ln t 0 , } .
(A3) 
There are Ω 1 -type functions g 1 , g 3 , an Ω 0 -type function g 2 and Ψ -type functions φ i , i = 1 , 2 , 3 with φ 1 φ 2 φ 3 such that:
f ( t , u , υ ) g 1 ( t ) φ 1 ln t t 0 1 α + α 1 u + g 2 ( t ) φ 2 υ , ( t , u , υ ) E 1 ,
h ( t , s , u ) g 3 ( s ) φ 3 ln s t 0 1 α + α 2 u , ( t , s , u ) E 2 ,
0 α 1 < α < 1 a n d 0 α 2 < α < 1 .
Functions satisfying the above hypotheses are given in Section 4.
In the next result, we shall show that the solution for Problem (1) satisfies the useful nonlinear inequality below.
Lemma 9. 
Suppose that u ( t ) is a C 1 α , ln α + 1 t 0 , -solution for Problem (1) and the functions f and h satisfy (A 1 ), (A 2 ) and (A 3 ). Then,
ln t t 0 1 α u ( t ) v ( t ) , t > t 0 ,
and
ln t t 0 1 α + α j H D t 0 α j u ( t ) v j ( t ) , t > t 0 ,
where
v ( t ) = u 1 Γ ( α ) + ln t t 0 Γ ( α + 1 ) u 2 + t 0 t g 1 ( s ) φ 1 ln s t 0 1 α + α 1 H D t 0 α 1 u ( s ) + g 2 ( s ) φ 2 t 0 τ g 3 ( τ ) φ 3 ln τ t 0 1 α + α 2 H D t 0 α 2 u τ d τ d s ,
and
v j ( t ) = u 1 Γ α α j + ln t t 0 Γ α α j + 1 u 2 + t 0 t g 1 ( s ) φ 1 ln s t 0 1 α + α 1 H D t 0 α 1 u ( s ) + g 2 ( s ) φ 2 t 0 τ g 3 ( τ ) φ 3 ln τ t 0 1 α + α 2 H D t 0 α 2 u τ d τ d s ,
for all t > t 0 and j = 1 , 2 .
Proof. 
Firstly, we recall the relations (12) and (13) from the proof of Lemma 5,
u ( t ) = u 1 ln t t 0 α 1 Γ ( α ) + u 2 ln t t 0 α Γ ( α + 1 ) + 1 Γ ( α ) t 0 t ln t s α 1 × t 0 s f τ , H D t 0 α 1 u ( τ ) , t 0 τ h ( τ , σ , H D t 0 α 2 u ( σ ) ) d σ d τ d s s , t > t 0 ,
u ( t ) = u 1 ln t t 0 α 1 Γ ( α ) + u 2 ln t t 0 α Γ ( α + 1 ) + 1 Γ ( α + 1 ) t 0 t ln t τ α × f τ , H D t 0 α 1 u ( τ ) , t 0 τ h ( τ , σ , H D t 0 α 2 u ( σ ) ) d σ d τ , t > t 0 .
Equation (36) leads to the following bound on u ( t ) ,
u ( t ) u 1 ln t t 0 α 1 Γ ( α ) + u 2 ln t t 0 α Γ ( α + 1 ) + ln t t 0 α Γ ( α + 1 ) × t 0 t f τ , H D t 0 α 1 u ( τ ) , t 0 τ h ( τ , σ , H D t 0 α 2 u ( σ ) ) d σ d τ , t > t 0 ,
or
ln t t 0 1 α u ( t ) u 1 Γ ( α ) + u 2 ln t t 0 Γ ( α + 1 ) + ln t t 0 Γ ( α + 1 ) × t 0 t f τ , H D t 0 α 1 u ( τ ) , t 0 τ h ( τ , σ , H D t 0 α 2 u ( σ ) ) d σ d τ ,
for all t > t 0 , which, in the light of the hypothesis (A 3 ), leads to (31) as desired.
Now, apply the Hadamard differentiation operator H D t 0 α j , j = 1 , 2 to (35), then employ Lemmas 1 and 2 (with β = α j and γ = 1 α , see (16)) to get:
H D t 0 α j u ( t ) = u 1 ln t t 0 α α j 1 Γ α α j + u 2 ln t t 0 α α j Γ α α j + 1 + 1 Γ α α j t 0 t ln t s α α j 1 × t 0 s f τ , H D t 0 α 1 u ( τ ) , t 0 τ h ( τ , σ , H D t 0 α 2 u ( σ ) ) d σ d τ d s s , t > t 0 .
In a similar way to that used to obtain (36), the relation (39) reduces to:
H D t 0 α j u ( t ) = u 1 ln t t 0 α α j 1 Γ α α j + u 2 ln t t 0 α α j Γ α α j + 1 + 1 Γ α α j + 1 t 0 t ln t τ α α j × f τ , H D t 0 α 1 u ( τ ) , t 0 τ h ( τ , σ , H D t 0 α 2 u ( σ ) ) d σ d τ , t > t 0 ,
and consequently,
H D t 0 α j u ( t ) u 1 ln t t 0 α α j 1 Γ α α j + u 2 ln t t 0 α α j Γ α α j + 1 + 1 Γ α α j + 1 t 0 t ln t τ α α j × f τ , H D t 0 α 1 u ( τ ) , t 0 τ h ( τ , σ , H D t 0 α 2 u ( σ ) ) d σ d τ u 1 ln t t 0 α α j 1 Γ α α j + ln t t 0 α α j Γ α α j + 1 u 2 + t 0 t f τ , H D t 0 α 1 u ( τ ) , t 0 τ h ( τ , σ , H D t 0 α 2 u ( σ ) ) d σ d τ , t > t 0 ,
which yields, in light of the hypothesis (A 3 ), the desired result in (32) and completes the proof. □
The main result of this section is given below.
Theorem 2. 
Suppose that the hypotheses (A 1 ), (A 2 ) and (A 3 ) are satisfied, then there exists a real number r such that any solution u C 1 α , ln α + 1 t 0 , of Problem (1), has the following asymptotic property lim t u t ln t t 0 α = r .
Proof. 
Let us start by recalling the relations (31) and (32) in Lemma 9,
ln t t 0 1 α u ( t ) v ( t ) , t > t 0 ,
and
ln t t 0 1 α + α j H D t 0 α j u ( t ) v j ( t ) , t > t 0 , j = 1 , 2 ,
where u is a C 1 α , ln α + 1 t 0 , -solution for Problem (1), v ( t ) and v j ( t ) , j = 1 , 2 are given in (33) and (34), respectively.
Let
a 1 = u 1 max j = 1 , 2 1 Γ ( α ) , 1 Γ α α j , a 2 = u 2 max j = 1 , 2 1 Γ α + 1 , 1 Γ α α j + 1 , a 3 = max j = 1 , 2 1 Γ α + 1 , 1 Γ α α j + 1 ,
and
w ( t ) = a 1 + ln t t 0 a 2 + a 3 t 0 t g 1 ( s ) φ 1 ln s t 0 1 α + α 1 H D t 0 α 1 u ( s ) + g 2 ( s ) φ 2 t 0 τ g 3 ( τ ) φ 3 ln τ t 0 1 α + α 2 H D t 0 α 2 u τ d τ d s , t > t 0 ,
then,
v ( t ) w t and v j ( t ) w t for all t > t 0 , j = 1 , 2 .
From the nondecreasingness property of the functions φ i , i = 1 , 2 , 3 , we have:
w ( t ) = a 1 + ln t t 0 a 2 + a 3 t 0 t g 1 ( s ) φ 1 v 1 ( s ) + g 2 ( s ) φ 2 t 0 τ g 3 ( τ ) φ 3 v 2 ( τ ) d τ d s a 1 + ln t t 0 a 2 + a 3 t 0 t g 1 ( s ) φ 1 w s + g 2 ( s ) φ 2 t 0 τ g 3 ( τ ) φ 3 w τ d τ d s ,
for all t > t 0 , which is Inequality (17). Therefore, we get from Theorem 1 that:
w ( t ) ln t t 0 Φ 3 1 K 1 for all t t 0 e ,
and as a consequence of the estimates (42) and (46), we find:
ln t t 0 1 α u ( t ) ln t t 0 Φ 3 1 K 1 for all t > t 0 e ,
or
u ( t ) ln t t 0 α K 3 : = Φ 3 1 K 1 for all t > t 0 e ,
where K 1 is as in (19).
Now, in the light of Hypothesis (A 3 ), and the inequalities (43) and (46), we deduce that:
J : = t 0 t f s , H D t 0 α 1 u ( s ) , t 0 s h ( s , τ , H D t 0 α 2 u ( τ ) ) d τ d s t 0 t f s , H D t 0 α 1 u ( s ) , t 0 s h ( s , τ , H D t 0 α 2 u ( τ ) ) d τ d s t 0 t g 1 ( s ) φ 1 ln s t 0 1 α + α 1 H D t 0 α 1 u ( s ) d s + t 0 t g 2 ( s ) × φ 2 t 0 s g 3 ( τ ) φ 3 ln τ t 0 1 α + α 2 H D t 0 α 2 u ( τ ) d τ d s t 0 t 0 e g 1 ( s ) φ 1 w ( s ) d s + t 0 t 0 e g 2 ( s ) φ 2 t 0 s g 3 ( τ ) φ 3 w ( τ ) d τ d s + t 0 e t g 1 ( s ) φ 1 w ( s ) d s + t 0 e t g 2 ( s ) φ 2 t 0 s g 3 ( τ ) φ 3 w ( τ ) d τ d s = J 1 + J 2 .
The first integral,
J 1 : = t 0 t 0 e g 1 ( s ) φ 1 w ( s ) d s + t 0 t 0 e g 2 ( s ) φ 2 t 0 s g 3 ( τ ) φ 3 w ( τ ) d τ d s t 0 t 0 e g 1 ( s ) φ 1 Φ 3 1 K d s + t 0 t 0 e g 2 ( s ) φ 2 t 0 s g 3 ( τ ) φ 3 Φ 3 1 K d τ d s φ 1 K 4 t 0 t 0 e g 1 ( s ) d s + t 0 t 0 e g 2 ( s ) φ 2 φ 3 K 4 t 0 s g 3 ( τ ) d τ d s
is finite by (21) from Theorem 1, K 4 : = Φ 3 1 K and K is given in (18).
The second integral,
J 2 : = t 0 e t g 1 ( s ) φ 1 w ( s ) d s + t 0 e t g 2 ( s ) φ 2 t 0 s g 3 ( τ ) φ 3 w ( τ ) d τ d s , t t 0 e ,
can be estimated in view of (47), Theorem 1 and the type of the functions φ i , i = 1 , 2 , 3 as follows:
J 2 t 0 e t ln s t 0 g 1 ( s ) φ 1 K 3 d s + t 0 e t g 2 ( s ) φ 2 t 0 t 0 e g 3 ( τ ) φ 3 K 4 d τ + t 0 e s ln τ t 0 g 3 ( τ ) φ 3 K 3 d τ d s φ 1 K 3 t 0 e t ln s t 0 g 1 ( s ) d s + t 0 e t g 2 ( s ) φ 2 φ 3 K 4 t 0 t 0 e g 3 ( τ ) d τ + φ 3 K 3 t 0 e s ln τ t 0 g 3 ( τ ) d τ d s , for all t t 0 e .
Since g 1 and g 3 are Ω 1 -type functions and g 2 is an Ω 0 -type function, we see that the integral J 2 is uniformly bounded and so is the integral J .
Therefore, the integral t 0 t f s , H D t 0 α 1 u ( s ) , t 0 s h ( s , τ , H D t 0 α 2 u ( τ ) ) d τ d s absolutely convergent and, consequently,
lim t t 0 t f s , H D t 0 α 1 u ( s ) , t 0 s h ( s , τ , H D t 0 α 2 u ( τ ) ) d τ d s < ,
Using Lemma 5, there exits a finite real number r such that:
lim t u ( t ) ln t t 0 α = 1 Γ ( α + 1 ) u 2 + lim t t 0 t f s , H D t 0 α 1 u ( s ) , t 0 s h ( s , τ , H D t 0 α 2 u ( τ ) ) d τ d s = r ,
as wanted. □
Remark 2. 
If the condition (A 3 ) is replaced by the condition:
(A4)
There are Ω 1 -type functions ξ 1 ,   ξ 2 and Ψ -type functions ϕ i , i = 1 , 2 , 3 with ϕ 1 ϕ 2     ϕ 3 such that
f ( t , u , υ ) ξ 1 ( t ) ϕ 1 ln t t 0 1 α + α 1 u ϕ 2 υ , ( t , u , υ ) E 1 , 0 α 1 < α < 1 ,
and
h ( t , s , u ) ξ 2 ( s ) ϕ 3 ln s t 0 1 α + α 2 u , ( t , s , u ) E 2 , 0 α 2 < α < 1 ,
then the conclusion of Theorem 2 is still valid. We state this fact below.
Theorem 3. 
Suppose that the functions f and h satisfy the conditions (A 1 ), (A 2 ) and (A 4 ). Then, there exists a real number r such that any solution u C 1 α , ln α + 1 t 0 , for Problem (1) satisfies lim t u t ln t t 0 α = r .
The proof is skipped as it can be shown in a similar manner to that used in the proof of Theorem 2.
Remark 3. 
The problem,
H D t 0 α u ( t ) = f t , u ( t ) , t 0 t h t , s , u ( s ) d s , t > t 0 > 0 , 0 < α < 1 , H I t 0 1 α u ( t 0 + ) = u 1 , H D t 0 α u ( t 0 + ) = u 2 , u 1 , u 2 R .
is a special case of Problem (1) when α 1 = α 2 = 0 . Therefore, we can conclude that there exists a constant r R such that any solution u C 1 α , ln α + 1 t 0 , for Problem (56) satisfies lim t u t ln t t 0 α = r .

4. Example

Consider the initial value problem:
H D t 0 α u ( t ) = t σ 1 ln t t 0 β 1 H D t 0 α 1 u ( t ) γ 1 + t σ 2 ln t t 0 β 2 t 0 t t + s σ 3 × ln s t 0 β 3 H D t 0 α 2 u ( s ) γ 3 d s γ 2 , t > t 0 > 0 , H I t 0 1 α u ( t 0 + ) = u 1 , H D t 0 α u ( t 0 + ) = u 2 , u 1 , u 2 R ,
where 0 α 1 < α < 1 , 0 α 2 < α < 1 , σ j > 1 , j = 1 , 2 , 3 , β 1 > 1 α + α 1 γ 1 1 , β 2 > 1 , β 3 > 1 α + α 2 γ 3 1 and 0 < γ 1 γ 2 γ 3 < 1 .
The source function,
f t , u , v = t σ 1 ln t t 0 β 1 u γ 1 + t σ 2 ln t t 0 β 2 v γ 2 , t > t 0 > 0 ,
satisfies
f t , H D t 0 α 1 u ( t ) , t 0 t h t , s , H D t 0 α 2 u ( s ) d s = t σ 1 ln t t 0 β 1 H D t 0 α 1 u ( t ) γ 1 + t σ 2 ln t t 0 β 2 t 0 t t + s σ 3 × ln s t 0 β 3 H D t 0 α 2 u ( s ) γ 3 d s γ 2 g 1 ( t ) φ 1 ln t t 0 1 α + α 1 H D t 0 α 1 u ( t ) + g 2 ( t ) φ 2 t 0 t g 3 ( s ) × φ 3 ln s t 0 1 α + α 2 H D t 0 α 2 u ( s ) d s ,
with
h t , s , u = t + s σ 3 ln s t 0 β 3 u γ 3 , g 1 ( t ) = t σ 1 ln t t 0 β 1 + α α 1 1 γ 1 , g 2 ( t ) = t σ 2 ln t t 0 β 2 , g 3 ( t ) = t σ 3 ln t t 0 β 3 + α α 2 1 γ 3 , φ j t = t γ j , t > t 0 > 0 ,
where β 1 > 1 α + α 1 γ 1 1 , β 2 > 1 , β 3 > 1 α + α 2 γ 3 1 , σ j > 1 , j = 1 , 2 , 3 and 0 < γ 1 γ 2 γ 3 1 .
The functions g j , j = 1 , 2 , 3 are continuous, nonnegative on t 0 , ,
t 0 t 0 e g 1 ( t ) d t = t 0 t 0 e ln t t 0 β 1 + α α 1 1 γ 1 d t t σ 1 1 t 0 σ 1 1 t 0 t 0 e ln t t 0 β 1 + α α 1 1 γ 1 d t t = 1 t 0 σ 1 1 0 1 s β 1 + α α 1 1 γ 1 d s = 1 β 1 + α α 1 1 γ 1 + 1 t 0 σ 1 1 < ,
t 0 t 0 e g 3 ( t ) d t = t 0 t 0 e ln t t 0 β 3 + α α 2 1 γ 1 d t t σ 3 1 β 3 + α α 3 1 γ 3 + 1 t 0 σ 3 1 < ,
and
t 0 t 0 e g 2 ( t ) d t = t 0 t 0 e ln t t 0 β 2 d t t σ 2 1 t 0 σ 2 1 t 0 t 0 e ln t t 0 β 2 d t t = 1 β 2 + 1 t 0 σ 2 1 < .
The function g 2 is an Ω 0 -type because:
t 0 e g 2 ( t ) d t = t 0 e ln t t 0 β 2 d t t σ 2 = 1 t 0 σ 2 1 t 0 e ln t t 0 β 2 t 0 t σ 2 1 d t t 1 t 0 σ 2 1 t 0 ln t t 0 β 2 e σ 2 1 ln t t 0 d t t = 1 t 0 σ 2 1 0 s β 2 e σ 2 1 s d s = 1 σ 2 1 β 2 + 1 t 0 σ 2 1 0 τ β 2 e τ d τ = Γ β 2 + 1 σ 2 1 β 2 + 1 t 0 σ 2 1 < ,
whereas the functions g 1 and g 3 are Ω 1 -type functions since
t 0 e ln t t 0 g 1 ( t ) d t = t 0 e ln t t 0 β 1 + α α 1 1 γ 1 + 1 d t t σ 1 1 t 0 σ 1 1 t 0 ln t t 0 β 1 + α α 1 1 γ 1 + 1 e σ 1 1 ln t t 0 d t t = 1 t 0 σ 1 1 0 s β 1 + α α 1 1 γ 1 + 1 e σ 1 1 s d s = Γ β 1 + α α 1 1 γ 1 + 2 σ 1 1 β 1 + α α 1 1 γ 1 + 2 t 0 σ 1 1 < ,
and
t 0 e ln t t 0 g 3 ( t ) d t = t 0 e ln t t 0 β 3 + α α 2 1 γ 3 + 1 d t t σ 3 Γ β 3 + α α 2 1 γ 3 + 2 σ 3 1 β 3 + α α 2 1 γ 3 + 2 t 0 σ 3 1 < .
The functions φ j t = t γ j , j = 1 , 2 , 3 , 0 < γ 1 γ 2 γ 3 1 , are Ψ -type functions with φ 1     φ 2     φ 3 . These three functions are continuous, positive, nondecreasing on 0 , ,
φ j a = a γ j b a b γ j = b φ j a b for all a 1 , b 1 , 0 < γ j 1 ,
and
ρ x d t φ j t = ρ x d t t γ j as x for any ρ > 0 .
As 0 < γ 1 γ 2 γ 3 1 , then φ 3 / φ 2 and φ 2 / φ 1 are nondecreasing functions and so φ 1     φ 2     φ 3 .
Obviously, Conditions (A 1 ), (A 2 ) and (A 3 ) are fulfilled with these functions and the hypotheses of Theorem 2 as well. Therefore, there exists a real constant r R such that solutions of Problem (57) satisfy lim t u t ln t t 0 α = r .

5. Conclusions

In this article, we considered a fractional integro-differential problem with Hadamard-Type fractional derivatives. The nonlinear source function depends on a lower-order Hadamard fractional derivative of the state and an integral involving another lower-order Hadamard fractional derivative of the state. We assumed the boundedness of the nonlinearities in question by some special kinds of functions in some appropriate spaces. Under these nonlinear growth conditions, we demonstrated that solutions of the initial value fractional problem under consideration are not only bounded by logarithmic functions but actually they converge asymptotically to logarithmic functions.

Author Contributions

Conceptualization, A.M. and N.-E.T.; methodology, A.M. and N.-E.T.; validation, A.M. and N.-E.T.; formal analysis, A.M. and N.-E.T.; investigation, A.M.; writing original draft preparation, A.M.; writing, review and editing, N.-E.T.; Funding acquisition, A.M.; project administration, A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC were funded by King Fahd University of Petroleum and Minerals through project number SB191023.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to acknowledge the support provided by King Fahd University of Petroleum and Minerals (KFUPM) through project number SB191023.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Halanay, A. On the asymptotic behavior of the solutions of an integro-differential equation. J. Math. Anal. Appl. 1965, 10, 319–324. [Google Scholar] [CrossRef] [Green Version]
  2. Cohen, D.S. The asymptotic behavior of a class of nonlinear differential equations. Proc. Am. Math. Soc. 1967, 18, 607–609. [Google Scholar] [CrossRef]
  3. Dannan, F.M. Integral inequalities of Gronwall-Bellman-Bihari type and asymptotic behavior of certain second order nonlinear differential equations. J. Math. Anal. Appl. 1985, 108, 151–164. [Google Scholar] [CrossRef] [Green Version]
  4. Constantin, A. On the asymptotic behavior of second order nonlinear differential equations. Rendiconti di Matematica e delle sue Applicazioni 1993, 13, 627–634. [Google Scholar]
  5. Philos, C.G.; Tsamatos, P.C. Solutions approaching polynomials at infinity to nonlinear ordinary differential equations. Electron. J. Differ. Equ. 2005, 2005, 1–25. [Google Scholar]
  6. Mustafa, O.G.; Rogovchenko, Y.V. Asymptotic integration of a class of nonlinear differential equations. Appl. Math. Lett. 2006, 19, 849–853. [Google Scholar] [CrossRef] [Green Version]
  7. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science Limited: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  8. Magin, R.L.; Abdullah, O.; Baleanu, D.; Zhou, X.J. Anomalous diffusion expressed through fractional order differential operators in the Bloch-Torrey equation. J. Magn. Reson. 2008, 190, 255–270. [Google Scholar] [CrossRef]
  9. Sabatier, J.; Agrawal, O.P.; Machado, J.T. Advances in Fractional Calculus: Theoretical Developments and Applications in Physics and Engineering; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar] [CrossRef]
  10. Shah, N.A.; Elnaqeeb, T.; Animasaun, I.; Mahsud, Y. Insight into the natural convection flow through a vertical cylinder using caputo time-fractional derivatives. Int. J. Appl. Comput. Math. 2018, 4, 1–18. [Google Scholar] [CrossRef]
  11. Ahmad, A.M.; Furati, K.M.; Tatar, N.E. Asymptotic power type behavior of solutions to a nonlinear fractional integro-differential equation. Electron. J. Differ. Equ. 2017, 2017, 1–16. [Google Scholar]
  12. Ahmad, A.M.; Furati, K.M.; Tatar, N.E. Asymptotic Behavior of Solutions for a Class of Fractional Integro-differential Equations. Mediterr. J. Math. 2018, 15, 188. [Google Scholar] [CrossRef]
  13. Băleanu, D.; Mustafa, O.G. On the asymptotic integration of a class of sublinear fractional differential equations. J. Math. Phys. 2009, 50, 123520. [Google Scholar] [CrossRef] [Green Version]
  14. Băleanu, D.; Mustafa, O.G.; Agarwal, R.P. Asymptotic integration of (1+ α)-order fractional differential equations. Comput. Math. Appl. 2011, 62, 1492–1500. [Google Scholar] [CrossRef] [Green Version]
  15. Băleanu, D.; Agarwal, R.P.; Mustafa, O.G.; Coşulschi, M. Asymptotic integration of some nonlinear differential equations with fractional time derivative. J. Phys. A Math. Theor. 2011, 44, 055203. [Google Scholar] [CrossRef]
  16. Brestovanska, E.; Medved’, M. Asymptotic behavior of solutions to second-order differential equations with fractional derivative perturbations. Electron. J. Differ. Equ. 2014, 2014, 1–10. [Google Scholar]
  17. Grace, S.R. On the oscillatory behavior of solutions of nonlinear fractional differential equations. Appl. Math. Comput. 2015, 266, 259–266. [Google Scholar] [CrossRef]
  18. Medveď, M. On the asymptotic behavior of solutions of nonlinear differential equations of integer and also of non-integer order. Electron. J. Qual. Theory Differ. Equ. 2012, 10, 1–9. [Google Scholar]
  19. Medveď, M. Asymptotic integration of some classes of fractional differential equations. Tatra Mt. Math. Publ. 2013, 54, 119–132. [Google Scholar] [CrossRef] [Green Version]
  20. Medved’, M.; Pospíšil, M. Asymptotic Integration of Fractional Differential Equations with Integrodifferential Right-Hand Side. Math. Model. Anal. 2015, 20, 471–489. [Google Scholar] [CrossRef] [Green Version]
  21. Ahmad, A.M.; Furati, K.M.; Tatar, N.E. Boundedness and power-type decay of solutions for a class of generalized fractional Langevin equations. Arab. J. Math. 2019, 8, 79–94. [Google Scholar] [CrossRef] [Green Version]
  22. Ahmad, A.M. On the Asymptotic Behavior of Solutions for a Fractional Differential Equation with a Singular Kernel. In Proceedings of the 24th International Conference on Circuits, Systems, Communications and Computers (CSCC), Chania, Greece, 19–22 July 2020; pp. 130–134. [Google Scholar] [CrossRef]
  23. Ahmad, A.M. The Asymptotic Behavior of Solutions of a Fractional Integro-differential Equation. WSEAS Trans. Syst. Control 2020, 15, 341–348. [Google Scholar] [CrossRef]
  24. Graef, J.R.; Grace, S.R.; Tunç, E. Asymptotic behavior of solutions of nonlinear fractional differential equations with Caputo-Type Hadamard derivatives. Fract. Calc. Appl. Anal. 2017, 20, 71–87. [Google Scholar] [CrossRef]
  25. Kassim, M.D.; Tatar, N.E. Stability of logarithmic type for a Hadamard fractional differential problem. J. Pseudo-Differ. Oper. Appl. 2020, 11, 447–466. [Google Scholar] [CrossRef]
  26. Kassim, M.D.; eddine Tatar, N. Asymptotic behavior of solutions of fractional differential equations with Hadamard fractional derivatives. Fract. Calc. Appl. Anal. 2021, 24, 483–508. [Google Scholar] [CrossRef]
  27. Li, M.; Wang, J. Analysis of nonlinear Hadamard fractional differential equations via properties of Mittag-Leffler functions. J. Appl. Math. Comput. 2016, 51, 487–508. [Google Scholar] [CrossRef]
  28. Li, C.; Li, Z. Stability and Logarithmic Decay of the Solution to Hadamard-Type Fractional Differential Equation. J. Nonlinear Sci. 2021, 31, 1–60. [Google Scholar] [CrossRef]
  29. Ahmad, B.; Alsaedi, A.; Ntouyas, S.K.; Tariboon, J. Hadamard-Type Fractional Differential Equations, Inclusions and Inequalities; Springer: Cham, Switzerland, 2017. [Google Scholar] [CrossRef]
  30. Băleanu, D.; Diethelm, K.; Scalas, E.; Trujillo, J.J. Fractional Calculus: Models and Numerical Methods; World Scientific: Singapore, 2012; Volume 3. [Google Scholar]
  31. De Barra, G. Measure Theory and Integration; Elsevier: Amsterdam, The Netherlands, 2003. [Google Scholar]
  32. Bihari, I. A generalization of a lemma of Bellman and its application to uniqueness problems of differential equations. Acta Math. Hung. 1956, 7, 81–94. [Google Scholar] [CrossRef]
  33. LaSalle, J. Uniqueness theorems and successive approximations. Ann. Math. 1949, 50, 722–730. [Google Scholar] [CrossRef]
  34. Pinto, M. Integral inequalities of Bihari-type and applications. Funkc. Ekvacioj 1990, 33, 387–403. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mugbil, A.; Tatar, N.-E. Hadamard-Type Fractional Integro-Differential Problem: A Note on Some Asymptotic Behavior of Solutions. Fractal Fract. 2022, 6, 267. https://doi.org/10.3390/fractalfract6050267

AMA Style

Mugbil A, Tatar N-E. Hadamard-Type Fractional Integro-Differential Problem: A Note on Some Asymptotic Behavior of Solutions. Fractal and Fractional. 2022; 6(5):267. https://doi.org/10.3390/fractalfract6050267

Chicago/Turabian Style

Mugbil, Ahmad, and Nasser-Eddine Tatar. 2022. "Hadamard-Type Fractional Integro-Differential Problem: A Note on Some Asymptotic Behavior of Solutions" Fractal and Fractional 6, no. 5: 267. https://doi.org/10.3390/fractalfract6050267

APA Style

Mugbil, A., & Tatar, N. -E. (2022). Hadamard-Type Fractional Integro-Differential Problem: A Note on Some Asymptotic Behavior of Solutions. Fractal and Fractional, 6(5), 267. https://doi.org/10.3390/fractalfract6050267

Article Metrics

Back to TopTop