Next Article in Journal
Three Decades of Fuzzy AHP: A Bibliometric Analysis
Previous Article in Journal
Patron–Prophet Artificial Bee Colony Approach for Solving Numerical Continuous Optimization Problems
Previous Article in Special Issue
On Implicit Time–Fractal–Fractional Differential Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Nonlinearity for the Wave Equation with Friction and Viscoelastic Damping

by
Abdelhamid Mohammed Djaouti
* and
Muhammad Amer Latif
Preparatory Year Deanship, King Faisal University, Hofuf 31982, Saudi Arabia
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(10), 524; https://doi.org/10.3390/axioms11100524
Submission received: 31 August 2022 / Revised: 27 September 2022 / Accepted: 28 September 2022 / Published: 1 October 2022
(This article belongs to the Special Issue Fractional Calculus and Differential Equations)

Abstract

:
In this paper we consider a fractional nonlinearity for the wave equation with friction and viscoelastic damping. Using Fixed point theorem a global in time existence of small data solutions to the Cauchy problem is investigated in this research. Our main interest is to show the influence of the fractional nonlinearity parameter to the admissible range of exponent ς comparing with power nonlinearity and also the generating of loss of decay.

1. Introduction and Tools

The left and right Riemann-Liouville fractional derivatives of order a > 0 for a function u are defined as:
Definition 1.
Let u H 1 0 , ϱ 2 , ϱ 2 > ϱ 1 , a > 0 then, the definition of the left and right Riemann-Liouville fractional derivatives are:
R L D ϱ 1 a u z = 1 Γ n a d n d z n ϱ 1 z z λ u λ d λ
and
R L D ϱ 2 a u z = 1 n Γ n a d n d z n z ϱ 2 λ z u λ d λ ,
where n is an integer which satisfies n 1 a < n and Γ · is the Euler’s gamma function.
The left and right Riemann-Liouville fractional integrals of order a > 0 for a function u are defined as:
Definition 2.
Let u H 1 0 , ϱ 2 , ϱ 2 > ϱ 1 , a > 0 then, the definition of the left and right Riemann-Liouville fractional integrals are:
R L T ϱ 1 a u z = 1 Γ a ϱ 1 z z λ u λ d λ
and
R L I ϱ 2 a u z = 1 Γ a z ϱ 2 λ z u λ d λ ,
where Γ · is the Euler’s gamma function.
An important concept, not long ago, that has emerged is the Caputo-Fabrizio integral operator, which has been established in the last few years. This is how it is defined:
Definition 3
([1]). Let u H 1 0 , ϱ 2 , ϱ 2 > ϱ 1 , a 0 , 1 then, the definition of the new Caputo fractional derivative is:
C F D a u z = M a 1 a ϱ 1 z u λ exp a z λ 1 a d λ ,
where M a is normalization function.
The integral formula for the Caputo-Fabrizio fractional derivative is as follows.
Definition 4
([2]). Let u H 1 0 , ϱ 2 , ϱ 2 > ϱ 1 , a 0 , 1 then, the definition of the left and right side of Caputo-Fabrizio fractional integrals are:
ϱ 1 C F I a u z = 1 a B a u z + a B a ϱ 1 z u λ d λ
and
C F I ϱ 2 a u z = 1 a B a u z + a B a z ϱ 2 u λ d λ ,
where B a is normalization function.
Atangana-Baleanu [3] has found a solution to the problem of the Caputo-Fabrizio operator not being reduced to the original function in a special case, despite the fact that the operator is an effective tool in the solution of many systems of differential equations. The features of the Caputo-Fabrizio operator are present in the normalization function.
Some fractional order derivative and integral operators include the power law in their kernel. Nature does not usually exhibit power law behavior. This novel derivative and integral operator incorporates the Mittag-Leffler function [3]. The Mittag-Leffler function is required to model nature. This improved the Atangana-Baleanu operator and piqued researchers’ interest. That the work uses the Atangana-Baleanu operator for Hermite-Hadamard inequalities is unusual. When the parameter is set to zero, the Atangana-Baleanu original function can be derived and compared to the Caputo-Fabrizio results.
Definition 5
([3]). Let u H 1 0 , ϱ 2 , ϱ 2 > ϱ 1 , a 0 , 1 then, the definition of the new fractional derivative is given below
ϱ 1 ABC D z a u z = B a 1 a ϱ 1 z u λ E a a z λ a 1 a d λ .
Definition 6
([3]). Let u H 1 0 , ϱ 2 , ϱ 2 > ϱ 1 , a 0 , 1 then, the definition of the new fractional derivative is given below:
ϱ 1 ABC D z a u z = B a 1 a d d z ϱ 1 z u λ E a a z λ a 1 a d λ .
Equations (8) and (9) have a non-local kernel. Also in (9) when the function is constant we get zero. Associated integral operator for Atangana-Baleanu fractional derivative has been defined as follows.
Definition 7
([3]). The fractional integral associate to the new fractional derivative with non-local kernel of a function u H 1 0 , ϱ 2 as defined:
ϱ 1 AB I a u z = 1 a B a u z + a B a Γ a ϱ 1 z u λ z λ a 1 d λ ,
where ϱ 2 > ϱ 1 , a 0 , 1 .
The authors of [4] described the integral operator’s right hand side as follows:
AB I ϱ 2 a u z = 1 a B a u z + a B a Γ a z ϱ 2 u λ λ z a 1 d λ ,
where ϱ 2 > ϱ 1 , a 0 , 1 . Differential equations of arbitrary real order γ > 0 are used to model various physical models arising in many branches of science and engineering. Applications of such mathematical models can be seen from statistical mechanics and Brownian motion, via visco-elasticity problems, to continuum and quantum mechanics, biosciences, chemical engineering, and control theory, see for instance [5]. Nonlinearity is very important in studying wave equations, evolution equations, damped equations etc. and in some problems fractional nonlinearity is considered as the Riemann-Liouville fractional operator to get a nonlinear memory term also know as nonlinearity. In this paper, we study the global (in time) existence of small data solutions to
ζ λ λ Δ ζ + ζ λ Δ ζ λ = I γ ( ζ ς ) for ( λ , z ) ( 0 , ) × R r , ζ ( 0 , z ) = ζ 0 ( z ) , ζ λ ( 0 , z ) = ζ 1 ( z ) for z R r ,
where
I γ ( ζ ) 0 λ ( λ s ) γ ζ ( s ) d s
denote the fractional Riemann-Liouville integral operator of order γ for some γ ( 0 , 1 ) . We also derive decay estimates for solutions to (12) and show the influence of fractional integral parameter γ . This problem important to the researchers working in field of differential equations and mathematical modeling. Over the last decade several papers have been devoted to the study of semilinear evolution model with the nonlinear term of memory type as in (13). In the pioneering paper [6] the authors determine the critical exponent for the semilinear heat equation with nonlinear memory term. Afterwards, this kind of result has been generalized for fractional (either in space or in time) heat equations [7,8,9] and for weakly coupled system of heat equations [8,10,11]. Mezadek et al. [12] considered the Cauchy problem for the semilinear wave equation with friction damping, visco-elastic damping and power nonlinearity
ζ λ λ Δ ζ + ζ λ Δ ζ λ = ζ ς for ( λ , z ) ( 0 , ) × R r , ζ ( 0 , z ) = ζ 0 ( z ) , ζ λ ( 0 , z ) = ζ 1 ( z ) for z R r ,
where the data ζ 0 and ζ 1 are known as Cauchy data. Mezadek et al. [12] defined a generalized diffusion phenomena and demonstrated that the long time asymptotic of solutions is a mixture of diffusion and wave equation solutions. A more general case also treated recently in [13].
For the semilinear classical damped wave equation with friction damping and no viscoelastic damping, many mathematicians have attempted to solve the following Cauchy problem:
ζ λ λ Δ ζ + ζ λ = ζ ς for ( λ , z ) ( 0 , ) × R r , ζ ( 0 , z ) = ζ 0 ( z ) , ζ λ ( 0 , z ) = ζ 1 ( z ) for z R r .
Energy solutions for compactly supported data ζ 0 , ζ 1 H 1 × L 2 and ς > 1 (and ς 1 + 2 r 2 if r 3 ) were found to exist locally (in time) by Nakao and Ono [14]. Todorova and Yordanov [15] established the global (in time) existence for ς > 1 + 2 r , where the data ζ 0 , ζ 1 H 1 × L 2 and that in terms of nonlinearities, the Fujita exponent is the vital exponent. ζ ς ς > 1 as can be seen from the results of Ikeda et al. [16] and Ikeda et al. [16]. In this case a blow-up result is proved for ς ( 1 , 1 + 2 r ] , even for small data from ( H 1 L 1 ) × ( L 2 L 1 ) , (resp. ( H σ , 0 H 0 , ρ ) × ( H σ 1 , 0 H 0 , ρ ) , where ρ satisfies r 1 a 1 2 < ρ < 2 ς 1 r 2 and a [ 1 , 2 ] . Ikehata and Ohta [17] studied (14) under the assumption ζ 0 , ζ 1 ( H 1 L k ) × ( L 2 L k ) for additional regularity L k , k [ 1 , 2 ) for the data. Ikehata and Ohta obtained a new critical exponent ς c r i λ = 1 + 2 k r for small data Sobolev solutions from both the global (in time) and blow-up perspectives.
Semilinear viscoelastic damped wave equations with viscoelastic damping and without friction have a large body of literature
ζ λ λ Δ ζ Δ ζ λ = ζ ς , ζ ( 0 , z ) = ζ 0 z , ζ λ ( 0 , z ) = ζ 1 z .
Ikehata et al. [18] addressed the related linear Cauchy problem with vanishing right-hand side and demonstrated global well posedness. Shibata [19] provided L ς L β estimations for Sobolev solutions and examined the diffusion phenomenon. Ikehata [20] investigated the asymptotic characteristic of Sobolev solutions as λ , assuming that ζ 0 , ζ 1 ( H 1 L 1 , 1 ) × ( L 2 L 1 , 1 ) , where L 1 , 1 represents a weighted L 1 space. D’Abbicco and Reissig [21] devoted their findings to determining the crucial exponent ς . The critical exponent denotes the existence of global (in time) Sobolev solutions for little data for ς > ς and simply the existence of local (in time) Sobolev solutions for large data for 1 < ς ς . Thus, a blow-up behavior can be expected in general. Since D’Abbicco and Reissig were unable to verify such a critical exponent so the critical exponent remains an open problem. It should be noted that the results of Theorem 2 in D’Abbicco and Reissig [21] are based on the usage of higher order regularity for the data, specifically, second order in space, and on the blending of multiple regularities for the data. The data, in particular, belong to the classical energy space.
Thus, it is feasible to investigate (13) using both friction and viscoelastic damping terms. The subject of the qualitative features of solutions to (13) emerges. It is an interesting point to understand the relationship between friction and viscoelastic damping. As a result, Ikehata and Sawada [22] demonstrated that the frictional damping effect is more prominent than the viscoelastic damping effect for the asymptotic profile as λ . D’Abbicco [23] recently addressed the Cauchy problem (13), in which the data ζ 0 , ζ 1 are assumed to belong to the energy space with an additional L 1 regularity, namely to ( H 1 L 1 ) × ( L 2 L 1 ) . The author established that small data energy solutions exist globally (in time) for admissible exponents ς 1 + 2 r , 1 + 2 r 2 for r > 1 . Ikehata and Takeda [24] investigated the identical Cauchy problem (13) with the following data assumptions ζ 0 , ζ 1 ( W r 2 + ε , 1 W r 2 + ε , ) × ( L 1 L ) , where ε > 0 . The authors established a result regarding the global (in time) existence of small data Sobolev solutions exclusively for r = 1 , 2 , 3 . They obtained this result for the range of permissible exponents ς meeting the condition ς > 1 + 2 r . Mezadek et al. [12] examined the Cauchy problem (13) under particular data and dimension r assumptions and investigated the effect of the regularity parameters σ 1 , σ 2 R + , and the additional regularity parameter k 1 , 2 on the data ζ 0 , ζ 1 ( H σ 1 L k ) × ( H σ 2 L k ) on the acceptable ranges of exponents ς , which enables the global (in time) existence of small data Sobolev or energy solutions with a sufficient degree of regularity. Additionally, the authors investigated the effect of σ 1 , σ 2 and k on solution regularity. Mezadek et al. [12] demonstrated the global (in time) existence of small data solutions to the semilinear Cauchy problem (13) in any space dimension r 1 by applying estimates of linear Cauchy problem solutions to the semilinear Cauchy problem with power nonlinearity ζ ς .
Small data solutions to the Cauchy problem for a semilinear wave equation with friction, viscoelastic damping and a fractional nonlinearity are the main goal of this paper. We want to find out if these solutions are globally available in time for small data. The main objectives of this paper is to show that fractional nonlinearity has an effect on the range of exponent ς that can be used and it also causes the decay rate to slow down with respect to the solution to the corresponding linear problem.

2. Main Results

2.1. Strategies

For the Cauchy problem (12), we will prove several results here. Our primary goal is to demonstrate the worldwide viability of small data solutions in the near future. As a result, we can assume that the zero solution is stable right away.
We introduce for σ 1 σ 2 0 and k [ 1 , 2 ) the function space
A k σ 1 , σ 2 : = ( H σ 1 L k ) × ( H σ 2 L k ) , A σ 1 , σ 2 = A 1 σ 1 , σ 2
with the norm
( ζ , υ ) A k σ 1 , σ 2 : = ζ H σ 1 + ζ L k + υ H σ 2 + υ L k .
In the first instance, we suppose that the low regular data come from A σ 1 , 0 . We will demonstrate that Sobolev solutions exist globally in time.
C ( [ 0 , ) , H σ 1 ) C 1 ( [ 0 , ) , L 2 )
where for large dimensions, an upper bound on the power nonlinearity exponent ς is necessary because it is larger than the modified Fujita exponent.
Case two outcomes are related to the highly regular data, which implies they belong to A k σ 1 , σ 2 and have a σ 2 > 0 . We use tools from harmonic analysis from [25] to prove the existence of a global (in time) solution (see Appendix A). All these results require the condition ς > σ 2 + 1 if σ 2 ( 0 , r 2 ] . Here, we denote by σ 2 : = min { a Z : σ 2 a } the ceiling function in σ 2 .
Finally, if σ 2 > r 2 , then using fractional powers the last condition ς > σ 2 will be weakened to ς > σ 2 .

2.2. Low Regular Data

2.2.1. Low Dimension

Theorem 1.
Let r < 2 ( 2 σ 1 ) and σ 1 ( 0 , 2 ) be a real number and the data ( ζ 0 , ζ 1 ) are in A σ 1 , 0 . Suppose that the exponent ς satisfies
ς > 2 + r r + 2 ( γ 1 )
and
2 ς i f σ 1 r 2 , 2 ς ς G N , σ 1 ( r ) : = r r 2 σ 1 i f σ 1 < r 2 ,
thus, there exists a small ϵ such that if
( ζ 0 , ζ 1 ) A σ 1 , 0 ϵ ,
then there exists a uniquely determined globally (in time) energy solution to (12) in
C ( [ 0 , ) , H σ 1 ) C 1 ( [ 0 , ) , L 2 ) .
According to estimates, the solution satisfies
ζ ( λ , · ) L 2 ( 1 + λ ) r 4 + ( 1 γ ) ( ζ 0 , ζ 1 ) A σ 1 , 0 , | D | σ 1 ζ ( λ , · ) L 2 ( 1 + λ ) r 4 σ 1 2 + ( 1 γ ) ( ζ 0 , ζ 1 ) A σ 1 , 0 , ζ λ ( λ , · ) L 2 ( 1 + λ ) r 4 1 + ( 2 γ ) ( ζ 0 , ζ 1 ) A σ 1 , 0 .
Example 1.
In order to determine the permissible range for the fractional integral parameter γ, we look at the conditions (16) and (17) in some particular situations of low dimension that are dependent on the parameters r , k , σ and γ.
rRegularity σ1Admissible range for ς
r = 1 σ 1 1 2 , 3 2 3 2 γ 1 < ς < .
σ 1 2 γ 3 , 1 2 3 2 γ 1 < ς 1 1 2 σ 1 .
r = 2 σ 1 1 , 2 2 γ < ς <
σ 1 1 γ 2 , 1 2 γ < ς 1 1 σ 1 .
Example 2.
In this example we treat the case of third dimension r = 3 for fixed γ = 1 2 . The model that we have is given by
ζ λ λ Δ ζ + ζ λ Δ ζ λ = I 1 2 ( ζ ς ) for ( λ , z ) ( 0 , ) × R 3 ,
ζ ( 0 , z ) = ζ 0 ( z ) , ζ λ ( 0 , z ) = ζ 1 ( z ) for z ( 0 , ) × R 3 .
Then, the admissible range for ς can be described as follows:
5 2 < ς i f σ 1 3 2 , 2 , 5 2 < ς 3 3 2 σ 1 i f σ 1 9 10 , 3 2 .

2.2.2. Higher Dimension

Theorem 2.
Let r > 4 and σ 1 2 be a real number and the data ( ζ 0 , ζ 1 ) are in A σ 1 , 0 , where k [ 1 , 2 ) . Suppose that the exponent ς satisfies
ς > 1 γ
and
2 ς i f σ 1 r 2 , 2 ς ς G N , σ 1 ( r ) i f σ 1 < r 2 ,
thus, there exists a small ϵ such that if
( ζ 0 , ζ 1 ) A σ 1 , 0 ϵ ,
then there exists a uniquely determined globally (in time) energy solution to (12) in
C ( [ 0 , ) , H σ 1 ) C 1 ( [ 0 , ) , L 2 ) .
According to estimates, the solution satisfies
ζ ( λ , · ) L 2 ( 1 + λ ) r 4 + A ζ ( ζ 0 , ζ 1 ) A σ 1 , 0 , | D | σ 1 ζ ( λ , · ) L 2 ( 1 + λ ) r 4 σ 1 2 + A | D | σ 1 v ( ζ 0 , ζ 1 ) A σ 1 , 0 , ζ λ ( λ , · ) L 2 ( 1 + λ ) r 4 1 + A ζ λ ( ζ 0 , ζ 1 ) A σ 1 , 0 ,
where A ζ r 4 γ , A | D | σ 1 ζ = A ζ + σ 1 2 and A ζ λ = A ζ + 1 define the loss of decay in comparison with the corresponding decay estimates for the solution ζ to the linear Cauchy problem with vanishing right-hand side.
In the following theorems we will take the additional regularity L k since the goal is the global existence where there is no blow-up results exist in the literature and also use the last case of estimates (22) in Proposition 1.

2.3. Data from Sobolev Spaces with Suitable Regularity

Theorem 3.
Let r < 2 k 2 k ( 2 σ 1 ) and the data ( ζ 0 , ζ 1 ) are supposed to belong to A k σ 1 , σ 2 , where k [ 1 , 2 ) , σ 2 r 2 . Suppose that for the exponent ς :
ς > max 1 + r 2 k r 2 k 1 + γ ; σ 2 + 1 ,
and
1 + 2 ( σ 1 σ 2 ) r ς i f σ 1 r 2 , r 2 + σ 2 , 1 + 2 ( σ 1 σ 2 ) r ς 1 + min 2 r 2 σ 1 ; 2 ( σ 1 σ 2 ) r 2 σ 1 i f σ 1 0 , r 2 .
there exists a small ϵ such that, if
( ζ 0 , ζ 1 ) A k σ 1 , σ 2 ϵ ,
and there exists a uniquely determined globally (in time) energy solution to (12) in
C ( [ 0 , ) , H σ 1 ) C 1 ( [ 0 , ) , H σ 2 ) .
Furthermore, the solution satisfies:
ζ ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 + 1 γ ( ζ 0 , ζ 1 ) A k σ 1 , 0 , | D | σ 1 ζ ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 σ 1 2 + 1 γ ( ζ 0 , ζ 1 ) A k σ 1 , 0 , ζ λ ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 + 1 γ ( ζ 0 , ζ 1 ) A k σ 1 , 0 , | D σ 2 ζ λ ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 σ 2 2 + 1 γ ( ζ 0 , ζ 1 ) A k σ 1 , 0 .

2.4. Large Regular Data

Theorem 4.
Let r < 2 k 2 k ( 2 σ 1 ) and the data ( ζ 0 , ζ 1 ) are supposed to belong to A k σ 1 , σ 2 , where k [ 1 , 2 ) , σ 2 > r 2 . Suppose that for the exponent ς : , we have
ς > max 1 + r 2 k r 2 k 1 + γ ; σ 2 ; 2 .
Then for small ϵ, we have
( ζ 0 , ζ 1 ) A k σ 1 , σ 2 ϵ ,
and there exists a uniquely determined globally (in time) energy solution to (12) in
C ( [ 0 , ) , H σ 1 ) C 1 ( [ 0 , ) , H σ 2 ) .
Moreover, the solution meets the estimates:
ζ ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 + 1 γ ( ζ 0 , ζ 1 ) A k σ 1 , 0 , | D | σ 1 ζ ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 σ 1 2 + 1 γ ( ζ 0 , ζ 1 ) A k σ 1 , 0 , ζ λ ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 + 1 γ ( ζ 0 , ζ 1 ) A k σ 1 , 0 , | D | σ 2 ζ λ ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 σ 2 2 + 1 γ ( ζ 0 , ζ 1 ) A k σ 1 , 0 .
Remark 1.
If we take in Theorems 3 and 4 the data belongs to A σ 1 , σ 2 , then we have to take r > 2 ( 2 σ 1 ) and the decay of the solution and their derivative will be define by γ as power similarly to Theorem 2.

3. Philosophy of Our Approach and Proofs

In order to proceed to our proofs, we first introduce some tools and previous results.
Proposition 1.
Let a R , b > 1 and γ ( 0 , 1 ) . Then,
0 λ ( 1 + λ τ ) a 0 τ ( τ s ) γ ( 1 + s ) b d s d τ ( 1 + λ ) γ i f a > 1 , ( 1 + λ ) γ log ( 2 + λ ) i f a = 1 , ( 1 + λ ) γ + 1 a i f a < 1 .
Proposition 2.
Let a R and b , γ ( 0 , 1 ) . Then,
0 λ ( 1 + λ τ ) a 0 τ ( τ s ) γ ( 1 + s ) b d s d τ ( 1 + λ ) γ + 1 b i f a > 1 , ( 1 + λ ) γ + 1 b log ( 2 + λ ) i f a = 1 , ( 1 + λ ) γ + 2 a b i f a < 1 .
Reference [26] contains the evidence supporting Propositions 1 and 2.
Proposition 3.
Let ( ζ 0 , ζ 1 ) ( H σ 1 L 1 ) × ( H σ 2 L 1 ) with σ 2 + 2 σ 1 σ 2 0 , k [ 1 , 2 ) . Then the solution to the Cauchy problem
ζ λ λ Δ ζ ζ λ Δ ζ λ = o , ζ ( 0 , z ) = ζ 0 ( z ) , ζ λ ( 0 , z ) = ζ 1 ( z )
satisfies the decay estimates
ζ ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 ( ζ 0 , ζ 1 ) A σ 1 , 0 ,
| D | σ 2 ζ λ ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 σ 2 2 1 ( ζ 0 , ζ 1 ) A σ 1 , 0 ,
ζ λ ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 1 ( ζ 0 , ζ 1 ) A σ 1 , 0 ,
| D | σ 1 ζ ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 σ 1 2 ( ζ 0 , ζ 1 ) A σ 1 , 0 .
Proof of Proposition 3 can be found in [12].
We define the space of solutions χ ( λ ) by
χ ( λ ) = C ( [ 0 , λ ] , H σ 1 ) C 1 ( [ 0 , λ ] , H σ 2 ) ,
where the norm of χ ( λ ) will be proposed separately for each theorem. We introduce the operator N by
N : ζ χ ( λ ) N ζ = N ζ ( λ , z ) : = ζ l n ( λ , z ) + ζ n l ( λ , z ) .
We denote by E 0 ( λ , 0 , z ) and E 1 ( λ , 0 , z ) the fundamental solutions to the linear equation, namely
ζ l n ( λ , z ) : = E 0 ( λ , 0 , z ) ( z ) ζ 0 ( z ) + E 1 ( λ , 0 , z ) ( z ) ζ 1 ( z )
is a solution to the Cauchy problem
ζ λ λ Δ ζ ζ λ Δ ζ λ = 0 , ζ ( 0 , z ) = ζ 0 ( z ) , ζ λ ( 0 , z ) = ζ 1 ( z ) ,
and
ζ n l ( λ , z ) : = 0 λ E 1 ( λ , τ , z ) ( z ) 0 τ ( τ s ) γ | ζ ( s , χ ) | ς d s d τ
is a solution to the Cauchy problem
ζ λ λ Δ ζ ζ λ Δ ζ λ = I γ ( ζ ς ) , ζ ( 0 , z ) = 0 , ζ λ ( 0 , z ) = 0 .
As a result of Proposition A4, we may prove the following inequalities:
N ζ χ ( λ ) C 0 ( λ ) ( ζ 0 , ζ 1 ) A σ 1 , σ 2 + C 1 ( λ ) ζ χ ( λ ) ς ,
N ζ N ζ χ ( λ ) C 2 ( λ ) ζ ζ χ ( λ ) ζ χ ( λ ) ς 1 + ζ χ ( λ ) ς 1 ,
where C 1 ( λ ) , C 2 ( λ ) 0 for λ 0 and C 1 ( λ ) , C 2 ( λ ) C for all λ [ 0 , ) .

3.1. Proof of Theorem 1

Let
ζ χ ( λ ) = sup τ [ 0 , λ ] | η | + J 1 ( 1 + τ ) r 4 + | η | 2 σ 1 + J ( 1 + J γ ) z η λ J ζ ( τ , · ) L 2 ,
where J + | η | = 0 , 1 . We remark that if ζ χ ( λ ) , then ζ χ ( τ ) ζ χ ( λ ) for any 0 τ λ .
We begin the proof of (29). From the estimates (24) to (27) of Theorem 3 and the definition of the norm of solutions space χ ( λ ) we have
ζ l n χ ( λ ) = sup τ [ 0 , λ ] { ( 1 + τ ) r 4 1 + γ ζ l n ( τ , · ) L 2 + ( 1 + τ ) r 4 + σ 1 2 1 + γ | D | σ 1 ζ l n ( τ , · ) L 2 ( R r ) + ( 1 + τ ) r 4 1 + γ ζ λ l n ( τ , · ) L 2 ( R r ) } sup τ [ 0 , λ ] { ( 1 + τ ) r 4 1 + γ ( 1 + τ ) r 4 ( ζ 0 , ζ 1 ) A σ 1 , 0 + ( 1 + τ ) r 4 + σ 1 2 1 + γ ( 1 + τ ) r 4 σ 1 2 ( ζ 0 , ζ 1 ) A σ 1 , 0 + ( 1 + τ ) r 4 1 + γ ( 1 + τ ) r 4 1 ( ζ 0 , ζ 1 ) A σ 1 , 0 } ( ζ 0 , ζ 1 ) A σ 1 , 0 .
Consequently,
ζ l n χ ( λ ) ( ζ 0 , ζ 1 ) A σ 1 , 0 .
Simply proving (29), it suffices to prove:
ζ n l χ ( λ ) ζ χ ( λ ) ς .
For the term ζ n l L 2 we have
ζ n l L 2 0 λ ( 1 + λ τ ) r 4 0 τ ( τ s ) γ | ζ ( τ , · ) | ς ] L 1 L 2 d s d τ .
Thus, we have
| ζ ( τ , · ) | ς L 1 ( R r ) L 2 ( R r ) | ζ ( τ , z ) | ς L 1 + | ζ ( τ , z ) | ς L 2 .
Proposition A1 and Gagliardo-Nirenberg inequality can be used to estimate both of the right-hand side terms. For the first term, we are able to get
| ζ ( τ , · ) | ς L 1 = R r | ζ ( τ , z ) | ς d χ 1 ς ς = ζ ( τ , · ) L ς ( R r ) ς ζ ( τ , · ) L 2 ( R r ) ( 1 θ ) ς | D | σ ζ ( τ , · ) L 2 θ ς ,
where
θ = r s 1 2 1 ς [ 0 , 1 ]
which is due to the condition (17) for ς .
By using the norm of solutions space χ ( λ ) for 0 τ λ , we get
ζ ( τ , · ) L k ς ( R r ) ς ( 1 + τ ) ( 1 θ ) ς r 4 + 1 γ + θ ς r 4 s 2 + 1 γ ζ χ ( λ ) ς .
Then
ζ ( τ , · ) L k ς ( R r ) ς ( 1 + τ ) r 2 ς + r 2 + ( 1 γ ) ς ζ χ ( λ ) ς .
Estimating | ζ ( τ , · ) | ς L 2 follows the same principles. This is done by employing the Gagliardo-Nirenberg inequalities and the defining of the norm of the solution space χ ( λ ) .
| ζ ( τ , · ) | ς L 2 = R r | ζ ( τ , z ) | 2 ς d χ 1 2 ς ς = ζ ( τ , · ) L 2 ς ( R r ) ς ζ ( τ , · ) L 2 ( R r ) ς ( 1 θ ˜ ) | D | σ ζ ( τ , · ) L 2 ς θ ˜ ( 1 + τ ) ( 1 θ ˜ ) ς r 4 + 1 γ + θ ˜ ς r 4 σ 2 + 1 γ ζ χ ( λ ) ς ( 1 + τ ) r 2 ς + r 4 + ( 1 γ ) ς ζ χ ( λ ) ς ,
where
θ ˜ = r σ 1 2 1 2 ς [ 0 , 1 ]
this is from condition (17) for ς . Hence, we may conclude for 0 τ λ the following estimate:
u ( τ , · ) L 2 ς ς ( 1 + τ ) r 2 ς + r 4 + ( 1 γ ) ς ζ χ ( λ ) ς .
All together leads to
ζ n l L 2 0 λ ( 1 + λ τ ) r 4 0 τ ( τ s ) γ ( 1 + τ ) r 2 ς + r 2 + ( 1 γ ) ς ζ χ ( λ ) ς d s d τ ζ χ ( λ ) ς ( 1 + λ ) r 4 + 1 γ ,
Moreover, we have to assume the condition r 2 ς + r 2 + ( 1 γ ) ς < 1 which generate a Fujita-like upper bound (16). So, it follows the desired estimate
ζ n l L 2 ζ χ ( λ ) ς ( 1 + λ ) r 4 + 1 γ .
In the same way one can get
| D | σ 1 ζ n l L 2 ζ χ ( λ ) ς ( 1 + λ ) r 4 + σ 1 2 + 1 γ ,
and
ζ λ n l L 2 ζ χ ( λ ) ς ( 1 + λ ) r 4 + 1 γ .
Taking into consideration (31) and (35) to (37) the estimate (29) is proved.
To prove (30) we assume that ζ and ξ belong to χ ( λ ) . Then
N ζ N ξ = 0 λ E 1 ( λ , τ , z ) ( z ) 0 τ ( τ s ) γ | ζ ( τ , z ) | ς | ξ ( τ , z ) | ς d s d τ .
We control all norms appearing in N ζ N ξ χ ( λ ) . We have
N ζ N ξ L 2 0 λ ( 1 + λ τ ) r 4 0 τ ( τ s ) γ | ζ ( τ , z ) | ς | ξ ( τ , z ) | ς L 1 L 2 d s d τ .
Hölder’s inequality implies
| ζ ( τ , · ) | ς | ξ ( τ , · ) | ς L 2 ζ ( τ , · ) ξ ( τ · ) L 2 ς ζ ( τ , · ) L 2 ς ς 1 + ξ ( τ , · ) L 2 ς ς 1 ,
and
| ζ ( τ , · ) | ς | ξ ( τ , z ) | ς L 1 ζ ( τ , · ) ξ ( τ · ) L ς ζ ( τ , · ) L ς ς 1 + ξ ( τ , · ) L ς ς 1 .
By using the norm of solution space χ ( λ ) and after applying the classical Gagliardo-Nirenberg inequality as we did for (33) and (34) we obtain the following estimates for 0 τ λ :
ζ ( τ , · ) ξ ( τ · ) L 2 ς ( 1 + τ ) r 2 + r 4 ς ζ ζ χ ( λ ) , ζ ( τ , · ) L 2 ς ς 1 ( 1 + τ ) r 2 + r 4 ς ( ς 1 ) ζ χ ( λ ) ς 1 , ξ ( τ , · ) L 2 ς ς 1 ( 1 + τ ) r 2 + r 4 ς ( ς 1 ) ξ χ ( λ ) ς 1 , ζ ( τ , · ) ξ ( τ , · ) L ς ( 1 + τ ) r 2 + r 2 ς ζ ξ χ ( λ ) , ζ ( τ , · ) L ς ς 1 ( 1 + τ ) r 2 + r 2 ς ( ς 1 ) u χ ( λ ) ς 1 , ξ ( τ , · ) L ς ς 1 ( 1 + τ ) r 2 + r 2 ς ( ς 1 ) ξ χ ( λ ) ς 1 .
Then we get
| ζ ( τ , · ) | ς | ξ ( τ , z ) | ς L 2 ( 1 + τ ) r 2 ς + r 4 ζ ξ χ ( λ ) ζ χ ( λ ) ς 1 + ξ χ ( λ ) ς 1 ,
and
| ζ ( τ , · ) | ς | ξ ( τ , z ) | ς L 1 ( 1 + τ ) r 2 ς + r 2 ζ ξ χ ( λ ) ζ χ ( λ ) ς 1 + ξ χ ( λ ) ς 1 .
Applying the same ideas as we did to estimate | D | σ ζ n l ( λ , · ) L 2 , this means, after plugging (41) and (42) into (38) one can get after using (16) the following estimates:
| D | σ ( N ζ N ξ ) ( λ , · ) L 2 ( 1 + λ ) r 4 σ 2 ζ ξ χ ( λ ) ζ χ ( λ ) ς 1 + ξ χ ( λ ) ς 1 ,
and
( N ζ N ξ ) ( λ , · ) L 2 ( 1 + λ ) r 4 ζ ξ χ ( λ ) ζ χ ( λ ) ς 1 + ξ χ ( λ ) ς 1 .
Then from the definition of χ ( λ ) , the proof of (30) is completed.

3.2. Proof of Theorem 2

χ ( λ ) = C ( [ 0 , λ ] , H σ 1 ) C 1 ( [ 0 , λ ] , L 2 ) ,
with the norm
ζ χ ( λ ) = sup τ [ 0 , λ ] | η | + J 1 ( 1 + τ ) r 4 A ζ + | η | 2 σ 1 + J z η λ J ζ ( τ , · ) L 2 ,
where J + | η | = 0 , 1 . As a result of the introduced norm for the solution space’s elements χ ( λ ) and the estimates of solutions to the homogeneous damped wave equation, we can derive (31). By demonstrating the inequality, we complete the prior theorem’s proofs.
ζ n l χ ( λ ) ζ χ ( λ ) ς
  • For ζ n l L 2 we have
    ζ n l L 2 0 λ ( 1 + λ τ ) r 0 τ ( τ s ) γ | ζ ( τ , · ) | ς L 1 L 2 d s d τ .
    Similarly to (33) and under the same conditions described in the theorem we get
    | ζ ( τ , · ) | ς L 1 L 2 ζ χ ( λ ) ς ( 1 + τ ) r 2 ς + r + A ζ ς r 4 ς .
    Using the last estimate we obtain
    ζ n l L 2 ζ χ ( λ ) ς 0 λ ( 1 + λ τ ) r 4 0 τ ( τ s ) γ ( 1 + τ ) r 2 ς + r + A ζ ς r 4 ς d s d τ ζ χ ( λ ) ς ( 1 + λ ) γ ζ χ ( λ ) ς ( 1 + λ ) r 4 + A ζ .
    Hence
    ζ n l L 2 ζ χ ( λ ) ς ( 1 + λ ) r 4 + A ζ .
  • For | D | σ 1 ζ n l L 2 we have
    | D | σ 1 ζ n l L 2 0 λ ( 1 + λ τ ) r 4 σ 1 2 0 τ ( τ s ) γ | ζ ( τ , · ) | ς L 1 L 2 d s d τ .
    Using (46) we get
    | D | σ 1 ζ n l L 2 ζ χ ( λ ) ς 0 λ ( 1 + λ τ ) r 4 σ 1 2 0 τ ( τ s ) r 2 ς + r + A ζ ς r 4 ς d s d τ ζ χ ( λ ) ς ( 1 + λ ) γ ζ χ ( λ ) ς ( 1 + λ ) r 4 + A ζ ζ χ ( λ ) ς ( 1 + λ ) r 4 σ 1 2 + A | D | σ 1 ζ .
    Hence
    | D | σ 1 ζ n l L 2 ζ χ ( λ ) ς ( 1 + λ ) r 4 σ 1 2 + A | D | σ 1 ζ .
  • For ζ λ n l L 2 analogously to (47) and (48) we get
    ζ λ n l L 2 ζ χ ( λ ) ς ( 1 + λ ) r 4 1 + A ζ λ .
From (47)–(49) the prove of (45) is completed.
To control all norms appearing in N ζ N ξ χ ( λ ) we follow exactly the same steps used in the proof of previous theorem.

3.3. Proof of Theorem 3

We define the space of solutions χ ( λ ) by
χ ( λ ) = C ( [ 0 , ) , H σ 1 ) C 1 ( [ 0 , ) , H σ 2 ) .
Our goal is to prove (29) and (30). Let
ζ χ ( λ ) = sup τ [ 0 , λ ] ( 1 + τ ) r 2 1 k 1 2 ( 1 γ ) u ( τ , · ) L 2 + ( 1 + τ ) r 2 1 k 1 2 + 1 ( 2 γ ) ζ λ ( τ , · ) L 2 + ( 1 + τ ) r 2 1 k 1 2 + 1 + σ 2 2 ( 2 γ ) ζ λ ( τ , · ) H σ 2 + ( 1 + τ ) r 2 1 k 1 2 + σ 1 2 ( 1 γ ) ζ ( τ , · ) H σ 1 .
From the defined norm and the estimates (24) to (27) we get
ζ L n χ ( λ ) ( ζ 0 , ζ 1 ) A k σ 1 , σ 2 .
Then, its remains to prove (32) to get (29).
For ζ n L L 2 we have
ζ n L L 2 0 λ ( 1 + λ τ ) r 2 1 k 1 2 0 τ ( τ s ) γ | ζ ( τ , · ) | ς L k H σ 2 d s d τ .
Similarly to (33) and (34) we obtain
ζ ( τ , · ) L k ς ς ( 1 + τ ) r 2 k ς + r 2 k + ( 1 γ ) ς ζ χ ( λ ) ς
and
ζ ( τ , · ) L 2 ς ς ( 1 + τ ) r 2 k ς + r 4 + ( 1 γ ) ς ζ χ ( λ ) ς ,
provided that
2 k ς i f σ 1 r 2 , 2 k ς r r 2 σ 1 i f σ 1 < r 2 .
Still it remains to estimate the norm
| ζ ( τ , · ) | ς H σ 2 .
The fractional chain rule, which the reader can find in citation [25] or Appendix A, is therefore used as a harmonic analysis tool in this paper. We can estimate ς > σ 2 and 0 τ λ by considering the Propositions A1 and A3, in particular Formula (A2).
| ζ ( τ , · ) | ς H σ 2 ζ ( τ , · ) L q 1 ς 1 | D | σ 2 ζ ( τ , · ) L q 2 ζ ( τ , · ) L 2 ( ς 1 ) ( 1 θ 1 ) | D | σ 1 ζ ( τ , · ) L 2 ( ς 1 ) θ 1 ζ ( τ , · ) L 2 1 θ 2 | D | σ 1 ζ ( τ , · ) L 2 θ 2 ( 1 + τ ) r 2 k ς + r 4 σ 2 2 + ( 1 γ ) ς ζ χ ( λ ) ς .
Then
| ζ ( τ , · ) | ς H σ 2 ( 1 + τ ) r 2 k ς + r 4 σ 2 2 + ( 1 γ ) ς ζ χ ( λ ) ς ,
where
ς 1 q 1 + 1 q 2 = 1 2 , θ 1 = r σ 1 1 2 1 q 1 [ 0 , 1 ] , θ 2 = r σ 1 1 2 1 q 2 + σ 2 σ 1 σ 2 σ 1 , 1 .
Last conditions implies
2 q 1 and 2 q 2 if σ 1 r 2 + σ 2 , , 2 q 1 and 2 q 2 2 r r 2 ( σ 1 σ 2 ) if σ 1 r 2 , r 2 + σ 2 , 2 q 1 2 r r 2 σ 1 and 2 q 2 2 r r 2 ( σ 1 σ 2 ) if σ 1 0 , r 2 .
Then we get the following bounds:
1 + 2 ( σ 1 σ 2 ) r ς if σ 1 r 2 , r 2 + σ 2 , 1 + 2 ( σ 1 σ 2 ) r ς 1 + 2 r 2 σ 1 if σ 1 0 , r 2 .
Using (51)–(53) with the last estimates of (22) in (50), we get
ζ n L L 2 0 λ ( 1 + λ τ ) r 2 1 k 1 2 0 τ ( τ s ) γ ( 1 + τ ) r 2 k ς + r 2 k + ( 1 γ ) ς ζ χ ( λ ) ς d s d τ ζ χ ( λ ) ς ( 1 + λ ) r 2 1 k 1 2 + 1 γ ,
where r 2 1 k 1 2 < 1 from the restriction of dimension described in the theorem and using the conditions r 2 k ς + r 2 k + ( 1 γ ) ς > 1 which generate a Fujita-like upper bound (16).
Finally, we get
ζ n L L 2 ζ χ ( λ ) ς ( 1 + λ ) r 2 1 k 1 2 + 1 γ .
Similarly, we can get
ζ λ n L L 2 ζ χ ( λ ) ς ( 1 + λ ) n 2 1 k 1 2 + 1 γ .
| D | σ 2 ζ λ n L L 2 ζ χ ( λ ) ς ( 1 + λ ) r 2 1 k 1 2 + σ 2 2 + 1 γ .
| D | σ 1 ζ n L L 2 ζ χ ( λ ) ς ( 1 + λ ) r 2 1 k 1 2 + σ 1 2 + 1 γ .
Last estimates complete the proof of (32). Now we prove the second inequality (30). Let ζ and ξ belong to χ ( λ ) . Then
N ζ N ξ L 2 0 λ ( 1 + λ τ ) r 2 1 k 1 2 0 τ ( τ s ) γ | ζ ( τ , z ) | ς | ξ ( τ , z ) | ς L k H σ 2 d s d τ .
Similar to (41) and (42), we get
| ζ ( τ , · ) | ς | ξ ( τ , z ) | ς L 2 ( 1 + τ ) r 2 k ς + r 4 ζ ξ χ ( λ ) ζ χ ( λ ) ς 1 + ξ χ ( λ ) ς 1 ,
| ζ ( τ , · ) | ς | ξ ( τ , z ) | ς L k ( 1 + τ ) r 2 k ς + r 2 k ζ ξ χ ( λ ) ζ χ ( λ ) ς 1 + ξ χ ( λ ) ς 1 .
In the next step we may control | ζ ( τ , · ) | ς | ξ ( τ , z ) | ς H σ 2 ( R r ) . Indeed, using the fractional Leibniz rule from Proposition A2 we get
| ζ ( τ , · ) | ς | ξ ( τ , · ) | ς H σ 2 0 1 | D | σ 2 ( ζ ξ ) ( ζ ι ( ζ ξ ) ) | ζ ι ( ζ ξ ) | ς 2 L 2 d r 0 1 | D | σ 2 ( ζ ξ ) L q 1 ( ζ ι ( ζ ξ ) ) | ζ ι ( ζ ξ ) | ς 2 L q 2 d r + 0 1 ζ ξ L q 3 | D | σ 2 ( ζ ι ( ζ ξ ) ) | ζ ι ( ζ ξ ) | ς 2 L q 4 d r ,
where
1 2 = 1 q 1 + 1 q 2 = 1 q 3 + 1 q 4 .
For the first integral we use the classical Gagliardo-Nirenberg inequality and obtain for 0 τ λ
| D | σ 2 ( ζ ξ ) L q 1 ζ ξ L 2 1 θ 1 | D | σ 1 ( ζ ξ ) L 2 θ 1 ( 1 + τ ) r 2 k σ 2 2 + r 2 q 1 + ( 1 γ ) ζ ξ χ ( λ ) ,
and
( ζ ι ( ζ ξ ) ) | ζ ι ( ζ ξ ) | ς 2 L q 2 ζ ι ( ζ ξ ) L 2 ( 1 θ 2 ) ( ς 1 ) | D | σ 1 ζ ι ( ζ ξ ) L 2 θ 2 ( ς 1 ) ( 1 + τ ) r 2 k ( ς 1 ) + r 2 q 2 ζ ι ( ζ ξ ) χ ( λ ) ς 1
for
θ 1 = r σ 1 1 2 1 q 1 + σ 2 σ 1 σ 2 σ 1 , 1 , θ 2 = r σ 1 1 2 1 q 2 ( ς 1 ) [ 0 , 1 ] .
The last conditions implies that
2 q 1 and 2 ς 1 q 2 if σ 1 r 2 + σ 2 , , 2 q 1 and 2 ς 1 q 2 2 r r 2 ( σ 1 σ 2 ) if σ 1 r 2 , r 2 + σ 2 , 2 q 1 2 r r 2 ( σ 1 σ 2 ) and 2 ς 1 q 2 2 r r 2 ( σ 1 σ 2 ) if σ 1 0 , r 2 .
Thus we get the following bounds:
To estimate the first term in the second integral we use again the Gagliardo-Nirenberg inequality.
In this way we, obtain
ζ ξ L q 3 ζ ξ L 2 1 θ 3 | D | σ 1 ( ζ ξ ) L 2 θ 3 ( 1 + τ ) r 2 k + r 2 q 3 + 1 γ ζ ξ χ ( λ ) ,
where
θ 3 = r σ 1 1 2 1 q 3 [ 0 , 1 ] .
To estimate the second term we use the fractional chain rule from Proposition A3.
Hence, we get
| D | σ 2 ( ζ ι ( ζ ξ ) ) | ζ ι ( ζ ξ ) | ς 2 L q 4 ζ ι ( ζ ξ ) L q 5 ς 2 | D | σ 2 ( ζ ι ( ζ ξ ) ) L q 6 ,
where
1 q 4 = ς 2 q 5 + 1 q 6 .
Using Gagliardo-Nirenberg inequality to estimate the last two norms, we get
ζ ι ( ζ ξ ) L q 5 ς 2 ζ ι ( ζ ξ ) L 2 ( p 2 ) ( 1 θ 5 ) | D | σ 1 ( ζ ι ( ζ ξ ) ) L 2 ( p 2 ) θ 5 ( 1 + τ ) r 2 k + r 2 q 5 + 1 γ ( p 2 ) ζ ι ( ζ ξ ) χ ( λ ) ς 2 ,
and
| D | σ 2 ( ζ ι ( ζ ξ ) ) L q 6 ζ ι ( ζ ξ ) L 2 1 θ 6 | D | σ 1 ( ζ ι ( ζ ξ ) ) L 2 θ 6 ( 1 + τ ) r 2 k + r 2 q 6 σ 2 2 + 1 γ ζ ι ( ζ v ) χ ( λ )
for
θ 5 = r σ 1 1 2 1 q 5 [ 0 , 1 ] , θ 6 = r σ 1 1 2 1 q 6 + σ 2 σ 1 σ 2 σ 1 , 1 .
The last conditions imply that
2 q 3 , q 5 and 2 q 6 if σ 1 r 2 + σ 2 , , 2 q 3 , q 5 and 2 q 6 2 r r 2 ( σ 1 σ 2 ) if σ 1 r 2 , r 2 + σ 2 , 2 q 3 , q 5 2 r r 2 σ 1 and 2 q 6 2 r r 2 ( σ 1 σ 2 ) if σ 1 0 , r 2 .
One possibility to choose the parameters q 3 , q 4 , q 5 and q 6 satisfying the last conditions is
q 3 = r ( ς 1 ) σ 1 σ 2 , q 4 = 2 n ( ς 1 ) r ( ς 1 ) 2 ( σ 1 σ 2 ) , q 5 = r ( ς 1 ) σ 1 σ 2 , q 6 = 2 r r 2 ( σ 1 σ 2 ) .
These choices imply the condition
1 + 2 ( σ 1 σ 2 ) r ς if σ 1 r 2 , r 2 + σ 2 , 1 + 2 ( σ 1 σ 2 ) r ς 1 + 2 ( σ 1 σ 2 ) r 2 σ 1 if σ 1 0 , r 2 .
Consequently, we obtain for 0 τ λ the estimate
| ζ ( τ , · ) | ς | ξ ( τ , z ) | ς H σ 2 ( R r ) 0 1 ( 1 + τ ) r 2 k ς + r 4 σ 2 2 + ( 1 γ ) ς ζ ξ χ ( λ ) ζ ι ( ζ ξ ) χ ( λ ) ς 1 d r 0 1 ( 1 + τ ) r 2 k ς + r 4 σ 2 2 ( 1 γ ) ς ζ ξ χ ( λ ) u χ ( λ ) ς 1 + v χ ( λ ) ς 1 d r ( 1 + τ ) r 2 k ς + r 4 σ 2 2 ( 1 γ ) ς ζ ξ χ ( λ ) u χ ( λ ) ς 1 + v χ ( λ ) ς 1 .
All together similar to the control of N ζ N ζ χ ( λ ) in Theorem 1 leads to
| D | σ ( N ζ N ξ ) ( λ , · ) L 2 ( 1 + λ ) r 2 1 k 1 2 σ 2 ζ ξ χ ( λ ) ζ χ ( λ ) ς 1 + ξ χ ( λ ) ς 1 ,
( N ζ N ζ ) ( λ , · ) L 2 ( R r ) ( 1 + λ ) r 2 1 k 1 2 ζ ξ χ ( λ ) ζ χ ( λ ) ς 1 + ξ χ ( λ ) ς 1 .
The proof is completed.

3.4. Proof of Theorem 4

As we did in the proof to Theorem 3 we need just to modify the estimate of the norm | D | σ 2 ζ λ n L ( λ , · ) L 2 by using the Lemma A1 introduced by D’Abbicco in [27]. Then we get for σ < r 2 the following estimate:
| ζ ( τ , · ) | ς H σ 2 ( R r ) ζ ( τ , · ) H σ 2 ( R r ) ζ ( τ , · ) L ( R r ) ς 1 ζ ( τ , · ) H σ 2 ( R r ) ζ ( τ , · ) H σ ( R r ) + ζ ( τ , · ) H σ 2 ( R r ) ς 1 ζ ( τ , · ) H σ 2 ( R r ) ς + ζ ( τ , · ) H σ 2 ( R r ) ζ ( τ , · ) H σ ( R r ) ς 1 .
Gagliardo-Nirenberg inequality and the definition of the solution space χ ( λ ) lead us to conclusion
| ζ ( τ , · ) | ς H σ 2 ( 1 + B ( τ , 0 ) ) r 2 1 k 1 2 p σ 2 2 σ 2 ( ς 1 ) + ( 1 γ ) ς ζ χ ( λ ) ς .
If we choose σ = r 2 ε , then we get
r 2 1 k 1 2 ς σ 2 2 σ 2 ( ς 1 ) r 2 k ς + r 2 k .
Thus, we obtain
| ζ ( τ , · ) | ς L k H σ 2 ( 1 + τ ) r 2 k ς + r 2 k + ( 1 γ ) ς ζ χ ( λ ) ς .
Using the most recent estimate, the same steps as in the proof of Theorem 3 can be used to conclude the proof.

4. Concluding Remarks

There are a number of Semiliear Cauchy problems in the literature which have the same decay estimates for the homogeneous problem as that of the Cauchy problem (12) with friction, viscoelastic damping. The results that we have acquired also hold for the semilinear Cauchy problems that already exist in mathematical literature. In this study, we have discussed the global existence in time of small data solutions to the Cauchy problem (12) with friction, viscoelastic damping and a fractional nonlinearity, where the data are supposed to belong to different classes of regularity and I γ denote the Caputo fractional integral of order γ defined by I γ ( ζ ) 0 λ ( λ s ) γ ζ ( s ) d s for some γ ( 0 , 1 ) . We have also shown the influence of the fractional nonlinearity to the admissible range of exponent ς comparing with power nonlinearity and also the generating of loss of decay. Indeed the Cauchy problem studied in this paper is more general than the Cauchy problems (13)–(15) since it contains the global existence in time of small data solutions when the data are supposed to belong to different classes of regularity.

Author Contributions

Conceptualization, A.M.D.; Formal analysis, A.M.D.; Writing—original draft, A.M.D.; Writing—review & editing, M.A.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia [Project No. GRANT492], (The Number of the old project that is transferred to the new project is No. AN000441).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no competing interest.

Appendix A

Here we state some inequalities which come into play in our proofs.
Proposition A1.
Let 1 < ς , ς 0 , ς 1 < , ϱ > 0 and s [ 0 , ϱ ) . Then the following fractional Gagliardo-Nirenberg inequality holds for all ζ L ς 0 H ς 1 ϱ :
ζ H ς ϱ ζ L ς 0 1 θ ζ H ς 1 ϱ θ ,
where
θ = θ s , ϱ : = 1 ς 0 1 ς + ϱ r 1 ς 0 1 ς 1 + ϱ r a n d ϱ ϱ θ 1 .
For the proof see [28,29,30,31,32,33,34].
Proposition A2.
Let us assume σ > 0 and 1 m , 1 < ς 1 , ς 2 , q 1 , q 2 satisfying the relation
1 m = 1 ς 1 + 1 ς 2 = 1 q 1 + 1 q 2 .
Then the following fractional Leibniz rule holds:
| D | σ ( f g ) L m | D | σ f L ς 1 g L ς 2 + f L q 1 | D | σ g L q 2 ,
for all f H ς 1 σ L q 1 and g H q 2 σ L ς 2 .
For more details concerning fractional Leibniz rule see [29].
Proposition A3.
Let us choose σ > 0 , ς > σ and 1 < m , m 1 , m 2 < satisfying
1 m = ς 1 m 1 + 1 m 2 .
Let us denote by F ( u ) one of the functions | ζ | ς , ± | ζ | ς 1 ζ . Then the following fractional chain rule holds:
| D | σ F ( u ) L m ζ L m 1 ς 1 | D | σ ζ L m 2 ,
For the proof see [25].
Lemma A1.
Let 0 < 2 σ < r < 2 σ . Then for any function f H σ H σ one has the estimate
f L f H σ + f H σ .
For the proof see [6].
Proposition A4.
The operator N maps χ ( λ ) into itself and has one and only one fixed point ζ χ ( λ ) if the following inequalities hold:
N ζ χ ( λ ) C 0 ( λ ) ( ζ 0 , ζ 1 ) A m , s + C 1 ( λ ) ζ χ ( λ ) ς ,
N ζ N ζ χ ( λ ) C 2 ( λ ) ζ ζ χ ( λ ) ζ χ ( λ ) ς 1 + ζ χ ( λ ) ς 1 ,
where C 1 ( λ ) , C 2 ( λ ) 0 for λ + 0 and C 0 ( λ ) , C 1 ( λ ) , C 2 ( λ ) C for all λ [ 0 , ) .
For the proof see [35].

References

  1. Caputo, M.; Fabrizio, M.A. A new definition of fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar]
  2. Abdeljawad, T.; Baleanu, D. On fractional derivatives with exponential kernel and their discrete versions. Rep. Math. Phys. 2017, 80, 11–27. [Google Scholar] [CrossRef] [Green Version]
  3. Atangana, A.; Baleanu, D. New fractional derivatives with non-local and non-singular kernel, theory and application to heat transfer model. Therm. Sci. 2016, 20, 763–769. [Google Scholar] [CrossRef] [Green Version]
  4. Abdeljawad, T.; Baleanu, D. Integration by parts and its applications of a new nonlocal fractional derivative with Mittag-Leffler nonsingular kernel. J. Nonlinear Sci. Appl. 2017, 10, 1098–1107. [Google Scholar] [CrossRef] [Green Version]
  5. Podlubny, I. Fractional differential equations: An introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications. In Mathematics in Science and Engineering; Elsevier: New York, NY, USA, 1999; Volume 198. [Google Scholar]
  6. D’Abbicco, M. The threshold of effective damping for semilinear wave equations. Math. Meth. Appl. Sci. 2015, 38, 1032–1045. [Google Scholar] [CrossRef] [Green Version]
  7. Cazenave, T.; Dickstein, F.; Weissler, F.B. An equation whose Fujita critical exponent is not given by scaling. Nonlinear Anal. 2008, 68, 862–874. [Google Scholar] [CrossRef]
  8. Lai, N.A.; Zhou, Y. An elementary proof of Strauss conjecture. J. Funct. Anal. 2014, 267, 1364–1381. [Google Scholar] [CrossRef]
  9. Wakasa, K.; Yordanov, B. Blow-up of solutions to critical semilinear wave equations with variable coefficients. J. Differ. Equ. 2019, 266, 5360–5376. [Google Scholar] [CrossRef] [Green Version]
  10. Fino, A.Z.; Kirane, M. Qualitative properties of solutions to a nonlocal evolution system. Math. Methods Appl. Sci. 2011, 34, 1125–1143. [Google Scholar] [CrossRef]
  11. Xu, Y.Q.; Tan, Z. Blow-up of solutions for a time-space fractional evolution system. Acta Math. Sin.-Engl. Ser. 2013, 29, 1067–1074. [Google Scholar] [CrossRef]
  12. Mezadek, M.K.; Mezadek, M.K.; Reissig, M. Semilinear wave models with friction and viscoelastic damping. Math. Meth. Appl. Sci. 2020, 47, 3117–3147. [Google Scholar]
  13. Chen, W.; D’Abbicco, M.; Girardi, G. Global small data solutions for semilinear waves with two dissipative terms. Ann. Mat. 2022, 201, 529–560. [Google Scholar] [CrossRef]
  14. Nakao, M.; Ono, K. Existence of global solutions to the Cauchy problem for the semilinear dissipative wave equations. Math. Z. 1993, 214, 325–342. [Google Scholar] [CrossRef]
  15. Todorova, M.; Yordanov, B. Critical exponent for a nonlinear wave equation with damping. J. Diff. Eq. 2001, 174, 464–489. [Google Scholar] [CrossRef]
  16. Ikehata, R.; Miyaoka, Y.; Nakatake, T. Decay estimates of solutions for dissipative wave equations in R r with lower power nonlinearities. J. Math. Soc. Jpn. 2004, 56, 365–373. [Google Scholar] [CrossRef]
  17. Ikehata, R.; Ohta, M. Critical exponents for semilinear dissipative wave equations in R N . J. Math. Anal. Appl. 2002, 269, 87–97. [Google Scholar] [CrossRef] [Green Version]
  18. Ikehata, R.; Tanizawa, K. Global existence of solutions for semilinear damped wave equations in R r with noncompactly supported initial data. Nonlinear Anal. 2005, 61, 1189–1208. [Google Scholar] [CrossRef]
  19. Shibata, Y. On the rate of decay of solutions to linear viscoelastic equation. Math. Meth. Appl. Sci. 2000, 23, 203–226. [Google Scholar] [CrossRef]
  20. Ikehata, R. Asymptotic profiles for wave equations with strong damping. J. Diff. Equ. 2014, 257, 2159–2177. [Google Scholar] [CrossRef]
  21. D’Abbicco, M.; Reissig, M. Semilinear structural damped waves. Math. Methods Appl. Sci. 2014, 37, 1570–1592. [Google Scholar] [CrossRef] [Green Version]
  22. Ikehata, R.; Sawada, A. Asymptotic profile of solutions for wave equations with frictional and viscoelastic damping terms. Asymptot Anal. 2016, 98, 59–77. [Google Scholar] [CrossRef]
  23. D’Abbicco, M. L1-L1 estimates for a doubly dissipative semilinear wave equation. Nonlinear Differ. Equ. Appl. (NoDEA) 2017, 24, 1–23. [Google Scholar]
  24. Ikehata, R.; Takeda, H. Exponent for nonlinear wave equations with frictional and viscoelastic damping terms. Nonlinear Anal. 2017, 148, 228–253. [Google Scholar] [CrossRef]
  25. Palmieri, A.; Reissig, M. Semi-linear wave models with power non-linearity and scale invariant time-dependent mass and dissipation II. Math. Nachr. 2018, 291, 1859–1892. [Google Scholar] [CrossRef]
  26. D’Abbicco, M.; Girardi, G. A structurally damped σ-evolution equation with nonlinear memory. Math. Meth. Appl. Sci. 2020. [Google Scholar] [CrossRef]
  27. D’Abbicco, M.; Ebert, M.R.; Lucente, S. Self-similar asymptotic profile of the solution to a nonlinear evolution equation with critical dissipation. Math. Meth. Appl. Sci. 2017, 40, 6480. [Google Scholar] [CrossRef]
  28. Christ, F.; Weinstein, M. Dispersion of Small-Amplitude Solutions of the Generalized Korteweg-de Vries Equation. J. Funct. Anal. 1991, 100, 87–109. [Google Scholar] [CrossRef] [Green Version]
  29. Grafakos, L. Classical and Modern Fourier Analysis; Prentice Hall: Upper Saddle River, NJ, USA, 2004. [Google Scholar]
  30. Grafakos, L.; Oh, S. The Kato Ponce inequality. Commun. Partial. Differ. Equ. 2014, 39, 1128–1157. [Google Scholar] [CrossRef]
  31. Gulisashvili, A.; Kon, M. Exact smoothing properties of Schrndinger semigroups. Am. J. Math. 1996, 118, 1215–1248. [Google Scholar]
  32. Hajaiej, H.; Molinet, L.; Ozawa, T.; Wang, B. Necessary and sufficient conditions for the fractional Gagliardo-Nirenberg inequalities and applications to Navier-Stokes and generalized boson equations. In Harmonic Analysis and Nonlinear Partial Differential Equations B26; 2011; pp. 159–175. Available online: https://arxiv.org/pdf/1004.4287v2.pdf (accessed on 30 August 2022).
  33. Kato, T.; Ponce, G. Well-posedness and scattering results for the generalized Kortewegde-Vries equation via the contraction principle. Commun. Pure Appl. Math. 1993, 46, 527–620. [Google Scholar]
  34. Kenig, C.E.; Ponce, G.; Vega, L. Commutator estimates and the Euler and Navier Stokes equations. Commun. Pure Appl. Math. 1988, 41, 891–907. [Google Scholar]
  35. Djaouti, A.M. Semilinear Systems of Weakly Coupled Damped Waves. Ph.D. Thesis, TU Bergakademie Freiberg, Freiberg, Germany, 2018. [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

Mohammed Djaouti, A.; Latif, M.A. Fractional Nonlinearity for the Wave Equation with Friction and Viscoelastic Damping. Axioms 2022, 11, 524. https://doi.org/10.3390/axioms11100524

AMA Style

Mohammed Djaouti A, Latif MA. Fractional Nonlinearity for the Wave Equation with Friction and Viscoelastic Damping. Axioms. 2022; 11(10):524. https://doi.org/10.3390/axioms11100524

Chicago/Turabian Style

Mohammed Djaouti, Abdelhamid, and Muhammad Amer Latif. 2022. "Fractional Nonlinearity for the Wave Equation with Friction and Viscoelastic Damping" Axioms 11, no. 10: 524. https://doi.org/10.3390/axioms11100524

APA Style

Mohammed Djaouti, A., & Latif, M. A. (2022). Fractional Nonlinearity for the Wave Equation with Friction and Viscoelastic Damping. Axioms, 11(10), 524. https://doi.org/10.3390/axioms11100524

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop