the green s function and large time behavior of solutions for the one dimensional boltzmann equation we study the pointwise behavior of the green s function of the boltzmann equation our results reveal the particle and fluid aspects of the equation the particle aspect is represented by singular waves these waves are carried by transport equations and dominate the short time behavior of the solution we devise a picard type iteration for constructing the increasingly regular particlelike waves the fluidlike waves reveal the dissipative behavior of the type of navier stokes equations as usually seen by the chapman enskog expansion these waves are constructed as part of the long wave expansion in the spectrum of the fourier mode for the space variable the fluidlike waves represent the long time behavior of the solution as an application we obtain the pointwise description of the large time behavior of the convergence to the global maxwellian when the initial perturbation is not necessarily smooth in our analysis of the exchanges of the microscopic velocity decay and space decay we make essential use of the hard sphere models c 2004 wiley periodicals inc contents

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

wavelet like bases for the fast solutions of second kind integral equations a class of vector space bases is introduced for the sparse representation of discretizations of integral operators an operator with a smooth nonoscillatory kernel possessing a finite number of singularities in each row or column is represented in these bases as a sparse matrix to high precision a method is presented that employs these bases for the numerical solution of second kind integral equations in time bounded by o n log 2 n where n is the number of points in the discretization numerical results are given which demonstrate the effectiveness of the approach and several generalizations and applications of the method are discussed

a maxwellian lower bound for solutions to the boltzmann equation we prove that the solution of the spatially homogeneous boltzmann equation is bounded pointwise from below by a maxwellian i e a function of the formc1 exp c2v2 this holds for any initial data with bounded mass energy and entropy and for any positive timet t0 the constantsc1 andc2 depend on the mass energy and entropy of the initial data and ont0 0 only

general mapping discrete time models of a descriptor system with an arbitrary initial condition abstract the definition of discretization proposed previously by the authors for continuous time signals is extended to include an impulse and its derivatives and then applied to descriptor systems that produce such signals this definition provides bases for the continuous time impulses to have their discrete time counterparts by taking values at multiple sampling points into account rather than a single point that has traditionally been used for ordinary functions theorems are developed based on the definition one of which provides a useful formula for choosing a sufficient initial condition for a discrete time descriptor system such that its initial response corresponds in the sense of the proposed discretization definition to the continuous time original this solves the long standing issue of ambiguity surrounding discontinuous initial responses and that of inconsistent initial conditions that have obscured the validity of discrete time descriptor models and of their numerical computations

conservative numerical schemes for the vlasov equation abstract a new scheme for solving the vlasov equation using a phase space grid is proposed the algorithm is based on the conservation of the flux of particles and the distribution function is reconstructed using various techniques that allow control of spurious oscillations or preservation of the positivity several numerical results are presented in two and four dimensional phase space and the scheme is compared with the semilagrangian method this method is almost as accurate as the semilagrangian one and the local reconstruction technique is well suited for parallel computation

on the trend to global equilibrium for spatially inhomogeneous kinetic systems the boltzmann equation as part of our study of convergence to equilibrium for spatially inhomogeneous kinetic equations started in 21 we derive estimates on the rate of convergence to equilibrium for solutions of the boltzmann equation like o t our results hold conditionally to some strong but natural estimates of smoothness decay at large velocities and strict positivity which at the moment have only been established in certain particular cases among the most important steps in our proof are 1 quantitative variants of boltzmann s h theorem as proven in 52 60 based on symmetry features hypercontractivity and information theoretical tools 2 a new quantitative version of the instability of the hydrodynamic description for non small knudsen number 3 some functional inequalities with geometrical content in particular the korn type inequality which we established in 22 and 4 the study of a system of coupled differential inequalities of second order by a treatment inspired from 21 we also briefly point out the particular role of conformal velocity fields when they are allowed by the geometry of the problem

a parametrized sum of fuzzy numbers with applications to fuzzy initial value problems abstract this paper presents a new parametrized family of joint possibility distributions that can be used to control the pompeiu hausdorff norm of the sum resulting from the extension principle we are able to generate sums with norms covering the entire interval ranging from the smallest possible value to the largest one finally we apply our approach to numerical methods for solving first order fuzzy initial value problems

a consistency result for a discrete velocity model of the boltzmann equation in this paper we study the link between a certain class of discrete velocity models dvms and the boltzmann equation those models possess an infinite number of velocities laying on a regular grid of step h in r3 our aim is to prove that it is possible to construct models consistent with the boltzmann equation i e such that the discrete collision term can be seen as an approximation of the collision integral of the boltzmann equation

a numerical method for the accurate solution of the fokker planck landau equation in the nonhomogeneous case abstract a new approach for the accurate numerical solution of the fokker planck landau fpl equation in the nonhomogeneous case is presented the method couples through a time splitting algorithm a finite volume scheme for the transport with a fast spectral solver for the efficient solution of the collision operator recently introduced the scheme allows the use of different grids in the velocity space for the transport and the collision phases the use of a suitable explicit runge kutta solver for the time integration of the collision phase permits avoid once of excessive small time steps induced by the stiffness of the diffusive collision operator numerical results for both space homogeneous and space nonhomogeneous situations show the efficiency and the accuracy of the present method

the vlasov maxwell boltzmann system near maxwellians perhaps the most fundamental model for dynamics of dilute charged particles is described by the vlasov maxwell boltzmann system in which particles interact with themselves through collisions and with their self consistent electromagnetic field despite its importance no global in time solutions weak or strong have been constructed so far it is shown in this article that any initially smooth periodic small perturbation of a given global maxwellian which preserves the same mass total momentum and reduced total energy 22 leads to a unique global in time classical solution for such a master system the construction is based on a recent nonlinear energy method with a new a priori estimate for the dissipation the linear collision operator l not its time integration is positive definite for any solutionf t x v with small amplitude to the vlasov maxwell boltzmann system 8 and 12 as a by product such an estimate also yields an exponential decay for the simpler vlasov poisson boltzmann system 24

brownian approximation and monte carlo simulation of the non cutoff kac equation the non cutoff boltzmann equation can be simulated using the dsmc method by a truncation of the collision term however even for computing stationary solutions this may be very time consuming in particular in situations far from equilibrium by adding an appropriate diffusion to the dsmc method the rate of convergence when the truncation is removed may be greatly improved we illustrate the technique on a toy model the kac equation as well as on the full boltzmann equation in a special case

solving the boltzmann equation in n log 2 n in c mouhot and l pareschi fast algorithms for computing the boltzmann collision operator math comp to appear c mouhot and l pareschi c r math acad sci paris 339 2004 pp 71 76 fast deterministic algorithms based on spectral methods were derived for the boltzmann collision operator for a class of interactions including the hard spheres model in three dimensions these algorithms are implemented for the solution of the boltzmann equation in two and three dimensions first for homogeneous solutions then for general nonhomogeneous solutions the results are compared to explicit solutions when available and to monte carlo methods in particular the computational cost and accuracy are compared to those of monte carlo methods as well as to those of previous spectral methods finally for inhomogeneous solutions we take advantage of the great computational efficiency of the method to show an oscillation phenomenon of the entropy functional in the trend to equilibrium which was suggested in the work l desvillettes and c villani invent math 159 2005 pp 245 316

a class of asymptotic preserving schemes for kinetic equations and related problems with stiff sources in this paper we propose a general time discrete framework to design asymptotic preserving schemes for initial value problem of the boltzmann kinetic and related equations numerically solving these equations are challenging due to the nonlinear stiff collision source terms induced by small mean free or relaxation time we propose to penalize the nonlinear collision term by a bgk type relaxation term which can be solved explicitly even if discretized implicitly in time moreover the bgk type relaxation operator helps to drive the density distribution toward the local maxwellian thus naturally imposes an asymptotic preserving scheme in the euler limit the scheme so designed does not need any nonlinear iterative solver or the use of wild sum it is uniformly stable in terms of the possibly small knudsen number and can capture the macroscopic fluid dynamic euler limit even if the small scale determined by the knudsen number is not numerically resolved it is also consistent to the compressible navier stokes equations if the viscosity and heat conductivity are numerically resolved the method is applicable to many other related problems such as hyperbolic systems with stiff relaxation and high order parabolic equations

a successive penalty based asymptotic preserving scheme for kinetic equations we propose an asymptotic preserving ap scheme for kinetic equations that is efficient also in the hydrodynamic regimes this scheme is based on the bhantnagar gross krook bgk penalty method introduced by filbet and jin j comput phys 229 2010 pp 7625 7648 but uses the penalization successively to achieve the desired asymptotic property this method possesses a stronger ap property than the original method of filbet and jin with the additional feature of being positivity preserving when applied on the boltzmann equation it is also general enough to be applicable to several important classes of kinetic equations including the boltzmann equation and the landau equation numerical experiments verify these properties

a stochastic galerkin method for the boltzmann equation with uncertainty we develop a stochastic galerkin method for the boltzmann equation with uncertainty the method is based on the generalized polynomial chaos gpc approximation in the stochastic galerkin framework and can handle random inputs from collision kernel initial data or boundary data we show that a simple singular value decomposition of gpc related coefficients combined with the fast fourier spectral method in velocity space allows one to compute the high dimensional collision operator very efficiently in the spatially homogeneous case we first prove that the analytical solution preserves the regularity of the initial data in the random space and then use it to establish the spectral accuracy of the proposed stochastic galerkin method several numerical examples are presented to illustrate the validity of the proposed scheme

a hybrid method for hydrodynamic kinetic flow part i a particle grid method for reducing stochastic noise in kinetic regimes in this work we present a hybrid particle grid monte carlo method for the boltzmann equation which is characterized by a significant reduction of the stochastic noise in the kinetic regime the hybrid method is based on a first order splitting in time to separate the transport from the relaxation step the transport step is solved by a deterministic scheme while a hybrid dsmc based method is used to solve the collision step such a hybrid scheme is based on splitting the solution in a collisional and a non collisional part at the beginning of the collision step and the dsmc method is used to solve the relaxation step for the collisional part of the solution only this is accomplished by sampling only the fraction of particles candidate for collisions from the collisional part of the solution performing collisions as in a standard dsmc method and then projecting the particles back onto a velocity grid to compute a piecewise constant reconstruction for the collisional part of the solution the latter is added to a piecewise constant reconstruction of the non collisional part of the solution which in fact remains unchanged during the relaxation step numerical results show that the stochastic noise is significantly reduced at large knudsen numbers with respect to the standard dsmc method indeed in this algorithm the particle scheme is applied only on the collisional part of the solution so only this fraction of the solution is affected by stochastic fluctuations but since the collisional part of the solution reduces as the knudsen number increases stochastic noise reduces as well at large knudsen numbers

fast algorithms for computing the boltzmann collision operator the development of accurate and fast numerical schemes for the five fold boltzmann collision integral represents a challenging problem in scientific computing for a particular class of interactions including the so called hard spheres model in dimension three we are able to derive spectral methods that can be evaluated through fast algorithms these algorithms are based on a suitable representation and approximation of the collision operator explicit expressions for the errors in the schemes are given and spectral accuracy is proved parallelization properties and adaptivity of the algorithms are also discussed

fast multipole method for multivariable integrals we give a fast numerical algorithm to evaluate a class of multivariable integrals a direct numerical evaluation of these integrals costs n m where m is the number of variables and n is the number of the quadrature points for each variable for m 2 and m 3 and for only one dimensional variables we present an algorithm which is able to reduce this cost from n m to a cost of the order of o log epsilon mu m n where epsilon is the desired accuracy and mu m is a constant that depends only on m then we make some comments about possible extensions of such algorithms to number of variables m geq 4 and to higher dimensions this recursive algorithm can be viewed as an extension of fast multipole methods to situations where the interactions between particles are more complex than the standard case of binary interactions numerical tests illustrating the efficiency and the limitation of this method are presented

a unified gas kinetic scheme for continuum and rarefied flows iv fluid dynamic equations are valid in their respective modeling scales such as the particle mean free path scale of the boltzmann equation and the hydrodynamic scale of the navier stokes ns equations with a variation of the modeling scales theoretically there should have a continuous spectrum of fluid dynamic equations even though the boltzmann equation is claimed to be valid in all scales many boltzmann solvers including direct simulation monte carlo method require the cell resolution to the order of particle mean free path scale therefore they are still single scale methods in order to study multiscale flow evolution efficiently the dynamics in the computational fluid has to be changed with the scales a direct modeling of flow physics with a changeable scale may become an appropriate approach the unified gas kinetic scheme ugks is a direct modeling method in the mesh size scale and its underlying flow physics depends on the resolution of the cell size relative to the particle mean free path the cell size of ugks is not limited by the particle mean free path with the variation of the ratio between the numerical cell size and local particle mean free path the ugks recovers the flow dynamics from the particle transport and collision in the kinetic scale to the wave propagation in the hydrodynamic scale the previous ugks is mostly constructed from the evolution solution of kinetic model equations even though the ugks is very accurate and effective in the low transition and continuum flow regimes with the time step being much larger than the particle mean free time it still has space to develop more accurate flow solver in the region where the time step is comparable with the local particle mean free time in such a scale there is dynamic difference from the full boltzmann collision term and the model equations this work is about the further development of the ugks with the implementation of the full boltzmann collision term in the region where it is needed the central ingredient of the ugks is the coupled treatment of particle transport and collision in the flux evaluation across a cell interface where a continuous flow dynamics from kinetic to hydrodynamic scales is modeled the newly developed ugks has the asymptotic preserving ap property of recovering the ns solutions in the continuum flow regime and the full boltzmann solution in the rarefied regime in the mostly unexplored transition regime the ugks itself provides a valuable tool for the non equilibrium flow study the mathematical properties of the scheme such as stability accuracy and the asymptotic preserving will be analyzed in this paper as well

numerical solution of the boltzmann equation i spectrally accurate approximation of the collision operator in this paper we show that the use of spectral galerkin methods for the approximation of the boltzmann equation in the velocity space permits one to obtain spectrally accurate numerical solutions at a reduced computational cost we prove that the spectral algorithm preserves the total mass and approximates with infinite order accuracy momentum and energy consistency of the method is also proved and a stability result for a smoothed positive scheme is given we demonstrate that the fourier coefficients associated with the collision kernel of the equation have a very simple structure and in some cases can be computed explicitly numerical examples for homogeneous test problems in two and three dimensions confirm the advantages of the method

high order numerical methods for the space non homogeneous boltzmann equation in this paper we present accurate methods for the numerical solution of the boltzmann equation of rarefied gas the methods are based on a time splitting technique the transport is solved by a third order accurate in space positive and flux conservative pfc method the collision step is treated by a fourier approximation of the collision integral which guarantees spectral accuracy in velocity coupled with several high order integrators in time strang splitting is used to achieve second order accuracy in space and time several numerical tests illustrate the properties of the methods

relaxation schemes for nonlinear kinetic equations a class of numerical schemes for nonlinear kinetic equations of boltzmann type is described following wild s approach the solution is represented as a power series with parameter depending exponentially on the knudsen number this permits us to derive accurate and stable time discretizations for all ranges of the mean free path these schemes preserve the main physical properties positivity conservation of mass momentum and energy moreover for some particular models the entropy property is also shown to hold

accurate numerical methods for the collisional motion of heated granular flows in this paper we extend the spectral method developed in l pareschi b perthame a fourier spectral method for homogeneous boltzmann equations trans theo stat phys 25 1996 369 383 l pareschi g russo numerical solution of the boltzmann equation i spectrally accurate approximation of the collision operator siam j numer anal 37 2000 1217 1245 to the case of the inelastic boltzmann equation describing the collisional motion of a granular gas with and without a heating source the schemes are based on a fourier representation of the equation in the velocity space and provide a very accurate description of the time evolution of the distribution function several numerical results in dimension one to three show the efficiency and accuracy of the proposed algorithms some mathematical and physical conjectures are also addressed with the aid of the numerical simulations

existence stability and convergence of solutions of discrete velocity models to the boltzmann equation we prove the convergence of finite difference approximations to solutions of the boltzmann equation an essential step is the proof of convergence of discrete approximations to the collision integral this proof relies on our previous results on the consistency of this approximation for the space homogeneous problem we prove strong convergence of our discrete approximation to the strong solution of the boltzmann equation in the space dependent case we prove weak convergence to diperna lions solutions