No Arabic abstract
The coefficient of restitution of colliding viscoelastic spheres is analytically known as a complete series expansion in terms of the impact velocity where all (infinitely many) coefficients are known. While beeing analytically exact, this result is not suitable for applications in efficient event-driven Molecular Dynamics (eMD) or Monte Carlo (MC) simulations. Based on the analytic result, here we derive expressions for the coefficient of restitution which allow for an application in efficient eMD and MC simulations of granular Systems.
The coefficient of restitution of a spherical particle in contact with a flat plate is investigated as a function of the impact velocity. As an experimental observation we notice non-trivial (non-Gaussian) fluctuations of the measured values. For a fixed impact velocity, the probability density of the coefficient of restitution, $p(epsilon)$, is formed by two exponential functions (one increasing, one decreasing) of different slope. This behavior may be explained by a certain roughness of the particle which leads to energy transfer between the linear and rotational degrees of freedom.
The current model of planet formation lacks a good understanding of the growth of dust particles inside the protoplanetary disk beyond mm sizes. In order to investigate the low-velocity collisions between this type of particles, the NanoRocks experiment was flown on the International Space Station (ISS) between September 2014 and March 2016. We present the results of this experiment. We quantify the damping of energy in systems of multiple particles in the 0.1 to 1 mm size range while they are in the bouncing regime, and study the formation of clusters through sticking collisions between particles. We developed statistical methods for the analysis of the large quantity of collision data collected by the experiment. We measured the average motion of particles, the moment of clustering, and the cluster size formed. In addition, we ran simple numerical simulations in order to validate our measurements. We computed the average coefficient of restitution (COR) of collisions and find values ranging from 0.55 for systems including a population of fine grains to 0.94 for systems of denser particles. We also measured the sticking threshold velocities and find values around 1 cm/s, consistent with the current dust collision models based on independently collected experimental data. Our findings have the following implications that can be useful for the simulation of particles in PPDs and planetary rings: (1) The average COR of collisions between same-sized free-floating particles at low speeds (< 2 cm/s) is not dependent on the collision velocity; (2) The simplified approach of using a constant COR value will accurately reproduce the average behavior of a particle system during collisional cooling; (3) At speeds below 5 mm/s, the influence of particle rotation becomes apparent on the collision behavior; (4) Current dust collision models predicting sticking thresholds are robust.
The accuracy of calculation of spectral line shapes in one-dimensional approximation is studied analytically in several limiting cases for arbitrary collision kernel and numerically in the rigid spheres model. It is shown that the deviation of the line profile is maximal in the center of the line in case of large perturber mass and intermediate values of collision frequency. For moderate masses of buffer molecules the error of one-dimensional approximation is found not to exceed 5%.
We pursue here the development of models for complex (viscoelastic) fluids in shallow free-surface gravity flows which was initiated by [Bouchut-Boyaval, M3AS (23) 2013] for 1D (translation invariant) cases. The models we propose are hyperbolic quasilinear systems that generalize Saint-Venant shallow-water equations to incompressible Maxwell fluids. The models are compatible with a formulation of the thermo-dynamics second principle. In comparison with Saint-Venant standard shallow-water model, the momentum balance includes extra-stresses associated with an elastic potential energy in addition to a hydrostatic pressure. The extra-stresses are determined by an additional tensor variable solution to a differential equation with various possible time rates. For the numerical evaluation of solutions to Cauchy problems, we also propose explicit schemes discretizing our generalized Saint-Venant systems with Finite-Volume approximations that are entropy-consistent (under a CFL constraint) in addition to satisfy exact (discrete) mass and momentum conservation laws. In comparison with most standard viscoelastic numerical models, our discrete models can be used for any retardation-time values (i.e. in the vanishing solvent-viscosity limit). We finally illustrate our hyperbolic viscoelastic flow models numerically using computer simulations in benchmark test cases. On extending to Maxwell fluids some free-shear flow testcases that are standard benchmarks for Newtonian fluids, we first show that our (numerical) models reproduce well the viscoelastic physics, phenomenologically at least, with zero retardation-time. Moreover, with a view to quantitative evaluations, numerical results in the lid-driven cavity testcase show that, in fact, our models can be compared with standard viscoelastic flow models in sheared-flow benchmarks on adequately choosing the physical parameters of our models. Analyzing our models asymptotics should therefore shed new light on the famous High-Weissenberg Number Problem (HWNP), which is a limit for all the existing viscoelastic numerical models.
We present a fully general derivation of the Laplace--Young formula and discuss the interplay between the intrinsic surface geometry and the extrinsic one ensuing from the immersion of the surface in the ordinary euclidean three-dimensional space. We prove that the (reversible) work done in a general surface deformation can be expressed in terms of the surface stress tensor and the variation of the intrinsic surface metric.