Lattice Boltzmann method for direct numerical simulation of turbulent flows

We present three-dimensional numerical simulations (DNS) of the Kida vortex flow, a prototypical turbulent flow, using a novel high-order lattice Boltzmann model. Extensive comparisons of various global and local statistical quantities obtained with an incompressible flow spectral element solver are reported. It is demonstrated that the lattice Boltzmann method is a promising alternative for DNS as it quantitatively capturesall the computed statistics of fluid turbulence

Lattice Boltzmann method for direct numerical simulation of turbulent flows

We present three-dimensional direct numerical simulations (DNS) of the Kida vortex flow, a prototypical turbulent flow, using a novel high-order lattice Boltzmann (LB) model. Extensive comparisons of various global and local statistical quantities obtained with an incompressible-flow spectral element solver are reported. It is demonstrated that the LB method is a promising alternative for DNS as it quantitatively captures all the computed statistics of fluid turbulenc

The development of thermal lattice Boltzmann models in incompressible limit

In this paper, an incompressible two-dimensional (2-D) and three-dimensional (3-D) thermohydrodynamics for the lattice Boltzmann scheme are developed. The basic idea is to solve the velocity field and the temperature field using two different distribution functions. A derivation of the lattice Boltzmann scheme from the continuous Boltzmann equation for 2-D is discussed in detail. By using the same procedure as in the derivation of the discretised density distribution function, we found that new lattice of four-velocity (2-D) and eight-velocity(3-D) models for internal energy density distribution function can be developed where the viscous and compressive heating effects are negligible. These models are validated by the numerical simulation of the 2-D porous plate Couette flow problem where the analytical solution exists and the natural convection flows in a cubic cavity

Link-wise Artificial Compressibility Method

The Artificial Compressibility Method (ACM) for the incompressible Navier-Stokes equations is (link-wise) reformulated (referred to as LW-ACM) by a finite set of discrete directions (links) on a regular Cartesian mesh, in analogy with the Lattice Boltzmann Method (LBM). The main advantage is the possibility of exploiting well established technologies originally developed for LBM and classical computational fluid dynamics, with special emphasis on finite differences (at least in the present paper), at the cost of minor changes. For instance, wall boundaries not aligned with the background Cartesian mesh can be taken into account by tracing the intersections of each link with the wall (analogously to LBM technology). LW-ACM requires no high-order moments beyond hydrodynamics (often referred to as ghost moments) and no kinetic expansion. Like finite difference schemes, only standard Taylor expansion is needed for analyzing consistency. Preliminary efforts towards optimal implementations have shown that LW-ACM is capable of similar computational speed as optimized (BGK-) LBM. In addition, the memory demand is significantly smaller than (BGK-) LBM. Importantly, with an efficient implementation, this algorithm may be one of the few which is compute-bound and not memory-bound. Two- and three-dimensional benchmarks are investigated, and an extensive comparative study between the present approach and state of the art methods from the literature is carried out. Numerical evidences suggest that LW-ACM represents an excellent alternative in terms of simplicity, stability and accuracy.Comment: 62 pages, 20 figure

Simulation of incompressible viscous flows around moving objects by a variant of immersed boundary-Lattice Boltzmann method

A variant of immersed boundary-lattice Boltzmann method (IB-LBM) is presented in this paper to simulate incompressible viscous flows around moving objects. As compared with the conventional IB-LBM where the force density is computed explicitly by Hook's law or the direct forcing method and the non-slip condition is only approximately satisfied, in the present work, the force density term is considered as the velocity correction which is determined by enforcing the non-slip condition at the boundary. The lift and drag forces on the moving object can be easily calculated via the velocity correction on the boundary points. The capability of the present method for moving objects is well demonstrated through its application to simulate flows around a moving circular cylinder, a rotationally oscillating cylinder, and an elliptic flapping wing. Furthermore, the simulation of flows around a flapping flexible airfoil is carried out to exhibit the ability of the present method for implementing the elastic boundary condition. It was found that under certain conditions, the flapping flexible airfoil can generate larger propulsive force than the flapping rigid airfoil

Fluid-Structure Interaction Simulation of a Coriolis Mass Flowmeter using a Lattice Boltzmann Method

In this paper we use a fluid-structure interaction (FSI) approach to simulate a Coriolis mass flowmeter (CMF). The fluid dynamics are calculated by the open source framework OpenLB, based on the lattice Boltzmann method (LBM). For the structural dynamics we employ the open source software Elmer, an implementation of the finite element method (FEM). A staggered coupling approach between the two software packages is presented. The finite element mesh is created by the mesh generator Gmsh to ensure a complete open source workflow. The Eigenmodes of the CMF, which are calculated by modal analysis are compared with measurement data. Using the estimated excitation frequency, a fully coupled, partitioned, FSI simulation is applied to simulate the phase shift of the investigated CMF design. The calculated phaseshift values are in good agreement to the measurement data and verify the suitability of the model to numerically describe the working principle of a CMF