A Relaxation Scheme for Solving the Boltzmann Equation Based on the Chapman-Enskog Expansion

Abstract In [16] a visco-elastic relaxation system, called the relaxed Burnett system, was proposed by Jinand Slemrod as a moment approximation to the Boltzmann equation. The relaxed Burnett system is weaklyparabolic, has a linearly hyperbolic convection part, and is endowed with a generalized eotropy inequality. Itagrees with the solution of the Boltzmann equation up to the Burnett order via the Chapman-Enskog expansion. We develop a one-dimensional non-oscillatory numerical scheme based on the relaxed Burnett system forthe Boltzmann equation. We compare numerical results for stationary shocks based on this relaxation scheme,and those obtained by the DSMC (Direct Simulation Monte Carlo), by the Navier-Stokes equations and bythe extended thermodynamics with thirteen moments (the Grad equations). Our numerical experiments showthat the relaxed Burnett gives more accurate approximations to the shock profiles of the Boltzmann equationobtained by the DSMC, for a range of Mach numbers for hypersonic flows, th     

A higher-order moment method of the lattice Boltzmann model for the conservation law equation

In this paper, we proposed a higher-order moment method in the lattice Boltzmann model for the conservation law equation. In contrast to the lattice Bhatnagar–Gross–Krook (BGK) model, the higher-order moment method has a wide flexibility to select equilibrium distribution function. This method is based on so-called a series of partial differential equations obtained by using multi-scale technique and Chapman–Enskog expansion. According to Hirt’s heuristic stability theory, the stability of the scheme can be controlled by modulating some special moments to design the third-order dispersion term and the fourth-order dissipation term. As results, the conservation law equation is recovered with higher-order truncation error. The numerical examples show the higher-order moment method can be used to raise the accuracy of the truncation error of the lattice Boltzmann scheme for the conservation law equation.     

Discontinuous Shock Structure in a Reacting Mixture Modelled by Grad 13 Moment Approximation

A preliminary analysis on the possible occurrence of sub-shocks into a gas mixture is carried out. The mixture, undergoing a reversible bimolecular reaction, is described by macroscopic equations obtained by Grad 13 moment approximation of the reactive Boltzmann equation.     

A new lattice Boltzmann model for the laplace equation

In this paper, a new lattice Boltzmann equation which is independent of time is proposed. Based on the new lattice Boltzmann equation, some steady problems can be modeled by the lattice Boltzmann method. In the further study, the Laplace equation is investigated with the method of the higher-order moment of equilibrium distribution functions and a series of partial differential equations in different space scales. The numerical results show that the new method is effective.     

Diffusive semiconductor moment equations using Fermi�CDirac statistics

Diffusive moment equations with an arbitrary number of moments are formally derived from the semiconductor Boltzmann equation employing a moment method and a Chapman?CEnskog expansion. The moment equations are closed by employing a generalized Fermi?CDirac distribution function obtained from entropy maximization. The current densities allow for a drift-diffusion-type formulation or a ??symmetrized?? formulation, using dual-entropy variables from nonequilibrium thermodynamics. Furthermore, drift-diffusion and new energy-transport equations based on Fermi?CDirac statistics are obtained and their degeneracy limit is studied.     

Center manifolds for a class of degenerate evolution equations and existence of small-amplitude kinetic shocks

We construct center manifolds for a class of degenerate evolution equations including the steady Boltzmann equation and related kinetic models, establishing in the process existence and behavior of small-amplitude kinetic shock and boundary layers. Notably, for Boltzmann's equation, we show that elements of the center manifold decay in velocity at near-Maxwellian rate, in accord with the formal Chapman–Enskog picture of near-equilibrium flow as evolution along the manifold of Maxwellian states, or Grad moment approximation via Hermite polynomials in velocity. Our analysis is from a classical dynamical systems point of view, with a number of interesting modifications to accommodate ill-posedness of the underlying evolution equation.     

Propagation of and Maxwellian weighted bounds for derivatives of solutions to the homogeneous elastic Boltzmann equation

We consider the n-dimensional space homogeneous Boltzmann equation for elastic collisions for variable hard potentials with Grad (angular) cutoff. We prove sharp moment inequalities, the propagation of L¹-Maxwellian weighted estimates, and consequently, the propagation L^∞-Maxwellian weighted estimates to all derivatives of the initial value problem associated to the afore mentioned equation. More specifically, we extend to all derivatives of the initial value problem associated to this class of Boltzmann equations corresponding sharp moment (Povzner) inequalities and time propagation of L¹-Maxwellian weighted estimates as originally developed Bobylev [A.V. Bobylev, Moment inequalities for the Boltzmann equation and applications to spatially homogeneous problems, J. Statist. Phys. 88 (1997) 1183–1214] in the case of hard spheres in 3 dimensions. To achieve this goal we implement the program presented in Bobylev–Gamba–Panferov [A.V. Bobylev, I.M. Gamba, V. Panferov, Moment inequalities and high-energy tails for Boltzmann equation with inelastic interactions, J. Statist. Phys. 116 (5–6) (2004) 1651–1682], which includes a full analysis of the moments by means of sharp moment inequalities and the control of L¹-exponential bounds, in the case of stationary states for different inelastic Boltzmann related problems with ‘heating’ sources where high energy tail decay rates depend on the inelasticity coefficient and the type of ‘heating’ source. More recently, this work was extended to variable hard potentials with angular cutoff by Gamba–Panferov–Villani [I.M. Gamba, V. Panferov, C. Villani, Upper Maxwellian bounds for the spatially homogeneous Boltzmann equation, ARMA (2008), in press] in the elastic case collision case where the L¹-Maxwellian weighted norm was shown to propagate if initial states have such property. In addition, we also extend to all derivatives the propagation of L^∞-Maxwellian weighted estimates, proven in [I.M. Gamba, V. Panferov, C. Villani, Upper Maxwellian bounds for the spatially homogeneous Boltzmann equation, ARMA (2008), in press], to solutions of the initial value problem to the Boltzmann equations for elastic collisions for variable hard potentials with Grad (angular) cutoff.     

Truncating the hierarchy of moment equations based on point distribution—application to innovation diffusion

The phenomenon of innovation diffusion, modeled as a nonlinear birth process, leads to a hierarchy of moment equations. For gaining insight into the evolution of moments of a number of adopters in innovation diffusion, truncation procedures based on point distributions are proposed for truncating the hierarchy of moment equations. It is found that the results obtained by employing new truncation procedures are in close agreement with analytical results based on system size expansion for large population size.     

从离散速度模型到矩方法

为了求解动理学方程,我们通过研究一维情形下的离散速度模型,发现通过对离散速度点使用自适应技术可以直接得到Grad矩方程组.作为一个统一的认识,矩方程组可以看作是对离散速度点自适应的离散速度模型,而离散速度模型可以看作是取特别形式的"矩"的矩方程组.这使得我们可以在一致的框架下来理解离散速度模型和矩方法,而不是将它们对立起来.为了建立这样的一致框架,最近在[2]中发展的正则化理论是根本性的.     

Lattice Boltzmann model for the bimolecular autocatalytic reaction–diffusion equation

A lattice Boltzmann model for the bimolecular autocatalytic reaction–diffusion equation is proposed. By using multi-scale technique and the Chapman–Enskog expansion on complex lattice Boltzmann equation, we obtain a series of complex partial differential equations, complex equilibrium distribution function and its complex moments. Then, the complex reaction–diffusion equation is recovered with higher-order accuracy of the truncation error. This equation can be used to describe the bimolecular autocatalytic reaction–diffusion systems, in which a rich variety of behaviors have been observed. Based on this model, the Fitzhugh–Nagumo model and the Gray–Scott model are simulated. The comparisons between the LBM results and the Alternative Direction Implicit results are given in detail. The numerical examples show that assumptions of source term can be used to raise the accuracy of the truncation error of the lattice Boltzmann scheme for the complex reaction–diffusion equation.     

The fluid dynamic limit of the nonlinear boltzmann equation

Solutions of the nonlinear Boltzmann equation are constructed up to the first appearance of shocks in the corresponding fluid dynamics. This construction assumes the knowledge of solutions of the Euler equations for compressible gas flow. The Boltzmann solution is found as a truncated Hilbert expansion with a remainder, and the remainder term solves a weakly nonlinear equation which is solved by iteration. The solutions found have special initial values. They should serve as “outer expansions” to which initial layers, boundary layers and shock layers can be matched.     

Global well-posedness of the velocity–vorticity-Voigt model of the 3D Navier–Stokes equations

The velocity–vorticity formulation of the 3D Navier–Stokes equations was recently found to give excellent numerical results for flows with strong rotation. In this work, we propose a new regularization of the 3D Navier–Stokes equations, which we call the 3D velocity–vorticity-Voigt (VVV) model, with a Voigt regularization term added to momentum equation in velocity–vorticity form, but with no regularizing term in the vorticity equation. We prove global well-posedness and regularity of this model under periodic boundary conditions. We prove convergence of the model's velocity and vorticity to their counterparts in the 3D Navier–Stokes equations as the Voigt modeling parameter tends to zero. We prove that the curl of the model's velocity converges to the model vorticity (which is solved for directly), as the Voigt modeling parameter tends to zero. Finally, we provide a criterion for finite-time blow-up of the 3D Navier–Stokes equations based on this inviscid regularization.     

Asymptotic Stability of Solutions of the Cauchy Problem for Models of Nonequilibrium Thermodynamics. Stable Hyperbolic Pencils

The paper distinguishes a class of stable polynomial pencils reproduced by first moments in the Grad hierarchy for the kinetic Boltzmann and Fokker-Planck equations. __________ Translated from Sovremennaya Matematika i Ee Prilozheniya (Contemporary Mathematics and Its Applications), Vol. 12, Partial Differential Equations, 2004.     

Derivation of the Boltzmann equation and entropy production in functional mechanics

A derivation of the Boltzmann equation from the Liouville equation by the use of the Grad limiting procedure in a finite volume is proposed. We introduce two scales of space-time: macro- and microscale and use the BBGKY hierarchy and the functional formulation of classical mechanics. According to the functional approach to mechanics, a state of a system of particles is formed from the measurements, which are rational numbers. Hence, one can speak about the accuracy of the initial probability density function in the Liouville equation. We assume that the initial data for the microscopic density functions are assigned by the macroscopic one (so one can say about a kind of hierarchy and subordination of the microscale to the macroscale) and derive the Boltzmann equation, which leads to the entropy production.     

用格子Boltzmann方法模拟非线性热传导方程

对具有非线性源项和非线性扩散项的热传导方程建立格子Boltzmann求解模型.在演化方程中增加了两个关于源项分布函数的微分算子,对演化方程实施Chapman-Enskog展开.通过对演化方程的进一步改进,恢复出具有高阶截断误差的宏观方程.对不同参数选取下的非线性热传导方程进行了数值模拟,数值解与精确解吻合得很好.该模型也可以用于同类型的其他偏微分方程的数值计算中.     

Semiclassical lattice hydrodynamics of rarefied channel flows

The two-dimensional channel flows of gas of arbitrary statistics in the slip and transition regimes as characterized by the Knudsen number are studied using a newly developed semiclassical lattice Boltzmann method. The method is directly derived by projecting the Uehling-Uhlenbeck Boltzmann-BGK equations onto the tensor Hermite polynomials using moment expansion method. The intrinsic discrete nodes of the Gauss-Hermite quadrature provide the natural lattice velocities for the semiclassical lattice Boltzmann method. The mass flow rates and the velocity profiles are calculated for the three particle statistics over wide range of Knudsen numbers and the Knudsen minimum can be captured. The results indicate distinct characteristics of the effects of quantum statistics.     

A robust multigrid approach for variational image registration models

Variational registration models are non-rigid and deformable imaging techniques for accurate registration of two images. As with other models for inverse problems using the Tikhonov regularization, they must have a suitably chosen regularization term as well as a data fitting term. One distinct feature of registration models is that their fitting term is always highly nonlinear and this nonlinearity restricts the class of numerical methods that are applicable. This paper first reviews the current state-of-the-art numerical methods for such models and observes that the nonlinear fitting term is mostly ‘avoided’ in developing fast multigrid methods. It then proposes a unified approach for designing fixed point type smoothers for multigrid methods. The diffusion registration model (second-order equations) and a curvature model (fourth-order equations) are used to illustrate our robust methodology. Analysis of the proposed smoothers and comparisons to other methods are given. As expected of a multigrid method, being many orders of magnitude faster than the unilevel gradient descent approach, the proposed numerical approach delivers fast and accurate results for a range of synthetic and real test images.     

The Boltzmann–Grad limit for a particle system in a slab with stochastic walls

It is shown that a stochastic system of N interacting particles in a slab approximates, in the Boltzmann–Grad limit, a one-dimensional Boltzmann equation with diffusive boundary conditions. © 1998 B. G. Teubner Stuttgart–John Wiley & Sons, Ltd.     

Global existence of weak solutions for a system of non-linear Boltzmann equations in semiconductor physics

In this paper we give a proof of the global existence of weak solutions for the semiconductor Boltzmann equation. This equation rules the evolution of the distribution function of carriers in the kinetic model of semiconductors. The main tool for the proof consists of a recent compactness result on velocity averages of solutions of transport equations. This result needs a L²-estimate of the electric field, which is obtained from the energy estimate, using the original regularization procedure of the problem, proposed in this paper.     

A robust multigrid approach for variational image registration models

Variational registration models are non-rigid and deformable imaging techniques for accurate registration of two images. As with other models for inverse problems using the Tikhonov regularization, they must have a suitably chosen regularization term as well as a data fitting term. One distinct feature of registration models is that their fitting term is always highly nonlinear and this nonlinearity restricts the class of numerical methods that are applicable. This paper first reviews the current state-of-the-art numerical methods for such models and observes that the nonlinear fitting term is mostly ‘avoided’ in developing fast multigrid methods. It then proposes a unified approach for designing fixed point type smoothers for multigrid methods. The diffusion registration model (second-order equations) and a curvature model (fourth-order equations) are used to illustrate our robust methodology. Analysis of the proposed smoothers and comparisons to other methods are given. As expected of a multigrid method, being many orders of magnitude faster than the unilevel gradient descent approach, the proposed numerical approach delivers fast and accurate results for a range of synthetic and real test images.