A low-Mach methodology for efficient direct numerical simulations of variable property thermally driven flows

被引:8
|
作者
Demou, A. D. [1 ]
Frantzis, C. [1 ]
Grigoriadis, D. G. E. [1 ]
机构
[1] Univ Cyprus, Dept Mech & Mfg Engn, UCY CompSci, 1 Panepistimiou Ave, CY-2109 Nicosia, Cyprus
关键词
Low-Math approximation; Variable properties; Natural convection; Direct numerical simulation; Differentially heated cavity; NATURAL-CONVECTION FLOWS; LARGE TEMPERATURE DIFFERENCES; LARGE-EDDY SIMULATIONS; POISSONS EQUATION; INCOMPRESSIBLE FLOWS; BOUNDARY-CONDITIONS; BENCHMARK PROBLEM; NUMBER SOLVERS; ALGORITHM; VALIDITY;
D O I
10.1016/j.ijheatmasstransfer.2018.12.018
中图分类号
O414.1 [热力学];
学科分类号
摘要
Thermally driven flows can be approximated as constant property flows only when temperature differences are relatively small. In that case, the Oberbeck-Boussinesq approximation holds and the set of governing equations is simplified. For larger temperature differences, the variation of the fluid properties with temperature cannot be ignored and in the case of gasses the governing equations take the low Mach form. This significantly complicates the numerical solution of variable property against constant property flows because of the emergence of a variable coefficient Poisson equation for the pressure. In the present study, a numerical methodology for the direct numerical simulation of thermally driven low-Mach flows is presented. In this framework, a pressure-splitting scheme is utilised to transform the variable coefficient Poisson equation for the pressure into a constant coefficient Poisson equation, improving the efficiency of the numerical solution. The consistent boundary conditions for the pressure are derived and presented. Furthermore, the proposed methodology is validated for the natural convection of air inside a differentially heated cavity, for a wide range of temperature differences, both within and outside the limits of applicability of the Oberbeck-Boussinesq approximation. Finally, to demonstrate the potential of this methodology, the three-dimensional natural convection of air inside a differentially heated cavity is simulated for a Rayleigh number Ra = 2.0 x 10(9) and for temperature differences Delta T = 50 K, 100 K, and 200 K. It is found that, with increasing temperature difference, the symmetry around the centre of the cavity that characterises the Oberbeck-Boussinesq solution is lost. In addition, the laminar-turbulent transition point on the heated and cooled walls changes position, moving to a more upstream position on the heated wall and a more downstream position on the cooled wall. (C) 2018 Elsevier Ltd. All rights reserved.
引用
收藏
页码:539 / 549
页数:11
相关论文
共 29 条