🔢 mathematics

On decay of solutions to the anisotropic Boussinesq equations near the hydrostatic balance in half space $\mathbb{R}_+^3$

This paper establishes the global stability and derives decay rates for perturbations of the 3D incompressible anisotropic Boussinesq system with horizontal dissipation near hydrostatic balance in the half-space $\mathbb{R}_+^3$ by employing energy methods, bootstrapping arguments, and specific Fourier transforms.

Original authors: Wangrong Yang, Aibin Zang

Published 2026-01-27

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: Wangrong Yang, Aibin Zang

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

Imagine a vast, invisible ocean of air or water stretching infinitely to the sides and upwards, but with a solid floor at the bottom. This is the setting for the research paper you provided. The scientists, Wanrong Yang and Aibin Zang, are studying how this fluid moves and changes temperature when it's slightly disturbed from a state of perfect calm.

Here is a breakdown of their work using simple analogies:

1. The Setup: A Calm Lake with a Twist

The paper focuses on the Boussinesq equations. Think of these as the "rulebook" for how fluids (like air in the atmosphere or water in the ocean) move when heat is involved.

The "Hydrostatic Balance": Imagine a perfectly still lake where the water is slightly warmer at the bottom and cooler at the top. It's stable; nothing is moving. This is the "hydrostatic balance."
The Disturbance: Now, imagine dropping a pebble into this lake or giving the water a tiny nudge. The paper asks: If we nudge this fluid just a little bit, will it eventually settle back down, or will it go wild and crash?

2. The Special Rules of the Game

This isn't just any fluid; it has two special constraints that make the math tricky:

One-Way Friction (Anisotropic Dissipation): Usually, fluids slow down because of friction in all directions (like stirring honey). In this model, the fluid only feels "friction" or resistance when moving side-to-side (horizontally). It's like a car that has brakes on its left and right wheels but no brakes on its front and back. The fluid can slide up and down freely without slowing down, which makes predicting its behavior much harder.
The Floor Rules (Boundary Conditions): The fluid is sitting on a floor (the ground). The rules say the fluid can't go through the floor, and it can't spin wildly right against the floor. It's a "slippery" floor where the fluid can glide along but not stick or penetrate.

3. The Main Discovery: Stability

The authors proved a very important thing: If the initial nudge (the pebble) is small enough, the fluid will always calm down.

The Guarantee: They showed that no matter how long you wait, the fluid won't explode into chaos. It will return to a stable state.
The "Energy" Check: They used a mathematical tool called "energy methods." Imagine the fluid has a certain amount of "wiggle energy." They proved that even though the fluid is slippery in the vertical direction, the horizontal friction is strong enough to drain that wiggle energy away over time, eventually bringing the system to rest.

4. The Speed of Calming Down (Decay)

The second part of the paper is like a stopwatch. They didn't just prove the fluid would calm down; they calculated how fast it happens.

The Slow Fade: They found that the speed of the fluid and the temperature differences don't disappear instantly. Instead, they fade away slowly, like a sound echoing in a large hall that gets quieter and quieter.
The Formula: They created a specific formula (a decay rate) that tells you exactly how quiet the fluid will be after a certain amount of time.
The Catch: Because the fluid is "slippery" vertically (no vertical friction), it takes a bit longer to settle than it would in a normal fluid. However, the interaction between the movement and the temperature actually helps smooth things out a little bit faster than expected in some specific ways.

5. How They Solved It: The "Magic Mirror"

How did they figure this out? They didn't just guess. They used a mathematical technique called Fourier transforms.

The Analogy: Imagine the fluid's movement is a complex song made of many different notes (frequencies). The authors used a "magic mirror" (the Fourier transform) to break that complex song down into individual notes.
The Process: They looked at how each individual note behaves over time. They found that the "horizontal notes" die out quickly because of the friction, while the "vertical notes" behave in a specific, predictable pattern. By putting all these notes back together, they could reconstruct the whole picture and prove that the entire song eventually fades to silence.

Summary

In short, Yang and Zang proved that for a specific type of fluid moving on a slippery floor with heat involved:

Stability: If you start with a small disturbance, the system is safe and will not go out of control.
Decay: The disturbance will fade away over time, and they calculated the exact speed of that fading.

They did this by breaking the complex movement of the fluid into simple, manageable pieces (like musical notes) and showing that the "friction" in the horizontal direction is strong enough to eventually quiet the whole system down.

Technical Summary: Decay of Solutions to the Anisotropic Boussinesq Equations near Hydrostatic Balance in Half Space

Problem Statement
This paper investigates the initial-boundary value problem for the three-dimensional incompressible anisotropic Boussinesq system in the half-space $\mathbb{R}^3_+$ . The system models geophysical fluids and incorporates horizontal dissipation (viscosity and thermal diffusivity acting only in the horizontal directions, $\Delta_h = \partial_{x_1}^2 + \partial_{x_2}^2$ ) while lacking vertical dissipation. The study focuses on the stability and large-time behavior of perturbations near a specific hydrostatic equilibrium state defined by $u^{(0)} = (0,0,0)$ , $\Theta^{(0)} = x_3$ , and $P^{(0)} = \frac{1}{2}x_3^2$ . The system is subject to Navier slip boundary conditions for the velocity ( $u_3=0, \partial_3 u_1 = \partial_3 u_2 = 0$ ) and Dirichlet boundary conditions for the temperature ( $\Theta=0$ ) at the boundary $x_3=0$ .

Methodology
The authors employ a combination of energy methods, bootstrapping arguments, and integral representations via Fourier analysis to address the global stability and decay rates.

Nonlinear Stability (Energy Methods): To prove global existence and stability in the Sobolev space $H^3(\mathbb{R}^3_+)$ , the authors utilize a bootstrapping argument. They define energy functionals involving $H^3$ norms and horizontal dissipation terms. A critical component of this proof is handling the boundary terms arising from integration by parts. Specifically, the authors establish that under the given boundary conditions and initial data constraints (specifically $\partial_3^2 \theta_0 = 0$ at $x_3=0$ ), the third-order vertical derivative of the pressure vanishes at the boundary ( $\partial_3^3 p = 0$ ), which allows them to close the energy estimates without losing regularity. They rely on calculus inequalities in Sobolev spaces and anisotropic interpolation inequalities (similar to those in Lemma 2.2) to bound the nonlinear terms.
Decay Rates (Fourier Analysis and Integral Representation): To derive the decay rates for the nonlinear system, the authors move beyond direct energy estimates, which they note are insufficient for capturing optimal decay. Instead, they:
- Linearize the system around the hydrostatic balance.
- Apply a Fourier transform in the horizontal variables ( $x_h = (x_1, x_2)$ ) and Fourier cosine/sine transforms in the vertical variable ( $x_3$ ) to the linearized equations.
- Derive an explicit solution formula for the linearized problem involving seven kernel functions in the frequency space. These kernels exhibit degeneracy and anisotropy.
- Utilize Duhamel's principle to express the solution of the nonlinear system as a sum of the linear evolution and integral terms involving the nonlinearities.
- Estimate these integral terms by combining the decay properties of the linear kernels with the assumed decay rates of the solution (bootstrapping the decay).

Key Contributions and Results

Global Stability (Theorem 1.1): The paper establishes the global existence and uniqueness of solutions in $H^3(\mathbb{R}^3_+)$ for sufficiently small initial data. The solutions satisfy a uniform bound in $H^3$ and exhibit dissipation in the horizontal gradients ( $\nabla_h u, \nabla_h \theta$ ). This result confirms the stability of the hydrostatic balance under the specified anisotropic dissipation and boundary conditions.
Decay Rates (Theorem 1.2): The authors derive explicit algebraic decay rates for the global solution and its derivatives. Assuming the initial data possesses certain low-frequency regularity (measured by the operator $\Lambda_h^{-\sigma}$ with $9/10 < \sigma < 1$ ), the solution decays as follows for large time $t$ :
- The $L^2$ norms of velocity and temperature decay as $(1+t)^{-\sigma/2 + \delta}$ .
- The vertical derivatives ( $\partial_3 u, \partial_3 \theta$ ) decay as $(1+t)^{-\sigma/2 + 3\delta}$ .
- The horizontal gradients ( $\nabla_h u, \nabla_h \theta$ ) and the vertical derivative of the vertical velocity component ( $\partial_3 u_3$ ) decay as $(1+t)^{-\sigma + 1/2 + \delta}$ .
  Here, $\delta$ is a small positive parameter constrained by the regularity of the initial data.

Significance and Claims
The paper claims that while the decay rates for the linearized equations match those found in previous Cauchy problem studies (specifically referencing Ji, Yan, and Wu [10]), the large-time behavior for the nonlinear initial-boundary value problem is distinct. The authors emphasize that in the presence of boundaries and anisotropic dissipation, the decay rates depend heavily on the high regularity of the solutions ( $H^3$ ) and the specific coupling within the system.

The work highlights the challenges posed by anisotropic dissipation in bounded domains, where standard energy methods are insufficient to capture the precise decay rates. By constructing an integral representation with degenerate and anisotropic kernels, the authors demonstrate that the coupling between velocity and temperature provides extra smoothing, influencing the decay structure. The results extend the understanding of Boussinesq systems from whole-space Cauchy problems to half-space boundary value problems with partial dissipation, a setting relevant to geophysical fluid dynamics where vertical mixing is often negligible compared to horizontal mixing.