![Arxiv:1407.2172V1 [Math.AP] 8 Jul 2014 E Od N Phrases](http://data.docslib.org/img/c44e26c0e5dd10b340b40ad27809c387-1.webp)
arXiv:1407.2172v1 [math.AP] 8 Jul 2014 e od n phrases. and words Key 2010 ..Etnint o-isptv systems non-dissipative References to Extension equation beam 3.2. Euler-Bernoulli Damped result 3.1. applications Main Some of proof the 3. of End 2.3. basis Riesz of spectrum the 2.2. of Description result main the 2.1. of Proof 2. result Main Result Main 1.1. and Introduction 1. AEO EA FSM ERWK-IEDISSIPATIVE PETROWSKY-LIKE SOME OF DECAY OF RATE ahmtc ujc Classification. Subject Mathematics oeseta a odto.W ieas oeapiain oillustra the to applications of some part also give real We the condition. corresp of gap the if spectral supremum semigroup, some the underlying the by of generator given infinitesimal is systems dissipative Abstract. KA nti ae,w hwta h ats ea aefrsm Petr some for rate decay fastest the that show we paper, this In ¨ SAMR,MUZDMSI N AE ZERZERI MAHER AND DIMASSI, MOUEZ AMMARI, IS aeo ea,Ptosylk ytm,Seta bcsa Ries Abscissa, Spectral systems, Petrowsky-like Decay, of Rate 41 3Q23B034L20). 35B40 (35Q72 74K10 SYSTEMS Contents A B 1 nigoeao satisfied operator onding eorsetting. our te pcrmo the of spectrum owsky-like Basis. z 14 12 12 11 11 10 5 5 4 2 2 K.AMMARI,M.DIMASSI,ANDM.ZERZERI
1. Introduction and Main Result
The determination of optimal decay rate is difficult and has not a complete answer in the general case. In the 1-d case, it was performed mostly, see [1, 2, 5, 6, 7, 8, 9], and to references therein. For higher dimension, G. Lebeau gives in [13] the explicit (and optimal) value of the best decay rate in terms of the spectral abscissa of the generator of the semigroup and the mean value of a damping cofficient along the rays of geometrical optics. In this paper, we describe, in abstract setting, the optimal decay rate for some Petrowsky- like dissipative systems in terms of spectral quantity of the corresponding infinitesimal generator of the underlying dynamic. The main idea is to identify the optimal energy decay rate with the supremum of the real part of the associated dissipative operator. To do this, it is enough to show that the set of the corresponding generalized eigenvectors forms a Riesz basis of the energy space. The approach used here is based on some resolvent estimates which is obtained by a perturbative method. In case of the damped wave operator, Cox and Zuazua ([7]) adopt the shooting method based on an ansatz of Horn. This approach consists in constructing an explicit approxima- tion of the characteristic equation of the underlying system. Under the assumption that the damping is of bounded variation, they obtained high frequency asymptotic expansions of the spectrum. The shooting method can be used only for one-dimensional boundary value problems. In the both cases, we require precise knowledge of the spectrum of the corresponding non self-adjoint operators, more precisely, the behavior of the high frequency set. The advantage of our approach is that it works in any dimension and in a very general setting (see [15] and also [12]). Let us introduce the abstract setting. Let H be a Hilbert space equipped with the norm H . Let A be an unbounded operator on H, self-adjoint, positive and with compact k·kinverse. We denote its domain by (A). D Let B be a bounded operator from U to H, where U, U is another Hilbert space which will be identified with its dual. k·k We consider the following system:
x¨(t)+ Ax(t)+ BB∗x˙(t)=0, (1.1) x(0), x˙(0) =(x0, x1) H 1 H, ∈ 2 ×
1 where t [0, ) is the time and H 1 = (A 2 ) the scaled Hilbert space with the norm 2 ∈ 1 ∞ D z 1 = A 2 z H , z H 1 . k k 2 k k ∀ ∈ 2 From now on, we set := H 1 H. We endowe this space with the inner product: H 2 ×
1 1 [f,g] , [u, v] := A 2 f, A 2 u H + g, v H, for all [f,g], [u, v] in . H h i h i H D E RATE OF DECAY OF THE ENERGY 3
We can rewrite the system (1.1) as a first order differential equation, by putting Y (t)= T x(t), x˙(t) : Y˙ (t)+ Y (t)=0, B (1.2) Y (0) =AT (x , x ) , 0 1 ∈ H where B := : ( B)= ( ) , with A A0 − B D A D A0 ⊂H→H 0 I 0 = : ( 0)= (A) H 1 , A A 0 D A D × 2 ⊂H→H − 0 0 and = ( ). B 0 BB∗ ∈L H The operator is skew-adjoint on hence it generates a strongly continuous group A0 H of unitary operators on , denoted by S (t) . Since B is dissipative and onto, it H 0 t∈R A generates a contraction semi-group on , denoted by S (t) . The system (1.1) is B t∈R+ well-posed. More precisely, the followingH classical result holds. Proposition 1.1. Suppose that (x0, x1) . Then the problem (1.1) admits a unique ∈ H 1 solution t x(t) in the space C [0, + ); H 1 C [0, + ); H . Moreover the solution 7→ ∞ 2 ∩ ∞ t x(t) satisfies the following energy identity: 7→ t E x(0) E x(t) = B∗x˙(s) 2 ds, for all t 0, (1.3) − U ≥ Z0 1 2 where E x(t) = x(t), x˙(t) . 2 H
2 From (1.3) it follows that the mapping t x(t), x˙(t) is non-increasing. In many 7−→ H applications it is important to know if this mapping decays exponentially when t + , → ∞ i.e., if the system (1.1) is exponentially stable. One of the methods currently used for proving such exponential stability results is based on an observability inequality for the conservative system associated to the initial value problem φ¨(t)+ Aφ(t)=0, φ(0), φ˙(0) =(x , x ) . (1.4) 0 1 ∈ H It is well-known that (1.4) is well-posed in . The result below, proved in [11] (see also [3]), shows that the exponential stability of (1.1)H is equivalent to an observability inequality for (1.4). Proposition 1.2. The system described by (1.1) is exponentially stable in if and only H if there exist T > 0, and CT > 0 such that T 2 C B∗ t Y dt Y 2 , Y . T 0 S0( ) 0 0 H for all 0 (1.5) 0 U ≥ ∈ H Z The spectrum of A is given by 0 < µ1 µ2 µ3 µn + and the ≤ ≤ ≤ ···≤ ≤ ···→ ∞ family (vn)n≥1 of corresponding normalized eigenvectors of A is an orthonormal basis of H. Now, we can describe the spectrum of the skew-adjoint operator by the following: A0 4 K.AMMARI,M.DIMASSI,ANDM.ZERZERI
Lemma 1.3. The eigenvalues of and the corresponding eigenvectors are given by: A0 vk 1 ∗ 0V±k = i√µk V±k, where V±k = , i , for all k N . (1.6) A ± √2 √µk ± ∈ h i Moreover, the family V k ∗ is an orthonormal basis of the energy space . ± k∈N H 1.1. Main result. As mentionned in the above, we give the value of the fastest decay rate of solution of the equation (1.1), in terms of the spectral abscissa of the generator B. A Let µ(B) be the spectral abscissa of B given by: A µ( B) = sup Re(λ); λ σ( B) . (1.7) A ∈ A Here σ( B) denotes the spectrum of B. In order to state the result on the optimal decay rate, weA define the decay rate, dependingA on B, as ω(B) = inf ω; there exists C = C(ω) > 0 such that
E(x(t)) C(ω) e2ωtE(x(0)) for every solution of (1.1) with initial data in . (1.8) ≤ H According to (1.3), ω(B) 0 (see [11] and also [3]). It follows easily that, ≤ µ( B) ω(B). (1.9) A ≤
The following assumption concern the high frequencies of 0. More precisely, it deals with the behavior of the gap between two consecutive high frequenA cies of . For k N∗, A0 ∈ we define δ±k := i(√µk+1 √µk) = √µk+1 √µk. We assume that | ± − | − (A1) lim δk =+ , k→+∞ ∞ and δ (A2) k+1 l2(N∗), where l2(N∗) is the space of square integrable sequences. δ2 ∈ k k≥1 Remarks (i) The assumption (A1) implies that the high frequencies of are simple. A0 (ii) Assumption (A2) implies δ lim k+1 =0. (1.10) k→+∞ δ2 k (iii) Note that, in general, assumption (A2) does not imply hypothesis (A1). Now, our main result on the optimal decay rate is: Theorem 1.4. Assume (A1) and (A2). Then,
ω(B)= µ( B). (1.11) A RATE OF DECAY OF THE ENERGY 5
In other words if all finite energy solutions of (1.1) are exponentially stable then the fastest decay rate of the solution of (1.1) satisfies (1.11). Outline of the proof: In the following, we give an idea of the proof of the main result.
First, for B = 0, the operator 0 is skew-adjoint with compact resolvent in . From general operator theory, all its eigenvaluesA lie on the imaginary axis and the geometricH and algebraic multiplicity of each eigenvalue are the same. Moreover, there is a sequence of eigenvectors of which forms a Riesz (orthonormal, actually) basis for . A0 H In our setting, i.e., B (U, H), we give in Proposition 2.1 rough preliminary bounds ∈ L on the spectrum of B. Moreover, since B is a bounded perturbation of skew-adjoint operator with compactA resolvent it followsA from [10, Chapter 5, Theorem 10.1] that the generalized eigenvectors of B are complete in . For instance, these results are not enough A H to prove Theorem 1.4. We need to study the high frequency of B, and in particular their algebraic multiplicities. Using the fact that the distance betweeA n two consecutive eigenvalues tends to infinity at infinity, as well as the fact that the dissipation is bounded, (k) we construct in Subsection 2.1 a closed curves (Γ )|k|>N0 (for some integer N0 sufficiently large) in the complex plane such that: ∗ (±n) (i) For all n N , Γ is centered in ( i√µn). ∈ ± (n) (ii) Inside each Γ there exits exactly one simple eigenvalue of B. A (k) (iii) The operator B has exactly 2N0 eigenvalues including multiplicity in C ( Γ ). A \ |k|>N∪ 0 B 0 2 B 0 (iv) P (k) P (k) < , where P k (resp. P k ) denotes the Riesz projection k Γ − Γ kL(H) ∞ Γ( ) Γ( ) |k|>N0 X (k) associated to B (resp. ) corresponding to Γ . A A0 The proof of the above statements are based on some resolvent estimates of the operators B and . Since the generalized eigenvectors of B are complete and the systems of A A0 A (generalized)-eigenvectors of B and 0 are quadratically close in (see (iv) above), it follows from [14, Appendix D,A TheoremA 3] that the system of generaHlized eigenvectors of B constitutes a Riesz basis in . Now, by a standard argument, we identify the optimal A H energy decay rate with the supremum of the real part of B, which complete the proof of Theorem 1.4. A
2. Proof of the main result
As indicated in the introduction, we will establish Theorem 1.4 by proving that the system of generalized eigenvectors of the operator B constitutes a Riesz basis in the energy A space , and that all eigenvalues of B with sufficiently large modulus are algebraically simple.H A
2.1. Description of the spectrum of B. The operator B is a bounded perturbation A A of a skew-adjoint operator 0 then, according to [10, Chapter 5, Theorem 10.1], we have the following spectral result:A 6 K.AMMARI,M.DIMASSI,ANDM.ZERZERI
Proposition 2.1. The following properties hold:
i) The resolvent of B is compact. In particular the spectrum of B is discrete, i.e., A A B has a discrete eigenvalues of finite algebraic multiplicity. A ii) The spectrum of B is symmetric about the real axis and is contained in , where A C ∪ I = λ C; λ √µ , β Re(λ) 0 (2.1) C ∈ | | ≥ 1 − ≤ ≤ n 1 1 o = β β2 µ 2 , β2 µ 2 . (2.2) I − − − 1 + − 1 + 1 ∗ 2 Here β := B h< + , µ1 > 0, is the first eigenvaluei of A and (γ)+ = 2 k kL(H,U) ∞ max(γ, 0). iii) The root vectors of B are complete in . A H
Proof. We give only the proof of the second point of the proposition. Let λk := λk( B) A be an eigenvalue of B. We denote by W ( ; λk) the corresponding eigenvector. Then T A · W ( ; λk)= u( ,λk) (1,λk), where u( ; λk) satisfies · · · 2 ∗ λ u( ; λk)+ λkBB u( ; λk)+ Au( ; λk) = 0 with u( ,λk) H . (2.3) k · · · · ∈ 1 Since B is real it follows that W ( ; λk)= W ; λk is an eigenvector of B corresponding A · · A to the eigenvalue λk. We take the scalar product of the equation (2.3) with u( ,λk), we obtain: · 1 2 4 2 2 1 ∗ u( ; λk) 1 ∗ u( ; λk) 1 u( ; λk) λ±k = B · B · A 2 · . −2 u( ; λk) H U ± 4 u( ; λk) H U − u( ; λk) H H k · k k · k k · k Hence, if λ is a non-real eigenvalue, we find k
u λ 2 1 u λ 2 u λ 4 1 ∗ ( ; k) ( ; k) 1 ∗ ( ; k) λ±k = B · i A 2 · B · , −2 u( ; λk) H U ± s u( ; λk) H H − 4 u( ; λk) H U k · k k · k k · k ∗ which implies that, since B is bounded from H to U, 2 1 ∗ u( ; λk) 0 < β Re(λ±k)= B · 0, − ≤ −2 u( ; λk) H U ≤ k · k where β := 1 B∗ 2 < + , and 2 k kL(H,U) ∞ 2 2 1 u( ; λk) λ±k = A 2 · µ1. | | u( ; λk) H H ≥ k · k If λ is real we observe that k
4 2 1 1 u( ; λk) 1 u( ; λk) 2 B∗ A 2 β µ 2 . · · 1 + s4 u( ; λk) H U − u( ; λk) H H ≤ − k · k k · k Here µ1 > 0, is the first eigenvalue of A. RATE OF DECAY OF THE ENERGY 7
To state the principal result of this subsection (see Theorem 2.2), we need to introduce some notations. For n N∗, we define the three complex numbers: ∈ 1 1 1 an = √µn + δn , bn = δn + i√µn and dn = δn + i√µn, (2.4) −1 2 −1 2 −2 ∗ (n) where δk := √µk+1 √µk for k N . Let Int(Γ ) denote the rectangle with sides (n) (n) (n) (n−) ∈ γ1 ,γ2 ,γ3 and γ4 , (see Figure 1), where
(n) δn γ := λ C; Im(λ)= an and Re(λ) < , 1 ∈ | | 2 (n) δn γ := λ C; Re(λ)= and an Im(λ) an , 2 ∈ 2 ≤ ≤ +1
(n) δn δn γ := λ C; Im(λ)= an and Re(λ) goes from to , 3 ∈ +1 2 − 2 and (n) δn γ := λ C;Re(λ)= and Im(λ) goes from an to an . 4 ∈ − 2 +1 For n =1, 2, ...., we set Γ(n) = γ(n) γ(n) γ(n) γ(n), Γ(−n) := z C; z Γ(n) (2.5) 1 ∪ 2 ∪ 3 ∪ 4 { ∈ ∈ } and (n) δn−1 δn−1 C = z C; Im(z) < √µn + and Re(z) < . ∈ | | −1 2 | | 2 Note that by construction Int(Γ(k)) Int(Γ(n))= for all k, n Z∗ such that k = n. Here we denote the interior of Γ(k) by Int(Γ∩ (k)). Moreover,∅ for all∈N N∗ we have6 C(N) ( Int(Γ(k))), where and are given by (2.1) and (2.2).∈ C∪I ⊂ |k|≥∪N C I TheoremS 2.2. We assume (A1) and that (1.10) is satisfied. Then, there exists N N∗ 0 ∈ large enough such that the operator B has exactly 2N0 eigenvalues, including multiplicity, A (k) in CN and one simple eigenvalue in Int(Γ ) for each k with k > N . This exhausts the 0 | | 0 spectrum of B. A We have divided the proof into a sequence of lemmas. Lemma 2.3. Assume (A1). Then, there exists C > 0 and N N (large enough) such 0 ∈ that for n > N0, the following properties hold: (±n) (n) (i) Γ ∂C C σ( B) σ( ) . ∪ ⊂ \ A ∪ A0 (ii)
−1 −1 C (±n) (n) (λ B) (λ 0) L(H) 2 , uniformly on λ Γ ∂C , (2.6) k −A − −A k ≤ δn−1 ∈ ∪ where ∂C(n) is the boundary of the rectangle C(n). 8 K.AMMARI,M.DIMASSI,ANDM.ZERZERI
ipµn+1
ia (n) n+1 Γ
dn iõn bn
ian
ipµn 1 −
δn 0 δn − 2 2
Re(z) = B 2 −k ∗k (U,H) L
Figure 1. Location of σ( B) A
Proof. Since is skew-adjoint, it follows that A0 1 (λ )−1 . (2.7) k −A0 kL(H) ≤ dist(λ, σ( )) A0 By construction of Γ(±n) and C(n), we have:
(±n) dist(Γ , σ( )) = min bn i√µn , dn i√µn , ian i√µn , ian i√µn A0 | − | | − | | +1 − +1| | − −1| δ = n−1 , 2 and dist(∂C(n), σ( )) δn−1 , which together with (2.7) yields A0 ≥ 2 −1 2 (±n) (n) (λ 0) L(H) , uniformly on λ Γ ∂C . (2.8) k −A k ≤ δn−1 ∈ ∪
Recalling that = B is a bounded linear operator on defined by B A0 −A H 0 0 = . B 0 BB∗ From (2.8), we have 2 B∗ 2 −1 L(H,U) (±n) (n) (λ 0) L(H) k k uniformly on λ Γ ∂C . (2.9) kB −A k ≤ δn−1 ∈ ∪ RATE OF DECAY OF THE ENERGY 9
By (A1), we choose N such that for n N : 0 ≥ 0 ∗ 2 2 B H,U k kL( ) κ< 1 . δn−1 ≤ Now the first statement of the lemma follows from (2.8), (2.9) and the following obvious equality: −1 λ B = Id + (λ ) (λ ). (2.10) −A B −A0 −A0 On the other hand (2.10) yields h i
−1 −1 −1 −1 p (λ B) =(λ ) +(λ ) (λ ) , −A −A0 −A0 − B −A0 p≥1 X which together with (2.8) and (2.9) imply (2.6).
According to Lemma 2.3, for n N the following Riesz projections are well defined: ≥ 0 B 1 −1 0 1 −1 PΓ(±n) := (λ B) dλ, PΓ(±n) := (λ 0) dλ, (2.11) 2πi (±n) −A 2πi (±n) −A ZΓ ZΓ B 1 −1 0 1 −1 P∂C(n) := (λ B) dλ and P∂C(n) := (λ 0) dλ. 2πi (n) −A 2πi (n) −A Z∂C Z∂C The following result is a simple consequence of (A1), (1.10), (2.6) and the estimate of the measure on ∂C(n), Γ(±n).
Lemma 2.4. We assume (A1) and that (1.10) is satisfied. Then, there exists C > 0 (independent of n) and N N such that for n N , we have 0 ∈ ≥ 0 B 0 δn PΓ(±n) PΓ(±n) L(H) C 2 < 1, (2.12) k − k ≤ δn−1
B 0 δn P∂C(n) P∂C(n) L(H) C 2 < 1. (2.13) k − k ≤ δn−1 End of the proof of Theorem 2.2: First, recalling that if P and Q are two projectors with P Q < 1, then rank(P ) = rank(Q) (see Lemma 3.1 in [10]). Thus, in the notation of kLemma− k 3.3, we have B 0 B 0 rank(P (n) ) = rank(P (n) ), rank(P (±n) ) = rank(P (±n) ), for n N . ∂C ∂C Γ Γ ≥ 0 Next, we conclude from (2.1) and (2.2) that C(N0) Int(Γ(k)) , hence that C∪I⊂ |k|≥∪N0 (N0) (k) S σ( B) is a subset of C Int(Γ ) . Now Theorem 2.2 follows from the fact A |k|≥∪N0 that S 0 0 N ±n rank(P∂C( 0) )=2N0 and rank(PΓ( ) )=1. 10 K.AMMARI,M.DIMASSI,ANDM.ZERZERI
Remark 2.5. In the proofs of Lemmas 2.3-2.4, we have used only the fact that the distance between two consecutive eigenvalues of 0 tends to infinity at infinity and the fact that 0 is a skew-adjoint operator. Similar generalA result are well-known (see Theorem 4.15aA in [12]).
2.2. Riesz basis. We start this subsection by constructing the eigenvectors associated to ∗ the high frequencies of B. Since the high frequencies of B are simple then for all k N , A A ∈ k > N0 (N0 given by Theorem 2.2), we define B ϕ±k = P±kV±k, (2.14)
where V±k is the eigenvector of 0 associated to the eigenvalue i√µk given by (1.6), and B A 0 ∗ ± P ±k is given by (2.11). Note that P Vn = Vn for all n Z . Γ( ) n ∈ ∗ For n Z , we denote the eigenvalue of B by λn(B). We have the following proposition: ∈ A Proposition 2.6. We assume (A1) and that (1.10) is satisfied. For k Z∗, such that ∈ k > N , the function ϕk is an eigenvector of B associated to the eigenvalue λk := | | 0 A λk( B). Moreover, there exists C > 0 such that A δn ϕn Vn H C 2 , for all n such that n > N0. (2.15) k − k ≤ δn−1 | |
Here N0 is given by Theorem 2.2. ∗ In particular, ϕn =1+ o(1) uniformly for n Z , n > N . k kH ∈ | | 0 Z∗ B B Proof. For all m , m > N0, we have Bϕm = BPΓ(m) Vm = λm(B)PΓ(m) Vm = ∈ | | A 0 A λm(B)ϕm. Using Lemma 2.4 and the fact that P n Vn = Vn with Vn 2 = 1, we get: Γ( ) k kV ×L B 0 B 0 δm ϕm Vm H = (PΓ(m) PΓ(m) )Vm H PΓ(m) PΓ(m) L(H) C 2 , k − k k − k ≤k − k ≤ δm−1 ∗ for all m Z , m > N0,(C independent of m). In particular, parallelogram inequality and ∈ | | ∗ recalling that Vm = 1 give that ϕm =1+ o(1) uniformly for m Z , m > N . k kH k kH ∈ | | 0
Now, we complete the sequence (ϕk)|k|>N0 of the eigenvectors associated to the high frequencies of B by considering the generalized eigenvectors associated to the low fre- A quencies of B. Note that the number of these generalized eigenvectors associated to A ∗ the low frequencies of B is finite, at most 2N by Theorem 2.2. For k Z such that A 0 ∈ k N0, we denote by mk the algebraic multiplicity of λk := λk( B) and we associated | | ≤ mk−1 A W to it the Jordan chain of generalized eigenvectors, k,p p=0 , i.e., a Jordan basis of the mk root subspace k := W ; B λk W =0 , E ∈ H A − n o ′ BWk,l = λkWk,l, Wk,l, Wk,l′ =0, l Here pk is the dimension of the eigenspace Ek := W ; B λk W =0 , Ek k. ∈ H A − ⊂E Now, we take the family of generalized eigenvectorsn of B : o A B := Wk,p ϕn . |k|≤N0,0≤p≤mk−1 ∪ |n|>N0 Since Vect(B) = (see Proposition 2.1, (iii)) and by assumption 1.10 the family B is H quadratically close to the orthonormal basis Vk ∗ of eigenvectors of the operator k∈Z A0 (see (2.15)). Then it follows from the Fredholm Alternative, see e.g., [14, Appendix D, Theorem 3], the following result: Theorem 2.7. Assume (A1) and (A2). Then the set B is a Riesz basis for the energy space . Moreover, there exists a linear isomorphism Φ of such that for all n Z∗, H H ∈ n > N , ΦVn = ϕn and Φ Vect(Vn, n N ) = Vect(Wk,p, k N , 0 p mk 1). | | 0 | | ≤ 0 | | ≤ 0 ≤ ≤ − 2.3. End of the proof of Main result. Using Theorem 2.7, we may expand the initial data as mk−1 0 0 u , v = ck,pWk,p + cnϕn . |k|≤N p=0 |n|>N X 0 X X 0 Then the solution of (1.2) is given by mk p −1 tp−l u, ∂ u = exp(λ t) c W + c exp(λ t)ϕ . (2.18) t k k,p (p l)! k,l n n n |k|≤N0 p=0 l=0 − |n|>N0 X X X X Recalling from Theorem 2.2 that at most 2N0 eigenvalues may be of algebraic multi- plicity greater than one and that 2N0 is the maximum of such multiplicity, and the family V ±k k∈N∗ is an orthonormal basis of the energy space (see Lemma 1.3), then, by the linear isomorphism Φ, we get H 2 2 −1 2 2N0 E u(t) = u, ∂tu Φ Φ 1+ t exp 2µ(B)t E u(0) . H ≤k k k k Then ω(B) µ( B), this with inequality (1.9) we have established our main result. ≤ A Remark 2.8. Note that, we talk about under (resp. over) damping if 1 B∗ 2 is less 2 k kL(H,U) (resp. greater) than √µ1, see [7]. Recall that µ1 > 0, is the first eigenvalue of A. 3. Some applications Firstly, we give examples of dissipative systems which satisfy Assumption 1.10 and we deduce the main result for these samples. In the second part, we extend our result to some non-dissipative systems and we give an example that illustrates this situation. 12 K.AMMARI,M.DIMASSI,ANDM.ZERZERI 3.1. Damped Euler-Bernoulli beam equation. We consider the following system: 2 4 ∂t u(x, t)+ ∂xu(x, t)+2a(x)∂tu(x, t)=0, 0 2 2 2 1 U = L (0, 1), H = L (0, 1), H 1 = H (0, 1) H0 (0, 1), 2 ∩ d2u d2u (A)= u H4(0, 1) H1(0, 1); (0) = (1) = 0 , D ∈ ∩ 0 dx2 dx2 = [H2(0, 1) H1(0, 1)] L2(0, 1), H ∩ 0 × d4 A = , Bφ = B∗φ = 2a(x)φ, φ L2(0, 1). dx4 ∀ ∈ So, p 0 I 0 I 0 = d4 , B = d4 . A 4 0 A 4 2a(x) − dx − dx − The operator 0 is skew-adjoint and with compact inverse and the spectrum is given • A 2 2 ∗ by σ( 0)= ik π ,k N , then Assumptions (A1) and (A2) are satisfied. NoteA that the{± inequality∈ (1.5)} is satisfied according to [11], if a satisfies (3.4). So, • ω(B) < 0. As a direct implication of Theorem 1.4, we have the following result (this result • was proved in [4]): Proposition 3.1. The fastest decay rate is given by the spectral abscissa, i.e., ω(B)= µ( B). A 3.2. Extension to non-dissipative systems. We consider the system described by: x¨(t)+ Ax(t)+ Kx(t)=0, x(0), x˙(0) =(x0, x1) H 1 H, (3.6) ∈ 2 × where A is the same operator as above and K (H1 ,H). ∈L 2 We can rewrite the system (3.6) as a first order differential equation, by putting Y (t)= T x(t), x˙(t) : T Y˙ (t)+ KY (t)=0,Y (0) = (x , x ) , (3.7) A 0 1 ∈ H RATE OF DECAY OF THE ENERGY 13 where K := : ( )= ( ) , with A A0 −K D AK D A0 ⊂H→H 0 I 0 = : ( 0)= (A) H 1 , A A 0 D A D × 2 ⊂H→H − 0 0 and = ( ). K K 0 ∈L H − The system (3.6) is well-posed. More precisely, the following classical result holds. Proposition 3.2. Suppose that (x0, x1) . Then the problem (1.1) admits a unique ∈ H 1 solution x in the following space C [0, + ); H 1 C [0, + ); H . ∞ 2 ∩ ∞ We denote, 1 2 E x(t) = x(t), x˙(t) . 2 H Let µ(K) be the spectral abscissa of K given by: A µ(K) = sup Re(λ); λ σ( K) . (3.8) ∈ A Here σ( K) denotes the spectrum of K. A A We define the growth bound, depending on K, as ω(K) = inf ω; there exists C = C(ω) > 0 such that E(x(t)) C(ω) e2ωtE(x(0)) for every solution of (3.6) with initial data in . (3.9) ≤ H As above we can prove the following result. Theorem 3.3. Assume (A1) and (A2). Then, (i) The eigenvectors of the associated operator K form a Riesz basis in the energy space A (ii) H ω(K)= µ( K). (3.10) A Example : Euler-Bernoulli equation with force term We consider the following initial and boundary value problem: 2 4 2 ∂t u(x, t)+ ∂xu(x, t)+ p ∂xu(x, t)=0, 0 2 2 du H = L (0, 1), H 1 = u H (0, 1); u(0) = 0, (0) = 0 , 2 ∈ dx 14 K.AMMARI,M.DIMASSI,ANDM.ZERZERI and the operators are defined 0 Id 0 = d4 d2 , A 4 p 2 0 − dx − dx 2 3 4 d u d u ( 0)= (u, v) H (0, 1) H 1 H 1 , (1) = 0, (1) = 0 , D A ∈ ∩ 2 × 2 dx2 dx3 d2 and K = p 2 (H 1 ,H). dx ∈L 2 0 1 2 We have the for all (u ,u ) H 1 L (0, 1) the problem (3.11)-(3.13) admits a unique ∈ 2 × solution 1 2 u C([0, + ); H 1 ) C ([0, + ); L (0, 1)). ∈ ∞ 2 ∩ ∞ 2 2 The spectrum of 0 is given by ( ik π )k∈N∗ . Then Assumptions (A1) and (A2) are satisfied. Therefore,A according to Theorem± 3.3, we obtain Proposition 3.4. (i) The generalized eigenvectors of the associated operator K form a Riesz basis in the energy space . A H (ii) ω(K)= µ( K). A References [1] K. Ammari, A. Henrot, and M. Tucsnak, Optimal location of the actuator for the pointwise stabilization of a string, C. R. Acad. Sci. Paris S´er. I Math. 330 (2000), no. 4, 275–280. [2] K. Ammari, A. Henrot and M. Tucsnak, Asymptotic behaviour of the solutions and optimal location of the actuator for the pointwise stabilization of a string, Asymptot. Anal. 28 (2001), no. 3-4, 215–240. [3] K. Ammari and M. Tucsnak, Stabilization of second order evolution equations by a class of unbounded feedback, ESAIM Control Optim. Calc. Var. 6 (2001), 361–386. [4] K. Ammari, M. Dimassi and M. Zerzeri, The rate at which energy decays in a viscously damped hinged Euler-Bernoulli beam, accepted for publication in Journal of Differential Equations, Juin 2014. [5] A. Benaddi and B. Rao, Energy decay rate of wave equations with indefinite damping, J. Differential Equations 161 (2000), no. 2, 337–357. [6] C. Castro and S. Cox, Achieving arbitrarily large decay in the damped wave equation, SIAM J. Control Optim. 39 (2001), no. 6, 1748–1755. [7] S. Cox and E. Zuazua, The rate at which energy decays in a damped string., Comm. Partial Differential Equations 19 (1994), no. 1-2, 213–243. [8] , The rate at which energy decays in a string damped at one end, Indiana Univ. Math. J. 44 (1995), no. 2, 545–573. [9] P. Freitas, Optimizing the rate of decay of solutions of the wave equation using genetic algorithms: a counterexample to the constant damping conjecture, SIAM J. Control Optim. 37 (1999), no. 2, 376–387. [10] I.C. Gohberg and M.G. Kre˘ın, Introduction to the theory of linear nonselfadjoint operators, American Mathematical Society,, vol. 18, Providence, R.I., 1969. [11] A. Haraux, Une remarque sur la stabilisation de certains syst`emes du deuxi`eme ordre en temps, Portugal. Math. 46 (1989), no. 3, 245–258. [12] T. Kato, Perturbation theory for linear operators, Springer-Verlag, 1966. [13] G. Lebeau, Equation´ des ondes amorties. algebraic and geometric methods in mathematical physics (kaciveli, 1993), Math. Phys. Stud. 19 (1996), 73–109. RATE OF DECAY OF THE ENERGY 15 [14] J. P¨oschel and E. Trubowitz, Inverse spectral theory, Pure and Applied Mathematics, vol. 130, Aca- demic Press, Inc., Boston, MA, 1987. [15] A-G. Ramm, On the basis property for root vectors of some nonselfadjoint operators, Journal of Mathematical Analysis and Applications, 80 (1981), 57–66. Ka¨ıs Ammari, UR Analyse et Controleˆ des Edp, UR 13ES64, Dep.´ de Mathematiques,´ Faculte´ des Sciences de Monastir, Universite´ de Monastir, 5019 Monastir, Tunisie E-mail address: [email protected] Mouez Dimassi, Universite´ Bordeaux 1, CNRS, UMR 5251 IMB, 351, cours de la Libration 33405 Talence cedex, France E-mail address: [email protected] Maher Zerzeri, Universite´ Paris XIII, CNRS, UMR 7539 LAGA, 99, ave Jean-Baptiste Clement´ F-93430 Villetaneuse, France E-mail address: [email protected]