Consistent second-order boundary implementations for convection-diffusion lattice Boltzmann method

Novel Schemes of No-Slip Boundary Conditions for the Discrete Unified Gas Kinetic Scheme Based on the Moment Constraints

The boundary conditions are crucial for numerical methods. This study aims to contribute to this growing area of research by exploring boundary conditions for the discrete unified gas kinetic scheme (DUGKS). The importance and originality of this study are that it assesses and validates the novel schemes of the bounce back (BB), non-equilibrium bounce back (NEBB), and Moment-based boundary conditions for the DUGKS, which translate boundary conditions into constraints on the transformed distribution functions at a half time step based on the moment constraints. A theoretical assessment shows that both present NEBB and Moment-based schemes for the DUGKS can implement a no-slip condition at the wall boundary without slip error. The present schemes are validated by numerical simulations of Couette flow, Poiseuille flow, Lid-driven cavity flow, dipole-wall collision, and Rayleigh-Taylor instability. The present schemes of second-order accuracy are more accurate than the original schemes. Both present NEBB and Moment-based schemes are more accurate than the present BB scheme in most cases and have higher computational efficiency than the present BB scheme in the simulation of Couette flow at high Re. The present Moment-based scheme is more accurate than the present BB, NEBB schemes, and reference schemes in the simulation of Poiseuille flow and dipole-wall collision, compared to the analytical solution and reference data. Good agreement with reference data in the numerical simulation of Rayleigh-Taylor instability shows that they are also of use to the multiphase flow. The present Moment-based scheme is more competitive in boundary conditions for the DUGKS.

Unified directional parabolic-accurate lattice Boltzmann boundary schemes for grid-rotated narrow gaps and curved walls in creeping and inertial fluid flows

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. On this basis, we propose three groups of LI^{+} nonequilibrium local corrections: (1) the LI_{1}^{+} is parametrized, meaning that its steady-state solution is physically consistent: the momentum accuracy is viscosity-independent in Stokes flow, and it is fixed by the Reynolds number (Re) in inertial flow; (2) the LI_{3}^{+} is parametrized, exact for arbitrary grid-rotated Poiseuille force-driven Stokes flow and thus most accurate in porous flow; and (3) the LI_{4}^{+} is parametrized, exact for pressure and inertial term gradients, and hence advantageous in very narrow porous gaps and at higher Reynolds range. The directional, two-relaxation-time collision operator plays a crucial role for all these features, but also for efficiency and robustness of the boundary schemes due to a proposed nonequilibrium linear stability criterion which reliably delineates their suitable coefficients and relaxation space. Our methodology allows one to improve any directional rule for Stokes or Navier-Stokes accuracy, but their parametrization is not guaranteed. In this context, the parametrized two-node EMR class enlarges the single-node schemes to match exactness in a grid-rotated linear Couette flow modeled with an equilibrium distribution designed for the Navier-Stokes equation (NSE). However, exactness of a grid-rotated Poiseuille NSE flow requires us to perform (1) the modification of the standard NSE term for exact bulk solvability and (2) the EMR extension towards the third neighbor node. A unique relaxation and equilibrium exact configuration for grid-rotated Poiseuille NSE flow allows us to classify the Galilean invariance characteristics of the boundary schemes without any bulk interference; in turn, its truncated solution suggests how, when increasing the Reynolds number, to avoid a deterioration of the mass-leakage rate and momentum accuracy due to a specific Reynolds scaling of the kinetic relaxation collision rate. The optimal schemes and strategies for creeping and inertial regimes are then singled out through a series of numerical tests, such as grid-rotated channels and rotated Couette flow with wall-normal injection, cylindrical porous array, and Couette flow between concentric cylinders, also comparing them against circular-shape fitted FEM solutions.

Three-Dimensional Simulations of Anisotropic Slip Microflows Using the Discrete Unified Gas Kinetic Scheme

The specific objective of the present work study is to propose an anisotropic slip boundary condition for three-dimensional (3D) simulations with adjustable streamwise and spanwise slip length by the discrete unified gas kinetic scheme (DUGKS). The present boundary condition is proposed based on the assumption of nonlinear velocity profiles near the wall instead of linear velocity profiles in a unidirectional steady flow. Moreover, a 3D corner boundary condition is introduced to the DUGKS to reduce the singularities. Numerical tests validate the effectiveness of the present method, which is more accurate than the bounce-back and specular reflection slip boundary condition in the lattice Boltzmann method. It is of significance to study the lid-driven cavity flow due to its applications and its capability in exhibiting important phenomena. Then, the present work explores, for the first time, the effects of anisotropic slip on the two-sided orthogonal oscillating micro-lid-driven cavity flow by adopting the present method. This work will generate fresh insight into the effects of anisotropic slip on the 3D flow in a two-sided orthogonal oscillating micro-lid-driven cavity. Some findings are obtained: The oscillating velocity of the wall has a weaker influence on the normal velocity component than on the tangential velocity component. In most cases, large slip length has a more significant influence on velocity profiles than small slip length. Compared with pure slip in both top and bottom walls, anisotropic slip on the top wall has a greater influence on flow, increasing the 3D mixing of flow. In short, the influence of slip on the flow field depends not only on slip length but also on the relative direction of the wall motion and the slip velocity. The findings can help in better understanding the anisotropic slip effect on the unsteady microflow and the design of microdevices.

Lattice Boltzmann method for simulation of diffusion magnetic resonance imaging physics in multiphase tissue models

We report an implementation of the lattice Boltzmann method (LBM) to integrate the Bloch-Torrey equation, which describes the evolution of the transverse magnetization vector and the fate of the signal of diffusion magnetic resonance imaging (dMRI). Motivated by the need to interpret dMRI experiments in biological tissues, and to offset the small time-step limitation of classical LBM, a hybrid LBM scheme is introduced and implemented to solve the Bloch-Torrey equation. A membrane boundary condition is presented which is able to accurately represent the effects of thin curvilinear membranes typically found in biological tissues. As implemented, the hybrid LBM scheme accommodates piece-wise uniform transport, dMRI parameters, periodic and mirroring outer boundary conditions, and finite membrane permeabilities on non-boundary-conforming inner boundaries. By comparing with analytical solutions of limiting cases, we demonstrate that the hybrid LBM scheme is more accurate than the classical LBM scheme. The proposed explicit LBM scheme maintains second-order spatial accuracy, stability, and first-order temporal accuracy for a wide range of parameters. The parallel implementation of the hybrid LBM code in a multi-CPU computer system, as well as on GPUs, is straightforward and efficient. Along with offering certain advantages over finite element or Monte Carlo schemes, the proposed hybrid LBM constitutes a flexible scheme that can by easily adapted to model more complex interfacial conditions and physics in heterogeneous multiphase tissue models and to accommodate sophisticated dMRI sequences.

Analysis and assessment of the no-slip and slip boundary conditions for the discrete unified gas kinetic scheme

The discrete unified gas kinetic scheme (DUGKS) with a force term is a finite volume solver for the Boltzmann equation. Unlike the standard lattice Boltzmann method (LBM), DUGKS can be applied on nonuniform grids. For both the LBM and DUGKS, the boundary conditions need to be processed through the density distribution function. So researchers introduced the boundary conditions from the LBM frame into the DUGKS. However, the accuracy of these boundary conditions in the DUGKS has not been studied thoroughly. Through strict theoretical deduction, we find that the bounce-back (BB) scheme leads to a different dependence of the numerical error term in the DUGKS as compared to the LBM. The error term is influenced by the relaxation time and the body force. And it can be reduced by lowering the kinetic viscosity. Unlike the BB scheme, the nonequilibrium bounce-back scheme has the ability to implement real no-slip boundary condition. Furthermore, two slip boundary conditions incorporated with Navier's slip model are introduced from the LBM framework into the DUGKS. The tangential momentum change-based (TMAC) scheme can be used directly in the DUGKS because it generates no numerical error term in the DUGKS. For the combination of the bounce-back and specular reflection schemes (BSR), the relation between the slip length and the combination parameter should be modified in accordance with the numerical error term. Analysis shows that the TMAC scheme can simulate a wider range of slip length than the BSR scheme. Numerical simulations of the Couette flow and the Poiseuille flow confirm our theoretical analysis.

Mixed bounce-back boundary scheme of the general propagation lattice Boltzmann method for advection-diffusion equations

In this work, a mixed bounce-back boundary scheme of general propagation lattice Boltzmann (GPLB) model is proposed for isotropic advection-diffusion equations (ADEs) with Robin boundary condition, and a detailed asymptotic analysis is also conducted to show that the present boundary scheme for the straight walls has a second-order accuracy in space. In addition, several numerical examples, including the Helmholtz equation in a square domain, the diffusion equation with sinusoidal concentration gradient, one-dimensional transient ADE with Robin boundary and an ADE with a source term, are also considered. The results indicate that the numerical solutions agree well with the analytical ones, and the convergence rate is close to 2.0. Furthermore, through adjusting the two parameters in the GPLB model properly, the present boundary scheme can be more accurate than some existing lattice Boltzmann boundary schemes.