Lattice Boltzmann method for oscillatory Stokes flow with applications to micro- and nanodevices

Quartz Crystal Microbalance Frequency Response to Discrete Adsorbates in Liquids

Quartz crystal microbalance with dissipation monitoring (QCM-D) has become a major tool enabling accurate investigation of the adsorption kinetics of nanometric objects such as DNA fragments, polypeptides, proteins, viruses, liposomes, polymer, and metal nanoparticles. However, in liquids, a quantitative analysis of the experimental results is often intricate because of the complex interplay of hydrodynamic and adhesion forces varying with the physicochemical properties of adsorbates and functionalized QCM-D sensors. In the present paper, we dissect the role of hydrodynamics for the analytically tractable case of stiff contact, whereas the adsorbed rigid particles oscillate with the resonator without rotation. Under the assumption of the low surface coverage, we theoretically study the excess shear force exerted on the resonator, which has two contributions: (i) the fluid-mediated force due to flow disturbance created by the particle and (ii) the force exerted on the particle by the fluid and transmitted to the sensor via contact. The theoretical analysis enables an accurate interpretation of the QCM-D impedance measurements. It is demonstrated inter alia that for particles of the size comparable with protein molecules, the hydrodynamic force dominates over the inertial force and that the apparent mass derived from QCM independently of the overtone is about 10 times the Sauerbrey (inertial) mass. The theoretical results show excellent agreement with the results of experiments and advanced numerical simulations for a wide range of particle sizes and oscillation frequencies.

Variational solution to the lattice Boltzmann method for Couette flow

Literature studies of the lattice Boltzmann method (LBM) demonstrate hydrodynamics beyond the continuum limit. This includes exact analytical solutions to the LBM, for the bulk velocity and shear stress of Couette flow under diffuse reflection at the walls through the solution of equivalent moment equations. We prove that the bulk velocity and shear stress of Couette flow with Maxwell-type boundary conditions at the walls, as specified by two-dimensional isothermal lattice Boltzmann models, are inherently linear in Mach number. Our finding enables a systematic variational approach to be formulated that exhibits superior computational efficiency than the previously reported moment method. Specifically, the number of partial differential equations (PDEs) in the variational method grows linearly with quadrature order while the number of moment method PDEs grows quadratically. The variational method directly yields a system of linear PDEs that provide exact analytical solutions to the LBM bulk velocity field and shear stress for Couette flow with Maxwell-type boundary conditions. It is anticipated that this variational approach will find utility in calculating analytical solutions for novel lattice Boltzmann quadrature schemes and other flows.

Origin of spurious oscillations in lattice Boltzmann simulations of oscillatory noncontinuum gas flows

Oscillatory noncontinuum gas flows at the micro and nanoscales are characterized by two dimensionless groups: a dimensionless molecular length scale, the Knudsen number Kn, and a dimensionless frequency θ, relating the oscillatory frequency to the molecular collision frequency. In a recent study [Shi et al., Phys. Rev. E 89, 033305 (2014)10.1103/PhysRevE.89.033305], the accuracy of the lattice Boltzmann (LB) method for simulating these flows at moderate-to-large Kn and θ was examined. In these cases, the LB method exhibits spurious numerical oscillations that cannot be removed through the use of discrete particle velocities drawn from higher-order Gauss-Hermite quadrature. Here, we identify the origin of these spurious effects and formulate a method to minimize their presence. This proposed method splits the linearized Boltzmann Bhatnagar-Gross-Krook (BGK) equation into two equations: (1) a homogeneous "gain-free equation" that can be solved directly, containing terms responsible for the spurious oscillations; and (2) an inhomogeneous "remainder equation" with homogeneous boundary conditions (i.e., stationary boundaries) that is solved using the conventional LB algorithm. This proposed "splitting method" is validated using published high-accuracy numerical solutions to the linearized Boltzmann BGK equation where excellent agreement is observed.

Separate-phase model and its lattice Boltzmann algorithm for liquid-vapor two-phase flows in porous media

In this article, we formulate a set of the separate-phase governing equations at the representative-elementary-volume scale and develop its double-distribution function lattice Boltzmann (LB) algorithm to describe liquid-vapor two-phase flows with or without phase change in porous media. Different from those previous studies, the mathematical description in this article involves the Darcy force, viscous force, and pressure gradient, and the resulting LB simulations can well describe two-phase flows and mass transfer throughout porous media under the compounding effects of these forces. The LB algorithm was validated by simulating single-phase flows in porous media. Its results are in good agreement with those available analytical solutions. We also applied it to model water flows through a semi-infinite porous region bounded by a heated solid wall, where liquid-vapor phase change takes place. The numerical simulations recover the previous results in the limit of the zero Darcy number. Significantly, it reveals much richer two-phase flow and mass transfer characteristics in porous media adjacent to solid walls. The separate-phase model and its lattice Boltzmann algorithm in this article are effective means to gain more profound and clearer understandings of complex two-phase transport processes in a porous system.

Performance assessment of bio-inspired systems: flow sensing MEMS hairs

Despite vigorous growth in biomimetic design, the performance of man-made devices relative to their natural templates is still seldom quantified, a procedure which would however significantly increase the rigour of the biomimetic approach. We applied the ubiquitous engineering concept of a figure of merit (FoM) to MEMS flow sensors inspired by cricket filiform hairs. A well known mechanical model of a hair is refined and tailored to this task. Five criteria of varying importance in the biological and engineering fields are computed: responsivity, power transfer, power efficiency, response time and detection threshold. We selected the metrics response time and detection threshold for building the FoM to capture the performance in a single number. Crickets outperform actual MEMS on all criteria for a large range of flow frequencies. Our approach enables us to propose several improvements for MEMS hair-sensor design.

Lattice Boltzmann method for linear oscillatory noncontinuum flows

Oscillatory gas flows are commonly generated by micro- and nanoelectromechanical systems. Due to their small size and high operating frequencies, these devices often produce noncontinuum gas flows. Theoretical analysis of such flows requires solution of the unsteady Boltzmann equation, which can present a formidable challenge. In this article, we explore the applicability of the lattice Boltzmann (LB) method to such linear oscillatory noncontinuum flows; this method is derived from the linearized Boltzmann Bhatnagar-Gross-Krook (BGK) equation. We formulate four linearized LB models in the frequency domain, based on Gaussian-Hermite quadratures of different algebraic precision (AP). The performance of each model is assessed by comparison to high-accuracy numerical solutions to the linearized Boltzmann-BGK equation for oscillatory Couette flow. The numerical results demonstrate that high even-order LB models provide superior performance over the greatest noncontinuum range. Our results also highlight intrinsic deficiencies in the current LB framework, which is incapable of capturing noncontinuum behavior at high oscillation frequencies, regardless of quadrature AP and the Knudsen number.

Transverse harmonic oscillations of laminae in viscous fluids: a lattice Boltzmann study

In this paper, we use the lattice Boltzmann method with the Bhatnagar-Gross-Krook linear collision operator to study the flow physics induced by a rigid lamina undergoing moderately large harmonic oscillations in a viscous fluid. We propose a refill procedure for the hydrodynamic quantities in the lattice sites that are in the vicinity of the oscillating lamina. The numerically estimated flow field is used to compute the complex hydrodynamic function that describes the added mass and hydrodynamic damping experienced by the lamina. Results of the numerical simulations are validated against theoretical predictions for small amplitude vibrations and experimental and numerical findings for moderately large oscillations.