Large eddy simulation of rotating turbulent flows and heat transfer by the lattice Boltzmann method

• Published in: Physics of fluids (1994) 2018-01, Vol.30 (1)
• Main Authors: Liou, Tong-Miin, Wang, Chun-Sheng
• Format: Article
• Language: English
• Subjects: Aerodynamics; Centrifugal compressors; Channel flow; Computational fluid dynamics; Computer simulation; Direct numerical simulation; Distribution functions; Fluid dynamics; Fluid flow; Gas turbine engines; Heat flux; Heat transfer; Induction heating; Large eddy simulation; Mathematical analysis; Mathematical models; Physics; Prandtl number; Relaxation time; Reynolds number; Reynolds stress; Rotation; Scale models; Simulation; Turbulent flow; Velocity; Working fluids
• ISSN: 1070-6631; 1089-7666
• DOI: 10.1063/1.5005901

Due to its advantage in parallel efficiency and wall treatment over conventional Navier-Stokes equation-based methods, the lattice Boltzmann method (LBM) has emerged as an efficient tool in simulating turbulent heat and fluid flows. To properly simulate the rotating turbulent flow and heat transfer, which plays a pivotal role in tremendous engineering devices such as gas turbines, wind turbines, centrifugal compressors, and rotary machines, the lattice Boltzmann equations must be reformulated in a rotating coordinate. In this study, a single-rotating reference frame (SRF) formulation of the Boltzmann equations is newly proposed combined with a subgrid scale model for the large eddy simulation of rotating turbulent flows and heat transfer. The subgrid scale closure is modeled by a shear-improved Smagorinsky model. Since the strain rates are also locally determined by the non-equilibrium part of the distribution function, the calculation process is entirely local. The pressure-driven turbulent channel flow with spanwise rotation and heat transfer is used for validating the approach. The Reynolds number characterized by the friction velocity and channel half height is fixed at 194, whereas the rotation number in terms of the friction velocity and channel height ranges from 0 to 3.0. A working fluid of air is chosen, which corresponds to a Prandtl number of 0.71. Calculated results are demonstrated in terms of mean velocity, Reynolds stress, root mean square (RMS) velocity fluctuations, mean temperature, RMS temperature fluctuations, and turbulent heat flux. Good agreement is found between the present LBM predictions and previous direct numerical simulation data obtained by solving the conventional Navier-Stokes equations, which confirms the capability of the proposed SRF LBM and subgrid scale relaxation time formulation for the computation of rotating turbulent flows and heat transfer.