A numerical methodology for efficient simulations of non-Oberbeck-Boussinesq flows

• Published in: International journal of heat and mass transfer 2018-10, Vol.125, p.1156-1168
• Main Authors: Demou, A.D., Frantzis, C., Grigoriadis, D.G.E.
• Format: Article
• Language: English
• Subjects: Boussinesq equations; Coefficients; Computational fluid dynamics; Computer simulation; Conservation equations; Density; Direct numerical simulation; Direct numerical simulations; Energy dissipation; Fluid flow; Free convection; Heat transfer; Incompressible flow; Material properties; Methodology; Natural convection; Navier-Stokes equations; Non-Oberbeck-Boussinesq; Poisson equation; Rayleigh-Bénard; Simulation; Temperature dependence; Temperature gradients; Thermophysical properties; Variable properties
• ISSN: 0017-9310; 1879-2189
• DOI: 10.1016/j.ijheatmasstransfer.2018.04.135

•An efficient methodology for non-Oberbeck-Boussinesq flows is proposed.•A pressure correction scheme allows the utilisation of Fast Direct Solvers.•The methodology is validated for a wide range of Prandtl and Rayleigh numbers.•The potential of the methodology is demonstrated with a series of 3D simulations.•Significant differences between 2D and 3D results are revealed and discussed. The Oberbeck-Boussinesq (OB) form of the Navier-Stokes equations can provide reliable solutions only for problems where the density differences are relatively small and all other material properties can be considered constant. For heat transfer problems at large temperature differences, the dependence of the material properties on temperature is the main source of non-Oberbeck-Boussinesq (NOB) effects. When a numerical solution is attempted for this category of problems, due to the density variations of the heated medium, a variable coefficient Poisson equation for the pressure emerges which is difficult to solve in a computationally efficient manner. In the present study, an efficient methodology to treat incompressible flows with variable properties outside the range of validity of the OB approximation is proposed. In the context of this methodology, the variable coefficient Poisson equation for the pressure is transformed into a constant coefficient Poisson equation using a pressure-correction scheme. In addition, all thermophysical properties are considered to be temperature dependent for all terms of the conservation equations. The proposed methodology is validated against results provided by previous studies on the natural convection flow of air, water and glycerol. Moreover, the potential of this methodology is demonstrated with the direct numerical simulation (DNS) of the NOB Rayleigh-Bénard convection of water inside a three-dimensional (3D) cavity. The comparison between two-dimensional (2D) and 3D results reveals significant differences, highlighting the need for efficient methodologies capable to accurately simulate NOB problems in 3D.