微尺度振荡Couette流的格子Boltzmann模拟

采用有效多松弛时间-格子Boltzmann方法(Effective MRT-LBM)数值模拟了微尺度条件下的振荡Couette和Poiseuille流动. 在微流动LBM中引入Knudsen边界层模型,对松弛时间进行修正. 模拟时平板或外力以正弦周期振动,Couette流中考虑了单平板振动、上下板同相振动这两类情况. 研究结果表明,修正后的MRT-LBM模型能有效用于这类非平衡的微尺度流动模拟;对于Couette流,随着Kn数的增大,壁面滑移效应变得越明显. St越大,板间速度剖面的非线性特性越剧烈;两板同相振荡时,若Kn,St均较小,板间流体受到平板拖动剪切的影响很小,板间速度几乎重叠在一起;在振荡Poiseuille流动中,St数增大到一定值时,相位滞后现象减弱;相对于Kn数,St数对振荡Couette 和Poiseuille流中不同位置处速度相位差的产生有较大影响. 关键词： 格子Boltzmann方法 有效MRT模型 Knudsen层 振荡流     

Slip, immiscibility, and boundary conditions at the liquid-liquid interface

The conventional boundary conditions at the interface between two flowing liquids include continuity of the tangential velocity. We have tested this assumption with molecular dynamics simulations of Couette and Poiseuille flows of two-layered liquid systems, with various molecular structures and interactions. When the total liquid density near the interface drops significantly compared to the bulk values, the tangential velocity varies very rapidly there, and would appear discontinuous at continuum resolution. The value of this apparent slip is given by a Navier boundary condition.     

Modified relaxation time Monte Carlo method for continuum-transition gas flows

Gas flows in the continuum-transition regime often occur in micro-electro-mechanical systems. The relaxation time Monte Carlo (RTMC) method was modified by using an ellipsoid statistical model and a multiple translational temperature model in the BGK model equation to simulate continuum-transition gas flows. The modified RTMC method uses a simplified form of the generalized relaxation time, which is related to the macro velocity and the local Knudsen number. The results for Couette flow and Poiseuille flow in microchannels predicted using the modified RTMC and the DSMC are in good agreement with the modified RTMC being much faster than the DSMC for continuum-transition gas flow simulations.     

Discrete Boltzmann Method with Maxwell-Type Boundary Condition for Slip Flow

The rarefied effect of gas flow in microchannel is significant and cannot be well described by traditional hydrodynamic models. It has been known that discrete Boltzmann model(DBM) has the potential to investigate flows in a relatively wider range of Knudsen number because of its intrinsic kinetic nature inherited from Boltzmann equation.It is crucial to have a proper kinetic boundary condition for DBM to capture the velocity slip and the flow characteristics in the Knudsen layer. In this paper, we present a DBM combined with Maxwell-type boundary condition model for slip flow. The tangential momentum accommodation coefficient is introduced to implement a gas-surface interaction model.Both the velocity slip and the Knudsen layer under various Knudsen numbers and accommodation coefficients can be well described. Two kinds of slip flows, including Couette flow and Poiseuille flow, are simulated to verify the model.To dynamically compare results from different models, the relation between the definition of Knudsen number in hard sphere model and that in BGK model is clarified.     

Boundary Slip and Surface Interaction： A Lattice Boltzmann Simulation

The factors affecting slip length in Couette geometry flows are analysed by means of a two-phase mesoscopic lattice Boltzmann model including non-ideal fluid-fluid and fluid-wall interactions. The main factors influencing the boundary slip are the strength of interactions between fluid-fluid and fluid-wall particles. Other factors, such as fluid viscosity, bulk pressure may also change the slip length. We find that boundary slip only occurs under a certain density （bulk pressure）. If the density is large enough, the slip length will tend to zero. In our simulations, a low density layer near the wall does not need to be postulated a priori but emerges naturally from the underlying non-ideal mesoscopic dynamics. It is the low density layer that induces the boundary slip. The results may be helpful to understand recent experimental observations on the slippage of micro flows.     

纳米通道滑移流动的分子动力学模拟研究

本文采用非平衡分子动力学方法对平板纳米通道滑移流动进行了非平衡分子动力学模拟,获得了不同壁面势能和不同温度时流体的速度分布及密度分布。研究结果表明滑移速度在很大程度上决定于流体温度和壁面吸引力作用强度的大小。由于不同壁面吸引力时流体的密度分布受温度的影响规律不同,使得不同壁面吸引力时流体的滑移速度受温度影响规律也不一致。而且,流体结构受壁面流速的影响要受到壁面势能的制约。     

Exact solutions for the flow of non-Newtonian fluid with fractional derivative in an annular pipe

Generally speaking, rheological properties of materials are specified by their so-called constitutive equations. The simplest constitutive equation for a fluid is a Newtonian one, on which the classical Navier-Stokes theory is based. The mechanical behavior of many fluids is well described by this theory. However, there are many rheologically compli- cated fluids such as polymer solutions, blood and heavy oils which are inadequately de- scribed by a Newtonian constitutive equation that does …     

Implementation of diffuse reflection boundary conditions in a thermal lattice Boltzmann model with flux limiters

We discuss three new implementation versions of diffuse reflection boundary conditions in a thermal lattice Boltzmann model. Their accuracy is investigated in the case of Couette flow by considering the slip regime. The best results are recovered with versions 2 and 3, which rely on outgoing fluxes to express the particle distribution functions in the ghost nodes outside the flow domain. Version 2 is found to be more economical since it involves no interpolation procedure. This version was thereafter used to investigate the temperature profile in Couette flow for various values of Prandtl number, as well as the capability of the thermal LB model to capture the Knudsen minimum in Poiseuille flow.     

Slip flow over structured surfaces with entrapped microbubbles

On hydrophobic surfaces, roughness may lead to a transition to a superhydrophobic state, where gas bubbles at the surface can have a strong impact on a detected slip. We present two-phase lattice Boltzmann simulations of a Couette flow over structured surfaces with attached gas bubbles. Even though the bubbles add slippery surfaces to the channel, they can cause negative slip to appear due to the increased roughness. The simulation method used allows the bubbles to deform due to viscous stresses. We find a decrease of the detected slip with increasing shear rate which is in contrast to some recent experimental results implicating that bubble deformation cannot account for these experiments. Possible applications of bubble surfaces in microfluidic devices are discussed.     

Simulation of micro- and nano-scale flows via the lattice Boltzmann method

We use the lattice Boltzmann method (LBM) for analysis of high and moderate Knudsen number phenomena. Simulation results are presented for microscale Couette and Poiseuille flows. The slip velocity, nonlinear pressure drop, and mass flow rate are compared with previous numerical results and/or experimental data. The Knudsen minimum is successfully predicted for the first time within the LBM framework. These results validate the usage of the LBM based commercial, arbitrary geometry code PowerFLOW for simulating nanoscale problems.     

Analytic Solutions of Linearized Lattice Boltzmann Equation for Simple Flows

A general procedure to obtain analytic solutions of the linearized lattice Boltzmann equation for simple flows is developed. As examples, the solutions for the Poiseuille and the plane Couette flows in two-dimensional space are obtained and studied in detail. The solutions not only have a component which is the solution of the Navier–Stokes equation, they also include a kinetic component which cannot be obtained by the Navier–Stokes equation. The kinetic component of the solutions is due to the finite-mean-free-path effect. Comparison between the analytic results and the numerical results of lattice-gas simulations is made, and they are found to be in accurate agreement.     

Cross-stream migration of flexible molecules in a nanochannel

Molecular dynamics simulations are used to examine the cross-stream chain migration phenomenon in dilute polymer solutions that are flowing in nanochannels. In particular, both uniform planar shear (Couette) and pressure driven (Poiseuille) flows of a dilute polymer solution are studied using a bead-spring representation of polymer chains and a coarse grained model for the solvent. Our results show that three mechanisms govern the migration of deformable molecules in a nanochannel: (1) chain-wall hydrodynamic interactions, (2) thermal diffusion, and (3) gradient in chain mobility. These results are discussed in the context of recent experimental, numerical and theoretical work.     

非等温非牛顿黏性流体流动问题的修正光滑粒子动力学方法模拟

为准确、有效地模拟非等温非牛顿黏性流体的流动问题,本文基于一种不含核导数计算的核梯度修正格式和不可压缩条件给出了一种改进光滑粒子动力学(SPH)离散格式,它较传统SPH离散格式具有较高精度和较好稳定性.同时,为准确地描述温度场的演化过程,建立了非牛顿黏性的SPH温度离散模型.通过对等温Poiseuille流、喷射流和非等温Couette流、4:1收缩流进行模拟,并与其他数值结果作对比,分别验证了改进SPH方法模拟非牛顿黏性流动问题的可靠性和提出的SPH温度离散模型求解非等温流动问题的有效性和准确性.随后,运用改进SPH方法结合SPH温度离散模型对环形腔和C形腔内非等温非牛顿黏性流体的充模过程进行了试探性模拟研究,分析了数值模拟的收敛性,讨论了不同位置处热流参数对温度和流动的影响.     

Non-Maxwell slippage induced by surface roughness for microscale gas flow: a molecular dynamics simulation

Rarefied gas flows in rough microchannels are investigated by non-equilibrium molecular dynamics simulations. The surface roughness is modelled by an array of triangular modules. The Maxwell slip model is found to break down due to the surface roughness for gas flows in microchannels with large surface roughness. Non-Maxwell slippage shows that the slip length is smaller than that predicted by the Maxwell model and is nonlinearly related to the mean free path. For larger surface roughness and smaller Knudsen number, the non-Maxwell effect becomes more pronounced. The boundary conditions, generally including velocity slip, no-slip and negative slip, depend not only on the Knudsen number but also on the surface roughness. Simulation results show that A/λ?≈?1 is a good criterion to validate the no-slip boundary condition and A/λ?>?0.3 can be a criterion to judge the occurrence of non-Maxwell slippage, where A is the surface roughness size and?λ?is the mean free path of gas molecules. The permeability enhanced by the surface roughness may be responsible for the roughness-induced non-Maxwell slippage.     

微通道疏水表面滑移的耗散粒子动力学研究

了解疏水表面的滑移规律对其在流动减阻方面的应用至关重要.利用耗散粒子动力学(dissipative particle dynamics, DPD)方法研究了微通道疏水表面的滑移现象.采用固定住的粒子并配合修正的向前反弹机制,构建了DPD固体壁面边界模型,利用该边界模型模拟了平板间的Couette流动.研究结果表明,通过调整壁面与流体间排斥作用强度,壁面能实现从无滑移到滑移的转变,壁面与流体间排斥作用越强,即疏水性越强,壁面滑移越明显,并且滑移长度与接触角之间存在近似的二次函数关系.无滑移时壁面附近密度分布均匀,有滑移时壁面附近存在低密度区域,低密度区域阻碍了动量传递,致使壁面产生滑移.     

Density waves and the effect of wall roughness in granular Poiseuille flow: Simulationand linear stability

The formation of density waves and the effect of wall roughness on them are studied using molecular dynamics simulations of gravity-driven granular Poiseuille flow. Three basic types of structures are found in moderately dense flows: a plug, a sinuous wave and a slug; a new varicose wave mode has been identified in dense flows with channels of large widths at moderate dissipations; only clump-like structures appear in dilute flows. The simulation results are contrasted with the predictions of a linear stability analysis of the kinetic-theory continuum equations for granular Poiseuille flow. The theoretical predictions on the form of density waves are in qualitative agreement with simulations in denser flows, however, there are discrepancies between simulation and theory in dilute flows.     

H theorem, regularization, and boundary conditions for linearized 13 moment equations

An H theorem for the linearized Grad 13 moment equations leads to regularizing constitutive equations for higher fluxes and to a complete set of boundary conditions. Solutions for Couette and Poiseuille flows show good agreement with direct simulation Monte Carlo calculations. The Knudsen minimum for the relative mass flow rate is reproduced.     

A Lattice Boltzmann Kinetic Model for Microflow and Heat Transfer

In this paper, we propose a lattice Boltzmann BGK model for simulation of micro flows with heat transfer based on kinetic theory and the thermal lattice Boltzmann method (He et al., J. Comp. Phys. 146:282, 1998). The relaxation times are redefined in terms of the Knudsen number and a diffuse scattering boundary condition (DSBC) is adopted to consider the velocity slip and temperature jump at wall boundaries. To check validity and potential of the present model in modelling the micro flows, two two-dimensional micro flows including thermal Couette flow and thermal developing channel flow are simulated and numerical results obtained compare well with previous studies of the direct simulation Monte Carlo (DSMC), molecular dynamics (MD) approaches and the Maxwell theoretical analysis     

Phénomènes de glissement à l'interface liquide–solide

The problem of liquid–solid slip is described here, in a simplified manner. Today, several experiments have shown that substantial slip appears when a non-wetting liquid flows along a surface which is smooth on an atomic scale. This phenomena is characterised by a length, called the slip length, or Navier length, generally denoted by L_s. A number of experiments indicate that this quantity may be as large as several hundreds of nanometers. Numerical simulations also show the existence of slip in non-wetting conditions, but the corresponding lengths found here are much smaller than those found experimentally. A theory, based on the existence of a gas film of nanometre thickness has been proposed, but has not yet been experimentally confirmed. Experiments on this are difficult, and sometimes controversial. To cite this article: P. Tabeling, C. R. Physique 5 (2004).     

Improving the Stability of the Multiple-Relaxation-Time Lattice Boltzmann Method by a Viscosity Counteracting Approach

Numerical instability may occur when simulating high Reynolds number flows by the lattice Boltzmann method (LBM). The multiple-relaxation-time (MRT) model of the LBM can improve the accuracy and stability, but is still subject to numerical instability when simulating flows with large single-grid Reynolds number (Reynolds number/grid number). The viscosity counteracting approach proposed recently is a method of enhancing the stability of the LBM. However, its effectiveness was only verified in the single-relaxation-time model of the LBM (SRT-LBM). This paper aims to propose the viscosity counteracting approach for the multiple-relaxation-time model (MRT-LBM) and analyze its numerical characteristics. The verification is conducted by simulating some benchmark cases: the two-dimensional (2D) lid-driven cavity flow, Poiseuille flow, Taylor-Green vortex flow and Couette flow, and three-dimensional (3D) rectangular jet. Qualitative and Quantitative comparisons show that the viscosity counteracting approach for the MRT-LBM has better accuracy and stability than that for the SRT-LBM.