Fractional Fokker-Planck equations for subdiffusion with space- and time-dependent forces

We derive a fractional Fokker-Planck equation for subdiffusion in a general space- and time-dependent force field from power law waiting time continuous time random walks biased by Boltzmann weights. The governing equation is derived from a generalized master equation and is shown to be equivalent to a subordinated stochastic Langevin equation.     

Convergence from Boltzmann to Landau Processes with Soft Potential and Particle Approximations

Our aim in this paper is to show how a probabilistic interpretation of the Boltzmann and Landau equations gives a microscopic understanding of these equations. We firstly associate stochastic jump processes with the Boltzmann equations we consider. Then we renormalize these equations following asymptotics which make prevail the grazing collisions, and prove the convergence of the associated Boltzmann jump processes to a diffusion process related to the Landau equation. The convergence is pathwise and also implies a convergence at the level of the partial differential equations. The best feature of this approach is the microscopic understanding of the transition between the Boltzmann and the Landau equations, by an accumulation of very small jumps. We deduce from this interpretation an approximation result for a solution of the Landau equation via colliding stochastic particle systems. This result leads to a Monte-Carlo algorithm for the simulation of solutions by a conservative particle method which enables to observe the transition from Boltzmann to Landau equations. Numerical results are given.     

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 …     

Dissipative dynamics for hard spheres

The dynamics for a system of hard spheres with dissipative collisions is described at the levels of statistical mechanics, kinetic theory, and simulation. The Liouville operator(s) and associated binary scattering operators are defined as the generators for time evolution in phase space. The BBGKY hierarchy for reduced distribution functions is given, and an approximate kinetic equation is obtained that extends the revised Enskog theory to dissipative dynamics. A Monte Carlo simulation method to solve this equation is described, extending the Bird method to the dense, dissipative hard-sphere system. A practical kinetic model for theoretical analysis of this equation also is proposed. As an illustration of these results, the kinetic theory and the Monte Carlo simulations are applied to the homogeneous cooling state of rapid granular flow.     

A discontinuous Galerkin finite-element method for a 1D prototype of the Boltzmann equation

To develop and analyze new computational techniques for the Boltzmann equation based on model or approximation adaptivity, it is imperative to have disposal of a compliant model problem that displays the essential characteristics of the Boltzmann equation and that admits the extraction of highly accurate reference solutions. For standard collision processes, the Boltzmann equation itself fails to meet the second requirement for d = 2, 3 spatial dimensions, on account of its setting in 2d, while for d = 1 the first requirement is violated because the Boltzmann equation then lacks the convergence-to-equilibrium property that forms the substructure of simplified models. In this article we present a numerical investigation of a new one-dimensional prototype of the Boltzmann equation. The underlying molecular model is endowed with random collisions, which have been fabricated such that the corresponding Boltzmann equation exhibits convergence to Maxwell–Boltzmann equilibrium solutions. The new Boltzmann model is approximated by means of a discontinuous Galerkin (DG) finite-element method. To validate the one-dimensional Boltzmann model, we conduct numerical experiments and compare the results with Monte-Carlo simulations of equivalent molecular-dynamics models.     

Quantum Lattice-Gas Model for the Burgers Equation

A quantum algorithm is presented for modeling the time evolution of a continuous field governed by the nonlinear Burgers equation in one spatial dimension. It is a microscopic-scale algorithm for a type-II quantum computer, a large lattice of small quantum computers interconnected in nearest neighbor fashion by classical communication channels. A formula for quantum state preparation is presented. The unitary evolution is governed by a conservative quantum gate applied to each node of the lattice independently. Following each quantum gate operation, ensemble measurements over independent microscopic realizations are made resulting in a finite-difference Boltzmann equation at the mesoscopic scale. The measured values are then used to re-prepare the quantum state and one time step is completed. The procedure of state preparation, quantum gate application, and ensemble measurement is continued ad infinitum. The Burgers equation is derived as an effective field theory governing the behavior of the quantum computer at its macroscopic scale where both the lattice cell size and the time step interval become infinitesimal. A numerical simulation of shock formation is carried out and agrees with the exact analytical solution.     

基于莱维飞行的鸟群优化算法

针对鸟群优化算法(BSA)在求解高维多极值优化问题时容易陷入局部最优解和出现早熟收敛的情况,在原始鸟群算法的基础上,在模拟鸟群飞行行为的过程中引入莱维飞行,提出了一种基于莱维飞行的改进算法——莱维-鸟群算法(LBSA)。这种算法替换了原算法中随机的飞行位置跳变,而采用莱维飞行更新鸟群飞行后的位置,大幅提高了鸟群的位置变化活力,提高了算法的有效性。仿真结果表明,在求解高维多极值优化问题时,该算法性能优于原始鸟群算法。     

Lattice Boltzmann method on quadtree grids for simulating fluid flow through porous media: A new automatic algorithm

During the past two decades, the lattice Boltzmann (LB) method has been introduced as a class of computational fluid dynamic methods for fluid flow simulations. In this method, instead of solving the Navier Stocks equation, the Boltzmann equation is solved to simulate the flow of a fluid. This method was originally developed based on uniform grids. However, in order to model complex geometries such as porous media, it can be very slow in comparison with other techniques such as finite differences and finite elements. To eliminate this limitation, a number of studies have aimed to formulate the lattice Boltzmann on the unstructured grids. This paper deals with simulating fluid flow through a synthetic porous medium using the LB method and on the quadtree grid structure. To this end, the LB method was used on nonuniform grids coupled with a technique for image reconstruction which resulted in the quadtree grids for simulation of fluid flow through porous media. Accuracy and efficiency of this algorithm is compared against the conventional LB method based on uniform grids. While the decrease in computational time in the proposed LB method on nonuniform grids is found to be significant regarding the size of the initial and reconstructed images, the same level of accuracy is obtained when compared with the conventional LB method on uniform grids.     

An exact energy conservation property of the quantum lattice Boltzmann algorithm

The quantum lattice Boltzmann algorithm offers a unitary and readily parallelisable discretisation of the Dirac equation that is free of the fermion-doubling problem. The expectation of the discrete time-advance operator is an exact invariant of the algorithm. Its imaginary part determines the expectation of the Hamiltonian operator, the energy of the solution, with an accuracy that is consistent with the overall accuracy of the algorithm. In the one-dimensional case, this accuracy may be increased from first to second order using a variable transformation. The three-dimensional quantum lattice Boltzmann algorithm uses operator splitting to approximate evolution under the three-dimensional Dirac equation by a sequence of solutions of one-dimensional Dirac equations. The three-dimensional algorithm thus inherits the energy conservation property of the one-dimensional algorithm, although the implementation shown remains only first-order accurate due to the splitting error.     

Long Time Energy Transfer in the Random Schrödinger Equation

We consider the long time behavior of solutions of the d-dimensional linear Boltzmann equation that arises in the weak coupling limit for the Schrödinger equation with a time-dependent random potential. We show that the intermediate mesoscopic time limit satisfies a Fokker–Planck type equation with the wave vector performing a Brownian motion on the (d ? 1)-dimensional sphere of constant energy, as in the case of a time-independent Schrödinger equation. However, the long time limit of the solution with an isotropic initial data satisfies an equation corresponding to the energy being the square root of a Bessel process of dimension d/2.     

A high-order Monte Carlo algorithm for the direct simulation of Boltzmann equation

A high-order algorithm of the direct simulation Monte Carlo (DSMC) method, H-DSMC, has been developed to simulate rarefied flow regimes. The mth order Taylor series expansion has been employed to obtain a more generalized form of the time discretization for the collision part of Boltzmann equation. In the purposed algorithm, the higher order collision terms are introduced as well as higher order terms in the time step of the probabilistic coefficients. These newly implemented higher order terms improve the accuracy and efficiency of the solution and enhance the convergence rate quite significantly. Comparison between results of the classic DSMC method and the H-DSMC method shows the promising performance of the introduced technique.     

气体运动论统一算法在跨流域转动非平衡效应模拟中的应用

从考虑转动松弛变化特性的Rykov模型出发,基于Boltzmann模型方程求解跨流域气体运动论统一算法原理与计算规则,采用转动惯量描述气体分子自旋运动,研究考虑转动非平衡影响的Boltzmann模型方程数值求解方法.通过对氮气激波结构、二维钝头体和三维尖双锥跨流域绕流的模拟分析,验证该算法的跨流域一致适用性.     

Boltzmann Limit and Quasifreeness for a Homogenous Fermi Gas in a Weakly Disordered Random Medium

We discuss some basic aspects of the dynamics of a homogenous Fermi gas in a weak random potential, under negligence of the particle pair interactions. We derive the kinetic scaling limit for the momentum distribution function with a translation invariant initial state and prove that it is determined by a linear Boltzmann equation. Moreover, we prove that if the initial state is quasifree, then the time evolved state, averaged over the randomness, has a quasifree kinetic limit. We show that the momentum distributions determined by the Gibbs states of a free fermion field are stationary solutions of the linear Boltzmann equation; this includes the limit of zero temperature.     

Simulation of euclidean quantum field theories by a random walk process

A new Monte Carlo method for euclidean lattice field theory is introduced by writing the Boltzmann distribution e^?s as a solution of a diffusion type equation and constructing the associated random walk process. It is practically tested for a quantum mechanical model and a non-compact version of lattice QCD. It is explained where the main interest in this algorithm lies: the diffusion process coming from an action that can be generalized to include non-conservative forces. This possibility is exploited in our QCD version to implement gauge fixing without Faddeev-Popov ghosts.     

Development of Lattice Boltzmann Flux Solver for Simulation of Incompressible Flows

A lattice Boltzmann flux solver (LBFS) is presented in this work for simulation of incompressible viscous and inviscid flows. The new solver is based on Chapman-Enskog expansion analysis, which is the bridge to link Navier-Stokes (N-S) equations and lattice Boltzmann equation (LBE). The macroscopic differential equations are discretized by the finite volume method, where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers. The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. LBFS is validated by its application to simulate the viscous decaying vortex flow, the driven cavity flow, the viscous flow past a circular cylinder, and the inviscid flow past a circular cylinder. The obtained numerical results compare very well with available data in the literature, which show that LBFS has the second order of accuracy in space, and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.     

含振动能激发Boltzmann模型方程气体动理论统一算法验证与分析

为模拟研究高温高马赫数下多原子气体内能激发对跨流域非平衡流动的影响,将转动能、振动能分别作为气体分子速度分布函数的自变量,把转动能和振动能处理为连续分布的能量模式,将Boltzmann方程的碰撞项分解成弹性碰撞项和非弹性碰撞项,同时将非弹性碰撞按一定松弛速率分解为平动-转动能松弛过程和平动-转动-振动能松弛过程,构造了一类考虑振动能激发的Boltzmann模型方程,并证明了其守恒性和H定理.基于内部能量变量对分布函数无穷积分,引入三个约化速度分布函数,得到一组考虑振动能激发的约化速度分布函数控制方程组,使用离散速度坐标法,基于LU-SGS隐式格式和有限体积法求解离散速度分布函数,建立含振动能激发的气体动理论统一算法.通过开展高稀薄流到连续流圆柱绕流问题统一算法与直接模拟蒙特卡罗法模拟结果对比分析,特别是过渡流区平动、转动、振动非平衡效应对绕流流场与物面力热特性的影响机制,证实了所建立的含振动能激发的Boltzmann模型方程及气体动理论统一算法的准确可靠性.     

A Simplified Linearized Lattice Boltzmann Method for Acoustic Propagation Simulation

A simplified linearized lattice Boltzmann method (SLLBM) suitable for the simulation of acoustic waves propagation in fluids was proposed herein. Through Chapman–Enskog expansion analysis, the linearized lattice Boltzmann equation (LLBE) was first recovered to linearized macroscopic equations. Then, using the fractional-step calculation technique, the solution of these linearized equations was divided into two steps: a predictor step and corrector step. Next, the evolution of the perturbation distribution function was transformed into the evolution of the perturbation equilibrium distribution function using second-order interpolation approximation of the latter at other positions and times to represent the nonequilibrium part of the former; additionally, the calculation formulas of SLLBM were deduced. SLLBM inherits the advantages of the linearized lattice Boltzmann method (LLBM), calculating acoustic disturbance and the mean flow separately so that macroscopic variables of the mean flow do not affect the calculation of acoustic disturbance. At the same time, it has other advantages: the calculation process is simpler, and the cost of computing memory is reduced. In addition, to simulate the acoustic scattering problem caused by the acoustic waves encountering objects, the immersed boundary method (IBM) and SLLBM were further combined so that the method can simulate the influence of complex geometries. Several cases were used to validate the feasibility of SLLBM for simulation of acoustic wave propagation under the mean flow.     

Generalized Boltzmann equation for lattice gas automata

In this paper a theory is formulated that predicts velocity and spatial correlations between occupation numbers that occur in lattice gas automata violating semi-detailed balance. Starting from a coupled BBGKY hierarchy for then-particle distribution functions, cluster expansion techniques are used to derive approximate kinetic equations. In zeroth approximation the standard nonlnear Boltzmann equation is obtained; the next approximation yields the ring kinetic equation, similar to that for hard-sphere systems, describing the time evolution of pair correlations. The ring equation is solved to determine the (nonvanishing) pair correlation functions in equilibrium for two models that violate semidetailed balance. One is a model of interacting random walkers on a line, the other one is a two-dimensional fluid-type model on a triangular lattice. The numerical predictions agree very well with computer simulations.     

Cooling Process for Inelastic Boltzmann Equations for Hard Spheres, Part I: The Cauchy Problem

We develop the Cauchy theory of the spatially homogeneous inelastic Boltzmann equation for hard spheres, for a general form of collision rate which includes in particular variable restitution coefficients depending on the kinetic energy and the relative velocity as well as the sticky particles model. We prove (local in time) non-concentration estimates in Orlicz spaces, from which we deduce weak stability and existence theorem. Strong stability together with uniqueness and instantaneous appearance of exponential moments are proved under additional smoothness assumption on the initial datum, for a restricted class of collision rates. Concerning the long-time behaviour, we give conditions for the cooling process to occur or not in finite time. Mathematics Subject Classification (2000): 76P05 Rarefied gas flows, Boltzmann equation [See also 82B40, 82C40, 82D05].     

小波分析在层析图像重构中的应用研究

小波分析作为一种非平稳信号分析方法,具有良好的时（ 空）、频局部化特性和多分辨率特性。介绍了小波分析的基本原理和应用,引入小波分析进行图像重构,利用小波分解后得到的多分辨率的稀疏矩阵表示,设计了一迭代重构算法。通过计算机仿真试验,验证了小波分析在计算层析（ＣＴ） 成像光谱技术中能够应用于图像重构,并证实了小波分析的迭代重构算法是稳定、多分辨率和快速收敛的。