Stability of contact discontinuity for the Boltzmann equation

The Boltzmann equation which describes the time evolution of a large number of particles through the binary collision in statistics physics has close relation to the systems of fluid dynamics, that is, Euler equations and Navier-Stokes equations. As for a basic wave pattern to Euler equations, we consider the nonlinear stability of contact discontinuities to the Boltzmann equation. Even though the stability of the other two nonlinear waves, i.e., shocks and rarefaction waves has been extensively studied, there are few stability results on the contact discontinuity because unlike shock waves and rarefaction waves, its derivative has no definite sign, and decays slower than a rarefaction wave. Moreover, it behaves like a linear wave in a nonlinear setting so that its coupling with other nonlinear waves reveals a complicated interaction mechanism. Based on the new definition of contact waves to the Boltzmann equation corresponding to the contact discontinuities for the Euler equations, we succeed in obtaining the time asymptotic stability of this wave pattern with a convergence rate. In our analysis, an intrinsic dissipative mechanism associated with this profile is found and used for closing the energy estimates.     

A lattice Boltzmann model for Maxwell’s equations

In this paper, a lattice Boltzmann model for the Maxwell’s equations is proposed by taking separate sets of distribution functions for the electric and magnetic fields, and a lattice Boltzmann model for the Maxwell vorticity equations with third order truncation error is proposed by using the higher-order moment method. At the same time, the expressions of the equilibrium distribution function and the stability conditions for this model are given. As numerical examples, some classical electromagnetic phenomena, such as the electric and magnetic fields around a line current source, the electric field and equipotential lines around an electrostatic dipole, the electric and magnetic fields around oscillating dipoles are given. These numerical results agree well with classical ones.     

Macroscopic regularity for the Boltzmann equation

The regularity of solutions to the Boltzmann equation is a fundamental problem in the kinetic theory. In this paper, the case with angular cut-off is investigated. It is shown that the macroscopic parts of solutions to the Boltzmann equation, i.e., the density, momentum and total energy are continuous functions of(x, t) in the region R3×(0, +∞). More precisely, these macroscopic quantities immediately become continuous in any positive time even though they are initially discontinuous and the discontinuities of solutions propagate only in the microscopic level. It should be noted that such kind of phenomenon can not happen for the compressible Navier-Stokes equations in which the initial discontinuities of the density never vanish in any finite time, see [22]. This hints that the Boltzmann equation has better regularity effect in the macroscopic level than compressible Navier-Stokes equations.     

Uncertainty principle and kinetic equations

A large number of mathematical studies on the Boltzmann equation are based on the Grad's angular cutoff assumption. However, for particle interaction with inverse power law potentials, the associated cross-sections have a non-integrable singularity corresponding to the grazing collisions. Smoothing properties of solutions are then expected. On the other hand, the uncertainty principle, established by Heisenberg in 1927, has been developed so far in various situations, and it has been applied to the study of the existence and smoothness of solutions to partial differential equations. This paper is the first one to apply this celebrated principle to the study of the singularity in the cross-sections for kinetic equations. Precisely, we will first prove a generalized version of the uncertainty principle and then apply it to justify rigorously the smoothing properties of solutions to some kinetic equations. In particular, we give some estimates on the regularity of solutions in Sobolev spaces w.r.t. all variables for both linearized and nonlinear space inhomogeneous Boltzmann equations without angular cutoff, as well as the linearized space inhomogeneous Landau equation.     

Conservative numerical method for solving the averaged Boltzmann equation

A method is proposed for averaging the Boltzmann kinetic equation with respect to transverse velocities. A system of two integro-differential equations for two desired functions depending only on the longitudinal velocity is derived. The gas particles are assumed to interact as absolutely hard spheres. The integrals in the equations are double. The reduction in the number of variables in the desired functions and the low multiplicity of the integrals ensure the computational efficiency of the averaged equations. A numerical method of discrete ordinates is developed that effectively solves the gas relaxation problem based on the averaged equations. The method is conservative, and the number of particles, momentum, and energy are conserved automatically at every time step.     

Decay of Dissipative Equations and Negative Sobolev Spaces

We develop a general energy method for proving the optimal time decay rates of the solutions to the dissipative equations in the whole space. Our method is applied to classical examples such as the heat equation, the compressible Navier-Stokes equations and the Boltzmann equation. In particular, the optimal decay rates of the higher-order spatial derivatives of solutions are obtained. The negative Sobolev norms are shown to be preserved along time evolution and enhance the decay rates. We use a family of scaled energy estimates with minimum derivative counts and interpolations among them without linear decay analysis.     

High order moment closure for Vlasov-Maxwell equations

The extended magnetohydrodynamic models are derived based on the moment closure of the Vlasov-Maxwell (VM) equations. We adopt the Grad type moment expansion which was firstly proposed for the Boltzmann equation. A new regularization method for the Grad’s moment system was recently proposed to achieve the globally hyperbolicity so that the local well-posedness of the moment system is attained. For the VM equations, the moment expansion of the convection term is exactly the same as that in the Boltzmann equation, thus the new developed regularization applies. The moment expansion of the electromagnetic force term in the VM equations turns out to be a linear source term, which can preserve the conservative properties of the distribution function in the VM equations perfectly.     

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     

高分辨KFVS有限体积方法及其CFD应用

1．引言文中研究三维Euler方程组的数值求解·儿1）中p,（。。,。。,。z）,p和E分别表示流体密度,流体速度矢量,压力和总能．方程组（1．1）是不封闭的,除非增加一个额外的方程一状态方程p一pk句,e表示单位质量内能．本文仅限于理想气体,此时状态方程为p一（、－1加e．队2）近H十年来,涌现了许多求解方程组（1．1）的无振荡、高分辨格式,例如TVD格式问,＊＊O格式问等,它们在一定程度上促进了航空航天和造船事业的发展．其中有一类根据双曲方程组（1．l）特征值的符号建立的迎风格式尤为突出,与中心格式相比,迎风格式的耗…     

一类非线性偏微分方程的四阶格子Boltzmann模型

通过Chapman-Enskog展开技术和多尺度分析,建立了一种新的D1Q4带修正项的四阶格子Boltzmann模型,一类非线性偏微分方程从连续的Boltzmann方程得到正确恢复．统一了KdV和Burgers等已知方程类型的格子BGK模型,还首次给出了组合KdV-Burgers,广义Burgers—Huxley等方程...     

结合气体分子动力学方法的RKDG有限元格式求解一维Euler可压方程组

In this paper, two numerical methods are developed for solving one-dimensionl compressible ELder equations by the RKDG finite element method.The schemesare obtained based on an important relation between the Boltzmann equation andthe ELder equations.The schemes have the TVD-like property under the uniform meshes.Several numerical results also present the performance of the schemes.     

Singular behavior of vacuum states for compressible fluids

In this survey paper, we will present the recent work on the study of the compressible fluids with vacuum states by illustrating its interesting and singular behavior through some systems of fluid dynamics, that is, Euler equations, Euler–Poisson equations and Navier–Stokes equations. The main concern is the well-posedness of the problem when vacuum presents and the singular behavior of the solution near the interface separating the vacuum and the gas. Furthermore, the relation of the solutions for the gas dynamics with vacuum to those of the Boltzmann equation will also be discussed. In fact, the results obtained so far for vacuum states are far from being complete and satisfactory. Therefore, this paper can only be served as an introduction to this interesting field which has many open and challenging mathematical problems. Moreover, the problems considered here are limited to the author's interest and knowledge in this area.     

Classical and quantum transport in random media

We study in this article the transport of particles in time-dependent random media, in the so-called weak coupling limit. We show the convergence of a Liouville equation to a Fokker–Planck equation. We also obtain the semi-classical limit of Schrödinger equations. This limit is described by a linear Boltzmann equation. In both cases, the ratio between a typical time scale and the scale of the media determines whether the limit diffusion and the collision process are elastic or not.     

Modifications of weighted Monte Carlo algorithms for nonlinear kinetic equations

Test problems for the nonlinear Boltzmann and Smoluchowski kinetic equations are used to analyze the efficiency of various versions of weighted importance modeling as applied to the evolution of multiparticle ensembles. For coagulation problems, a considerable gain in computational costs is achieved via the approximate importance modeling of the “free path” of the ensemble combined with the importance modeling of the index of a pair of interacting particles. A weighted modification of the modeling of the initial velocity distribution was found to be the most efficient for model solutions to the Boltzmann equation. The technique developed can be useful as applied to real-life coagulation and relaxation problems for which the model problems considered give approximate solutions.     

外部流动的Oseen耦合方法,I:Oseen耦合逼近

这篇文章考虑了具有非齐次边界条件的二维非定常外部Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程．通过将内部区域的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和外部区域的Ｏｓｅｅｎ方程相耦合,得到了Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ问题的逼近问题： Ｏｓｅｅｎ耦合问题,此外,我们证明了 Ｏｓｅｅｎ耦合方程弱解的存在性,唯一性和收敛精度．     

The Passage of Weakly Coupled Nonlinear Oscillators through Internal Resonance

An asymptotic procedure for deriving equations governing the passage of a weakly coupled nonlinear system of oscillators is discussed. The procedure avoids an inner-outer-matching technique and is valid when the small coupling and detuning parameters are arbitrary. Resonance is permitted to occur at one or several instances of time or to last for a finite length of time. Numerical results are discussed.     

Oseen Coupling Method for the Exterior Flow Part I: Oseen Coupling Approximation

Abstract In this paper, we consider the bidimensional exterior unsteady Navier-Stokes equations with nonhomogeneous boundary conditions and present an Oseen coupling problem which approximates the Navier-Stokes problem, obtained by coupling the Navier-Stokes equations in the inner region and the Oseen equations in the outer region. Moreover, we prove the existence, uniqueness and the approximate accuracy of the weak solution of the Oseen coupling equations. Project supported by NSF of China & State Major Key Project of Basic Research     

A new algebraic procedure to construct exact solutions of nonlinear differential-difference equations

This paper presents a new algebraic procedure to construct exact solutions of selected nonlinear differential-difference equations. The discrete sine-Gordon equation and differential-difference asymmetric Nizhnik-Novikov-Veselov equations are chosen as examples to illustrate the efficiency and effectiveness of the new procedure, where various types of exact travelling wave solutions for these nonlinear differential-difference equations have been constructed. It is anticipated that the new procedure can also be used to produce solutions for other nonlinear differential-difference equations.     

Convergence to rarefaction waves for the nonlinear Boltzmann equation and compressible Navier-Stokes equations

The main purpose of this paper is to study the asymptotic equivalence of the Boltzmann equation for the hard-sphere collision model to its corresponding Euler equations of compressible gas dynamics in the limit of small mean free path. When the fluid flow is a smooth rarefaction (or centered rarefaction) wave with finite strength, the corresponding Boltzmann solution exists globally in time, and the solution converges to the rarefaction wave uniformly for all time (or away from t=0) as ?→0. A decomposition of a Boltzmann solution into its macroscopic (fluid) part and microscopic (kinetic) part is adopted to rewrite the Boltzmann equation in a form of compressible Navier-Stokes equations with source terms. In this setting, the same asymptotic equivalence of the full compressible Navier-Stokes equations to its corresponding Euler equations in the limit of small viscosity and heat conductivity (depending on the viscosity) is also obtained.     

KINETIC FLUX VECTOR SPLITTING FOR THE EULER EQUATIONS WITH GENERAL PRESSURE LAWS

This paper attempts to develop kinetic flux vector splitting(KFVS)for the Euler equa-tions with general pressure laws.It is well known that the gas distribution function forthe local equilibrium state plays an important role in the construction of the gas-kineticschemes.To recover the Euler equations with a general equation of state(EOS),a newlocal equilibrium distribution is introduced with two parameters of temperature approx-imation decided uniquely by macroscopic variables.Utilizing the well-known connectionthat the Euler equations of motion are the moments of the Boltzmann equation wheneverthe velocity distribution function is a local equilibrium state,a class of high resolutionMUSCL-type KFVS schemes are presented to approximate the Euler equations of gas dy-namics with a general EOS.The schemes are finally applied to several test problems for ageneral EOS.In comparison with the exact solutions,our schemes give correct location andmore accurate resolution of discontinuities.The extension of our idea to multidimensionalcase is natural.