Next Article in Journal
Least-Squares Estimators of Drift Parameter for Discretely Observed Fractional Ornstein–Uhlenbeck Processes
Next Article in Special Issue
Reflection-Like Maps in High-Dimensional Euclidean Space
Previous Article in Journal
General Homogenised Formulation for Thick Viscoelastic Layered Structures for Finite Element Applications
Previous Article in Special Issue
Boundary Value Problems for a Class of First-Order Fuzzy Delay Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Logarithmic Decay of Wave Equation with Kelvin-Voigt Damping

1
Laboratoire de Mathématiques Appliquées, UMR 8100 du CNRS, Université Paris–Saclay (site UVSQ), 45 avenue des Etats Unis, 78035 Versailles, France
2
School of Mathematics and Statistics, Beijing Key Laboratory on MCAACI, Beijing Institute of Technology, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(5), 715; https://doi.org/10.3390/math8050715
Submission received: 2 April 2020 / Revised: 22 April 2020 / Accepted: 25 April 2020 / Published: 3 May 2020
(This article belongs to the Special Issue Modern Analysis and Partial Differential Equation)

Abstract

:
In this paper, we analyze the longtime behavior of the wave equation with local Kelvin-Voigt Damping. Through introducing proper class symbol and pseudo-diff-calculus, we obtain a Carleman estimate, and then establish an estimate on the corresponding resolvent operator. As a result, we show the logarithmic decay rate for energy of the system without any geometric assumption on the subdomain on which the damping is effective.
MSC:
Primary 93B05; Secondary 93B07 and 35B37

1. Introduction

In this paper, a wave equation with local Kelvin-Voigt damping is considered. More precisely, we assume that the wave propagates through two segments consisting of an elastic and a Kelvin-Voigt medium. The latter material is a viscoelastic material having the properties both of elasticity and viscosity. We analyze long time behaviour of energy of solution of the system. Let Ω R d be a bounded domain with smooth boundary Γ = Ω . Denote by n the unit outward normal vector on boundary Γ . The PDE model is as follows.
y t t ( t , x ) div [ y ( t , x ) + a ( x ) y t ( t , x ) ] = 0 in ( 0 , ) × Ω , y ( t , x ) = 0 on ( 0 , ) × Γ , y ( 0 , x ) = y 0 , y t ( 0 , x ) = y 1 in Ω ,
where the coefficient function a ( · ) L 1 ( Ω ) is non-negative and not identically null.
The natural energy of system (1) is
E ( t ) = 1 2 Ω | y ( t ) | 2 + | y t ( t ) | 2 d x .
A direct computation gives that
d d t E ( t ) = supp a a ( x ) | y t ( t ) | 2 d x .
Formula (3) shows that the only dissipative mechanism acting on the system is the viscoelastic damping div [ a y t ] , which is only effective on supp a .
To rewrite the system as an evolution equation, we set the energy space as
H = H 0 1 ( Ω ) × L 2 ( Ω ) ,
with norm
Y H = y 1 L 2 ( Ω ) 2 + y 2 L 2 ( Ω ) 2 , Y = y 1 , y 2 H ,
Define an unbounded operator A : D ( A ) H H by
A Y = y 2 , div ( y 1 + a y 2 ) , Y = y 1 , y 2 D ( A ) ,
and
D ( A ) = y 1 , y 2 H : y 2 H 0 1 ( Ω ) , div ( y 1 + a y 2 ) L 2 ( Ω ) } .
Let Y ( t ) = ( y ( t ) , y t ( t ) ) . Then system (1) can be written as
d d t Y ( t ) = A Y ( t ) , t 0 , Y ( 0 ) = ( y 0 , y 1 ) .
It is known from Reference [1,2,3] that if supp a is non-empty, the operator A generates a contractive C 0 semigroup e t A on H and i R ρ ( A ) , the resolvent of A . Consequently, the semigroup e t A is strongly stable. Moreover, if the entire medium is of the viscoelastic type (i.e., supp a = Ω ¯ ), the damping for the wave equation not only induces exponential energy decay, but also restricts the spectrum of the associated semigroup generator to a sector in the left half plane, and the associated semigroup is analytic [4]. However, when the Kelvin-Voigt damping is local and the material coefficient a ( · ) is a positive constant on supp a , the energy of system (1) does not decay exponentially for any geometry of Ω and supp a [1,5]. The reason is that the strong damping leads to the reflection of waves at the interface γ ( supp a ) \ Γ , which then fails to be effectively damped because they do not enter the region of damping [6,7,8]. It turns out that the viscoelastic damping does not follow the assumption that the “geometric optics” condition implies exponential stability [9].
On the other hand, it has been proved that the properties of regularity and the stability of the 1-d system (1) depend on the continuousness of coefficient function a ( · ) . More precisely, assume that Ω = ( 1 , 1 ) and a ( x ) behaviours like x α with α > 0 in supp a = [ 0 , 1 ] . Then the solution of (1) is eventually differentiable for α > 1 , exponentially stable for α 1 , polynomially stable of order 1 1 α for 0 < α < 1 , and polynomially stable of optimal decay rate 2 for α = 0 (see [10,11,12,13]).
For the higher dimensional system, the corresponding semigroup is exponentially stable when a ( · ) C 2 ( Ω ) and supp a Γ , polynomially stable of order 1 2 when a ( · ) a 0 > 0 on supp a and supp a satisfies certain geometry conditions [14,15]. Then, a natural problem is—how about the decay rate when Ω and supp a is arbitrary? In this paper, we analyze the logarithmic decay properties of the solution to (1) and obtain that the decay rate of system energy stays at a rate 1 [ log ( t + 1 ) ] 8 / 5 . This result was recently improved by Burq [16].
The main result reads as follows.
Theorem 1.
Suppose that the coefficient function a ( · ) C 0 ( Ω ) is non-negative and supp a Ω is non-empty. Then the energy of the solution of (1) decays at logarithmic speed. More precisely, one has that there exists a positive constant C such that
e t A Y 0 H C [ log ( t + 1 ) ] 4 k / 5 Y 0 D ( A k ) , t > 0 , Y 0 = ( y 0 , y 1 ) D ( A k ) .
We shall prove Theorem 1 through the resolvent estimate [17]. The main idea is to introduce proper operators, class of symbol and pseudo-diff-calculus. Then, a Carleman estimate on the subdomain far away from the boundary will be proven. Combining these with a classical Carleman estimation up to the boundary, one can obtain the desired Carleman estimate and resolvent estimate. This method was developed in References [18,19,20,21,22,23] and the references cited therein.
The rest of the paper is organized as follows—we present some preliminaries in Section 2. Section 3 and Section 4 are devoted to the proof of the Carleman estimate and the logarithmic stability of system (1). Finally, some proofs of the classic result are given in Appendix A to complete the paper.
Throughout this paper, we use · , ( · | · ) to denote the norm and inner product in L 2 ( Ω ) if there are no comments. When writing f g (or f g ), we mean that there exists a positive constant C such that f C g (or f C g ).

2. Preliminaries

We shall prove Theorem 1 by Weyl-H o ¨ rmander calculus, which was introduced H o ¨ rmander [24,25]. First, let’s define some definitions and results on the class of symbol and pseudo-diff-calculus.

2.1. Symbol and Symbolic Calculus

Let V be a bounded open set in R d . For any ( x , ξ ) V ¯ × R d , λ R and τ > 0 , we introduce the metric
g = g x , ξ = λ d x 2 + μ 2 d ξ 2 , where μ 2 = μ ( τ , ξ ) 2 = τ 2 + | ξ | 2 ,
and the weight
ν = ν ( x , λ ) = 1 + λ 2 a ( x ) 2 .
Note that g x , ξ ( X , Ξ ) = λ | X | 2 + μ 2 ( τ , ξ ) | Ξ | 2 for all X , Ξ R d . Then we have the following results. It’s proof is given in Appendix A for the sake of completeness.
Lemma 1.
Assume that there exist positive constants C and λ 0 such that λ λ 0 and τ max { C λ , 1 } . It holds
(i)
The metric g = g x , ξ defined by (8) is admissible, i.e., it is slowly varying and temperate.
(ii)
The weight ν = ν ( x , λ ) defined by (9) is admissible, that is, it is g-continuous and g-temperate.
Definition 1
(Section 18.4.2 in [24]). Let the weight m ( x , ξ ) be admissible and the metric g be defined by (8). Assume q ( x , ξ , λ , τ ) C ( R d × R d ) , where parameters λ , τ satisfy conditions in Lemma 1. q ( x , ξ , λ , τ ) is a symbol in class S ( m , g ) if for all α , β N d there exist C α , β independent of τ and λ such that
| x α ξ β q ( x , ξ , λ , τ ) | C α , β m ( x , ξ ) λ | α | / 2 μ ( τ , ξ ) | β | .
Remark 1.
(i)
It is clear that μ = τ 2 + | ξ | 2 S ( μ , g ) since | ξ β μ ( τ , ξ ) | μ 1 | β | for all β R d .
(ii)
Let ν be the weight defined by (9). It is easy to get that λ a S ( ν , g ) . In fact, if | α | 2 , it holds that | x α ( λ a ( x ) ) | C α λ C α λ | α | / 2 ν ( x ) , where C α > 0 . For the case | α | = 1 , it follows from (A1) that
| x α ( λ a ( x ) ) | 2 a L ( Ω ) 1 2 λ 1 2 λ a ( x ) 1 2
Note that | λ a ( x ) | < C ν 2 ( x ) for some C > 0 . This together with the above inequality, we have that there exists a positive constant C such that
| x α ( λ a ( x ) ) | < 2 C a L ( Ω ) 1 2 λ 1 2 ν ( x ) .
(iii)
It is known from Lemma 18.4.3 of [24] that if the metric g and weights m 1 , m 2 are admissible, symbols a S ( m 1 , g ) and b S ( m 2 , g ) , then a b S ( m 1 m 2 , g ) . In particular, ( λ a ) j μ k S ( ν j μ k , g ) for all j , k N { 0 } .
Definition 2.
Let b S ( m , g ) be a symbol and u S ( R d ) , we set
b ( x , D , τ ) u ( x ) = Op ( b ) u ( x ) : = ( 2 π ) d R d e i x · ξ b ( x , ξ , τ ) u ^ ( ξ ) d ξ .
It is known that Op ( b ) : S ( R d ) S ( R d ) is continuous and Op ( b ) can be uniquely extended to S ( R d ) continuously. The following two lemmas are consequences of Theorem 18.5.4 and 18.5.10 in Reference [24].
Lemma 2.
Let b S ( m , g ) where m is an admissible weight and g is defined by (8). Then there exists c S ( m , g ) such that Op ( b ) * = Op ( c ) and c ( x , ξ ) = b ( x , ξ ) ¯ + r ( x , ξ ) where the remainder r S ( λ 1 2 μ 1 m , g ) .
Lemma 3.
Let b S ( m 1 , g ) and c S ( m 2 , g ) where m j are admissible weights for j = 1 , 2 and g is defined by (8). Then,
(i)
there exists d S ( m 1 m 2 , g ) such that Op ( b ) Op ( c ) = Op ( d ) and d ( x , ξ ) = b ( x , ξ ) c ( x , ξ ) + r ( x , ξ ) where r S ( λ 1 2 μ 1 m 1 m 2 , g ) .
(ii)
for commutator i [ Op ( b ) , Op ( c ) ] = Op ( f ) , it holds that f S ( λ 1 2 μ 1 m 1 m 2 , g ) and f ( x , ξ ) = { b , c } ( x , ξ ) + r ( x , ξ ) where r S ( λ μ 2 m 1 m 2 , g ) .
The operators in S ( ν j μ k , g ) act on Sobolev spaces adapted to the class of symbol. Let b S ( ν j μ k , g ) , where μ and g are defined by (8). Then there exists C > 0 such that
Op ( b ) u L 2 ( R d ) C ν j Op ( μ k ) u L 2 ( R d ) , u S ( R d ) .
By symbolic calculus, the above estimate is equivalent to Op ( μ k ν j ) Op ( b ) acts on L 2 ( R d ) since the operators associated with symbol in S ( 1 , g ) act on L 2 ( R d ) . In particular, if b S ( ν j μ , g ) , then for any λ λ 0 , τ max { C λ , 1 } and u S ( R d ) , it holds that
Op ( b ) u L 2 ( R d ) C τ ν j u L 2 ( R d ) + C ν j D u L 2 ( R d ) ,
where C > 0 depends on λ 0 and C .

2.2. Commutator Estimate

In this subsection, we suppose that λ = 1 since the symbol does not depend on λ . The the metric in (8) becomes
g ˜ = d x 2 + μ 2 d ξ 2 , where μ is   defined   by   ( 8 ) .
To get the commutator estimate, we shall use the following Gårding inequality [24].
Lemma 4.
Let b S ( μ 2 k , g ˜ ) be real valued. μ and g ˜ are defined by (10). We assume there exists C > 0 such that b ( x , ξ , τ ) C μ 2 k . Then there exists C ˜ > 0 such that
Re ( Op ( b ) w | w ) C ˜ Op ( μ k ) w 2 , w S ( R d ) .
Define the operator P ( x , D , λ ) = D 2 + i λ D a ( x ) D λ 2 , where D 2 = j = 1 d D j 2 , D j = i x j and D a ( x ) D = j = 1 d D j a ( x ) D j . Introduce the weight function φ C ( R d ; R ) . The associated conjugate operator of P ( x , D , λ ) is P φ ( x , D , λ ) = e τ φ P ( x , D , λ ) e τ φ . Then,
P φ = ( D + i τ φ ( x ) ) 2 + i λ ( D + i τ φ ( x ) ) a ( x ) ( D + i τ φ ( x ) ) λ 2 .
By setting Q 2 = 1 2 ( P φ + P φ * ) and Q 1 = 1 2 i ( P φ P φ * ) , we have P φ = Q 2 + i Q 1 , where
Q 2 = D 2 τ 2 | φ | 2 λ τ D a ( x ) φ ( x ) λ τ φ ( x ) a ( x ) D λ 2 , Q 1 = τ D φ ( x ) + τ φ ( x ) D + λ D a ( x ) D λ τ 2 a ( x ) | φ | 2 .
Definition 3.
Let V be a bounded open set in R d . We say that the weight function φ C ( R d ; R ) satisfies the sub-ellipticitycondition in V ¯ if | d φ | > 0 in V ¯ and there exists constant C > 0 ,
p φ ( x , ξ , τ ) = 0 , ( x , ξ ) V ¯ × R d , τ > 0 { q 2 , q 1 } ( x , ξ , τ ) C ( | ξ | 2 + τ 2 ) 3 / 2 ,
where p φ ( x , ξ , τ ) = | ξ + i τ φ ( x ) | 2 = q 2 ( x , ξ , τ ) + i q 1 ( x , ξ , τ ) and q 1 , q 2 are real valued.
Lemma 5
([24]). Let V be a bounded open set in R d and ψ C ( R d ; R ) be such that | ψ | > 0 in V ¯ . Then, for γ > 0 sufficiently large, φ = e γ ψ fulfills the sub-ellipticity property in V.
Lemma 6.
Assume that φ satisfies the sub-ellipticity in Definition 3. For all w C 0 ( Ω ) , there exist C 1 , C 2 > 0 and τ 0 > 0 such that the following inequality holds for all τ τ 0 ,
C 1 τ 3 w 2 + C 1 τ D w 2 Op ( { ξ 2 τ 2 φ ( x ) 2 , 2 τ ξ · φ ( x ) } ) w | w + C 2 τ 1 Op ( ξ 2 τ 2 φ ( x ) 2 ) w 2 + C 2 τ 1 Op ( 2 τ ξ · φ ( x ) ) w 2 .
The proof of Lemma 6 can be found in Appendix A.

3. Carleman Estimate

In this section, we shall prove the Carleman estimate by the tool introduced by Hörmander [24,25] and called Weyl-Hörmander calculus. It allows us to define proper operators, class of symbol and pseudo-diff-calculus. Throughout this section, we denote by · V , ( · | · ) V the norm and inner product in L 2 ( V ) for V Ω , respectively.
Let the metric g and weight ν be defined by (8) and (9). Furthermore, due to the results in Section 2.1, we know that D + i τ φ ( x ) is an operator with a symbol in S ( μ , g ) class, λ a is in S ( ν , g ) , and ( 1 + i λ a ( x ) ) ( ξ + i τ φ ( x ) ) 2 , the principal symbol of P φ belongs to S ( ν μ 2 , g ) . It follows from Lemmas 2 and 3 that
P φ = Op ( 1 + i λ a ( x ) ) ( ξ + i τ φ ( x ) ) 2 λ 2 + R 3 Q 2 = Op ( q 2 ) λ 2 + R 2 Q 1 = Op ( q 1 ) + R 1 ,
where q 2 = | ξ | 2 τ 2 | φ ( x ) | 2 2 λ τ a ( x ) ξ · φ ( x ) , q 1 = 2 τ ξ · φ ( x ) + λ a ( x ) ( | ξ | 2 τ 2 | φ ( x ) | 2 ) and the symbols of R j are in S ( λ 1 2 ν μ , g ) for j = 1 , 2 , 3 . It is clear that
P φ v V 2 = Q 2 v V 2 + Q 1 v V 2 + 2 Re ( Q 2 v | i Q 1 v ) V .
In what follows, several Carleman estimates are introduced. First, we give an estimation on the subdomain which is far away from the boundary Γ .
Theorem 2.
Suppose φ satisfies sub-ellipticity condition on V Ω . Then, there exist positive constants C , K ˜ and λ 0 , such that for every u C 0 ( V ) , it holds
τ 3 ( 1 + λ 2 a ( x ) 2 ) 1 2 e τ φ u V 2 + τ ( 1 + λ 2 a ( x ) 2 ) 1 2 e τ φ D u V 2 C e τ φ P u V 2 ,
where λ λ 0 and τ max { K ˜ | λ | 5 / 4 , 1 } .
Proof. 
We shall prove that (17) is equivalent to
P φ v V 2 τ ν ( x ) D v V 2 + τ 3 ν ( x ) v V 2 .
First, assume (18) holds. Set v = e τ φ u . Then, D v = e τ φ ( D u i τ φ u ) and e τ φ D u = D v + i τ φ v . Then there exist positive constants c 1 , c 2 such that
c 1 ν ( x ) D v V + τ ν ( x ) v V ν ( x ) e τ φ D u V + τ ν ( x ) e τ φ u V c 2 ν ( x ) D v V + τ ν ( x ) v V .
Combining this with (18), we conclude that
τ 3 ν ( x ) e τ φ u V 2 + τ ν ( x ) e τ φ D u V 2 ( τ ν ( x ) D v V 2 + τ 3 ν ( x ) v V 2 ) e τ φ P u V 2 .
On the other hand, (17) implies that
τ 3 ν ( x ) e τ φ u V 2 + τ ν ( x ) e τ φ D u V 2 P φ v V 2 .
Then, we proved (18) from the above estimate and (19).
Now we are going to prove (18). Note that
2 Re ( Q 2 v | i Q 1 v ) V = ( Q 2 v | i Q 1 v ) V + ( i Q 1 v | Q 2 v ) V = ( i [ Q 2 , Q 1 ] v | v ) V ,
where we denote by [ O p ( a ) , O p ( b ) ] = O p ( a ) O p ( b ) O p ( b ) O p ( a ) and Poisson bracket { a , b } ( x , ξ , τ ) = 1 j d ( ξ j a x j b x j a ξ j b ) ( x , ξ , τ ) . From symbolic calculus, the principal symbol of i [ Q 2 , Q 1 ] is { q 2 , q 1 } . Due to the results in Section 2.1, we obtain that
i [ Q 2 , Q 1 ] = Op ( { q 2 , q 1 } ) + R 4 ,
where { q 2 , q 1 } S ( ν 2 μ 3 λ 1 2 , g ) and the symbol of R 4 is in S ( ν 2 μ 2 λ , g ) .
A direct computation gives that
{ q 2 , q 1 } = ( 1 + a 2 ( x ) λ 2 ) { ξ 2 τ 2 φ ( x ) 2 , 2 τ ξ . φ ( x ) } + { ξ 2 τ 2 φ ( x ) 2 , λ a ( x ) } λ a ( x ) { 2 τ ξ . φ ( x ) , λ a ( x ) } ( ξ 2 τ 2 φ ( x ) 2 ) 2 τ ξ . φ ( x ) λ a ( x ) { λ a ( x ) , ξ 2 τ 2 φ ( x ) 2 } + { λ a ( x ) , 2 τ ξ . φ ( x ) } .
Since
q 2 ( x , ξ ) + λ a ( x ) q 1 ( x , ξ ) = ( 1 + λ 2 a ( x ) 2 ) ( ξ 2 τ 2 φ ( x ) 2 ) , q 1 ( x , ξ ) λ a ( x ) q 2 ( x , ξ ) = 2 τ ξ . φ ( x ) ( 1 + λ 2 a 2 ( x ) ) ,
we have
{ q 2 , q 1 } = ν 2 ( x ) ξ 2 τ 2 φ ( x ) 2 , 2 τ ξ . φ ( x ) + ν 2 ( x ) ( q 2 ( x , ξ ) + λ a ( x ) q 1 ( x , ξ ) ) { ξ 2 τ 2 φ ( x ) 2 , λ a ( x ) } λ a ( x ) { 2 τ ξ . φ ( x ) , λ a ( x ) } ν 2 ( x ) ( q 1 ( x , ξ ) λ a ( x ) q 2 ( x , ξ ) ) λ a ( x ) { λ a ( x ) , ξ 2 τ 2 φ ( x ) 2 } + { λ a ( x ) , 2 τ ξ . φ ( x ) } .
Then, it follows from symbolic calculus that
i [ Q 2 , Q 1 ] = ν ( x ) Op { ξ 2 τ 2 φ ( x ) 2 , 2 τ ξ . φ ( x ) } ν ( x ) + B 1 ν 1 ( x ) Op q 2 ( x , ξ ) + λ a ( x ) q 1 ( x , ξ ) + λ 2 B 2 ν 1 ( x ) Op q 1 ( x , ξ ) λ a ( x ) q 2 ( x , ξ ) λ 3 a ( x ) + R 5 ,
where the symbol of B 1 is ν 1 ( x ) { ξ 2 τ 2 φ ( x ) 2 , λ a ( x ) } λ a ( x ) { 2 τ ξ . φ ( x ) , λ a ( x ) } S ( λ 1 2 ν μ , g ) , the symbol of B 2 is ν 1 ( x ) λ a ( x ) { λ a ( x ) , ξ 2 τ 2 φ ( x ) 2 } + { λ a ( x ) , 2 τ ξ . φ ( x ) } S ( λ 1 2 ν μ , g ) and the symbol of R 5 is in S ( λ ν 2 μ 2 , g ) . We refer to Section 2.1 where the rules on symbolic calculus are given and precise. Therefore, by the continuity of pseudo-differential operator, we have that for j = 1 , 2 , k = 0 , 1 and = 1 , 2 ,
| ( B j ν ( x ) 1 ( λ a ( x ) ) k Q v | v ) V | = | ( ν ( x ) 1 ( λ a ( x ) ) k Q v | B j * v ) | C Q v V B j * v V ( 1 / 10 ) Q v V 2 + C λ τ 2 ν ( x ) v V 2 + C λ ν ( x ) D v V 2 ,
| ( B j ν 1 ( x ) λ 2 ( λ a ( x ) ) k v | v ) V | C λ 5 / 2 τ ν ( x ) v V + ν ( x ) D v V ν ( x ) v V C λ 5 / 2 τ ν ( x ) v V 2 + C λ 5 / 2 τ 1 ν ( x ) D v V 2 ,
and
| ( R 5 v | v ) V | λ τ 2 ν ( x ) v V 2 + λ ν ( x ) D v V 2 .
Let w = ( 1 + λ 2 a ( x ) 2 ) 1 2 v = ν v in (14). We obtain
C 1 τ 3 ν v V 2 + C 1 τ D ( ν ( x ) v ) V 2 ( Op ( { ξ 2 τ 2 φ ( x ) 2 , 2 τ ξ . φ ( x ) } ) ν v | ν v ) V + C 2 τ 1 Op ( ξ 2 τ 2 φ ( x ) 2 ) ν v V 2 + C 2 τ 1 Op ( 2 τ ξ . φ ( x ) ) ν v V 2 , C 1 , C 2 > 0 .
Since the symbol of [ D , ν ( x ) ] is in S ( λ 1 2 ν , g ) , we have
ν ( x ) D v V 2 D ( ν ( x ) v ) V 2 + [ D , ν ( x ) ] v V 2 D ( ν ( x ) v ) V 2 + C λ ν ( x ) v V 2 .
Consequently, for τ C λ with C > 0 sufficiently large, it holds
τ 3 ν ( x ) v V 2 + τ ν ( x ) D v V 2 τ 3 ν ( x ) v V 2 + τ D ( ν ( x ) v ) V 2 .
It follows from (26) and (27) that
C 1 τ 3 ν v V 2 + τ ν ( x ) D v V 2 ( Op ( { ξ 2 τ 2 φ ( x ) 2 , 2 τ ξ . φ ( x ) } ) ν v | ν v ) V + C 2 τ 1 Op ( ξ 2 τ 2 φ ( x ) 2 ) ν v V 2 + C 2 τ 1 Op ( 2 τ ξ . φ ( x ) ) ν v V 2 .
Now we are going to estimate the terms Op ( ξ 2 τ 2 φ ( x ) 2 ) ν v V 2 and Op ( 2 τ ξ . φ ( x ) ) ν v V 2 . It follows from (15) and (21) that
Op ( ξ 2 τ 2 φ ( x ) 2 ) ν v V 2 + Op ( 2 τ ξ . φ ( x ) ) ν v V 2 = Op ( ν 2 ( q 2 + λ a ( x ) q 1 ) ν v V 2 + Op ( ν 2 ( q 1 λ a ( x ) q 2 ) ν v V 2 = ν 1 Q 2 + λ 2 R 2 + λ a ( x ) ( Q 1 R 1 ) v V 2 + ν 1 Q 1 R 1 λ a ( x ) ( Q 2 + λ 2 R 2 ) v V 2 .
Combining this with the fact that the symbols of R j are in S ( λ 1 2 ν μ , g ) for j = 1 , 2 , 3 , we have that there exists a positive constant C such that
Op ( ξ 2 τ 2 φ ( x ) 2 ) ν v V 2 + Op ( 2 τ ξ . φ ( x ) ) ν v V 2 j = 1 , 2 Q j v V 2 + R 1 v V 2 + R 2 v V 2 + λ 4 v V 2 j = 1 , 2 Q j v V 2 + λ τ 2 ν v V 2 + λ ν D v V 2 + λ 4 v V 2 .
From (28) and (30), we have that for τ C λ , there exist positive constants C 1 , C 2 such that
C 1 τ 3 ν ( x ) v V 2 + τ ν ( x ) D v V 2 ( Op ( { ξ 2 τ 2 φ ( x ) 2 , 2 τ ξ . φ ( x ) } ) ν v | ν v ) V + C 2 ( 1 / τ ) j = 1 , 2 Q j v V 2 + λ τ ν v V 2 + ν D v V 2 + λ 3 v V 2 .
On the other hand, due to (22)–(24), we have that there exists a positive constant C such that
Op ( { ξ 2 τ 2 φ ( x ) 2 , 2 τ ξ . φ ( x ) } ) ν v | ν v V ( i [ Q 2 , Q 1 ] v | v ) V + ( 1 / 4 ) j = 1 , 2 Q j v V 2 + C ( λ τ 2 ν ( x ) v V 2 + λ ν ( x ) D v V 2 + λ 5 / 2 τ ν ( x ) v V 2 + λ 5 / 2 τ 1 ν ( x ) D v V 2 ) .
Finally, by (31) and (32), one can choose τ K ˜ | λ | 5 / 4 with K ˜ sufficiently large and τ τ 0 with τ 0 sufficiently large, such that for some C > 0 ,
C τ 3 ν ( x ) v V 2 + τ ν ( x ) D v V 2 ( i [ Q 2 , Q 1 ] v | v ) V + ( 1 / 2 ) j = 1 , 2 Q j v V 2 . + ε τ 3 ν ( x ) v V 2 + ε τ ν ( x ) D v V 2 ,
where ε > 0 is arbitrary. Choosing ε small with respect C 1 , using (16), (20) and (33), we obtain
C 1 τ 3 ν ( x ) v V 2 + C 1 τ ν ( x ) D v V 2 + ( 1 / 2 ) Q 1 v V 2 + ( 1 / 2 ) Q 2 v V 2 ( i [ Q 2 , Q 1 ] v , v ) + Q 1 v V 2 + Q 2 v V 2 = P φ v V 2 ,
this implies (18). ☐
Remark 2.
The estimates on the previous terms impose the assumption τ K ˜ | λ | 5 / 4 . The other remainder terms only impose the condition τ C λ . This condition is related with the principal normal condition. Indeed for a complex operator, with symbol p 1 + i p 2 where p 1 , p 2 are both real valued, the Carleman estimate is only true if { p 1 , p 2 } = 0 on p 1 = p 2 = 0 . Here the symbol of operator before conjugaison by weight is ξ 2 λ 2 + i λ a ( x ) ξ 2 , and the Poisson bracket is { ξ 2 λ 2 , λ a ( x ) ξ 2 } = 2 λ ( ξ . a ( x ) ) ξ 2 . We can estimate this term, uniformly in a neighborhood of a ( x ) = 0 , by C λ a 1 2 ( x ) | ξ | 3 . This explanation does not justify the power | λ | 5 / 4 found at the end of computations but shows the difficulties.
Since there is higher order term d i v ( a ( x ) y t ) in system (1), it is necessary to deal with the term d i v ( a ( x ) f ) for f H 1 ( Ω ) when proving the resolvent estimate. The following result is analogue to the work by [26].
Theorem 3.
Suppose φ satisfies sub-ellipticity condition on V Ω . Then, there exist C , K ˜ > 0 , λ 0 > 0 , such that for all u C 0 ( V ) satisfying
P u = g 0 + j = 0 d x j g j , w h e r e g j L 2 ( V ) , j = 0 , 1 , , d ,
it holds
τ ( 1 + λ 2 a ( x ) 2 ) 1 2 e τ φ u L 2 ( V ) 2 + τ 1 ( 1 + λ 2 a ( x ) 2 ) 1 2 e τ φ D u L 2 ( V ) 2 C j = 0 d e τ φ g j L 2 ( V ) 2 ,
where λ λ 0 and τ max { K ˜ | λ | 5 / 4 , 1 } .
Proof. 
First, from (15), we have D 2 = Q 2 + S 2 and λ a D 2 = Q 1 + S 1 where S 1 and S 2 have symbols in S ( τ ν μ , g ) if τ λ . It follows that
D 2 v 2 Q 2 v 2 + τ 2 ( τ 2 ν ( x ) v 2 + ν ( x ) D v 2 ) , λ a D 2 v 2 Q 1 v 2 + τ 2 ( τ 2 ν ( x ) v 2 + ν ( x ) D v 2 ) .
Using the fact that ν ( x ) D 2 v 2 2 ( D 2 v 2 + λ a D 2 v 2 ) , (34), we obtain
τ 3 ν ( x ) v 2 + τ ν ( x ) D v 2 + τ 1 ν ( x ) D 2 v 2 P φ v 2 .
Let w and χ be in C 0 ( Ω ) such that χ = 1 on a neighborhood of supp w . From (27), we obtain
τ ν ( x ) w 2 + τ 1 ν ( x ) D w 2 τ ν ( x ) w 2 + τ 1 D ( ν ( x ) w ) 2 τ 3 Op ( μ 1 ) ν ( x ) χ w 2 + τ 1 D 2 Op ( μ 1 ) ν χ w 2 ,
where the last estimate is obtained by Fourier transform and by the inequality
τ + | ξ | 2 τ τ 3 τ 2 + | ξ | 2 + | ξ | 4 τ ( τ 2 + | ξ | 2 ) .
From the results in Section 2.1, we have Op ( μ 1 ) ν χ = ν χ Op ( μ 1 ) + R 1 , where R 1 has a symbol in S ( μ 2 ν λ 1 2 , g ) , and D 2 Op ( μ 1 ) ν χ = ν D 2 χ Op ( μ 1 ) + R 2 , where R 2 has a symbol in S ( ν λ 1 2 , g ) . Then, it follows from (39) that
τ ν ( x ) w 2 + τ 1 ν ( x ) D w 2 τ 3 ν ( x ) χ Op ( μ 1 ) w 2 + τ 1 ν D 2 χ Op ( μ 1 ) w 2 + λ ν ( x ) w 2 .
For τ max { C λ , 1 } with C large enough, one has the following result from the above inequality.
τ ν ( x ) w 2 + τ 1 ν ( x ) D w 2 τ 3 ν ( x ) χ Op ( μ 1 ) w 2 + τ 1 ν D 2 χ Op ( μ 1 ) w 2 .
Now, we apply (38) to v = χ Op ( μ 1 ) w to have
τ 3 ν ( x ) χ Op ( μ 1 ) w 2 + τ ν ( x ) D χ Op ( μ 1 ) w 2 + τ 1 ν ( x ) D 2 χ Op ( μ 1 ) w 2 P φ χ Op ( μ 1 ) w 2 .
Thus, combining this with (40) yields
τ ν ( x ) w 2 + τ 1 ν ( x ) D w 2 P φ χ Op ( μ 1 ) w 2 .
Finally, note that P φ has a symbol in S ( ν μ 2 , g ) . Consequently, P φ χ Op ( μ 1 ) = Op ( μ 1 ) P φ χ + R , where R has a symbol in S ( ν λ 1 2 , g ) . Then, we can deduce from (41) that
τ ν ( x ) w 2 + τ 1 ν ( x ) D w 2 Op ( μ 1 ) P φ w 2 + λ ν w 2 .
When τ C λ , with C large enough, the error term ν λ 1 2 w 2 can be absorbed by the left hand side. Hence,
τ ν ( x ) w 2 + τ 1 ν ( x ) D w 2 Op ( μ 1 ) P φ w 2 .
For w = e τ φ u , we have
P φ w = e τ φ P u = e τ φ g 0 + j = 1 d e τ φ x j g j = e τ φ g 0 + j = 1 d x j ( e τ φ g j ) τ g j x j φ .
Obviously, one has that Op ( μ 1 ) P φ w 2 j = 0 d e τ φ g j 2 . Combining this with (42), we obtain Theorem 3. ☐
Remark 3.
Since a ( · ) is non-negative and not identically null, there exists δ > 0 such that { x Ω : a ( x ) > δ } . We introduce several sets as follows.
W 1 = Ω supp a , W 2 = Ω { x Ω : a ( x ) > δ } O ( Γ ) , W 3 = Ω ¯ O { x Ω : a ( x ) 0 } , W 4 = Ω ¯ { x Ω : a ( x ) > δ / 2 } ,
where O ( Γ ) means the neighborhood of Γ .
It is known that there exists a function ψ C ( Ω ) such that [25]
(1)
ψ ( x ) = 0 for x Ω .
(2)
n ψ ( x ) < 0 for x Ω .
(3)
ψ ( x ) 0 for x W 1 .
Let φ = e γ ψ . It follows from Lemma 5 that φ satisfies the sub-ellipticity condition on x W 2 if γ > 0 is sufficiently large. Then, in Theorem 2 and 3, one can choose V as W 2 .
The following result can be obtained from the classical Carleman estimate up to the boundary [27] (Proposition 2). Note that this estimate corresponds to the Laplacian with Dirichlet boundary condition.
Lemma 7.
Suppose φ satisfies the sub-ellipticity condition on W 4 . Then, there exist C , K > 0 , λ 0 > 0 , such that for all u C ( Ω ¯ ) satisfying supp u W 3 and u = 0 on Γ, it holds that
τ 3 e τ φ u 2 + τ e τ φ D u 2 C e τ φ ( D 2 λ 2 ) u 2 ,
where λ λ 0 and τ max { K λ , 1 } .
Theorem 4.
Suppose φ satisfies sub-ellipticity condition on W 4 . Let u C ( Ω ¯ ) satisfy
P u = f 0 + j = 1 d x j f j i n Ω , u = 0 o n Γ ,
where f j L 2 ( Ω ) , supp f 0 Ω and supp f j { x Ω : a ( x ) 0 } for j = 1 , , d . Then, there exist C , K ˜ > 0 , λ 0 > 0 , such that for all λ λ 0 and τ max { K ˜ | λ | 5 / 4 , 1 } , it holds
τ e τ φ u 2 + τ 1 e τ φ D u 2 C j = 0 d e τ φ f j 2 + λ e τ φ u L 2 ( { x Ω , a ( x ) δ / 2 } ) 2 ,
where the positive constant δ is defined as in Remark 3.
Proof. 
Let χ 1 , χ 2 C ( Ω ) satisfy χ 1 and χ 2 are supported on W 3 and W 2 , respectively. We assume that χ 1 + χ 2 = 1 on W 3 . In particular χ 1 = 1 on a neighborhood of Ω , and χ 2 = 1 in a neighborhood of the frontier of supp a .
Since χ 1 x j f j = 0 for j = 1 , , d , we have P χ 1 u = χ 1 f 0 + [ P , χ 1 ] u , where [ P , χ 1 ] is a first order operator and supported on supp χ 1 . By using χ 1 u instead of u in Lemma 7, we obtain
τ 3 e τ φ χ 1 u 2 + τ e τ φ D χ 1 u 2 e τ φ χ 1 f 0 2 + e τ φ [ P , χ 1 ] u 2 .
On the other hand, it is clear that P χ 2 u = χ 2 f 0 + j = 1 d x j ( χ 2 f j ) f j x j χ 2 + [ P , χ 2 ] u . Since [ P , χ 2 ] is a first order operator, we have that there exist a 0 , a 1 , , a d and b 0 , b 1 , , b d such that [ P , χ 2 ] u = a 0 u + j = 1 d x j ( a j u ) + λ j = 1 d x j ( b j u ) , where a j and b j ( j = 0 , 1 , , d ) are supported on the set { x Ω : δ a ( x ) δ / 2 } ( Ω supp a ) . Then, applying Theorem 3 with χ 2 u instead of u, we obtain
τ ( 1 + λ 2 a ( x ) 2 ) 1 2 e τ φ χ 2 u 2 + τ 1 ( 1 + λ 2 a ( x ) 2 ) 1 2 e τ φ D χ 2 u 2 j = 0 d e τ φ f j 2 + e τ φ u L 2 ( Ω supp a ) 2 + λ e τ φ u L 2 ( { x Ω : a ( x ) δ / 2 } ) 2 .
Summing this estimate with (44) multiply by τ 2 , using ( 1 + λ 2 a ( x ) 2 ) 1 , we obtain
τ e τ φ ( χ 1 + χ 2 ) u 2 + τ 1 e τ φ D ( χ 1 + χ 2 ) u 2 j = 0 d e τ φ f j 2 + e τ φ u L 2 ( Ω supp a ) 2 + λ e τ φ u L 2 ( { x Ω , a ( x ) δ / 2 } ) 2 + τ 2 e τ φ [ P , χ 1 ] u 2 .
As D ( χ 1 + χ 2 ) is supported on { x Ω : a ( x ) δ / 2 } , χ 1 + χ 2 = 1 on Ω supp a , and [ P , χ 1 ] supported where χ 1 + χ 2 = 1 , we have, for τ sufficiently large
τ e τ φ ( χ 1 + χ 2 ) u 2 + τ 1 e τ φ ( χ 1 + χ 2 ) D u 2 j = 0 d e τ φ f j 2 + λ e τ φ u L 2 ( { x Ω : a ( x ) δ / 2 } ) 2 .
Then as χ 1 + χ 2 = 1 on Ω { x Ω : a ( x ) > δ / 2 } , we obtain the statement of Theorem 4. ☐

4. Resolvent Estimate

In this section, we prove the main result. From Batty-Duyckaerts [17], the estimate of the energy decay in Theorem 1 can be obtained through the following result.
Theorem 5.
There exists C > 0 , such that for every ν R with λ large, we have
( A i λ ) 1 L ( H ) C e C | λ | 5 / 4 .
Let λ be a real number such that | λ | is large. Consider the resolvent equation:
F = ( A i λ ) X , where X = y 1 , y 2 D ( A ) , F = f 1 , f 2 H ,
which yields
div ( y 1 + a y 2 ) + λ 2 y 1 = i λ f 1 + f 2 , in Ω , y 2 = i λ y 1 + f 1 in Ω , y 1 | Γ = 0 .
In what follows, we shall prove the solution to (46) satisfies Theorem 5. As a result, the logarithmic decay of (7) is arrived at. The proof is divided into two lemmas.
Lemma 8.
Let η > 0 and χ C 0 ( { x Ω : a ( x ) > η } ) be real valued functions, there exists C > 0 such that
χ y 1 C ( y 1 L 2 ( { x Ω : a ( x ) > η } ) + f 1 H 1 ( Ω ) + f 2 L 2 ( Ω ) ) ,
where y 1 , satisfies (47).
Proof. 
By (47), one has that y 1 satisfies
div ( y 1 + i λ a y 1 ) + λ 2 y 1 = i λ f 1 + f 2 div a f 1 .
Multiplying this equation by χ 2 y ¯ 1 , and integrating on Ω , we obtain,
Ω y 1 ( χ 2 y ¯ 1 ) d x i λ Ω a ( x ) y 1 ( χ 2 y ¯ 1 ) d x + λ 2 χ y 1 2 = i λ Ω f 1 χ 2 y ¯ 1 d x + Ω f 2 χ 2 y ¯ 1 d x + Ω a ( x ) f 1 ( χ 2 y ¯ 1 ) d x .
Since ( χ 2 y ¯ 1 ) = χ 2 y ¯ 1 + 2 y ¯ 1 χ χ , we have
| Ω a ( x ) y 1 ( χ 2 y ¯ 1 ) d x | D y 1 L 2 ( { x Ω : a ( x ) > η } ) 2 + χ y 1 D y 1 L 2 ( { x Ω : a ( x ) > η } ) 2 D y 1 L 2 ( { x Ω : a ( x ) > η } ) 2 + 1 / 8 χ y 1 2 .
Substituting (50) into (49), we obtain
λ 2 χ y 1 2 ( 2 | λ | + 1 ) D y 1 L 2 ( { x Ω : a ( x ) > η } ) 2 + λ 2 / 8 χ y 1 2 + λ 2 / 4 χ y 1 2 + f 1 2 + λ 2 f 2 2 + 2 ( 1 + λ 2 ) f 1 H 1 ( Ω ) 2 + λ 2 / 8 χ y 1 2 + D y 1 L 2 ( { x Ω : a ( x ) > η } ) 2 .
Then, the proof is finished. ☐
Lemma 9.
For every y 1 , f 1 and f 2 satisfying (47), there exists C > 0 such that
Ω a | y 1 | 2 d x C f 1 H 1 ( Ω ) 2 + ( f 1 + f 2 ) y 1 .
Proof. 
From (47), multiplying by y ¯ 1 , we have
Ω div ( y 1 + a y 2 ) + λ 2 y 1 y ¯ 1 d x = Ω ( i λ f 1 + f 2 ) y ¯ 1 d x .
Replacing y 2 by f 1 + i λ y 1 and integrating by parts, we have
Ω ( y 1 · y ¯ 1 + λ 2 | y 1 | 2 a ( i λ y 1 + f 1 ) · y ¯ 1 ) d x = Ω ( i λ f 1 + f 2 ) y ¯ 1 d x .
Consequently,
Ω [ ( i λ a + 1 ) | y 1 | 2 + λ 2 | y 1 | 2 ] d x = Ω [ a f 1 · y ¯ 1 + ( i λ f 1 + f 2 ) y ¯ 1 ] d x .
Taking the imaginary part, we obtain
Ω a | y 1 | 2 d x = λ 1 Im Ω [ a f 1 · y 1 ¯ + ( i λ f 1 + f 2 ) y 1 ¯ ] d x .
Using the Cauchy-Schwarz inequality in the above equality, one has
Ω a | y 1 | 2 d x | λ | 1 Ω a | y 1 | 2 d x 1 2 f 1 H 1 ( Ω ) + | λ | 1 ( | λ | f 1 + f 2 ) y 1 1 / 4 Ω a | y 1 | 2 d x + | λ | 2 f 1 H 1 ( Ω ) 2 + | λ | 1 ( | λ | f 1 + f 2 ) y 1 .
Thus, we obtain the result. ☐
Proof of Theorem 5.
In Lemma 8, taking χ ( x ) = 1 is a ( x ) 2 η , we have
y 1 L 2 ( { x Ω : a ( x ) > 2 η } ) 2 C ( y 1 L 2 ( { x Ω : a ( x ) > η } ) 2 + f 1 H 1 ( Ω ) 2 + f 2 L 2 ( Ω ) 2 ) .
It follows from Lemma 9 that
y 1 L 2 ( { x Ω : a ( x ) > η } ) 2 C f 1 H 1 ( Ω ) 2 + ( f 1 + f 2 ) y 1 .
Both inequalities imply
y 1 H 1 ( { x Ω : a ( x ) > 2 η } ) 2 C f 1 H 1 ( Ω ) 2 + f 2 L 2 ( Ω ) 2 + ( f 1 + f 2 ) y 1 .
Applying Theorem 4 to y 1 satisfying (48), we obtain
τ e τ φ y 1 2 + τ 1 e τ φ D y 1 2 λ 2 e τ φ f 1 2 + e τ φ f 2 2 + j = 1 d e τ φ a x j f 1 2 + λ e τ φ y 1 L 2 ( { x Ω : a ( x ) δ / 2 } ) 2 .
Let c 1 = min x Ω φ ( x ) and c 2 = max x Ω φ ( x ) , we conclude from the above inequality that
e 2 c 1 τ τ y 1 2 + τ 1 e 2 c 1 τ D y 1 2 λ 2 e 2 c 2 τ f 1 H 1 ( Ω ) 2 + e 2 c 2 τ f 2 2 + λ e 2 c 2 τ y 1 H 1 ( { x Ω : a ( x ) δ / 2 } ) 2 .
Substituting (52) with η = δ / 2 into (53), we obtain
e 2 c 1 τ τ y 1 2 + τ 1 e 2 c 1 τ D y 1 2 λ 2 e 2 c 2 τ f 1 H 1 ( Ω ) 2 + e 2 c 2 τ f 2 2 + e 2 c 2 τ ( f 1 + f 2 ) y 1 .
Let c 3 = 2 ( c 2 c 1 ) + 1 . For τ max { K ˜ | λ | 5 / 4 , 1 } , One has
y 1 2 + D y 1 2 e c 3 τ f 1 H 1 ( Ω ) 2 + e c 3 τ f 2 2 + e c 3 τ ( f 1 + f 2 ) y 1 .
Using e c 3 τ ( f 1 + f 2 ) y 1 ε y 1 2 + ε 1 e 2 c 3 τ ( f 1 + f 2 ) 2 in the above estimate, we conclude that
y 1 2 + D y 1 2 e 2 c 3 τ f 1 H 1 ( Ω ) 2 + e 2 c 3 τ f 2 2 ,
which gives the desired result taking τ = K | λ | 5 / 4 with K large enough. ☐

5. Conclusions

It is known that the stability property of the wave equation system with local viscoelastic damping depends on both continuousness and geometry of the support set of the damping function. In this paper, we obtain the logarithmic decay of a wave equation system with local Kelvin-Voigt damping, where the damping function is smooth and its support is an arbitrary non-empty subset. The approach is based on Batty-Duyckaerts’ result that the resolvent estimate (45) implies the logarithmic decay of the semigroup. Through introducing proper operators, class of symbol and pseudo-diff-calculus, we obtain a Carleman estimate on the subdomain far away from the boundary. Combining these with a classical Carleman estimate up to the boundary, we arrive at the desired Carleman estimate and resolvent estimate.

Author Contributions

Investigation, Q.Z.; Methodology, L.R. The contributions of all authors are equal. Both authors have read and agreed to the published version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Beijing Natural Science Foundation (grant No. 4182059) and National Natural Science Foundation of China (grants No. 61873036).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this Appendix, we shall prove Lemmas 1 and 6. First, we claim that for a compactly supported and nonnegative function a C 2 ( Ω ) , the following inequality holds:
| a ( x ) | 2 2 a ( x ) a , x Ω .
In fact, from the following identity
a ( x + h ) = a ( x ) + a ( x ) h + 0 1 ( 1 t ) a ( x + t h ) h 2 d t , h R ,
one can get
a ( x ) + a ( x ) h + 1 2 a | h | 2 0 .
Let h = y a ( x ) , where y R and x Ω are arbitrary. It follows from the above inequality that
a ( x ) + | a ( x ) | 2 y + 1 2 a | a ( x ) | 2 y 2 0 , x Ω , y R .
Then,
| a ( x ) | 4 2 a ( x ) a | a ( x ) | 2 0 ,
and (A1) is proved.
Proof of Lemma 1.
(i) From Definition 18.4.1 in [24], the metric g x , ξ defined by (8) is slowly varying if there exist δ > 0 and C > 0 such that
g x , ξ ( y x , η ξ ) δ implies g y , η ( X , Ξ ) C g x , ξ ( X , Ξ ) , x , y , ξ , η , X , Ξ R d ,
where the constants δ and C are independent on the parameters λ and τ .
Suppose 0 < δ 1 / 4 and
g x , ξ ( y x , η ξ ) = λ | y x | 2 + ( τ 2 + | ξ | 2 ) 1 | η ξ | 2 δ .
Then, we have
τ 2 + | ξ | 2 τ 2 + 2 | ξ η | 2 + 2 | η | 2 τ 2 + 2 δ ( τ 2 + | ξ | 2 ) + 2 | η | 2 .
This implies that τ 2 + | ξ | 2 4 ( τ 2 + | η | 2 ) . Consequently,
g y , η ( X , Ξ ) = λ | X | 2 + ( τ 2 + | η | 2 ) 1 | Ξ | 2 λ | X | 2 + 1 4 ( τ 2 + | ξ | 2 ) 1 | Ξ | 2 g x , ξ ( X , Ξ ) .
Therefore, g is slowly varying.
For a given metric g x , ξ , the associated metric g x , ξ σ is defined by g x , ξ σ = ( τ 2 + | ξ | 2 ) d x 2 + λ 1 d ξ 2 . The metric g x , ξ is temperate if there exist C > 0 and N > 0 , such that
g x , ξ ( X , Ξ ) C g y , η ( X , Ξ ) 1 + g x , ξ σ ( x y , ξ η ) N , x , y , ξ , η , X , Ξ R d ,
where the constants C and N are independent on the parameters λ and τ (Definition 18.5.1 in [24]).
For the metric g = g x , ξ defined by (8), (A2) is equivalent to
λ | X | 2 + ( τ 2 + | ξ | 2 ) 1 | Ξ | 2 C λ | X | 2 + ( τ 2 + | η | 2 ) 1 | Ξ | 2 1 + ( τ 2 + | ξ | 2 ) | x y | 2 + λ 1 | ξ η | 2 N .
First, assume that τ 2 + | η | 2 4 ( τ 2 + | ξ | 2 ) . It follows that
( τ 2 + | η | 2 ) C ( τ 2 + | ξ | 2 ) 1 + λ 1 | ξ η | 2 N , C > 0 , N > 0 .
Then it is easy to obtain (A3) from (A4).
Secondly, consider the case τ 2 + | η | 2 > 4 ( τ 2 + | ξ | 2 ) . Then
| η | > 2 | ξ | , | η | > 3 τ ,
and
| ξ η | > 1 2 | η | > 3 2 τ > 3 2 C λ .
It follows from (A5) and (A6) that
λ 1 | ξ η | 2 > 3 2 C | ξ η | > 3 4 C | η | .
Consequently,
1 + λ 1 | ξ η | 2 2 > 3 16 C 2 | η | 2 > 3 32 C 2 ( | η | 2 + 3 τ 2 ) .
This together with τ 2 + | ξ | 2 1 yields that there exists a positive constant C such that (A4) holds with N = 2 .
(ii) It is known from Definition 18.4.2 in [24] that a weight ν ( x ) is g-continuous if there exist δ > 0 and C > 0 such that
g x , ξ ( y x , η ξ ) δ implies C 1 ν ( x ) ν ( y ) C ν ( x ) , x , y , ξ , η R d .
where the constants δ and C are independent on the parameters λ and τ . Since the weight ν ( x ) defined by (9) does not depend on ξ , the above condition is reduced to
λ | x y | 2 δ implies C 1 ν ( x ) ν ( y ) C ν ( x ) , x , y R d .
The weight ν ( x ) is g-temperate if there exist C > 0 and N > 0 such that
ν ( y ) C ν ( x ) 1 + g y , η σ ( x y , ξ η ) N , x , y , ξ , η R d ,
where the constants C and N do not depend on the parameters λ and τ (18.5.1 in [24]). The weight ν ( x ) is admissible if it is g-continuous and g-temperate. When a weight is admissible, all the powers of this weight are g-continuous and g-temperate. Therefore, it suffice to prove that 1 + λ a ( x ) is admissible.
Let s [ 0 , t ] and t [ 0 , 1 ] . Define f ( s ) = λ a ( x + s ( y x ) ) and F ( t ) = sup s [ 0 , t ] f ( s ) where x , y Ω satisfying λ | x y | 2 δ . It is clear that f ( s ) = λ a ( x + s ( y x ) ) ( y x ) . Combining this with (A1) yields
| f ( s ) | λ | a ( x + s ( y x ) ) | | y x | 2 λ a 1 2 [ a ( x + s ( y x ) ) ] 1 2 | y x | .
Consequently,
sup s [ 0 , t ] | f ( s ) | 2 λ 1 2 a 1 2 F ( t ) 1 2 | y x | .
Since f ( t ) f ( 0 ) + t sup s [ 0 , t ] | f ( s ) | , F is non-decreasing and λ | x y | 2 δ , we obtain that for all t [ 0 , 1 ] ,
f ( t ) f ( 0 ) + C λ 1 2 F ( t ) 1 2 | y x | f ( 0 ) + C δ F ( t ) 1 2 f ( 0 ) + C δ F ( α ) 1 2 ,
where C = 2 a 1 2 and α [ t , 1 ] . Note that f ( 0 ) = F ( 0 ) . It follows that
F ( α ) = sup t [ 0 , α ] f ( t ) F ( 0 ) + C δ F ( α ) 1 2 .
This yields
1 + F ( α ) 1 + F ( 0 ) + C δ 1 + F ( α ) 1 2 1 + F ( 0 ) + C δ 1 + F ( α ) .
By choosing δ sufficiently small such that C δ 1 / 2 , one can deduce from (A8) that
1 + F ( α ) 2 ( 1 + F ( 0 ) ) , α [ t , 1 ] .
In particular, we have
1 + λ a ( y ) 2 ( 1 + λ a ( x ) ) .
The above inequality remains true if we exchange x and y. Therefore, the weight 1 + λ a ( x ) is g-continuous.
On the other hand, note that 1 + λ a ( x ) is independent to ξ . Then, to obtain the weight 1 + λ a ( x ) is σ -temperate, it is sufficient to prove that
1 + λ a ( y ) C ( 1 + λ a ( x ) ) ( 1 + τ 2 | x y | 2 ) N .
In fact, it is clear that 1 + λ a ( y ) 1 + λ ( a ( x ) + C | x y | ) where C = a . Therefore, there exists positive constant C = C C 1 such that
1 + λ a ( y ) ( 1 + λ a ( x ) ) ( 1 + C τ | x y | ) ( 1 + λ a ( x ) ) ( 2 + 2 ( C τ | x y | ) 2 ) 1 2 .
Thus, we obtain (A9) with N = 1 2 , C = 2 max { 1 , C } . ☐
Proof of Lemma 6.
In what follows, we use the symbolic calculus with λ = 1 since the symbol does not depend on λ . First, by homogeneity in ( ξ , τ ) , compactness arguments and sub-ellipticity condition, we claim that there exist constants C , δ > 0 such that
C | 2 ξ · φ ( x ) | 2 + μ 2 ( | ξ | 2 τ 2 | φ ( x ) | 2 ) 2 + | ξ | 2 τ 2 | φ ( x ) | 2 , 2 ξ · φ ( x ) δ μ 2 .
The proof of (A10) is classical. In fact, set
K = { ( x , ξ , τ ) R d × R d × R : x Ω ¯ , | ξ | 2 + τ 2 = 1 , τ 0 } ,
and for ( x , ξ , τ ) K , κ > 0 ,
G ( x , ξ , τ , κ ) = κ | 2 ξ · φ ( x ) | 2 + μ 2 ( | ξ | 2 τ 2 | φ ( x ) | 2 ) 2 + { ξ 2 τ 2 | φ ( x ) | 2 , 2 ξ · φ ( x ) } .
If p φ = 0 for ( x , ξ , τ ) K , then | 2 ξ · φ ( x ) | 2 + μ 2 ( | ξ | 2 τ 2 | φ ( x ) | 2 ) 2 = 0 . it is clear that there exists a positive constant δ such that (A10) holds due to the fact that ϕ is sub-elliptic. When | 2 ξ · φ ( x ) | 2 + μ 2 ( | ξ | 2 τ 2 | φ ( x ) | 2 ) 2 > 0 , there exists a positive constant κ x , ξ , τ such that G ( x , ξ , τ , κ ) > 0 for every κ κ x , ξ , τ since { ξ 2 τ 2 | φ ( x ) | 2 , 2 ξ · φ ( x ) } is bounded on K . By continuity of G ( x , ξ , τ , κ ) , there exists a neighborhood of ( x , ξ , τ ) , denoted by V x , ξ , τ , such that G ( x , ξ , τ , κ ) > 0 for all ( x , ξ , τ ) V x , ξ , τ and κ κ x , ξ , τ . Since K is compact, there exist finite sets V j = V x j , ξ j , τ j and corresponding constants κ j = κ x j , ξ j , τ j ( i = 1 , 2 , , n ) , such that K i = 1 n V j and G ( x , ξ , τ , κ ) > 0 for all ( x , ξ , τ ) V j and κ > κ j . Let κ ˜ = max { κ j : j = 1 , 2 , , n } . It follows that G ( x , ξ , τ , κ ˜ ) > 0 for all ( x , ξ , τ ) K and κ κ ˜ . Finally, using the compactness of K again, we conclude that there exists δ > 0 such that G ( x , ξ , τ , κ ˜ ) δ . Thus, (A10) is reached since g is a homogeneous function of degree 2 with respect to variables ( ξ , τ ) .
By Gårding inequality (11), there exists a constant C ˜ > 0 such that, for τ τ 0 where τ 0 sufficiently large,
C ˜ Op ( μ ) w 2 Re ( Op ( | 2 ξ · φ ( x ) | 2 + μ 2 ( | ξ | 2 τ 2 | φ ( x ) | 2 ) 2 + | ξ | 2 τ 2 | φ ( x ) | 2 , 2 ξ · φ ( x ) ) w | w ) .
Now we are going to estimate the terms τ Op ( | 2 ξ · φ ( x ) | 2 ) and τ Op ( μ 2 ( | ξ | 2 τ 2 | φ ( x ) | 2 ) 2 ) . Firstly, it follows from Lemma 3 that
τ 1 Op ( | 2 τ ξ · φ ( x ) | 2 ) = τ 1 Op ( 2 τ ξ · φ ( x ) ) * Op ( 2 τ ξ · φ ( x ) ) + τ Op ( r 1 ) ,
where r 1 S ( μ , g ˜ ) and g ˜ is defined by (10). Therefore, for any ε > 0 , there exists a positive constant C ε such that
| τ 1 Op ( | 2 τ ξ · φ ( x ) | 2 ) w | w | τ 1 Op ( 2 τ ξ · φ ( x ) ) w 2 + τ | ( Op ( r 1 ) w | w ) | τ 1 Op ( 2 τ ξ · φ ( x ) ) w 2 + ε τ Op ( μ ) w 2 + C ε τ w 2 .
Substituting (A13) into (A11) and choosing ε small enough, we have
C τ Op ( μ ) w 2 Re Op τ μ 2 ( | ξ | 2 τ 2 | φ ( x ) | 2 ) 2 + τ | ξ | 2 τ 2 | φ ( x ) | 2 , 2 ξ · φ ( x ) w | w + τ 1 Op ( 2 τ ξ · φ ( x ) ) w 2 + C ε τ w 2 .
Secondly, by symbolic calculus, we have that
τ Op ( μ 2 ( | ξ | 2 τ 2 | φ ( x ) | 2 ) 2 ) = τ Op ( r 0 ) Op ( | ξ | 2 τ 2 | φ ( x ) | 2 ) + τ Op ( r 2 ) ,
where r 0 ( x , ξ ) = μ 2 ( | ξ | 2 τ 2 | φ ( x ) | 2 ) S ( 1 , g ˜ ) and r 2 S ( μ , g ˜ ) . Therefore, for all ε > 0 , there exists C ε > 0 such that
| ( τ Op ( r 0 ) Op ( | ξ | 2 τ 2 | φ ( x ) | 2 ) w | w ) | C ε τ 1 Op ( r 0 ) Op ( | ξ | 2 τ 2 | φ ( x ) | 2 ) w 2 + ε τ 3 w 2 C ε τ 1 Op ( | ξ | 2 τ 2 | φ ( x ) | 2 ) w 2 + ε τ 3 w 2 .
We choose ε small enough and combine (A15) and (A16) with (A14) to get
C τ Op ( μ ) w 2 Re τ | ξ | 2 τ 2 | φ ( x ) | 2 , 2 ξ · φ ( x ) w | w ) + C ε τ 1 Op ( | ξ | 2 τ 2 | φ ( x ) | 2 ) w 2 + τ 1 Op ( 2 τ ξ · φ ( x ) ) w 2 + C ε ( τ + ε τ 3 ) w 2 .
Finally, it is clear that there exist positive constant C such that
τ 3 w 2 + τ D w 2 C τ Op ( μ ) w 2 .
Thus, we obtain (14) by using (A17) and (A18), choosing ε small enough and letting τ > τ 0 big enough. ☐

References

  1. Chen, S.; Liu, K.; Liu, Z. Spectrum and stability for elastic systems with global or local Kelvin-Voigt damping. SIAM J. Appl. Math. 1998, 59, 651–668. [Google Scholar]
  2. Liu, K.; Rao, B. Exponential stability for the wave equation with local Kelvin-Voigt damping. Z. Angew. Math. Phys. 2006, 57, 419–432. [Google Scholar] [CrossRef]
  3. Pazy, A. Semigroups of Linear Operators and Applications to Partial Differential Equations; Springer: New York, NY, USA, 1983. [Google Scholar]
  4. Huang, F. On the mathematical model for linear elastic systems with analytic damping. SIAM J. Control Optim. 1988, 26, 714–724. [Google Scholar] [CrossRef]
  5. Zhang, Q. On the lack of exponential stability for an elastic-viscoelastic waves interaction system. Nonlinear Anal. Real World Appl. 2017, 37, 387–411. [Google Scholar] [CrossRef]
  6. Duyckaerts, T. Optimal decay rates of the energy of a hyperbolic-parabolic system coupled by an interface. Asymptot. Anal. 2007, 51, 17–45. [Google Scholar]
  7. Rauch, J.; Zhang, X.; Zuazua, E. Polynomial decay for a hyperbolic-parabolic coupled system. J. Math. Pures Appl. 2005, 84, 407–470. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, X.; Zuazua, E. Long time behavior of a coupled heat-wave system arising in fluid-structure interaction. Arch. Rat. Mech. Anal. 2007, 184, 49–120. [Google Scholar] [CrossRef]
  9. Bardos, C.; Lebeau, G.; Rauch, J. Sharp sufficient conditions for the observation, control and stabilization of waves from the boundary. SIAM J. Control Optim. 1992, 30, 1024–1065. [Google Scholar] [CrossRef]
  10. Liu, K.; Liu, Z.; Zhang, Q. Eventual differentiability of a string with local Kelvin-Voigt damping. ESAIM Control. Optim. Caculus Var. 2017, 23, 443–454. [Google Scholar] [CrossRef]
  11. Liu, Z.; Zhang, Q. Stability of a string with local Kelvin-Voigt damping and non-smooth coefficient at interface. SIAM J. Control. Optim. 2016, 54, 1859–1871. [Google Scholar] [CrossRef]
  12. Renardy, M. On localized Kelvin-Voigt damping. Z. Angew. Math. Mech. 2004, 84, 280–283. [Google Scholar] [CrossRef]
  13. Zhang, Q. Exponential stability of an elastic string with local Kelvin-Voigt damping. Z. Angew. Math. Phys. 2010, 6, 1009–1015. [Google Scholar] [CrossRef]
  14. Liu, Z.; Rao, B. Frequency domain characterization of rational decay rate for solution of linear evolution euqations. Z. Angew. Math. Phys. 2005, 56, 630–644. [Google Scholar] [CrossRef] [Green Version]
  15. Zhang, Q. Polynomial decay of an elastic/viscoelastic waves interaction system. Z. Angew. Math. Phys. 2018, 69, 88. [Google Scholar] [CrossRef]
  16. Burq, N. Decays for Kelvin-Voigt damped wave equations I: The black box perturbative method. submitted.
  17. Batty, C.; Duyckaerts, T. Non-uniform stability for bounded semi-groups on Banach spaces. J. Evol. Equ. 2008, 8, 765–780. [Google Scholar] [CrossRef]
  18. Bellassoued, M. Distribution of resonances and decay rate of the local energy for the elastic wave equation. Commun. Math. Phys. 2000, 215, 375–408. [Google Scholar] [CrossRef]
  19. Bellassoued, M. Carleman estimates and distribution of resonances for the transparent obstacle and application to the stabilization. Asymptot. Anal. 2003, 35, 257–279. [Google Scholar]
  20. Bellassoued, M. Decay of solutions of the elastic wave equation with a localized dissipation. Ann. Fac. Sci. Toulouse Math. 2003, 12, 267–301. [Google Scholar] [CrossRef]
  21. Lebeau, G. Équation des ondes amortizes. In Algebraic and Geometric Methods in Mathematical Physics (Kaciveli, 1993), Volume 19 of Mathematical Physics Studies; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1996. [Google Scholar]
  22. Robbiano, L. Fonction de coût et contrôle des solutions des équations hyperboliques. Asymptotic Anal. 1995, 10, 95–115. [Google Scholar] [CrossRef]
  23. Lebeau, G.; Robbiano, L. Stabilisation de l’équation des ondes par le bord. Duke Math. J. 1997, 86, 465–491. [Google Scholar] [CrossRef]
  24. Hörmander, L. The Analysis of Linear Partial Differential Operators; Springer: Berlin, Germany, 2003. [Google Scholar]
  25. Le Rousseau, J.; Robbiano, L. Carleman estimate for elliptic operators with coefficients with jumps at an interface in arbitrary dimension and application to the null controllability of linear parabolic equations. Arch. Ration. Mech. Anal. 2010, 195, 953–990. [Google Scholar] [CrossRef] [Green Version]
  26. Imanuvilov, O.Y.; Puel, J.-P. Global Carleman estimates for weak solutions of elliptic nonhomogeneous Dirichlet problems. Int. Math. Res. Not. 2003, 16, 883–913. [Google Scholar] [CrossRef]
  27. Lebeau, G.; Robbiano, L. Contrôle exact de l’équation de la chaleur. Comm. Partial. Differ. Equ. 1995, 20, 335–356. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Robbiano, L.; Zhang, Q. Logarithmic Decay of Wave Equation with Kelvin-Voigt Damping. Mathematics 2020, 8, 715. https://doi.org/10.3390/math8050715

AMA Style

Robbiano L, Zhang Q. Logarithmic Decay of Wave Equation with Kelvin-Voigt Damping. Mathematics. 2020; 8(5):715. https://doi.org/10.3390/math8050715

Chicago/Turabian Style

Robbiano, Luc, and Qiong Zhang. 2020. "Logarithmic Decay of Wave Equation with Kelvin-Voigt Damping" Mathematics 8, no. 5: 715. https://doi.org/10.3390/math8050715

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