🔢 mathematics

Optimal Polynomial Stabilization of the Linearized Periodic Whitham--Boussinesq System

This paper establishes the well-posedness and optimal polynomial stability of the linearized periodic Whitham--Boussinesq system under localized damping, proving a $t^{-2}$ energy decay rate that is shown to be sharp through the construction of high-frequency quasimodes.

Original authors: Roberto de A. Capistrano Filho (DMat/UFPE), William Artiles Roqueta (DMat/UFPE)

Published 2026-08-27

📖 4 min read🧠 Deep dive

CC0 1.0

Original authors: Roberto de A. Capistrano Filho (DMat/UFPE), William Artiles Roqueta (DMat/UFPE)

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Water waves are a familiar sight, from the gentle ripples in a pond to the powerful swells of the ocean. When scientists try to predict how these waves move, they often rely on mathematical models that simplify the complex physics of fluid motion. One such model, known as the Whitham–Boussinesq system, is particularly valuable because it captures the exact way water waves spread out and disperse, a feature that older, simpler models often miss. This system describes how waves travel in both directions along a channel of finite depth, balancing the push of the water's weight against the fluid's tendency to spread. While these equations are excellent for describing how waves form and travel, a different question arises when we consider how to stop them. In the real world, waves eventually lose energy due to friction or obstacles, but in a mathematical model, we must introduce a mechanism to represent this loss. The challenge is to understand exactly how fast a wave system will calm down if we apply a damping force, and whether that force needs to be applied everywhere or if a small, localized patch is enough to stop the motion entirely.

In a recent study, researchers tackled this problem by examining a simplified, linear version of the Whitham–Boussinesq system on a repeating loop, essentially a circular channel where waves can travel forever without hitting a wall. Their goal was to determine the precise rate at which energy dissipates when a damping force is applied to the system. They found that the system does not settle down quickly in an exponential fashion, where the energy would drop by a fixed percentage every second. Instead, the energy fades away much more slowly, following a polynomial decay. This means the energy decreases at a rate proportional to the inverse square of time; if you wait twice as long, the energy is only one-quarter of what it was, rather than a tiny fraction. The researchers proved that this specific rate of decay is the best possible outcome for this type of system when the damping is applied only to a specific, limited region of the channel.

The key to this discovery lies in the unique way high-frequency waves behave in this model. In many physical systems, waves of different frequencies are spaced out in a regular pattern, which allows for predictable and rapid stabilization. However, in the Whitham–Boussinesq system, the frequencies of these waves grow in a sublinear way, meaning that as the waves get smaller and faster, their frequencies do not increase as rapidly as one might expect. This creates a situation where the high-frequency waves are harder to control. The researchers showed that if the damping force is not applied everywhere, these fast, small waves can hide in the undamped regions, slipping through the net of the control mechanism. By constructing specific mathematical examples of these elusive waves, they demonstrated that the system's ability to stabilize is fundamentally limited by this high-frequency behavior.

To confirm that their calculated decay rate was not just a guess but the absolute best possible, the team used a sophisticated mathematical tool that links the behavior of the system's waves to the speed at which it settles. They proved that the system's response to external forces grows linearly with frequency, a property that directly translates to the slow, polynomial decay of energy. This result is significant because it establishes a hard limit: no matter how cleverly one designs the damping mechanism, as long as it is confined to a specific area, the energy cannot disappear faster than this polynomial rate. The study rules out the possibility of exponential stability for this system under localized damping, showing that the slow fade is an intrinsic feature of the physics, not a flaw in the method.

The work also highlights the delicate balance between the geometry of the damping and the nature of the waves. If the damping were applied everywhere on the loop, the outcome might be different, potentially allowing for faster stabilization. However, under the realistic assumption that damping is only present in a specific zone, the system is forced to rely on the slow process of waves eventually drifting into the damped area to lose their energy. The researchers' findings provide a clear, quantitative answer to how long one must wait for the waves to die down, offering a precise benchmark for engineers and scientists working with similar fluid dynamics problems. Ultimately, the study reveals that for this class of water wave models, patience is a mathematical necessity; the energy will vanish, but it will do so at a steady, predictable, and relatively slow pace.

Technical Summary: Optimal Polynomial Stabilization of the Linearized Periodic Whitham–Boussinesq System

Problem Formulation
This study investigates the long-time stabilization properties of the linearized periodic Whitham–Boussinesq system on the one-dimensional torus $\mathbb{T} = \mathbb{R}/(2\pi\mathbb{Z})$ . The system models bidirectional surface water waves, retaining the exact linear dispersion relation of the full water-wave problem while employing a weakly nonlinear Boussinesq-type approximation. The authors focus exclusively on the linearized dynamics around the zero equilibrium, governed by the system:
$\begin{cases} h_t - T_\beta u = 0, \\ u_t + h_x = 0, \end{cases}$
where $h$ and $u$ represent the surface elevation and velocity, respectively, and $T_\beta$ is a nonlocal Fourier multiplier defined by the symbol $-i \frac{\tanh(\beta k)}{\beta}$ . The natural energy space is identified as $\mathcal{H} = L^2(\mathbb{T}) \times H^{-1/2}(\mathbb{T})$ .

To induce stabilization, the authors introduce spatially dependent damping terms, leading to the closed-loop system:
$\begin{cases} h_t - T_\beta u + a^2(x)h = 0, \\ u_t + h_x + b^2(x)P_\beta u = 0, \end{cases}$
where $a, b \in C^\infty(\mathbb{T})$ are non-negative, non-trivial damping coefficients, and $P_\beta$ is a positive self-adjoint Fourier multiplier of order $-1$. The central question is whether this dissipation forces the energy to decay to zero and, if so, to determine the optimal rate of decay.

Methodology and Proof Strategy
The analysis proceeds through a rigorous spectral and resolvent-based approach, avoiding classical non-harmonic Fourier methods due to the specific high-frequency behavior of the system.

Well-posedness and Functional Framework: The authors first establish that the conservative operator $A$ is skew-adjoint in $\mathcal{H}$ , generating a unitary $C_0$ -group. The damped operator $A_d = A - D$ (where $D$ is the damping operator) is shown to be maximal dissipative, generating a contraction semigroup $S_d(t)$ .
Spectral Analysis: The spectrum of the conservative operator is characterized as purely discrete, with eigenvalues $\pm i\omega_k$ where the frequencies satisfy $\omega_k \sim \beta^{-1/2}|k|^{1/2}$ as $|k| \to \infty$ . This sublinear growth ( $|k|^{1/2}$ ) is critical, as it violates the uniform gap condition required for classical exponential stability proofs.
Strong Stability: Using the Arendt–Batty–Lyubich–Vũ theorem, the authors prove that the damped semigroup is strongly stable (energy decays to zero) provided the damping coefficients are non-trivial. This relies on showing that the damped generator has no eigenvalues on the imaginary axis.
Resolvent Estimates: The core of the analysis involves determining the growth of the resolvent $(i\tau I - A_d)^{-1}$ $(i τ I - A_{d})^{- 1}$ as $|\tau| \to \infty$ $∣ τ ∣ \to \infty$ .
- Upper Bound: By reducing the coupled system to a scalar operator via the transformation $v = P_\beta^{1/2}u$ , the authors adapt high-frequency stationary estimates for damped fractional wave equations (specifically referencing Alazard, Marzuola, and Wang). They prove that $\|(i\tau I - A_d)^{-1}\|_{\mathcal{L}(\mathcal{H})} \lesssim |\tau|$ .
- Lower Bound: Under the assumption that the damping on the first component ( $a$ ) is "genuinely localized" (i.e., vanishes on a non-empty open set), the authors construct high-frequency quasimodes concentrated in the undamped region. These quasimodes demonstrate that the resolvent growth is at least linear, i.e., $\|(i\tau I - A_d)^{-1}\|_{\mathcal{L}(\mathcal{H})} \gtrsim |\tau|$ .
Decay Rates: The linear resolvent growth estimate is combined with the Borichev–Tomilov theorem to deduce polynomial decay rates for the semigroup and the energy.

Key Results
The main result (Theorem 1.1) establishes that for initial data in the domain of the generator $D(A_d)$ :

The solution decays as $\|Z(t)\|_{\mathcal{H}} \leq C t^{-1} \|Z_0\|_{D(A_d)}$ .
The energy decays as $E(t) \leq C t^{-2} \|Z_0\|_{D(A_d)}^2$ .

Furthermore, if the damping $a$ is genuinely localized, these polynomial rates are optimal. The paper proves that no faster polynomial decay rate (e.g., $t^{-p}$ with $p > 1$ ) is possible, as this would contradict the linear lower bound on the resolvent derived from the quasimode construction. Consequently, exponential stability is impossible under genuine localization.

Significance and Claims
The paper claims to provide the first optimal polynomial stabilization result for the linearized periodic Whitham–Boussinesq system. Its significance lies in:

Handling Sublinear Dispersion: It addresses the specific challenge posed by the $|k|^{1/2}$ dispersion relation, which prevents exponential decay even with localized damping, a phenomenon distinct from classical wave equations where geometric control conditions often yield exponential stability.
Optimality: It rigorously demonstrates that the $t^{-2}$ energy decay rate is the best possible polynomial rate under the given damping assumptions, linking the decay rate directly to the high-frequency spectral structure of the nonlocal operator.
Methodological Adaptation: The work successfully adapts techniques from fractional wave equations (stationary resolvent estimates) to a coupled, nonlocal dispersive system, bridging the gap between full-dispersion models and control theory.

The authors conclude that while the linear analysis is complete, the persistence of these polynomial rates in the nonlinear regime and the influence of damping geometry (e.g., non-vanishing damping everywhere) remain open questions for future investigation.

More like this