Spectral-Lagrangian methods for collisional models of non-equilibrium statistical states

We propose a new spectral Lagrangian based deterministic solver for the non-linear Boltzmann transport equation (BTE) in d-dimensions for variable hard sphere (VHS) collision kernels with conservative or non-conservative binary interactions. The method is based on symmetries of the Fourier transform of the collision integral, where the complexity in its computation is reduced to a separate integral over the unit sphere S^d^-^1. The conservation of moments is enforced by Lagrangian constraints. The resulting scheme, implemented in free space, is very versatile and adjusts in a very simple manner to several cases that involve energy dissipation due to local micro-reversibility (inelastic interactions) or elastic models of slowing down process. Our simulations are benchmarked with available exact self-similar solutions, exact moment equations and analytical estimates for the homogeneous Boltzmann equation, both for elastic and inelastic VHS interactions. Benchmarking of the simulations involves the selection of a time self-similar rescaling of the numerical distribution function which is performed using the continuous spectrum of the equation for Maxwell molecules as studied first in Bobylev et al. [A.V. Bobylev, C. Cercignani, G. Toscani, Proof of an asymptotic property of self-similar solutions of the Boltzmann equation for granular materials, Journal of Statistical Physics 111 (2003) 403-417] and generalized to a wide range of related models in Bobylev et al. [A.V. Bobylev, C. Cercignani, I.M. Gamba, On the self-similar asymptotics for generalized non-linear kinetic Maxwell models, Communication in Mathematical Physics, in press. URL: ]. The method also produces accurate results in the case of inelastic diffusive Boltzmann equations for hard spheres (inelastic collisions under thermal bath), where overpopulated non-Gaussian exponential tails have been conjectured in computations by stochastic methods [T.V. Noije, M. Ernst, Velocity distributions in homogeneously cooling and heated granular fluids, Granular Matter 1(57) (1998); M.H. Ernst, R. Brito, Scaling solutions of inelastic Boltzmann equations with over-populated high energy tails, Journal of Statistical Physics 109 (2002) 407-432; S.J. Moon, M.D. Shattuck, J. Swift, Velocity distributions and correlations in homogeneously heated granular media, Physical Review E 64 (2001) 031303; I.M. Gamba, S. Rjasanow, W. Wagner, Direct simulation of the uniformly heated granular Boltzmann equation, Mathematical and Computer Modelling 42 (2005) 683-700] and rigorously proven in Gamba et al. [I.M. Gamba, V. Panferov, C. Villani, On the Boltzmann equation for diffusively excited granular media, Communications in Mathematical Physics 246 (2004) 503-541(39)] and [A.V. Bobylev, I.M. Gamba, V. Panferov, Moment inequalities and high-energy tails for Boltzmann equations with inelastic interactions, Journal of Statistical Physics 116 (2004) 1651-1682].

[1]  W. Wagner.  A convergence proof for Bird's direct simulation Monte Carlo method for the Boltzmann equation , 1992 .

[2]  Leslie Greengard,et al.  Spectral Approximation of the Free-Space Heat Kernel , 2000 .

[3]  Lorenzo Pareschi,et al.  Numerical Solution of the Boltzmann Equation I: Spectrally Accurate Approximation of the Collision Operator , 2000, SIAM J. Numer. Anal..

[4]  Katsuhisa Koura,et al.  Comment on “Direct Simulation Scheme Derived from the Boltzmann Equation. I. Monocomponent Gases” , 1981 .

[5]  Irene M. Gamba,et al.  On Some Properties of Kinetic and Hydrodynamic Equations for Inelastic Interactions , 2000 .

[6]  G. Toscani,et al.  Self-Similarity and Power-Like Tails in Nonconservative Kinetic Models , 2010, 1009.2760.

[7]  C. Villani,et al.  Upper Maxwellian Bounds for the Spatially Homogeneous Boltzmann Equation , 2007 .

[8]  Lorenzo Pareschi,et al.  A Fourier spectral method for homogeneous boltzmann equations , 1996 .

[9]  J. Broadwell,et al.  Study of rarefied shear flow by the discrete velocity method , 1964, Journal of Fluid Mechanics.

[10]  Carlo Cercignani,et al.  Shear Flow of a Granular Material , 2001 .

[11]  Steven G. Johnson,et al.  The Fastest Fourier Transform in the West , 1997 .

[12]  M D Shattuck,et al.  Velocity distributions and correlations in homogeneously heated granular media. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  S. Rjasanow,et al.  Difference scheme for the Boltzmann equation based on the Fast Fourier Transform , 1997 .

[14]  Sergej Rjasanow,et al.  Numerical solution of the Boltzmann Equation using fully conservative difference scheme based on the Fast Fourier Transform , 2000 .

[15]  Lorenzo Pareschi,et al.  Fast algorithms for computing the Boltzmann collision operator , 2006, Math. Comput..

[16]  C. Cercignani,et al.  On the Self-Similar Asymptotics for Generalized Nonlinear Kinetic Maxwell Models , 2006, math-ph/0608035.

[17]  Alexander V. Bobylev,et al.  The inverse laplace transform of some analytic functions with an application to the eternal solutions of the Boltzmann equation , 2002, Appl. Math. Lett..

[18]  Mackintosh,et al.  Driven granular media in one dimension: Correlations and equation of state. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[19]  G. Toscani,et al.  Relaxation Schemes for Nonlinear Kinetic Equations , 1997 .

[20]  Product measure steady states of generalized zero range processes , 2005, cond-mat/0506776.

[21]  Giuseppe Toscani,et al.  Proof of an Asymptotic Property of Self-Similar Solutions of the Boltzmann Equation for Granular Materials , 2003 .

[22]  Stéphane Mischler,et al.  A spatially homogeneous Boltzmann equation for elastic, inelastic and coalescing collisions , 2005 .

[23]  Carlo Cercignani,et al.  Exact Eternal Solutions of the Boltzmann Equation , 2002 .

[24]  Henri Cabannes.  Global solution of the initial value problem for the discrete Boltzmann equation , 1978 .

[25]  A. Bobylev,et al.  Exact solutions of the nonlinear Boltzmann equation and the theory of relaxation of a Maxwellian gas , 1984 .

[26]  Irene M. Gamba,et al.  On the Boltzmann Equation for Diffusively Excited Granular Media , 2004 .

[27]  J. Kestin,et al.  Handbook of fluid dynamics , 1948 .

[28]  Nicolas,et al.  A Boltzmann equation for elastic, inelastic and coalescing collisions , 2003 .

[29]  R Brito,et al.  Driven inelastic Maxwell models with high energy tails. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[30]  T. Wu,et al.  Formation of Maxwellian Tails , 1976 .

[31]  Carlo Cercignani,et al.  Moment Equations for a Granular Material in a Thermal Bath , 2002 .

[32]  Francis Filbet,et al.  High order numerical methods for the space non-homogeneous Boltzmann equation , 2003 .

[33]  W. Wagner,et al.  Stochastic Numerics for the Boltzmann Equation , 2005 .

[34]  Yoshio Sone,et al.  Molecular gas dynamics , 2007 .

[35]  L. Mieussens.  Discrete-Velocity Models and Numerical Schemes for the Boltzmann-BGK Equation in Plane and Axisymmetric Geometries , 2000 .

[36]  Carlo Cercignani,et al.  Self-Similar Asymptotics for the Boltzmann Equation with Inelastic and Elastic Interactions , 2003 .

[37]  S. Rjasanow,et al.  Fast deterministic method of solving the Boltzmann equation for hard spheres , 1999 .

[38]  W. Wagner,et al.  A Stochastic Weighted Particle Method for the Boltzmann Equation , 1996 .

[39]  J. M. Montanero,et al.  Computer simulation of uniformly heated granular fluids , 2000 .

[40]  L. C. Pitchford,et al.  A Numerical Solution of the Boltzmann Equation , 1983 .

[41]  Mohammed Seaïd,et al.  Discrete-Velocity Models and Relaxation Schemes for Traffic Flows , 2006, SIAM J. Sci. Comput..

[42]  Irene M. Gamba,et al.  Moment Inequalities and High-Energy Tails for Boltzmann Equations with Inelastic Interactions , 2004 .

[43]  Sergej Rjasanow,et al.  Numerical solution of the Boltzmann equation on the uniform grid , 2002, Computing.

[44]  Henning Struchtrup,et al.  A linearization of Mieussens's discrete velocity model for kinetic equations , 2007 .

[45]  Lorenzo Pareschi,et al.  Solving the Boltzmann Equation in N log2N , 2006, SIAM J. Sci. Comput..

[46]  Shuichi Kawashima.  Global solution of the initial value problem for a discrete velocity model of the Boltzmann equation , 1981 .

[47]  Sergej Rjasanow,et al.  Direct simulation of the uniformly heated granular boltzmann equation, , 2005, Math. Comput. Model..

[48]  Carlo Cercignani,et al.  Shock Waves for a Discrete Velocity Gas Mixture , 2000 .

[49]  Maria Groppi,et al.  Approximate solutions to the problem of stationary shear flow of smooth granular materials , 2002 .

[50]  A. Bobylev.  Exact solutions of the Boltzmann equation , 1975 .

[51]  Irene M. Gamba,et al.  Propagation of L1 and L∞ Maxwellian weighted bounds for derivatives of solutions to the homogeneous elastic Boltzmann equation , 2007 .

[52]  M. H. Ernst,et al.  Velocity Distributions in Homogeneously Cooling and Heated Granular Fluids , 1998 .

[53]  Carlo Cercignani,et al.  Discrete Velocity Models Without Nonphysical Invariants , 1999 .

[54]  A. V. Bobylev,et al.  Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions and Power Like Tails , 2006 .

[55]  M. H. Ernst,et al.  Scaling Solutions of Inelastic Boltzmann Equations with Over-Populated High Energy Tails , 2002 .