On the role of discrete mass conservation for non-Boussinesq flow simulations in enclosures. (January 2017)

• Record Type: Journal Article
• Authors: Kumar, Mukesh; Natarajan, Ganesh
• Is Part Of: International journal of heat and mass transfer. Volume 104(2017:Jan.)
• Page Start: 1283
• Page End: 1299
• Publication Date: 2017-01
• Subjects: Thermobuoyant flows -- Discrete consistency -- Discrete conservation -- Quasi-incompressible -- Non-Boussinesq convection -- Low Mach number; Heat -- Transmission -- Periodicals; Mass transfer -- Periodicals; Chaleur -- Transmission -- Périodiques; Transfert de masse -- Périodiques; Electronic journals; 621.4022
• DOI: 10.1016/j.ijheatmasstransfer.2016.09.073
• Languages: English
• ISSNs: 0017-9310

Abstract

Highlights: Comparison of two low Mach number algorithms for thermobuoyant flows. Both algorithms perform similarly in the Boussinesq limit. Discrete mass conservation errors lead to numerical instabilities for low Prandtl non-Boussinesq convection. Discrete mass conservation errors affect flow physics for unsteady non-Boussinesq convection. Abstract: We discuss the role of discrete mass conservation on the simulations of thermobuoyant flows in enclosures through a detailed numerical investigation in the Boussinesq and non-Boussinesq regimes. Employing a low Mach number formulation in conjunction with a pressure-based approach, we devise two unified algorithms that are capable of handling a wide range of temperature differences. These algorithms can be broadly categorized as discretely conservative and discretely consistent algorithms based on their implementation. While the discretely conservative algorithm computes the density by solving the mass conservation equation in each control volume (along with other conservation laws), the discretely consistent algorithm employs the algebraic equation of state to evaluate the local density. We show that the discretely consistent algorithm is prone to a discrete mass conservation error which can be detrimental to numerical stability in steady non-Boussinesq flows and can significantly influence the computed flowfield in case of unsteady non-Boussinesq convection. We also show that although the discretely conservative algorithm does not satisfy the discrete equation of state it introduces no significant errors in thermodynamic pressure for closed systems. Several numerical experiments over a range of temperature differences, Rayleigh and Prandtl numbers are carried out and the results whilst supporting the theoretical arguments also indicate that the discrete mass conservation errors are strongly linked respectively to numerical stability and solution accuracy in steady and unsteady non-Boussinesq flows.