Large-Time Behavior of Solutions to Darcy–Boussinesq Equations with Non-Vanishing Scalar Acceleration Coefficient

We study the large-time behavior of solutions to Darcy–Boussinesq equations with a non-vanishing scalar acceleration coefficient, which model buoyancy-driven flows in porous media with spatially varying gravity. First, we show that the system admits steady-state solutions of the form ( u , ρ , p ) = ( 0 , ρ s , p s ) , where ρ s is characterised by the hydrostatic balance ∇ p s = − ρ s ∇ Ψ . Second, we prove that the steady-state solution satisfying ∇ ρ s = δ ( x , y ) ∇ Ψ is linearly stable provided that δ ( x , y ) < δ 0 < 0 , while the system exhibits Rayleigh–Taylor instability if Ψ = g y , ρ s = δ 0 g and δ 0 > 0 . Finally, despite the inherent Rayleigh–Taylor instability that may trigger exponential growth in time, we prove that for any sufficiently regular initial data, the solutions of the system asymptotically converge towards the vicinity of a steady-state solution, where the velocity field is zero, and the new state is determined by hydrostatic balance. This work advances porous media modeling for geophysical and engineering applications, emphasizing the critical interplay of gravity, inertia, and boundary conditions.