We investigate the instability and stability of specific steady-state solutions of the two-dimensional non-homogeneous, incompressible, and viscous Navier-Stokes equations under the influence of a general potential f. This potential is commonly used to model fluid motions in celestial bodies. First, we demonstrate that the system admits only steady-state solutions of the form (, V, p) = (₀, 0, P₀), where P₀ and ₀ satisfy the hydrostatic balance condition P₀=-₀ f. Additionally, the relationship between ₀ and the potential function f is constrained by the condition (ₘ₀, ₗ₀) (ₗf, ₘf) =0, which allows us to express ₀ as h (x, y) f. Second, when there exists a point (x₀, y₀) such that h (x₀, y₀) >0, we establish the linear instability of these solutions. Furthermore, we demonstrate their nonlinear instability in both the Lipschitz and Hadamard senses through detailed nonlinear energy estimates. This instability aligns with the well-known Rayleigh-Taylor instability. Our study signficantly extends and generalizes the existing mathematical results, which have predominantly focused on the scenarios involving a uniform gravitational field characterized by f= (0, g). Finally, we show that these steady states are linearly stable provided that h (x, y) <0 holds throughout the domain. Moreover, they exhibit nonlinear stability when h (x, y) is a negative constant.
Li et al. (Mon,) studied this question.