过渡区微尺度流动的有效黏性多松弛系数格子Boltzmann模拟

为提高采用二维九速离散速度模型的格子Boltzmann方法 (LBM)模拟微尺度流动中非线性现象的精度和效率,引入Dongari等提出的有效平均分子自由程对黏性进行修正(Dongari N,Zhang Y H,Reese J M2011 J.Fluids Eng.133 071101);并针对以往研究微尺度流动时采用边界处理格式含有离散误差的问题,采用多松弛系数格子Boltzmann方法结合二阶滑移边界条件,对微尺度Couette流动和周期性Poiseuille流动进行模拟,并将速度分布以及质量流量等模拟结果与直接模拟蒙特卡罗方法模拟数据、线性Boltzmann方程的数值解以及现有的LBM模型模拟结果进行对比.结果表明,相对于现有的LBM模型,引入新的修正函数所建立的有效黏性多松弛系数LBM模型有效提高了LBM模拟过渡区的微尺度流动中的非线性现象的能力.     

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.     

Gas flow through a slit into a vacuum in a wide range of rarefaction

Gas flow through a two-dimensional slit into a vacuum is investigated by a direct simulation Monte Carlo method. Results for the mass flow rate are obtained as a function of the rarefaction parameter, which is inversely proportional to the Knudsen number. The distributions of density, temperature and mass velocity, and streamlines are presented. In the free molecular flow regime and in the hydrodynamic limit, our results agree with theoretical asymptotes, and in the transition regime, they compare well with numerical simulations by other authors. The text was submitted by the author in English.     

Boundary scheme for lattice Boltzmann modeling of micro-scale gas flow in organic-rich pores considering surface diffusion

We propose a boundary scheme for addressing multi-mechanism flow in a porous medium in slip and early transition flow regimes, which is frequently encountered in shale gas reservoirs. Micro-gaseous flow in organic-rich shale involves a complex flow mechanism. A self-developed boundary scheme that combines the non-equilibrium extrapolation scheme and the combined diffusive reflection and bounce-back scheme(half-way DBB) to embed the Langmuir slip boundary into the single-relaxation-time lattice Boltzmann method(SRT-LBM) enables us to describe this process, namely, the coupling effect of micro-gaseous flow and surface diffusion in organic-rich nanoscale pores. The present LBM model comes with the careful consideration of the local Knudsen number, local pressure gradient, viscosity correction model, and regularization procedure to account for the rarefied gas flows in irregular pores. Its validity and accuracy are verified by several benchmarking cases, and the calculated results by this boundary scheme accord well with our analytical solutions.This boundary scheme shows a higher accuracy than the existing studies. Additionally, a subiteration strategy is presented to tackle the coupled micro-gaseous flow and surface diffusion, which necessitates the iteration process matching of these two mechanisms. The multi-mechanism flow in the self-developed irregular pores is also numerically investigated and analyzed over a wide range of parameters. The results indicate that the present model can effectively capture the coupling effect of micro-gaseous flow and surface diffusion in a tree-like porous medium.     

Knudsen diffusion in silicon nanochannels

Measurements on helium and argon gas flow through an array of parallel, linear channels of 12 nm diameter and 200 microm length in a single crystalline silicon membrane reveal a Knudsen diffusion type transport from 10(2) to 10(7) in Knudsen number Kn. The classic scaling prediction for the transport diffusion coefficient on temperature and mass of diffusing species, D(He) is proportional to square root T, is confirmed over a T range from 40 K to 300 K for He and for the ratio of D(He)/D(Ar) is proportional to square root (m(Ar)/m(He)). Deviations of the channels from a cylindrical form, resolved with electron microscopy down to subnanometer scales, quantitatively account for a reduced diffusivity as compared to Knudsen diffusion in ideal tubular channels. The membrane permeation experiments are described over 10 orders of magnitude in Kn, encompassing the transition flow regime, by the unified flow model of Beskok and Karniadakis.     

微通道内气体流动的格子Boltzmann方法研究

采用格子Boltzmann方法模拟了微通道在滑移区内不同Knudsen数下的微气体Poiseuille流,分析了微气体流动的速度分布以及流量与压降的关系,并给出了相对滑移长度和Poiseuille数随Knudsen数的变化特性。研究结果表明,微气体Poiseuille流的速度轮廓呈抛物线分布,但是边界速度大于0,出现速...     

Gas Flow in Microchannels – A Lattice Boltzmann Method Approach

Gas flow in microchannels can often encounter tangential slip motion at the solid surface even under creeping flow conditions. To simulate low speed gas flows with Knudsen numbers extending into the transition regime, alternative methods to both the Navier–Stokes and direct simulation Monte Carlo approaches are needed that balance computational efficiency and simulation accuracy. The lattice Boltzmann method offers an approach that is particularly suitable for mesoscopic simulation where details of the molecular motion are not required. In this paper, the lattice Boltzmann method has been applied to gas flows with finite Knudsen number and the tangential momentum accommodation coefficient has been implemented to describe the gas-surface interactions. For fully-developed channel flows, the results of the present method are in excellent agreement with the analytical slip-flow solution of the Navier–Stokes equations, which are valid for Knudsen numbers less than 0.1. The present paper demonstrates that the lattice Boltzmann approach is a promising alternative simulation tool for the design of microfluidic devices.     

基于UGKS的过渡区方柱绕流数值研究

本文基于统一气体动理学格式(Unified Gas Kinetic Scheme,UGKS),对微尺度过渡区气体绕流方柱开展数值模拟研究,计算分析了Kn数对气体流动传热过程的影响规律。研究发现,随着Kn数增加,方柱壁面气体速度滑移和温度跳跃增大,壁面上压力、剪切力和热流也相应增大,方柱壁面-气体之间的换热得到强化;方柱对气体流动的阻碍作用减小,方柱前滞止区影响范围相对增大,方柱温度对柱后区域气体温度影响相对减小。     

A Comparison and Unification of Ellipsoidal Statistical and Shakhov BGK Models

The Ellipsoidal Statistical model (ES-model) and the Shakhov model (S-model) were constructed to correct the Prandtl number of the original BGK model through the modification of stress and heat flux. With the introduction of a new parameter to combine the ES-model and S-model, a generalized kinetic model can be developed. This new model can give the correct Navier-Stokes equations in the continuum flow regime. Through the adjustment of the new parameter, it provides abundant dynamic effect beyond the ES-model and S-model. Changing the free parameter, the physical performance of the new model has been tested numerically. The unified gas kinetic scheme (UGKS) is employed for the study of the new model. In transition flow regime, many physical problems, i.e., the shock structure and micro-flows, have been studied using the generalized model. With a careful choice of the free parameter, good results can be achieved for most test cases. Due to the property of the Boltzmann collision integral, the new parameter in the generalized kinetic model cannot be fully determined. It depends on the specific problem. Generally speaking, the S-model predicts more accurate numerical solutions in most test cases presented in this paper than the ES-model, while ES-model performs better in the cases where the flow is mostly driven by temperature gradient, such as a channel flow with large boundary temperature variation at high Knudsen number.     

跨流域高超声速绕流环境Boltzmann模型方程统一算法研究

如何准确可靠地模拟从外层空间高稀薄流到近地面连续流的航天器高超声速绕流环境与复杂流动变化机理是流体物理的前沿基础科学问题. 基于对Boltzmann方程碰撞积分的物理分析与可计算建模, 确立了可描述自由分子流到连续流区各流域不同马赫数复杂流动输运现象统一的Boltzmann模型速度分布函数方程, 发展了适于高、低不同马赫数绕流问题的离散速度坐标法和直接求解分子速度分布函数演化更新的气体动理论数值格式, 建立了模拟复杂飞行器跨流域高超声速飞行热环境绕流问题的气体动理论统一算法. 对稀薄流到连续流不同Knudsen数0.002 ≤Kn_∞ ≤1.618、不同马赫数下可重复使用卫星体再入过程(110–70 km)中高超声速绕流问题进行算法验证分析, 计算结果与典型文献的Monte Carlo直接模拟值及相关理论分析符合得较好. 研究揭示了飞行器跨流域不同高度高超声速复杂流动机理、绕流现象与气动力/热变化规律, 提出了一个通过数值求解介观Boltzmann模型方程, 可靠模拟高稀薄自由分子流到连续流跨流域高超声速气动力/热绕流特性统一算法.     

致密多孔介质中气体渗流的格子Boltzmann模拟

为了能用格子Boltzmann方法来正确地刻画致密多孔介质中微尺度流动问题,对单通道模型进行推广,并将其应用于孔隙群里气体渗流的数值模拟.通过对部分具有代表性的多孔介质中的实际流动问题进行模拟,研究渗透率与平均压力和Knudsen数之间的相互关系.基于理论分析及相关文献中的试验结果,验证模拟结果,为用格子Boltzmann方法深入研究气体渗流问题奠定基础.     

Tracer dispersion in a percolation network with spatial correlations

We analyze the transport properties of a neutral tracer in a carrier fluid flowing through percolationlike porous media with spatial correlations. We model convection in the mass transport process using the velocity field obtained by the numerical solution of the Navier-Stokes and continuity equations in the pore space. We find that the resulting statistical properties of the tracer show a transition from a subdiffusion regime at low Peclet number to an enhanced diffusion regime at high Peclet number.     

A unified gas-kinetic scheme for continuum and rarefied flows

With discretized particle velocity space, a multiscale unified gas-kinetic scheme for entire Knudsen number flows is constructed based on the BGK model. The current scheme couples closely the update of macroscopic conservative variables with the update of microscopic gas distribution function within a time step. In comparison with many existing kinetic schemes for the Boltzmann equation, the current method has no difficulty to get accurate Navier–Stokes (NS) solutions in the continuum flow regime with a time step being much larger than the particle collision time. At the same time, the rarefied flow solution, even in the free molecule limit, can be captured accurately. The unified scheme is an extension of the gas-kinetic BGK-NS scheme from the continuum flow to the rarefied regime with the discretization of particle velocity space. The success of the method is due to the un-splitting treatment of the particle transport and collision in the evaluation of local solution of the gas distribution function. For these methods which use operator splitting technique to solve the transport and collision separately, it is usually required that the time step is less than the particle collision time. This constraint basically makes these methods useless in the continuum flow regime, especially in the high Reynolds number flow simulations. Theoretically, once the physical process of particle transport and collision is modeled statistically by the kinetic Boltzmann equation, the transport and collision become continuous operators in space and time, and their numerical discretization should be done consistently. Due to its multiscale nature of the unified scheme, in the update of macroscopic flow variables, the corresponding heat flux can be modified according to any realistic Prandtl number. Subsequently, this modification effects the equilibrium state in the next time level and the update of microscopic distribution function. Therefore, instead of modifying the collision term of the BGK model, such as ES-BGK and BGK–Shakhov, the unified scheme can achieve the same goal on the numerical level directly. Many numerical tests will be used to validate the unified method.     

Dispersion measurements using time-of-flight remote detection MRI

Remote detection nuclear magnetic resonance and magnetic resonance imaging can be used to study fluid flow and dispersion in a porous medium from a purely Eulerian point of view (i.e., in a laboratory frame of reference). Information about fluid displacement is obtained on a macroscopic scale in a long-time regime, while local velocity distributions are averaged out. It is shown how these experiments can be described using the common flow propagator formalism and how experimental data can be analyzed to obtain effective porosity, flow velocity inside the porous medium, fluid dispersion and flow tracing of fluid.     

饱和含混多孔介质中环绕水平圆柱的自然对流

本文运用达西流模型，对饱和含湿多孔介质中环绕水平等温圆柱发生的自然对流进行了数值模拟和分析。流函数和温度分布采用傅里叶级数的形式，以联立求解二维流动的控制方程组，傅里叶级数项的系数由一维常微分方程的解得到。文中同时给出局部和平均努谢尔数随修正瑞利数的变化关系以及修正瑞利数处在０．１～１０范围内时二维流场和温度场分布。     

随机多孔介质的蒙特卡罗重构与渗流特性模拟

根据多孔介质微观结构的分形尺度标度特征,采用蒙特卡罗方法分别重构随机多孔介质的微观颗粒和孔隙结构,并基于分形毛管束模型研究多尺度多孔介质的气体渗流特性,建立多孔介质微观结构和宏观渗流特性的定量关系。结果表明：分形蒙特卡罗重构的多孔介质微细结构接近真实介质结构,气体渗流特性的计算结果与格子玻尔兹曼模拟数据较为吻合; 多孔介质气体渗透率随着克努森数的增加而增大,孔隙分形维数对于气体渗流的微尺度效应具有显著影响,而迂曲度分形维数对于表观渗透率和固有渗透率的比值影响可以忽略。提出的分形蒙特卡罗方法具有收敛速度快且计算误差与维数无关的优点,有利于深入理解多尺度多孔介质的渗流机理。     

A hybrid method for hydrodynamic-kinetic flow – Part II – Coupling of hydrodynamic and kinetic models

In this work we present a non stationary domain decomposition algorithm for multiscale hydrodynamic-kinetic problems, in which the Knudsen number may span from equilibrium to highly rarefied regimes. Our approach is characterized by using the full Boltzmann equation for the kinetic regime, the Compressible Euler equations for equilibrium, with a buffer zone in which the BGK-ES equation is used to represent the transition between fully kinetic to equilibrium flows.In this fashion, the Boltzmann solver is used only when the collision integral is non-stiff, and the mean free path is of the same order as the mesh size needed to capture variations in macroscopic quantities. Thus, in principle, the same mesh size and time steps can be used in the whole computation. Moreover, the time step is limited only by convective terms.Since the Boltzmann solver is applied only in wholly kinetic regimes, we use the reduced noise DSMC scheme we have proposed in Part I of the present work. This ensures a smooth exchange of information across the different domains, with a natural way to construct interface numerical fluxes. Several tests comparing our hybrid scheme with full Boltzmann DSMC computations show the good agreement between the two solutions, on a wide range of Knudsen numbers.     

Instabilities of uniform filtration flows with phase transition

New mechanisms of instability are described for vertical flows with phase transition through horizontally extended two-dimensional regions of a porous medium. A plane surface of phase transition becomes unstable at an infinitely large wavenumber and at zero wavenumber. In the latter case, the unstable flow undergoes reversible subcritical bifurcations leading to the development of secondary flows (which may not be horizontally uniform). The evolution of subcritical modes near the instability threshold is governed by the Kolmogorov-Petrovskii-Piskunov equation. Two examples of flow through a porous medium are considered. One is the unstable flow across a water-bearing layer above a layer that carries a vapor-air mixture under isothermal conditions in the presence of capillary forces at the phase transition interface. The other is the vertical flow with phase transition in a high-temperature geothermal reservoir consisting of two high-permeability regions separated by a low-permeability stratum.     

Vortex dynamics of a trapezoidal bluff body placed inside a circular pipe

The influence of Reynolds number and blockage ratio on the vortex dynamics of a trapezoidal bluff body placed inside a circular pipe is studied experimentally and numerically. Low aspect ratio, high blockage ratio, curved end conditions (junction of pipe and bluff body), axisymmetric upstream flow with shear and turbulence are some of the intrinsic features of this class of bluff body flows which have been scarcely addressed in the literature. A large range (200:200,000) of Reynolds number (Re_D) is covered in this study, encompassing all the three pipe flow regimes (laminar, transition and turbulent). Four different flow regimes are defined based on the distinct features of Strouhal number (St)–Re_D relation: steady, laminar irregular, transition and turbulent. The wake in the steady regime is stationary with no oscillations in the shear layer. The laminar regime is termed as irregular owing to irregular vortex shedding. The vortex shedding in this regime is observed to be symmetric. The emergence of separation bubble downstream of the bluff body on either side is another interesting feature of this regime, which is further observed to be symmetric. Two pairs of mean streamwise vortices are noticed in the near-wake regime, which are termed as reverse dipole-type wake topology. Beyond the irregular laminar regime, the Strouhal number falls gradually and vortex shedding becomes more periodic. This regime is named transition and occurs close to the Reynolds number at which transition to turbulence takes place in a fully developed pipe. The turbulent regime is characterised by a nearly constant Strouhal number. Typical Karman-type vortex shedding is noticed in this regime. The convection velocity, wake width formation length and irrecoverable pressure loss are quantified to highlight the influence of blockage ratio. These results will be useful to develop basic understanding of vortex dynamics of confined bluff body flow for several practical applications.