Boundary stabilization of a vibrating string with variable length
Abstract.
We study small vibrations of a string with time-dependent length and boundary damping. The vibrations are described by a 1-d wave equation in an interval with one moving endpoint at a speed slower than the speed of propagation of the wave c=1. With no damping, the energy of the solution decays if the interval is expanding and increases if the interval is shrinking. The energy decays faster when the interval is expanding and a constant damping is applied at the moving end. However, to ensure the energy decay in a shrinking interval, the damping factor must be close enough to the optimal value , corresponding to the transparent condition. In all cases, we establish lower and upper estimates for the energy with explicit constants.
Key words and phrases:
Wave equation, time-dependent domains, generalized Fourier series, boundary stabilization.2010 Mathematics Subject Classification:
35L05, 35C10, 93D15.1. Introduction
We consider small transversal vibrations of a uniform string, with a time dependent length. The mechanical setting is sketched in Figure 1 where the left end of the string is fixed while the right moving end is also allowed to move transversely and attached to a damping device (a dash-pot with a damping factor ).

Denoting the displacement function by , depending on the position along the string and the time , the model can be stated as follows
(WP) |
The subscripts and in (WP) stand for the derivatives in time and space variables respectively. The functions and represents the initial shape and the initial transverse speed of the string, respectively. The initial length of the string is denoted by .
We assume that and that
(1.1) |
which means that the speed of variation of the length of the string is strictly less then the speed of propagation of the wave (here simplified to ).
In this paper, we are mainly interested in the asymptotic behaviour in time of the energy of the solution, defined as
(1.2) |
For the time-independent interval, i.e. when for , it is well known that:
-
•
If , then it is the energy is constant, i.e. for .
- •
The question of behaviour in time is more delicate if the interval depends on time. For instance, for the Dirichlet boundary conditions (which corresponds to in (WP)), the energy decays if the interval is expanding and increases if the interval is shrinking, see [2]. See also [12, 13, 14] when the variation of the length is uniform in time.
Regarding the case with a velocity feedback at the moving endpoint :
- •
- •
For the special case , the boundary condition at reads
This is a transparent condition, i.e. there is no reflections of waves from the moving endpoint and consequently all the initial disturbances leave the interval at most after a time
where see Figure 2. Hence, wether the interval is expanding or shrinking, the linear velocity feedback steers the solution to the zero state in the finite time . See for instance [6, 7], and for the particular case see [15, 4]. In the remaining of this paper, we will assume that and

The approach used in the paper at hands is based on generalized Fourier series. This is possible since the closed form for the solution of (WP) is given by the series
(1.4) |
where the coefficients can be explicitly computed in function of the initial data We denoted by a sequence of complex number given by
(1.5) |
Observe that since for every the real part of is nonpositive.
By we denoted a real function satisfying the functional equation
(1.6) |
This equation often called Moor’s equation following his paper [10]. Some pairs of solutions can be found in [17].
We will assume that is differentiable and increasing. More precisely,
(1.7) |
The assumption is needed in the sequel since will serve as a weight for an space. Besides, a deacreasing can not satisfy (1.6) since . See [7] for further discussions on the regularity of the solutions of (1.6).
In this work, we demonstrate how the series formulas (1.4) can be used to achieve the following results:
-
•
For the undamped case, i.e. , the energy of the solution satisfies
(1.8) where
(1.9) See Theorem 2 and its corollaries for sharper estimates.
- •
After the present introduction, we derive the exact solution and the expression for the coefficients of the series formula (1.4). In Section 3, we establish upper and lower estimates for the energy of the undamped equation. In Section 4, we deal with the damped case. Some examples will be included as a last section.
2. Exact solution
Let us introduce the following family of Hilbert spaces
and assume that the initial data satisfies
(2.1) |
Let Then, we have the following existence result for Problem (WP).
Theorem 1.
Proof.
The exact solution: This part of solution is slightly different from the approach in [16] where the author considered instead of in the boundary condition at . We include it here for the sake clarity. The general solution of (WP) is given by D’Alembert’s formula
(2.5) |
where and are arbitrary continuous functions. The boundary conditions at the endpoint , we have
The condition at implies that
hence
(2.6) |
where and . Then, noting that
(2.7) |
we can rewrite (2.6) as
(2.8) |
By integration, it follows that
(2.9) |
Let us assume for the moment that Then, it is convenient to search for in the form , for a constant and some function . Substituting in (2.9), we get
Assuming that satisfies (1.6), we are led to the following cases:
-
-
If , then and we get
Solving this equation for , we obtain a sequence of values where
-
-
If , we have and we obtain this time
Thus, if and , we always have and given by (1.5).
Due to the superposition principal, it follows that can be written as
where are complex coefficients to be determined later. Since then D’Alembert’s formula for the solution yields the series
(2.10) |
If in (2.9), then we can check that
solves (2.9). However, this will not affect the solution of (WP) since and thus the constant parts of and will be cancelled in the expression (2.10).
Computing the coefficients : We extend to an odd function on the interval This ensures in particular that the boundary condition at is satisfied for It follows also that and are respectively an odd and an even function. Going back to (2.10), we infer that
for and . Using the definition of we get
(2.11) |
which implies that
(2.12) |
Taking into account that is an orthonormal basis of the weighted space we deduce that
for . Whether or in both cases, we have
(2.13) |
Taking , we obtain (2.3) as claimed.
3. The undamped case
In this section, we show some results for the undamped case, i.e. in Problem (WP). To know wether the energy is increasing or in creasing, we compute Thus, using Leibniz’s rule for differentiation under the integral sign, we get
Since , then and it follows
(3.1) |
Then, we have the following result.
Lemma 1.
The energy of solution of Problem (WP), with , satisfies
Proof.
Remark 1.
The above lemma means that
(3.2) |
The same result holds for the multidimensional wave equation, in a time-dependent domain, with homogenous Dirichlet boundary conditions, see [2].
The next theorem show that the asymptotic behaviour of is dictated by
Theorem 2.
Proof.
If then and . The identity (2.14) becomes
(3.5) |
Since is an even function of and that is an odd one, then changing by in the last formula, we also obtain
(3.6) |
Remark 2.
If is monotone then the asymptotic behaviour is dictated by and we have the following refinements.
Corollary 1.
Under the assumption of Theorem 2, assume that
(3.9) |
-
•
If on , then and are nondeacreasing and
(3.10) -
•
If on , then and are nonincreasing and
(3.11)
Proof.
We already have (3.2). The derivation of the identity (1.6) yields
(3.12) |
Taking into consideration the variation of the function on the interval , it follows that
hence
(3.13) | |||
(3.14) |
Of course, if is monotone and then (3.13) means that is necessarily nondecreasing, and so is and The same argument can be made when . This shows the corollary. ∎
Remark 3.
The assumption (3.9) is satisfied in all the examples of the last section.
Recall that can be computed using the initial data by setting in (3.4). If one needs to compare with for then we have the next result.
-
•
If on , then and are nondecreasing and satisfy
(3.16) -
•
If on , then and are nonincreasing and
(3.17)
Proof.
4. The damped case
In this section, we investigate the case with boundary damping when
in Problem (WP). Let us recall that we still have (3.1), i.e.
but now the boundary condition at is
Let us discuss the sign of for different values of
-
•
If the interval is shrinking for where then we have the following cases:
- If then the boundary condition at reads and thus
- If then and thus
- Assume that then after some computation we can rewrite (3.1) as
(4.1) The sign of is opposite to the sign of
Due to (1.1), the polynomial has a discriminant and thus has two real roots
(4.2) Figure 3. Variation of and in function of when . -
–
If for then and by consequence
-
–
If or for then i.e. the energy is constant.
-
–
If for then and we have
-
–
-
•
If the interval is independent of time, i.e. for then and thus
- •
Remark 4.
To summarize, under the assumption
(4.3) |
Let us now estimate is using and .
Theorem 3.
Proof.
We argue as in the proof of Theorem 2. Noting that is an even function of and that is an odd one, then the identity (2.14) yields
(4.6) |
Hence, summing, we get
Expanding squares, we get
(4.7) |
As the function under the integral sign is even, then (4.4) follows.
For and let us denote
Then, we can rewrite (4.4) as
Using the algebraic inequality
we get
Thanks to (3.8), the precedent estimation yields
for By consequence
This implies (4.5). ∎
Since for and is nondecreasing, we have the following immediate corollary.
If is monotone, then (4.8) can be replaced by more explicit estimation.
-
•
If on , then is nondecreasing and satisfies
(4.9) -
•
If on , then and are nonincreasing and
(4.10)
Proof.
It suffices to argue as in the proof of Corollary 1. ∎
Remark 5.
To compare with the energy for we have the following result.
Corollary 5.
-
•
If on , then is nondecreasing and satisfies
(4.12) -
•
If on , then and are nonincreasing and
(4.13)
Proof.
Acknowledgements
The authors have been supported by the General Direction of Scientific Research and Technological Development (Algerian Ministry of Higher Education and Scientific Research) PRFU # C00L03UN280120220010.
ORCID
Abdelmouhcene Sengouga https://orcid.org/0000-0003-3183-7973.
References
- [1] K. Ammari, A. Bchatnia, and K. El Mufti. Stabilization of the wave equation with moving boundary. Eur. J. Control, 39:35–38, 2018. ISSN 0947-3580.
- [2] C. Bardos and G. Chen. Control and stabilization for the wave equation. III: Domain with moving boundary. SIAM J. Control Optim., 19:123–138, 1981.
- [3] M. Cherkaoui. Estimation optimale du taux de décroissance de l’énergie pour une équation des ondes avec contrôle frontière. Rapport de recherche N°2328 - INRIA, 1994.
- [4] S. Cox and E. Zuazua. The rate at which energy decays in a string damped at one end. Indiana Univ. Math. J., pages 545–573, 1995.
- [5] S. E. Ghenimi and A. Sengouga. Energy decay estimates for an axially travelling string damped at one end. Submitted. URL http://arxiv.org/abs/2211.10537. 14 pages.
- [6] M. Gugat. Optimal boundary feedback stabilization of a string with moving boundary. IMA J. Math. Control Inform., 25(1):111–121, 2008.
- [7] B. H. Haak and D.-T. Hoang. Exact observability of a 1-dimensional wave equation on a noncylindrical domain. SIAM J. Control Optim., 57(1):570–589, 2019.
- [8] L. Lu and Y. Feng. Observability and stabilization of wave equations with moving boundary feedback. Acta Appl. Math., 170(1):731–753, 2020.
- [9] Y. Mokhtari. Boundary controllability and boundary time-varying feedback stabilization of the 1d wave equation in non-cylindrical domains. Evol. Equ. Control Theory, 11(2):373–397, 2022.
- [10] G. T. Moore. Quantum theory of electromagnetic field in a variable-length one-dimensional cavity. J. Math. Phys., 11(9):2679–2691, 1970.
- [11] J. P. Quinn and D. L. Russell. Asymptotic stability and energy decay rates for solutions of hyperbolic equations with boundary damping. Proc. R. Soc. Edinb. A: Math., 77(1-2):97–127, 1977.
- [12] A. Sengouga. Observability of the 1-D wave equation with mixed boundary conditions in a non-cylindrical domain. Mediterr. J. Math., 15(62):1–22, 2018.
- [13] A. Sengouga. Observability and controllability of the 1-D wave cquation in domains with moving boundary. Acta Appli. Math., 157:117–128, 2018.
- [14] H. Sun, H. Li, and L. Lu. Exact controllability for a string equation in domains with moving boundary in one dimension. Electron. J. Diff. Equations, 2015(98):1–7, 2015.
- [15] K. Veselić. On linear vibrational systems with one dimensional damping. Appl. Anal., 29(1-2):1–18, 1988.
- [16] A. I. Vesnitskii. The inverse problem for a one-dimensional resonator the dimensions of which vary with time. Radiophysics and Quantum Electronics, 14(10):1209–1215, 1971.
- [17] A. I. Vesnitskii and A. I. Potapov. Some general properties of wave processes in one-dimensional mechanical systems of variable length. Sov. App. Mechanics, 11(4):422–426, 1975. ISSN 1573-8582.