Two‐dimensional flow past circular cylinders using finite volume lattice Boltzmann formulation

In this article, an extension to the total variation diminishing finite volume formulation of the lattice Boltzmann equation method on unstructured meshes was presented. The quadratic least squares procedure is used for the estimation of first‐order and second‐order spatial gradients of the particle distribution functions. The distribution functions were extrapolated quadratically to the virtual upwind node. The time integration was performed using the fourth‐order Runge–Kutta procedure. A grid convergence study was performed in order to demonstrate the order of accuracy of the present scheme. The formulation was validated for the benchmark two‐dimensional, laminar, and unsteady flow past a single circular cylinder. These computations were then investigated for the low Mach number simulations. Further validation was performed for flow past two circular cylinders arranged in tandem and side‐by‐side. Results of these simulations were extensively compared with the previous numerical data. Copyright © 2011 John Wiley & Sons, Ltd.     

Time‐dependent finite element simulations of a shear‐thinning viscoelastic fluid with application to blood flow

A new finite element method is developed to simulate time‐dependent viscoelastic shear‐thinning flows characterized by the generalized Oldroyd‐B model. The focus of the algorithm is improved stability through a free‐energy dissipative scheme by using low‐order piecewise‐constant finite element approximations for stress. The algorithm is further modified by incorporating a pressure‐projection method, a DG‐upwinding scheme, a symmetric interior penalty DG method to solve the elliptic pressure‐update equation and a geometric multigrid preconditioner. The improved stability and cost to accuracy is compared when using higher order discontinuous bilinear approximation, where in addition, we consider the influence of a slope limiter for these elements. The algorithm is applied to the 2D start‐up‐driven cavity problem, and the stability of the free energy is illustrated and compared between element choices. An application of the model to modelling blood in small arterioles and channels is considered by simulating pulsatile blood flow through a stenotic arteriole. The individual influences of viscoelasticity and shear‐thinning within the generalized Oldroyd‐B model are investigated by comparing results to the Newtonian, generalized Newtonian and Oldroyd‐B models. Copyright © 2014 John Wiley & Sons, Ltd.     

Comparative study of lattice‐Boltzmann and finite volume methods for the simulation of laminar flow through a 4:1 planar contraction

In the present paper, a comparative study of numerical solutions for Newtonian fluids based on the lattice‐Boltzmann method (LBM) and the classical finite volume method (FVM) is presented for the laminar flow through a 4:1 planar contraction at a Reynolds number of value one, Re=1. In this study, the stress field for LBM is directly obtained from the distribution function. The calculations of the stress based on the FVM‐data use the evaluations of velocity gradients with finite differences. The stress field for both LBM and FVM is expressed in the present study in terms of the shear stress and the first normal stress difference. The lateral and axial profiles of the velocity, the shear stress and the first normal stress difference for both methods are investigated. It is shown that the LBM results for the velocity and the stresses are in excellent agreement with the FVM results. Copyright © 2004 John Wiley & Sons, Ltd.     

Lattice Boltzmann simulations for flow and heat/mass transfer problems in a three‐dimensional porous structure

The lattice Boltzmann method (LBM) for a binary miscible fluid mixture is applied to problems of transport phenomena in a three‐dimensional porous structure. Boundary conditions for the particle distribution function of a diffusing component are described in detail. Flow characteristics and concentration profiles of diffusing species at a pore scale in the structure are obtained at various Reynolds numbers. At high Reynolds numbers, the concentration profiles are highly affected by the flow convection and become completely different from those at low Reynolds numbers. The Sherwood numbers are calculated and compared in good agreement with available experimental data. The results indicate that the present method is useful for the investigation of transport phenomena in porous structures. Copyright © 2003 John Wiley & Sons, Ltd.     

Discontinuous Galerkin spectral element lattice Boltzmann method on triangular element

Discontinuous Galerkin spectral element method is used to solve the lattice Boltzmann equation (LBE) in the discrete velocity space. The triangular elements are adopted because of their flexibility to deal with complex geometries. The flow past a circular cylinder is simulated by the proposed scheme. The results are consistent with those obtained from the previous numerical methods and experiments. Copyright © 2003 John Wiley & Sons, Ltd.     

An application of the finite difference‐based lattice Boltzmann model to simulating flow‐induced noise

This paper presents a novel approach to simulate aerodynamically generated sounds by modifying the finite difference‐based lattice BGK compressible fluid model for the purpose of speeding up the calculation and also stabilizing the numerical scheme. With the model, aerodynamic sounds generated by a uniform flow around a two‐dimensional circular cylinder at Re = 150 are simulated. The third‐order‐accurate up‐wind scheme is used for the spatial derivatives, and the second‐order‐accurate Runge–Kutta method is applied for the time marching. The results show that we successively capture very small acoustic pressure fluctuations, with the same frequency of the Karman vortex street, much smaller than the whole pressure fluctuation around a circular cylinder. The propagation velocity of the acoustic waves shows that the points of peak pressure are biased upstream owing to the Doppler effect in the uniform flow. For the downstream, on the other hand, it is faster. It is also apparent that the amplitude of sound pressure is proportional to r^?1/2, r being the distance from the centre of the circular cylinder. Moreover, the edgetone generated by a two‐dimensional jet impinging on a wedge to predict the frequency characteristics of the discrete oscillations of a jet‐edge feedback cycle is investigated. The jet is chosen long enough to guarantee the parabolic velocity profile of the jet at the outlet, and the edge is of an angle of α = 23°. At a stand‐off distance w, the edge is inserted along the centreline of the jet, and a sinuous instability wave with real frequency is assumed to be created in the vicinity of the nozzle exit and to propagate towards the downstream. We have succeeded in capturing small pressure fluctuations resulting from periodic oscillation of jet around the edge. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical simulations of droplet formation in a cross‐junction microchannel by the lattice Boltzmann method

An immiscible liquid–liquid multiphase flow in a cross‐junction microchannel was numerically studied using the lattice Boltzmann method. An improved, immiscible lattice BGK model was proposed by introducing surface tension force based on the continuum surface force (CSF) method. Recoloring step was replaced by the anti‐diffusion scheme in the mixed region to reduce the side‐effect and control the thickness of the interface. The present method was tested by the simulation of a static bubble. Laplace's law and spurious velocities were examined. The results show that our model is more advantageous for simulations of immiscible fluids than the existing immiscible lattice BGK models. Computational results of multiphase flow in a cross‐junction microchannel were obtained and analyzed based on dimensionless numbers. It is found that the flow pattern is decided mostly by the capillary number at a small inlet flux. However, at the same capillary number, a large inlet flux will lead to much smaller droplet generation. For this case, the flow is determined by both the capillary number and the Weber number. Copyright © 2007 John Wiley & Sons, Ltd.     

Lattice‐Boltzmann and meshless point collocation solvers for fluid flow and conjugate heat transfer

The applicability and performance of the lattice‐Boltzmann (LB) and meshless point collocation methods as CFD solvers in flow and conjugate heat transfer processes are investigated in this work. Lid‐driven cavity flow and flow in a slit with an obstacle including heat transfer are considered as case studies. A comparison of the computational efficiency accuracy of the two methods with that of a finite volume method as implemented in a commercial package (ANSYS CFX, ANSYS Inc., Canonsburg, PA) is made. Utilizing the analogy between heat and mass transfer, an advection–diffusion LB model was adopted to simulate the heat transfer part of the slit flow problem followed by a rigorous mapping of the mass transfer variables to the heat transfer quantities of interest, thus circumventing the need for a thermal LB model. Direct comparison among the results of the three methods revealed excellent agreement over a wide range of Reynolds and Prandtl number values. Furthermore, an integrated computational scheme is proposed, utilizing the rapid convergence of the LB model in the flow part of the conjugate heat transfer problem with that of the meshless collocation method for the heat transfer part. The meshless treatment remains sufficiently rapid even for conduction‐controlled processes in contrast to the LB method, which is very rapid in the convection‐controlled case only. A single, common computational grid, composed of regularly distributed nodes is used, saving significant computational and coding time and ensuring convergence of the discrete Laplacian operator in the heat transfer part of the computations. Copyright © 2012 John Wiley & Sons, Ltd.     

Development of a lattice Boltzmann method for two‐layered shallow‐water flow

In this paper the dynamics of a two‐layered liquid, made of two immiscible shallow‐layers of different density, has been investigated within the framework of the lattice Boltzmann method (LBM). The LBM developed in this paper for the two‐layered, shallow‐water flow has been obtained considering two separate sets of LBM equations, one for each layer. The coupling terms between the two sets have been defined as external forces, acted on one layer by the other. Results obtained from the LBM developed in this paper are compared with numerical results obtained solving the two‐layered, shallow‐water equations, with experimental and other numerical results published in literature. The results are interesting. First, the numerical results obtained by the LBM and by the shallow‐water model can be considered as equivalent. Second, the LBM developed in this paper is able to simulate motion conditions on nonflat topography. Third, the agreement between the LBM (and also shallow‐water model) numerical results and the experimental results is good when the evolution of the flow does not depend on the viscosity, that is, during the initial phase of the flow, dominated by gravity and inertia forces. When the viscous forces dominate the evolution of the flow the agreement between numerical and experimental results depends strongly on the viscosity; it is good if the numerical LBM viscosity has the same order of magnitude of the liquid's kinematic viscosity. Copyright © 2012 John Wiley & Sons, Ltd.     

A hybrid FVM–LBM method for single and multi‐fluid compressible flow problems

The lattice Boltzmann method (LBM) has established itself as an alternative approach to solve the fluid flow equations. In this work we combine LBM with the conventional finite volume method (FVM), and propose a non‐iterative hybrid method for the simulation of compressible flows. LBM is used to calculate the inter‐cell face fluxes and FVM is used to calculate the node parameters. The hybrid method is benchmarked for several one‐dimensional and two‐dimensional test cases. The results obtained by the hybrid method show a steeper and more accurate shock profile as compared with the results obtained by the widely used Godunov scheme or by a representative flux vector splitting scheme. Additional features of the proposed scheme are that it can be implemented on a non‐uniform grid, study of multi‐fluid problems is possible, and it is easily extendable to multi‐dimensions. These features have been demonstrated in this work. The proposed method is therefore robust and can possibly be applied to a variety of compressible flow situations. Copyright © 2009 John Wiley & Sons, Ltd.     

Mobility‐dependent bifurcations in capillarity‐driven two‐phase fluid systems by using a lattice Boltzmann phase‐field model

Bifurcations in capillarity‐driven two‐phase fluid systems, due to different mobilities in phase‐field models for such systems, are studied by using a lattice Boltzmann method (LBM). Specifically, two‐dimensional (2D) and three‐dimensional (3D) droplets on a flat wall with given wettability variations are investigated. It is found that the mobility controls the rate of diffusive relaxation of the phase field from non‐equilibrium toward equilibrium, and similar to previous findings on mechanically driven two‐phase systems, the mobility is closely related to the contact line velocity. For the cases investigated, different mobilities across a critical value result in fundamentally different system evolution routes and final stable equilibrium states. These results may provide some implications for phase‐field study of droplet manipulations by surface wettability adjustments in microfluidics. Copyright © 2008 John Wiley & Sons, Ltd.     

Lattice Boltzmann modeling of pore‐scale fluid flow through idealized porous media

In order to capture the hydro‐mechanical impacts on the solid skeleton imposed by the fluid flowing through porous media at the pore‐scale, the flow in the pore space has to be modeled at a resolution finer than the pores, and the no‐slip condition needs to be enforced at the grain–fluid interface. In this paper, the lattice Boltzmann method (LBM), a mesoscopic Navier–Stokes solver, is shown to be an appropriate pore‐scale fluid flow model. The accuracy and lattice sensitivity of LBM as a fluid dynamics solver is demonstrated in the Poiseuille channel flow problem (2‐D) and duct flow problem (3‐D). Well‐studied problems of fluid creeping through idealized 2‐D and 3‐D porous media (J. Fluid Mech. 1959; 5 (2):317–328, J. Fluid Mech. 1982; 115 :13–26, Int. J. Multiphase Flow 1982; 8 (4):343–360, Phys. Fluids A 1989; 1 (1):38–46, Int. J. Numer. Anal. Meth. Geomech. 1999; 23 :881–904, Int. J. Numer. Anal. Meth. Geomech. 2010; DOI: 10.1002/nag.898, Int. J. Multiphase Flow 1982; 8 (3):193–206) are then simulated using LBM to measure the friction coefficient for various pore throats. The simulation results agree well with the data reported in the literature. The lattice sensitivity of the frictional coefficient is also investigated. Copyright © 2010 John Wiley & Sons, Ltd.     

Application of the lattice‐Boltzmann method to study flow and dispersion in channels with and without expansion and contraction geometry

The lattice‐Boltzmann (LB) method, derived from lattice gas automata, is a relatively new technique for studying transport problems. The LB method is investigated for its accuracy to study fluid dynamics and dispersion problems. Two problems of relevance to flow and dispersion in porous media are addressed: (i) Poiseuille flow between parallel plates (which is analogous to flow in pore throats in two‐dimensional porous networks), and (ii) flow through an expansion–contraction geometry (which is analogous to flow in pore bodies in two‐dimensional porous networks). The results obtained from the LB simulations are compared with analytical solutions when available, and with solutions obtained from a finite element code (FIDAP) when analytical results are not available. Excellent agreement is found between the LB results and the analytical/FIDAP solutions in most cases, indicating the utility of the lattice‐Boltzmann method for solving fluid dynamics and dispersion problems. Copyright © 1999 John Wiley & Sons, Ltd.     

A Hermite finite element method for incompressible fluid flow

We describe some Hermite stream function and velocity finite elements and a divergence‐free finite element method for the computation of incompressible flow. Divergence‐free velocity bases defined on (but not limited to) rectangles are presented, which produce pointwise divergence‐free flow fields (∇· u _h≡0). The discrete velocity satisfies a flow equation that does not involve pressure. The pressure can be recovered as a function of the velocity if needed. The method is formulated in primitive variables and applied to the stationary lid‐driven cavity and backward‐facing step test problems. Copyright © 2009 John Wiley & Sons, Ltd.     

A high‐order compact finite‐difference lattice Boltzmann method for simulation of steady and unsteady incompressible flows

A high‐order compact finite‐difference lattice Boltzmann method (CFDLBM) is proposed and applied to accurately compute steady and unsteady incompressible flows. Herein, the spatial derivatives in the lattice Boltzmann equation are discretized by using the fourth‐order compact FD scheme, and the temporal term is discretized with the fourth‐order Runge–Kutta scheme to provide an accurate and efficient incompressible flow solver. A high‐order spectral‐type low‐pass compact filter is used to stabilize the numerical solution. An iterative initialization procedure is presented and applied to generate consistent initial conditions for the simulation of unsteady flows. A sensitivity study is also conducted to evaluate the effects of grid size, filtering, and procedure of boundary conditions implementation on accuracy and convergence rate of the solution. The accuracy and efficiency of the proposed solution procedure based on the CFDLBM method are also examined by comparison with the classical LBM for different flow conditions. Two test cases considered herein for validating the results of the incompressible steady flows are a two‐dimensional (2‐D) backward‐facing step and a 2‐D cavity at different Reynolds numbers. Results of these steady solutions computed by the CFDLBM are thoroughly compared with those of a compact FD Navier–Stokes flow solver. Three other test cases, namely, a 2‐D Couette flow, the Taylor's vortex problem, and the doubly periodic shear layers, are simulated to investigate the accuracy of the proposed scheme in solving unsteady incompressible flows. Results obtained for these test cases are in good agreement with the analytical solutions and also with the available numerical and experimental results. The study shows that the present solution methodology is robust, efficient, and accurate for solving steady and unsteady incompressible flow problems even at high Reynolds numbers. Copyright © 2014 John Wiley & Sons, Ltd.     

A spectral‐element discontinuous Galerkin thermal lattice Boltzmann method for conjugate heat transfer applications

We present a spectral‐element discontinuous Galerkin thermal lattice Boltzmann method for fluid–solid conjugate heat transfer applications. Using the discrete Boltzmann equation, we propose a numerical scheme for conjugate heat transfer applications on unstructured, non‐uniform grids. We employ a double‐distribution thermal lattice Boltzmann model to resolve flows with variable Prandtl (Pr) number. Based upon its finite element heritage, the spectral‐element discontinuous Galerkin discretization provides an effective means to model and investigate thermal transport in applications with complex geometries. Our solutions are represented by the tensor product basis of the one‐dimensional Legendre–Lagrange interpolation polynomials. A high‐order discretization is employed on body‐conforming hexahedral elements with Gauss–Lobatto–Legendre quadrature nodes. Thermal and hydrodynamic bounce‐back boundary conditions are imposed via the numerical flux formulation that arises because of the discontinuous Galerkin approach. As a result, our scheme does not require tedious extrapolation at the boundaries, which may cause loss of mass conservation. We compare solutions of the proposed scheme with an analytical solution for a solid–solid conjugate heat transfer problem in a 2D annulus and illustrate the capture of temperature continuities across interfaces for conductivity ratio γ > 1. We also investigate the effect of Reynolds (Re) and Grashof (Gr) number on the conjugate heat transfer between a heat‐generating solid and a surrounding fluid. Steady‐state results are presented for Re = 5?40 and Gr = 10⁵?10⁶. In each case, we discuss the effect of Re and Gr on the heat flux (i.e. Nusselt number Nu) at the fluid–solid interface. Our results are validated against previous studies that employ finite‐difference and continuous spectral‐element methods to solve the Navier–Stokes equations. Copyright © 2016 John Wiley & Sons, Ltd.     

Efficient coupling of direct forcing immersed boundary-lattice Boltzmann method and finite element method to simulate fluid-structure interactions

Fluid-structure interactions (FSI) of rigid and flexible bodies are simulated in this article. For the fluid flow, multidirect forcing immersed boundary method (IBM) is adopted to capture the moving boundary, and lattice Boltzmann method (LBM) is used to evolve the flow field. Compared with our previous no-penetration IBM, less iterations are required in this work. In addition, larger velocity in lattice units can be used and the nonphysical force oscillations are suppressed due to the C³ 6-point kernel. Multi-relaxation-time collision operator and local grid refinement are also adopted in LBM to enhance the numerical stability and efficiency. To improve the efficiency of the FSI coupling algorithm, the mesh of the deformable structure can be coarser than the Lagrangian mesh using Newton-Cotes formulas to integrate the traction on the structure surface. A variety of benchmarks, including flow around a circular cylinder with Reynold number ranging from 20 to 200, forced oscillation of a circular cylinder, vortex-induced vibration (VIV) of an elastically mounted circular cylinder in two degrees of freedom, and VIV of an elastic cantilever beam attached to a circular cylinder, are carried out to evaluate the accuracy and stability of the present coupling algorithm. For the benchmark FSI problem considered in this article, a reduction of 54% of the calculation time is achieved using coarser structure mesh. As an application of the FSI coupling algorithm, the mechanism of an elastic beam in the wake of a circular cylinder is discussed.     

Application of the lattice Boltzmann method to flow in aneurysm with ring‐shaped stent obstacles

To resolve the characteristics of a highly complex flow, a lattice Boltzmann method with an extrapolation boundary technique was used in aneurysms with and without transverse objects on the upper wall, and results were compared with the non‐stented aneurysm. The extrapolation boundary concept allows the use of Cartesian grids even when the boundaries do not conform to Cartesian coordinates. To ease the code development and facilitate the incorporation of new physics, a new scientific programming strategy based on object‐oriented concepts was developed. The reduced flow, smaller vorticity magnitude and wall shear stress, and smaller du/dy near the dome of the aneurysm were observed when the proposed stent obstacles were used. The height of the stent obstacles was more effective to reduce the vorticity near the dome of the aneurysm than the width of the stent. The rectangular stent with 20% height‐of‐vessel radius was observed to be optimal and decreased the magnitude of the vorticity by 21% near the dome of the aneurysm. Copyright © 2008 John Wiley & Sons, Ltd.