Next Article in Journal
Statistical Modeling of High Frequency Datasets Using the ARIMA-ANN Hybrid
Previous Article in Journal
Peak-to-Peak Stabilization of Sampled-Data Systems Subject to Actuator Saturation and Its Practical Application to an Inverted Pendulum
Previous Article in Special Issue
Lump-Type Solutions, Lump Solutions, and Mixed Rogue Waves for Coupled Nonlinear Generalized Zakharov Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

General Stability for the Viscoelastic Wave Equation with Nonlinear Time-Varying Delay, Nonlinear Damping and Acoustic Boundary Conditions

1
Department of Mathematics, Pusan National University, Busan 46241, Republic of Korea
2
Department of Applied Mathematics, Pukyong National University, Busan 48513, Republic of Korea
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(22), 4593; https://doi.org/10.3390/math11224593
Submission received: 4 October 2023 / Revised: 26 October 2023 / Accepted: 6 November 2023 / Published: 9 November 2023
(This article belongs to the Special Issue Applications of Partial Differential Equations)

Abstract

:
This paper is focused on energy decay rates for the viscoelastic wave equation that includes nonlinear time-varying delay, nonlinear damping at the boundary, and acoustic boundary conditions. We derive general decay rate results without requiring the condition a 2 > 0 and without imposing any restrictive growth assumption on the damping term f 1 , using the multiplier method and some properties of the convex functions. Here we investigate the relaxation function ψ , namely ψ ( t ) μ ( t ) G ( ψ ( t ) ) , where G is a convex and increasing function near the origin, and μ is a positive nonincreasing function. Moreover, the energy decay rates depend on the functions μ and G, as well as the function F defined by f 0 , which characterizes the growth behavior of f 1 at the origin.

1. Introduction

In this paper, we study the energy decay rates for the viscoelastic wave equation with nonlinear time-varying delay, nonlinear damping at the boundary, and acoustic boundary conditions
u t t ( x , t ) Δ u ( x , t ) + 0 t ψ ( t s ) Δ u ( x , s ) d s = 0 , in Ω × ( 0 , ) ,
u ( x , t ) = 0 , on Γ 0 × ( 0 , ) ,
u ν ( x , t ) 0 t ψ ( t s ) u ν ( x , s ) d s + a 1 f 1 ( u t ( x , t ) ) + a 2 f 2 ( u t ( x , t ϱ ( t ) ) ) = w t ( x , t ) , on Γ 1 × ( 0 , ) ,
u t ( x , t ) + h ( x ) w t ( x , t ) + m ( x ) w ( x , t ) = 0 , on Γ 1 × ( 0 , ) ,
u ( x , 0 ) = u 0 ( x ) , u t ( x , 0 ) = u 1 ( x ) , in Ω ,
u t ( x , t ) = j 0 ( x , t ) , in Γ 1 × ( ϱ ( 0 ) , 0 ) ,
where Ω is a bounded domain in R n ( n 1 ) with smooth boundary Γ of class C 2 ; Γ = Γ 0 Γ 1 , where Γ 0 and Γ 1 are closed and disjoint; w ( x , t ) is the normal displacement into the domain of a point x Γ 1 at time t; and h , m : Γ 1 R are essential bounded functions that represent resistivity and spring constant per unit area, respectively. f 1 , f 2 : R R are given functions, and f 1 represents the nonlinear frictional damping. a 1 , a 2 are real numbers with a 1 > 0 , a 2 0 . The integral term is the memory responsible for the viscoelastic damping. The functions ψ and ϱ ( t ) represent the kernel of the memory term and the time-varying delay, respectively. ν is the outward unit normal vector to Γ . The initial data ( u 0 , u 1 , j 0 ) belong to a suitable space. Boundary conditions (3) and (4) are called acoustic boundary conditions.
In the past decades, the non-delayed wave equation with a viscoelastic term has garnered significant attention in the field of partial differential equations. Research on the energy decay rate of the solution to the viscoelastic wave equation is vital in various fields, contributing to technological advancements, safety assurance, environmental protection, energy efficiency, and academic exploration. The stability of solutions for such equations has recently been studied by many authors (see [1,2,3] and references therein). When a 1 = a 2 = 0 , models (1)–(5) are pertinent to noise control and suppression in practical applications. The noise propagates through some acoustic medium, like air, in a room that is defined by a bounded domain Ω and whose floor, walls, and ceiling are determined by the boundary conditions [4,5]. Under the conditions that 0 ψ ( s ) d s < 1 2 and ψ ( t ) μ ( t ) ψ ( t ) , for t 0 , Park and Park [6] considered the general decay for problems (1)–(5). Liu [7] improved the research of [6] by achieving arbitrary rates of decay, which may not necessarily be an exponential or a polynomial one. Recently, Yoon et al. [8] generalized the work of [6,7] without the assumption condition 0 ψ ( s ) d s < 1 2 . The assumption on relaxation function ψ has been weakened compared to the conditions assumed in previous literature [6,7].
Numerous phenomena are influenced by both the current state and the previous occurrences of the system. There has been a notable increase in the research on the equation with delay effects, which frequently arise in various physical, biological, chemical, medical, and economic problems [9,10,11]. However, the delay effects can generally be considered a cause of instability. In order to stabilize a system containing delay terms, additional control terms will be necessary. Kirane and Said-Houari [12] showed the global existence and asymptotic stability for the following wave equation with memory and constant delay,
u t t ( x , t ) Δ u ( x , t ) + 0 t ψ ( t s ) Δ u ( x , s ) d s + a 1 u t ( x , t ) + a 2 u t ( x , t ϱ ) = 0 ,
where a 1 , a 2 , and ϱ are positive constants. They used the damping term a 1 u t ( x , t ) to control the delay term in obtaining the decay estimate of the energy. They proved that its energy was exponentially decaying when a 2 a 1 . Dai and Yang [13] investigated the exponential decay of an unsolved problem proposed by Kirane and Said-Houari [12], namely, the problem with a 1 = 0 . In the case of constant weight and constant delay, the delay term typically considers the past history of strain, only up to some finite time ϱ ( t ) ϱ . Nicaise and Pignotti [14] investigated the following wave equation with internal time-varying delay instead of constant delay,
u t t ( x , t ) Δ u ( x , t ) + a 1 u t ( x , t ) + a 2 u t ( x , t ϱ ( t ) ) = 0 ,
where ϱ ( t ) > 0 ,   a 1 , and a 2 are real numbers with a 1 > 0 . They proved the exponential stability result for the wave equation under the condition | a 2 | < 1 ζ 0 a 1 , where the constant ζ 0 satisfies ϱ ( t ) ζ 0 < 1 , t > 0 . Liu [15] studied the following wave equation involving memory and time-varying delay:
u t t ( x , t ) Δ u ( x , t ) + α ( t ) 0 t ψ ( t s ) Δ u ( x , s ) d s + a 1 u t ( x , t ) + a 2 u t ( x , t ϱ ( t ) ) = 0 .
Systems with time-varying delays have been extensively considered by many authors (see [16,17,18,19,20,21,22] and references therein). Recently, Zennir [23] considered the stability for solutions of plate equations with a time-varying delay and weak viscoelasticity in R n . Moreover, Benaissa et al. [24] proved the global existence and stability for solutions of the following wave equation with a time-varying delay in the weakly nonlinear feedback,
u t t ( x , t ) Δ u ( x , t ) + a 1 σ ( t ) f 1 ( u t ( x , t ) ) + a 2 σ ( t ) f 2 ( u t ( x , t ϱ ( t ) ) ) = 0 ,
where ϱ ( t ) > 0 ,   a 1 , and a 2 are positive real numbers, and f 1 , f 2 satisfy some conditions. This result extended the previous work [10,14]. Park [25] investigated the decay result of the energy for a von Karman equation with time-varying delay by dropping the restriction a 2 > 0 under the same conditions as ϱ ,   f 1 , and f 2 in [24]. For the viscoelastic problem with time-varying delay in the nonlinear internal or boundary feedback, we also refer to [26,27]. As far as we know, there are few results for the viscoelastic wave equation with a nonlinear time-varying delay. Recently, Djeradi et al. [28] and Mukiawa et al. [29] showed the stability of the thermoelastic laminated beam and thermoelastic Timoshenko beam with nonlinear time-varying delay, respectively. The papers introduced so far have studied the energy decay rate of the solution for the equation with nonlinear time-varying delay in the Dirichlet boundary condition.
Motivated by these results, we study the general decay rates of the solution for problems (1)–(6) with a nonlinear time-varying delay term, nonlinear damping at the boundary, and acoustic boundary conditions. Research on the energy decay rate of solutions for the viscoelastic wave equation with nonlinear time-delay terms plays a critical role in various application areas, including stability assessment, understanding complex behaviors, advancing neuroscience, disaster preparedness, and improving energy efficiency. We consider the general assumption on the relaxation function ψ ,
ψ ( t ) μ ( t ) G ( ψ ( t ) ) ,
where μ : R + R + is a positive nonincreasing function, and G is linear or is a strictly increasing and strictly convex function. We derive the general decay rate results without requiring the condition a 2 > 0 and without imposing any restrictive growth assumption on the damping term f 1 . The energy decay rates depend on the functions μ and G , as well as the function F defined by f 0 , which represents the growth f 1 at the origin. Our result improves upon previous work [6,7,8].
This paper is composed of the following. In Section 2, we prepare some notations and materials needed for our work. In Section 3, we introduce some technical lemmas to prove our stability result. In Section 4, we state and prove the general energy decay.

2. Preliminaries

In this section, we present some materials required for our results. Throughout this paper, we use the notation
V = { u H 1 ( Ω ) : u = 0 on Γ 0 } .
For simplicity, we denote · L 2 ( Ω ) and · L 2 ( Γ 1 ) by · and · Γ 1 , respectively.
The Poincaré inequality holds in V; that is, there exist the positive constants λ 0 and λ 1 such that
u 2 λ 0 u 2 and u Γ 1 2 λ 1 u 2 for all u V .
As in [1,3,8,26,30], we consider the following assumptions for ψ , f 1 , f 2 , ϱ , h , and m .
(H1) ψ : [ 0 , ) R + is a differentiable function satisfying
1 0 ψ ( s ) d s = l > 0 ,
and there exists a C 1 function G : R + R + that is linear or is a strictly convex and strictly increasing C 2 function on ( 0 , r 0 ] , r 0 ψ ( 0 ) such that
ψ ( t ) μ ( t ) G ( ψ ( t ) ) , t 0 ,
where G ( 0 ) = G ( 0 ) = 0 , and μ is a positive nonincreasing differentiable function. The function G was first introduced in [31]. These are weaker conditions on G than those introduced in [31].
(H2) f 1 : R R is a nondecreasing C 0 function such that there exists a strictly increasing function f 0 C 1 ( R + ) , with f 0 ( 0 ) = 0 , and positive constants c 0 , c 1 , and ε such that
f 0 ( | s | ) | f 1 ( s ) | f 0 1 ( | s | ) for all | s | ε ,
c 0 | s | | f 1 ( s ) | c 1 | s | for all | s | ε .
Moreover, we assume that the function F, defined by F ( s ) = s f 0 ( s ) , is a strictly convex C 2 function on ( 0 , r 1 ] , for some r 1 > 0 , when f 0 is nonlinear.
(H3) f 2 : R R is an odd nondecreasing C 1 function such that there exist positive constants c 2 , c 3 , and c 4 that satisfy
| f 2 ( s ) | c 2 , c 3 s f 2 ( s ) F 2 ( s ) c 4 s f 1 ( s ) , for s R ,
where F 2 ( s ) = 0 s f 2 ( t ) d t .
(H4) ϱ W 2 , ( [ 0 , T ] ) is a function such that
0 < ϱ 1 ϱ ( t ) ϱ 2 and ϱ ( t ) ϱ 3 < 1 for all t > 0 ,
where T ,   ϱ 1 , and ϱ 2 are positive constants. Moreover, the weight of dissipation and the delay satisfy
0 < | a 2 | < c 3 ( 1 ϱ 3 ) c 4 ( 1 c 3 ϱ 3 ) a 1 .
(H5) We assume that h , m C ( Γ 1 ) , h ( x ) > 0 , and m ( x ) > 0 for all x Γ 1 . Then, there exist positive constants h i and m i ( i = 1 , 2 ) such that
h 1 h ( x ) h 2 , m 1 m ( x ) m 2 for all x Γ 1 .
Remark 1. 
1. The assumption (H2) implies that s f 1 ( s ) > 0 , for all s 0 .
2. The assumption (11) of function f 1 has been weakened compared to the condition assumed in [24,25].
3. Since f 2 is an odd nondecreasing function, F 2 is an even and convex function. Furthermore, it is satisfied that F 2 ( s ) = 0 s f 2 ( t ) d t s f 2 ( s ) . From (13), we find that c 3 1 .
Remark 2 
([3]). 1. By (H1), we obtain lim t + ψ ( t ) = 0 . Then, there exists t 0 0 large enough that
ψ ( t 0 ) = r 0 ψ ( t ) r 0 , t t 0 .
Given ψ and μ are positive nonincreasing continuous functions, G is a positive continuous function, and for (10), we have, for some positive constant c 5 ,
ψ ( t ) μ ( t ) G ( ψ ( t ) ) c 5 ψ ( t ) , t [ 0 , t 0 ] .
2. If G is a strictly convex and strictly increasing C 2 function on ( 0 , r 0 ] , with G ( 0 ) = G ( 0 ) = 0 , then it has an extension G ¯ , which is a strictly convex and strictly increasing C 2 function on ( 0 , ) . The same remark can be established for F ¯ .
We recall the well-known Jensen inequality, which plays a pivotal role in proving our main result. If ϕ is a convex function on [ a , b ] , p : Ω [ a , b ] and k represents integrable functions on Ω such that k ( x ) 0 and Ω k ( x ) d x = k 0 > 0 , then Jensen’s inequality holds:
ϕ 1 k 0 Ω p ( x ) k ( x ) d x 1 k 0 Ω ϕ [ p ( x ) ] k ( x ) d x .
Let H * be the conjugate of the convex function H defined by H * ( s ) = sup r 0 ( s r H ( r ) ) , then
s r H * ( s ) + H ( r ) , s , r 0 .
Moreover, due to the argument provided in [32], it holds that
H * ( s ) = s ( H ) 1 ( s ) H ( H ) 1 ( s ) , s 0 .
As in [10,14], we introduce the following new function:
v ( x , κ , t ) = u t ( x , t κ ϱ ( t ) ) , for ( x , κ , t ) Γ 1 × ( 0 , 1 ) × ( 0 , ) .
Then, problems (1)–(6) can be expressed as follows:
u t t ( x , t ) Δ u ( x , t ) + 0 t ψ ( t s ) Δ u ( x , s ) d s = 0 , in Ω × ( 0 , ) ,
ϱ ( t ) v t ( x , κ , t ) + ( 1 κ ϱ ( t ) ) v κ ( x , κ , t ) = 0 , in Γ 1 × ( 0 , 1 ) × ( 0 , ) ,
u ( x , t ) = 0 , in Γ 0 × ( 0 , ) ,
u ν ( x , t ) 0 t ψ ( t s ) u ν ( x , s ) d s + a 1 f 1 ( u t ( x , t ) ) + a 2 f 2 ( v ( x , 1 , t ) ) = w t ( x , t ) , on Γ 1 × ( 0 , ) ,
u t ( x , t ) + h ( x ) w t ( x , t ) + m ( x ) w ( x , t ) = 0 , on Γ 1 × ( 0 , ) ,
u ( x , 0 ) = u 0 ( x ) , u t ( x , 0 ) = u 1 ( x ) , in Ω ,
v ( x , κ , 0 ) = j 0 ( x , κ ϱ ( 0 ) ) , in Γ 1 × ( 0 , 1 ) .
We state the global existence result that can be established by the arguments of [24,33].
Theorem 1. 
Let initial data ( u 0 , u 1 ) ( V H 2 ( Ω ) ) × V and j 0 L 2 ( Γ 1 × ( 0 , 1 ) ) . Suppose that (H1)–(H5) hold. Then, for any T > 0 , there exists a unique pair of functions ( u , w , v ) that are the solution to problems (22)–(28) in the class
u L ( 0 , T ; V H 2 ( Ω ) ) , u t L ( 0 , T ; V ) , u t t L ( 0 , T ; L 2 ( Ω ) ) , v L ( 0 , T ; L 2 ( Γ 1 × ( 0 , 1 ) ) ) , w , w t L 2 ( 0 , ; L 2 ( Γ 1 ) ) .
As in [6,25], we introduce the energy for problems (22)–(28),
E ( t ) = 1 2 u t ( t ) 2 + 1 2 1 0 t ψ ( s ) d s u ( t ) 2 + 1 2 ( ψ u ) ( t ) + 1 2 Γ 1 m ( x ) w 2 ( t ) d Γ + ζ ϱ ( t ) 2 Γ 1 0 1 F 2 ( v ( x , κ , t ) ) d κ d Γ ,
where
( ψ u ) ( t ) = 0 t ψ ( t s ) u ( t ) u ( s ) 2 d s
and
2 | a 2 | ( 1 c 3 ) c 3 ( 1 ϱ 3 ) < ζ < 2 ( a 1 | a 2 | c 4 ) c 4 .
Thanks to (15), this makes sense.
To show the main results of this paper, we need the following lemma.
Lemma 1. 
Assume that (H3)–(H5) hold. Then, there exist positive constants γ 0 and γ 1 satisfying
E ( t ) 1 2 ( ψ u ) ( t ) 1 2 ψ ( t ) u ( t ) 2 h 1 | | w t ( t ) | | Γ 1 2 γ 0 Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ γ 1 Γ 1 f 2 ( v ( x , 1 , t ) ) v ( x , 1 , t ) d Γ .
Proof. 
Multiplying by u t ( t ) in (22), using Green’s formula, (25), and (26), we have
1 2 d d t u t ( t ) 2 + 1 0 t ψ ( s ) d s u ( t ) 2 + ( ψ u ) ( t ) + Γ 1 m ( x ) w 2 ( t ) d Γ = 1 2 ( ψ u ) ( t ) 1 2 ψ ( t ) u ( t ) 2 Γ 1 h ( x ) w t 2 ( t ) d Γ a 1 Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ a 2 Γ 1 f 2 ( v ( x , 1 , t ) ) u t ( t ) d Γ ,
where we used the relation
Ω u t ( t ) 0 t ψ ( t s ) u ( s ) d s d x = d d t 1 2 ( ψ u ) ( t ) 1 2 0 t ψ ( s ) d s u ( t ) 2 1 2 ( ψ u ) ( t ) + 1 2 ψ ( t ) u ( t ) 2 .
From (29) and (32), we have
E ( t ) = 1 2 ( ψ u ) ( t ) 1 2 ψ ( t ) u ( t ) 2 Γ 1 h ( x ) w t 2 ( t ) d Γ a 1 Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ a 2 Γ 1 f 2 ( v ( x , 1 , t ) ) u t ( t ) d Γ + ζ ϱ ( t ) 2 Γ 1 0 1 F 2 ( v ( x , κ , t ) ) d κ d Γ + ζ ϱ ( t ) 2 Γ 1 0 1 f 2 ( v ( x , κ , t ) ) v t ( x , κ , t ) d κ d Γ ,
where F 2 ( t ) = 0 t f 2 ( s ) d s . In (23), we multiply by f 2 ( v ( x , κ , t ) ) and integrate over Γ 1 × ( 0 , 1 ) to obtain
ζ ϱ ( t ) 2 Γ 1 0 1 f 2 ( v ( x , κ , t ) ) v t ( x , κ , t ) d κ d Γ = ζ 2 Γ 1 ( 1 ϱ ( t ) ) F 2 ( v ( x , 1 , t ) ) F 2 ( v ( x , 0 , t ) ) + 0 1 ϱ ( t ) F 2 ( v ( x , κ , t ) ) d κ d Γ .
Applying this to (33) and noting that v ( x , 0 , t ) = u t ( x , t ) , it follows that
E ( t ) = 1 2 ( ψ u ) ( t ) 1 2 ψ ( t ) u ( t ) 2 Γ 1 h ( x ) w t 2 ( t ) d Γ a 1 Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ a 2 Γ 1 f 2 ( v ( x , 1 , t ) ) u t ( t ) d Γ ζ 2 Γ 1 ( 1 ϱ ( t ) ) F 2 ( v ( x , 1 , t ) ) F 2 ( u t ( x , t ) ) d Γ .
From (13) and (14), we obtain
ζ 2 Γ 1 ( 1 ϱ ( t ) ) F 2 ( v ( x , 1 , t ) ) F 2 ( u t ( x , t ) ) d Γ ζ c 3 2 ( 1 ϱ 3 ) Γ 1 f 2 ( v ( x , 1 , t ) ) v ( x , 1 , t ) d Γ + ζ c 4 2 Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ .
Substituting (35) into (34), we obtain
E ( t ) 1 2 ( ψ u ) ( t ) 1 2 ψ ( t ) u ( t ) 2 Γ 1 h ( x ) w t 2 ( t ) d Γ a 1 ζ c 4 2 Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ ζ c 3 2 ( 1 ϱ 3 ) Γ 1 f 2 ( v ( x , 1 , t ) ) v ( x , 1 , t ) d Γ a 2 Γ 1 f 2 ( v ( x , 1 , t ) ) u t ( t ) d Γ .
Now, we estimate the last term in the right-hand side of (36). The definition of F 2 and (21) give
F 2 * ( s ) = s f 2 1 ( s ) F 2 ( f 2 1 ( s ) ) , for s 0 .
When f 2 ( v ( x , 1 , t ) ) < 0 and u t ( t ) 0 , using (20) and (37) with s = f 2 ( v ( x , 1 , t ) ) and r = u t ( t ) , we obtain (see details in [25])
a 2 Γ 1 ( f 2 ( v ( x , 1 , t ) ) ) u t ( t ) d Γ | a 2 | Γ 1 f 2 ( v ( x , 1 , t ) ) ( v ( x , 1 , t ) ) F 2 ( v ( x , 1 , t ) ) + F 2 ( u t ( t ) ) d Γ = | a 2 | Γ 1 f 2 ( v ( x , 1 , t ) ) v ( x , 1 , t ) F 2 ( v ( x , 1 , t ) ) + F 2 ( u t ( t ) ) d Γ ,
where we used the fact that f 2 is odd and F 2 is even. When f 2 ( v ( x , 1 , t ) ) 0 and u t ( t ) < 0 , with s = f 2 ( v ( x , 1 , t ) ) and r = u t ( t ) , we obtain
a 2 Γ 1 f 2 ( v ( x , 1 , t ) ) ( u t ( t ) ) d Γ | a 2 | Γ 1 f 2 ( v ( x , 1 , t ) ) ( v ( x , 1 , t ) ) F 2 ( v ( x , 1 , t ) ) + F 2 ( u t ( t ) ) d Γ = | a 2 | Γ 1 f 2 ( v ( x , 1 , t ) ) v ( x , 1 , t ) F 2 ( v ( x , 1 , t ) ) + F 2 ( u t ( t ) ) d Γ .
From (38) and (39), for the case f 2 ( v ( x , 1 , t ) ) u t ( t ) 0 , we have
a 2 Γ 1 f 2 ( v ( x , 1 , t ) ) u t ( t ) d Γ | a 2 | Γ 1 f 2 ( v ( x , 1 , t ) ) v ( x , 1 , t ) F 2 ( v ( x , 1 , t ) ) + F 2 ( u t ( t ) ) d Γ .
Similarly, (40) holds when f 2 ( v ( x , 1 , t ) ) u t ( t ) 0 . Hence, using (13) and (40), we see that
a 2 Γ 1 f 2 ( v ( x , 1 , t ) ) u t ( t ) d Γ | a 2 | ( 1 c 3 ) Γ 1 f 2 ( v ( x , 1 , t ) ) v ( x , 1 , t ) d Γ + c 4 Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ .
By using (16), (36), and (41), and by selecting ζ satisfying (30), we obtain the desired inequality (31) where γ 0 = a 1 ζ c 4 2 | a 2 | c 4 > 0 and γ 1 = ζ c 3 2 ( 1 ϱ 3 ) | a 2 | ( 1 c 3 ) > 0 .

3. Technical Lemmas

In this section, we prove the following lemmas to obtain the general decay rates of the solution to problems (22)–(28).
Lemma 2. 
Under the assumption (H1), the functional Φ 1 defined by
Φ 1 ( t ) = Ω u ( t ) u t ( t ) d x + Γ 1 u ( t ) w ( t ) d Γ + 1 2 Γ 1 h ( x ) w 2 ( t ) d Γ
satisfies
Φ 1 ( t ) u t ( t ) 2 l 2 u ( t ) 2 + 2 C ( ξ ) l ( i u ) ( t ) + 8 λ 1 l w t ( t ) Γ 1 2 + a 1 a 3 l Γ 1 f 1 2 ( u t ( t ) ) d Γ + | a 2 | a 3 l Γ 1 f 2 2 ( v ( x , 1 , t ) ) d Γ Γ 1 m ( x ) w 2 ( t ) d Γ ,
for any 0 < ξ < 1 , where
i ( t ) = ξ ψ ( t ) ψ ( t ) and C ( ξ ) = 0 ψ 2 ( s ) i ( s ) d s .
Proof. 
Using Equations (22) and (24)–(26), and utilizing (9) and Young’s inequality, we obtain
Φ 1 ( t ) = u t ( t ) 2 1 0 t ψ ( s ) d s u ( t ) 2 + 0 t ψ ( t s ) ( u ( s ) u ( t ) , u ( t ) ) d s a 1 Γ 1 f 1 ( u t ( t ) ) u ( t ) d Γ a 2 Γ 1 f 2 ( v ( x , 1 , t ) ) u ( t ) d Γ + 2 Γ 1 u ( t ) w t ( t ) d Γ Γ 1 m ( x ) w 2 ( t ) d Γ u t ( t ) 2 7 l 8 u ( t ) 2 + 2 l Ω 0 t ψ ( t s ) | u ( s ) u ( t ) | d s 2 d x a 1 Γ 1 f 1 ( u t ( t ) ) u ( t ) d Γ a 2 Γ 1 f 2 ( v ( x , 1 , t ) ) u ( t ) d Γ + 2 Γ 1 u ( t ) w t ( t ) d Γ Γ 1 m ( x ) w 2 ( t ) d Γ .
Using the Cauchy–Schwarz inequality and (43), we have (see [3,34])
Ω 0 t ψ ( t s ) | u ( s ) u ( t ) | d s 2 d x 0 t ψ 2 ( s ) i ( s ) d s ( i u ) ( t ) C ( ξ ) ( i u ) ( t ) .
Applying Young’s inequality and (8), we obtain, for η > 0 ,
| a 1 Γ 1 f 1 ( u t ( t ) ) u ( t ) d Γ | η a 1 λ 1 u ( t ) 2 + a 1 4 η Γ 1 f 1 2 ( u t ( t ) ) d Γ ,
| a 2 Γ 1 f 2 ( v ( x , 1 , t ) ) u ( t ) d Γ | η | a 2 | λ 1 u ( t ) 2 + | a 2 | 4 η Γ 1 f 2 2 ( v ( x , 1 , t ) ) d Γ ,
and
2 Γ 1 u ( t ) w t ( t ) d Γ l 8 u ( t ) 2 + 8 λ 1 l w t ( t ) Γ 1 2 .
Combining estimates (44)–(47), we see that
Φ 1 ( t ) u t ( t ) 2 ( 3 l 4 η a 1 λ 1 η | a 2 | λ 1 ) u ( t ) 2 + 2 C ( ξ ) l ( i u ) ( t ) + 8 λ 1 l w t ( t ) Γ 1 2 + a 1 4 η Γ 1 f 1 2 ( u t ( t ) ) d Γ + | a 2 | 4 η Γ 1 f 2 2 ( v ( x , 1 , t ) ) d Γ Γ 1 m ( x ) w 2 ( t ) d Γ .
Setting a 3 = ( a 1 + | a 2 | ) λ 1 and choosing η = l 4 a 3 leads to (42). □
Lemma 3. 
Under the assumption (H1), the functional Φ 2 defined by
Φ 2 ( t ) = Ω u t ( t ) 0 t ψ ( t s ) ( u ( t ) u ( s ) ) d s d x
satisfies
Φ 2 ( t ) 0 t ψ ( s ) d s δ u t ( t ) 2 + δ u ( t ) 2 + C 1 ( 1 + C ( ξ ) ) δ ( i u ) ( t ) + δ λ 1 w t ( t ) Γ 1 2 + δ a 1 λ 1 Γ 1 f 1 2 ( u t ( t ) ) d Γ + δ | a 2 | λ 1 Γ 1 f 2 2 ( v ( x , 1 , t ) ) d Γ ,
for any 0 < δ < 1 .
Proof. 
Using Equations (22), (24), and (25), we obtain
Φ 2 ( t ) = 1 0 t ψ ( s ) d s Ω u · 0 t ψ ( t s ) ( u ( t ) u ( s ) ) d s d x + Ω 0 t ψ ( t s ) ( u ( t ) u ( s ) ) d s 2 d x Γ 1 w t ( t ) 0 t ψ ( t s ) ( u ( t ) u ( s ) ) d s d Γ + a 1 Γ 1 f 1 ( u t ( t ) ) 0 t ψ ( t s ) ( u ( t ) u ( s ) ) d s d Γ + a 2 Γ 1 f 2 ( v ( x , 1 , t ) ) 0 t ψ ( t s ) ( u ( t ) u ( s ) ) d s d Γ Ω u t ( t ) 0 t ψ ( t s ) ( u ( t ) u ( s ) ) d s d x 0 t ψ ( s ) d s u t ( t ) 2 = ϑ 1 + ϑ 2 + + ϑ 6 0 t ψ ( s ) d s u t ( t ) 2 .
By Young’s inequality, (8), and (44), we obtain, for δ > 0 ,
ϑ 1 δ u ( t ) 2 + C ( ξ ) 4 δ ( i u ) ( t ) , ϑ 2 C ( ξ ) ( i u ) ( t ) , | ϑ 3 | δ λ 1 w t ( t ) Γ 1 2 + C ( ξ ) 4 δ ( i u ) ( t ) , | ϑ 4 | δ a 1 λ 1 Γ 1 f 1 2 ( u t ( t ) ) d Γ + a 1 C ( ξ ) 4 δ ( i u ) ( t ) , | ϑ 5 | δ | a 2 | λ 1 Γ 1 f 2 2 ( v ( x , 1 , t ) ) d Γ + | a 2 | C ( ξ ) 4 δ ( i u ) ( t ) .
Using Young’s inequality, (8), (9), (43), and (44), we see that
ϑ 6 = Ω u t ( t ) 0 t i ( t s ) ( u ( t ) u ( s ) ) d s d x ξ Ω u t ( t ) 0 t ψ ( t s ) ( u ( t ) u ( s ) ) d s d x δ u t ( t ) 2 + 1 2 δ Ω 0 t i ( t s ) | u ( s ) u ( t ) | d s 2 d x + ξ 2 2 δ Ω 0 t ψ ( t s ) | u ( t ) u ( s ) | d s 2 d x δ u t ( t ) 2 + λ 0 ( ψ ( 0 ) + ξ ) 2 δ ( i u ) ( t ) + λ 0 ξ 2 C ( ξ ) 2 δ ( i u ) ( t ) .
Combining all above estimates and taking C 1 = max { λ 0 ( ψ ( 0 ) + ξ ) 2 , δ + 1 + λ 0 ξ 2 2 + a 1 + | a 2 | 4 } , the desired inequality (48) is established. □
Lemma 4. 
Under the assumptions (H3) and (H4), the functional Φ 3 defined by
Φ 3 ( t ) = ϱ ( t ) Γ 1 0 1 e κ ϱ ( t ) F 2 ( v ( x , κ , t ) ) d κ d Γ
satisfies
Φ 3 ( t ) e ϱ 2 ϱ ( t ) Γ 1 0 1 F 2 ( v ( x , κ , t ) ) d κ d Γ c 3 ( 1 ϱ 3 ) e ϱ 2 Γ 1 f 2 ( v ( x , 1 , t ) ) v ( x , 1 , t ) d Γ + c 4 Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ .
Proof. 
Using Equation (23), integration by parts, (13), and (14), we obtain (see [26])
Φ 3 ( t ) = ϱ ( t ) Γ 1 0 1 e κ ϱ ( t ) F 2 ( v ( x , κ , t ) ) d κ d Γ ϱ ( t ) Γ 1 0 1 κ ϱ ( t ) e κ ϱ ( t ) F 2 ( v ( x , κ , t ) ) d κ d Γ Γ 1 0 1 e κ ϱ ( t ) ( 1 κ ϱ ( t ) ) d d κ F 2 ( v ( x , κ , t ) ) d κ d Γ = Φ 3 ( t ) e ϱ ( t ) Γ 1 ( 1 ϱ ( t ) ) F 2 ( v ( x , 1 , t ) ) d Γ + Γ 1 F 2 ( u t ( x , t ) ) d Γ e ϱ 2 ϱ ( t ) Γ 1 0 1 F 2 ( v ( x , κ , t ) ) d κ d Γ c 3 ( 1 ϱ 3 ) e ϱ 2 Γ 1 f 2 ( v ( x , 1 , t ) ) v ( x , 1 , t ) d Γ + c 4 Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ .
Lemma 5 
([3]). Under the assumption (H1), the functional Φ 4 defined by
Φ 4 ( t ) = Ω 0 t G 2 ( t s ) | u ( s ) | 2 d s d x
satisfies
Φ 4 ( t ) 3 ( 1 l ) u ( t ) 2 1 2 ( ψ u ) ( t ) ,
where G 2 ( t ) = t ψ ( s ) d s .
Next, let us define the perturbed modified energy by
L ( t ) = N E ( t ) + N 1 Φ 1 ( t ) + N 2 Φ 2 ( t ) + Φ 3 ( t ) + b 1 E ( t ) ,
where N , N 1 , N 2 , and b 1 are some positive constants.
As in [6,26], for a large enough N > 0 , there exist positive constants β 1 and β 2 such that
β 1 E ( t ) L ( t ) β 2 E ( t ) .
Lemma 6. 
Assume that (H1) and (H3)–(H5) hold. Then, there exist positive constants β 3 , β 4 , and β 5 such that
L ( t ) β 3 E ( t ) + β 4 t 0 t ψ ( s ) Ω | u ( t ) u ( t s ) | 2 d x d s + β 5 Γ 1 f 1 2 ( u t ( t ) ) d Γ , t t 0 ,
where t 0 was introduced in (17).
Proof. 
Let ψ 0 = 0 t 0 ψ ( s ) d s . Using the fact that i ( t ) = ξ ψ ( t ) ψ ( t ) and combining (31), (42), (48), (49), and (51), we obtain, for all t t 0 ,
L ( t ) ξ N 2 ( ψ u ) ( t ) l N 1 2 δ N 2 u ( t ) 2 ψ 0 N 2 δ N 2 N 1 u t ( t ) 2 N 2 2 C ( ξ ) N 1 l C 1 ( 1 + C ( ξ ) ) N 2 δ ( i u ) ( t ) N 1 Γ 1 m ( x ) w 2 ( t ) d Γ + b 1 E ( t ) h 1 N 8 λ 1 N 1 l δ λ 1 N 2 w t ( t ) Γ 1 2 e ϱ 2 ϱ ( t ) Γ 1 0 1 F 2 ( v ( x , κ , t ) ) d κ d Γ γ 0 N c 4 Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ γ 1 N + c 3 ( 1 ϱ 3 ) e ϱ 2 Γ 1 f 2 ( v ( x , 1 , t ) ) v ( x , 1 , t ) d Γ + a 1 a 3 N 1 l + δ a 1 λ 1 N 2 Γ 1 f 1 2 ( u t ( t ) ) d Γ + | a 2 | a 3 N 1 l + δ | a 2 | λ 1 N 2 Γ 1 f 2 2 ( v ( x , 1 , t ) ) d Γ .
From (13), we find that
Γ 1 f 2 2 ( v ( x , 1 , t ) ) d Γ c 2 Γ 1 f 2 ( v ( x , 1 , t ) ) v ( x , 1 , t ) d Γ .
Applying (54) to (53) and taking δ = l 4 N 2 , we obtain, for all t t 0 ,
L ( t ) ξ N 2 ( ψ u ) ( t ) l N 1 2 l 4 u ( t ) 2 ψ 0 N 2 N 1 l 4 u t ( t ) 2 N 2 4 C 1 N 2 2 l C ( ξ ) 2 N 1 l + 4 C 1 N 2 2 l ( i u ) ( t ) N 1 Γ 1 m ( x ) w 2 ( t ) d Γ h 1 N 8 λ 1 N 1 l l λ 1 4 w t ( t ) Γ 1 2 e ϱ 2 ϱ ( t ) Γ 1 0 1 F 2 ( v ( x , κ , t ) ) d κ d Γ γ 0 N c 4 Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ + a 1 a 3 N 1 l + a 1 l λ 1 4 Γ 1 f 1 2 ( u t ( t ) ) d Γ + b 1 E ( t ) γ 1 N + c 3 ( 1 ϱ 3 ) e ϱ 2 | a 2 | a 3 c 2 N 1 l | a 2 | c 2 l λ 1 4 Γ 1 f 2 ( v ( x , 1 , t ) ) v ( x , 1 , t ) d Γ .
We choose N 1 large enough so that
l N 1 2 l 4 > 4 ( 1 l ) ,
then N 2 large enough so that
ψ 0 N 2 N 1 l 4 > 1 .
Using the fact that ξ ψ 2 ( s ) i ( s ) < ψ ( s ) and the Lebesgue dominated convergence theorem, we deduce that
ξ C ( ξ ) = 0 ξ ψ 2 ( s ) i ( s ) d s 0 as ξ 0 .
Hence, there is 0 < ξ 0 < 1 such that if ξ < ξ 0 , then
ξ C ( ξ ) 2 N 1 l + 4 C 1 N 2 2 l < 1 8 .
Finally, selecting ξ = 1 2 N and choosing N large enough so that
N > max 16 C 1 N 2 2 l , 1 h 1 8 λ 1 N 1 l + l λ 1 4 , c 4 γ 0 , 1 γ 1 | a 2 | a 3 c 2 N 1 l + | a 2 | c 2 l λ 1 4 c 3 ( 1 ϱ 3 ) e ϱ 2 ,
we obtain
L ( t ) u t ( t ) 2 4 ( 1 l ) u ( t ) 2 + 1 4 ( ψ u ) ( t ) N 1 Γ 1 m ( x ) w 2 ( t ) d Γ e ϱ 2 ϱ ( t ) Γ 1 0 1 F 2 ( v ( x , κ , t ) ) d κ d Γ + β 5 Γ 1 f 1 2 ( u t ( t ) ) d Γ + b 1 E ( t ) , t t 0 ,
where β 5 = a 1 a 3 N 1 l + a 1 l λ 1 4 . Using (18) and (31), we find that, for any t t 0 ,
0 t 0 ψ ( s ) Ω | u ( t ) u ( t s ) | 2 d x d s 1 c 5 0 t 0 ψ ( s ) Ω | u ( t ) u ( t s ) | 2 d x d s 2 c 5 E ( t ) .
Combining (29), (55), and (56) and making a suitable choice for b 1 , we obtain the estimate (52). □
To evaluate the two terms on the right side of (52), we establish the following lemmas.
Lemma 7 
([1]). Assume that (H2) holds and max { r 1 , f 0 ( r 1 ) } < ε , where ε was introduced in (11). Then, there exist positive constants C 2 , C 3 , and C 4 such that
Γ 1 f 1 2 ( u t ( t ) ) d Γ C 2 Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ , if f 0 is linear ,
Γ 1 f 1 2 ( u t ( t ) ) d Γ C 3 F 1 ( χ ( t ) ) C 3 E ( t ) , if f 0 is nonlinear ,
where
χ ( t ) = 1 | Γ 11 | Γ 11 f 1 ( u t ( t ) ) u t ( t ) d Γ C 4 E ( t ) ,
Γ 11 = { x Γ 1 : | u t ( t ) | ε 1 } and 0 < ε 1 = min { r 1 , f 0 ( r 1 ) } .
Lemma 8. 
Assume that (H1) and (H3)–(H5) hold and that f 0 is linear. Then, the energy functional satisfies
0 E ( s ) d s < .
Proof. 
We introduce the functional
L ( t ) = L ( t ) + Φ 4 ( t ) ,
which is nonnegative. From (50) and (55), we see that, for all t t 0 ,
L ( t ) u t ( t ) 2 ( 1 l ) u ( t ) 2 1 4 ( ψ u ) ( t ) N 1 Γ 1 m ( x ) w 2 ( t ) d Γ e ϱ 2 ϱ ( t ) Γ 1 0 1 F 2 ( v ( x , κ , t ) ) d κ d Γ + β 5 Γ 1 f 1 2 ( u t ( t ) ) d Γ + b 1 E ( t ) .
Applying (29), (31), and (57), we have
L ( t ) d 1 E ( t ) + b 1 β 5 C 2 γ 0 E ( t ) ,
where d 1 is some positive constant. Selecting a suitable choice for b 1 , we obtain
L ( t ) d 1 E ( t ) .
This implies that
d 1 t 0 t E ( s ) d s L ( t 0 ) L ( t ) L ( t 0 ) < .
Next, we define Υ ( t ) by
Υ ( t ) : = t 0 t ψ ( s ) Ω | u ( t ) u ( t s ) | 2 d x d s 2 E ( t ) .
Lemma 9. 
Assume that (H1) and (H2) hold and that G is nonlinear. Then, the solution to (22)–(28) satisfies the estimates
t 0 t ψ ( s ) Ω | u ( t ) u ( t s ) | 2 d x d s 1 θ G ¯ 1 θ Υ ( t ) μ ( t ) , t t 0 , if f 0 is linear ,
t 0 t ψ ( s ) Ω | u ( t ) u ( t s ) | 2 d x d s t t 0 θ G ¯ 1 θ Υ ( t ) ( t t 0 ) μ ( t ) , t > t 0 , if f 0 is nonlinear ,
where θ ( 0 , 1 ) , and G ¯ is an extension of G such that G ¯ is a strictly convex and strictly increasing C 2 function on ( 0 , ) .
Proof. 
First, we prove the estimate (62) when f 0 is linear. For 0 < θ < 1 , we define I ( t ) by
I ( t ) : = θ t 0 t Ω | u ( t ) u ( t s ) | 2 d x d s .
By (60), θ is taken so small that, for all t t 0 ,
I ( t ) < 1 .
Since G is strictly convex on ( 0 , r 0 ] , then
G ( q ζ ) q G ( ζ ) ,
where 0 q 1 and ζ ( 0 , r 0 ] . Using the fact that μ is a positive nonincreasing function and applying (10), (64), (65), and Jensen’s inequality (19), we find that (see details in [1,3])
Υ ( t ) μ ( t ) θ I ( t ) t 0 t I ( t ) G ( ψ ( s ) ) Ω θ | u ( t ) u ( t s ) | 2 d x d s μ ( t ) θ I ( t ) t 0 t G ( I ( t ) ψ ( s ) ) Ω θ | u ( t ) u ( t s ) | 2 d x d s μ ( t ) θ G ¯ θ t 0 t ψ ( s ) Ω | u ( t ) u ( t s ) | 2 d x d s .
Since G ¯ is strictly increasing, we obtain
t 0 t ψ ( s ) Ω | u ( t ) u ( t s ) | 2 d x d s 1 θ G ¯ 1 θ Υ ( t ) μ ( t ) .
Now, we show the estimate (63) when f 0 is nonlinear. Since we cannot guarantee (60), we define the following function:
I 1 ( t ) : = θ t t 0 t 0 t Ω | u ( t ) u ( t s ) | 2 d x d s , t > t 0 .
Using the fact that E ( t ) 0 and (29), we have
I 1 ( t ) 2 θ t t 0 t 0 t ( | | u ( t ) | | 2 + | | u ( t s ) | | 2 ) d s 8 θ E ( 0 ) l .
Choose θ small enough so that, for all t > t 0 ,
I 1 ( t ) < 1 .
Similar to (67), using (10), (65), (67), and Jensen’s inequality (19), we obtain
Υ ( t ) = t t 0 θ I 1 ( t ) t 0 t I 1 ( t ) ( ψ ( s ) ) Ω θ t t 0 | u ( t ) u ( t s ) | 2 d x d s ( t t 0 ) μ ( t ) θ I 1 ( t ) t 0 t G ( I 1 ( t ) ψ ( s ) ) Ω θ t t 0 | u ( t ) u ( t s ) | 2 d x d s ( t t 0 ) μ ( t ) θ G ¯ θ t t 0 t 0 t ψ ( s ) Ω | u ( t ) u ( t s ) | 2 d x d s .
This implies that
t 0 t ψ ( s ) Ω | u ( t ) u ( t s ) | 2 d x d s t t 0 θ G ¯ 1 θ Υ ( t ) ( t t 0 ) μ ( t ) .

4. General Decay of the Energy

In this section, we state and prove the main result of our work.
Theorem 2. 
Assume that (H1)–(H5) hold and that f 0 is linear. Then, there exist positive constants k 1 , k 2 , k 3 , and k 4 such that the energy functional satisfies, for all t t 0 ,
E ( t ) k 2 e k 1 t 0 t μ ( s ) d s , if G is linear ,
E ( t ) k 4 G 1 1 k 3 t 0 t μ ( s ) d s , if G is nonlinear ,
where G 1 ( t ) = t r 0 1 s G ( s ) d s is strictly decreasing and convex on ( 0 , r 0 ] .
Proof. 
Now, we consider the following two cases.
Case 1: G ( t ) is linear. Multiplying (52) by the positive nonincreasing function μ ( t ) and using (10), (31), and (57), we obtain
μ ( t ) L ( t ) β 3 μ ( t ) E ( t ) + β 4 t 0 t μ ( s ) ψ ( s ) Ω | u ( t ) u ( t s ) | 2 d x d s + β 5 μ ( t ) Γ 1 f 1 2 ( u t ( t ) ) d Γ β 3 μ ( t ) E ( t ) β 4 t 0 t ψ ( s ) Ω | u ( t ) u ( t s ) | 2 d x d s + β 5 C 2 μ ( 0 ) Γ 1 f 1 ( u t ( t ) ) u t ( t ) d Γ β 3 μ ( t ) E ( t ) C 5 E ( t ) ,
where C 5 = 2 β 4 + β 5 C 2 μ ( 0 ) γ 0 is a positive constant. Since μ ( t ) is nonincreasing, we have
( μ L + C 5 E ) ( t ) β 3 μ ( t ) E ( t ) , t t 0 .
Since μ ( t ) L ( t ) + C 5 E ( t ) E ( t ) , for some positive constants k 1 and k 2 , we obtain
E ( t ) k 2 e k 1 t 0 t μ ( s ) d s .
Case 2: G ( t ) is nonlinear. This case is obtained through the ideas presented in [3] as follows.
Using (31), (52), (57), and (62), we obtain
L ( t ) β 3 E ( t ) + β 4 θ G ¯ 1 θ Υ ( t ) μ ( t ) β 5 C 2 γ 0 E ( t ) , t t 0 .
Let L 1 ( t ) = L ( t ) + β 5 C 2 γ 0 E ( t ) E ( t ) , and then (70) becomes
L 1 ( t ) β 3 E ( t ) + β 4 θ G ¯ 1 θ Υ ( t ) μ ( t ) , t t 0 .
For 0 < ε 0 < r 0 , using (71) and the fact that E 0 , G ¯ > 0 and G ¯ > 0 , we find that the functional L 2 , defined by
L 2 ( t ) : = G ¯ ε 0 E ( t ) E ( 0 ) L 1 ( t ) E ( t ) ,
satisfies
L 2 ( t ) β 3 E ( t ) G ¯ ε 0 E ( t ) E ( 0 ) + β 4 θ G ¯ ε 0 E ( t ) E ( 0 ) G ¯ 1 θ Υ ( t ) μ ( t ) , t t 0 .
With s = G ¯ ε 0 E ( t ) E ( 0 ) and r = G ¯ 1 θ Υ ( t ) μ ( t ) , using (20), (21), and (72), we obtain
L 2 ( t ) β 3 E ( t ) G ε 0 E ( t ) E ( 0 ) + ε 0 β 4 θ E ( t ) E ( 0 ) G ε 0 E ( t ) E ( 0 ) + β 4 Υ ( t ) μ ( t ) ,
where we have used that ε 0 E ( t ) E ( 0 ) < r 0 and G ¯ = G on ( 0 , r 0 ] . Multiplying this by μ ( t ) and using (61), we obtain
μ ( t ) L 2 ( t ) β 3 E ( 0 ) ε 0 β 4 θ μ ( t ) E ( t ) E ( 0 ) G ε 0 E ( t ) E ( 0 ) 2 β 4 E ( t ) .
By defining L 3 ( t ) = μ ( t ) L 2 ( t ) + 2 β 4 E ( t ) , we see that, for some positive constants γ 2 and γ 3 ,
γ 2 L 3 ( t ) E ( t ) γ 3 L 3 ( t ) .
With a suitable choice of ε 0 , we obtain, for some positive constant d 2 ,
L 3 ( t ) d 2 μ ( t ) E ( t ) E ( 0 ) G ε 0 E ( t ) E ( 0 ) = d 2 μ ( t ) G 2 E ( t ) E ( 0 ) , t t 0 ,
where G 2 ( t ) = t G ( ε 0 t ) . Applying the strict convexity of G on ( 0 , r 0 ] and G 2 ( t ) = G ( ε 0 t ) + ε 0 t G ( ε 0 t ) , we see that G 2 ( t ) , G 2 ( t ) > 0 on ( 0 , 1 ] . Finally, defining
Q ( t ) = γ 2 L 3 ( t ) E ( 0 )
and using (73), we have
Q ( t ) E ( t ) E ( 0 ) 1 and Q ( t ) E ( t ) .
From (74), (75), and the fact that G 2 ( t ) > 0 on ( 0 , 1 ] , we arrive at
Q ( t ) k 3 μ ( t ) G 2 ( Q ( t ) ) , t t 0 ,
where k 3 = d 2 γ 2 E ( 0 ) is a positive constant. Integrating this over ( t 0 , t ) and using variable transformation, we find that (see details in [3])
t t 0 ε 0 Q ( s ) ε 0 Q ( s ) G ( ε 0 Q ( s ) ) d s k 3 t 0 t μ ( s ) d s ε 0 Q ( t ) ε 0 Q ( t 0 ) 1 s G ( s ) d s k 3 t 0 t μ ( s ) d s .
Since ε 0 < r 0 and Q ( t ) 1 , for all t t 0 , we have
G 1 ( ε 0 Q ( t ) ) = ε 0 Q ( t ) r 0 1 s G ( s ) d s k 3 t 0 t μ ( s ) d s Q ( t ) 1 ε 0 G 1 1 k 3 t 0 t μ ( s ) d s ,
where G 1 ( t ) = t r 0 1 s G ( s ) d s . Here, we have used the fact that G 1 is a strictly decreasing function on ( 0 , r 0 ] . Therefore, using (75) and (76), the estimate (69) is established. □
Theorem 3. 
Assume that (H1)–(H5) hold and that f 0 is nonlinear. Then, there exist positive constants α 1 , α 2 , α 3 , and α 4 such that the energy functional satisfies
E ( t ) α 2 F 1 1 α 1 t 0 t μ ( s ) d s , t t 0 , if G is linear ,
where F 1 ( t ) = t r 1 1 s F ( s ) d s and
E ( t ) α 4 ( t t 0 ) K 1 1 α 3 ( t t 0 ) t 1 t μ ( s ) d s , t t 1 , if G is nonlinear ,
where K 1 ( t ) = t K ( ε 2 t ) , 0 < ε 2 < r 2 = min { r 0 , r 1 } and K = G ¯ 1 + F ¯ 1 1 .
Proof. Case 1: G ( t ) is linear. Multiplying (52) by the positive nonincreasing function μ ( t ) and using (10), (31), and (58), we obtain
μ ( t ) L ( t ) β 3 μ ( t ) E ( t ) + β 5 C 3 μ ( t ) F 1 ( χ ( t ) ) C 6 E ( t ) ,
where C 6 = 2 β 4 + β 5 C 3 μ ( 0 ) is a positive constant. Since μ ( t ) is nonincreasing, (79) becomes
F 3 ( t ) β 3 μ ( t ) E ( t ) + β 5 C 3 μ ( t ) F 1 ( χ ( t ) ) , t t 0 ,
where F 3 ( t ) = μ ( t ) L ( t ) + C 6 E ( t ) E ( t ) . For 0 < ε 1 < r 1 , using (80) and the fact that E 0 , F > 0 and F > 0 on ( 0 , r 1 ] , the functional F 4 , defined by
F 4 ( t ) : = F ε 1 E ( t ) E ( 0 ) F 3 ( t ) E ( t ) ,
satisfies
F 4 ( t ) β 3 μ ( t ) E ( t ) F ε 1 E ( t ) E ( 0 ) + β 5 C 3 μ ( t ) F ε 1 E ( t ) E ( 0 ) F 1 ( χ ( t ) ) .
Given (20) and (21) with s = F ε 1 E ( t ) E ( 0 ) and r = F 1 ( χ ( t ) ) , using (59), we obtain that
F 4 ( t ) β 3 μ ( t ) E ( t ) F ε 1 E ( t ) E ( 0 ) + ε 1 β 5 C 3 μ ( t ) E ( t ) E ( 0 ) F ε 1 E ( t ) E ( 0 ) + β 5 C 3 μ ( 0 ) χ ( t ) β 3 E ( 0 ) ε 1 β 5 C 3 μ ( t ) E ( t ) E ( 0 ) F ε 1 E ( t ) E ( 0 ) β 5 C 3 C 4 μ ( 0 ) E ( t ) , t t 0 .
Let F 5 ( t ) = F 4 ( t ) + β 5 C 3 C 4 μ ( 0 ) E ( t ) ; then it satisfies, for positive constants γ 4 and γ 5 ,
γ 4 F 5 ( t ) E ( t ) γ 5 F 5 ( t ) .
Consequently, with a suitable choice of ε 1 , we have, for some positive constant d 3 ,
F 5 ( t ) d 3 μ ( t ) E ( t ) E ( 0 ) F ε 1 E ( t ) E ( 0 ) = d 3 μ ( t ) F 0 E ( t ) E ( 0 ) , t t 0 ,
where F 0 ( t ) = t F ( ε 1 t ) . From the strict convexity of F on ( 0 , r 1 ] , we obtain F 0 ( t ) , F 0 ( t ) > 0 on ( 0 , 1 ] . Let
J ( t ) = γ 4 F 5 ( t ) E ( 0 ) ,
and from (81) and (82), we obtain
J ( t ) E ( t ) E ( 0 ) 1 and J ( t ) α 1 μ ( t ) F 0 ( J ( t ) ) , t t 0 ,
where α 1 = d 3 γ 4 E ( 0 ) is a positive constant. Then, similar to (76), the integration over ( t 0 , t ) and variable transformation yield
J ( t ) 1 ε 1 F 1 1 α 1 t 0 t μ ( s ) d s ,
where F 1 ( t ) = t r 1 1 s F ( s ) d s , which is a strictly decreasing function on ( 0 , r 1 ] . Combining (81) and (83), the estimate (77) is proved.
Case 2: G ( t ) is nonlinear. This case is obtained by the arguments presented in [1] as follows. Using (52), (58), and (63), we obtain
L ( t ) β 3 E ( t ) + β 4 ( t t 0 ) θ G ¯ 1 θ Υ ( t ) ( t t 0 ) μ ( t ) + β 5 C 3 F 1 ( χ ( t ) ) β 5 C 3 E ( t ) , t > t 0 .
Since lim t 1 t t 0 = 0 , there exists t 1 > t 0 such that
1 t t 0 < 1 , t t 1 .
Using the strictly convex and strictly increasing function of F ¯ and (65) with q = 1 t t 0 , we see that
F ¯ 1 ( χ ( t ) ) ( t t 0 ) F ¯ 1 χ ( t ) t t 0 , t t 1 .
Combining (84) and (86), we arrive at
R 1 ( t ) β 3 E ( t ) + β 4 ( t t 0 ) θ G ¯ 1 θ Υ ( t ) ( t t 0 ) μ ( t ) + β 5 C 3 ( t t 0 ) F ¯ 1 χ ( t ) t t 0 , t t 1 ,
where R 1 ( t ) = L ( t ) + β 5 C 3 E ( t ) E ( t ) . Let
r 2 = min { r 0 , r 1 } , φ ( t ) = max θ Υ ( t ) ( t t 0 ) μ ( t ) , χ ( t ) t t 0 and K = G ¯ 1 + F ¯ 1 1 , t t 1 .
Therefore, (87) reduces to
R 1 ( t ) β 3 E ( t ) + C 7 ( t t 0 ) K 1 ( φ ( t ) ) , t t 1 ,
where C 7 = max { β 4 θ , β 5 C 3 } . The strictly increasing and strictly convex properties of G ¯ and F ¯ imply that
K = G ¯ F ¯ G ¯ + F ¯ > 0 and K = G ¯ ( F ¯ ) 2 + ( G ¯ ) 2 F ¯ ( G ¯ + F ¯ ) 2 > 0 ,
on ( 0 , r 2 ] .
Now, for 0 < ε 2 < r 2 , using (85), we see that ε 2 t t 0 E ( t ) E ( 0 ) < r 2 . Defining
R 2 ( t ) = K ε 2 t t 0 E ( t ) E ( 0 ) R 1 ( t ) , t t 1 ,
and using (89) and (90), we find that
R 2 ( t ) = ε 2 ( t t 0 ) 2 E ( t ) E ( 0 ) + ε 2 t t 0 E ( t ) E ( 0 ) K ε 2 t t 0 E ( t ) E ( 0 ) R 1 ( t ) + K ε 2 t t 0 E ( t ) E ( 0 ) R 1 ( t ) β 3 E ( t ) K ε 2 t t 0 E ( t ) E ( 0 ) + C 7 ( t t 0 ) K ε 2 t t 0 E ( t ) E ( 0 ) K 1 ( φ ( t ) ) , t t 1 .
Using (20) and (21) with s = K ε 2 t t 0 E ( t ) E ( 0 ) and r = K 1 ( φ ( t ) ) and applying (91), we obtain
R 2 ( t ) β 3 E ( t ) K ε 2 t t 0 E ( t ) E ( 0 ) + ε 2 C 7 E ( t ) E ( 0 ) K ε 2 t t 0 E ( t ) E ( 0 ) + C 7 ( t t 0 ) φ ( t ) .
From (59), (61), and (88), we obtain
( t t 0 ) μ ( t ) φ ( t ) C 8 E ( t ) ,
where C 8 = min { 2 θ , C 4 μ ( 0 ) } . Multiplying (92) by the positive nonincreasing function μ ( t ) and using (93), we have
R 3 ( t ) β 3 E ( 0 ) ε 2 C 7 μ ( t ) E ( t ) E ( 0 ) K ε 2 t t 0 E ( t ) E ( 0 ) , t t 1 ,
where R 3 ( t ) = μ ( t ) R 2 ( t ) + C 7 C 8 E ( t ) E ( t ) . For a suitable choice of ε 2 , we find that
R 3 ( t ) d 4 μ ( t ) E ( t ) E ( 0 ) K ε 2 t t 0 E ( t ) E ( 0 ) , t t 1 ,
where d 4 is a positive constant. An integration of (94) yields
d 4 E ( 0 ) t 1 t E ( s ) K ε 2 s t 0 E ( s ) E ( 0 ) μ ( s ) d s t t 1 R 3 ( s ) d s R 3 ( t 1 ) .
Using (90) and the non-increasing property of E , we see that the map t E ( t ) K ε 2 t t 0 E ( t ) E ( 0 ) is non-increasing and, consequently, we obtain
d 4 E ( t ) E ( 0 ) K ε 2 t t 0 E ( t ) E ( 0 ) t 1 t μ ( s ) d s R 3 ( t 1 ) , t t 1 .
Multiplying (95) by 1 t t 0 , we obtain
d 4 K 1 1 t t 0 E ( t ) E ( 0 ) t 1 t μ ( s ) d s R 3 ( t 1 ) t t 0 , t t 1 ,
where K 1 ( s ) = s K ( ε 2 s ) , which is strictly increasing. Therefore, we deduce that
E ( t ) α 4 ( t t 0 ) K 1 1 α 3 ( t t 0 ) t 1 t μ ( s ) d s , t t 1 ,
where α 3 and α 4 are positive constants. This completes the proof. □
Examples. We provide examples to explain the decay of energy (see [1]).
  • Case: f 0 and G are linear.
Let ψ ( t ) = a e b ( 1 + t ) , μ ( t ) = b , and G ( t ) = t , where b > 0 , and a > 0 is small enough. Assume that f 0 ( t ) = c t and F ( t ) = t f 0 ( t ) = c t . Then, we can obtain
E ( t ) k 2 e k 1 t , for all t t 0 .
2.
Case: f 0 is linear and G is nonlinear.
Let ψ ( t ) = a e t p , μ ( t ) = 1 , and G ( t ) = p t ( ln ( a t ) ) 1 p 1 , where 0 < p < 1 , and a > 0 is small enough. Assume that f 0 ( t ) = c t and F ( t ) = t f 0 ( t ) = c t . Then, G satisfies the condition (H1) on ( 0 , r 0 ] for any 0 < r 0 < a .
G 1 ( t ) = t r 0 1 s G ( s ) d s = t r 0 [ ln a s ] 1 p s [ 1 p + p ln a s ] d s = ln a r 0 ln a t u 1 p 1 p + p u d u ln a t 1 p .
Then, we can have
E ( t ) k 4 e k 3 t p , for all t t 0 .
3.
Case: f 0 is nonlinear and G is linear.
Let ψ ( t ) = a e b ( 1 + t ) , μ ( t ) = b , and G ( t ) = t , where b > 0 , and a > 0 is small enough. Assume that f 0 ( t ) = c t p , where p > 1 and F ( t ) = t f 0 ( t ) = c t p + 1 2 . Then,
F 1 ( t ) = t r 1 1 s F ( s ) d s = t r 1 2 c ( p + 1 ) s p + 1 2 d s = α 0 r 1 p 1 2 t p 1 2
and
F 1 1 ( t ) = ( r 1 p 1 2 + 1 α 0 t ) 2 p 1 ,
where α 0 = 4 c ( p + 1 ) ( p 1 ) . Therefore, we find that
E ( t ) ( α 1 t + α 2 ) 2 p 1 , for all t t 0 .
4.
Case: f 0 is nonlinear and G is nonlinear.
Let ψ ( t ) = a ( 1 + t ) 2 , μ ( t ) = b , and G ( t ) = t 3 2 , where b > 0 , and a > 0 is taken so that (9) remains valid. Assume that f 0 ( t ) = t 5 and F ( t ) = t 3 . Then,
K ( s ) = ( G 1 + F 1 ) 1 ( s ) = 1 + 1 + 4 s 2 3 .
Therefore, we see that
E ( t ) α 3 ( t t 0 ) 1 3 , for all t t 1 ,
where t 1 > t 0 .

5. Conclusions

Numerous phenomena are influenced by both the current state and the previous occurrences of the system. There has been a notable increase in the research on the equation with delay effects, which frequently arise in various physical, biological, chemical, medical, and economic problems. In this paper, we study the energy decay rates for the viscoelastic wave equation with nonlinear time-varying delay, nonlinear damping at the boundary, and acoustic boundary conditions. We consider the relaxation function ψ , namely ψ ( t ) μ ( t ) G ( ψ ( t ) ) , where G is an increasing and convex function near the origin, and μ is a positive nonincreasing function. We establish general decay rate results without the need for the condition a 2 > 0 and without imposing any limiting growth assumption on the damping term f 1 , using the multiplier method and some properties of the convex functions. Moreover, the energy decay rates depend on the functions μ and G, as well as the function F defined by f 0 , which characterizes the growth behavior of f 1 at the origin.

Author Contributions

M.J.L. and J.-R.K. contributed equally to the writing of this paper and reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Pukyong National University Industry–University Cooperation Research Fund in 2023 (202311540001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors would like to thank the handling editor and the referees for their relevant remarks and corrections in order to improve the final version.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Al-Gharabli, M.M.; Al-Mahdi, A.M.; Messaoudi, S.A. General and optimal decay result for a viscoelastic problem with nonlinear boundary feedback. J. Dyn. Control Syst. 2019, 25, 551–572. [Google Scholar] [CrossRef]
  2. Messaoudi, S.A. General decay of solutions of a viscoelastic equation. J. Math. Anal. Appl. 2008, 341, 1457–1467. [Google Scholar] [CrossRef]
  3. Mustafa, M.I. Optimal decay rates for the viscoelastic wave equation. Math. Meth. Appl. Sci. 2018, 41, 192–204. [Google Scholar] [CrossRef]
  4. Beale, J.T.; Rosencrans, S.I. Acoustic boundary conditions. Bull. Am. Math. Soc. 1974, 80, 1276–1278. [Google Scholar] [CrossRef]
  5. Munoz Rivera, J.E.; Qin, Y.M. Polynomial decay for the energy with an acoustic boundary condition. Appl. Math. Lett. 2003, 16, 249–256. [Google Scholar] [CrossRef]
  6. Park, J.Y.; Park, S.H. Decay rate estimates for wave equation of memory type with acoustic boundary conditions. Nonlinear Anal. Theory Methods Appl. 2011, 74, 993–998. [Google Scholar] [CrossRef]
  7. Liu, W.J. Arbitrary rate of decay for a viscoelastic equation with acoustic boundary coditions. Appl. Math. Lett. 2014, 38, 155–161. [Google Scholar] [CrossRef]
  8. Yoon, M.; Lee, M.J.; Kang, J.R. General decay result for the wave equation with memory and acoustic boundary conditions. Appl. Math. Lett. 2023, 135, 108385. [Google Scholar] [CrossRef]
  9. Feng, B.W. Long-time dynamics of a plate equation with memory and time delay. Bull. Braz. Math. Soc. 2018, 49, 395–418. [Google Scholar] [CrossRef]
  10. Nicaise, S.; Pignotti, D. Stability and instability results of the wave equation with a delay term in the boundary or internal feedbacks. SIAM J. Control Optim. 2006, 45, 1561–1585. [Google Scholar] [CrossRef]
  11. Nicaise, S.; Pignotti, C. Stability of the wave equation with boundary or internal distributed delay. Differ. Integral Equ. 2008, 21, 935–958. [Google Scholar]
  12. Kirane, M.; Said-Houari, B. Existence and asymptotic stability of a viscoelastic wave equation with a delay. Z. Angew. Math. Phys. 2011, 62, 1065–1082. [Google Scholar] [CrossRef]
  13. Dai, Q.; Yang, Z.F. Global existence and exponential decay of the solution for a viscoelastic wave equation with a delay. Z. Angew. Math. Phys. 2014, 65, 885–903. [Google Scholar] [CrossRef]
  14. Nicaise, S.; Pignotti, C. Interior feedback stabilization of wave equations with time dependent delay. Electron. J. Differ. Equ. 2011, 2011, 1–20. [Google Scholar]
  15. Liu, W.J. General decay rate estimate for the energy of a weak viscoelastic equation with an internal time-varying delay term. Taiwanese J. Math. 2013, 17, 2101–2115. [Google Scholar] [CrossRef]
  16. Feng, B.W. Well-posedness and exponential stability for a plate equation with time-varying delay and past history. Z. Angew. Math. Phys. 2017, 68, 1–24. [Google Scholar] [CrossRef]
  17. Lee, M.J.; Kim, D.W.; Park, J.Y. General decay of solutions for Kirchhoff type containing Balakrishnan-Taylor damping with a delay and acoustic boundary conditions. Bound. Value Probl. 2016, 2016, 173. [Google Scholar] [CrossRef]
  18. Liu, G.W.; Diao, L. Energy decay of the solution for a weak viscoelastic equation with a time-varying delay. Acta Appl. Math. 2018, 155, 9–19. [Google Scholar] [CrossRef]
  19. Mustafa, M.I. Asymptotic behavior of second sound thermoelasticity with internal time-varying delay. Z. Angew. Math. Phys. 2013, 64, 1353–1362. [Google Scholar] [CrossRef]
  20. Park, S.H. Decay rate estimates for a weak viscoelastic beam equation with time-varying delay. Appl. Math. Lett. 2014, 31, 46–51. [Google Scholar] [CrossRef]
  21. Park, S.H.; Kang, J.R. General decay for weak viscoelastic Kirchhoff plate equations with delay boundary conditions. Bound. Value Probl. 2017, 2017, 96. [Google Scholar] [CrossRef]
  22. Zitouni, S.; Zennir, K.; Bouzettouta, L. Uniform decay for a viscoelastic wave equation with density and time-varying delay in Rn. Filomat. 2019, 33, 961–970. [Google Scholar] [CrossRef]
  23. 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]
  24. Benaissa, A.; Benaissa, A.; Messaoudi, S.A. Global existence and energy decay of solutions for the wave equation with a time varying delay term in the weakly nonlinear internal feedbacks. J. Math. Phys. 2012, 53, 123514. [Google Scholar] [CrossRef]
  25. Park, S.H. Energy decay for a von Karman equation with time-varying delay. Appl. Math. Lett. 2016, 55, 10–17. [Google Scholar] [CrossRef]
  26. Kang, J.R.; Lee, M.J.; Park, S.H. Asymptotic stability of a viscoelastic problem with Balakrishnan-Taylor damping and time-varying delay. Comput. Math. Appl. 2017, 74, 1506–1515. [Google Scholar] [CrossRef]
  27. Lee, M.J.; Park, J.Y.; Park, S.H. General decay of solutions of quasilinear wave equation with time-varying delay in the boundary feedback and acoustic boundary conditions. Math. Meth. Appl. Sci. 2017, 40, 4560–4576. [Google Scholar] [CrossRef]
  28. Djeradi, F.S.; Yazid, F.; Georgiev, S.G.; Hajjej, Z.; Zennir, K. On the time decay for a thermoelastic laminated beam with microtemperature effects, nonlinear weight, and nonlinear time-varying delay. AIMS Math. 2023, 8, 26096–26114. [Google Scholar] [CrossRef]
  29. Mukiawa, S.E.; Enyi, C.D.; Messaoudi, S.A. Stability of thermoelastic Timoshenko beam with suspenders and time-varying feedback. Adv. Contin. Disc. Models. 2023, 2023, 7. [Google Scholar] [CrossRef]
  30. Al-Gharabli, M.M.; Balegh, M.; Feng, B.W.; Hajjej, Z.; Messaoudi, S.A. Existence and general decay of Balakrishnan-Taylor viscoelastic equation with nonlinear frictional damping and logarithmic source term. Evol. Equ. Control Theory. 2022, 11, 1149–1173. [Google Scholar] [CrossRef]
  31. Alabau-Boussouira, F.; Cannarsa, P. A general method for proving sharp energy decay rates for memory dissipative evolution equations. Comptes Rendus Math. 2009, 347, 867–872. [Google Scholar] [CrossRef]
  32. Arnold, V.I. Mathematical Methods of Classical Mechanics; Springer: New York, NY, USA, 1989. [Google Scholar]
  33. Park, J.Y.; Ha, T.G. Well-posedness and uniform decay rates for the Klein–Gordon equation with damping term and acoustic boundary conditions. J. Math. Phys. 2009, 50, 013506. [Google Scholar] [CrossRef]
  34. Jin, K.P.; Liang, J.; Xiao, T.J. Coupled second order evolution equations with fading memory: Optimal energy decay rate. J. Differ. Equ. 2014, 257, 1501–1528. [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

Lee, M.J.; Kang, J.-R. General Stability for the Viscoelastic Wave Equation with Nonlinear Time-Varying Delay, Nonlinear Damping and Acoustic Boundary Conditions. Mathematics 2023, 11, 4593. https://doi.org/10.3390/math11224593

AMA Style

Lee MJ, Kang J-R. General Stability for the Viscoelastic Wave Equation with Nonlinear Time-Varying Delay, Nonlinear Damping and Acoustic Boundary Conditions. Mathematics. 2023; 11(22):4593. https://doi.org/10.3390/math11224593

Chicago/Turabian Style

Lee, Mi Jin, and Jum-Ran Kang. 2023. "General Stability for the Viscoelastic Wave Equation with Nonlinear Time-Varying Delay, Nonlinear Damping and Acoustic Boundary Conditions" Mathematics 11, no. 22: 4593. https://doi.org/10.3390/math11224593

APA Style

Lee, M. J., & Kang, J. -R. (2023). General Stability for the Viscoelastic Wave Equation with Nonlinear Time-Varying Delay, Nonlinear Damping and Acoustic Boundary Conditions. Mathematics, 11(22), 4593. https://doi.org/10.3390/math11224593

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