Numerical simulation of laminar mixing flow by develop of compressible algorithm

Abstract

The use of the classical Boussinesq approximation is a straightforward strategy for taking into account the buoyancy effect in incompressible solvers. This strategy is highly effective if density variation is low. Whenever the density variation is high, this can cause considerable deviation from the correct prediction of fluid flow behavior and the accurate estimation of heat transfer rate. In this study, an incompressible algorithm is suitably extended to solve high-density-variation fields caused by strong natural-convection with mixing of oxygen and nitrogen in unsteady laminar compressible flow in a cavity. The continuity, momentum, energy and species equations are discritized based on finite volume methods and then numerically solved with extended algorithm with SIMPEL method. This new algorithm is capable of solving both Boussinesq and non-Boussinesq regimes. The fluid is assumed to be calorically an ideal gas and its thermodynamic properties depend on temperature and pressure. The extended algorithm is then verified by solving the benchmark convecting cavity problem at Rayleigh 106 and a temperature range of ε = 0.01–0.6. The results show that the method can vigorously solve unsteady mixing flow fields with extreme density variation.

Keywords


[1] Lessani B, and Papalexandris MV (2006) Time-accurate calculation of variable density flows with strong temperature gradients and combustion. J Comput. Phys 212: 218–246.
[2] Venkateswaran S, Merkle L (1999) Analysis of preconditioning methods for the Euler and Navier–Stokes equations. Von Karman Lecture Series: 1999–03.
[3] Chenoweth DR and Paolucci S (1986) Natural convection in an enclosed vertical air layer with large horizontal temperature differences. J Fluid Mech: 173–210.
[4] Darbandi M and Hosseinizadeh SF (2006) Numerical simulation of thermobuoyant flow with large temperature variation. J Thermophys Heat Transfer: 285–296.
[5] Darbandi M and Hosseinizadeh SF (2007) Numerical study of natural convection in vertical enclosures using a novel non-Boussinesq algorithm. Numer Heat Tr, Part A 52: 849–873.
[6] Ibrahim A and Lemonnier D (2009) Numerical study of coupled double-diffusive natural convection and radiation in a square cavity filled with a N2-CO2 mixture. Int Commun Heat Mass 36: 197–202.
[7] Ebrahimi Kebria HR, Darbandi M, Hosseinizadeh SF (2011) Numerical simulation of low Mach number laminar reacting and non-reacting flow using a novel dual purpose algorithm. Numerical Heat Transfer, Part A, 59: 595–514.
[8] Versteeg HK and Malalasekera W (1995) An introduction to computational fluid dynamics, the finite volume method. Longman, Malaysia.
[9] Patankar, SV, Spalding DB (1972) A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows. Int J Heat Mass Transfer: 1787–1806.
[10] Paillere, H, Le Quere P, Weisman C, Vierendeels J, Dick E, Braack M, Dabbene F, Beccantini A, Studer E, Kloczko T, Corre C, Darbandi M, Hosseinizadeh SF (2005) Modeling of natural convection flow with large temperature differences: A benchmark problem for low Mach number solvers. Part 2. Contributions to the June 2004 conference, ESAIM: Mathematical Modeling and Numerical Analysis: 617–621.
[11] De Vahl Davis G (1983) Natural convection of air in a square cavity: A bench mark numerical solution. Int J for Nume. Methods in Fluids: 249–264.