Influence of mesh non-orthogonality on numerical simulation of buoyant jet flows

Abstract In the present research, we discuss the influence of mesh non-orthogonality on numerical solution of a type of buoyant flow. Buoyant jet flows are simulated numerically with hexahedral and prismatic mesh elements in an open source Computational Fluid Dynamics (CFD) code called “OpenFOAM”. Buoyant jet instability obtained with the prismatic meshes may be overestimated compared to that obtained with the hexahedral meshes when non-orthogonal correction is not applied in the code. Although the non-orthogonal correction method can improve the instability generated by mesh non-orthogonality, it may increase computation time required to reach a convergent solution. Thus, we propose modified solvers that can reduce the influence of mesh non-orthogonality and reduce the computation time compared to the existing solvers in OpenFOAM. It is demonstrated that calculations for a buoyant jet with a large temperature difference are performed faster by the modified solver.

[1]  G. D. Davis.  Natural convection of air in a square cavity: A bench mark numerical solution , 1983 .

[2]  Weeratunge Malalasekera,et al.  An introduction to computational fluid dynamics - the finite volume method , 2007 .

[3]  Christophe Eric Corre,et al.  Numerical simulations of a transient injection flow at low Mach number regime , 2008 .

[4]  S. Muzaferija,et al.  Adaptive finite volume method for flow prediction using unstructured meshes and multigrid approach , 1994 .

[5]  M. Darwish,et al.  The Finite Volume Method in Computational Fluid Dynamics: An Advanced Introduction with OpenFOAM® and Matlab , 2015 .

[6]  J. Remacle,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[7]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[8]  Hrvoje Jasak,et al.  Error analysis and estimation for the finite volume method with applications to fluid flows , 1996 .

[9]  Kai H. Luo,et al.  Large-eddy simulation of a turbulent forced plume , 2001 .

[10]  Jean-Marc Hérard,et al.  A finite volume method to solve the Navier-Stokes equations for incompressible flows on unstructured meshes , 2000 .

[11]  Anupam Dewan,et al.  URANS computations with buoyancy corrected turbulence models for turbulent thermal plume , 2014 .

[12]  Simon J. Rees,et al.  Simulation of steady-state natural convection using CFD , 2005 .

[13]  V. Babu,et al.  Simulation of high Rayleigh number natural convection in a square cavity using the lattice Boltzmann method , 2006 .

[14]  Taisuke Yonomoto,et al.  RANS analyses on erosion behavior of density stratification consisted of helium–air mixture gas by a low momentum vertical buoyant jet in the PANDA test facility, the third international benchmark exercise (IBE-3) , 2015 .

[15]  Michele Andreani,et al.  Synthesis of the OECD/NEA-PSI CFD benchmark exercise , 2016 .

[16]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[17]  Zhenghua Yan,et al.  A two-equation turbulence model and its application to a buoyant diffusion flame , 1999 .

[18]  Soonil Nam,et al.  Numerical simulation of thermal plumes , 1993 .

[19]  R. I. Issa,et al.  AN IMPROVED PISO ALGORITHM FOR THE COMPUTATION OF BUOYANCY-DRIVEN FLOWS , 2001 .