A Hybrid Lattice Boltzmann Flux Solver for Simulation of Viscous Compressible Flows

In this paper, a hybrid lattice Boltzmann flux solver (LBFS) is proposed for simulation of viscous compressible flows. In the solver, the finite volume method is applied to solve the Navier-Stokes equations. Different from conventional Navier-Stokes solvers, in this work, the inviscid flux across the cell interface is evaluated by local reconstruction of solution using one-dimensional lattice Boltzmann model, while the viscous flux is still approximated by conventional smooth function approximation. The present work overcomes the two major drawbacks of existing LBFS [28–31], which is used for simulation of inviscid flows. The first one is its ability to simulate viscous flows by including evaluation of viscous flux. The second one is its ability to effectively capture both strong shock waves and thin boundary layers through introduction of a switch function for evaluation of inviscid flux, which takes a value close to zero in the boundary layer and one around the strong shock wave. Numerical experiments demonstrate that the present solver can accurately and effectively simulate hypersonic viscous flows.     

基于浸入边界-多松弛时间格子玻尔兹曼通量求解法的流固耦合算法研究

针对流固耦合问题,发展了基于浸入边界-多松弛时间格子玻尔兹曼通量求解法(immersed boundary method multi-relaxation-time lattice Boltzmann flux solver,IB-MRT-LBFS)的弱耦合算法.依据多尺度Chapman-Enskog展开,建立不可压宏观方程状态变量和通量与格子玻尔兹曼方程中粒子密度分布函数之间的关系;采用强制浸入边界法处理流固界面使固壁表面满足无滑移边界条件,根据修正的速度求解动量方程力源项;结构运动方程采用四阶龙格-库塔法求解.格子模型与浸入边界法的引入使流固耦合计算可以在笛卡尔网格下进行,无需生成贴体网格及运用动网格技术,简化了计算过程.数值模拟了单圆柱横向涡激振动、单圆柱及串列双圆柱双自由度涡激振动问题.结果表明,IB-MRT-LBFS能够准确预测圆柱涡激振动的锁定区间、振动响应、受力情况以及捕捉尾流场结构形态,验证了该算法在求解流固耦合问题的有效性和可行性.     

High-resolution viscous flow simulations at arbitrary Mach number

A characteristic-based unsteady viscous flow solver is developed with preconditioning that is uniformly applicable for Mach numbers ranging from essentially incompressible to supersonic. A preconditioned flux-difference formulation for nondimensional primitive variables is a key element of the present approach. The simple primitive-variable numerical flux is related to Roe’s flux-difference scheme and preserves contact discontinuities using primitive variables, with or without preconditioning. Preconditioning by a single-parameter diagonal matrix conditions the system eigenvalues in terms of nondimensional local velocity and local temperature. An iterative implicit solution algorithm is given for the preconditioned formulation and is used for several simple test and validation cases. These include an inviscid shock-tube case, flat-plate boundary layer flow at low Mach number, viscous flow past a circular cylinder at low Reynolds number and with different thermal boundary conditions, and validation cases for incompressible and transonic flows.     

任意马赫数非定常流动数值模拟的统一算法

发展适用于从低速到高速任意马赫数非定常流动数值模拟的统一算法.通过引入一个伪时间导数项和一个新的预处理矩阵,得到双时间非定常预处理可压缩Navier-Stokes方程.方程的对流项采用三阶Roe通量近似差分格式离散,粘性项采用二阶中心差分格式离散.基于数值通量的线性化技术,实现伪时间步的隐式ADI-LU格式迭代,进而获得物理时间步的二阶推进精度.重点以低马赫数流动为例,求解了圆柱绕流和NACA0015翼型等速上仰动态失速问题.计算结果表明该统一算法能够较好地模拟低马赫数乃至任意马赫数非定常流动.     

Nested Cartesian grid method in incompressible viscous fluid flow

In this work, the local grid refinement procedure is focused by using a nested Cartesian grid formulation. The method is developed for simulating unsteady viscous incompressible flows with complex immersed boundaries. A finite-volume formulation based on globally second-order accurate central-difference schemes is adopted here in conjunction with a two-step fractional-step procedure. The key aspects that needed to be considered in developing such a nested grid solver are proper imposition of interface conditions on the nested-block boundaries, and accurate discretization of the governing equations in cells that are with block-interface as a control-surface. The interpolation procedure adopted in the study allows systematic development of a discretization scheme that preserves global second-order spatial accuracy of the underlying solver, and as a result high efficiency/accuracy nested grid discretization method is developed. Herein the proposed nested grid method has been widely tested through effective simulation of four different classes of unsteady incompressible viscous flows, thereby demonstrating its performance in the solution of various complex flow–structure interactions. The numerical examples include a lid-driven cavity flow and Pearson vortex problems, flow past a circular cylinder symmetrically installed in a channel, flow past an elliptic cylinder at an angle of attack, and flow past two tandem circular cylinders of unequal diameters. For the numerical simulations of flows past bluff bodies an immersed boundary (IB) method has been implemented in which the solid object is represented by a distributed body force in the Navier–Stokes equations. The main advantages of the implemented immersed boundary method are that the simulations could be performed on a regular Cartesian grid and applied to multiple nested-block (Cartesian) structured grids without any difficulty. Through the numerical experiments the strength of the solver in effectively/accurately simulating various complex flows past different forms of immersed boundaries is extensively demonstrated, in which the nested Cartesian grid method was suitably combined together with the fractional-step algorithm to speed up the solution procedure.     

A high-order gas-kinetic Navier–Stokes flow solver

The foundation for the development of modern compressible flow solver is based on the Riemann solution of the inviscid Euler equations. The high-order schemes are basically related to high-order spatial interpolation or reconstruction. In order to overcome the low-order wave interaction mechanism due to the Riemann solution, the temporal accuracy of the scheme can be improved through the Runge–Kutta method, where the dynamic deficiencies in the first-order Riemann solution is alleviated through the sub-step spatial reconstruction in the Runge–Kutta process. The close coupling between the spatial and temporal evolution in the original nonlinear governing equations seems weakened due to its spatial and temporal decoupling. Many recently developed high-order methods require a Navier–Stokes flux function under piece-wise discontinuous high-order initial reconstruction. However, the piece-wise discontinuous initial data and the hyperbolic-parabolic nature of the Navier–Stokes equations seem inconsistent mathematically, such as the divergence of the viscous and heat conducting terms due to initial discontinuity. In this paper, based on the Boltzmann equation, we are going to present a time-dependent flux function from a high-order discontinuous reconstruction. The theoretical basis for such an approach is due to the fact that the Boltzmann equation has no specific requirement on the smoothness of the initial data and the kinetic equation has the mechanism to construct a dissipative wave structure starting from an initially discontinuous flow condition on a time scale being larger than the particle collision time. The current high-order flux evaluation method is an extension of the second-order gas-kinetic BGK scheme for the Navier–Stokes equations (BGK-NS). The novelty for the easy extension from a second-order to a higher order is due to the simple particle transport and collision mechanism on the microscopic level. This paper will present a hierarchy to construct such a high-order method. The necessity to couple spatial and temporal evolution nonlinearly in the flux evaluation can be clearly observed through the numerical performance of the scheme for the viscous flow computations.     

Entropic lattice Boltzmann method for turbulent flow simulations: Boundary conditions

We introduce a new concept of boundary conditions for realization of the lattice Boltzmann simulations of turbulent flows. The key innovation is the use of a universal distribution function for particles, analogous to the Tamm–Mott-Smith solution for the shock wave in the classical Boltzmann kinetic equation. Turbulent channel flow simulations demonstrate that the new boundary enables accurate results even with severely under-resolved grids. Generalization to complex boundary is illustrated with an example of turbulent flow past a circular cylinder.     

Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier–Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.     

微通道中气体流动的格子Boltzmann数值模拟

采用可压缩格子Boltzmann模型及非平衡外推边界条件,数值模拟微通道中的气体在滑移区域(Kn≤0.1)内的流动,计算结果包括出口速度剖面、通道中心压力分布以及质量流率等,与理论结果及其他实验结果符合得很好.还模拟了180°弯曲通道中的气体流动.结果表明,滑移速度的存在抑制了边界层的分离,因此在弯曲处不存在漩涡.计算结果还表明,弯道的存在显著影响了气体的质量流率.     

A momentum exchange-based immersed boundary-lattice Boltzmann method for simulating incompressible viscous flows

A momentum exchange-based immersed boundary-lattice Boltzmann method is presented in this Letter for simulating incompressible viscous flows. This method combines the good features of the lattice Boltzmann method (LBM) and the immersed boundary method (IBM) by using two unrelated computational meshes, an Eulerian mesh for the flow domain and a Lagrangian mesh for the solid boundaries in the flow. In this method, the non-slip boundary condition is enforced by introducing a forcing term into the lattice Boltzmann equation (LBE). Unlike the conventional IBM using the penalty method with a user-defined parameter or the direct forcing scheme based on the Navier–Stokes (NS) equations, the forcing term is simply calculated by the momentum exchange of the boundary particle density distribution functions, which are interpolated by the Lagrangian polynomials from the underlying Eulerian mesh. Numerical examples show that the present method can provide very accurate numerical results.     

A Switch Function-Based Gas-Kinetic Scheme for Simulation of Inviscid and Viscous Compressible Flows

In this paper, a switch function-based gas-kinetic scheme (SF-GKS) is presented for the simulation of inviscid and viscous compressible flows. With the finite volume discretization, Euler and Navier-Stokes equations are solved and the SF-GKS is applied to evaluate the inviscid flux at cell interface. The viscous flux is obtained by the conventional smooth function approximation. Unlike the traditional gas-kinetic scheme in the calculation of inviscid flux such as Kinetic Flux Vector Splitting (KFVS), the numerical dissipation is controlled with a switch function in the present scheme. That is, the numerical dissipation is only introduced in the region around strong shock waves. As a consequence, the present SF-GKS can well capture strong shock waves and thin boundary layers simultaneously. The present SF-GKS is firstly validated by its application to the inviscid flow problems, including 1-D Euler shock tube, regular shock reflection and double Mach reflection. Then, SF-GKS is extended to solve viscous transonic and hypersonic flow problems. Good agreement between the present results and those in the literature verifies the accuracy and robustness of SF-GKS.     

基于格子Boltzmann方法的超声速非平衡流模拟

基于D1Q4可压缩格子Boltzmann模型,按照流通矢量分裂方法的思路,采用坐标旋转技术构造求解三维带化学反应Navier-Stokes方程对流通量求解器.结合有限体积法求解三维化学非平衡流Navier-Stokes方程,采用时间算子分裂算法解决化学反应刚性问题,数值模拟超声速化学非平衡流的三个经典算例.数值结果表明:在高马赫数下,采用D1Q4可压缩格子Boltzmann模型构造的三维对流通量求解器数值模拟中没有出现非物理解,同时在超声速化学非平衡流场中正确分辨激波、燃烧波等物理现象,精度和分辨率均较高,验证了本文构造的三维对流通量求解器的可靠性,拓宽了D1Q4可压缩格子Boltzmann模型的应用范围,为计算超声速化学非平衡流提供一种新方法.     

Lattice Boltzmann simulations of 2D laminar flows past two tandem cylinders

We apply the lattice Boltzmann equation (LBE) with multiple-relaxation-time (MRT) collision model to simulate laminar flows in two-dimensions (2D). In order to simulate flows in an unbounded domain with the LBE method, we need to address two issues: stretched non-uniform mesh and inflow and outflow boundary conditions. We use the interpolated grid stretching method to address the need of non-uniform mesh. We demonstrate that various inflow and outflow boundary conditions can be easily and consistently realized with the MRT-LBE. The MRT-LBE with non-uniform stretched grids is first validated with a number of test cases: the Poiseuille flow, the flow past a cylinder asymmetrically placed in a channel, and the flow past a cylinder in an unbounded domain. We use the LBE method to simulate the flow past two tandem cylinders in an unbounded domain with Re = 100. Our results agree well with existing ones. Through this work we demonstrate the effectiveness of the MRT-LBE method with grid stretching.     

通道中串列旋转圆柱的格子Boltzmann-虚拟区域方法的模拟

通过结合格子Boltzmann方法(LBM)和虚拟区域(Fictitiou sDomain)思想,建立格子Boltzmann-虚拟区域(LB-DF/FD)方法.采用两套网格系统,欧拉网格用于流体,拉格朗日网格用于固体.原有的LBM在计算运动固体的受力方面存在数据振荡,LB-DF/FD方法改进了此缺陷.为验证该方法,模拟圆柱绕流、圆形颗粒在无限长通道中平动及在无限大流场中转动三种情况,结果与其他数值解及理论解符合得很好.利用该方法模拟低雷诺数下通道中串列旋转圆柱周围的流场,分析圆柱间距(g)及雷诺数(Re)对流场结构的影响.给出Re=0.001,0.1和10下,0.2≤g≤8.0的流线结构、圆柱升力、阻力以及力矩等数值结果.结果表明,g对流场的结构及圆柱的受力有显著影响,Re对圆柱阻力及Stokes单元数目的影响较大.     

Lattice Boltzmann Simulation of Free-Surface Temperature Dispersion in Shallow Water Flows

We develop a lattice Boltzmann method for modeling free-surface temperature dispersion in the shallow water flows. The governing equations are derived from the incompressible Navier-Stokes equations with assumptions of shallow water flows including bed frictions, eddy viscosity, wind shear stresses and Coriolis forces. The thermal effects are incorporated in the momentum equation by using a Boussinesq approximation. The dispersion of free-surface temperature is modelled by an advection-diffusion equation. Two distribution functions are used in the lattice Boltzmann method to recover the flow and temperature variables using the same lattice structure. Neither upwind discretization procedures nor Riemann problem solvers are needed in discretizing the shallow water equations. In addition, the source terms are straightforwardly included in the model without relying on well-balanced techniques to treat flux gradients and source terms. We validate the model for a class of problems with known analytical solutions and we also present numerical results for sea-surface temperature distribution in the Strait of Gibraltar.     

Kinetic theory based lattice Boltzmann equation with viscous dissipation and pressure work for axisymmetric thermal flows

A lattice Boltzmann equation (LBE) for axisymmetric thermal flows is proposed. The model is derived from the kinetic theory which exhibits several features that distinguish it from other previous LBE models. First, the present thermal LBE model is derived from the continuous Boltzmann equation, which has a solid foundation and clear physical significance; Second, the model can recover the energy equation with the viscous dissipation term and work of pressure which are usually ignored by traditional methods and the existing thermal LBE models; Finally, unlike the existing thermal LBE models, no velocity and temperature gradients appear in the force terms which are easy to realize in the present model. The model is validated by thermal flow in a pipe, thermal buoyancy-driven flow, and swirling flow in vertical cylinder by rotating the top and bottom walls. It is found that the numerical results agreed excellently with analytical solution or other numerical results.     

A lattice Boltzmann method for immiscible multiphase flow simulations using the level set method

We consider the lattice Boltzmann method for immiscible multiphase flow simulations. Classical lattice Boltzmann methods for this problem, e.g. the colour gradient method or the free energy approach, can only be applied when density and viscosity ratios are small. Moreover, they use additional fields defined on the whole domain to describe the different phases and model phase separation by special interactions at each node. In contrast, our approach simulates the flow using a single field and separates the fluid phases by a free moving interface. The scheme is based on the lattice Boltzmann method and uses the level set method to compute the evolution of the interface. To couple the fluid phases, we develop new boundary conditions which realise the macroscopic jump conditions at the interface and incorporate surface tension in the lattice Boltzmann framework. Various simulations are presented to validate the numerical scheme, e.g. two-phase channel flows, the Young–Laplace law for a bubble and viscous fingering in a Hele-Shaw cell. The results show that the method is feasible over a wide range of density and viscosity differences.     

实数型格子气模型及方柱绕流的模拟

在4Bits的HPP格子气模型基础上提出了用实数表示节点每一速度方向的粒产密度的实数型HPP模型,建立了粒子碰撞的运算关系式(代替碰撞规则),用此方法对方柱绕流问题进行了计算并介绍了边界条件的处理方法,计算结果较好地反映了真实流动特性。     

A time-reversal lattice Boltzmann method

In this paper we address the time-reversed simulation of viscous flows by the lattice Boltzmann method (LB). The theoretical derivation of the reversed LB from the Boltzmann equation is detailed, and the method implemented for weakly compressible flows using the D2Q9 scheme. The implementation of boundary conditions is also discussed. The accuracy and stability are illustrated by four test cases, namely the propagation of an acoustic wave in a medium at rest and in an uniform mean flow, the Taylor–Green vortex decay and the vortex pair–wall collision.     

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.