过渡区微尺度流动的有效黏性多松弛系数格子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模拟过渡区的微尺度流动中的非线性现象的能力.     

Effect of viscosity on stability and accuracy of the two-component lattice Boltzmann method with a multiple-relaxation-time collision operator investigated by the acoustic attenuation model

A two-component lattice Boltzmann method (LBM) with a multiple-relaxation-time (MRT) collision operator is presented to improve the numerical stability of the single relaxation time (SRT) model. The macroscopic and the momentum conservation equations can be retrieved through the Chapman—Enskog (C-E) expansion analysis. The equilibrium moment with the diffusion term is calculated, a diffusion phenomenon is simulated by utilizing the developed model, and the numerical stability is verified. Furthermore, the binary mixture channel model is designed to simulate the sound attenuation phenomenon, and the obtained simulation results are found to be consistent with the analytical solutions. The sound attenuation model is used to study the numerical stability and calculation accuracy of the LBM model. The simulation results show the stability and accuracy of the MRT model and the SRT model under different viscosity conditions. Finally, we study the influence of the error between the macroscopic equation of the MRT model and the standard incompressible Navier—Stokes equation on the calculation accuracy of the model to demonstrate the general applicability of the conclusions drawn by the sound attenuation model in the present study.     

多孔介质内溶解与沉淀过程的格子Boltzmann方法模拟

本文采用格子Boltzmann方法模拟了多孔介质内的溶解和沉淀现象, 并分析了雷诺数、施密特数、达姆科勒数对多孔介质孔隙结构及浓度分布的影响. 结果表明: 对于多孔介质内的溶解(沉淀)过程, 当雷诺数越大时, 孔隙率越大(小), 平均浓度值越小(大); 当达姆科勒数或施密特数较小时, 溶解和沉淀过程均受反应控制, 此时反应在多孔介质的固体表面较为均匀的发生; 当达姆科勒数或施密特数较大时, 溶解和沉淀过程均受扩散控制, 此时反应主要发生在上游及大孔隙区域.     

Lattice Boltzmann Simulation of the Cross Flow Over a Cantilevered and Longitudinally Vibrating Circular Cylinder

The multiple-relaxation-time lattice Boltzmann method （MRT-LBM） is implemented to numerically simulate the cross flow over a longitudinal vibrating circular cylinder. This research is carried out on a three-dimensional （3D） finite cantilevered cylinder to investigate the effect of forced vibration on the wake characteristics and the 3D effect of a cantilevered cylinder. To meet the accuracy of this method, the present calculation is carried out at a low Reynolds number Re =100, as well as to make the vibration obvious, we make the vibration strong enough. The calculation results indicate that the vibration has significant influence on the wake characteristics. When the vibrating is big enough, our early works show that the 2D vortex shedding would be locked up by vibration. Contrarily, this phenomenon would not appear in the present 3D case because of the end effect of the cantilevered cylinder.     

Optimal low-dispersion low-dissipation LBM schemes for computational aeroacoustics

Lattice Boltmzann Methods (LBM) have been proved to be very effective methods for computational aeroacoustics (CAA), which have been used to capture the dynamics of weak acoustic fluctuations. In this paper, we propose a strategy to reduce the dispersive and disspative errors of the two-dimensional (2D) multi-relaxation-time lattice Boltzmann method (MRT-LBM). By presenting an effective algorithm, we obtain a uniform form of the linearized Navier–Stokes equations corresponding to the MRT-LBM in wave-number space. Using the matrix perturbation theory and the equivalent modified equation approach for finite difference methods, we propose a class of minimization problems to optimize the free-parameters in the MRT-LBM. We obtain this way a dispersion-relation-preserving LBM (DRP-LBM) to circumvent the minimized dispersion error of the MRT-LBM. The dissipation relation precision is also improved. And the stability of the MRT-LBM with the small bulk viscosity is guaranteed. Von Neuman analysis of the linearized MRT-LBM is performed to validate the optimized dispersion/dissipation relations considering monochromatic wave solutions. Meanwhile, dispersion and dissipation errors of the optimized MRT-LBM are quantitatively compared with the original MRT-LBM. Finally, some numerical simulations are carried out to assess the new optimized MRT-LBM schemes.     

多孔介质中流体流动及扩散的耦合格子Boltzmann模型

多孔介质中高Péclet数和大黏性比下混溶流体的流动和扩散广泛存在于二氧化碳驱油、化工生产等工业过程中.用数值方法对该问题进行研究时,关键在于如何正确描述高Péclet数和大黏性比下多孔介质内流体的行为.为此,提出了一种基于多松弛模型和格子动理模型的耦合格子Boltzmann模型.通过Chapman-Enskog分析,证明该模型能有效求解不可压Navier-Stokes方程和对流扩散方程.数值结果表明,该模型不仅具有二阶精度和良好的稳健性,而且对于高Péclet数和大黏性比的问题具有良好的数值稳定性,为模拟此类问题提供了有效工具.     

微尺度振荡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层 振荡流     

Stabilized finite element method for incompressible flows with high Reynolds number

In the following paper, we discuss the exhaustive use and implementation of stabilization finite element methods for the resolution of the 3D time-dependent incompressible Navier–Stokes equations. The proposed method starts by the use of a finite element variational multiscale (VMS) method, which consists in here of a decomposition for both the velocity and the pressure fields into coarse/resolved scales and fine/unresolved scales. This choice of decomposition is shown to be favorable for simulating flows at high Reynolds number. We explore the behaviour and accuracy of the proposed approximation on three test cases. First, the lid-driven square cavity at Reynolds number up to 50,000 is compared with the highly resolved numerical simulations and second, the lid-driven cubic cavity up to Re = 12,000 is compared with the experimental data. Finally, we study the flow over a 2D backward-facing step at Re = 42,000. Results show that the present implementation is able to exhibit good stability and accuracy properties for high Reynolds number flows with unstructured meshes.     

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.     

An interaction potential based lattice Boltzmann method with adaptive mesh refinement (AMR) for two-phase flow simulation

The lattice Boltzmann method (LBM) for two-phase flow simulation is often hindered by insufficient resolution at the interface. As a result, the LBM simulation of bubbles in bubbling flows is commonly limited to spherical or slightly deformed bubble shapes. In this study, the adaptive mesh refinement method for the LBM is developed to overcome such a problem. The approach for this new method is based on the improved interaction potential model, which is able to maintain grid-independent fluid properties in the two-fluid phases and at the interface. The LBM–AMR algorithm is described, especially concerning the LBM operation on a non-uniform mesh and the improved interaction potential model. Numerical simulations have been performed to validate the method in both single phase and multiphase flows. The 2D and 3D simulations of the buoyant rise of bubbles are conducted under various conditions. The agreement between the simulated bubble shape and velocity with experiments illustrates the capability of the LBM–AMR approach in predicting bubble dynamics even under the large bubble deformation conditions. Further, the LBM–AMR technique is capable of simulating a complex topology change of the interface. Integration of LBM with AMR can significantly improve the accuracy and reduce computation cost. The method developed in this study may appreciably enhance the capability of LBM in the simulation of complex multiphase flows under realistic conditions.     

Finite-Difference Lattice Boltzmann Scheme for High-Speed Compressible Flow: Two-Dimensional Case

Lattice Boltzmann (LB) modeling of high-speed compressible flows has long been attempted by various authors. One common weakness of most of previous models is the instability problem when the Mach number of the flow is large. In this paper we present a finite-difference LB model, which works for flows with flexible ratios of specific heats and a wide range of Mach number, from $0$ to 30 or higher. Besides the discrete-velocity-model by Watari [Physica A 382 (2007) 502], a modified Lax--Wendroff finite difference scheme and an artificial viscosity are introduced. The combination of the finite-difference scheme and the adding of artificial viscosity must find a balance of numerical stability versus accuracy. The proposed model is validated by recovering results of some well-known benchmark tests: shock tubes and shock reflections. The new model may be used to track shock waves and/or to study the non-equilibrium procedure in the transition between the regular and Mach reflections of shock waves, etc.     

Dynamic subgrid scale modeling of turbulent flows using lattice-Boltzmann method

In this paper, we discuss the incorporation of dynamic subgrid scale (SGS) models in the lattice-Boltzmann method (LBM) for large-eddy simulation (LES) of turbulent flows. The use of a dynamic procedure, which involves sampling or test-filtering of super-grid turbulence dynamics and subsequent use of scale-invariance for two levels, circumvents the need for empiricism in determining the magnitude of the model coefficient of the SGS models. We employ the multiple relaxation times (MRT) formulation of LBM with a forcing term, which has improved physical fidelity and numerical stability achieved by proper separation of relaxation time scales of hydrodynamic and non-hydrodynamic modes, for simulation of the grid-filtered dynamics of large-eddies. The dynamic procedure is illustrated for use with the common Smagorinsky eddy-viscosity SGS model, and incorporated in the LBM kinetic approach through effective relaxation time scales. The strain rate tensor in the SGS model is locally computed by means of non-equilibrium moments of the MRT-LBM. We also discuss proper sampling techniques or test-filters that facilitate implementation of dynamic models in the LBM. For accommodating variable resolutions, we employ conservative, locally refined grids in this framework. As examples, we consider the canonical anisotropic and inhomogeneous turbulent flow problem, i.e. fully-developed turbulent channel flow at two different shear Reynolds numbers Re_∗ of 180 and 395. The approach is able to automatically and self-consistently compute the values of the Smagorinsky coefficient, C_S. In particular, the computed value in the outer or bulk flow region, where turbulence is generally more isotropic, is about 0.155 (or the model coefficient ) which is in good agreement with prior data. It is also shown that the model coefficient becomes smaller and approaches towards zero near walls, reflecting the dampening of turbulent length scales near walls. The computed turbulence statistics at these Reynolds numbers are also in good agreement with prior data. The paper also discusses a procedure for incorporation of more general scale-similarity based SGS stress models.     

Three-dimensional finite-difference lattice Boltzmann model and its application to inviscid compressible flows with shock waves

In this paper, a three-dimensional (3D) finite-difference lattice Boltzmann model for simulating compressible flows with shock waves is developed in the framework of the double-distribution-function approach. In the model, a density distribution function is adopted to model the flow field, while a total energy distribution function is adopted to model the temperature field. The discrete equilibrium density and total energy distribution functions are derived from the Hermite expansions of the continuous equilibrium distribution functions. The discrete velocity set is obtained by choosing the abscissae of a suitable Gauss–Hermite quadrature with sufficient accuracy. In order to capture the shock waves in compressible flows and improve the numerical accuracy and stability, an implicit–explicit finite-difference numerical technique based on the total variation diminishing flux limitation is introduced to solve the discrete kinetic equations. The model is tested by numerical simulations of some typical compressible flows with shock waves ranging from 1D to 3D. The numerical results are found to be in good agreement with the analytical solutions and/or other numerical results reported in the literature.     

Assessment of localized artificial diffusivity scheme for large-eddy simulation of compressible turbulent flows

The localized artificial diffusivity method is investigated in the context of large-eddy simulation of compressible turbulent flows. The performance of different artificial bulk viscosity models are evaluated through detailed results from the evolution of decaying compressible isotropic turbulence with eddy shocklets and supersonic turbulent boundary layer. Effects of subgrid-scale (SGS) models and implicit time-integration scheme/time-step size are also investigated within the framework of the numerical scheme used. The use of a shock sensor along with artificial bulk viscosity significantly improves the scheme for simulating turbulent flows involving shocks while retaining the shock-capturing capability. The proposed combination of Ducros-type sensor with a negative dilatation sensor removes unnecessary bulk viscosity within expansion and weakly compressible turbulence regions without shocks and allows it to localize near the shocks. It also eliminates the need for a wall-damping function for the bulk viscosity while simulating wall-bounded turbulent flows. For the numerical schemes used, better results are obtained without adding an explicit SGS model than with SGS model at moderate Reynolds number. Inclusion of a SGS model in addition to the low-pass filtering and artificial bulk viscosity results in additional damping of the resolved turbulence. However, investigations at higher Reynolds numbers suggest the need for an explicit SGS model. The flow statistics obtained using the second-order implicit time-integration scheme with three sub-iterations closely agrees with the explicit scheme if the maximum Courant–Friedrichs–Lewy is kept near unity.     

格子Boltzmann亚格子模型的研究

为了将格子Boltzmann法应用于大雷诺数流动的模拟,本文将Smagorinsky亚格子模型和LBGK模型相结合,并对该亚格子LBM模型进行了研究。利用该亚格子LBM模型,对二维顶盖驱动流进行了模拟,得到了若干大雷诺数下流线图和方腔中心线上无量纲速度分布。计算结果与基准解进行比较,两者相互吻合。     

二维顶盖驱动半圆腔内流动MRT-LBM研究

本文采用多弛豫时间格子玻尔兹曼方法(multi relaxation time lattice Boltzmann method—MRT LBM)对二维顶盖驱动半圆腔内流动进行了数值模拟,得到了雷诺数为500～50000范围内半圆腔内流场分布情况。在二维顶盖驱动半圆腔流场中,随雷诺数的增大,流场内旋涡的数目逐渐增加,且流动依次呈现出稳定流、周期流、混沌流等状态。本文计算结果表明,MRT-LBM模型可显著提高计算的稳定性,适用于大范围的雷诺数流动情况。     

改进的物理粘性SPH方法及其在溃坝问题中的应用

在低雷诺数物理粘性SPH方法基础上引入再生核粒子法进行密度重构,既避免了用人工粘性所导致的数值耗散问题,又提高了低雷诺数物理粘性SPH方法的数值稳定性;以溃坝问题为例,对比分析低雷诺数物理粘性SPH方法和本文方法的仿真结果表明,本文方法可有效消除数值不稳定,压强和速度分布更加光滑,粒子秩序更好,可应用于雷诺数较高或粘性不可忽略的流动问题.     

Existence of dual solutions and three-dimensional instability in helical pipe flow

Three-dimensional (3D) direct numerical simulations (DNS) of the viscous incompressible fluid flow through a helical pipe with circular cross section were performed. The flow is governed by three parameters: the Dean number (or the Reynolds number), curvature, and torsion. First, we obtained steady solutions by steady 3D calculations, where dual solutions were found, one was uniform in the pipe axial direction and the other varied very slowly, if torsion exceeded a critical value. Then, the instability of the steady solutions obtained was studied by unsteady 3D calculations. We obtained critical Reynolds numbers of steady to unsteady transition by observing the behaviors of the unsteady solutions. The present results of the critical Reynolds number nearly agreed with those by the 2D linear stability analysis (Yamamoto et al. [9]) except for the lowest critical Reynolds number region, where the present study gave the critical Reynolds number much less than that obtained by the 2D linear stability analysis. We found the vortical structures in the form of a standing wave slightly above the marginal instability state, which is a trigger of explosive 3D instability.     

An Immersed Boundary Method with Formal Second-Order Accuracy and Reduced Numerical Viscosity

A formally second-order accurate immersed boundary method is presented and tested in this paper. We apply this new scheme to simulate the flow past a circular cylinder and study the effect of numerical viscosity on the accuracy of the computation by comparing the numerical results with those of a first-order method. The numerical evidence shows that the new scheme has less numerical viscosity and is therefore a better choice for the simulation of high Reynolds number flows with immersed boundaries.     

Lattice Boltzmann Solver for Multiphase Flows: Application to High Weber and Reynolds Numbers

The lattice Boltzmann method, now widely used for a variety of applications, has also been extended to model multiphase flows through different formulations. While already applied to many different configurations in low Weber and Reynolds number regimes, applications to higher Weber/Reynolds numbers or larger density/viscosity ratios are still the topic of active research. In this study, through a combination of a decoupled phase-field formulation—the conservative Allen–Cahn equation—and a cumulant-based collision operator for a low-Mach pressure-based flow solver, we present an algorithm that can be used for higher Reynolds/Weber numbers. The algorithm was validated through a variety of test cases, starting with the Rayleigh–Taylor instability in both 2D and 3D, followed by the impact of a droplet on a liquid sheet. In all simulations, the solver correctly captured the flow dynamics andmatched reference results very well. As the final test case, the solver was used to model droplet splashing on a thin liquid sheet in 3D with a density ratio of 1000 and kinematic viscosity ratio of 15, matching the water/air system at We = 8000 and Re = 1000. Results showed that the solver correctly captured the fingering instabilities at the crown rim and their subsequent breakup, in agreement with experimental and numerical observations reported in the literature.