Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Q-Learning-Based Dumbo Octopus Algorithm for Parameter Tuning of Fractional-Order PID Controller for AVR Systems
Previous Article in Journal
Analytical C2 Continuous Surface Blending
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Existence and Stability Results for Thermodiffusion Laminated Beam System with Delay Feedback

1
Laboratory of Pure and Applied Mathematics, Amar Teledji Laghouat University, Laghouat 03000, Algeria
2
Department of Mathematics, University of Khemis Miliana, Road of Theniet El-Had, Khemis Miliana 44225, Algeria
3
Energy and Intelligent Systems Laboratory, University of Khemis Miliana, Khemis Miliana 44225, Algeria
4
Computer Science and Mathematics Laboratory, University of Laghouat, Laghouat 03000, Algeria
5
Department of Mathematics, College of Sciences, Qassim University, Buraydah 51452, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(19), 3097; https://doi.org/10.3390/math12193097
Submission received: 2 September 2024 / Revised: 25 September 2024 / Accepted: 28 September 2024 / Published: 3 October 2024

Abstract

:
In this paper, a one-dimensional thermodiffusion laminated beam system with delay feedback is studied. The existence of a solution for our system is discussed within the context of the semigroup approach. In addition, under different boundary conditions, two results of stability properties independent of initial data are investigated.

1. Introduction

The aim of this study is to investigate the following system, for ( x , t ) ( 0 , 1 ) × ( 0 , ) :
ρ ϱ t t + G ϖ ϱ x x + μ ϱ t ( x , t τ ) = 0 , I ρ ( 3 s ϖ ) t t D ( 3 s ϖ ) x x G ϖ ϱ x γ 1 θ x γ 2 P x = 0 , I ρ s t t D s x x + G ϖ ϱ x + 4 3 γ s + 4 3 β s t = 0 , c θ t + d P t κ θ x x γ 1 ( 3 s ϖ ) x t = 0 , d θ t + r P t P x x γ 2 ( 3 s ϖ ) x t = 0 ,
along with the initial conditions ( x ( 0 , 1 ) ) :
ϱ ( x , t = 0 ) = ϱ 0 ( x ) , ϖ ( x , t = 0 ) = ϖ 0 ( x ) , s ( x , t = 0 ) = s 0 ( x ) , θ ( x , t = 0 ) = θ 0 ( x ) , P ( x , t = 0 ) = P 0 ( x ) , ϱ t ( x , t = 0 ) = ϱ 1 ( x ) , ϖ t ( x , t = 0 ) = ϖ 1 ( x ) , s t ( x , t = 0 ) = s 1 ( x ) ϱ t ( x , t ) = f 0 ( x , t ) ,
where τ > 0 is the time delay, and the non-zero real number μ > 0 is the weight of delay. With some restrictions on the parameters, we take into account the two boundary conditions listed below:
Boundary I ( t > 0 ):
ϱ x ( x = 0 , t ) = ϖ ( x = 0 , t ) = s ( x = 0 , t ) = θ ( x = 0 , t ) = P ( x = 0 , t ) = 0 , ϱ ( x = 1 , t ) = ϖ x ( x = 1 , t ) = s x ( x = 1 , t ) = θ ( x = 1 , t ) = P ( x = 1 , t ) = 0 ,
and Boundary II ( t > 0 ):
ϱ ( x = 0 , t ) = ϖ ( x = 0 , t ) = s ( x = 0 , t ) = θ ( x = 0 , t ) = P ( x = 0 , t ) = 0 , ϱ ( x = 1 , t ) = ϖ ( x = 1 , t ) = s ( x = 1 , t ) = θ ( x = 1 , t ) = P ( x = 1 , t ) = 0 .
γ 1 , γ 2 , r , c and d are positive constants. We assume that the wave speeds are identical:
ρ G = I ρ D .
We suppose also that the symmetric matrix
Γ = c d d r
is positive definite. Then, we deduce that
δ = c r d 2 > 0 .
Thus, for all θ and P, we have
c θ 2 + 2 d θ P + r P 2 > 0 .
It might be assumed that dissipation can only be explained through thermal conductivity in Timoshenko beams. However, research stemming from the development of advanced technologies after World War II has shown that the role of diffusion in solids cannot be ignored. This leads to a key question: what occurs when both diffusion and thermal effects are considered in laminated beams? Diffusion is the random movement of particles from areas of higher concentration to lower concentration. In elastic solids, thermodiffusion results from the interaction between strain, temperature, and mass diffusion fields. Heat and mass diffusion processes are crucial in various engineering applications, such as satellite systems, returning spacecraft, and aircraft landings on water or land. Currently, there is a strong interest in diffusion processes within the production of integrated circuits, resistors, semiconductor substrates, and MOS transistors. Oil companies are also focusing on this phenomenon to optimize oil extraction conditions.
First, we evaluate similar findings in the literature. Apalara [1] studied
ρ φ t t + G ( ψ φ x ) x = 0 , in 0 , 1 × 0 , + , I ρ ( ψ + 3 s ) t t D ( ψ + 3 s ) x x G ( ψ φ x ) = 0 , in 0 , 1 × 0 , + , I ρ s t t D s x x + 3 G ( ψ φ x ) + 4 γ s + δ θ x = 0 , in 0 , 1 × 0 , + , δ s x t λ θ x x + ρ 3 θ t = 0 , in 0 , 1 × 0 , + ,
Its soultion is exponential stable when G ρ = D I ρ . Aouadi et al. [2] investigated the following Timoshenko beam with thermodiffusion effects:
ρ 1 φ t t k φ x + ψ x = 0 , in 0 , L × 0 , + , ρ 2 ψ t t α ψ x x + k φ x + ψ γ 1 θ x γ 2 P x = 0 , in 0 , L × 0 , + , c θ t + d P t k θ x x γ 1 ψ x t = 0 , in 0 , L × 0 , + , d θ t + r P t h P x x γ 2 ψ x t = 0 , in 0 , L × 0 , + ,
They proved the existence and exponential stability result. We refer the reader to [3,4,5,6] and the references cited therein for more related results.
This work is new due to the idea of mass diffusion in laminated beams, which can have important physical consequences in addition to deformations of the body. In micro-beam resonators, for example, the current research has concentrated on the effect of mass diffusion on the damping ratio (see, e.g., [7]). Moreover, in addition to the traditional critical thickness linked to thermoelastic damping, mass diffusion provides a new one. The aforementioned arguments imply that mass diffusion will be essential to understanding laminated beam models’ thermomechanical behavior. According to the literature review, no numerical or theoretical studies investigating the impact of mass diffusion on the thermal vibrations of laminated beams have been carried out. This work therefore aims to explore the impact of mass diffusion in addition to temperature impacts on laminated beam behavior.
The structure of the paper is as follows: Section 2 provides our first main result about well posedness and the procedure used to achieve this result is a semigroup approach. For the two types of boundary conditions, we shall state and prove the exponential stability in Section 3.

2. Well-Posedness Results of the System

In order to show the well-posedness result of system (1), (2), and (4) or problems (1)–(3), we denote ξ = 3 s ϖ ; then, we introduce this new variable as follows:
z x , σ , t = ϱ t x , t τ σ , x 0 , 1 , σ 0 , 1 , t > 0 .
Then, we obtain the following equation:
τ z t x , σ , t + z σ x , σ , t = 0 , x , σ , t 0 , 1 × 0 , 1 × 0 , + .
Then, system (1) can be transformed as follows:
ρ ϱ t t + G 3 s ξ ϱ x x + μ z ( x , 1 , t ) = 0 , I ρ ξ t t D ξ x x G 3 s ξ ϱ x γ 1 θ x γ 2 P x = 0 , I ρ s t t D s x x + G 3 s ξ ϱ x + 4 3 γ s + 4 3 β s t = 0 , c θ t + d P t κ θ x x γ 1 ξ x t = 0 , d θ t + r P t P x x γ 2 ξ x t = 0 , τ z t x , σ , t + z σ x , σ , t = 0 ,
With the common initial data ( x ( 0 , 1 ) ), have the following:
ϱ ( x , t = 0 ) = ϱ 0 ( x ) , ξ ( x , t = 0 ) = ξ 0 ( x ) , s ( x , t = 0 ) = s 0 ( x ) , θ ( x , 0 ) = θ 0 ( x ) , P ( x , t = 0 ) = P 0 ( x ) , ϱ t ( x , t = 0 ) = ϱ 1 ( x ) , ξ t ( x , t = 0 ) = ξ 1 ( x ) , s t ( x , 0 ) = s 1 ( x ) , z ( x , 0 , t ) = ϱ t ( x , t ) , z ( x , σ , 0 ) = f 0 ( x , σ t ) , ( σ , t ) 0 , 1 × 0 , + .
The boundary conditions are taken in two different cases as
Boundary I ( t > 0 ) :
ϱ x ( x = 0 , t ) = ξ ( x = 0 , t ) = s ( x = 0 , t ) = θ ( x = 0 , t ) = P ( x = 0 , t ) = 0 , ϱ ( x = 1 , t ) = ξ x ( x = 1 , t ) = s x ( x = 1 , t ) = θ ( x = 1 , t ) = P ( x = 1 , t ) = 0 ,
and Boundary II ( t > 0 ) :
ϱ ( x = 0 , t ) = ξ ( x = 0 , t ) = s ( x = 0 , t ) = θ ( x = 0 , t ) = P ( x = 0 , t ) = 0 , ϱ ( x = 1 , t ) = ξ ( x = 1 , t ) = s ( x = 1 , t ) = θ ( x = 1 , t ) = P ( x = 1 , t ) = 0 ,
Let u = ϱ t , v = ξ t and y = s t , and introducing the vector function
U ( t ) = ϱ , u , ξ , v , s , y , θ , P , z T ,
we define the spaces
H a 1 ( 0 , 1 ) = y H 1 ( 0 , 1 ) : y ( x = 0 ) = 0 , H b 1 ( 0 , 1 ) = v H 1 ( 0 , 1 ) : v ( x = 1 ) = 0 , H b 2 ( 0 , 1 ) = y H 2 ( 0 , 1 ) : y x ( x = 0 ) = 0 , H a 2 ( 0 , 1 ) = v H 2 ( 0 , 1 ) : v x ( x = 1 ) = 0 ,
as well as the phase spaces H 1 and H 2 with respect to Boundary I and Boundary II as
H 1 = H b 1 ( 0 , 1 ) × L 2 ( 0 , 1 ) × H a 1 ( 0 , 1 ) × L 2 ( 0 , 1 ) 2 × L 2 ( 0 , 1 ) 2 × L 2 ( 0 , 1 ) × ( 0 , 1 ) ,
and
H 2 = H 0 1 ( 0 , 1 ) × L 2 ( 0 , 1 ) 3 × L 2 ( 0 , 1 ) 2 × L 2 ( 0 , 1 ) × ( 0 , 1 ) ,
Both state spaces are equipped with the same inner product, which is given by the following:
U , U ˜ H i = 0 1 [ ρ u u ˜ + G 3 s ξ ϱ x 3 s ˜ ξ ˜ ϱ ˜ x + I ρ v v ˜ + D ξ x ξ ˜ x + 3 I ρ y y ˜ + 3 D s x s ˜ x + 4 γ s s ˜ + c θ θ ˜ + d θ ˜ P + d θ P ˜ + r P P ˜ ] d x + τ | μ | 0 1 0 1 z ( x , σ ) z ˜ ( x , σ ) d σ d x .
This is defined for all U ( t ) = ( ϱ , u , ξ , v , s , y , θ , P , z ) , U ˜ ( t ) = ( ϱ ˜ , u ˜ , ξ ˜ , v ˜ , s ˜ , y ˜ , θ ˜ , P ˜ , z ˜ ) H i , i = 1 , 2 .
Then, system (8) can be transformed into a first-order Cauchy problem:
d d t U ( t ) = A i U ( t ) , i = 1 , 2 , U ( 0 ) = U 0 = ϱ 0 , ϱ 1 , ξ 0 , ξ 1 , s 0 , s 1 , θ 0 , P 0 , f 0 T ,
Therein, the operator A i , i = 1 , 2 is given by
A i U = u 1 ρ G 3 s ξ ϱ x x + μ z ( x , 1 , t ) v 1 I ρ D ξ x x + G 3 s ξ ϱ x + γ 1 θ x + γ 2 P x y 1 I ρ D s x x G 3 s ξ ϱ x 4 3 γ s 4 β 3 y κ r δ θ x x d δ P x x + γ 1 r γ 2 d δ v x c δ P x x κ d δ θ x x + γ 2 c γ 1 d δ v x 1 τ z σ ( x , σ ) .
Therein, we have
D A 1 = { U H 1 / ϱ H b 2 ( 0 , 1 ) H b 1 ( 0 , 1 ) , ξ , s H a 2 ( 0 , 1 ) H a 1 ( 0 , 1 ) , u H b 1 ( 0 , 1 ) v , y H a 1 ( 0 , 1 ) , θ , P H 2 ( 0 , 1 ) H 0 1 ( 0 , 1 ) , z , z σ L 2 ( 0 , 1 ) × ( 0 , 1 ) } .
D A 2 = { U H 2 / ϱ , ξ , s H 2 ( 0 , 1 ) H 0 1 ( 0 , 1 ) , u , v , y H 0 1 ( 0 , 1 ) , θ , P H 2 ( 0 , 1 ) H 0 1 ( 0 , 1 ) , z , z σ L 2 ( 0 , 1 ) × ( 0 , 1 ) } .
The first result related to the well posedness in both boundary conditions is presented in the following theorem.
Theorem 1.
For U 0 H 1 (resp. U 0 H 2 ), the problem (1)−(3) (resp. (1)−(2)) and (4) has a unique weak solution:
U ( t ) C [ 0 , ) , H 1 r e s p . U ( t ) C [ 0 , ) , H 2 .
In addition, if U 0 D A 1 (resp. U 0 D A 2 ), then
U ( t ) C [ 0 , ) , D A 1 C 1 [ 0 , ) , H 1 r e s p . U ( t ) C [ 0 , ) , D A 2 C 1 [ 0 , ) , H 2 .
Proof. 
In order to prove Theorem 1, we use the semigroup approach. That is, this is done to show that the operator A i generates a C 0 semigroup in H i . At this stage, we first prove the dissipativity of the operator
B i = A i c 1 I , i = 1 , 2 .
For
U = ( ϱ , u , ξ , v , s , y , θ , P , z ) D A i , i = 1 , 2 ,
we have
A i U ( t ) , U ( t ) H i = 0 1 ρ u G ρ 3 s ξ ϱ x x μ ρ z ( x , 1 ) d x + 0 1 G 3 s ξ ϱ x 3 y v u x d x + 0 1 I ρ v D I ρ ξ x x + G I ρ 3 s ξ ϱ x + γ 1 I ρ θ x + γ 2 I ρ P x d x + D 0 1 ξ x v x d x + 3 I ρ 0 1 y D I ρ s x x G I ρ 3 s ξ ϱ x 4 γ 3 I ρ s 4 β 3 I ρ y d x + 3 D 0 1 s x y x d x + 4 γ 0 1 s y d x + c 0 1 θ κ r δ θ x x d δ P x x + γ 1 r γ 2 d δ v x d x + r 0 1 P c δ P x x κ d δ θ x x + γ 2 c γ 1 d δ v x d x + d 0 1 P κ r δ θ x x d δ P x x + γ 1 r γ 2 d δ v x d x + d 0 1 θ c δ P x x κ d δ θ x x + γ 2 c γ 1 d δ v x d x | μ | 0 1 0 1 z ( x , σ ) z σ ( x , σ ) d σ d x .
Using integration by parts and taking into account the boundary conditions, we obtain
A i U ( t ) , U ( t ) H i κ 0 1 θ x 2 d x 0 1 P x 2 d x 4 β 0 1 y 2 d x + | μ | 0 1 u 2 d x , | μ | 0 1 u 2 d x , c 1 U ( t ) , U ( t ) H i ,
where c 1 = max ( | μ | , 1 ) min ( ϱ , G , 3 I ρ , 3 D , 4 δ , c , 2 d , r , τ | μ | ) , which shows that the operator A i c 1 I ( i = 1 , 2 ) is dissipative.
Let
F = f 1 , , f 9 H i ,
and we seek a solution U D A i of
I d A i U = F .
Then, we obtain
ϱ u = f 1 , ρ u + G 3 s ξ ϱ x x + μ z ( x , 1 ) = ρ f 2 , ξ v = f 3 , I ρ v D ξ x x G 3 s ξ ϱ x γ 1 θ x γ 2 P x = I ρ f 4 , s y = f 5 , 4 β + 3 I ρ y 3 D s x x + 3 G 3 s ξ ϱ x + 4 γ s = 3 I ρ f 6 , δ θ κ r θ x x + d P x x γ 1 r γ 2 d v x = δ f 7 , δ P c P x x + κ d θ x x γ 2 c γ 1 d v x = δ f 8 τ z + z σ = τ f 9 .
Suppose that ϱ is found with the appropriate regularity. Therefore, the first equation in (10) yields
u = ϱ f 1 .
It is clear that u H b 1 0 , 1 if U H 1 and u H 0 1 0 , 1 if U H 2 . Furthermore, we can find z as
z x , 0 = u x , for x 0 , 1 .
We obtain using (10)9
z x , σ = u x e σ τ + τ e σ τ 0 σ f 9 x , η e η τ d η .
From (11), we have
z x , σ = ϱ x e σ τ f 1 e σ τ + τ e σ τ 0 σ f 9 x , η e η τ d η .
From (13), we have
z x , 1 = ϱ x e τ + z 0 x ,
where x 0 , 1 , and
z 0 x = f 1 e τ + τ e τ 0 1 f 9 x , η e η τ d η .
It is clear from the above formula that z 0 depends only on f i , i = 1 , 9 .
By using (14), (12) in (10), and substituting (10)1 into (10)2, (10)3 into (10)4, (10)5 into (10)6, (10)7 into (10)8, we obtain
ρ ϱ + G 3 s ξ ϱ x x + μ ϱ e τ = ρ ( f 1 + f 2 ) + μ f 1 e τ μ τ e τ 0 1 f 9 x , η e η τ d η , I ρ ξ D ξ x x G 3 s ξ ϱ x γ 1 θ x γ 2 P x = I ρ ( f 4 + f 3 ) , 4 β + 3 I ρ s 3 D s x x + 3 G 3 s ξ ϱ x + 4 γ s = 3 I ρ f 6 + 4 β + 3 I ρ f 5 , δ θ κ r θ x x + d P x x γ 1 r γ 2 d ξ x = δ f 7 γ 1 r γ 2 d f 3 x , δ P c P x x + κ d θ x x γ 2 c γ 1 d ξ x = δ f 8 γ 2 c γ 1 d f 3 x .
Multiplying (16)1 by ϱ ¯ , (16)2 by ξ ¯ , (16)3 by s ¯ , (16)4 by c δ θ ¯ , (16)5 by r δ P ¯ , (16)4 by d δ P ¯ and (16)5 by d δ θ ¯ and integrating their sum over ( 0 , 1 ) to obtain the subsequent variations formulation, we have
B ( ϱ , ξ , s , θ , P ) , ( ϱ ¯ , ξ ¯ , s ¯ , θ ¯ , P ¯ ) = L ( ϱ ¯ , ξ ¯ , s ¯ , θ ¯ , P ¯ ) ,
where the bilinear form
B : H b 1 ( 0 , 1 ) × H a 1 ( 0 , 1 ) × H a 1 ( 0 , 1 ) × L 2 ( 0 , 1 ) × L 2 ( 0 , 1 ) 2 R ,
or
B : H 0 1 ( 0 , 1 ) × H 0 1 ( 0 , 1 ) × H 0 1 ( 0 , 1 ) × L 2 ( 0 , 1 ) × L 2 ( 0 , 1 ) 2 R ,
is given by
B ( ( ϱ , ξ , s , θ , P ) , ( ϱ ¯ , ξ ¯ , s ¯ , θ ¯ , P ¯ ) ) = ( ρ + μ e τ ) 0 1 ϱ ϱ ¯ d x + G 0 1 3 s ξ ϱ x 3 s ¯ ξ ¯ ϱ ¯ x d x + I ρ 0 1 ξ ξ ¯ d x + D 0 1 ξ x ξ ¯ x d x + 4 β + 3 I ρ + 4 γ 0 1 s s ¯ d x + 3 D 0 1 s x s ¯ x d x + γ 1 0 1 θ ξ ¯ x d x γ 1 0 1 ξ x θ ¯ d x + γ 2 0 1 P ξ ¯ x d x γ 2 0 1 ξ x P ¯ d x + c 0 1 θ θ ¯ d x + κ 0 1 θ x θ ¯ x d x + r 0 1 P P ¯ d x + 0 1 P x P ¯ x d x + d 0 1 θ P ¯ d x + d 0 1 P θ ¯ d x ,
and the linear form
L : H b 1 ( 0 , 1 ) × H a 1 ( 0 , 1 ) × H a 1 ( 0 , 1 ) × L 2 ( 0 , 1 ) × L 2 ( 0 , 1 ) R ,
or
L : H 0 1 ( 0 , 1 ) × H 0 1 ( 0 , 1 ) × H 0 1 ( 0 , 1 ) × L 2 ( 0 , 1 ) × L 2 ( 0 , 1 ) R ,
is defined by
L ( ϱ ¯ , ξ ¯ , s ¯ , θ ¯ , P ¯ ) = ρ 0 1 f 1 + f 2 μ ρ z 0 ϱ ¯ d x + I ρ 0 1 f 3 + f 4 ξ ¯ d x + 4 β + 3 I ρ 0 1 f 5 s ¯ d x + 3 I ρ 0 1 f 6 s ¯ d x + c 0 1 f 7 θ ¯ d x γ 1 0 1 f 3 x θ ¯ d x + r 0 1 f 8 P ¯ d x γ 2 0 1 f 3 x P ¯ d x + d 0 1 f 7 P ¯ d x + d 0 1 f 8 θ ¯ d x .
It is easy to verify that B is continuous and coercive, and L is continuous. So, applying the Lax–Milgram theorem, we deduce that there exists a positive constant C such that
B ( ϱ , ξ , s , θ , p ) , ( ϱ , ¯ ξ ¯ , s ¯ , θ ¯ , p ¯ ) C | | ( ϱ , ξ , s , θ , p ) | | H i 2 , i = 1 , 2 .
This shows that B is coercive on H i × H i , i = 1 , 2 .
Owing to the Lax–Milgram Theorem, we find that (17) has a one solution:
ϱ H b 1 ( 0 , 1 ) , ξ , s H a 1 ( 0 , 1 ) , θ , P L 2 ( 0 , 1 ) ,
or
ϱ , ξ , s H 0 1 ( 0 , 1 ) , θ , P L 2 ( 0 , 1 ) .
To achieve more regularity, we should take
( ξ ¯ , s ¯ , θ ¯ , P ¯ ) = ( 0 , 0 , 0 , 0 ) H a 1 ( 0 , 1 ) × H a 1 ( 0 , 1 ) × L 2 ( 0 , 1 ) × L 2 ( 0 , 1 ) .
From (17), we have
( ρ + μ e τ ) 0 1 ϱ ϱ ¯ d x G 0 1 3 s ξ ϱ x ϱ ¯ x = ρ 0 1 f 1 + f 2 μ ρ z 0 ϱ ¯ d x , ϱ ¯ H b 1 ,
which implies that
G ϱ x x = ( ρ + μ e τ ) ϱ + G ( 3 s ξ ) x ρ f 1 + f 2 μ ρ z 0 L 2 ( 0 , 1 ) .
Then, we have
ϱ H 2 ( 0 , 1 ) .
As opposed to that, for any φ C 1 [ 0 , 1 ] with φ ( x = 1 ) = 0 , we find that (19) also holds true. Then,
( ρ + μ e τ ) 0 1 ϱ φ d x G 0 1 3 s ξ ϱ x φ x ρ 0 1 f 1 + f 2 μ ρ z 0 φ d x = 0 , φ C 1 [ 0 , 1 ] .
Utilizing integration by parts, we obtain
G ( 3 s ξ ϱ x ) ( 1 ) φ ( 1 ) G ϱ x ( 0 ) φ ( 0 ) + ( ρ + μ e τ ) 0 1 ϱ φ d x + G 0 1 ( 3 s ξ ϱ x ) x φ d x ρ 0 1 f 1 + f 2 μ ρ z 0 φ d x = 0 .
Thus,
G ϱ x ( x = 0 ) φ ( x = 0 ) = 0 , φ C 1 [ 0 , 1 ] .
Since φ is arbitrary, we can obtain ϱ x ( x = 0 ) = 0 ; thus,
ϱ H b 2 ( 0 , 1 ) H b 1 ( 0 , 1 ) .
By using the same approach, we can arrive at
ξ , s H a 2 ( 0 , 1 ) H a 1 ( 0 , 1 ) , θ , P H 2 ( 0 , 1 ) H 0 1 ( 0 , 1 ) .
So, by applying the Lax–Milgram theorem, our problem has a unique solution. Consequently, ϱ , ξ , s , θ , P D A 1 and A 1 are maximal likewise, so we may acquire ( ϱ , ξ , s , θ , P ) D A 2 and A 2 as maximal. Therefore, the conclusion of the proof follows from the Lumer–Philips Theorem (see [8] ). □

3. Exponential Stability

The stability of problem (1) arises very concretely during the qualitative study, since the presence of a delay most often leads to oscillations or even instabilities. If the mathematical nature of the delay phenomenon is not taken into account, the only thermal diffusion effect is to provide the stability behavior to permanently accelerate dynamic performances. Various methods are available. The most practical method is the Lyapunov approach. We will use it here because it generally serves as a basis for suitable methods of investigating stability. To begin with, we define the energy of the system as
E ( t ) = 1 2 0 1 ρ ϱ t 2 + G 3 s ξ ϱ x 2 + I ρ ξ t 2 + D ξ x 2 + 3 I ρ s t 2 + 3 D s x 2 + 4 γ s 2 + c θ 2 + 2 d θ P + r P 2 d x + τ | μ | 2 0 1 0 1 z 2 ( x , σ , t ) d σ d x .
Then, according to (7), the energy E is positive.
It is not hard to see for the two boundary conditions that
d d t E ( t ) κ 0 1 θ x 2 d x 0 1 P x 2 d x 4 β 0 1 s t 2 d x + | μ | 0 1 ϱ t 2 d x .
The stability result is stated in the following theorem.
Theorem 2.
Let U be a regular solution of problem (1)–(3) or (1), (2), and (4). Let (5) hold. Then, there exist positive constants μ ¯ , a 1 and a 2 , independent of initial data, such that | μ | < μ ¯ and
E ( t ) a 2 E ( 0 ) e a 1 t , t > 0 .
To prove Theorem 2, with the two boundary conditions, we shall divide into the following two subsections.

3.1. Exponential Decay with Boundary I

We need the following several lemmas.
Lemma 1.
Let the functional F 1 be defined by
F 1 ( t ) = I ρ 0 1 ξ t φ d x ,
where
γ 1 φ x = c θ + d P , φ ( x = 0 ) = φ ( x = 1 ) = 0 .
Then, the there exist constants c 1 > 0 and c 2 > 0 such that for any ε 1 > 0 and ε 2 > 0 , we have
d d t F 1 ( t ) I ρ 2 0 1 ξ t 2 d x + G ε 1 2 0 1 3 s ξ ϱ x 2 d x + D ε 2 2 0 1 ξ x 2 d x + c 1 0 1 θ x 2 d x + c 2 0 1 P x 2 d x .
Proof. 
Differentiate F 1 and use (8)2 to obtain
d d t F 1 ( t ) = 0 1 D ξ x x + G 3 s ξ ϱ x + γ 1 θ x + γ 2 P x φ d x + I ρ 0 1 ξ t φ t d x = G 0 1 3 s ξ ϱ x φ d x D 0 1 ξ x φ x d x + 0 1 γ 1 θ x + γ 2 P x φ d x + I ρ 0 1 ξ t φ t d x .
Using Cauchy–Schwarz and Poincare’s inequalities, we obtain ε 1 > 0
G 0 1 3 s ξ ϱ x φ d x ε 1 G 2 0 1 3 s ξ ϱ x 2 d x + G 2 ε 1 0 1 φ x 2 d x ε 1 G 2 0 1 3 s ξ ϱ x 2 d x + G 2 ε 1 0 L c γ 1 θ + d γ 1 P 2 d x ε 1 G 2 0 1 3 s ξ ϱ x 2 d x + G d 2 β 2 ε 1 γ 1 2 0 1 P x 2 d x
+ G c 2 β 1 ε 1 γ 1 2 0 1 θ x 2 d x ,
and for any ε 2 > 0
D 0 1 ξ x φ x d x D ε 2 2 0 1 ξ x 2 d x + D c 2 β 3 ε 2 γ 1 2 0 1 θ x 2 d x + D d 2 β 4 ε 2 γ 1 2 0 1 P x 2 d x ,
and
0 1 γ 1 θ x + γ 2 P x φ d x 1 2 0 1 γ 1 θ x + γ 2 P x 2 d x + 1 2 0 1 φ 2 d x γ 1 2 + c 2 β 5 γ 1 2 0 1 θ x 2 d x + γ 2 2 + β 6 d 2 γ 1 2 0 1 P x 2 d x .
Using (8)4 and Young’s inequality, we have
I ρ 0 1 ξ t φ t d x = I ρ γ 1 0 1 ξ t 0 x κ θ y y + γ 1 ξ t y d y d x = I ρ 0 1 ξ t 2 d x I ρ κ γ 1 0 1 ξ t θ x d x I ρ 2 0 1 ξ t 2 d x + I ρ κ 2 2 γ 1 2 0 1 θ x 2 d x .
Inserting (24)–(28) into (23), we establish (22) with
c 1 = I ρ κ 2 2 γ 1 2 + G c 2 β 1 ε 1 γ 1 2 + D c 2 β 3 ε 2 γ 1 2 + γ 1 2 + c 2 β 5 γ 1 2 ,
and
c 2 = G d 2 β 2 ε 1 γ 1 2 + D d 2 β 4 ε 2 γ 1 2 + γ 2 2 + β 6 d 2 γ 1 2 .
Lemma 2.
The functional F 2 defined by
F 2 ( t ) = I ρ 0 1 ξ t ξ d x ,
satisfies, for any t > 0 ,
d d t F 2 ( t ) D 2 0 1 ξ x 2 d x + I ρ 0 1 ξ t 2 d x + G ϵ 3 2 0 1 3 s ξ ϱ x 2 d x + 2 γ 1 2 D 0 1 θ x 2 d x + 2 γ 2 2 D 0 1 P x 2 d x .
Proof. 
By using (8)2, we obtain
d d t F 2 ( t ) = I ρ 0 1 ξ t 2 d x + G 0 1 3 s ξ ϱ x ξ d x D 0 1 ξ x 2 d x + 0 1 γ 1 θ x + γ 2 P x ξ d x .
Making use of Young’s and Poincare’s inequalities, we obtain
G 0 1 3 s ξ ϱ x ξ d x G 2 ϵ 3 0 1 ξ 2 d x + G ϵ 3 2 0 1 3 s ξ ϱ x 2 d x G c 3 2 ϵ 3 0 1 ξ x 2 d x + G ϵ 3 2 0 1 3 s ξ ϱ x 2 d x ,
and
0 1 γ 1 θ x + γ 2 P x ξ d x D 4 0 1 ξ x 2 d x + 2 γ 1 2 D 0 1 θ x 2 d x + 2 γ 2 2 D 0 1 P x 2 d x .
By using (30)–(31) and choosing c 3 = D ϵ 3 2 G , we obtain (29). □
Lemma 3.
Let F 3 a functional defined by
F 3 ( t ) = I ρ 0 1 ξ t 3 s ξ ϱ x d x + ρ D G 0 1 ξ x ϱ t d x .
Assume that the condition (5) holds; then, F 3 satisfies, for any t > 0 ,
d d t F 3 ( t ) G 2 0 1 3 s ξ ϱ x 2 d x + 2 I ρ 0 1 ξ t 2 d x + 3 I ρ 4 0 1 s t 2 d x + γ 1 2 G 0 1 θ x 2 d x + γ 2 2 G 0 1 P x 2 d x + D 2 2 G 2 0 1 ξ x 2 d x + μ 2 2 0 1 z 2 ( x , 1 , t ) d x .
Proof. 
Using the first and the second equations of (8) and taking the derivative of F 3 , we have
d d t F 3 ( t ) = 0 1 G 3 s ξ ϱ x D ξ x x γ 1 θ x γ 2 P x 3 s ξ ϱ x d x I ρ 0 1 ξ t 3 s ξ ϱ x t d x + ρ D G 0 1 ξ x t ϱ t d x + D G 0 1 ξ x G 3 s ξ ϱ x x μ z ( x , 1 , t ) d x = G 0 1 3 s ξ ϱ x 2 d x 0 1 γ 1 θ x + γ 2 P x 3 s ξ ϱ x d x I ρ 0 1 ξ t ( 3 s ξ ) t d x + ρ D G I ρ 0 1 ξ x t ϱ t d x D μ G 0 1 ξ x z ( x , 1 , t ) d x .
Young’s inequality gives us
0 1 γ 1 θ x + γ 2 P x 3 s ξ ϱ x d x G 2 0 1 3 s ξ ϱ x 2 d x + 1 2 G 0 1 γ 1 θ x + γ 2 P x 2 d x G 2 0 1 3 s ξ ϱ x 2 d x + γ 1 2 G 0 1 θ x 2 d x + γ 2 2 G 0 1 P x 2 d x ,
and
D μ G 0 1 ξ x z ( x , 1 , t ) d x D 2 2 G 2 0 1 ξ x 2 d x + μ 2 2 0 1 z 2 ( x , 1 , t ) d x ,
and
I ρ 0 1 ξ t ( 3 s ξ ) t = I ρ 0 1 ξ t 2 d x 3 I ρ 0 1 ξ t s t d x 2 I ρ 0 1 ξ t 2 d x + 3 I ρ 4 0 1 s t 2 d x .
Replacing (34)–(36) into (33) and noting the condition (5), we obtain (32). □
Lemma 4.
Let the functional F 4 be defined by
F 4 ( t ) = ρ 0 1 ϱ ϱ t d x ρ 0 1 ϱ t x 1 ( 3 s ξ ) ( y ) d y d x .
Then, for all t > 0 , we have
d d t F 4 ( t ) ρ 2 0 1 ϱ t 2 d x + G + 1 0 1 3 s ξ ϱ x 2 d x + ρ 0 1 ξ t 2 d x + 9 ρ 0 1 s t 2 d x + μ 2 4 0 1 z 2 ( x , 1 , t ) d x .
Proof. 
It follows from (8) that
d d t F 4 ( t ) = ρ 0 1 ϱ t 2 d x + G 0 1 3 s ξ ϱ x x ϱ d x + μ 0 1 x 1 ( 3 s ξ ) ( y ) d y + ϱ z ( x , 1 , t ) d x + G 0 1 3 s ξ ϱ x x x 1 ( 3 s ξ ) ( y ) d y d x ρ 0 1 ϱ t x 1 ( 3 s ξ ) t ( y ) d y d x = ρ 0 1 ϱ t 2 d x + G 0 1 3 s ξ ϱ x 2 d x + ρ 0 1 ϱ t x 1 ξ t ( y ) d y d x 3 ρ 0 1 ϱ t x 1 s t ( y ) d y d x + μ 0 1 x 1 ( 3 s ξ ) ( y ) d y + ϱ z ( x , 1 , t ) d x ,
thereby yielding the following estimates
ρ 0 1 ϱ t x 1 ξ t ( y ) d y d x ρ 4 0 1 ϱ t 2 d x + ρ 0 1 ξ t 2 d x ,
and
3 ρ 0 1 ϱ t x 1 s t ( y ) d y d x ρ 4 0 1 ϱ t 2 d x + 9 ρ 0 1 s t 2 d x ,
and
μ 0 1 x 1 ( 3 s ξ ) ( y ) d y + ϱ z ( x , 1 , t ) d x 0 1 x 1 ( 3 s ξ ) ( y ) d y + ϱ 2 d x + μ 2 4 0 1 z 2 ( x , 1 , t ) d x 0 1 3 s ξ ϱ x 2 d x + μ 2 4 0 1 z 2 ( x , 1 , t ) d x ,
which imply (37). □
Lemma 5.
Let the functional F 5 be defined by
F 5 ( t ) = 3 I ρ 0 1 s t s d x + 2 β 0 1 s 2 d x 3 ρ 0 1 ϱ t x 1 s ( y ) d y d x .
Then, for any t > 0 , we have
d d t F 5 ( t ) 3 D 0 1 s x 2 d x 3 γ 0 1 s 2 d x + 3 ρ 2 0 1 ϱ t 2 d x + 3 I ρ + 1 2 ρ 0 1 s t 2 d x + 9 μ 2 4 γ 0 1 z 2 ( x , 1 , t ) d x .
Proof. 
By taking the derivative of F 5 and using (8)3, we obtain
d d t F 5 ( t ) = 3 I ρ 0 1 s t 2 d x 3 D 0 1 s x 2 d x 4 γ 0 1 s 2 d x 3 ρ 0 1 ϱ t x 1 s t ( y ) d y d x + 3 μ 0 1 z ( x , 1 , t ) x 1 s ( y ) d y d x ,
along with the estimates
3 ρ 0 1 ϱ t x 1 s t ( y ) d y d x 3 ρ 2 0 1 ϱ t 2 d x + 3 ρ 2 0 1 s t 2 d x ,
and
3 μ 0 1 z ( x , 1 , t ) x 1 s ( y ) d y d x 9 μ 2 4 γ 0 1 z 2 ( x , 1 , t ) d x + γ 0 1 s 2 d x ,
so we establish (38). □
Lemma 6.
The functional F 6 defined by
F 6 ( t ) = τ 0 1 0 1 e σ τ z 2 ( x , σ , t ) d σ d x ,
satisfies for m 1 > 0 the estimate
d d t F 6 ( t ) m 1 0 1 z 2 ( x , 1 , t ) d x m 1 τ 0 1 0 1 z 2 ( x , σ , t ) d σ d x + 0 1 ϱ t 2 d x .
Proof. 
Differentiate F 6 and use the last equation in (8) and z ( x , 0 , t ) = ϱ t as follows:
d d t F 6 ( t ) = 2 0 1 0 1 e τ σ z ( x , σ , t ) z σ ( x , σ , t ) d σ d x = 0 1 0 1 d d σ e τ σ z 2 ( x , σ , t ) d σ d x τ 0 1 0 1 e τ σ z 2 ( x , σ , t ) d σ d x = 0 1 e τ z 2 ( x , 1 , t ) z 2 ( x , 0 , t ) d x τ 0 1 0 1 e τ σ z 2 ( x , σ , t ) d σ d x = 0 1 e τ z 2 ( x , 1 , t ) d x + 0 1 ϱ t 2 d x τ 0 1 0 1 e τ σ z 2 ( x , σ , t ) d σ d x .
Observe that, σ ( 0 , 1 ) , the relation e τ e σ τ 1 holds. Therefore, for some m 1 = e τ , we arrive at (39). □
We are now in position to introduce the Lyapunov functional L in the following:
L ( t ) = N E ( t ) + N 1 F 1 ( t ) + N 2 F 2 ( t ) + N 3 F 3 ( t ) + F 4 ( t ) + N 4 F 5 ( t ) + N 5 F 6 ( t ) ,
where N and N i , i = 1 , , 5 will be determined later.
Lemma 7.
There exist two positive constants ν 1 and ν 2 such that
ν 1 E ( t ) L ( t ) ν 2 E ( t ) .
Proof. 
Owing to the Young’s inequality, there is a constant N 0 > 0 such that
| L ( t ) N E ( t ) | N 0 E ( t ) .
Then, we can take N > 0 large (if needed) such that N N 0 > 0 to obtain (40) with ν 1 = N N 0 and ν 2 = N + N 0 . □
Proof. 
(Of Theorem 2 for Boundary I.) By using Lemma 1–Lemma 6 and taking
ε 1 = N 3 4 N 1 , ε 2 = 2 D N 1 , ε 3 = N 3 4 N 2 , N 5 = | μ | , N 4 = 1 6 ,
we have
d d t L ( t ) κ N c 1 N 1 2 γ 1 2 D N 2 γ 1 2 G N 3 0 1 θ x 2 d x D 2 N 2 D 2 2 G 2 N 3 1 0 1 ξ x 2 d x N c 2 N 1 2 γ 2 2 D N 2 γ 2 2 G N 3 0 L P x 2 d x ρ 4 ( N + 1 ) | μ | 0 1 ϱ t 2 d x 4 β N 3 I ρ 4 N 3 9 ρ I ρ 2 + ρ 4 0 1 s t 2 d x D 2 0 1 s x 2 d x I ρ 2 N 1 I ρ N 2 2 I ρ N 3 ρ 0 1 ξ t 2 d x γ 2 0 1 s 2 d x G 4 N 3 G + 1 0 1 3 s ξ ϱ x 2 d x m 1 | μ | τ 0 1 0 1 e σ τ z 2 ( x , σ , t ) d σ d x | μ | m 1 | μ | 3 8 γ + 1 4 + N 3 2 0 1 z 2 ( x , 1 , t ) d x .
We first take N 3 > 0 such that
G 4 N 3 G + 1 > 0 .
For fixed N 3 , we take N 2 > 0 such that
D 2 N 2 D 2 2 G 2 N 3 1 > 0 .
For fixed N 2 and N 3 , we choose N 1 > 0 large enough such that
I ρ 2 N 1 I ρ N 2 2 I ρ N 3 ρ > 0 .
Then, we choose N so large such that
κ N c 1 N 1 2 γ 1 2 D N 2 γ 1 2 G N 3 > 0 , N c 2 N 1 2 γ 2 2 D N 2 γ 2 2 G N 3 > 0 ,
and
4 β N 3 I ρ 4 N 3 9 ρ I ρ 2 + ρ 4 > 0 .
Finally, we pick
μ ¯ = min m 1 3 8 γ + 1 4 + N 3 2 , ρ 4 ( N + 1 ) .
Therefore, there exists a positive constant ν 0 such that
d d t L ( t ) ν 0 E ( t ) .
According to (40) and (41), we see that
E ( t ) ν 2 ν 1 E ( 0 ) e ν 0 ν 2 t .
This completes the proof. □

3.2. Exponential Decay with Boundary I I

The exponential decay of system (1) with Boundary I I is established in this subsection. The functional F 1 ( t ) still fulfills the estimates (29) for Boundary II, and this is simple to confirm. We first establish the following lemmas in order to circumvent the boundary’s difficulty.
Lemma 8.
Let the functional H 1 be defined by
H 1 ( t ) = ρ 0 1 ϱ ϱ t d x .
Then,
d d t H 1 ( t ) ρ 0 1 ϱ t 2 d x + 4 0 1 ξ x 2 d x + 36 0 1 s x 2 d x + 1 + G + G 2 4 0 1 3 s ξ ϱ x 2 d x + μ 2 2 0 1 z 2 ( x , 1 , t ) d x .
Proof. 
Differentiating H 1 with respect to t and using (8)1, we have
d d t H 1 ( t ) = ρ 0 1 ϱ t 2 d x + G 0 1 3 s ξ ϱ x 2 d x G 0 1 3 s ξ ϱ x ( 3 s ξ ) d x + μ 0 1 ϱ z ( x , 1 , t ) d x .
It follows from Young’s and Poincare’s inequalities that
G 0 1 3 s ξ ϱ x ( 3 s ξ ) d x 0 1 ( 3 s ξ ) x 2 d x + G 2 4 0 1 3 s ξ ϱ x 2 d x ,
μ 0 1 ϱ z ( x , 1 , t ) d x μ 2 2 0 1 z 2 ( x , 1 , t ) d x + 1 2 0 1 ϱ 2 d x ,
and
1 2 0 1 ϱ 2 d x 1 2 0 1 ϱ x 2 d x = 1 2 0 1 3 s ξ ( 3 s ξ ) ϱ x 2 d x 0 1 3 s ξ ϱ x 2 d x + 0 1 ( 3 s ξ ) x 2 d x .
On the other hand,
0 1 ( 3 s ξ ) x 2 d x 2 0 1 ξ x 2 d x + 18 0 1 s x 2 d x ,
Together with (43)–(45), this yields (42). □
Lemma 9.
Let the functional H 2 be defined by
H 2 ( t ) = ρ 0 1 ϱ t f d x + I ρ 0 1 ξ t ξ d x ,
where
f x x = ξ x , f ( x = 0 ) = f ( x = 1 ) = 0 .
Then, H 2 satisfies, for any ε 3 > 0 ,
d d t H 2 ( t ) D 2 0 1 ξ x 2 d x + ε 3 0 1 ϱ t 2 d x + I ρ + ρ 2 4 ε 3 0 1 ξ t 2 d x + 2 γ 1 2 D 0 1 θ x 2 d x + 2 γ 2 2 D 0 1 P x 2 d x + μ 2 2 0 1 z 2 ( x , 1 , t ) d x .
Proof. 
It is easy to obtain
d d t H 2 ( t ) = G 0 1 3 s ξ ϱ x f x d x μ 0 1 f z ( x , 1 , t ) d x + ρ 0 1 ϱ t f t d x + I ρ 0 1 ξ t 2 d x D 0 1 ξ x 2 d x + G 0 1 3 s ξ ϱ x ξ d x + 0 1 γ 1 θ x + γ 2 P x ξ d x = ρ 0 1 ϱ t f t d x + I ρ 0 1 ξ t 2 d x D 0 1 ξ x 2 d x + 0 1 γ 1 θ x + γ 2 P x ξ d x μ 0 1 f z ( x , 1 , t ) d x .
It follows from Young’s inequality that
0 1 γ 1 θ x + γ 2 P x ξ d x D 4 0 1 ξ x 2 d x + 2 γ 1 2 D 0 1 θ x 2 d x + 2 γ 2 2 D 0 1 P x 2 d x .
and
μ 0 1 f z ( x , 1 , t ) d x μ 2 2 0 1 z 2 ( x , 1 , t ) d x + 1 2 0 1 f 2 d x ,
and
ρ 0 1 ϱ t f t d x ε 3 0 1 ϱ t 2 d x + ρ 2 4 ε 3 0 1 ξ t 2 d x , ε 3 > 0 ,
where we used the fact that
1 2 0 1 f 2 d x 1 2 0 1 f x 2 d x D 4 0 1 ξ x 2 d x ,
and
0 1 f t 2 d x 0 1 f x t 2 d x 0 1 ξ t 2 d x .
The proof is completed. □
The functional F 3 in the following lemma is still defined as in Lemma (3), but because the boundary conditions have changed, we now obtain the new estimate shown below.
Lemma 10.
Assume that the condition (5) holds; then, the functional F 3 is satisfied for any ε 4 > 0 and ε 5 > 0
d d t F 3 ( t ) G 2 6 + 2 G ε 4 D 4 ε 4 ε 5 0 1 3 s ξ ϱ x 2 d x + 36 G + 126 ε 4 0 1 s x 2 d x + 2 ρ G ε 4 0 1 ϱ t 2 d x + 2 I ρ + I ρ D 2 ε 4 0 1 ξ t 2 d x + 3 I ρ 4 0 1 s t 2 d x + 14 + 4 G ε 4 + D 4 ε 4 3 D + G 2 ε 5 0 1 ξ x 2 d x + γ 1 2 G + γ 1 2 2 ε 4 0 1 θ x 2 d x + μ 2 2 1 + 2 ε 4 G 0 1 z 2 ( x , 1 , t ) d x + γ 2 2 G + γ 2 2 2 ε 4 0 1 P x 2 d x ε 4 G d d t 0 1 ρ g ( x ) ϱ t ϱ x d x D 4 ε 4 d d t 0 1 I ρ g ( x ) ξ t ξ x d x ,
where
g ( x ) = 2 4 x , x [ 0 , 1 ] .
Proof. 
Given the boundary conditions (4) and the identical arguments in (33), we have
d d t F 3 ( t ) = G 0 1 3 s ξ ϱ x 2 d x + D ξ x ϱ x | x = 0 x = 1 0 1 γ 1 θ x + γ 2 P x 3 s ξ ϱ x d x I ρ 0 1 ξ t ( 3 s ξ ) t d x + ρ D G I ρ 0 1 ξ t x ϱ t d x μ D G 0 1 ξ x z ( x , 1 , t ) d x .
According to (34) and (36), we obtain
d d t F 3 ( t ) G 2 0 1 3 s ξ ϱ x 2 d x + D ξ x ϱ x | x = 0 x = 1 + γ 1 2 G 0 1 θ x 2 d x + γ 2 2 G 0 1 P x 2 d x + ρ D G I ρ 0 1 ξ x t ϱ t d x + 2 I ρ 0 1 ξ t 2 d x + 3 I ρ 4 0 1 s t 2 d x + D 2 2 G 2 0 1 ξ x 2 d x + μ 2 2 0 1 z 2 ( x , 1 , t ) d x .
We will now estimate the term D ξ x ϱ x | x = 0 x = 1 . Using Young’s inequality, we have ε 4 > 0
D ξ x ϱ x | x = 0 x = 1 ε 4 ϱ x 2 ( 1 ) + ϱ x 2 ( 0 ) + D 2 4 ε 4 ξ x 2 ( 1 ) + ξ x 2 ( 0 ) .
It follows from (8)2 that
d d t 0 1 I ρ g ( x ) ξ t ξ x d x = 0 1 g ( x ) D ξ x x + G 3 s ξ ϱ x + γ 1 θ x + γ 2 P x ξ x d x + I ρ 0 1 g ( x ) ξ t ξ x t d x = D 2 g ( x ) ξ x 2 | x = 0 x = 1 D 2 0 1 g ( x ) ξ x 2 d x + G 0 1 g ( x ) 3 s ξ ϱ x ξ x d x + 0 1 g ( x ) γ 1 θ x + γ 2 P x ξ x d x I ρ 2 0 1 g ( x ) ξ t 2 d x ,
Notice that
g ( 0 ) = 2 , g ( 1 ) = 2 , g ( x ) = 4 ,
yields
d d t 0 1 I ρ g ( x ) ξ t ξ x d x = D ξ x 2 ( 1 ) + ξ x 2 ( 0 ) + 2 D 0 1 ξ x 2 d x + G 0 1 g ( x ) 3 s ξ ϱ x ξ x + 2 I ρ 0 1 ξ t 2 d x + 0 1 g ( x ) γ 1 θ x + γ 2 P x ξ x d x .
From Young’s inequality, we obtain for any ε 5 > 0 the following:
G 0 1 g ( x ) 3 s ξ ϱ x ξ x d x ε 5 0 1 3 s ξ ϱ x 2 d x + G 2 ε 5 0 1 ξ x 2 d x ,
and
0 1 g ( x ) γ 1 θ x + γ 2 P x ξ x d x D 0 1 ξ x 2 d x + 2 γ 1 2 D 0 1 θ x 2 d x + 2 γ 2 2 D 0 L P x 2 d x .
Replacing (52) and (53) in (51), we see that
D ξ x 2 ( 1 ) + ξ x 2 ( 0 ) d d t 0 1 I ρ g ( x ) ξ t ξ x d x + ε 5 0 1 3 s ξ ϱ x 2 d x + 3 D + G 2 ε 5 0 1 ξ x 2 d x + 2 γ 1 2 D 0 1 θ x 2 d x + 2 γ 2 2 D 0 1 P x 2 d x + 2 I ρ 0 1 ξ t 2 d x .
On the other hand,
d d t 0 1 ρ g ( x ) ϱ t ϱ x d x = G 0 1 g ( x ) ϱ x x ϱ x d x G 0 1 g ( x ) ( 3 s ξ ) x ϱ x d x + ρ 0 1 g ( x ) ϱ t ϱ x t d x μ 0 1 g ( x ) ϱ x z ( x , 1 , t ) d x = G ϱ x 2 ( 1 ) + ϱ x 2 ( 0 ) + 2 G 0 1 ϱ x 2 d x G 0 1 g ( x ) ( 3 s ξ ) x ϱ x d x + 2 ρ 0 1 ϱ t 2 d x μ 0 1 g ( x ) ϱ x z ( x , 1 , t ) d x .
Then, the following estimates
G 0 1 g ( x ) ( 3 s ξ ) x ϱ x d x G 0 1 ϱ x 2 d x + G 0 1 ( 3 s ξ ) x 2 d x ,
and
μ 0 1 g ( x ) ϱ x z ( x , 1 , t ) d x 0 1 ϱ x 2 d x + μ 2 0 1 z 2 ( x , 1 , t ) d x ,
give us
d d t 0 1 ρ g ( x ) ϱ t ϱ x d x G ϱ x 2 ( 1 ) + ϱ x 2 ( 0 ) + 3 G 0 1 ϱ x 2 d x + G 0 1 ( 3 s ξ ) x 2 d x + 2 ρ 0 1 ϱ t 2 d x + 0 1 ϱ x 2 d x + μ 2 0 1 z 2 ( x , 1 , t ) d x .
Notice that
0 1 ( 3 s ξ ) x 2 d x 2 0 1 ξ x 2 d x + 18 0 1 s x 2 d x ,
and
0 1 ϱ x 2 d x = 0 1 3 s ξ ϱ x ( 3 s ξ ) 2 d x 2 0 1 3 s ξ ϱ x 2 d x + 2 0 1 ( 3 s ξ ) x 2 d x 2 0 1 3 s ξ ϱ x 2 d x + 4 0 1 ξ x 2 d x + 36 0 1 s x 2 d x .
Therefore, we conclude from (55) that
ϱ x 2 ( 1 ) + ϱ x 2 ( 0 ) 1 G d d t 0 1 ρ g ( x ) ϱ t ϱ x d x + 2 G + 6 0 1 3 s ξ ϱ x 2 d x + μ 2 G 0 1 z 2 ( x , 1 , t ) d x + 4 G + 14 0 1 ξ x 2 d x + 36 G + 126 0 1 s x 2 d x + 2 ρ G 0 1 ϱ t 2 d x .
Combining (54) and (56) with (50), we have for any ε 4 > 0 and ε 5 > 0 that
D ξ x ϱ x | x = 0 x = 1 ( 6 + 2 G ) ε 4 + D 4 ε 4 ε 5 0 1 3 s ξ ϱ x 2 d x + 36 G + 126 ε 4 0 1 s x 2 d x + 2 ρ G ε 4 0 1 ϱ t 2 d x + I ρ D 2 ε 4 0 1 ξ t 2 d x + γ 1 2 2 ε 4 0 1 θ x 2 d x + γ 2 2 2 ε 4 0 L P x 2 d x + 4 G + 14 ε 4 + D 4 ε 4 3 D + G 2 ε 5 0 1 ξ x 2 d x + μ 2 ε 4 G 0 1 z 2 ( x , 1 , t ) d x ε 4 G d d t 0 1 ρ g ( x ) ϱ t ϱ x d x D 4 ε 4 d d t 0 1 I ρ g ( x ) ξ t ξ x d x .
Inserting (57) into (49) and using (5), then (48) follows. □
Lemma 11.
Let the functional H 3 be defined by
H 3 ( t ) = 3 I ρ 0 1 s s t d x + 2 β 0 1 s 2 d x .
Then, H 3 satisfies
d d t H 3 ( t ) 3 D 0 1 s x 2 d x 3 γ 0 1 s 2 d x + 3 I ρ 0 1 s t 2 d x + 9 G 2 4 γ 0 1 3 s ξ ϱ x 2 d x .
Proof. 
Clearly,
d d t H 3 ( t ) = 3 I ρ 0 1 s t 2 d x 3 G 0 1 3 s ξ ϱ x s d x 3 D 0 1 s x 2 d x 4 γ 0 1 s 2 d x .
Afterward, using Young’s inequality, we arrive at (58). □
In the following, we define the functional L ^ by
L ^ ( t ) = M E ( t ) + M 1 F 1 ( t ) + H 1 ( t ) + M 3 H 2 ( t ) + M 4 F 3 ( t ) + M 5 H 3 ( t ) + M 6 F 6 + ε 4 G M 4 0 1 ρ g ( x ) ϱ t ϱ x d x + D 4 ε 4 M 4 0 1 I p g ( x ) ξ t ξ x d x ,
where M and M i , i = 1 , 3 , 4 , 5 , 6 will be determined later. As in (40), there exist two positive constants ν ^ 1 and ν ^ 2 such that
ν ^ 1 E ( t ) L ^ ( t ) ν ^ 2 E ( t ) .
Proof. 
(Of Theorem 2 for Boundary I I )
Combining (22), (42), (46), (48), and (58), as well as taking
ε 1 = 1 M 1 , ε 2 = 1 D M 1 , ε 3 = ρ 4 M 3 , ε 4 = G 8 M 4 , ε 5 = 1 M 4 2 , M 6 = | μ | ,
we have
d d t L ^ ( t ) κ M c 1 M 1 2 γ 1 2 D M 3 γ 1 2 G M 4 4 γ 1 2 G M 4 2 0 1 θ x 2 d x M c 2 M 1 2 γ 2 2 D M 3 γ 2 2 G M 4 4 γ 2 2 G M 4 2 0 1 P x 2 d x 4 β M 3 I ρ 4 M 4 3 I ρ M 5 0 1 s t 2 d x ρ 2 | μ | M + 1 0 1 ϱ t 2 d x I ρ 2 M 1 I ρ + ρ M 3 M 3 2 I ρ M 4 4 D I ρ G M 4 2 0 1 ξ t 2 d x D 2 M 3 5 7 G 4 2 D G M 4 2 3 D + G 2 M 4 2 0 1 ξ x 2 d x G 2 M 4 9 G 4 5 4 2 D G G 2 4 9 G 2 4 γ M 5 0 1 3 s ξ ϱ x 2 d x 3 D M 5 81 2 63 G 4 0 1 s x 2 d x 3 γ M 5 0 1 s 2 d x m 1 | μ | τ 0 1 0 1 e σ τ z 2 ( x , σ , t ) d σ d x | μ | 2 m 1 | μ | M 3 + M 4 + 9 8 0 1 z 2 ( x , 1 , t ) d x .
We first take M 5 > 0 such that
3 D M 5 81 2 63 G 4 > 0 .
For fixed M 5 , we choose M 4 > 0 large enough such that
G 2 M 4 9 G 4 5 4 2 D G G 2 4 9 G 2 4 γ M 5 > 0 .
For fixed M 4 , we choose M 3 > 0 large enough such that
D 2 M 3 5 7 G 4 2 D G M 4 2 3 D + G 2 M 4 2 > 0 .
For fixed M 4 and M 3 , we choose M 1 > 0 large enough such that
I ρ 2 M 1 I ρ + ρ M 3 M 3 2 I ρ M 4 4 D I ρ G M 4 2 > 0 .
For fixed M i , i = 1 , 3 , 4 , 5 , we at last take M large enough such that
κ M c 1 M 1 2 γ 1 2 D M 3 γ 1 2 G M 4 4 γ 1 2 G M 4 2 > 0 ,
M c 2 M 1 2 γ 2 2 D M 3 γ 2 2 G M 4 4 γ 2 2 G M 4 2 > 0 ,
and
4 β M 3 I ρ 4 M 4 3 I ρ M 5 > 0 .
Finally, we pick
μ ¯ = min m 1 M 3 + M 4 + 9 8 , ρ 2 ( M + 1 ) .
We find that there exists a positive constant ν ^ 0 so that
d d t L ^ ( t ) ν ^ 0 E ( t ) ,
together with (59), gives
E ( t ) ν ^ 2 ν ^ 1 E ( 0 ) e ν ^ 0 ν ^ 2 t .

4. Conclusions

In this work, we have investigated the existence of a solution for our system, which is discussed within the context of the semigroup approach. In addition, under different boundary conditions, two results of stability properties independent of initial data are studied (see [3,9,10,11,12,13]). In future work, we will follow the same study but with a time-varying delay term, and the numerical examples will be investigated.

Author Contributions

Z.K. and A.K.: writing original draft, Methodology, Resources, Methodology, formal analysis, Conceptualization; S.B. and A.H.A.: conceptualized, investigated, analyzed and validated the research while, I.L. and S.B.: formulated, investigated, reviewed, S.B.: Corresponding author, D.O. and S.B.: Supervision. All authors have read and agreed to the published version of the manuscript.

Funding

The researchers would like to thank the Deanship of Graduate Studies and Scientific Research at Qassim University for financial support (QU-APC-2024-9/1).

Institutional Review Board Statement

There is no ethical issue in this work. All the authors actively participated in this research and approved it for publication.

Data Availability Statement

There are no data associated with the current study.

Acknowledgments

The researchers would like to thank the Deanship of Graduate Studies and Scientific Research at Qassim University for financial support (QU-APC-2024-9/1).

Conflicts of Interest

There are no competing interests regarding this research work.

References

  1. Apalara, T.A. On the stability of a thermoelastic laminated beam. Acta Math. Sci. 2019, 39B, 1517–1524. [Google Scholar] [CrossRef]
  2. Aouadi, M.; Campo, M.; Copetti, M.I.M.; Fernandez, J.R. Existence, stability and numerical results for a Timoshenko beam with thermodiffusion effects. Z. Angew. Math. Phys. 2019, 70, 117. [Google Scholar] [CrossRef]
  3. Khalili, Z.; Ouchenane, D.; El-hamidi, A. Stability results for Laminated beam with Thermo-visco-Elastic effects and localized nonlinear damping. Nonlinear Anal. Real World Appl. 2022, 67, 103555. [Google Scholar] [CrossRef]
  4. Khalili, Z.; Ouchenane, D. A stability result for a Timoshenko system with infinite history and distributed delay term. Kragujev. J. Math. 2023, 47, 281–295. [Google Scholar] [CrossRef]
  5. Khalili, Z.; Ouchenane, D.; Choucha, A. Well-posedness and Stability result of a nonlinear damping Porous-elastic system in thermoelasticity of second sound with infinite memory and distributed delay terms. Dyn. Contin. Discret. Impuls. Syst. Ser. Math. Anal. 2023, 30, 21–53. [Google Scholar]
  6. Khalili, Z.; Ouchenane, D. Stability results for Laminated beam with thermodiffusion effects and distributed delay. Stud. Eng. Exact Sci. 2024, 5, e7725. [Google Scholar]
  7. Khanchehgardan, A.; Rezazadeh, G.; Shabani, R. Effect of mass diffusion on the damping ratio in micro-beam resonators. Int. J. Solids Struct. 2014, 51, 3147–3155. [Google Scholar] [CrossRef]
  8. Pazy, A. Semigroups of Linear Operators and Applications to Partial Differential Equations; Applied Mathematical Sciences; Springer: New York, NY, USA, 1983; Volume 44. [Google Scholar]
  9. Arnold, V.I. Mathematical Methods of Classical Mechanics; Springer: New York, NY, USA, 1989. [Google Scholar]
  10. Hansen, S.W.; Spies, R. Structural damping in a laminated beams due to interfacial slip. J. Sound Vib. 1997, 204, 183–202. [Google Scholar] [CrossRef]
  11. Li, G.; Kong, X.; Wenjun, L. General decay for a laminated beam with structural damping and memory, the case of non-equal-wave speeds. J. Integral Equ. Appl. 2018, 30, 95–116. [Google Scholar] [CrossRef]
  12. Mustafa, M.I. General decay result for nonlinear viscoelastic equations. J. Math. Anal. Appl. 2018, 457, 134–152. [Google Scholar] [CrossRef]
  13. Mustafa, M.I. Laminated Timoshenko beams with viscoelastic damping. J. Math. Anal. Appl. 2018, 466, 619–641. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Khalili, Z.; Ouchenane, D.; Krelifa, A.; Laribi, I.; Boulaaras, S.; Ahmed, A.H. Existence and Stability Results for Thermodiffusion Laminated Beam System with Delay Feedback. Mathematics 2024, 12, 3097. https://doi.org/10.3390/math12193097

AMA Style

Khalili Z, Ouchenane D, Krelifa A, Laribi I, Boulaaras S, Ahmed AH. Existence and Stability Results for Thermodiffusion Laminated Beam System with Delay Feedback. Mathematics. 2024; 12(19):3097. https://doi.org/10.3390/math12193097

Chicago/Turabian Style

Khalili, Zineb, Djamel Ouchenane, Ali Krelifa, Imene Laribi, Salah Boulaaras, and Ahmed Himadan Ahmed. 2024. "Existence and Stability Results for Thermodiffusion Laminated Beam System with Delay Feedback" Mathematics 12, no. 19: 3097. https://doi.org/10.3390/math12193097

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