A modified lattice Bhatnagar-Gross-Krook model for convection heat transfer in porous media

The lattice Bhatnagar-Gross-Krook (LBGK) model has become the most popular one in the lattice Boltzmann method for simulating the convection heat transfer in porous media. However, the LBGK model generally suffers from numerical instability at low fluid viscosities and effective thermal diffusivities. In this paper, a modified LBGK model is developed for incompressible thermal flows in porous media at the representative elementary volume scale, in which the shear rate and temperature gradient are incorporated into the equilibrium distribution functions. With two additional parameters, the relaxation times in the collision process can be fixed at a proper value invariable to the viscosity and the effective thermal diffusivity. In addition, by constructing a modified equilibrium distribution function and a source term in the evolution equation of temperature field, the present model can recover the macroscopic equations correctly through the Chapman-Enskog analysis, which is another key point different from previous LBGK models. Several benchmark problems are simulated to validate the present model with the proposed local computing scheme for the shear rate and temperature gradient, and the numerical results agree well with analytical solutions and/or those well-documented data in previous studies. It is also shown that the present model and the computational schemes for the gradient operators have a second-order accuracy in space, and better numerical stability of the present modified LBGK model than previous LBGK models is demonstrated.Comment: 38pages,50figure

Comparison of multiphase SPH and LBM approaches for the simulation of intermittent flows

Smoothed Particle Hydrodynamics (SPH) and Lattice Boltzmann Method (LBM) are increasingly popular and attractive methods that propose efficient multiphase formulations, each one with its own strengths and weaknesses. In this context, when it comes to study a given multi-fluid problem, it is helpful to rely on a quantitative comparison to decide which approach should be used and in which context. In particular, the simulation of intermittent two-phase flows in pipes such as slug flows is a complex problem involving moving and intersecting interfaces for which both SPH and LBM could be considered. It is a problem of interest in petroleum applications since the formation of slug flows that can occur in submarine pipelines connecting the wells to the production facility can cause undesired behaviors with hazardous consequences. In this work, we compare SPH and LBM multiphase formulations where surface tension effects are modeled respectively using the continuum surface force and the color gradient approaches on a collection of standard test cases, and on the simulation of intermittent flows in 2D. This paper aims to highlight the contributions and limitations of SPH and LBM when applied to these problems. First, we compare our implementations on static bubble problems with different density and viscosity ratios. Then, we focus on gravity driven simulations of slug flows in pipes for several Reynolds numbers. Finally, we conclude with simulations of slug flows with inlet/outlet boundary conditions. According to the results presented in this study, we confirm that the SPH approach is more robust and versatile whereas the LBM formulation is more accurate and faster

Large-scale grid-enabled lattice-Boltzmann simulations of complex fluid flow in porous media and under shear

Well designed lattice-Boltzmann codes exploit the essentially embarrassingly parallel features of the algorithm and so can be run with considerable efficiency on modern supercomputers. Such scalable codes permit us to simulate the behaviour of increasingly large quantities of complex condensed matter systems. In the present paper, we present some preliminary results on the large scale three-dimensional lattice-Boltzmann simulation of binary immiscible fluid flows through a porous medium derived from digitised x-ray microtomographic data of Bentheimer sandstone, and from the study of the same fluids under shear. Simulations on such scales can benefit considerably from the use of computational steering and we describe our implementation of steering within the lattice-Boltzmann code, called LB3D, making use of the RealityGrid steering library. Our large scale simulations benefit from the new concept of capability computing, designed to prioritise the execution of big jobs on major supercomputing resources. The advent of persistent computational grids promises to provide an optimal environment in which to deploy these mesoscale simulation methods, which can exploit the distributed nature of compute, visualisation and storage resources to reach scientific results rapidly; we discuss our work on the grid-enablement of lattice-Boltzmann methods in this context.Comment: 17 pages, 6 figures, accepted for publication in Phil.Trans.R.Soc.Lond.