Lattice Boltzmann Model for The Volume-Averaged Navier-Stokes Equations

A numerical method, based on the discrete lattice Boltzmann equation, is presented for solving the volume-averaged Navier-Stokes equations. With a modified equilibrium distribution and an additional forcing term, the volume-averaged Navier-Stokes equations can be recovered from the lattice Boltzmann equation in the limit of small Mach number by the Chapman-Enskog analysis and Taylor expansion. Due to its advantages such as explicit solver and inherent parallelism, the method appears to be more competitive with traditional numerical techniques. Numerical simulations show that the proposed model can accurately reproduce both the linear and nonlinear drag effects of porosity in the fluid flow through porous media.Comment: 9 pages, 2 figure

Microchannel fluid flow and heat transfer by lattice boltzmann method

This paper was presented at the 4th Micro and Nano Flows Conference (MNF2014), which was held at University College, London, UK. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute, ASME Press, LCN London Centre for Nanotechnology, UCL University College London, UCL Engineering, the International NanoScience Community, www.nanopaprika.eu.Micro flow has become a popular field of interest due to the advent of micro electromechanical systems (MEMS). In this work, the lattice Boltzmann method, a particle-based approach, is applied to simulate the two-dimensional micro channel fluid flow. We simulated fluid flow and heat transfer inside microchannel, the prototype application of this study is micro-heat exchangers. The main incentive to look at fluidic behaviour at micron scale is that micro devices tend to behave much differently from the objects we are used to handling in daily life. The choice of using LBM for micro flow simulation is a good one owing to the fact that it is based on the Boltzmann equation which is valid for the whole range of the Knudsen number. Slip velocity and temperature jump boundary conditions are used for the microchannel simulations with Knudsen number values covering the slip flow. The lattice Bhatnagar-Gross-Krook single relaxation time approximation was used. The results found are compared with the Navier-Stokes analytical and numerical results available in the literature and good matches are observed

A volume-preserving sharpening approach for the propagation of sharp phase boundaries in multiphase lattice Boltzmann simulations

Lattice Boltzmann models that recover a macroscopic description of multiphase flow of immiscible liquids typically represent the boundaries between phases using a scalar function, the phase field, that varies smoothly over several grid points. Attempts to tune the model parameters to minimise the thicknesses of these interfaces typically lead to the interfaces becoming fixed to the underlying grid instead of advecting with the fluid velocity. This phenomenon, known as lattice pinning, is strikingly similar to that associated with the numerical simulation of conservation laws coupled to stiff algebraic source terms. We present a lattice Boltzmann formulation of the model problem proposed by LeVeque and Yee [J. Comput. Phys. 86, 187] to study the latter phenomenon in the context of computational combustion, and offer a volume-conserving extension in multiple space dimensions. Inspired by the random projection method of Bao and Jin [J. Comput. Phys. 163, 216] we further generalise this formulation by introducing a uniformly distributed quasi-random variable into the term responsible for the sharpening of phase boundaries. This method is mass conserving and the statistical average of this method is shown to significantly delay the onset of pinning

Steady State Convergence Acceleration of the Generalized Lattice Boltzmann Equation with Forcing Term through Preconditioning

Several applications exist in which lattice Boltzmann methods (LBM) are used to compute stationary states of fluid motions, particularly those driven or modulated by external forces. Standard LBM, being explicit time-marching in nature, requires a long time to attain steady state convergence, particularly at low Mach numbers due to the disparity in characteristic speeds of propagation of different quantities. In this paper, we present a preconditioned generalized lattice Boltzmann equation (GLBE) with forcing term to accelerate steady state convergence to flows driven by external forces. The use of multiple relaxation times in the GLBE allows enhancement of the numerical stability. Particular focus is given in preconditioning external forces, which can be spatially and temporally dependent. In particular, correct forms of moment-projections of source/forcing terms are derived such that they recover preconditioned Navier-Stokes equations with non-uniform external forces. As an illustration, we solve an extended system with a preconditioned lattice kinetic equation for magnetic induction field at low magnetic Prandtl numbers, which imposes Lorentz forces on the flow of conducting fluids. Computational studies, particularly in three-dimensions, for canonical problems show that the number of time steps needed to reach steady state is reduced by orders of magnitude with preconditioning. In addition, the preconditioning approach resulted in significantly improved stability characteristics when compared with the corresponding single relaxation time formulation.Comment: 47 pages, 21 figures, for publication in Journal of Computational Physic