Small data existence for the Enskog equation inL 1

An existence theorem for the Enskog equation with small initial data is proved in anL ¹ setting. This type of result is not available for the Boltzmann equation.     

Global existence in L1 for the Enskog equation and convergence of the solutions to solutions of the Boltzmann equation

For the Enskog equation in a box an existence theorem is proved for initial data with finite mass, energy, and entropy. Then, by letting the diameter of the molecules go to zero, the weak convergence of solutions of the Enskog equation to solutions of the Boltzmann equation is proved.     

Lattice Boltzmann model for interface capturing of multiphase flows based on Allen–Cahn equation

A phase-field-based lattice Boltzmann model is proposed for the interface capturing of multi-phase flows based on the conservative Allen–Cahn equation(ACE).By adopting the improved form of a relaxation matrix and an equilibrium distribution function,the time derivative?t(φ)induced by recovering the diffusion term in ACE is eliminated.The conducted Chapman–Enskog analysis demonstrates that the correct conservative ACE is recovered.Four benchmark cases including Zalesak’s disk rotation,vortex droplet,droplet impact on thin film,and Rayleigh–Taylor instability are investigated to validate the proposed model.The numerical results indicate that the proposed model can accurately describe the complex interface deformation.     

Global existence inL 1 for the generalized Enskog equation

Various existence theorems are given for the generalized Enskog equation inR ³ and in a bounded spatial domain with periodic boundary conditions. A very general form of the geometric factorY is allowed, including an explicit space, velocity, and time dependence. The method is based on the existence of a Liapunov functional, an analog of theH-function in the Boltzmann equation, and utilizes a weak compactness argument inL ¹.     

Physical modeling of multiphase flow via lattice Boltzmann method: Numerical effects, equation of state and boundary conditions

The aims of the present paper are threefold. First, we further study the fast Fourier transform thermal lattice Boltzmann (FFT-TLB) model for van der Waals (VDW) fluids proposed in Phys. Rev. E, 2011, 84(4): 046715. We analyze the merits of the FFT approach over the traditional finite difference scheme and investigate the effects of smoothing factors on accuracy and stability in detail. Second, we incorporate the VDW equation of state with flexible parameters into the FFT-TLB model. As a result, the revised model may be used to handle multiphase flows with various critical densities and temperatures. Third, we design appropriate boundary conditions for systems with solid walls. These improvements, from the views of numerics and physics, significantly extend the application scope of the model in science and engineering.     

Thermodynamic Foundations of Kinetic Theory and Lattice Boltzmann Models for Multiphase Flows

This paper demonstrates that thermodynamically consistent lattice Boltzmann models for single-component multiphase flows can be derived from a kinetic equation using both Enskog's theory for dense fluids and mean-field theory for long-range molecular interaction. The lattice Boltzmann models derived this way satisfy the correct mass, momentum, and energy conservation equations. All the thermodynamic variables in these LBM models are consistent. The strengths and weaknesses of previous lattice Boltzmann multiphase models are analyzed.     

Nonequilibrium effects of reactive flow based on gas kinetic theory

How to accurately probe chemically reactive flows with essential thermodynamic nonequilibrium effects is an open issue. Via the Chapman–Enskog analysis, the local nonequilibrium particle velocity distribution function is derived from the gas kinetic theory. It is demonstrated theoretically and numerically that the distribution function depends on the physical quantities and derivatives, and is independent of the chemical reactions directly as the chemical time scale is longer than the molecular relaxation time. Based on the simulation results of the discrete Boltzmann model, the departure between equilibrium and nonequilibrium distribution functions is obtained and analyzed around the detonation wave. In addition, it has been verified for the first time that the kinetic moments calculated by summations of the discrete distribution functions are close to those calculated by integrals of their original forms.     

Kinetic Models for Granular Flow

The generalization of the Boltzmann and Enskog kinetic equations to allow inelastic collisions provides a basis for studies of granular media at a fundamental level. For elastic collisions the significant technical challenges presented in solving these equations have been circumvented by the use of corresponding model kinetic equations. The objective here is to discuss the formulation of model kinetic equations for the case of inelastic collisions. To illustrate the qualitative changes resulting from inelastic collisions the dynamics of a heavy particle in a gas of much lighter particles is considered first. The Boltzmann–Lorentz equation is reduced to a Fokker–Planck equation and its exact solution is obtained. Qualitative differences from the elastic case arise primarily from the cooling of the surrounding gas. The excitations, or physical spectrum, are no longer determined simply from the Fokker–Planck operator, but rather from a related operator incorporating the cooling effects. Nevertheless, it is shown that a diffusion mode dominates for long times just as in the elastic case. From the spectral analysis of the Fokker–Planck equation an associated kinetic model is obtained. In appropriate dimensionless variables it has the same form as the BGK kinetic model for elastic collisions, known to be an accurate representation of the Fokker–Planck equation. On the basis of these considerations, a kinetic model for the Boltzmann equation is derived. The exact solution for states near the homogeneous cooling state is obtained and the transport properties are discussed, including the relaxation toward hydrodynamics. As a second application of this model, it is shown that the exact solution for uniform shear flow arbitrarily far from equilibrium can be obtained from the corresponding known solution for elastic collisions. Finally, the kinetic model for the dense fluid Enskog equation is described.     

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.     

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

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

燃烧系统的离散Boltzmann建模与模拟研究进展

燃烧系统的诸多模拟依托于流体建模, 离散Boltzmann方法(discrete Boltzmann method, DBM) 是近年来发展起来的一种新的流体介观建模方法. 本文简要评述DBM发展的两个方向——Navier-Stokes等偏微分方程的数值逼近解法和复杂系统的微介观动理学建模. 主要介绍在燃烧系统模拟方面DBM已有的工作、新近的思路、与传统流体建模的异同以及近期的研究成果. 本文重点传递的信息为: 作为复杂系统微介观动理学建模出现的DBM在模拟过程中同时给出“流动”及其相伴随的、关系最密切的那部分“热动”非平衡效应; 它为燃烧等复杂系统中各类非平衡行为的描述、非平衡信息的提取、非平衡程度的度量提供了一种简洁、有效的方法; 它所提供的热动非平衡测量量有两类: 一类是直接比较分布函数和平衡态分布函数的动理学矩关系得到的, 一类是来自于Chapman-Enskog多尺度分析给出的热传导和黏性项. 基于第二类DBM, 可以实现(燃烧等)一大类复杂流体系统的多尺度物理建模.     

Comment on symmetry properties of the linear Enskog kinetic operators

The one-particle average consistent with the structure of the revised Enskog theory is introduced. Symmetry properties of the linear kinetic operators reflecting those of theN-particle pseudo-Liouville operators are derived, implying a recently proved symmetry of kinetic expressions for equilibrium time correlation functions.     

Computational modeling of MR flow imaging by the lattice Boltzmann method and Bloch equation

In this work, a computational model of magnetic resonance (MR) flow imaging is proposed. The first model component provides fluid dynamics maps by applying the lattice Boltzmann method. The second one uses the flow maps and couples MR imaging (MRI) modeling with a new magnetization transport algorithm based on the Eulerian coordinate approach. MRI modeling is based on the discrete time solution of the Bloch equation by analytical local magnetization transformations (exponential scaling and rotations).     

The Limiting Kinetic Equation of the Consistent Boltzmann Algorithm for Dense Gases

This paper establishes a theoretical foundation for the Consistent Boltzmann Algorithm (CBA) by deriving the limiting kinetic equation. The formulation is similar to the proof by one of the authors that the Boltzmann equation is the limiting kinetic equation for Direct Simulation Monte Carlo [W. Wagner, J. Statist. Phys. 66:1011 (1992)]. For a simplified model distilled from CBA, the limiting equation is solved numerically, and very good agreement with the predictions of the theory is found.     

Numerical analysis of a shock-wave solution of the Enskog equation obtained via a Monte Carlo method

In this paper a planar stationary shock-wave-like solution of the Enskog equation obtained via a Monte Carlo technique is studied; both the algorithm used to obtain the solution and the qualitative behavior of the macroscopic quantities are discussed in comparison with the corresponding solution of the Boltzmann equation.     

Approximation of the Boltzmann equation by discrete velocity models

Two convergence results related to the approximation of the Boltzmann equation by discrete velocity models are presented. First we construct a sequence of deterministic discrete velocity models and prove convergence (as the number of discrete velocities tends to infinity) of their solutions to the solution of a spatially homogeneous Boltzmann equation. Second we introduce a sequence of Markov jump processes (interpreted as random discrete velocity models) and prove convergence (as the intensity of jumps tends to infinity) of these processes to the solution of a deterministic discrete velocity model.     

Some comments on modeling the linearized Boltzmann equation

Some exact solutions of the homogeneous and the inhomogeneous linearized Boltzmann equation (LBE) for rigid-sphere collisions are used to define two model equations in the general area of rarefied-gas dynamics. These equations are obtained from a systematic development of two synthetic scattering kernels that yield model equations that have as exact solutions certain known exact solutions of the homogeneous and of the inhomogeneous LBE. The first model established is defined in terms of the collisional invariants and the Chapman-Enskog integral equations for viscosity and for heat conduction. An extended model is defined also in terms of the collisional invariants and the Chapman-Enskog functions for viscosity and heat conduction, but the first and second Burnett functions are also included in the model. The variable collision frequency or generalized BGK model is also obtained as a special case. In addition, the exact mean-free paths defined, for rigid-sphere collisions and the LBE, in terms of viscosity or heat conduction are employed to define approximations of these quantities that are consistent with the use of the variable collision frequency model.     

Gas kinetic algorithm for flows in Poiseuille-like microchannels using Boltzmann model equation

~~Gas kinetic algorithm for flows in Poiseuille-like microchannels using Boltzmann model equation1. Feynman, R., There's plenty of room at the bottom, Journal of Microelectromechanical Systems, 1992, 1: 60 -66. 2. Piekos, E. S., Breuer, K. S., Numerical modeling of micromechanical devices using the direct simulation Monte Carlo method, Transactions of the ASME, Journal of Fluids Engineering, 1996, 118: 464-469. 3. Beskok, A., Karniadakis, G. E., Trimmer, W., Rarefaction and …     

Study on the unified algorithm for three-dimensional complex problems covering various flow regimes using Boltzmann model equation

The Boltzmann simplified velocity distribution function equation describing the gas transfer phenomena from various flow regimes will be explored and solved numerically in this study. The discrete velocity ordinate method of the gas kinetic theory is studied and applied to simulate the complex multi-scale flows. Based on the uncoupling technique on molecular movement and colliding in the DSMC method, the gas-kinetic finite difference scheme is constructed to directly solve the discrete velocity distribution functions by extending and applying the unsteady time-splitting method from computational fluid dynamics. The Gauss-type discrete velocity numerical quadrature technique for different Mach number flows is developed to evaluate the macroscopic flow parameters in the physical space. As a result, the gas-kinetic numerical algorithm is established to study the three-dimensional complex flows from rarefied transition to continuum regimes. The parallel strategy adapted to the gas-kinetic numerical algorithm is investigated by analyzing the inner parallel degree of the algorithm, and then the HPF parallel processing program is developed. To test the reliability of the present gas-kinetic numerical method, the three-dimensional complex flows around sphere and spacecraft shape with various Knudsen numbers are simulated by HPF parallel computing. The computational results are found in high resolution of the flow fields and good agreement with the theoretical and experimental data. The computing practice has confirmed that the present gas-kinetic algorithm probably provides a promising approach to resolve the hypersonic aerothermodynamic problems with the complete spectrum of flow regimes from the gas-kinetic point of view of solving the Boltzmann model equation. Supported by the National Natural Science Foundation of China (Grant Nos. 90205009 and 10321002) and the National Parallel Computing Center     

介观尺度通道内多相流动的耗散粒子动力学模拟

采用四次方光滑函数构造了包含远程吸引近距排斥的保守力势函数, 应用该势函数, 采用耗散粒子动力学方法, 对十字型微通道内的多相流动进行了计算. 计算结果表明, 该势函数能有效模拟十字型微通道内的流动过程及模式. 关键词： 多相流 耗散粒子动力学 保守力势函数