Abstract
In this paper we present new results on the preservation of polynomial stability of damped wave equations under addition of perturbing terms. We in particular introduce sufficient conditions for the stability of perturbed two-dimensional wave equations on rectangular domains, a one-dimensional weakly damped Webster’s equation, and a wave equation with an acoustic boundary condition. In the case of Webster’s equation, we use our results to compute explicit numerical bounds that guarantee the polynomial stability of the perturbed equation.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
In this paper we study the stability properties of damped wave equations and abstract second-order differential equations of the form [17, 29]
on a Hilbert space \(X_0\). Here \(L: \mathrm{dom\,}L \subset X_0\rightarrow X_0\) is a negative self-adjoint operator with a bounded inverse and \(D_0\in {{\mathfrak {L}}}(U,X_0)\) for some Hilbert space U. Our main interest is in the preservation of stability under bounded perturbations in the situation where the unperturbed differential Eq. (1.1) is only polynomially stable [7, 10] (as opposed to being uniformly exponentially stable). The polynomial stability of (1.1) means that there exist constants \(\alpha ,M>0\) such that for all initial conditions \(w_0\in \mathrm{dom\,}L \) and \(w_1\in \mathrm{dom\,}(-L)^{1/2}\) the solutions of (1.1) satisfy [10]
Polynomial stability has been investigated in detail in the literature for damped wave equations on multi-dimensional domains [5, 11, 20], coupled partial differential equations [6, 16, 27, 30], as well as abstract damped second-order systems of the form (1.1) [3, 4, 13, 15, 21].
Polynomial stability is a strictly weaker concept than exponential stability, and it can in particular be destroyed under addition of arbitrarily small lower order terms in the partial differential equation. In this paper we employ and refine the general framework introduced in [23, 24] to present conditions for preservation of the polynomial stability of the abstract differential Eq. (1.1) under finite-rank and Hilbert–Schmidt perturbations. Moreover, we study preservation of polynomial stability for selected partial differential equation models, namely, a damped two-dimensional wave equation on a rectangular domain, a weakly damped Webster’s equation, and a one-dimensional wave equation with a dynamic boundary condition.
As our first main results we present general conditions for the polynomial stability of perturbed second-order systems of the form
Here the operators \(B_2\in {{\mathfrak {L}}}(Y,X_0)\), \(C_1\in {{\mathfrak {L}}}(\mathrm{dom\,}(-L)^{1/2},Y)\), and \(C_2\in {{\mathfrak {L}}}(X_0,Y)\) for some Hilbert space Y describe the perturbations to the nominal polynomially stable Eq. (1.1). As our first main results we adapt and improve the main results in [23, 24] to make them more easily verifiable for second-order systems of the form (1.3). Our results show that if the unperturbed equation is polynomially stable so that (1.2) is satisfied with some \(\alpha \in (0,2]\), then (1.3) is polynomially stable provided that for exponents \(\beta ,\gamma \in [0,1]\) satisfying \(\beta +\gamma \ge \alpha \) the graph norms \(\Vert (-L)^{\beta /2}B_2\Vert \), \(\Vert (-L)^{(\gamma -1)/2}C_1^*\Vert \), and \(\Vert (-L)^{\gamma /2}C_2^*\Vert \) are finite and sufficiently small. Our new results also provide concrete bounds for the required sizes of these graph norms based on lower bounds for the operator \(D_0^*\) restricted to the spectral subspaces of L. The results are applicable in the situations where Y is either finite-dimensional or where \(B_2\), \(C_1\), and \(C_2\) are Hilbert–Schmidt operators.
As the first concrete partial differential equation we study a wave equation with viscous damping on a rectangle \(\Omega =(0,a)\times (0,b)\),
We assume the damping coefficient \(d(\cdot ,\cdot )\ge 0\) is strictly positive on some non-empty open subset of \(\Omega \) which does not satisfy the Geometric Control Condition (see, e.g., [5]). We apply our abstract results to present conditions for the polynomial stability of perturbed wave equations of the form
where \(b_{k,2},c_{k,1},c_{k,2}\in L^2(\Omega )\). In particular, our results show that the perturbed wave equation is polynomially stable provided that the coefficient functions \(b_{k,2}\), \(c_{k,1}\), and \(c_{k,2}\) have sufficient smoothness properties in the sense that these functions belong to fractional domains of \(-\Delta \) (the Dirichlet Laplacian on \(\Omega \)), and the associated fractional graph norms are sufficiently small. We present also analogous results of Hilbert–Schmidt perturbations of (1.4). Finally, we analyse the stability of (1.4) in a situation where the damping term is perturbed in a non-dissipative way with a rank one operator.
Our second concrete partial differential equation is a Webster’s equation with a weak damping on (0, 1),
where \(a>0\) and \(d(\cdot )\in L^2(0,1)\) is the damping coefficient. In this article we focus on a special case where \(d(x)=1-x\). We begin by proving that the Webster’s equation is polynomially stable with this particular damping coefficient. Using our abstract results we then present conditions for the preservation of the Webster’s equation under addition of a perturbation term. We also present a numerical example where we compute numerical bounds for the coefficient functions in the perturbation to guarantee the preservation of polynomial stability of the Webster’s equation.
Finally, as our third partial differential equation we consider a one-dimensional wave equation with a dynamic boundary condition. The polynomial stability of this model was shown in [1, 22], and in this paper we present conditions for the preservation of the stability under addition of perturbation terms to the differential equation.
We use the following notation. Given a closed operator A on a Hilbert space X, which will be assumed to be complex, we denote its domain by \(\mathrm{dom\,}A\), its kernel by \({\mathrm{ker\,}\,}A\), and its range by \(\mathrm{ran\,}A\). The spectrum of A is denoted by \(\sigma (A)\), and given \(\lambda \in \rho (A):={\mathbb {C}}{\setminus }\sigma (A)\) we write \(R(\lambda ,A)\) for the resolvent operator \((\lambda -A)^{-1}\). The space of bounded linear operators on X is denoted by \({{\mathfrak {L}}}(X)\). Given two functions \(f,g:(0,\infty )\rightarrow {\mathbb {R}}_+\), we write \(f(t)=O(g(t))\) to indicate that \(f(t)\le Cg(t)\) for some constant \(C>0\) and for all sufficiently large \(t>0\).
2 Robustness of Stability for Generalized Wave Equations
2.1 Polynomial Stability of Strongly Continuous Semigroups
The second-order differential Eq. (1.1) with a negative and boundedly invertible operator \(L: \mathrm{dom\,}L \subset X_0\rightarrow X_0\) and \(D_0\in {{\mathfrak {L}}}(U,X_0)\) can be represented as a first-order abstract Cauchy problem with state \(u(t)= ( w(t), w_t(t) )^\top \) as
with the initial condition \(u(0)=(w_0,w_1)^\top \). We choose the state space of this linear system as
The space \({{\mathcal {H}}}\) is a Hilbert space with inner product defined by
for all \(u=(u_1,u_2)^\top ,v=(v_1,v_2)^\top \in {{\mathcal {H}}}\). The domain of A is \(\mathrm{dom\,}A=\mathrm{dom\,}L\times \mathrm{dom\,}(-L)^{1/2}\). The operator A has the form \(A=A_0-DD^*\) where
Here \(A_0\) is a skew-adjoint operator and A generates a strongly continuous semigroup T(t) on \({{\mathcal {H}}}\) by the Lumer-Phillips theorem [17, Sec. VI.3].
Definition 2.1
( [10]) A strongly continuous semigroup T(t) generated by a linear operator A is said to be polynomially stable with \(\alpha >0\) if it is uniformly bounded, i.e. \(\sup _{t\ge 0}\Vert T(t)\Vert <\infty \), if \(i{\mathbb {R}}\subset \rho (A)\), and if
for some constant \(M>0\).
2.2 Polynomial Stability of Perturbed Semigroups
We are interested in robustness of the polynomial stability of (1.1) under perturbations of the form
where \(B_2\in {{{\mathfrak {L}}}}(Y,X_0)\), \(C_1:\mathrm{dom\,}(-L)^{1/2}\subset X_0\rightarrow Y\), and \(C_2\in {{{\mathfrak {L}}}}(X_0, Y)\) for some Hilbert space Y are such that \(C_1(-L)^{-1/2}\in {{{\mathfrak {L}}}}(X_0,Y)\). If we define \(B:=(0,B_2)^\top \in {{\mathfrak {L}}}(Y,{{\mathcal {H}}})\) and \(C:=(C_1,C_2)\in {{\mathfrak {L}}}({{\mathcal {H}}},Y)\), the perturbed system can be represented as an abstract Cauchy problem \(\frac{d u}{d t}=(A+BC)u\). The following theorem presented in [24] provides general conditions for the preservation of the polynomial stability of the semigroup \(T_{A+BC}(t)\) generated by \(A+BC\).
Theorem 2.2
([24, Thm. 6]) Assume T(t) generated by A is polynomially stable with \(\alpha >0\), let \(\beta ,\gamma \ge 0\) be such that \(\beta +\gamma \ge \alpha \), and let \(\kappa >0\) satisfy
If \(B\in {{\mathfrak {L}}}(Y,{{\mathcal {H}}})\) and \(C\in {{\mathfrak {L}}}({{\mathcal {H}}},Y)\) are such that
if \((-A)^\beta B\) and \((-A^*)^\gamma C^*\) are Hilbert–Schmidt operators, and if
then the semigroup generated by \(A+BC\) is polynomially stable with the same \(\alpha \).
The following theorem introduces a concrete bound \(\kappa >0\) for the norms of the perturbations in Theorem 2.2 for \(A=A_0-DD^*\) with a skew-adjoint operator \(A_0\) in the important special case \(\beta ,\gamma \ge 0\) are chosen so that \(\beta +\gamma =\lceil \alpha \rceil \) (here \(\lceil \alpha \rceil \in {\mathbb {N}}\) denotes the ceiling of \(\alpha >0\)). The first part of the result is a special case of [13, Thm. 3.5] with a proof which has been modified in a trivial manner to yield an explicit constant \(M_R>0\).
Theorem 2.3
Let X and U be Hilbert spaces, and assume \(A=A_0-DD^*\) where \(A_0: \mathrm{dom\,}A_0 \subset X\rightarrow X\) is skew-adjoint and \(D\in {{{\mathfrak {L}}}}(U,X)\). Let \(P_{(a,b)}\in {{{\mathfrak {L}}}}(X)\) be the spectral projection of \(A_0\) corresponding to the interval \((ia,ib)\subset i\mathbb {R}\). Assume there exist \(\eta _0,\delta _0>0\) and functions \(\eta :{\mathbb {R}}\rightarrow (0,\eta _0] \) and \(\delta :\mathbb {R}\rightarrow (0,\delta _0]\) such that
Then
where
If there exists \(M_0>0\) such that \(\eta (s)^{-2}\delta (s)^{-2}\le M_0(1+|s|^{\alpha }) \) for all \(s\in {\mathbb {R}}\), then for \(\beta ,\gamma \ge 0\) with \(\beta +\gamma =\lceil \alpha \rceil \) in Theorem 2.2 it is possible to choose any \(\kappa >0\) such that
where \(M_C>0\) is defined with an arbitrary \(s_0>0\) by
Proof
Assume that the functions \(\eta \) and \(\delta \) satisfy the assumptions of the theorem. Let \(y\in X\) and \(s\in {\mathbb {R}}\) be arbitrary and write \(x=R(is,A)y\in \mathrm{dom\,}A \). We then have \((is-A_0+DD^*)x=y\), and thus
Denote \(P_s:=P_{(s-\delta (s),s+\delta (s))}\) for brevity and write \(X=X_s\oplus _\perp X_\infty \) where \(X_s=P_sX\) and \(X_\infty = (I-P_s)X\). If we write \(x=x_0+x_\infty \) and \(y=y_0+y_\infty \) according to this decomposition, then
since the restriction \((is-A_0)\vert _{X_\infty }\) of \(is-A_0\) to \(X_\infty \) is boundedly invertible. Since \(A_0\) is skew-adjoint and \(\sigma ((is-A_0)\vert _{X_\infty })\subset i{\mathbb {R}}{\setminus } (-i\delta (s),i\delta (s))\), we have
By assumption we have \(\Vert x_0\Vert \le \eta (s)^{-1}\Vert D^*x_0\Vert \le \eta (s)^{-1}(\Vert D^*x\Vert +\Vert D^*x_\infty \Vert )\). If we denote \(q(s)= 1+ 2\eta (s)^{-2}\Vert D\Vert ^2\), we can use \(\Vert D^*x\Vert ^2\le \Vert x\Vert \Vert y\Vert \) and the Young’s inequality to estimate
This estimate implies
Recall that \(1\le \eta _0^2\eta (s)^{-2}\) and \(1\le \delta _0^2\delta (s)^{-2} \). We have
The estimate (2.6) implies
This completes the first part of the proof.
Assume now that there exists \(M_0>0\) such that \(\eta (s)^{-2}\delta (s)^{-2}\le M_0(1+|s|^\alpha ) \) for all \(s\in {\mathbb {R}}\) and denote \(n_\alpha = \lceil \alpha \rceil \in {\mathbb {N}}\). Then \(\Vert R(is,A)\Vert \le M_0M_R(1+|s|^\alpha )\) for all \(s\in {\mathbb {R}}\), and Theorem 2.2 implies that if \(\beta ,\gamma \ge 0\) are such that \(\beta +\gamma =n_\alpha \), then the constant \(\kappa >0\) is required to satisfy \(\kappa <(\sup _{{\lambda }\in \mathbb {C}_+}\Vert R({\lambda },A)A^{-n_\alpha }\Vert )^{-1/2}\). The approach in the proof of [8, Lem. 5.3] can be used to show that
This implies that \(\kappa >0\) in Theorem 2.2 can be chosen to have any value \(\kappa <1/ \sqrt{2M_C}\) provided that the constant \(M_C>0\) in the statement of the theorem is such that \(\Vert R(is,A)A^{-n_\alpha }\Vert \le M_C\) for all \(s\in {\mathbb {R}}\). In order to show this, let \(s_0>0\) be arbitrary and fixed. For any \(s\in {\mathbb {R}}\) with \(|s|\le s_0\) we have
On the other hand, if \(|s|\ge s_0\), then using the resolvent identity \(R(is,A)A^{-1} = (is)^{-1} (R(is,A)+A^{-1})\) repeatedly shows that
Combining the above two estimates shows that \(\sup _{s\in {\mathbb {R}}}\Vert R(is,A)A^{-\beta -\gamma }\Vert \le M_C\) for the constant \(M_C>0\) in the statement of the theorem, and thus the proof is complete. \(\square \)
Remark 2.4
In the case where \(-L\) has a complete set of orthonormal eigenvectors \(-L\phi _n=\mu _n\phi _n\) with \(0<\mu _1\le \mu _2\le \cdots \), the operator \(A_0\) in (2.1) has eigenvalues \(\lambda _n = \mathrm{sign\,}(n)i\sqrt{\mu _n}\) and a complete set of orthonormal eigenvectors \(\{\psi _n\}_{n\in \mathbb {Z}{\setminus } \{0\}}\) such that
In this situation for every \(s\in {\mathbb {R}}\) the spectral subspace \(\mathrm{ran\,}P_{(s-\delta (s),s+\delta (s))}\) of \(A_0\) consists of linear combinations of the eigenvectors \(\psi _n\) with every \(n\in \mathbb {Z}{\setminus } \{0\}\) for which \(s-\delta (s)<\mathrm{sign\,}(n)\sqrt{\mu _n}<s+\delta (s)\). The functions \(\eta :\mathbb {R} \rightarrow (0,\eta _0]\) and \(\delta :\mathbb {R}\rightarrow (0,\delta _0]\) in Theorem 2.3 should then be chosen so that
for all \(s\in {\mathbb {R}}\) and \(x=(x_1,x_2)^\top \in \mathrm{ran\,}P_{(s-\delta (s),s+\delta (s))}\). In particular, if \(\delta :{\mathbb {R}}\rightarrow (0,\delta _0]\) is chosen in such a way that \(\delta (-s)=\delta (s)\) and every interval \((i(s-\delta (s)),i(s+\delta (s)))\) contains at most one eigenvalue \({\lambda }_n = \mathrm{sign\,}(n)i\sqrt{\mu _n}\), then \(\eta :\mathbb {R}\rightarrow (0,\eta _0]\) in Theorem 2.3 can be chosen to be an even function satisfying
Remark 2.5
The second part of the proof of Theorem 2.3 can be extended in a straightforward manner to the more general case where \(\beta ,\gamma \ge 0\) are any exponents satisfying \(\beta +\gamma \ge \alpha \). Indeed, if we denote \(n_{\beta \gamma }=\lceil \beta +\gamma \rceil \), the moment inequality [17, Thm. II.5.34] with \(\theta =(n_{\beta \gamma }-\beta -\gamma )/n_{\beta \gamma }\) and a constant \(M_{\beta ,\gamma }>0\) can first be used estimate
and \(\Vert R(is,A)A^{-n_{\beta \gamma }}\Vert \) can be estimated using the resolvent identity similarly as before. However, in this case the constant \(M_C\) in the bound for \(\kappa >0\) has a more complicated formula.
2.3 Robustness Results for Wave Equations
The structure of the operator A allows us to improve the assumptions of Theorem 2.2 to overcome the difficulty of computing the graph norms of the fractional powers of the damped generators \(-(A_0-DD^*)\) and \(-(A_0-DD^*)^*\). Instead, the conditions are given in terms of the graph norms of the fractional powers of the positive operator \(-L\). Throughout this section \(C_1^*\) denotes the adjoint of \(C_1\) as an operator \(C_1:\mathrm{dom\,}C_1\subset X_0\rightarrow Y\).
Theorem 2.6
Assume that the strongly continuous semigroup T(t) generated by
is polynomially stable with \(\alpha \le 2\), that \(0\le \beta ,\gamma \le 1\) are such that \(\beta +\gamma \ge \alpha \), and that \(\kappa >0\) is as in Theorem 2.2. If the perturbation operators \(B=(0,B_2)^\top \in {{\mathfrak {L}}}(Y,{{\mathcal {H}}})\) and \(C=(C_1,C_2)\in {{\mathfrak {L}}}({{\mathcal {H}}},Y)\) satisfy
if \((-L)^{\beta /2} B_2\), \((-L)^\frac{\gamma -1}{2} C^*_1\), and \((-L)^{\gamma /2} C^*_2\) are Hilbert–Schmidt operators, and if
then the semigroup generated by \(A+BC\) is polynomially stable with the same \(\alpha \). Here \(K_\theta = e^{\frac{1}{2}{\pi ^2\theta (1-\theta )}} M^\theta \) with \(\theta \in [0,1]\) and \(M=1+\Vert D_0\Vert ^2\Vert (-L)^{-1/2}\Vert \).
For proving this result we use the following theorem from [19].
Theorem 2.7
( [19, Thm. 1]) If \({\mathcal {A}}_1\), \({\mathcal {A}}_2\) are closed maximal accretive operators on a Hilbert space \({{\mathfrak {H}}}\) such that \(\mathrm{dom\,}{\mathcal {A}}_1 \subset \mathrm{dom\,}{\mathcal {A}}_2\) and \(\Vert {\mathcal {A}}_2 u\Vert \le M\Vert {\mathcal {A}}_1 u\Vert \) for some constant \(M>0\) and for all \(u\in \mathrm{dom\,}{\mathcal {A}}_1\), then \(\mathrm{dom\,}{\mathcal {A}}_1^\theta \subset \mathrm{dom\,}{\mathcal {A}}_2^\theta \) and
where \(K_\theta = e^{\frac{1}{2}{\pi ^2\theta (1-\theta )}} M^\theta \).
Proof
(Proof of Theorem 2.6) Let \(0\le \beta ,\gamma \le 1\) be such that \(\beta +\gamma \ge \alpha \). Our aim is to show that if \(B_2\), \(C_1\), and \(C_2\) satisfy the given assumptions, then \(B=(0,B_2)^\top \) and \(C=(C_1,C_2)\) satisfy (2.4) with the same \(\kappa >0\). The stability of the semigroup generated by \(A+BC\) then follows directly from Theorem 2.2. To this end let \(\kappa >0\) be as in Theorem 2.2 and suppose that (2.7) hold. Define \(A_d: \mathrm{dom\,}A_d\subset {{\mathcal {H}}}\rightarrow {{\mathcal {H}}}\) and \(A_0: \mathrm{dom\,}A_0\subset {{\mathcal {H}}}\rightarrow {{\mathcal {H}}}\) with domains \(\mathrm{dom\,}A_d = \mathrm{dom\,}A_0 = \mathrm{dom\,}A = \mathrm{dom\,}L\times \mathrm{dom\,}(-L)^{1/2}\) by
and let \(D=(0,D_0)^\top \in {{{\mathfrak {L}}}}(U,{{\mathcal {H}}})\). Clearly
for all \(0\le \theta \le 1\). Since \(C\in {{{\mathfrak {L}}}}({{\mathcal {H}}},Y)\) and \(C^*= (-L^{-1} C_1^*,C_2^*)^\top \), the assumptions on \(B_2\), \(C_1\), and \(C_2\) imply that \(\mathrm{ran\,}B\subset \mathrm{dom\,}A_d^\beta \) and \(\mathrm{ran\,}C^*\subset \mathrm{dom\,}A_d^\gamma \). For every \(u=(u_1,u_2)^\top \in \mathrm{dom\,}A_d \) we have \(\Vert A_0u\Vert _{{\mathcal {H}}}^2 = \Vert (-L)^{1/2}u_2\Vert ^2+\Vert Lu_1\Vert ^2 = \Vert A_d u\Vert _{{\mathcal {H}}}^2\), and
where \(\Vert DD^*A_0^{-1}\Vert \le \Vert D_0\Vert ^2 \Vert (-L)^{-1/2}\Vert \). An analogous argument shows that we have \(\Vert A^*u\Vert _{{\mathcal {H}}}= \Vert (A_0+DD^*) u\Vert \le (1+\Vert D_0\Vert ^2 \Vert (-L)^{-1/2}\Vert )\Vert A_d u\Vert _{{\mathcal {H}}}\) for all \(u\in \mathrm{dom\,}A_d\). Since \(-A\), \(-A^*\) and \(A_d\) are closed and maximally accretive operators and \(\mathrm{dom\,}A = \mathrm{dom\,}A^*= \mathrm{dom\,}A_d\), Theorem 2.7 implies that \(\mathrm{ran\,}B\subset \mathrm{dom\,}(-A)^\beta \) and \(\mathrm{ran\,}C^*\subset \mathrm{dom\,}(-A^*)^\gamma \), and for all \(y\in Y\) with \(\Vert y\Vert =1\) we have
By Theorem 2.2 the semigroup generated by \(A+BC\) is polynomially stable with \(\alpha \). \(\square \)
If the operator L is diagonalizable [29, Sec. 2.6], then for \(\theta \in \mathbb {R}\) the spaces \(\mathrm{dom\,}(-L)^\theta \) and the graph norms of \((-L)^\theta \) have the forms
where \(\mu _{k}\) are the eigenvalues of \(-L\) and \(\phi _{k}\) are the corresponding orthonormal eigenvectors. With these definitions the space \(H_{-\theta }(L)\) is the dual of \(H_\theta (L)\) with respect to the pivot space \(X_0\) [29, Sec. 2.9].
Corollary 2.8
Assume that L is diagonalisable, that the strongly continuous semigroup T(t) generated by A is polynomially stable with \(\alpha \le 2\), \(0\le \beta ,\gamma \le 1\) satisfy \(\beta +\gamma \ge \alpha \), and \(\kappa >0\) is as in Theorem 2.2. If the perturbation operators \(B=(0,B_2)^\top \in {{\mathfrak {L}}}(Y,{{\mathcal {H}}})\) and \(C=(C_1,C_2)\in {{\mathfrak {L}}}({{\mathcal {H}}},Y)\) satisfy
if \((-L)^{\beta /2} B_2\), \((-L)^\frac{\gamma -1}{2} C^*_1\), and \((-L)^{\gamma /2} C^*_2\) are Hilbert–Schmidt operators and if
then the semigroup generated by \(A+BC\) is polynomially stable with the same \(\alpha \). Here \(K_\theta = e^{\frac{1}{2}{\pi ^2\theta (1-\theta )}} M^\theta \) with \(\theta \in [0,1]\) and \(M=1+\Vert D_0\Vert ^2\Vert (-L)^{-1/2}\Vert \).
3 Perturbations of Damped Two-Dimensional Wave Equations
In this section we consider damped wave equations on rectangular domains with different damping functions. We use Theorem 2.6 to derive concrete conditions for preservation the polynomial stability of perturbed wave equations with finite rank and Hilbert-Schmidt perturbations. We consider the damped wave Eq. (1.4) on \(\Omega = (0,a)\times (0,b)\), \(a,b>0\), with a damping coefficient \(d(\cdot ,\cdot )\in L^\infty (\Omega )\). The equation is of the form (1.1) on \(X_0=L^2(\Omega )\) with the choice \(L=\Delta \) and domain \(\mathrm{dom\,}\Delta =H^2(\Omega )\cap H_0^1(\Omega )\), and with \(D_0\in {{{\mathfrak {L}}}}(L^2(\Omega ))\) defined as the multiplication operator such that \(D_0u=\sqrt{d(\cdot ,\cdot )}u(\cdot ,\cdot )\) for all \(u\in L^2(\Omega )\).
We suppose that the set \(\omega =\{d(x,y)>0\}\) contains an open, nonempty subset and does not satisfy Geometric Control Condition (GCC) (see a definition of GCC for example in [5, Sec. 1]). It was shown in [18] that for such damping the Schrödinger group is observable, i.e., the pair \((D_0^*, i(-\Delta ))\) is exactly observable [29, Def. 6.1.1] (see also [12]). In this case the damped wave equation (1.4) is polynomially stable with \(\alpha = 2\) by [5, Thm. 2.3].
Our assumptions together with the results in [18] and [13, Prop. 3.9] also imply that the condition (2.5) is satisfied for some functions \(\eta :{\mathbb {R}}\rightarrow (0,\eta _0] \) and \(\delta :\mathbb {R}\rightarrow (0,\delta _0]\) satisfying \(\eta (s)^{-2}\delta (s)^{-2}\le M_0(1+s^2) \) for all \(s\in {\mathbb {R}}\). Because of this, Theorem 2.3 could in principle be used to derive numerical values for \(\kappa >0\) for particular damping functions \(d(\cdot ,\cdot )\). In practice, however, finding suitable concrete functions \(\eta \) and \(\delta \) can be challenging, and in the case of the two-dimensional wave equation this is an important topic for further research.
Remark 3.1
In some cases of damping functions the estimate for the exponent of polynomial stability can be improved. For example, in [28] the exponent of polynomial stability for the damping function
was shown to be \(\alpha =3/2\). Moreover, additional differentiability assumptions on \(d(\cdot ,\cdot )\) improve the rate of polynomial decay, as shown in [5, 11, 14].
3.1 Rank One Perturbations
We begin by considering perturbed wave equations of the form
with \(b_{2},c_{2}\in L^2(\Omega )\) and \(c_{1}\in H_{-1/2}(\Delta )\). The following theorem presents sufficient conditions for the polynomial stability of (3.2).
Theorem 3.2
Assume that damped wave equation (1.4) is polynomially stable with \(\alpha \le 2\), \(0\le \beta ,\gamma \le 1\) are such that \(\beta +\gamma \ge \alpha \), and \(\kappa >0\) is as in Theorem 2.6. If \(b_2\in H_{\beta /2}(\Delta )\), \(c_1\in H_{{(\gamma -1)/2}}(\Delta )\), \(c_2\in H_{\gamma /2}(\Delta )\) satisfy
then perturbed wave equation (3.2) is polynomially stable with the same \(\alpha \). Here \(K_\theta = e^{\frac{1}{2}{\pi ^2\theta (1-\theta )}} M^\theta \) with \(M = 1+\frac{ab \Vert d\Vert _{L^\infty }}{\pi \sqrt{a^2+b^2}}\) and \(\theta \in [0,1]\). For such perturbations there exists \(M_T>0\) such that the solutions of (3.2) corresponding to initial conditions \(w_0\in H^2(\Omega )\cap H_0^1(\Omega )\) and \(w_1\in H_0^1(\Omega )\) satisfy
Proof
In this case the perturbed wave equation has the form (2.2) with \(Y=\mathbb {C}\) and
where \(\langle \cdot ,\cdot \rangle _{H_{1/2},H_{-1/2}}\) denotes the dual pairing between \(H_{1/2}(\Delta )\) and \(H_{-1/2}(\Delta )\). Since \(B_2=b_2\), \(C_1^*= c_1\) and \(C_2^*= c_2\), the claim follows from Corollary 2.8 and a suitable upper bound for \(\Vert D_0\Vert ^2\left\| (-\Delta )^{-1/2}\right\| \). Since \(-\Delta \) is a positive self-adjoint operator with compact resolvent and its smallest eigenvalue is \(\pi ^2 (a^2+b^2)/(a^2b^2)\), we have \(\left\| (-\Delta )^{-1/2}\right\| \le \frac{ab}{\pi \sqrt{a^2+b^2}}\). Since \(\Vert D_0\Vert = \Vert \sqrt{d(\cdot ,\cdot )}\Vert _{L^\infty }=\sqrt{\Vert d(\cdot ,\cdot )\Vert _{L^\infty }}\), the claim holds for \(M = 1+\frac{ ab\Vert d\Vert _{L^\infty }}{\pi \sqrt{a^2+b^2}}\). \(\square \)
Remark 3.3
Conditions (3.3) have the simplest form if we choose \(\beta =\gamma =1\)
where \(M = 1+\frac{ ab\Vert d\Vert _{L^\infty }}{\pi \sqrt{a^2+b^2}}\).
Remark 3.4
If the damping function is as in (3.1), then the exponent of polynomial stability is \(\alpha =3/2\). In this case \(\beta \), \(\gamma \) can be chosen as \(\beta =\frac{1}{2}\) and \(\gamma =1\) and conditions (3.3) take the form
where \(K_{1/2}=e^{\frac{\pi ^2}{8}}\sqrt{M}\) with \(M = 1+\frac{ ab}{\pi \sqrt{a^2+b^2}}\).
3.2 Finite Rank Perturbations
We now consider the wave Eq. (1.4) with a finite number of perturbation terms
with \(b_{k,2},c_{k,2}\in L^2(\Omega )\) and \(c_{k,1}\in H_{-1/2}(\Delta )\).
Theorem 3.5
Assume that damped wave equation (1.4) is polynomially stable with \(\alpha \le 2\), \(0\le \beta ,\gamma \le 1\) are such that \(\beta +\gamma \ge \alpha \), and \(\kappa >0\) is as in Theorem 2.6. If for all \(k\in \{1,\ldots ,m\}\) we have \(b_{k,2}\in H_{\beta /2}(\Delta )\), \(c_{k,1}\in H_{{(\gamma -1)/2}}(\Delta )\), and \(c_{k,2}\in H_{\gamma /2}(\Delta )\) and
then the perturbed wave equation (3.4) is polynomially stable with the same \(\alpha \). Here \(K_\theta = e^{\frac{1}{2}{\pi ^2\theta (1-\theta )}} M^\theta \) with \(M = 1+\frac{ ab\Vert d\Vert _{L^\infty }}{\pi \sqrt{a^2+b^2}}\) and \(\theta \in [0,1]\). For such perturbations there exists \(M_T>0\) such that the solutions of (3.4) corresponding to initial conditions \(w_0\in H^2(\Omega )\cap H_0^1(\Omega )\) and \(w_1\in H_0^1(\Omega )\) satisfy
Proof
The perturbation can be written in the form (2.2) with the choice \(Y=\mathbb {C}^m\) and defining \(B_2\in {{{\mathfrak {L}}}}(Y,X_0)\), \(C_1\in {{{\mathfrak {L}}}}(\mathrm{dom\,}(-L)^{1/2},Y)\), and \(C_2\in {{{\mathfrak {L}}}}(X_0,Y)\) so that
for all \( y=(y_k)_{k=1}^m \in Y\), \(x\in H_{1/2}(\Delta )\) and \(z\in L^2(\Omega )\). For any \(0\le \beta ,\gamma \le 1\) we have
and thus the claims follow from Corollary 2.8 as in the proof of Theorem 3.2. \(\square \)
3.3 Hilbert–Schmidt Perturbations
Now we consider a more general case of perturbations of the wave equation
where the functions \(b_{k,2},c_{k,2}\in L^2(\Omega )\) and \(c_{k,1}\in H_{-1/2}(\Delta )\) of the perturbation are assumed to satisfy
The stability of this perturbed wave equation can be studied using Corollary 2.8 for Hilbert–Schmidt perturbations.
Theorem 3.6
Assume that damped wave equation (1.4) is polynomially stable with \(\alpha \le 2\), \(0\le \beta ,\gamma \le 1\) are such that \(\beta +\gamma \ge \alpha \), and \(\kappa >0\) is as in Theorem 2.6. If for all \(k\in \mathbb {N}\) we have \(b_{k,2}\in H_{\beta /2}(\Delta )\), \(c_{k,1}\in H_{{(\gamma -1)/2}}(\Delta )\), and \(c_{k,2}\in H_{\gamma /2}(\Delta )\) and
then the perturbed wave equation (3.6) is polynomially stable with the same \(\alpha \). Here \(K_\theta = e^{\frac{1}{2}{\pi ^2\theta (1-\theta )}} M^\theta \) with \(M = 1+\frac{ ab\Vert d\Vert _{L^\infty }}{\pi \sqrt{a^2+b^2}}\) and \(\theta \in [0,1]\). For such perturbations there exists \(M_T>0\) such that the solutions of (3.6) corresponding to initial conditions \(w_0\in H^2(\Omega )\cap H_0^1(\Omega )\) and \(w_1\in H_0^1(\Omega )\) satisfy
Proof
The perturbations can be written in the form (2.2) with \(Y=\ell ^2({{\mathbb {C}}})\) if we define \(B_2\), \(C_1\), and \(C_2\) so that
The assumptions \((\Vert b_{k,2}\Vert _{L^2})_k\in \ell ^2\), \((\Vert c_{k,1}\Vert _{H_{-1/2}})_k\in \ell ^2\), and \((\Vert c_{k,2}\Vert _{L^2})_k\in \ell ^2\) imply that \(B_2\in {{{\mathfrak {L}}}}(Y,X_0)\), \(C_1\in {{{\mathfrak {L}}}}(H_{1/2},Y)\) and \(C_2\in {{{\mathfrak {L}}}}(X_0,Y)\). If we let \(0\le \beta ,\gamma \le 1\) be such that \(\beta +\gamma \ge \alpha \), then
imply that \((-\Delta )^{\beta /2} B_2\), \((-\Delta )^\frac{\gamma -1}{2} C^*_1\), and \((-\Delta )^{\gamma /2} C^*_2\) are Hilbert–Schmidt operators and that
Thus the claim follows from Corollary 2.8 as in the proof of Theorem 3.2. \(\square \)
3.4 Wave Equation with "Almost Dissipative" Damping
Finally, we consider the two-dimensional damped wave equation with a perturbed damping term, namely
with \(b_2,c\in L^2(\Omega )\). We also make an additional assumption that \(d\in C^2(\Omega )\). The structure of the perturbed semigroup generator is now \({{\widetilde{A}} }:=A_0-(D+B)D^*=A-BD^*\), where \(D=(0,\sqrt{d(\cdot )})^\top \) and \(B=(0,b_2\langle \cdot ,c\rangle _{L^2})^\top \). Because of this structure, the damping in the wave Eq. (3.7) can be thought to be “almost dissipative”.
Since we assumed that the damping coefficient is smooth, i.e. \(d(\cdot )\in C^2(\Omega )\), it is possible to characterise the higher order domain \(\mathrm{dom\,}A^2\) and the stability of (3.7) can be studied using Theorem 2.2 with the parameters \(\beta =2\) and \(\gamma =0\), as shown in the following theorem.
Theorem 3.7
Assume that damped wave equation (3.7) is polynomially stable with \(\alpha \le 2\) in the case where \(b_2=0\). There exists \(\kappa >0\) such that if \(b_2\in \mathrm{dom\,}\Delta \) and \(c\in L^2(\Omega )\) satisfy
where \(M=1+\frac{ab\Vert d\Vert ^2_{L^\infty }}{\pi \sqrt{a^2+b^2}}\), then (3.7) is polynomially stable with the same \(\alpha \). For such perturbations there exists \(M_T>0\) such that the solutions of (3.7) corresponding to initial conditions \(w_0\in H^2(\Omega )\cap H_0^1(\Omega )\) and \(w_1\in H_0^1(\Omega )\) satisfy
Proof
Let \(\kappa >0\) be as in Theorem 2.2 and suppose the assumptions on \(b_2\) and c are satisfied. We define \({{\widetilde{B}} }:=(0,b_2)^\top \) and \({{\widetilde{C}} }:=(0,\langle \cdot ,\sqrt{d}c\rangle )\). It is clear that \({{\widetilde{B}} }{{\widetilde{C}} }=BD^*\). Our aim is to verify that the conditions of Theorem 2.2 are satisfied for the perturbed operator for \(A-{{\widetilde{B}} }{{\widetilde{C}} }\) with parameters \(\beta =2\) and \(\gamma =0\). We have
and thus \(\mathrm{ran\,}B= \{0\}\times \mathrm{span}\{b_j\} \subset \mathrm{dom\,}A^2\) provided that \(b_2\in \mathrm{dom\,}\Delta \) and \(db_2\in \mathrm{dom\,}(-\Delta )^{1/2}\). Since \(d\in C^2(\Omega )\), the assumption \(b_2\in \mathrm{dom\,}\Delta \) also implies \(db_2\in \mathrm{dom\,}(-\Delta )^{1/2}\).
The norm of \(A^2B\) can be estimated by
where the last estimate is completed as in the proof of Theorem 3.2. Moreover,
Thus \(\Vert (-A)^2B\Vert <\kappa \). We also have that \(\Vert (-A^*)^0C^*\Vert =\Vert \sqrt{d}c\Vert _{L^2}<\kappa \). The polynomial stability of the semigroup generated by \(A-{{\widetilde{B}} }{{\widetilde{C}} }=A- BD^*\) follows from Theorem 2.2 and then wave Eq. (3.7) is polynomial stable with \(\alpha \). \(\square \)
4 Perturbations of Webster’s Equations
In this section we show the polynomial stability of weakly damped Webster’s equation and use Theorems 2.6 and 2.3 to derive sufficient conditions for the preservation of the stability under addition of perturbing terms. We begin by considering an undamped Webster’s equation on \(\Omega =(0,1)\) which has the form
We consider \(r(x)=e^{ax}\), where \(a\ge 0\). Then Webster’s equation takes the form
We denote by \(L_a^2(0,1)\) the Hilbert space \(L^2(0,1)\) with the inner product
Let us define the operator \(L=\frac{d^2 }{d x^2}+a\frac{d }{d x}\) from \(L^2_a(0,1)\) to \(L^2_a(0,1)\) with \(\mathrm{dom\,}L=\{h\in L^2_a(0,1):h,h' \text { are absolutely continuous }, h''\in L^2_a(0,1) \text { and } h(0)=h(1)=0\}\). In the next lemma we state some properties of L.
Lemma 4.1
The operator \(L=\frac{d^2 }{d x^2}+a\frac{d }{d x}\) from \(L^2_a(0,1)\) to \(L^2_a(0,1)\) is a negative self-adjoint operator with a bounded inverse. The eigenvalues and eigenvectors of L are
respectively, for \(n\in {\mathbb {N}}\).
Proof
We define a unitary mapping \(V:L^2(0,1)\rightarrow L^2_a(0,1)\) by the formula
Now we can consider an auxiliary operator \({{\widetilde{L}} }:L^2(0,1)\rightarrow L^2(0,1)\) defined by \({{\widetilde{L}} }= V^*LV\) with \(\mathrm{dom\,}{{\widetilde{L}} }=H^2(0,1)\cap H^1_0(0,1).\) Direct calculations yield that \({{\widetilde{L}} }f=\frac{d^2 f }{d x^2}-\frac{a^2}{4}f\). It is well known that \({{\widetilde{L}} }\) is a negative self-adjoint operator with a bounded inverse. Hence L is also a negative self-adjoint operator with a bounded inverse. The eigenvalues and the eigenvectors of the operator \({{\widetilde{L}} }\) are
Since \({{\widetilde{L}} }= V^*LV\), the operators \({{\widetilde{L}} }\) and L have the same eigenvalues and the eigenvectors of L are given by the formula \(\varphi _n(x)=(V{{\widetilde{\varphi }} }_n)(x)=e^{-\frac{ax}{2}}\sin (\pi nx)\). \(\square \)
Now we consider weakly damped Webster’s equation
where the damping coefficient is \(d\in L^2_a(0,1)\). This equation is of the form (1.1) on \(X_0=L_a^2(0,1)\) with L defined above and with a rank one operator \(D_0 = d(\cdot )\in {{{\mathfrak {L}}}}({{\mathbb {C}}},L_a^2(0,1))\) and \(D_0^*=\langle \cdot ,d\rangle _{L^2_a}\).
The polynomial stability of the weakly damped Webster’s equation can be analyzed using [25, Thm. 6.3]. The following result in particular shows that (4.1) is polynomially stable for the particular choice of damping \(d(x)=1-x\).
Proposition 4.2
The weakly damped Webster’s equation (4.1) with the damping function \(d(x)=1-x\) is polynomially stable with \(\alpha =2\).
Proof
We can write
with \(\mathrm{dom\,}A_0 = \mathrm{dom\,}A\) and \(D\in {{{\mathfrak {L}}}}({{\mathbb {C}}},X)\). We will use [25, Thm. 6.3] to show that \(\Vert R(is,A)\Vert \le M(1+s^2)\) for some \(M>0\). To this end, we need to estimate the quantities \(|D^*\psi _n|\) from below, where \(\psi _n\) are the normalized eigenvectors of \(A_0\). Since L has eigenvalues \(\mu _n=-\frac{a^2}{4}-\pi ^2n^2\) with the corresponding eigenvectors \(\varphi _n(x)=e^{-\frac{ax}{2}}\sin (\pi nx)\) for \(n\in {\mathbb {N}}\), the eigenvectors \(\psi _n\) and the corresponding eigenvalues \(\lambda _n\) of \(A_0\) are given by
For any \(n\in \mathbb {Z}{\setminus } \{0\}\) we thus have
for some constant \(c>0\) and for all sufficiently large |n|. By [25, Thm. 6.3] we have \(\Vert R(is,A)\Vert =O( s^2)\) for |s| large, and thus [10, Thm. 2.4] implies that the semigroup generated by \(A=A_0-DD^*\) is polynomially stable with \(\alpha =2\). \(\square \)
Remark 4.3
Note that if in this weakly damped Webster’s equation one takes \(a=0\) then we get a weakly damped wave equation on the interval (0, 1) with the same damping coefficient \(d(x)=1-x\) and such equation is also polynomially stable with \(\alpha = 2\).
We consider the weakly damped Webster’s equation with additional perturbing terms of the form
where \(b_2,c_2\in L^2(0,1)\) and \(c_1\in H_{-1/2}(L)\). The following theorem presents conditions for the polynomial stability of the perturbed Webster’s Eq. (4.2). The spaces \(H_{\theta }(L)\) and the corresponding norms are defined as in (2.8). The above perturbations correspond to rank one perturbation operators in the abstract wave equation. Addition of multiple perturbation terms can be treated similarly as in the case of the two-dimensional wave equation in Sect. 3.
Theorem 4.4
Assume that the weakly damped Webster’s equation (4.1) is polynomially stable with \(\alpha \le 2\), that \(0\le \beta ,\gamma \le 1\) such that \(\beta +\gamma \ge \alpha \), and that \(\kappa >0\) is as in Theorem 2.2. If \( b_2\in H_{\beta /2}(L)\), \(c_1\in H_{{(\gamma -1)/2}}(L)\), \(c_2\in H_{\gamma /2}(L)\) satisfy
then the perturbed Webster’s Eq. (4.2) is polynomially stable with the same \(\alpha \). Here \(K_\theta = e^{\frac{1}{2}{\pi ^2\theta (1-\theta )}} M^\theta \), \(\theta \in [0,1]\), and \(M=1+ \Vert d\Vert _{L_a^2}(\frac{a^2}{4}+\pi ^2)^{-1/2}.\)
Proof
The perturbed system is of the form (1.3) with
where \(\langle \cdot ,\cdot \rangle _{H_{1/2},H_{-1/2}}\) denotes the dual pairing between \(H_{1/2}(L)\) and \(H_{-1/2}(L)\). We have \(B_2=b_2\), \(C_1^*= c_1\) and \(C_2^*= c_2\), and \(\Vert D_0\Vert =\Vert d\Vert _{L_a^2}\). Since \(-L\) is positive and its smallest eigenvalue is given by \(-\mu _1 = a^2/4+\pi ^2\), we also have \(\Vert (-L)^{-1/2}\Vert =(-\mu _1)^{-1/2} = (a^2/4+\pi ^2)^{-1/2}\). Thus the claim follows from Corollary 2.8. \(\square \)
As shown in Proposition 4.2 the Webster’s equation with the damping function \(d(x)=1-x\) is polynomially stable with \(\alpha =2\). Since \(\Vert d\Vert _{L_a^2}^2=2a^{-3}(e^a-1-a-a^2/2)\) for this d, for the choices \(\beta =\gamma =1\) Theorem 4.4 has the following form.
Corollary 4.5
Let \(d(x)=1-x\). If \(\kappa >0\) is as in Theorem 2.2 with \(\beta =\gamma =1\) and if \( b_2,c_2\in H_{1/2}(L)\), \(c_1\in L_a^2(0,1)\) satisfy
where \(M=1+ 2a^{-3}(e^a-1-a-a^2/2)(\frac{a^2}{4}+\pi ^2)^{-1/2},\) then the perturbed Webster’s Eq. (4.2) is polynomially stable with \(\alpha =2\).
Example 4.6
We use Theorem 2.3 for computing an explicit numerical value of \(\kappa \) for the case \(d(x) = 1-x\) and \(a = 2\). To this end, we need to find functions \(\eta (\cdot )\) and \(\delta (\cdot )\) such that the condition (2.5) in Theorem 2.3 is satisfied. For \(a=2\) the eigenvalues and the corresponding eigenvectors of \(A_0\) are \(\lambda _{n}=\mathrm{sign\,}(n) i\sqrt{1+\pi ^2n^2}\),
For all \(n\in \mathbb {N}\) (using the inequality \(\sqrt{x + y}\le \sqrt{x}+ \sqrt{y}\))
since \(f(x) = \frac{\pi ^2x}{2 + \pi x}\) is increasing for \(x\in (1,\infty )\). If we choose \(\delta (s) \equiv \delta _0 = \frac{\pi ^2}{a + 3\pi }\), then every interval \((i(s-\delta _0),i(s+\delta _0))\) contains at most one eigenvalue and \(\mathrm{ran\,}P_{(s-\delta _0,s+\delta _0)}\) consists of the corresponding eigenvector. Similar computations as in the proof of Proposition 4.2 then show that
where \(\eta (s) = \frac{c}{s+\delta _0} \) with the \(c>0\) such that
To find a suitable \(c>0\), let us denote
and
It is obvious that \(F(n)\ge G(n)\) for \(n\ge 2\). The values G(n) for \(n\ge 2\) are increasing and \(G(n)\rightarrow 1\) as \(n\rightarrow \infty \) and therefore \(\min \limits _{n\ge 2} G(n) = G(2)\). Hence we can choose
Finally, the maximum of \(\eta (s)\) when \(s\ge 0\) is \(\eta _0 = c/\delta _0\).
In the next step we calculate \(M_R\). To this end, we need also \(\Vert D\Vert \) which is
We can now use Matlab to compute \(M_R = 5.451\).
Now we will find the constant \(M_0>0\). A direct estimate using \((x+y)^2\le 2(x^2+y^2)\) and \(\delta _0<1\) shows that
with \(M_0=\frac{2}{c^2\delta _0^2}\). To compute \(M_C>0\), we also need an estimate for \(\Vert A^{-1}\Vert \). We have
If we take \(s_0 = 2.8\) in the formula for \(M_C\), we obtain \(M_C = 17.0664\). This way, we finally see that \(\kappa >0\) in Theorem 2.2 can take any value such that \(\kappa < \frac{1}{\sqrt{2M_C}}= 0.1712\).
Now we are able to give explicit upper bounds for the norms of \(b_2\), \(c_1\), and \(c_2\) for the preserving of polynomial stability. From Corollary 4.5 we have that \( b_2,c_2\in H_{1/2}(L)\), \(c_1\in L_a^2(0,1)\) satisfy
where \(M = 1 + \frac{e^2-5}{4\sqrt{1+\pi ^2}}\), then the perturbed Webster’s Eq. (4.2) is polynomially stable with \(\alpha =2\).
5 Wave Equation with an Acoustic Boundary Condition
In this section we consider a one-dimensional wave equation with an "acoustic boundary condition" on the interval \(\Omega =(0,1)\),
with \(k,d>0\) [1, Sec. 6.1]. The spectral properties and polynomial stability of differential equations of this form (also on multidimensional spatial domains) have been studied in detail in [1, 2, 9, 22]. In particular, it was shown in [22, Thm. 1.3] that the energy of the classical solutions of (5.1) decays at a rational rate, and the optimality of this decay rate was proved in [1, Sec. 6.1]. This model is not of the form (1.1), but the preservation of its polynomial stability can be studied using Theorem 2.2.
Equation (5.1) can be formulated as an abstract Cauchy problem with state \(u(t) = ( w_x(t,\cdot ), w_t(t,\cdot ), a(t), a_t(t))^\top \) on the Hilbert space \( {{\mathcal {H}}}=L^2(0,1)\times L^2(0,1)\times \mathbb {C}^2 \) with inner product defined as
for all \(u=(u_1(\cdot ),u_2(\cdot ),u_3,u_4)^\top \), \(v=(v_1(\cdot ),v_2(\cdot ),v_3,v_4)^\top \in {{\mathcal {H}}}\). In this situation the semigroup generator is defined as
with domain
The operator A generates a contraction semigroup on \({{\mathcal {H}}}\), and it was shown in [22, Thm. 1.3] (see also [1, Sec. 6.1], [26, Sec. 4]) that this semigroup is polynomially stable with \(\alpha =2\). In the context of the wave Eq. (5.1) this means that there exists a constant \(M_T>0\) such that for all initial conditions \(w_0,w_1,a_0,a_1\) such that \((w_0',w_1,a_0,a_1)^\top \in \mathrm{dom\,}A \) the solutions of (5.1) satisfy
for all \(t>0\).
We can now study the stability of perturbed wave equations of the form
where \(b_2,c_1,c_2\in L^2(0,1)\) and \(c_3,c_4\in {{\mathbb {C}}}\). The following two theorems introduce conditions for the polynomial stability of (5.3).
Theorem 5.1
Assume \(\kappa >0\) is as in Theorem 2.2 with \(\beta =\gamma =1\). If \(b_2\in H^1_0(0,1)\), \(c_1,c_2\in H^1(0,1)\), and \(c_4=0\) satisfy
then the perturbed Eq. (5.3) is polynomially stable with \(\alpha =2\). For such perturbations there exists \(M_T>0\) such that the solutions of (5.3) corresponding to initial conditions \(w_0,w_1,a_0,a_1\) such that \((w_0',w_1,a_0,a_1)^\top \in \mathrm{dom\,}A \) satisfy (5.2) for all \(t>0\).
Proof
The perturbed system operator can be written as \(A+BC\) where \(B=(0,b_2(\cdot ),0,0)^\top \in {{{\mathfrak {L}}}}({{\mathbb {C}}},{{\mathcal {H}}})\) and \(C=(\langle \cdot ,c_1(\cdot )\rangle _{L^2},\langle \cdot ,c_2(\cdot )\rangle _{L^2},\overline{c_3},0)\in {{{\mathfrak {L}}}}({{\mathcal {H}}},{{\mathbb {C}}})\) with \(b_2\in H_0^1(0,1)\), \(c_1,c_2\in H^1(0,1)\) and \(c_3 \in {\mathbb {C}}\). A straightforward computation shows that the adjoint operator of A has the form
and that its domain \(\mathrm{dom\,}A^*\) contains the subspace
The assumptions therefore imply that \(\mathrm{ran\,}B\subset \mathrm{dom\,}A\) and \(\mathrm{ran\,}C^*\subset \mathrm{dom\,}A^*\), and
Here we have used the property \(|c_1(1)|\le \Vert c_1\Vert _{L^2}+\Vert c_1'\Vert _{L^2}\), which can be verified using the identity \(c_1(1)=\int _0^1 \frac{d}{dx}(xc_1(x))dx = \int _0^1(c_1(x)+xc_1'(x))dx\). Thus the claim follows from Theorem 2.2 with \(\beta =\gamma =1\). \(\square \)
Similarly, applying Theorem 2.2 with \(\beta =2\) and \(\gamma =0\) we obtain the following alternative conditions for the polynomial stability of (5.3).
Theorem 5.2
Assume \(\kappa >0\) is as in Theorem 2.2 with \(\beta =2\) and \(\gamma =0\). If \(b_2\in H_0^1(0,1)\cap H^2(0,1)\), \(c_1,c_2\in L^2(0,1)\), and \(c_3,c_4\in {{\mathbb {C}}}\) satisfy
then (5.3) is polynomially stable with \(\alpha =2\). For such perturbations there exists \(M_T>0\) such that the solutions of (5.3) corresponding to initial conditions \(w_0,w_1,a_0,a_1\) such that \((w_0',w_1,a_0,a_1)^\top \in \mathrm{dom\,}A \) satisfy (5.2) for all \(t>0\).
Proof
The perturbations have the same form as in the proof of Theorem 5.1. Since \(B=b:=(0,b_2(\cdot ),0,0)^\top \) with \(b_2\in H_0^1(0,1)\cap H^2(0,1)\), we have \((0,b_2(\cdot ),0,0)^\top \in \mathrm{dom\,}A\) and \(Ab=(b'_2(\cdot ),0,0,0)^\top \in \mathrm{dom\,}A\). Thus \(\mathrm{ran\,}B\subset \mathrm{dom\,}A^2\) and
implies \(\Vert A^2B\Vert ^2=\Vert b_2''\Vert _{L^2}^2+|b_2'(1)|^2\le 3\Vert b_2'\Vert _{H^1}^2\), since \(|b_2'(1)|^2\le 2\Vert b_2'\Vert _{H^1}^2\) similarly as in the proof of Theorem 5.1. The claim now follows from Theorem 2.2 with the choices \(\beta =2\) and \(\gamma =0\). \(\square \)
References
Abbas, Z., Nicaise, S.: Polynomial decay rate for a wave equation with general acoustic boundary feedback laws. SeMA J. 61, 19–47 (2013)
Abbas, Z., Nicaise, S.: The multidimensional wave equation with generalized acoustic boundary conditions II: polynomial stability. SIAM J. Control Optim. 53(4), 2582–2607 (2015)
Ammari, K., Bchatnia, A., El Mufti, K.: Non-uniform decay of the energy of some dissipative evolution systems. Z. Anal. Anwend. 36(2), 239–251 (2017)
Ammari, K., Tucsnak, M.: Stabilization of second order evolution equations by a class of unbounded feedbacks. ESAIM Control Optim. Calc. Var. 6, 361–386 (2001)
Anantharaman, N., Léautaud, M.: Sharp polynomial decay rates for the damped wave equation on the torus. Anal. PDE 7(1), 159–214 (2014). With an appendix by Stéphane Nonnenmacher
Avalos, G., Triggiani, R.: Rational decay rates for a PDE heat-structure interaction: a frequency domain approach. Evol. Equ. Control Theory 2(2), 233–253 (2013)
Bátkai, A., Engel, K.-J., Prüss, J., Schnaubelt, R.: Polynomial stability of operator semigroups. Math. Nachr. 279(13–14), 1425–1440 (2006)
Batty, C., Chill, R., Tomilov, Y.: Fine scales of decay of operator semigroups. J. Eur. Math. Soc. 18(4), 853–929 (2016)
Beale, J.T.: Spectral properties of an acoustic boundary condition. Indiana Univ. Math. J. 25(9), 895–917 (1976)
Borichev, A., Tomilov, Y.: Optimal polynomial decay of functions and operator semigroups. Math. Ann. 347(2), 455–478 (2010)
Burq, N., Hitrik, M.: Energy decay for damped wave equations on partially rectangular domains. Math. Res. Lett. 14(1), 35–47 (2007)
Burq, N., Zworski, M.: Rough controls for Schrödinger operators on 2-tori. Ann. H. Lebesgue 2, 331–347 (2019)
Chill, R., Paunonen, L., Seifert, D., Stahn, R., Tomilov, Y.: Non-uniform stability of damped contraction semigroups. Anal. PDE, accepted for publication. Preprint available at arXiv:1911.04804
Datchev, K., Kleinhenz, P.: Sharp polynomial decay rates for the damped wave equation with Hölder-like damping. Proc. Am. Math. Soc. 148(8), 3417–3425 (2020)
Dell’Oro, F., Pata, V.: Second order linear evolution equations with general dissipation. Appl. Math. Optim. 83(3), 1877–1917 (2021)
Duyckaerts, T.: Optimal decay rates of the energy of a hyperbolic-parabolic system coupled by an interface. Asymptot. Anal. 51(1), 17–45 (2007)
Engel, K.-J., Nagel, R.: One-Parameter Semigroups for Linear Evolution Equations. Springer, New York (2000)
Jaffard, S.: Contrôle interne exact des vibrations d’une plaque rectangulaire. Portugal. Math. 47(4), 423–429 (1990)
Kato, T.: A generalization of the Heinz inequality. Proc. Jpn. Acad. 37, 305–308 (1961)
Liu, Zhuangyi, Rao, Bopeng: Characterization of polynomial decay rate for the solution of linear evolution equation. Z. Angew. Math. Phys. 56(4), 630–644 (2005)
Liu, Zhuangyi, Zhang, Qiong: A note on the polynomial stability of a weakly damped elastic abstract system. Z. Angew. Math. Phys. 66(4), 1799–1804 (2015)
Muñoz Rivera, J.E., Qin, Yuming: Polynomial decay for the energy with an acoustic boundary condition. Appl. Math. Lett. 16(2), 249–256 (2003)
Paunonen, L.: Robustness of strong and polynomial stability of semigroups. J. Funct. Anal. 263, 2555–2583 (2012)
Paunonen, L.: Robustness of strong stability of semigroups. J. Differential Equations 257, 4403–4436 (2014)
Paunonen, L.: Robust controllers for regular linear systems with infinite-dimensional exosystems. SIAM J. Control Optim. 55(3), 1567–1597 (2017)
Paunonen, L.: On polynomial stability of coupled partial differential equations in 1D. In: Semigroups of Operators—Theory and Applications. Springer Proc. Math. Stat. 325, pp. 349–362 (2020)
Muñoz Rivera, J.E., Ávila, A.I.: Rates of decay to non homogeneous Timoshenko model with tip body. J. Differ. Equ. 258(10), 3468–3490 (2015)
Stahn, R.: Optimal decay rate for the wave equation on a square with constant damping on a strip. Z. Angew. Math. Phys. 68(2), 36 (2017)
Tucsnak, M., Weiss, G.: Observation and Control for Operator Semigroups. Birkhäuser Basel (2009)
Zhang, Xu., Zuazua, E.: Long-time behavior of a coupled heat-wave system arising in fluid-structure interaction. Arch. Ration. Mech. Anal. 184(1), 49–120 (2007)
Open Access
This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The research was funded by the Academy of Finland Grants 298182 and 310489 held by L. Paunonen.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Baidiuk, D., Paunonen, L. Robustness of Polynomial Stability of Damped Wave Equations. J Dyn Diff Equat 35, 1489–1510 (2023). https://doi.org/10.1007/s10884-021-10117-y
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10884-021-10117-y