Next Article in Journal
Predicting Rice Production in Central Thailand Using the WOFOST Model with ENSO Impact
Next Article in Special Issue
Numerical and Theoretical Stability Study of a Viscoelastic Plate Equation with Nonlinear Frictional Damping Term and a Logarithmic Source Term
Previous Article in Journal
Time to Critical Condition in Emergency Services
Previous Article in Special Issue
Simple Algebraic Expressions for the Prediction and Control of High-Temperature Annealed Structures by Linear Perturbation Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Well-Posedness and Stability Results for a Nonlinear Damped Porous–Elastic System with Infinite Memory and Distributed Delay Terms

1
Department of Mathematics, Faculty of Science, University of Hail, Hail 55425, Saudi Arabia
2
Laboratory of Pure and Applied Mathematics, University of Laghouat, Laghouat 03000, Algeria
3
Department of Mathematics, Faculty of Sciences, University 20 Aout 1955, Skikda 21000, Algeria
4
Department of Mathematics, College of Sciences and Arts, Qassim University, Ar-Rass 51921, Saudi Arabia
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2021, 26(4), 71; https://doi.org/10.3390/mca26040071
Submission received: 30 August 2021 / Revised: 7 October 2021 / Accepted: 13 October 2021 / Published: 16 October 2021

Abstract

:
In the present paper, we consider an important problem from the application perspective in science and engineering, namely, one-dimensional porous–elastic systems with nonlinear damping, infinite memory and distributed delay terms. A new minimal conditions, placed on the nonlinear term and the relationship between the weights of the different damping mechanisms, are used to show the well-posedness of the solution using the semigroup theory. The solution energy has an explicit and optimal decay for the cases of equal and nonequal speeds of wave propagation.

1. Introduction

As introduced in [1], the one-dimensional porous–elastic model constitutes a system of two partial differential equations with unknown ( u , φ ) given by
ρ 0 u t t = μ u x x + β φ x , in 0 , l   ×   0 , L , ρ 0 k φ t t = α φ x x β u x τ φ t ξ φ , in 0 , l   ×   0 , L ,
where l , L > 0 the constant ρ is the mass density, κ is the equilibrated inertia and the constants μ , α , β , τ , ξ are assumed to satisfy the appropriate conditions. This type of problem has been studied by many authors and a lot of results have been shown (please see [1,2,3,4,5,6,7,8,9]). The pioneering contribution was made by [10] for the problem (1). The basic evolution equations for one-dimensional theories of porous materials with memory effect are given by
ρ u t t = T x , J ϕ t t = H x + G ,
where T is the stress tensor, H is the equilibrated stress vector and G is the equilibrated body force. The variables u and ϕ are the displacement of the solid elastic material and the volume fraction, respectively. The constitutive equations are
T = μ u x + b ϕ , H = δ ϕ x 0 t g t s ϕ x s d s , G = b u x ξ ϕ .
A porous–elastic system was considered by [11] in the system
ρ u t t μ u x x b ϕ x = 0 , in 0 , 1 × 0 , , J ϕ t t δ ϕ x x + b u x + ξ ϕ + 0 t g t s ϕ x x x , s d s = 0 , in 0 , 1 × 0 , .
System (4) subjected Neumann–Dirichlet boundary conditions, where g is the relaxation function; the authors obtained a general decay result for the case of equal speeds of wave propagation (See [12,13]). In [14], the authors improved the case of non-equal speed of wave propagation. In [15] the authors considered the following system with memory and distributed delay terms
ρ u t t μ u x x b ϕ x = 0 J ϕ t t δ ϕ x x + b u x + ξ ϕ + 0 t g ( s ) ϕ x x ( t s ) d s + μ 1 ϕ t + τ 1 τ 2 | μ 2 ( ϱ ) | ϕ t ( x , t ϱ ) d ϱ = 0 .
The exponential stability results of systems with memory and distributed delay terms, for the case of equal speeds of wave propagation under a suitable assumptions, are proved. In [16], the following system was considered
ρ u t t μ u x x b ϕ x = 0 , J ϕ t t δ ϕ x x + b u x + a ϕ + 0 g ( s ) ϕ x x ( t s ) d s + α ( t ) f ( ϕ t ) = 0 .
The authors proved the global well-posedness and stability results of (6), which has been extended in [17] for the case of nonequal speeds of wave propagation. Very recently, one-dimensional equations of an homogeneous and isotropic porous–elastic solid with an interior time-dependent delay term feedbacks was treated by Borges Filho and M. Santos in [1].
The result in [10] for system (1) was improved by Apalara to exponential stability in [18]. For more papers related to our paper, please see [19,20,21,22].
Motivated by all the above papers, we investigate the well-posedness and stability results with distributed delay for the cases of equal and nonequal speeds of wave propagation, under additional conditions of the following system
ρ u t t μ u x x b ϕ x = 0 J ϕ t t δ ϕ x x + b u x + ξ ϕ + 0 g ( p ) ϕ x x ( t p ) d p + μ 1 ϕ t + τ 1 τ 2 | μ 2 ( ϱ ) | ϕ t ( x , t ϱ ) d ϱ + α ( t ) f ( ϕ t ) = 0 ,
where
( x , ϱ , t ) ( 0 , 1 ) × ( τ 1 , τ 2 ) × ( 0 , ) ,
with the Neumann–Dirichlet boundary conditions
u x ( 0 , t ) = u x 1 , t = ϕ 0 , t = ϕ 1 , t = 0 , t 0 ,
and the initial data
u x , 0 = u 0 x , u t x , 0 = u 1 x , x ( 0 , 1 ) ϕ x , 0 = ϕ 0 x , ϕ t x , 0 = ϕ 1 x , x ( 0 , 1 ) ϕ t ( x , t ) = f 0 ( x , t ) , ( x , t ) ( 0 , 1 ) × ( 0 , τ 2 ) .
Here, ρ , μ , J , b , δ , ξ and μ 1 are positive constants satisfying μ ξ > b 2 , the term α ( t ) f ( ϕ t ) , where the functions α and f are specified later, represent the nonlinear damping term. The term τ 1 τ 2 | μ 2 ( ϱ ) | ϕ t ( x , t ϱ ) d ϱ is a distributed delay that acts only on the porous equation and τ 1 , τ 2 are two real numbers with 0 τ 1 τ 2 , where μ 2 is an L function, and the function g is called the relaxation function. We first state the following assumptions:
Hypothesis 1 (H1).
g C 1 ( R + , R + ) satisfying
g ( 0 ) > 0 , δ 0 g ( p ) d p = l > 0 , 0 g ( p ) d p = g 0 .
Hypothesis 2 (H2).
There exists a non-increasing differentiable function α , η : R + R + such that
g ( t ) η ( t ) g ( t ) , t 0 ,
and
lim t α ( t ) α ( t ) = 0 .
Hypothesis 3 (H3).
f C 0 ( R , R ) is non-decreasing such that there exist v 1 , v 2 , ε > 0 and a strictly increasing function G C 1 ( [ 0 , ) ) , with G ( 0 ) = 0 and G is a linear or strictly convex C 2 -function on ( 0 , ε ] , such that
s 2 + f 2 ( s ) G 1 ( s f ( s ) ) , | s | < ε v 1 | s | | f ( s ) | v 2 | s | , | s | ε .
which implies that s f ( s ) > 0 for all s 0 . The function f satisfies
| f ( ψ 2 ) f ( ψ 1 ) | k 0 ( | ψ 2 | β + | ψ 1 | β ) | ψ 2 ψ 1 | , ψ 1 , ψ 2 R ,
where k 0 , β > 0 .
Hypothesis 4 (H4).
The bounded function μ 2 : [ τ 1 , τ 2 ] R , satisfying
τ 1 τ 2 | μ 2 ( ϱ ) | d ϱ < μ 1 .
Now, as in [23], taking the following new variable
y ( x , ρ , ϱ , t ) = ϕ t ( x , t ϱ ρ ) ,
then we obtain
ϱ y t ( x , ρ , ϱ , t ) + y ρ ( x , ρ , ϱ , t ) = 0 y ( x , 0 , ϱ , t ) = ϕ t ( x , t ) .
As in [24], we introduce the following new variable
η t ( x , s ) = ϕ ( x , t ) ϕ ( x , t s ) , ( x , t , s ) ( 0 , 1 ) × R + × R + ,
where η t is the relative history of ϕ satisfies
η t t + η s t = ϕ t ( x , t ) , ( x , t , s ) ( 0 , 1 ) × ( 0 , 1 ) × R + × R + .
Consequently, the problem (7) is equivalent to
ρ u t t μ u x x b ϕ x = 0 J ϕ t t l ϕ x x + b u x + ξ ϕ + 0 g ( p ) η x x t ( p ) d p + μ 1 ϕ t + τ 1 τ 2 | μ 2 ( ϱ ) | y ( x , 1 , ϱ , t ) d ϱ + α ( t ) f ( ϕ t ) = 0 ϱ y t ( x , ρ , ϱ , t ) + y ρ ( x , ρ , ϱ , t ) = 0 η t t + η s t = ϕ t ( x , t ) ,
where
( x , ρ , ϱ , t ) ( 0 , 1 ) × ( 0 , 1 ) × ( τ 1 , τ 2 ) × ( 0 , ) ,
with the following boundary and initial conditions
u x 0 , t = u x 1 , t = ϕ 0 , t = ϕ 1 , t = 0 , t 0 ,
u x , 0 = u 0 x , u t x , 0 = u 1 x , x ( 0 , 1 ) ϕ x , 0 = ϕ 0 x , ϕ t x , 0 = ϕ 1 x , x ( 0 , 1 ) y ( x , ρ , ϱ , 0 ) = f 0 ( x , ρ ϱ ) , x ( 0 , 1 ) , ρ ( 0 , 1 ) , ϱ ( 0 , τ 2 ) η t ( x , 0 ) = 0 , η 0 ( x , s ) = η 0 ( x , s ) , ( x , s ) ( 0 , 1 ) × R + .
Meanwhile, from (7)1 and (9), it follows that
d 2 d t 2 0 1 u x , t d x = 0 .
Therefire, by solving (18) and using the initial data of u, we get
0 1 u x , t d x = t 0 1 u 1 x d x + 0 1 u 0 x d x .
Consequently, if we let
u x , t = u x , t t 0 1 u 1 x d x 0 1 u 0 x d x ,
we get
0 1 u x , t d x = 0 , t 0 .
Therefore, the use of Poincaré’s inequality for u is justified. In addition, a simple substitution shows that ( u , ϕ , y , η t ) satisfies system (7). Hence, we work with u instead of u, but write u for simplicity of notation.
By imposing new appropriate conditions (H3), with the help of some special results, we obtain an unusual, weaker decay result using Lyaponov functiona, extending some earlier results known in the existing literature. The main results in this manuscript are as follows: Theorem 1 for the existence and uniqueness of solution and Theorem 2 for the general stability estimates.

2. Well-Posedness

In this section, we prove the existence and uniqueness result of the system (16)–(18) using the semigroup theory. To achieve our goal, we first introduce the vector function
U = ( u , u t , ϕ , ϕ t , y , η t ) T ,
and the new dependent variables v = u t , ψ = ϕ t , φ = η t ; then, the system (16) can be written as follows
U t = A U + Γ ( U ) U ( 0 ) = U 0 = ( u 0 , u 1 , ϕ 0 , ϕ 1 , f 0 , η 0 ) T ,
where A : D ( A ) H : H is the linear operator defined by
A U = v μ ρ u x x + b ρ ϕ x ψ l J ψ x x + b J u x ξ J ϕ x + 1 J 0 g ( p ) φ x x ( p ) d p μ 1 J ψ 1 J τ 1 τ 2 | μ 2 ( ϱ ) | y ( x , 1 , ϱ , t ) d ϱ 1 ϱ y ρ φ s + ψ ,
and
Γ ( U ) = 0 0 0 α ( t ) J f ( ψ ) 0 0 ,
and H is the energy space given by
H = H * 1 ( 0 , 1 ) × L * 2 ( 0 , 1 ) × H 0 1 ( 0 , 1 ) × L 2 ( 0 , 1 ) × L 2 ( ( 0 , 1 ) × ( 0 , 1 ) × ( τ 1 , τ 2 ) ) × L g ( 0 , 1 ) ,
where
L * 2 ( 0 , 1 ) = { Φ L 2 ( 0 , 1 ) / 0 1 Φ ( x ) d x = 0 } , H * 1 ( 0 , 1 ) = H 1 ( 0 , 1 ) L * 2 ( 0 , 1 ) , L g ( 0 , 1 ) = { Φ : R + H 0 1 ( 0 , 1 ) , 0 1 0 g ( s ) Φ x 2 ( p ) d p < } ,
where the space L g ( 0 , 1 ) is endowed with the following inner product
Φ 1 , Φ 2 L g ( 0 , 1 ) = 0 1 0 g ( p ) Φ 1 x ( p ) Φ 2 x ( p ) d p .
For any
U = ( u , v , ϕ , ψ , y , φ ) T H , U ^ = ( u ^ , v ^ , ϕ ^ , ψ ^ , y ^ , φ ^ ) T H .
The space H equipped with the inner product is defined by
U , U ^ H = ρ 0 1 v v ^ d x + μ 0 1 u x u ^ x d x + J 0 1 ψ ψ ^ d x + ξ 0 1 ϕ ϕ ^ d x + l 0 1 ϕ x ϕ ^ x d x + b 0 1 ( u x ϕ ^ + u ^ x ϕ ) d x + 0 1 0 1 τ 1 τ 2 ϱ | μ 2 ( ϱ ) | y y ^ d ϱ d ρ d x + φ , φ ^ L g ( 0 , 1 ) .
The domain of A is given by
D ( A ) = U H / u H * 2 H * 1 , ϕ H 2 H 0 1 , v H * 1 , ψ H 0 1 ( 0 , 1 ) , φ L g ( 0 , 1 ) , y , y ρ L 2 ( ( 0 , 1 ) × ( 0 , 1 ) × ( τ 1 , τ 2 ) ) , y ( x , 0 , ϱ , t ) = ψ ,
where
H * 2 ( 0 , 1 ) = Φ H 2 ( 0 , 1 ) / Φ x ( 1 ) = Φ x ( 0 ) = 0 .
Clearly, D ( A ) is dense in H . Now, we can state and prove the existence result.
Theorem 1.
Let U 0 H and assume that (10)–(15) hold. Then, there exists a unique solution U C ( R + , H ) of problem (20). Moreover, if U 0 D ( A ) , then
U C ( R + , D ( A ) ) C 1 ( R + , H ) .
Proof. 
First, we prove that the operator A is dissipative. For any U 0 D ( A ) and by using (23), we have
A U , U H = μ 1 0 1 ψ 2 d x 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | ψ y ( x , 1 , ϱ , t ) d ϱ d x 0 1 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y ρ y d ϱ d ρ d x 0 1 0 g ( p ) φ x p ( p ) φ x ( p ) d p d x .
For the third term of the RHS of (24), we have
0 1 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y ρ y d ϱ d ρ d x = 1 2 0 1 τ 1 τ 2 0 1 | μ 2 ( ϱ ) | d d ρ y 2 d ρ d ϱ d x = 1 2 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 ( x , 1 , ϱ , t ) d ϱ d x + 1 2 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 ( x , 0 , ϱ , t ) d ϱ d x .
Using Young’s inequality, we obtain
0 1 τ 1 τ 2 ϱ | μ 2 ( ϱ ) | ψ y ( x , 1 , ϱ , t ) d ϱ d x 1 2 ( τ 1 τ 2 | μ 2 ( ϱ ) | d ϱ ) 0 1 ψ 2 d x + 1 2 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 ( x , 1 , ϱ , t ) d ϱ d x .
By integrating the last term of the right-hand side of (24), we have
0 1 0 g ( p ) φ x p ( p ) φ x ( p ) d p d x = 1 2 0 1 0 g ( p ) φ x 2 ( p ) d p d x .
Substituting (25), (26) and (27) into (24), using the fact that y ( x , 0 , ϱ , t ) = ψ ( x , t ) and (15), we obtained
A U , U H μ 1 τ 1 τ 2 | μ 2 ( ϱ ) | d ϱ 0 1 ψ 2 d x + 1 2 0 1 0 g ( p ) φ x 2 ( p ) d p d x 0 .
Hence, the operator A is dissipative.
Next, we prove that the operator A is maximal. This is enough to show that the operator ( λ I A ) is surjective. Indeed, for any F = ( f 1 , f 2 , f 3 , f 4 , f 5 , f 6 ) T H , we prove that there is a unique V = ( u , v , ϕ , ψ , y , φ ) D ( A ) such that
( λ I A ) V = F .
That is
λ u v = f 1 H * 1 ( 0 , 1 ) ρ λ v μ u x x b ϕ x = ρ f 2 L * 2 ( 0 , 1 ) λ ϕ ψ = f 3 H 0 1 ( 0 , 1 ) J λ ψ l ϕ x x + b u x + ξ ϕ 0 g ( p ) φ x x ( p ) d p + μ 1 ψ + τ 1 τ 2 | μ 2 ( ϱ ) | y ( x , 1 , ϱ , t ) d ϱ = J f 4 L 2 ( 0 , 1 ) λ ϱ y t ( x , ρ , ϱ , t ) + y ρ ( x , ρ , ϱ , t ) = ϱ f 5 L 2 ( ( 0 , 1 ) × ( 0 , 1 ) × ( τ 1 , τ 2 ) ) λ φ + φ s ψ = f 6 L g ( 0 , 1 ) .
We note that the equation (30)5 with y ( x , 0 , ϱ , t ) = ψ ( x , t ) has a unique solution, given by
y ( x , ρ , ϱ , t ) = e λ ρ ϱ ψ + ϱ e λ ϱ ρ 0 ρ e λ ϱ σ f 5 ( x , σ , ϱ , t ) d σ ,
then
y ( x , 1 , ϱ , t ) = e λ ϱ ψ + ϱ e λ ϱ 0 1 e λ ϱ σ f 5 ( x , σ , ϱ , t ) d σ ,
and we infer from (30)6 that
φ = e λ s 0 s e τ ( ψ + f 6 ( τ ) ) d τ ,
and we have
v = λ u f 1 , ψ = λ ϕ f 3 .
Inserting (32), (33) and (34) in (30)2 and (30)4, we get
ρ λ 2 u μ u x x b ϕ x = h 1 L * 2 ( 0 , 1 ) μ 3 ϕ μ 4 ϕ x x + b u x = h 2 L 2 ( 0 , 1 ) ,
where
μ 3 = J λ 2 + ξ + λ μ 1 + λ τ 1 τ 2 | μ 2 ( ϱ ) | e λ ϱ d ϱ μ 4 = l + 0 g ( p ) ( 1 e λ p ) d p , h 1 = ρ ( λ f 1 + f 2 ) h 2 = ( J λ + μ 1 + τ 1 τ 2 | μ 2 ( ϱ ) | e λ ϱ d ϱ ) f 3 + J f 4 τ 1 τ 2 ϱ | μ 2 ( ϱ ) | e λ ϱ 0 1 e λ ϱ σ f 5 ( x , σ , ϱ , t ) d σ d ϱ + 0 g ( p ) e λ p 0 p e τ ( ψ + f 6 ( τ ) ) x x d τ d p .
We multiply (35) by u ^ , ϕ ^ , respectively and integrate their sum over ( 0 , 1 ) to obtain the following variational formulation
B ( ( u , ϕ ) , ( u ^ , ϕ ^ ) ) = Υ ( u ^ , ϕ ^ ) ,
where
B : ( H * 1 ( 0 , 1 ) × H 0 1 ( 0 , 1 ) ) 2 R ,
is the bilinear form defined by
B ( ( u , ϕ ) , ( u ^ , ϕ ^ ) ) = λ 2 ρ 0 1 u u ^ d x + μ 3 0 1 ϕ ϕ ^ d x + μ 0 1 u x u ^ x d x + μ 4 0 1 ϕ x ϕ ^ x d x + b 0 1 ( u x ϕ ^ + ϕ u ^ x ) d x ,
and
Υ : ( H * 1 ( 0 , 1 ) × H 0 1 ( 0 , 1 ) ) R ,
is the linear functional given by
Υ ( u ^ , ϕ ^ ) = 0 1 h 1 u ^ d x + 0 1 h 2 ϕ ^ d x
Now, for V = H * 1 ( 0 , 1 ) × H 0 1 ( 0 , 1 ) , equipped with the norm
( u , ϕ ) V 2 = u 2 2 + ϕ 2 2 + u x 2 2 + ϕ x 2 2 ,
we have
B ( ( u , ϕ ) , ( u , ϕ ) ) = λ 2 ρ 0 1 u 2 d x + μ 3 0 1 ϕ 2 d x + μ 0 1 u x 2 d x + 2 b 0 1 u x ϕ d x + μ 4 0 1 ϕ x 2 d x .
On the other hand, we can write
μ u x 2 + 2 b u x ϕ + μ 3 ϕ 2 = 1 2 [ μ u x + b μ ϕ 2 + μ 3 ϕ + b μ 3 u x 2 + μ b 2 μ 3 u x 2 + μ 3 b 2 μ ϕ 2 ] .
Since μ ξ > b 2 , we deduce that
μ u x 2 + 2 b u x ϕ + μ 3 ϕ 2 > 1 2 μ b 2 μ 3 u x 2 + μ 3 b 2 μ ϕ 2 ,
then, for some M 0 > 0
B ( ( u , ϕ ) , ( u , ϕ ) ) M 0 ( u , ϕ ) V 2 .
Thus, B is coercive, similarly,
Υ ( u ^ , ϕ ^ ) M 1 ( u ^ , ϕ ^ ) V 2 .
Consequently, using Lax–Milgram theorem, we conclude that (16) has a unique solution
( u , ϕ ) H * 1 ( 0 , 1 × H 0 1 ( 0 , 1 ) .
Substituting u , ϕ into (32), (33) and (34), respectively, we have
v H * 1 ( 0 , 1 ) , ψ H 0 1 ( 0 , 1 ) , φ L g ( 0 , 1 ) y , y ρ L 2 ( ( 0 , 1 ) × ( 0 , 1 ) × ( τ 1 , τ 2 ) ) .
Moreover, if we take u ^ = 0 H * 1 ( 0 , 1 ) in (37) to obtain
μ 3 0 1 ϕ ϕ ^ d x + b 0 1 u x ϕ ^ d x + μ 4 0 1 ϕ x ϕ ^ x d x = 0 1 h 2 ϕ ^ d ,
we get
μ 4 0 1 ϕ x ϕ ^ x d x = 0 1 ( h 2 μ 3 ϕ b u x ) ϕ ^ d x , ϕ ^ H 0 1 ( 0 , 1 ) ,
which yields
μ 4 ϕ x x = ( h 2 μ 3 ϕ b u x ) L 2 ( 0 , 1 ) .
Thus,
ϕ H 2 ( 0 , 1 ) H 0 1 ( 0 , 1 ) .
Consequently, (45) takes the following form
0 1 ( μ 4 ϕ x x h 2 + μ 3 ϕ + b u x ) ϕ ^ d x = 0 , ϕ ^ H 0 1 ( 0 , 1 ) .
Hence, we get
μ 4 ϕ x x + μ 3 ϕ + b u x = h 2 .
This give (35)2. Similarly, if we take ϕ ^ = 0 H 0 1 ( 0 , 1 ) in (36) to obtain
μ 0 1 u x u ^ x d x + b 0 1 ϕ u ^ x d x + λ 2 ρ 0 1 u u ^ d x = 0 1 h 1 u ^ d x ,
we obtain
μ 0 1 u x u ^ x d x = 0 1 ( h 1 + b ϕ x λ 2 ρ u ) u ^ d x , u ^ H * 1 ( 0 , 1 ) ,
which yields
μ u x x = ( h 1 + b ϕ x λ 2 ρ u ) L * 2 ( 0 , 1 ) .
Consequently, (48) takes the following form
0 1 ( μ u x x h 1 b ϕ x + λ 2 ρ u ) u ^ d x = 0 , u ^ H * 1 ( 0 , 1 ) .
Hence, we obtain
μ u x x b ϕ x + λ 2 ρ u = h 1 .
This give (35)1.
Moreover, (48) also holds for any Φ C 1 ( [ 0 , 1 ] ) . Then, by using integration by parts, we obtain
μ 0 1 u x Φ x d x + 0 1 ( h 1 b ϕ x + λ 2 ρ u ) Φ d x = 0 , Φ C 1 ( [ 0 , 1 ] ) .
Then, we obtain for any Φ C 1 ( [ 0 , 1 ] )
u x ( 1 ) Φ ( 1 ) u x ( 0 ) Φ ( 0 ) = 0 .
Since Φ is arbitrary, we obtain that u x ( 0 ) = u x ( 1 ) = 0 . Hence, u H * 2 ( 0 , 1 ) H * 1 ( 0 , 1 ) . Therefore, the application of regularity theory for the linear elliptic equations guarantees the existence of unique U D ( A ) such that (29) is satisfied. Consequently, we conclude that A is a maximal dissipative operator. Now, we prove that the operator Γ defined in
(22) is locally Lipschitz in H . Let
U = ( u , v , ϕ , ψ , y , φ ) T H , U ^ = ( u ^ , v ^ , ϕ ^ , ψ ^ , y ^ , φ ^ ) T H .
Then, we have
Γ ( U ) Γ ( U ^ ) H M 3 f ( ψ ) f ( ψ ^ ) L 2 ( 01 ) .
By using (14) and Holder and Poincaré’s inequalities, we can obtain
f ( ψ ) f ( ψ ^ ) L 2 ( 01 ) k 0 ( ψ 2 β β + ψ ^ 2 β β ) ψ ψ ^ k 1 ψ x ψ ^ x L 2 ( 01 ) ,
which gives us
Γ ( U ) Γ ( U ^ ) H M 4 U U ^ H .
Then, the operator Γ is locally Lipschitz in H . Consequently, the well-posedness result follows from the Hille–Yosida theorem. The proof is completed. □

3. Stability Result

In this section, we state and prove our decay result for the energy of the system
(16)–(18) using the multiplier technique. We need the following Lemmas.
Lemma 1.
The energy functional E , defined by
E t = 1 2 0 1 ρ u t 2 + μ u x 2 + J ϕ t 2 + l ϕ x 2 + ξ ϕ 2 + 2 b u x ϕ d x + 1 2 0 1 0 g ( p ) φ x 2 ( p ) d p d x + 1 2 0 1 0 1 τ 1 τ 2 ϱ | μ 2 ( ϱ ) | y 2 x , ρ , ϱ , t d ϱ d ρ d x ,
satisfies
E ( t ) η 0 0 1 ϕ t 2 d x + 1 2 0 1 0 g ( p ) φ x 2 ( p ) d p d x + α ( t ) 0 1 ϕ t f ( ϕ t ) d x 0 ,
where η 0 = μ 1 τ 1 τ 2 | μ 2 ( ϱ ) | d ϱ > 0 and φ ( s ) = η t = ϕ ( x , t ) ϕ ( x , t p ) .
Proof. 
Multiplying (16) 1 by u t and (16) 2 by ϕ t , then integration by parts over 0 , 1 and using (17), we get
1 2 d d t 0 1 ρ u t 2 + μ u x 2 + J ϕ t 2 + δ ϕ x 2 + ξ ϕ 2 + 2 b u x ϕ d x 0 1 ϕ x t 0 g ( p ) φ x ( p ) d p d x + μ 1 0 1 ϕ t 2 d x + 0 1 ϕ t τ 1 τ 2 | μ 2 ( ϱ ) | y x , 1 , ϱ , t d ϱ d x + α ( t ) 0 1 ϕ t f ( ϕ t ) d x = 0 .
The last term in the LHS of (54) is estimated as follows
0 1 ϕ t τ 1 τ 2 | μ 2 ( ϱ ) | y x , 1 , ϱ , t d ϱ d x 1 2 τ 1 τ 2 | μ 2 ( ϱ ) | d ϱ 0 1 ϕ t 2 d x + 1 2 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 x , 1 , ϱ , t d ϱ d x ,
and
0 1 ϕ x t 0 g ( p ) φ x p d p d x 1 2 d d t 0 1 0 g ( p ) φ x 2 p d p d x 1 2 0 1 0 g ( p ) φ x 2 p d p d x .
Now, multiplying the equation (16)3 by y | μ 2 ( ϱ ) | and integrating the result over ( 0 , 1 ) × ( 0 , 1 ) × ( τ 1 , τ 2 )
d d t 1 2 0 1 0 1 τ 1 τ 2 ϱ | μ 2 ( ϱ ) | y 2 ( x , ρ , ϱ , t ) d ϱ d ρ d x = 0 1 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y y ρ x , ρ , ϱ , t d ϱ d ρ d x = 1 2 0 1 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | d d ρ y 2 x , ρ , ϱ , t d ϱ d ρ d x = 1 2 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | ( y 2 x , 0 , ϱ , t y 2 ( x , 1 , ϱ , t ) ) d ϱ d x = 1 2 ( τ 1 τ 2 | μ 2 ( ϱ ) | d ϱ ) 0 1 ϕ t 2 d x 1 2 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 x , 1 , ϱ , t d ϱ d x .
Now, using (54), (55), (56) and (57), we have
E ( t ) μ 1 τ 1 τ 2 | μ 2 ( ϱ ) | d ϱ 0 1 ϕ t 2 d x + 1 2 0 1 0 g ( p ) φ x 2 p d p d x α ( t ) 0 1 ϕ t f ( ϕ t ) d x ,
then, by (10), there exists a positive constant η 0 , such that
E ( t ) η 0 0 1 ϕ t 2 d x + 1 2 0 1 0 g ( p ) φ x 2 p d p d x α ( t ) 0 1 ϕ t f ( ϕ t ) d x ,
hence, by (11)–(15) we obtain E is a non-increasing function. □
Remark 1.
Using ( μ ξ > b 2 ) , we conclude that the energy E ( t ) definie by (51) satisfies
E t > 1 2 0 1 ρ u t 2 + μ ^ u x 2 + J ϕ t 2 + l ϕ x 2 + ξ ^ ϕ 2 d x + 1 2 0 1 0 g ( p ) φ x 2 ( p ) d p d x + 1 2 0 1 0 1 τ 1 τ 2 ϱ | μ 2 ( ϱ ) | y 2 x , ρ , ϱ , t d ϱ d ρ d x ,
where
μ ^ = 1 2 ( μ b 2 ξ ) > 0 , ξ ^ = 1 2 ( ξ b 2 μ ) > 0 ,
then E ( t ) is positive function.
Lemma 2.
The functional
D 1 t : = J 0 1 ϕ t ϕ d x + b ρ μ 0 1 ϕ 0 x u t y d y d x + μ 1 2 0 1 ϕ 2 d x ,
satisfies
D 1 t l 2 0 1 ϕ x 2 d x μ ^ 0 1 ϕ 2 d x + ε 1 0 1 u t 2 d x + c ( 1 + 1 ε 1 ) 0 1 ϕ t 2 d x + c 0 1 0 g ( p ) φ x 2 ( p ) d p d x + c 0 1 f 2 ( ϕ t ) d x + c 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 x , 1 , ϱ , t d ϱ d x ,
where μ ^ = ξ b 2 μ > 0 .
Proof. 
Direct computation using integration by parts and Young’s inequality, for ε 1 > 0 , yields
D 1 t = l 0 1 ϕ x 2 d x ξ b 2 μ 0 1 ϕ 2 d x + b ρ μ 0 1 ϕ t 0 x u t y d y d x + 0 1 ϕ x 0 g p φ x p d p d x + α ( t ) 0 1 ϕ f ( ϕ t ) d x + J 0 1 ϕ t 2 d x 0 1 ϕ τ 1 τ 2 | μ 2 ( ϱ ) | y ( x , 1 , ϱ , t ) d ϱ d x l 0 1 ϕ x 2 d x ξ b 2 μ 0 1 ϕ 2 d x + c 1 + 1 ε 1 0 1 ϕ t 2 d x + ε 1 0 1 0 x u t y d y 2 d x + 0 1 ϕ x 0 g p φ x p d p d x 0 1 ϕ τ 1 τ 2 | μ 2 ( ϱ ) | y ( x , 1 , ϱ , t ) d ϱ d x + α ( t ) 0 1 ϕ f ( ϕ t ) d x .
According to Cauchy–Schwartz inequality, it is clear that
0 1 0 x u t y d y 2 d x 0 1 0 1 u t d x 2 d x 0 1 u t 2 d x .
Therefore, estimate (63) becomes
D 1 t δ 0 1 ϕ x 2 d x ξ b 2 μ 0 1 ϕ 2 d x + c 1 + 1 ε 1 0 1 ϕ t 2 d x + ε 1 0 1 u t 2 d x 0 1 ϕ τ 1 τ 2 | μ 2 ( ϱ ) | y ( x , 1 , ϱ , t ) d ϱ d x + 0 1 ϕ x 0 g p φ x p d p d x + α ( t ) 0 1 ϕ f ( ϕ t ) d x .
The last term in the RHS of (64) is estimated as follows
0 1 ϕ x 0 g p φ x p d p d x c δ 1 0 1 ϕ x 2 d x + c 4 δ 1 0 1 0 g p φ x 2 ( p ) d p d x ,
where we used Cauchy–Schwartz, Young and Poincaré’s inequalities, for δ 1 , δ 2 , δ 3 > 0 .
By substituting (65) into (63), we obtain
D 1 t l c δ 1 μ 1 c δ 2 c δ 3 0 1 ϕ x 2 d x ξ b 2 μ 0 1 ϕ 2 d x + ε 1 0 1 u t 2 d x + c ( 1 + 1 ε 1 ) 0 1 ϕ t 2 d x + c 4 δ 1 0 1 0 g p φ x 2 ( p ) d p d x + 1 4 δ 2 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 ( x , 1 , ϱ , t ) d ϱ d x + 1 4 δ 3 0 1 f 2 ( ϕ t ) d x .
Bearing in mind that μ ξ > b 2 and letting δ 1 = l 6 c , δ 2 = l 6 c μ 1 and δ 3 = l 6 c , we obtain estimate (62). □
Lemma 3.
Then, for any ε 2 > 0 the functional
D 2 t : = 0 1 ϕ x u t d x + 0 1 ϕ t u x d x ρ μ J 0 1 u t 0 g p ϕ x t p d p d x ,
satisfies
D 2 t b 2 J 0 1 u x 2 d x + c 0 1 ϕ x 2 d x + c ε 2 0 1 u t 2 d x + c 0 1 ϕ t 2 d x + c 0 1 0 g p φ x 2 p d p d x c ε 2 0 1 0 g p φ x 2 p d p d x + c 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 ( x , 1 , ϱ , t ) d ϱ d x + c 0 1 f 2 ( ϕ t ) d x + δ J μ ρ 0 1 u x ϕ x x d x .
Proof. 
By differentiating D 2 , then using (16), integration by parts and (17) we obtain
D 2 t = b J 0 1 u x 2 d x + l + g 0 J μ ρ 0 1 u x ϕ x x d x + ( b ρ b g 0 μ J ) 0 1 ϕ x 2 d x ξ J 0 1 u x ϕ d x b μ J 0 1 ϕ x 0 g ( p ) φ x p d p d x ρ μ J 0 1 u t 0 g p φ x p d p d x α ( t ) μ J 0 1 u x f ( ϕ t ) d x μ 1 J 0 1 ϕ t u x d x 1 J 0 1 u x τ 1 τ 2 | μ 2 ( ϱ ) | y 2 ( x , 1 , ϱ , t ) d ϱ d x .
In what follows, we estimate the last six terms in the RHS of (68), using Young, Cauchy–Schwartz and Poincaré’s inequalities. For δ 4 , δ 5 , ε 2 > 0 , we have
ξ J 0 1 u x ϕ d x ξ J δ 4 0 1 u x 2 d x + ξ 4 J δ 4 0 1 ϕ 2 d x .
By letting δ 4 = b 6 ξ , using Poincaré’s inequality, we get
ξ J 0 1 u x ϕ d x b 6 J 0 1 u x 2 d x + c 0 1 ϕ x 2 d x ,
and by Young and Chauchy–Schawrz’s inequalities, we get
b μ J 0 1 ϕ x 0 g p φ x p d p d x c δ 5 0 1 ϕ x 2 d x + c 4 δ 5 0 1 0 g p φ x 2 p d p d x .
By letting δ 5 = b 6 c J , we obtain
b μ J 0 1 ϕ x 0 g p φ x p d p d x b 6 J 0 1 ϕ x 2 d x + c 0 1 0 g p φ x 2 p d p d x .
Similarly, ε 2 > 0 we have
ρ μ J 0 1 u t 0 g p φ x p d p d x c ε 2 0 1 u t 2 d x + c ε 2 0 1 0 g p φ x 2 p d p d x ,
and
μ 1 J 0 1 ϕ t u x d x μ 1 δ 6 2 J 0 1 u x 2 d x + μ 1 2 J δ 6 0 1 ϕ t 2 d x ,
and
1 J 0 1 u x τ 1 τ 2 | μ 2 ( ϱ ) | y ( x , 1 , ϱ , t ) d ϱ d x δ 7 μ 1 2 J 0 1 u x 2 d x + 1 2 J δ 7 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 ( x , 1 , ϱ , t ) d ϱ ,
and
α ( t ) J 0 1 u x f ( ϕ t ) d x α ( 0 ) δ 8 2 J 0 1 u x 2 d x + α ( 0 ) 2 J δ 8 0 1 f 2 ( ϕ t ) d x .
Replacing (69)–(74) into (68) and letting δ 6 = δ 7 = b 6 μ 1 and δ 8 = b 6 α ( 0 ) , yields (67). □
Lemma 4.
The functional
D 3 t : = ρ 0 1 u t u d x ,
satisfies
D 3 t ρ 0 1 u t 2 d x + 3 μ 2 0 1 u x 2 d x + c 0 1 ϕ x 2 d x .
Proof. 
Direct computations give
D 3 t = ρ 0 1 u t 2 d x + μ 0 1 u x 2 d x + b 0 1 u x ϕ d x .
The estimate (75) easily follows using Young and Poincaré inequalities.
D 3 t ρ 0 1 u t 2 d x + μ 0 1 u x 2 d x + b ε 0 1 u x 2 d x + b 4 ε 0 1 ϕ 2 d x ρ 0 1 u t 2 d x + μ 0 1 u x 2 d x + b ε 0 1 u x 2 d x + b c 4 ε 0 1 ϕ x 2 d x ,
by taking ε = μ 2 b , we obtain (75). □
Lemma 5.
The functional
D 4 t : = 0 1 0 1 τ 1 τ 2 ϱ e ϱ ρ | μ 2 ( ϱ ) | y 2 x , ρ , ϱ , t d ϱ d ρ d x ,
satisfies
D 4 t η 1 0 1 0 1 τ 1 τ 2 ϱ | μ 2 ( ϱ ) | y 2 x , ρ , ϱ , t d ϱ d ρ d x + μ 1 0 1 ϕ t 2 d x η 1 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 x , 1 , ϱ , t d ϱ d x ,
where η 1 is a positive constant.
Proof. 
By differentiating D 4 , with respect to t, using the Equation (16)3, we have
D 4 t = 2 0 1 0 1 τ 1 τ 2 e ϱ ρ | μ 2 ( ϱ ) | y y ρ x , ρ , ϱ , t d ϱ d ρ d x = 0 1 0 1 τ 1 τ 2 ϱ e ϱ ρ | μ 2 ( ϱ ) | y 2 x , ρ , ϱ , t d ϱ d ρ d x 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | [ e ϱ y 2 x , 1 , ϱ , t y 2 x , 0 , ϱ , t ] d ϱ d x .
Using the fact that y ( x , 0 , ϱ , t ) = ϕ t ( x , t ) and e ϱ e ϱ ρ 1 , for all 0 < ρ < 1 , we obtain
D 4 t = η 1 0 1 0 1 τ 1 τ 2 ϱ | μ 2 ( ϱ ) | y 2 x , ρ , ϱ , t d ϱ d ρ d x 0 1 τ 1 τ 2 e ϱ | μ 2 ( ϱ ) | y 2 x , 1 , ϱ , t d ϱ d x + τ 1 τ 2 | μ 2 ( ϱ ) | d ϱ 0 1 ϕ t 2 d x .
Since e ϱ is an increasing function, we have e ϱ e τ 2 , for all ϱ [ τ 1 , τ 2 ] .
Finally, setting η 1 = e τ 2 and recalling (15), we obtain (76). We are now ready to prove the main result.
Theorem 2.
Assume (10)–(15) hold. Let h ( t ) = α ( t ) . η ( t ) be a positive non-increasing function. Then, for any U 0 D ( A ) , satisfying, for some c 0 > 0
m a x 0 1 ϕ 0 x 2 ( x , s ) d x , 0 1 ϕ 0 s x 2 ( x , s ) d x c 0 , s > 0 ,
there are positive constants β 1 , β 2 and β 3 such that the energy functional given by (52) satisfies
E t β 1 G 0 1 β 2 + β 3 0 t h ( p ) ϖ ( p ) d p 0 t h ( p ) d p ,
where
G 0 ( t ) = t G ( ε 0 t ) , ε 0 0 , a n d ϖ ( s ) = s g ( σ ) d σ .
Proof. 
We define a Lyapunov functional
L t : = N E t + N 1 D 1 t + N 2 D 2 t + D 3 t + N 4 D 4 t ,
where N , N 1 , N 2 , and N 4 are positive constants, to be chosen later. By differentiating (80) and using (53), (62), (67), (75), (76), we have
L t l N 1 2 c N 2 c 0 1 ϕ x 2 d x ρ N 1 ε 1 N 2 c ε 2 0 1 u t 2 d x b N 2 2 J 3 μ 2 0 1 u x 2 d x + c N 1 + N 2 0 1 0 g ( p ) φ x 2 ( p ) d p d x η 0 N c N 1 ( 1 + 1 ε 1 ) N 2 c μ 1 N 4 0 1 ϕ t 2 d x N 4 η 1 c N 1 c N 2 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 x , 1 , ϱ , t d ϱ d x N 1 μ ^ 0 1 ϕ 2 d x + N 2 c N 2 ε 2 0 1 0 g ( p ) φ x 2 ( p ) d p d x N 4 η 1 0 1 0 1 τ 1 τ 2 ϱ | μ 2 ( ϱ ) | y 2 x , ρ , ϱ , t d ϱ d ρ d x + c N 1 + N 2 0 1 f 2 ( ϕ t ) d x + N 2 χ 0 1 u x x ϕ x d x ,
where χ = ( μ ρ δ J ) and by setting
ε 1 = ρ 4 N 1 , ε 2 = ρ 4 c N 2 ,
we obtain
L t l N 1 2 c N 2 ( 1 + N 2 ) c 0 1 ϕ x 2 d x ρ 2 0 1 u t 2 d x b N 2 2 J 3 μ 2 0 1 u x 2 d x + c N 1 + N 2 0 1 0 g ( p ) φ x 2 ( p ) d p d x η 0 N c N 1 ( 1 + N 1 ) c N 2 μ 1 N 4 0 1 ϕ t 2 d x N 4 η 1 c N 1 c N 2 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 x , 1 , ϱ , t d ϱ d x N 1 μ ^ 0 1 ϕ 2 d x + N 2 c N 2 2 0 1 0 g ( p ) φ x 2 ( p ) d p d x N 4 η 1 0 1 0 1 τ 1 τ 2 ϱ | μ 2 ( s ) | y 2 x , ρ , ϱ , t d ϱ d ρ d x + c N 1 + N 2 0 1 f 2 ( ϕ t ) d x + N 2 χ 0 1 u x x ϕ x d x .
Next, we carefully choose our constants so that the terms inside the brackets are positive. We choose a N 2 that is large enough that
α 1 = b N 2 2 J 3 μ 2 > 0 ,
then, we choose a large enough N 1 that
α 2 = l N 1 4 c N 2 1 + N 2 c > 0 ,
then we choose a large enough N 4 that
α 3 = N 4 η 1 c N 1 c N 2 > 0 ,
thus, we arrive at
L t α 2 0 1 ϕ x 2 d x α 0 0 1 ϕ 2 d x ρ 2 0 1 u t 2 d x α 1 0 1 u x 2 d x [ η 0 N c ] 0 1 ϕ t 2 d x + N 2 c 0 1 0 g ( p ) φ x 2 ( p ) d p d x + c 0 1 0 g ( p ) φ x 2 ( p ) d p d x α 3 0 1 τ 1 τ 2 | μ 2 ( ϱ ) | y 2 x , 1 , ϱ , t d ϱ d x α 4 0 1 0 1 τ 1 τ 2 ϱ | μ 2 ( ϱ ) | y 2 x , ρ , ϱ , t d ϱ d ρ d x + c 0 1 f 2 ( ϕ t ) d x + α 5 0 1 u x x ϕ x d x .
where α 0 = μ ^ N 1 = ξ b 2 μ N 1 , a n d α 5 = N 2 χ = N 2 ( μ ρ δ J ) . On the other hand, if we let
L t = N 1 D 1 t + N 2 D 2 t + D 3 t + N 4 D 4 t ,
then
L t J N 1 0 1 ϕ ϕ t d x + b ρ N 1 μ 0 1 ϕ 0 x u t y d y d x + N 2 0 1 ϕ x u t + u x ϕ t ρ μ J u t 0 g p ϕ x t p d p d x + ρ 0 1 u t u d x + N 4 0 1 0 1 τ 1 τ 2 ϱ e ϱ ρ | μ 2 ( ϱ ) | y 2 x , ρ , ϱ , t d ϱ d ρ d x .
Exploiting Young, Cauchy–Schwartz and Poincaré inequalities, we obtain
L t c 0 1 u t 2 + ϕ t 2 + ϕ x 2 + u x 2 + ϕ 2 d x + c 0 1 0 g ( p ) φ x 2 ( p ) d p d x + c 0 1 0 1 τ 1 τ 2 ϱ | μ 2 ( s ) | y 2 x , ρ , ϱ , t d ϱ d ρ c E t .
Consequently, we obtain
L t = L t N E t c E t ,
that is
N c E t L t N + c E t .
Now, by choosing a large enough N that
N 2 c > 0 , N c > 0 , N η 0 c > 0 ,
and exploiting (52), estimates (81) and (82), respectively, we obtain
c 2 E t L t c 3 E t , t 0 ,
and
L t k 1 E t + k 2 0 1 0 g ( p ) φ x 2 ( p ) d p d x + k 3 0 1 ( ϕ t 2 + f 2 ( ϕ t ) ) d x + α 5 0 1 u x x ϕ x d x ,
for some k 1 , k 2 , k 3 , c 2 , c 3 > 0 .
Case 1.
If χ = ( μ ρ δ J ) = 0 , in this case, (84) takes the form
L t k 1 E t + k 2 0 1 0 g ( p ) φ x 2 ( p ) d p d x + k 3 0 1 ( ϕ t 2 + f 2 ( ϕ t ) ) d x .
By multiplying (85) by h ( t ) = α ( t ) . η ( t ) , we obtain
h ( t ) L t k 1 h ( t ) E t + k 2 h ( t ) 0 1 0 g ( p ) φ x 2 ( p ) d p d x + k 3 h ( t ) 0 1 ( ϕ t 2 + f 2 ( ϕ t ) ) d x .
We distinguish two cases
  • G is linear on [ 0 , ε ] . In this case, using the assumption (13)1 and (53), we can write
    k 3 h ( t ) 0 1 ( ϕ t 2 + f 2 ( ϕ t ) ) d x k 3 h ( t ) 0 1 ϕ t f ( ϕ t ) ) d x k 3 η ( t ) E ( t ) ,
    and by (11) we have
    h ( t ) 0 1 0 t g ( p ) φ x 2 ( p ) d p d x = α ( t ) 0 1 0 t η ( s ) g ( p ) φ x 2 ( p ) d p d x α ( t ) 0 1 0 t g ( p ) φ x 2 ( p ) d p d x α ( t ) 0 1 0 g ( p ) φ x 2 ( p ) d p d x 2 α ( t ) E t ,
    and by (77) we obtain
    0 1 φ x 2 ( s ) d x = 2 0 1 ϕ x 2 ( x , t ) d x + 2 0 1 ϕ x 2 ( x , t s ) d x 4 sup s > 0 0 1 ϕ x 2 ( x , s ) d x + 2 sup τ > 0 0 1 ϕ 0 x 2 ( x , τ ) d x 8 E ( 0 ) l + 2 c 0 ,
    then, we obtain
    h ( t ) 0 1 t g ( p ) φ x 2 ( p ) d p d x ( 8 E ( 0 ) l + 2 c 0 ) h ( t ) t g ( p ) d p .
    Hence,
    h ( t ) 0 1 0 g ( p ) φ x 2 ( p ) d p d x 2 α ( t ) E t + ( 8 E ( 0 ) l + 2 c 0 ) h ( t ) ϖ ( t ) .
    Inserting (87) and (91) in (86). Since h ( t ) 0 , α ( t ) 0 , η ( t ) 0 . Then, we have
    L 1 t k 1 h ( t ) E t + γ h ( t ) ϖ ( t ) ,
    and
    m 1 E ( t ) L 1 t m 2 E ( t ) ,
    with
    m 1 = τ 1 , m 2 = c 2 h ( 0 ) + k 3 η ( 0 ) + 2 k 2 α ( 0 ) + τ 1 ,
    where
    L 1 ( t ) = h ( t ) L ( t ) + ( k 3 η ( t ) + 2 k 2 α ( t ) + τ 1 ) E ( t ) E ( t ) , γ = ( 8 E ( 0 ) l + 2 c 0 ) , τ 1 > 0 a n d ϖ ( t ) = t g ( p ) d p .
    Since E ( t ) 0 , t 0 . By using (92), we have
    E ( T ) 0 T h ( t ) d t L 1 ( 0 ) k 1 + γ k 1 0 T h ( t ) ϖ ( t ) d t .
    Using the fact that G 0 1 is linear. Then,
    E ( T ) ζ G 0 1 L 1 ( 0 ) k 1 + γ k 1 0 T h ( t ) ϖ ( t ) d t 0 T h ( t ) d t .
    with β 1 = ζ , β 2 = L 1 ( 0 ) k 1 , β 3 = γ k 1 . This completes the proof.
  • G is nonlinear on [ 0 , ε ] , we choose 0 ε 1 ε and we consider
    I 1 ( t ) = { x ( 0 , 1 ) , | ϕ t | ε 1 } , I 2 = { x ( 0 , 1 ) , | ϕ t | > ε 1 } ,
    we define
    I = I 1 ϕ t f ( ϕ t ) d t .
    Using Jensen’s inequality and the assumption (13)1, we have
    k 3 h ( t ) 0 1 ( ϕ t 2 + f 2 ( ϕ t ) ) d x k 3 h ( t ) 0 1 ϕ t f ( ϕ t ) ) d x k 3 h ( t ) G 1 ( I ( t ) ) k 3 η ( t ) E ( t ) .
    Inserting (97) in (86), since α ( t ) 0 , η ( t ) 0 and E ( t ) 0 , we obtain
    L 2 t k 1 h ( t ) E t + γ h ( t ) ϖ ( t ) + k 3 h ( t ) G 1 ( I ( t ) ) .
    and
    m 3 E ( t ) L 2 t m 4 E ( t ) ,
    with
    m 3 = τ 1 , m 4 = c 2 h ( 0 ) + k 3 η ( 0 ) + 2 k 2 α ( 0 ) + τ 1 ,
    where
    L 2 ( t ) = h ( t ) L ( t ) + ( k 3 η ( t ) + 2 k 2 α ( t ) + τ 1 ) E ( t ) E ( t ) .
    Now, for ε 0 < ε 1 and by using E ( t ) 0 , G > 0 and G > 0 on ( 0 , ε ] , we define the functional L 3 ( t ) by,
    L 3 ( t ) = G ( ε 0 E ( t ) ) L 2 ( t ) + τ 2 E ( t ) E ( t ) , τ 2 > 0 ,
    satisfies
    L 3 ( t ) = E ( t ) ( ε 0 G ( ε 0 E ( t ) ) L 2 ( t ) + τ 2 ) + L 2 ( t ) G ( ε 0 E ( t ) ) k 1 h ( t ) G 0 ( E ( t ) ) + γ G ( ε 0 E ( t ) ) h ( t ) ϖ ( t ) + k 3 h ( t ) G ( ε 0 E ( t ) ) G 1 ( I ( t ) ) .
    To estimate the last term of (92), using the general Young’s inequality
    A B G * ( A ) + G ( B ) , i f A ( 0 , G ( ε ) ) , B ( 0 , ε ) ,
    where
    G * ( A ) = s ( G ) 1 ( s ) G ( ( G ) 1 ( s ) ) , i f s ( 0 , G ( ε ) ) ,
    satisfies
    k 3 h ( t ) G ( ε 0 E ( t ) ) G 1 ( I ( t ) ) k 3 ε 0 h ( t ) G 0 ( E ( t ) ) k 3 η ( t ) E ( t ) .
    Inserting (102) in (92) and letting ε 0 = k 1 2 k 3 we get
    L 3 ( t ) + k 3 η ( t ) E ( t ) k 1 h ( t ) G 0 ( E ( t ) ) + γ G ( ε 0 E ( t ) ) h ( t ) ϖ ( t ) .
    Since η ( t ) 0 , then
    L 4 ( t ) k 1 h ( t ) G 0 ( E ( t ) ) + γ G ( ε 0 E ( t ) ) h ( t ) ϖ ( t ) ,
    where
    L 4 ( t ) = L 3 ( t ) + k 3 η ( t ) E ( t ) E ( t ) .
    Since α ( t ) , G 0 ( E ( t ) ) , G ( ε 0 E ( t ) ) are non-increasing functions,
    then, for any T > 0
    k 1 G 0 ( E ( T ) ) 0 T h ( t ) d t k 1 0 T h ( t ) G 0 ( E ( t ) ) d t L 4 ( 0 ) + γ G ( ε 0 E ( 0 ) ) 0 T h ( t ) ϖ ( t ) d t ,
    which gives (78) with β 1 = 1 , β 2 = L 4 ( 0 ) k 1 and β 3 = γ G ( ε 0 E ( 0 ) ) k 1 .
    The proof is now completed.
Case 2.
If χ = ( μ ρ δ J ) 0 and
| χ | < k 1 μ 2 l 2 N 2 ( l ρ + b μ ) i f χ < 0 | χ | < k 1 μ 2 2 N 2 ρ i f χ > 0 .
This case is more important from the physical perspective, where waves are not necessarily of equal speeds. Let
E ( t ) = E ( u , ϕ , y , φ ) = E 1 ( t ) .
Denotes the first-order energy defined in (52) and
E 2 ( t ) = E ( u t , ϕ t , y t , φ t ) .
Denotes the second-order energy; then, we have
E 2 ( t ) η 0 0 1 ϕ t t 2 d x + 1 2 0 1 0 g ( p ) φ t x 2 ( p ) d p α ( t ) 0 1 ϕ t t f ( ϕ t ) d x α ( t ) 0 1 ϕ t t 2 f ( ϕ t ) d x = η 0 0 1 ϕ t t 2 d x + 1 2 0 1 0 g ( p ) φ t x 2 ( p ) d p + α ( t ) α ( t ) α ( t ) 0 1 ϕ t t f ( ϕ t ) d x 0 1 ϕ t t 2 f ( ϕ t ) d x .
Since f , g are non-decreasing functions, α ( t ) is a positive function and lim t α ( t ) α ( t ) = 0 , we deduce that
E 2 ( t ) η 0 0 1 ϕ t t 2 d x + 1 2 0 1 0 g ( s ) φ t x 2 η 0 0 1 ϕ t t 2 d x ,
where η 0 = μ 1 τ 1 τ 2 | μ 2 ( ϱ ) | d ϱ > 0 .
The last term in (84), by using (16)1, Young’s inequality and by setting K = χ N 2 ρ μ = α 5 ρ μ and α 5 = χ N 2 as follows
α 5 0 1 u x x ϕ x d x = α 5 ρ μ 0 1 ϕ x u t t d x b α 5 μ 0 1 ϕ x 2 d x = K d d t 0 1 ϕ t u x d x + 0 1 ϕ x u t d x K 0 1 u x ϕ t t 2 d x b α 5 μ 0 1 ϕ x 2 d x K d d t 0 1 ϕ t u x d x + 0 1 ϕ x u t d x + | K | 4 0 1 ϕ t t 2 d x + | K | 0 1 u x 2 d x .
Let
N ( t ) = 0 1 ϕ t u x d x + 0 1 ϕ x u t d x ,
then (84)
L ( t ) + K N ( t ) k 1 E 1 ( t ) + k 2 0 1 0 g ( p ) φ p 2 d p d x + | K | 4 0 1 ϕ t t 2 d x + | K | 0 1 u x 2 d x + k 3 0 1 ( ϕ t 2 + f 2 ( ϕ t ) ) d x k 4 E 1 ( t ) + k 2 0 1 0 g ( p ) φ p 2 d p d x + | K | 4 0 1 ϕ t t 2 d x + k 3 0 1 ( ϕ t 2 + f 2 ( ϕ t ) ) d x ,
where
k 4 = k 1 2 | K | μ > 0 .
Let
R ( t ) = L ( t ) + K N ( t ) + N 5 ( E 1 ( t ) + E 2 ( t ) ) .
Indeed, by using Young’s inequality, we obtain
| N ( t ) | = | 0 1 ϕ u x t d x | + | 0 1 ϕ t u x d x | 1 2 0 1 u t 2 d x + 1 2 0 1 ϕ t 2 d x + 1 2 0 1 ϕ x 2 d x + 1 2 0 1 u x 2 d x C 0 E 1 ( t ) ,
where C 0 = m a x { 1 J , 1 ξ , 1 ρ , 1 μ } .
By (83) and (109), we obtain
| R ( t ) N 5 ( E 1 ( t ) + E 2 ( t ) ) | ( c 3 + C 0 ) E 1 ( t ) c ( E 1 ( t ) + E 2 ( t ) ) ,
and
( N 5 c ) ( E 1 ( t ) + E 2 ( t ) ) R ( t ) ( N 5 + c ) ( E 1 ( t ) + E 2 ( t ) ) ,
and by using (105), (107) and (3), we obtain
R ( t ) = L ( t ) + K N ( t ) + N 5 ( E 1 ( t ) + E 2 ( t ) ) k 4 E 1 ( t ) + k 2 0 1 0 g ( p ) φ p 2 d p d x + k 3 0 1 ( ϕ t 2 + f 2 ( ϕ t ) ) d x ( η 0 N 5 | K | 4 ) 0 1 ϕ t t 2 d x .
We choose a large enough N 5 that
η 0 N 5 | K | 4 > 0 , N 5 c > 0 ,
we obtain
R ( t ) ( E 1 ( t ) + E 2 ( t ) ) ,
and
R ( t ) k 4 E 1 ( t ) + k 2 0 1 0 g ( p ) φ p 2 d p d x + k 3 0 1 ( ϕ t 2 + f 2 ( ϕ t ) ) d x .
By multiplying (114) by h ( t ) = α ( t ) . η ( t ) , we obtain
h ( t ) R t k 4 h ( t ) E t + k 2 h ( t ) 0 1 0 g ( p ) φ x 2 ( p ) d p d x + k 3 h ( t ) 0 1 ( ϕ t 2 + f 2 ( ϕ t ) ) d x .
We distinguish two cases
  • G is linear on [ 0 , ε ] . In the same way as in the previous case, we obtain
    R 1 t k 4 h ( t ) E t + γ h ( t ) ϖ ( t ) ,
    and
    m 1 ( E 1 ( t ) + E 2 ( t ) ) R 1 t m 2 ( E 1 ( t ) + E 2 ( t ) ) ,
    with
    m 1 = τ 1 , m 2 = c 2 h ( 0 ) + k 3 η ( 0 ) + 2 k 2 α ( 0 ) + τ 1 ,
    where
    R 1 ( t ) = h ( t ) R ( t ) + ( k 3 η ( t ) + 2 k 2 α ( t ) + τ 1 ) E ( t ) ( E 1 ( t ) + E 2 ( t ) ) γ = ( 8 E ( 0 ) l + 2 c 0 ) , τ 1 > 0 a n d ϖ ( t ) = t g ( p ) d p .
    Since E ( t ) 0 , t 0 . By using (116), we have
    E ( T ) 0 T h ( t ) d t R 1 ( 0 ) k 4 + γ k 4 0 T h ( t ) ϖ ( t ) d t .
    Using the fact that G 0 1 is linear. Then,
    E ( T ) ζ G 0 1 R 1 ( 0 ) k 4 + γ k 4 0 T h ( t ) ϖ ( t ) d t 0 T h ( t ) d t ,
    with β 1 = ζ , β 2 = R 1 ( 0 ) k 4 , β 3 = γ k 4 . This completes the proof.
  • G is nonlinear on [ 0 , ε ] , we choose 0 ε 1 ε . In a similar way to that in the previous case, we have
    R 2 t k 1 h ( t ) E t + γ h ( t ) ϖ ( t ) + k 3 h ( t ) G 1 ( I ( t ) ) ,
    and
    m 3 ( E 1 ( t ) + E 2 ( t ) ) R 2 t m 4 ( E 1 ( t ) + E 2 ( t ) ) ,
    with
    m 3 = τ 1 , m 4 = c 2 h ( 0 ) + k 3 η ( 0 ) + 2 k 2 α ( 0 ) + τ 1 ,
    where
    R 2 ( t ) = h ( t ) R ( t ) + ( k 3 η ( t ) + 2 k 2 α ( t ) + τ 1 ) E ( t ) ( E 1 ( t ) + E 2 ( t ) ) .
    Now, for ε 0 < ε 1 , and by using E ( t ) 0 , G > 0 and G > 0 on ( 0 , ε ] , we define the functional L 3 ( t ) by,
    R 3 ( t ) = G ( ε 0 E ( t ) ) R 2 ( t ) + τ 2 E ( t ) ( E 1 ( t ) + E 2 ( t ) ) , τ 2 > 0 ,
    satisfies
    R 3 ( t ) = E ( t ) ( ε 0 G ( ε 0 E ( t ) ) R 2 ( t ) + τ 2 ) + R 2 ( t ) G ( ε 0 E ( t ) ) k 4 h ( t ) G 0 ( E ( t ) ) + γ G ( ε 0 E ( t ) ) h ( t ) ϖ ( t ) + k 3 h ( t ) G ( ε 0 E ( t ) ) G 1 ( I ( t ) ) .
    To estimate the last term of (122), again using the general Young’s inequality (101).
    Inserting (122) in (121) and letting ε 0 = k 1 2 k 3 , we get
    R 3 ( t ) + k 3 η ( t ) E ( t ) k 4 h ( t ) G 0 ( E ( t ) ) + γ G ( ε 0 E ( t ) ) h ( t ) ϖ ( t ) .
    Since η ( t ) 0 , then
    R 4 ( t ) k 4 h ( t ) G 0 ( E ( t ) ) + γ G ( ε 0 E ( t ) ) h ( t ) ϖ ( t ) ,
    where
    R 4 ( t ) = R 3 ( t ) + k 3 η ( t ) E ( t ) ( E 1 ( t ) + E 2 ( t ) ) .
    Since α ( t ) , G 0 ( E ( t ) ) , G ( ε 0 E ( t ) ) are non-increasing functions, then, for any T > 0
    k 4 G 0 ( E ( T ) ) 0 T h ( t ) d t k 4 0 T h ( t ) G 0 ( E ( t ) ) d t R 4 ( 0 ) + γ G ( ε 0 E ( 0 ) ) 0 T h ( t ) ϖ ( t ) d t ,
    which gives (78) with β 1 = 1 , β 2 = R 4 ( 0 ) k 4 and β 3 = γ G ( ε 0 E ( 0 ) ) k 4 .
    The proof is completed.

Author Contributions

Writing—original draft preparation, D.O.; writing—review and editing, K.B.; visualization, A.M.; supervision, Y.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank the anonymous referees and the handling editor for their careful reading and for relevant remarks/suggestions to improve the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Santos, M.L.; Junior, D.A. On porous-elastic system with localized damping. Z. Angew. Math. Phys. 2016, 67, 63. [Google Scholar] [CrossRef]
  2. Bahri, N.; Beniani, A.; Zennir, K.; Zhang, H.W. Existence and exponential stability of solutions for laminated viscoelastic Timoshenko beams. Appl. Sci. 2020, 22, 1–16. [Google Scholar]
  3. Bouzettouta, L.; Zennir, K.; Zitouni, S. Uniform Decay for a Viscoelastic Wave Equation with Density and Time-Varying Delay in Rn. Filomat 2019, 33, 961–970. [Google Scholar]
  4. Djebabla, A.; Choucha, A.; Ouchenane, D.; Zennir, K. Explicit Stability for a Porous Thermoelastic System with Second Sound and Distributed Delay. Int. J. Appl. Comput. Math. 2021, 7, 50. [Google Scholar] [CrossRef]
  5. Dridi, H.; Zennir, K. Well-posedness and energy decay for some thermoelastic systems of Timoshenko type with Kelvin–Voigt damping. SeMA 2021, 78, 385–400. [Google Scholar] [CrossRef]
  6. Ouchenane, D.; Zennir, K.; Guidad, D. Well-posedness and a general decay for a nonlinear damped porous thermoelastic system with second sound and distributed delay terms. J. Appl. Nonlinear Dyn. in press.
  7. Ouchenane, D.; Zennir, K. General decay of solutions in one-dimensional porous-elastic system with memory and distributed delay term with second sound. Commun. Optim. Theory 2020, 2020, 18. [Google Scholar]
  8. Rivera, J.E.M.; Pamplona, P.X.; Quintanilla, R. On the decay of solutions for porous-elastic system with history. J. Math. Anal. Appl. 2011, 379, 682–705. [Google Scholar]
  9. Zennir, K. Stabilization for Solutions of Plate Equation with Time-Varying Delay and Weak-Viscoelasticity in Rn. Russ. Math. 2020, 64, 21–33. [Google Scholar] [CrossRef]
  10. Quintanilla, R. Slow decay in one-dimensional porous dissipation elasticity. Appl. Math. Lett. 2003, 16, 487–491. [Google Scholar] [CrossRef] [Green Version]
  11. Apalara, T.A. General decay of solution in one-dimensional porous-elastic system with memory. J. Math. Anal. Appl. 2017, 469, 457–471. [Google Scholar] [CrossRef]
  12. Zennir, K.; Ouchenane, D.; Choucha, A.; Biomy, M. Well-posedness and stability for Bresse–Timoshenko type systems with thermodiffusion effects and nonlinear damping. AIMS Math. 2021, 6, 2704–2721. [Google Scholar] [CrossRef]
  13. Hebhoub, F.; Zennir, K.; Miyasita, T.; Biomy, M. Blow up at well defined time for a coupled system of one spatial variable Emden–Fowler type in viscoelasticities with strong nonlinear sources. AIMS Math. 2021, 6, 442–455. [Google Scholar] [CrossRef]
  14. Feng, B. On the decay rates for a one-dimensional porous elasticity system with past history. Commun. Pure Appl. Anal. 2019, 19, 2905–2921. [Google Scholar] [CrossRef] [Green Version]
  15. Choucha, A.; Ouchenane, D.; Zennir, K. General Decay of Solutions in One-Dimensional Porous-Elastic with Memory and Distributed Delay Term. Tamkang J. Math. 2021, 52, 1–18. [Google Scholar] [CrossRef]
  16. Khochemane, H.; Bouzettouta, L.; Zitouni, S. General decay of a nonlinear damping porous-elastic system with past history. Ann. Univ. Ferrara 2019, 65, 249–275. [Google Scholar] [CrossRef]
  17. Khochemane, H.; Djebabla, A.; Zitouni, S.; Bouzettouta, L. Well-posedness and general decay of a nonlinear damping porous-elastic system with infinite memory. J. Math. Phys. 2020, 61, 021505. [Google Scholar] [CrossRef]
  18. Apalara, T.A. Exponential decay in one-dimensional porous dissipation elasticity. Q. J. Mech. Appl. Math. 2017, 70, 363–372, Correction in 2017, 70, 553–555. [Google Scholar] [CrossRef]
  19. Apalara, T.A.; Soufyane, A. Energy decay for a weakly nonlinear damped porous system with a nonlinear delay. Appl. Anal. 2021. [Google Scholar] [CrossRef]
  20. Apalara, T.A. A general decay for a weakly nonlinearly damped porous system. J. Dyn. Control Syst. 2019, 25, 311–322. [Google Scholar] [CrossRef]
  21. Dos Santos, M.J.; Feng, B.; Junior, D.S.A.; Santos, M.L. Global and exponential attractors for a nonlinear porous elastic system with delay term. Discrete Contin. Dyn. Syst. B 2021, 26, 2805. [Google Scholar]
  22. Freitas, M.M.; Santos, M.L.; Langa, J.A. Porous elastic system with nonlinear damping and sources terms. J. Diff. Equ. 2018, 264, 2970–3051. [Google Scholar] [CrossRef]
  23. Nicaise, A.S.; Pignotti, C. Stabilization of the wave equation with boundary or internal distributed delay. Differ. Integral Equ. 2008, 21, 935–958. [Google Scholar]
  24. Fernandez Sare, H.D.; Munnoz Rivera, J.E. Stability of Timoshenko systems with past history. J. Math. Anal. Appl. 2008, 339, 482–502. [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

Moumen, A.; Ouchenane, D.; Bouhali, K.; Altayeb, Y. Well-Posedness and Stability Results for a Nonlinear Damped Porous–Elastic System with Infinite Memory and Distributed Delay Terms. Math. Comput. Appl. 2021, 26, 71. https://doi.org/10.3390/mca26040071

AMA Style

Moumen A, Ouchenane D, Bouhali K, Altayeb Y. Well-Posedness and Stability Results for a Nonlinear Damped Porous–Elastic System with Infinite Memory and Distributed Delay Terms. Mathematical and Computational Applications. 2021; 26(4):71. https://doi.org/10.3390/mca26040071

Chicago/Turabian Style

Moumen, Abdelkader, Djamel Ouchenane, Keltoum Bouhali, and Yousif Altayeb. 2021. "Well-Posedness and Stability Results for a Nonlinear Damped Porous–Elastic System with Infinite Memory and Distributed Delay Terms" Mathematical and Computational Applications 26, no. 4: 71. https://doi.org/10.3390/mca26040071

APA Style

Moumen, A., Ouchenane, D., Bouhali, K., & Altayeb, Y. (2021). Well-Posedness and Stability Results for a Nonlinear Damped Porous–Elastic System with Infinite Memory and Distributed Delay Terms. Mathematical and Computational Applications, 26(4), 71. https://doi.org/10.3390/mca26040071

Article Metrics

Back to TopTop