The goal of this work is to advance the characteristics of existing lattice Boltzmann Dirichlet velocity boundary schemes in terms of the accuracy, locality, stability, and mass conservation for arbitrarily grid-inclined straight walls, curved surfaces, and narrow fluid gaps, for both creeping and inertial flow regimes. We reach this objective with two infinite-member boundary classes: (1) the single-node "Linear Plus" (LI^{+}) and (2) the two-node "Extended Multireflection" (EMR). The LI^{+} unifies all directional rules relying on the linear combinations of up to three pre- or postcollision populations, including their "ghost-node" interpolations and adjustable nonequilibrium approximations.
View Article and Find Full Text PDFWe propose a procedure to implement Dirichlet velocity boundary conditions for complex shapes that use data from a single node only, in the context of the lattice Boltzmann method. Two ideas are at the base of this approach. The first is to generalize the geometrical description of boundary conditions combining bounce-back rule with interpolations.
View Article and Find Full Text PDFPhilos Trans A Math Phys Eng Sci
July 2020
This work addresses the Dirichlet boundary condition for momentum in the lattice Boltzmann method (LBM), with focus on the steady-state Stokes flow modelling inside non-trivial shaped ducts. For this task, we revisit a local and highly accurate boundary scheme, called the local second-order boundary (LSOB) method. This work reformulates the LSOB within the two-relaxation-time (TRT) framework, which achieves a more standardized and easy to use algorithm due to the pivotal parametrization TRT properties.
View Article and Find Full Text PDFImpact of the unphysical tangential advective-diffusion constraint of the bounce-back (BB) reflection on the impermeable solid surface is examined for the first four moments of concentration. Despite the number of recent improvements for the Neumann condition in the lattice Boltzmann method-advection-diffusion equation, the BB rule remains the only known local mass-conserving no-flux condition suitable for staircase porous geometry. We examine the closure relation of the BB rule in straight channel and cylindrical capillary analytically, and show that it excites the Knudsen-type boundary layers in the nonequilibrium solution for full-weight equilibrium stencil.
View Article and Find Full Text PDFThe effect of the heterogeneity in the soil structure or the nonuniformity of the velocity field on the modeled resident time distribution (RTD) and breakthrough curves is quantified by their moments. While the first moment provides the effective velocity, the second moment is related to the longitudinal dispersion coefficient (k_{T}) in the developed Taylor regime; the third and fourth moments are characterized by their normalized values skewness (Sk) and kurtosis (Ku), respectively. The purpose of this investigation is to examine the role of the truncation corrections of the numerical scheme in k_{T}, Sk, and Ku because of their interference with the second moment, in the form of the numerical dispersion, and in the higher-order moments, by their definition.
View Article and Find Full Text PDFPhys Rev E Stat Nonlin Soft Matter Phys
February 2015
This work focuses on the numerical solution of the Stokes-Brinkman equation for a voxel-type porous-media grid, resolved by one to eight spacings per permeability contrast of 1 to 10 orders in magnitude. It is first analytically demonstrated that the lattice Boltzmann method (LBM) and the linear-finite-element method (FEM) both suffer from the viscosity correction induced by the linear variation of the resistance with the velocity. This numerical artefact may lead to an apparent negative viscosity in low-permeable blocks, inducing spurious velocity oscillations.
View Article and Find Full Text PDFPhys Rev E Stat Nonlin Soft Matter Phys
June 2008
We show that a consistent modeling of porous flows needs at least one free collision relaxation rate to avoid a nonlinear dependency of the numerical errors on the viscosity. This condition is necessary to get the viscosity-independent permeability from the Stokes flow and to parametrize properly (with nondimensional physical numbers) the lattice Boltzmann Brinkman schemes. The two-relaxation-time (TRT) collision operator controls all coefficients of the higher-order corrections in steady solutions with a specific combination of its two collision rates, a possibility lacking for the Bhatnagar-Gross-Krook (BGK)-based single-relaxation-time schemes.
View Article and Find Full Text PDFPhilos Trans A Math Phys Eng Sci
March 2002
This article provides a concise exposition of the multiple-relaxation-time lattice Boltzmann equation, with examples of 15-velocity and 19-velocity models in three dimensions. Simulation of a diagonally lid-driven cavity flow in three dimensions at Re = 500 and 2000 is performed. The results clearly demonstrate the superior numerical stability of the multiple-relaxation-time lattice Boltzmann equation over the popular lattice Bhatnagar-Gross-Krook equation.
View Article and Find Full Text PDFPhilos Trans A Math Phys Eng Sci
March 2002
The filling process of viscoplastic metal alloys and plastics in expanding cavities is modelled using the lattice Boltzmann method in two and three dimensions. These models combine the regularized Bingham model for viscoplastic fluids with a free-interface algorithm. The latter is based on a modified immiscible lattice Boltzmann model in which one species is the fluid and the other one is considered to be a vacuum.
View Article and Find Full Text PDFPhys Rev E Stat Nonlin Soft Matter Phys
December 2003
We present a general framework for several previously introduced boundary conditions for lattice Boltzmann models, such as the bounce-back rule and the linear and quadratic interpolations. The objectives are twofold: first to give theoretical tools to study the existing link-type boundary conditions and their corresponding accuracy; second to design boundary conditions for general flows which are third-order kinetic accurate. Using these new boundary conditions, Couette and Poiseuille flows are exact solutions of the lattice Boltzmann models for a Reynolds number Re=0 (Stokes limit) for arbitrary inclination with the lattice directions.
View Article and Find Full Text PDF