Consistent lattice Boltzmann modeling of low-speed isothermal flows at finite Knudsen numbers in slip-flow regime. II. Application to curved boundaries.

Abstract

Gaseous flows inside microfluidic devices often fall in the slip-flow regime. According to this theoretical description, the Navier-Stokes model remains applicable in bulk, while at solid walls a slip velocity boundary model shall be considered. Physically, it is well established that, to properly account for the wall curvature, the wall slip velocity must be determined by the shear stress, rather than the normal component of the velocity derivative alone, as commonly applied to planar surfaces. It follows that the numerical transcription of this type of boundary condition is generally a challenging task for standard computational fluid dynamics (CFD) techniques. This paper aims to show that the realization of the slip velocity condition on arbitrarily shaped boundaries can be accomplished in a natural way with the lattice Boltzmann method (LBM). To substantiate this conclusion, this work undertakes the following three studies. First, we examine the conditions under which the generic reflection-type boundary rules used by LBM become consistent models for the slip velocity boundary condition. This effort makes use of the second-order Chapman-Enskog expansion method, where we address both planar and curved boundaries. The analysis also clarifies the capabilities and limitations behind the considered reflection-type slip schemes. Second, we revisit the family of parabolic accurate LBM slip boundary schemes, originally formulated in [Phys. Rev. E 96, 013311 (2017)2470-004510.1103/PhysRevE.96.013311] on the basis of the multireflection framework, and discuss their characteristics when operating on curved boundaries as well as the limitations of other less accurate LBM slip boundary formulations, such as the linearly accurate slip schemes and the widely popular "kinetic-based" boundary schemes. In addition, we also discuss the numerical stability of the parabolic slip schemes previously developed, providing an heuristic strategy to improve their stable range of operation. Third, we evaluate the performance of the several slip boundary schemes debated in this paper. The numerical tests correspond to two classical 2D benchmark flow problems of slip over non-planar solid surfaces, namely: (i) the velocity profile of the cylindrical Couette flow, and (ii) the permeability of a slow rarefied gas over a periodic array of circular cylindrical obstacles. The obtained numerical results confirm the competitiveness of the LBM when equipped with slip boundary schemes of parabolic accuracy as CFD tool to simulate slippage phenomena over arbitrarily non-planar surfaces. Indeed, although operating on a simple uniform mesh discretization, the LBM yields a similar, or even superior, level of accuracy compared to state-of-the-art FEM simulations conducted on hardworking body-fitted meshes. This conclusion establishes the LBM as a very appealing CFD technique for simulating microfluidic flows in the slip-flow regime, a result that deserves further exploration in future studies.