On the global stability and large time behavior of solutions of the Boussinesq equations

We study the two-dimensional viscous Boussinesq equations, which model the motion of stratified flows in a circular domain influenced by a general gravitational potential f. First, we demonstrate that the Boussinesq equations admit steady-state solutions only in the form of hydrostatic equilibria, given by (u, , p) = (0, ₛ, pₛ), where the pressure gradient satisfies pₛ=-ₛ f. Subsequently, we establish that any hydrostatic equilibrium satisfying the condition ₛ= (x, y) f is linearly unstable if (x, y) is positive at some point (x, y) = (x₀, y₀), This instability corresponds to the well-known Rayleigh-Taylor instability. Thirdly, by employing a series of regularity estimates, we reveal that although the presence of the Rayleigh-Taylor instability increases the velocity, the system ultimately converges to a state of hydrostatic equilibrium. This result is achieved by analyzing perturbations around any state of hydrostatic equilibrium, including both stable and unstable configurations. Specifically, the state of hydrostatic equilibrium can be expressed as =- f+, where and are positive constants, provided that the global perturbation satisfies additional conditions. This highlights the system's tendency to stabilize into a hydrostatic state despite the presence of instabilities.