A moment method with closures based on Gaussian quadrature formulas is proposed to solve the Boltzmann kinetic equation with a hard-sphere collision kernel for mono-dispersed particles. Different orders of accuracy in terms of the moments of the velocity distribution function are considered, accounting for moments up to seventh order. Quadrature-based closures for four different models for inelastic collision-the Bhatnagar-Gross-Krook, ES-BGK, the Maxwell model for hard-sphere collisions, and the full Boltzmann hard-sphere collision integral-are derived and compared. The approach is validated studying a dilute non-isothermal granular flow of inelastic particles between two stationary Maxwellian walls. Results obtained from the kinetic models are compared with the predictions of molecular dynamics (MD) simulations of a nearly equivalent system with finite-size particles. The influence of the number of quadrature nodes used to approximate the velocity distribution function on the accuracy of the predictions is assessed. Results for constitutive quantities such as the stress tensor and the heat flux are provided, and show the capability of the quadrature-based approach to predict them in agreement with the MD simulations under dilute conditions.

The interplay between dissipation and long-range repulsive/attractive forces inhomogeneous, dilute, mono-disperse particle systems is studied. Thepseudo-Liouville operator formalism, originallyintroduced for hard-sphere interactions, is modified such that it provides very goodpredictions for systems with weak long-range forces at low densities, with the ratio ofpotential to fluctuation kinetic energy as control parameter. By numerical simulations,the theoretical results are generalized with empirical, density dependentcorrection-functions up to moderate densities.

The main result of this study on dissipative cooling is an analytical prediction for the reduced cooling rate due torepulsive forces and for the increased rate due to attractive forces. In the latter case,surprisingly, for intermediate densities, similar cooling behavior is observed as insystems without long-range interactions. In the attractive case, in general, dissipationleads to inhomogeneities earlier and faster than in the repulsive case.

The space missions Voyager and Cassini together with earthbound observations revealed a wealth of structures in Saturn’s rings. There are, for example, waves being excited at ring positions which are in orbital resonance with Saturn’s moons. Other structures can be assigned to embedded moons like empty gaps, moon induced wakes or S-shaped propeller features. Furthermore, irregular radial structures are observed in the range from 10 meters until kilometers. Here some of these structures will be discussed in the frame of hydrodynamical modeling of Saturn’s dense rings. For this purpose we will characterize the physical properties of the ring particle ensemble by mean field quantities and point to the special behavior of the transport coefficients. We show that unperturbed rings can become unstable and how diffusion acts in the rings. Additionally, the alternative streamline formalism is introduced to describe perturbed regions of dense rings with applications to the wake damping and the dispersion relation of the density waves.

A computer aided method using symbolic computations that enables the calculation of thesource terms (Boltzmann) in Grad’s method of moments is presented. The method is extremelypowerful, easy to program and allows the derivation of balance equations to very highmoments (limited only by computer resources). For sake of demonstration the method isapplied to a simple case: the one-dimensional stationary granular gas under gravity. Themethod should find applications in the field of rarefied gases, as well. Questions ofconvergence, closure are beyond the scope of this article.

We perform a numerical study of the fluctuations of the rescaled hydrodynamic transversevelocity field during the cooling state of a homogeneous granular gas. We are interestedin the role of Molecular Chaos for the amplitude of the hydrodynamic noise and itsrelaxation in time. For this purpose we compare the results of Molecular Dynamics (MD,deterministic dynamics) with those from Direct Simulation Monte Carlo (DSMC, randomprocess), where Molecular Chaos can be directly controlled. It is seen that the large timedecay of the fluctuation’s autocorrelation is always dictated by the viscosity coefficientpredicted by granular hydrodynamics, independently of the numerical scheme (MD or DSMC).On the other side, the noise amplitude in Molecular Dynamics, which is known toviolate the equilibrium Fluctuation-Dissipation relation, is not alwaysaccurately reproduced in a DSMC scheme. The agreement between the two models improves ifthe probability of recollision (controlling Molecular Chaos) is reduced by increasing thenumber of virtual particles per cells in the DSMC. This result suggests that DSMC is notnecessarily more efficient than MD, if the real number of particles is small(~103 ± 104) and if one is interested in accurately reproducefluctuations. An open question remains about the small-times behavior of theautocorrelation function in the DSMC, which in MD and in kinetic theory predictions is nota straight exponential.

Continuum mechanics (e.g., hydrodynamics, elasticity theory) is based on the assumptionthat a small set of fields provides a closed description on large space and time scales.Conditions governing the choice for these fields are discussed in the context of granularfluids and multi-component fluids. In the first case, the relevance of temperature orenergy as a hydrodynamic field is justified. For mixtures, the use of a total temperatureand single flow velocity is compared with the use of multiple species temperatures andvelocities.

The evolution of a force-free granular gas with a constant restitution coefficient isstudied by means of granular hydrodynamics. We numerically solve the hydrodynamicequations and analyze the mechanisms of cluster formation. According to our findings, thepresently accepted mode-enslaving mechanism may not be responsible for the latterphenomenon. On the contrary, we observe that the cluster formation is mainly driven byshock-waves, which spontaneously originate and develop in the system. This agrees with apreviously suggested mechanism of formation of density singularities in one-dimensionalgranular gases.

The kinetics of granular gases, including both freely cooling and steadily driven systems, is studied experimentally in quasi-two-dimensional cells. Under microgravity conditions achieved inside an aircraft flying parabolic trajectories, the frictional force is reduced. In both the freely cooling and steadily driven systems, we confirm that the velocity distribution function has the form exp(−α|v|β). The value of exponent β is close to 1.5 for the driven system in a highly excited case, which is consistent with theory derived under the assumption of the existence of the white-noise thermostat (van Noije & Ernst, Gran. Mat., vol. 1, 1998, p. 5764). In the freely cooling system, the value of β evolves from 1.5 to 1 as the cooling proceeds, and the system's energy decays algebraically (Tg = T0(1 + t/τ)−2), agreeing with Haff's law (Haff, J. Fluid Mech., vol. 134, 1983, p. 401430).