Simulation of Low-Speed Buoyant Flows with a Stabilized Compressible/Incompressible Formulation: The Full Navier-Stokes Approach versus the Boussinesq Model

This paper compares two strategies to compute buoyancy-driven flows using stabilized methods. Both formulations are based on a unified approach for solving compressible and incompressible flows, which solves the continuity, momentum, and total energy equations in a coupled entropy-consistent way. The first approach introduces the variable density thermodynamics of the liquid or gas without any artificial buoyancy terms, i.e., without applying any approximate models into the Navier-Stokes equations. Furthermore, this formulation holds for flows driven by high temperature differences. Further advantages of this formulation are seen in the fact that it conserves the total energy and it lacks the incompressibility inconsistencies due to volume changes induced by temperature variations. The second strategy uses the Boussinesq approximation to account for temperature-driven forces. This method models the thermal terms in the momentum equation through a temperature-dependent nonlinear source term. Computer examples show that the thermodynamic approach, which does not introduce any artificial terms into the Navier-Stokes equations, is conceptually simpler and, with the incompressible stabilization matrix, attains similar residual convergence with iteration count to methods based on the Boussinesq approximation. For the Boussinesq model, the SUPG and SGS methods are compared, displaying very similar computational behavior. Finally, the VMS a posteriori error estimator is applied to adapt the mesh, helping to achieve better accuracy for the same number of degrees of freedom.