-
Diffusion and velocity correlations of the phase transitions in a system of macroscopic rolling spheres
Authors:
F. Vega Reyes,
A. Rodríguez-Rivas,
J. F. González-Saavedra,
M. A. López-Castaño
Abstract:
We study an air-fluidized granular monolayer, composed of plastic spheres
which roll on a metallic grid. The air current is adjusted so that the spheres never loose contact
with the grid, so that the dynamics may be regarded as pseudo two-dimensional (or two-dimensional,
if the effects of sphere rolling are not taken into account). We find two surprising continuous
transitions, both of the…
▽ More
We study an air-fluidized granular monolayer, composed of plastic spheres
which roll on a metallic grid. The air current is adjusted so that the spheres never loose contact
with the grid, so that the dynamics may be regarded as pseudo two-dimensional (or two-dimensional,
if the effects of sphere rolling are not taken into account). We find two surprising continuous
transitions, both of them displaying two coexisting phases. Moreover, in all cases, we found the
coexisting phases display strong energy non-equipartition. In the first transition, at weak
fludization, a glassy phase coexists with a disordered fluid-like phase. In the second transition,
a hexagonal crystal coexists with the fluid phase. We analyze, for these two-phase systems, the
specific diffusive properties of each phase, as well as the velocity correlations. Surprisingly,
we find a glass phase at very low packing fraction and for a wide range of granular
temperatures. Both phases are characterized also by a strong anti-correlated velocities upon
collision. Thus, the dynamics observed for this quasi two-dimensional system unveils phase
transitions with peculiar properties, very different from the predicted behavior in well know
theories for their equilibrium counterparts.
△ Less
Submitted 11 November, 2022; v1 submitted 22 September, 2022;
originally announced September 2022.
-
Chiral flow in a binary mixture of two-dimensional active disks
Authors:
Miguel Ángel López-Castaño,
Álvaro Rodríguez-Rivas,
Francisco Vega Reyes
Abstract:
We study, experimentally, the dynamics of a binary mixture of air-fluidized disks. The disks are chiral since they incorporate a set of blades with constant tilt. Both species are identical except for their blades tilt angle, which is rotated by 180o in the second species. We analyze the phase behavior of the system. Our analysis reveals a wide range of different fluid dynamics, including chiral f…
▽ More
We study, experimentally, the dynamics of a binary mixture of air-fluidized disks. The disks are chiral since they incorporate a set of blades with constant tilt. Both species are identical except for their blades tilt angle, which is rotated by 180o in the second species. We analyze the phase behavior of the system. Our analysis reveals a wide range of different fluid dynamics, including chiral flow. This chiral flow features in its base state a large vortex. We report, for certain ranges of relative particle density of each species, inversion of the vorticity of this vortex. We discuss on the possible mechanisms behind these chiral flow transitions.
△ Less
Submitted 4 October, 2022; v1 submitted 21 June, 2022;
originally announced June 2022.
-
Diffusive regimes in a two-dimensional chiral fluid
Authors:
Francisco Vega Reyes,
Miguel A. López-Castaño,
Álvaro Rodríguez-Rivas
Abstract:
Diffusion is a fundamental aspect of transport processes in biological systems, and thus, in the development of life itself. And yet, the diffusive dynamics of active fluids with directed rotation, known as chiral fluids, has not been analyzed in detail so far. Here, we describe the diffusive regimes of a two-dimensional chiral fluid, composed in this case of a set of identical disk-shaped rotors.…
▽ More
Diffusion is a fundamental aspect of transport processes in biological systems, and thus, in the development of life itself. And yet, the diffusive dynamics of active fluids with directed rotation, known as chiral fluids, has not been analyzed in detail so far. Here, we describe the diffusive regimes of a two-dimensional chiral fluid, composed in this case of a set of identical disk-shaped rotors. We found strong experimental evidence of odd diffusion. This odd diffusion emerges in the form of a two-dimensional tensor with an antisymmetric part. In particular, we show that chiral diffusion is complex, featuring transitions between super, quasi-normal, and sub diffusion, and very slowly aging. Moreover, we show that the diffusion tensor elements, including off-diagonal elements; i.e., odd diffusion coefficient, change sign according to flow vorticity. Therefore, the chiral fluid has a self regulated diffusion, controlled by its vorticity.
△ Less
Submitted 5 October, 2022; v1 submitted 17 February, 2022;
originally announced February 2022.
-
Spin and velocity correlations in a confined two-dimensional fluid of disk-shaped active rotors
Authors:
Miguel Ángel López-Castaño,
Alejandro Márquez-Seco,
Alicia Márquez-Seco,
Álvaro Rodríguez-Rivas,
Francisco Vega Reyes
Abstract:
We study the velocity autocorrelations in an experimental configuration of confined two-dimensional active rotors (disks). We report persistent small scale oscillations in both rotational and translational velocity autocorrelations, with their characteristic frequency increasing as rotational activity increases. While these small oscillations are qualitatively similar in all experiments, we found…
▽ More
We study the velocity autocorrelations in an experimental configuration of confined two-dimensional active rotors (disks). We report persistent small scale oscillations in both rotational and translational velocity autocorrelations, with their characteristic frequency increasing as rotational activity increases. While these small oscillations are qualitatively similar in all experiments, we found that, at strong particle rotational activity, the large scale particle spin fluctuations tend to vanish, with the small oscillations around zero persisting in this case, and spins remain predominantly and strongly anti-correlated at longer times. For weaker rotational activity, however, spin fluctuations become increasingly larger and angular velocities remain de-correlated at longer times. We discuss in detail how the autocorrelation oscillations are related to the rotational activity and why this feature is generically a signal of the emergence of chirality in the dynamics of a particulate system.
△ Less
Submitted 28 February, 2023; v1 submitted 29 July, 2021;
originally announced July 2021.
-
Chirality transitions in a system of active flat spinners
Authors:
Miguel A. López-Castaño,
Alejandro Márquez Seco,
Alicia Márquez Seco,
Álvaro Rodríguez-Rivas,
Francisco Vega Reyes
Abstract:
We study in this work the 2D dynamics of an experimental system of disk-shaped rotors, fluidized by turbulent upflow. Contrary to previous knowledge, our experiments show the same particle chiral geometry can produce flows with different chiralities. In particular, we unveil a conspicuous complex chiral flow, which displays multiple persistent vortexes with either sign, located randomly in the sys…
▽ More
We study in this work the 2D dynamics of an experimental system of disk-shaped rotors, fluidized by turbulent upflow. Contrary to previous knowledge, our experiments show the same particle chiral geometry can produce flows with different chiralities. In particular, we unveil a conspicuous complex chiral flow, which displays multiple persistent vortexes with either sign, located randomly in the system. This peculiar phase mediates a continuous transition, which takes place as the kinetic energy input increases, from a flow with positive chirality (one vortex rotating in the same direction as particles spin) to a flow with negative chirality (one vortex in opposite sense to particle spin). We find that these surprising transitions are determined by the specific state of the statistical correlations between particle spin and translational velocity. We discuss how these correlations are determined in turn by the combined action of a series of mechanisms (for instance, heat dissipation at the boundaries, particle activity, average kinetic energy...), several of which are unveiled here.
△ Less
Submitted 21 September, 2022; v1 submitted 6 May, 2021;
originally announced May 2021.
-
The pseudo-two-dimensional dynamics in a system of macroscopic rolling spheres
Authors:
M. A. López-Castaño,
J. F. González-Saavedra,
A. Rodríguez-Rivas,
E. Abad,
S. B. Yuste,
F. Vega Reyes
Abstract:
We study in this work the dynamics of a collection of identical hollow spheres (ping-pong balls) that rest on a horizontal metallic grid. Fluidization is achieved by means of a turbulent air current coming from below. The upflow is adjusted so that the particles do not levitate over the grid, resulting in quasi-two-dimensional dynamics. We show that the behavior of diffusion and correlations in th…
▽ More
We study in this work the dynamics of a collection of identical hollow spheres (ping-pong balls) that rest on a horizontal metallic grid. Fluidization is achieved by means of a turbulent air current coming from below. The upflow is adjusted so that the particles do not levitate over the grid, resulting in quasi-two-dimensional dynamics. We show that the behavior of diffusion and correlations in this system is particularly rich. Noticeably as well (and related to the complex dynamical behavior), a variety of phases appear, with important peculiarities with respect to analogous setups. We observe gas, liquid, glass, and hexagonal crystal phases. Most notably, we show that the melting of the hexagonal crystal occurs in coexistence with a liquid phase. This strikingly differs from the corresponding transition in a purely two-dimensional systems of air-fluidized disks, for which no phase coexistence has been reported in the literature.
△ Less
Submitted 19 April, 2021; v1 submitted 26 June, 2020;
originally announced June 2020.
-
Memory effects in a gas of viscoelastic particles
Authors:
E. Mompó,
M. A. López-Castaño,
A. Torrente,
F. Vega Reyes,
A. Lasanta
Abstract:
We study a granular gas of viscoelastic particles (kinetic energy loss upon collision is a function of the particles' relative velocities at impact) subject to a stochastic thermostat. We show that the system displays anomalous cooling and heating rates during thermal relaxation processes, this causing the emergence of thermal memory. In particular, a significant \textit{Mpemba effect} is present;…
▽ More
We study a granular gas of viscoelastic particles (kinetic energy loss upon collision is a function of the particles' relative velocities at impact) subject to a stochastic thermostat. We show that the system displays anomalous cooling and heating rates during thermal relaxation processes, this causing the emergence of thermal memory. In particular, a significant \textit{Mpemba effect} is present; i.e., an initially hotter/cooler granular gas can cool down/heat up faster than an in comparison cooler/hotter granular gas. Moreover, a \textit{Kovacs effect} is also observed; i.e., a non-monotonic relaxation of the granular temperature --if the gas undergoes certain sudden temperature changes before fixing its value. Our results show that both memory effects have distinct features, very different and eventually opposed to those reported in theory for granular fluids under simpler collisional models. We study our system via three independent methods: approximate solution of the kinetic equation time evolution and computer simulations (both molecular dynamics simulations and Direct Simulation Monte Carlo method), finding good agreement between them.
△ Less
Submitted 15 June, 2021; v1 submitted 30 May, 2020;
originally announced June 2020.
-
Acoustic resonances in a confined set of macroscopic disks
Authors:
Juan Francisco González-Saavedra,
Álvaro Rodríguez-Rivas,
Miguel Ángel López-Castaño,
Francisco Vega Reyes
Abstract:
We study in this work a system of granular disks enclosed in a rectangular region. Granular disks are arranged, at high particle density, in a hexagonal lattice. Specifically, we are interested in the conditions for mechanical signal transmission. Two different potentials are considered. Our results are obtained by means of soft disks computer simulations. We analytically determine the eigenfreque…
▽ More
We study in this work a system of granular disks enclosed in a rectangular region. Granular disks are arranged, at high particle density, in a hexagonal lattice. Specifically, we are interested in the conditions for mechanical signal transmission. Two different potentials are considered. Our results are obtained by means of soft disks computer simulations. We analytically determine the eigenfrequencies spectra of the system from the corresponding Hessian matrix. Alternatively, the eigenfrequencies can also be obtained from the time evolution of molecular dynamics simulations. Previous works have shown that mechanical signal transmission in granular matter can be used to develop acoustic switches. Our results reveal that the working range of the granular switch changes dramatically for different interaction models. In particular, strong anisotropies in signal transmission may appear for certain interparticle forces.
△ Less
Submitted 14 November, 2019;
originally announced November 2019.
-
Statistical properties of a granular gas fluidized by turbulent air wakes
Authors:
Miguel Ángel López-Castaño,
Juan Francisco González-Saavedra,
Álvaro Rodríguez-Rivas,
Francisco Vega Reyes
Abstract:
We perform experiments with a granular system that consists of a collection of identical hollow spheres (ping-pong balls). Particles rest on a horizontal metallic grid and are confined within a circular region. Fluidization is achieved by means of a turbulent air current coming from below. Air flow is adjusted so that the balls do not elevate over the grid, as an approach to 2D dynamics. With a hi…
▽ More
We perform experiments with a granular system that consists of a collection of identical hollow spheres (ping-pong balls). Particles rest on a horizontal metallic grid and are confined within a circular region. Fluidization is achieved by means of a turbulent air current coming from below. Air flow is adjusted so that the balls do not elevate over the grid, as an approach to 2D dynamics. With a high-speed camera, we take images of the system. From these images we can infer horizontal particle positions and velocities by means of particle-tracking algorithms. With the obtained data we analyze: a) the systematic measurement error in the determination of positions and velocities from our digital images; b) the degree of homogeneity achieved in our experiments (which depends on possible deviations of the grid from the horizontal and on the homogeneity of turbulent air wakes). Interestingly, we have observed evidences of crystallization at high enough densities.
△ Less
Submitted 4 November, 2019;
originally announced November 2019.
-
Zero-gravity thermal convection in granular gases
Authors:
Álvaro Rodríguez-Rivas,
Miguel A. López-Castaño,
Francisco Vega Reyes
Abstract:
Previous experimental and theoretical evidence has shown that convective flow may appear in granular fluids if subjected to a thermal gradient and gravity (Rayleigh-Bénard-type convection). In contrast to this, we present here evidence of gravity-free thermal convection in a granular gas, with no presence of external thermal gradients either. Convection is here maintained steady by internal gradie…
▽ More
Previous experimental and theoretical evidence has shown that convective flow may appear in granular fluids if subjected to a thermal gradient and gravity (Rayleigh-Bénard-type convection). In contrast to this, we present here evidence of gravity-free thermal convection in a granular gas, with no presence of external thermal gradients either. Convection is here maintained steady by internal gradients due to dissipation and thermal sources at the same temperature. The granular gas is composed by identical disks and is enclosed in a rectangular region. Our results are obtained by means of an event-driven algorithm for inelastic hard disks.
△ Less
Submitted 18 July, 2020; v1 submitted 25 October, 2019;
originally announced October 2019.
-
Intruders in disguise: Mimicry effect in granular gases
Authors:
Antonio Lasanta,
Francisco Vega Reyes,
Vicente Garzó,
Andrés Santos
Abstract:
In general, the total kinetic energy in a multicomponent granular gas of inelastic and rough hard spheres is unequally partitioned among the different degrees of freedom. On the other hand, partial energy equipartition can be reached, in principle, under appropriate combinations of the mechanical parameters of the system. Assuming common values of the coefficients of restitution, we use kinetic-th…
▽ More
In general, the total kinetic energy in a multicomponent granular gas of inelastic and rough hard spheres is unequally partitioned among the different degrees of freedom. On the other hand, partial energy equipartition can be reached, in principle, under appropriate combinations of the mechanical parameters of the system. Assuming common values of the coefficients of restitution, we use kinetic-theory tools to determine the conditions under which the components of a granular mixture in the homogeneous cooling state have the same translational and rotational temperatures as those of a one-component granular gas (``mimicry'' effect). Given the values of the concentrations and the size ratios, the mimicry effect requires the mass ratios to take specific values, the smaller spheres having a larger particle mass density than the bigger spheres. The theoretical predictions for the case of an impurity immersed in a host granular gas are compared against both DSMC and molecular dynamics simulations with a good agreement.
△ Less
Submitted 17 June, 2019; v1 submitted 26 March, 2019;
originally announced March 2019.
-
Large Mpemba-like effect in a gas of inelastic rough hard spheres
Authors:
Aurora Torrente,
Miguel A. López-Castaño,
Antonio Lasanta,
Francisco Vega Reyes,
Antonio Prados,
Andrés Santos
Abstract:
We report the emergence of a giant Mpemba effect in the uniformly heated gas of inelastic rough hard spheres: The initially hotter sample may cool sooner than the colder one, even when the initial temperatures differ by more than one order of magnitude. In order to understand this behavior, it suffices to consider the simplest Maxwellian approximation for the velocity distribution in a kinetic app…
▽ More
We report the emergence of a giant Mpemba effect in the uniformly heated gas of inelastic rough hard spheres: The initially hotter sample may cool sooner than the colder one, even when the initial temperatures differ by more than one order of magnitude. In order to understand this behavior, it suffices to consider the simplest Maxwellian approximation for the velocity distribution in a kinetic approach. The largeness of the effect stems from the fact that the rotational and translational temperatures, which obey two coupled evolution equations, are comparable. Our theoretical predictions agree very well with molecular dynamics and direct simulation Monte Carlo data.
△ Less
Submitted 11 June, 2019; v1 submitted 25 February, 2019;
originally announced February 2019.
-
On the emergence of large and complex memory effects in nonequilibrium fluids
Authors:
Antonio Lasanta,
Francisco Vega Reyes,
Antonio Prados,
Andrés Santos
Abstract:
Control of cooling and heating processes is essential in many industrial and biological processes. In fact, the time evolution of an observable quantity may differ according to the previous history of the system. For example, a system that is being subject to cooling and then, at a given time $t_{w}$ for which the instantaneous temperature is $T(t_w)=T_{\mathrm{st}}$, is suddenly put in contact wi…
▽ More
Control of cooling and heating processes is essential in many industrial and biological processes. In fact, the time evolution of an observable quantity may differ according to the previous history of the system. For example, a system that is being subject to cooling and then, at a given time $t_{w}$ for which the instantaneous temperature is $T(t_w)=T_{\mathrm{st}}$, is suddenly put in contact with a temperature source at $T_{\mathrm{st}}$ may continue cooling down temporarily or, on the contrary, undergo a temperature rebound. According to current knowledge, there can be only one "spurious" and small peak/low. However, our results prove that, under certain conditions, more than one extremum may appear. Specifically, we have observed regions with two extrema and a critical point with three extrema. We have also detected cases where extraordinarily large extrema are observed, as large as the order of magnitude of the stationary value of the variable of interest. We show this by studying the thermal evolution of a low density set of macroscopic particles that do not preserve kinetic energy upon collision, i.e., a granular gas. We describe the mechanism that signals in this system the emergence of these complex and large memory effects, and explain why similar observations can be expected in a variety of systems.
△ Less
Submitted 11 April, 2019; v1 submitted 25 August, 2018;
originally announced August 2018.
-
Energy nonequipartition in gas mixtures of inelastic rough hard spheres: The tracer limit
Authors:
Francisco Vega Reyes,
Antonio Lasanta,
Andrés Santos,
Vicente Garzó
Abstract:
The dynamical properties of a tracer or impurity particle immersed in a host gas of inelastic and rough hard spheres in the homogeneous cooling state is studied. Specifically, the breakdown of energy equipartition as characterized by the tracer/host ratios of translational and rotational temperatures is analyzed by exploring a wide spectrum of values of the control parameters of the system (masses…
▽ More
The dynamical properties of a tracer or impurity particle immersed in a host gas of inelastic and rough hard spheres in the homogeneous cooling state is studied. Specifically, the breakdown of energy equipartition as characterized by the tracer/host ratios of translational and rotational temperatures is analyzed by exploring a wide spectrum of values of the control parameters of the system (masses, moments of inertia, sizes, and coefficients of restitution). Three complementary approaches are considered. On the theoretical side, the Boltzmann and Boltzmann--Lorentz equations (both assuming the molecular chaos ansatz) are solved by means of a multitemperature Maxwellian approximation for the velocity distribution functions. This allows us to obtain explicit analytical expressions for the temperature ratios. On the computational side, two different techniques are used. First, the kinetic equations are numerically solved by the direct simulation Monte Carlo (DSMC) method. Second, molecular dynamics simulations for dilute gases are performed. Comparison between theory and simulations shows a general good agreement. This means that (i) the impact of the molecular chaos ansatz on the temperature ratios is not significant (except at high inelasticities and/or big impurities) and (ii) the simple Maxwellian approximation yields quite reliable predictions.
△ Less
Submitted 17 September, 2019; v1 submitted 1 September, 2017;
originally announced September 2017.
-
Hydrodynamics of a Granular Gas in a Heterogeneous Environment
Authors:
Francisco Vega Reyes,
Antonio Lasanta
Abstract:
We analyze the transport properties of a low density ensemble of identical macroscopic particles immersed in an active fluid. The particles are modeled as inelastic hard spheres (granular gas). The non-homogeneous active fluid is modeled by means of a non-uniform stochastic thermostat. The theoretical results are validated with a numerical solution of the corresponding the kinetic equation (direct…
▽ More
We analyze the transport properties of a low density ensemble of identical macroscopic particles immersed in an active fluid. The particles are modeled as inelastic hard spheres (granular gas). The non-homogeneous active fluid is modeled by means of a non-uniform stochastic thermostat. The theoretical results are validated with a numerical solution of the corresponding the kinetic equation (direct simulation Monte Carlo method). We show a steady flow in the system that is accurately described by Navier-Stokes (NS) hydrodynamics, even for high inelasticity. Surprisingly, we find that the deviations from NS hydrodynamics for this flow are stronger as the inelasticity decreases. The active fluid action is modeled here with a non-uniform fluctuating volume force. This is a relevant result given that hydrodynamics of particles in complex environments, such as biological crowded environments, is still a question under intense debate.
△ Less
Submitted 14 October, 2017; v1 submitted 21 August, 2017;
originally announced August 2017.
-
Inherent thermal convection in a gas inside a box under a gravity field
Authors:
F. Vega Reyes,
A. Puglisi,
G. Pontuale,
A. Gnoli
Abstract:
We theoretically prove the existence in granular fluids of a thermal convection that is inherent, in the sense that is always present and has no thermal gradient threshold (convection occurs for all finite values of the Rayleigh number). More specifically, we study a gas of inelastic smooth hard disks enclosed in a rectangular region under a constant gravity field. The vertical walls act as energy…
▽ More
We theoretically prove the existence in granular fluids of a thermal convection that is inherent, in the sense that is always present and has no thermal gradient threshold (convection occurs for all finite values of the Rayleigh number). More specifically, we study a gas of inelastic smooth hard disks enclosed in a rectangular region under a constant gravity field. The vertical walls act as energy sinks (i.e., inelastic walls that are parallel to gravity) whereas the other two walls are perpendicular to gravity and act as energy sources. We show that this convection is due to the combined action of dissipative lateral walls and a volume force (in this case, gravitation). Hence, we call it \textit{dissipative lateral walls convection}, DLWC. Our theory, that describes also the limit case of elastic collisions, shows that inelastic particle collisions enhance the DLWC. We perform our study via numerical solutions (volume element method) of the corresponding hydrodynamic equations, in an extended Boussinesq approximation. We show our theory describes the essentials of the results for similar (but more complex) laboratory experiments.
△ Less
Submitted 28 September, 2018; v1 submitted 13 July, 2017;
originally announced July 2017.
-
Dissipative lateral walls are sufficient to trigger convection in vibrated granular gases
Authors:
G. Pontuale,
Andrea Gnoli,
Francisco Vega Reyes,
Andrea Puglisi
Abstract:
Buoyancy-driven (thermal) convection in dilute granular media, fluidized by a vibrating base, is known to appear without the need of lateral boundaries in a restricted region of parameters (inelasticity, gravity, intensity of energy injection). We have recently discovered a second buoyancy-driven convection effect which occurs at any value of the parameters, provided that the impact of particles w…
▽ More
Buoyancy-driven (thermal) convection in dilute granular media, fluidized by a vibrating base, is known to appear without the need of lateral boundaries in a restricted region of parameters (inelasticity, gravity, intensity of energy injection). We have recently discovered a second buoyancy-driven convection effect which occurs at any value of the parameters, provided that the impact of particles with the lateral walls is inelastic (Pontuale et al., Phys. Rev. Lett. 117, 098006 (2016)). It is understood that this novel convection effect is strictly correlated to the existence of perpendicular energy fluxes: a vertical one, induced by both bulk and wall inelasticity, and a horizontal one, induced only by dissipation at the walls. Here we first review those previous results, and then present new experimental and numerical data concerning the variations of box geometry, intensity of energy injection, number of particles and width of the box.
△ Less
Submitted 30 June, 2017;
originally announced June 2017.
-
When the Hotter Cools More Quickly: Mpemba Effect in Granular Fluids
Authors:
Antonio Lasanta,
Francisco Vega Reyes,
Antonio Prados,
Andrés Santos
Abstract:
Under certain conditions, two samples of fluid at different initial temperatures present a counterintuitive behavior known as the Mpemba effect: it is the hotter system that cools sooner. Here, we show that the Mpemba effect is present in granular fluids, both in uniformly heated and in freely cooling systems. In both cases, the system remains homogeneous, and no phase transition is present. Analy…
▽ More
Under certain conditions, two samples of fluid at different initial temperatures present a counterintuitive behavior known as the Mpemba effect: it is the hotter system that cools sooner. Here, we show that the Mpemba effect is present in granular fluids, both in uniformly heated and in freely cooling systems. In both cases, the system remains homogeneous, and no phase transition is present. Analytical quantitative predictions are given for how differently the system must be initially prepared to observe the Mpemba effect, the theoretical predictions being confirmed by both molecular dynamics and Monte Carlo simulations. Possible implications of our analysis for other systems are also discussed.
△ Less
Submitted 5 October, 2017; v1 submitted 15 November, 2016;
originally announced November 2016.
-
Thermal properties of an impurity immersed in a granular gas of rough hard spheres
Authors:
Francisco Vega Reyes,
Antonio Lasanta,
Andrés Santos,
Vicente Garzó
Abstract:
We study in this work the dynamics of a granular impurity immersed in a low-density granular gas of identical particles. For description of the kinetics of the granular gas and the impurity particles we use the rough hard sphere collisional model. We take into account the effects of non-conservation of energy upon particle collision. We find an (approximate) analytical solution of the pertinent ki…
▽ More
We study in this work the dynamics of a granular impurity immersed in a low-density granular gas of identical particles. For description of the kinetics of the granular gas and the impurity particles we use the rough hard sphere collisional model. We take into account the effects of non-conservation of energy upon particle collision. We find an (approximate) analytical solution of the pertinent kinetic equations for the single-particle velocity distribution functions that reproduces reasonably well the properties of translational/rotational energy non-equipartition. We assess the accuracy of the theoretical solution by comparing with computer simulations. For this, we use two independent computer data sets, from molecular dynamics (MD) and from Direct Simulation Monte Carlo method (DSMC). Both approach well, with different degrees, the kinetic theory within a reasonable range of parameter values.
△ Less
Submitted 12 July, 2017; v1 submitted 28 October, 2016;
originally announced October 2016.
-
Thermal convection in granular gases with dissipative lateral walls
Authors:
Giorgio Pontuale,
Andrea Gnoli,
Francisco Vega Reyes,
Andrea Puglisi
Abstract:
We consider a granular gas under the action of gravity, fluidized by a vibrating base. We show that a horizontal temperature gradient, here induced by limiting dissipative lateral walls (DLW), leads always to a granular thermal convection (DLW-TC) that is essentially different from ordinary buoyancy-driven convection (BD-TC). In an experiment where BD-TC is inhibited, by reducing gravity with an i…
▽ More
We consider a granular gas under the action of gravity, fluidized by a vibrating base. We show that a horizontal temperature gradient, here induced by limiting dissipative lateral walls (DLW), leads always to a granular thermal convection (DLW-TC) that is essentially different from ordinary buoyancy-driven convection (BD-TC). In an experiment where BD-TC is inhibited, by reducing gravity with an inclined plane, we always observe a DLW-TC cell next to each lateral wall. Such a cell squeezes towards the nearest wall as the gravity and/or the number of grains increase. Molecular dynamics simulations reproduce the experimental results and indicate that at large gravity or number of grains the DLW-TC is barely detectable.
△ Less
Submitted 23 August, 2016; v1 submitted 10 May, 2016;
originally announced May 2016.
-
Steady state in a gas of inelastic rough spheres heated by a uniform stochastic force
Authors:
Francisco Vega Reyes,
Andrés Santos
Abstract:
We study here the steady state attained in a granular gas of inelastic rough spheres that is subject to a spatially uniform random volume force. The stochastic force has the form of the so-called white noise and acts by adding impulse to the particle translational velocities. We work out an analytical solution of the corresponding velocity distribution function from a Sonine polynomial expansion t…
▽ More
We study here the steady state attained in a granular gas of inelastic rough spheres that is subject to a spatially uniform random volume force. The stochastic force has the form of the so-called white noise and acts by adding impulse to the particle translational velocities. We work out an analytical solution of the corresponding velocity distribution function from a Sonine polynomial expansion that displays energy non-equipartition between the translational and rotational modes, translational and rotational kurtoses, and translational-rotational velocity correlations. By comparison with a numerical solution of the Boltzmann kinetic equation (by means of the Direct Simulation Monte Carlo method) we show that our analytical solution provides a good description that is quantitatively very accurate in certain ranges of inelasticity and roughness. We also find three important features that make the forced granular gas steady state very different from the homogeneous cooling state (attained by an unforced granular gas). First, the marginal velocity distributions are always close to a Maxwellian. Second, there is a continuous transition to the purely smooth limit (where the effects of particle rotations are ignored). And third, the angular translational-rotational velocity correlations show a preference for a quasiperpendicular mutual orientation (which is called "lifted-tennis-ball" behavior).
△ Less
Submitted 4 November, 2015; v1 submitted 27 August, 2015;
originally announced August 2015.
-
Non-Newtonian hydrodynamics for a dilute granular suspension under uniform shear flow
Authors:
Moisés G. Chamorro,
F. Vega Reyes,
V. Garzó
Abstract:
We study in this work a steady shearing laminar flow with null heat flux (usually called "uniform shear flow") in a gas-solid suspension at low density. The solid particles are modeled as a gas of smooth hard spheres with inelastic collisions while the influence of the surrounding interstitial fluid on the dynamics of grains is modeled by means of a volume drag force, in the context of a rheologic…
▽ More
We study in this work a steady shearing laminar flow with null heat flux (usually called "uniform shear flow") in a gas-solid suspension at low density. The solid particles are modeled as a gas of smooth hard spheres with inelastic collisions while the influence of the surrounding interstitial fluid on the dynamics of grains is modeled by means of a volume drag force, in the context of a rheological model for suspensions. The model is solved by means of three different but complementary routes, two of them being theoretical (Grad's moment method applied to the corresponding Boltzmann equation and an exact solution of a kinetic model adapted to granular suspensions) and the other being computational (Monte Carlo simulations of the Boltzmann equation). Unlike in previous studies on granular sheared suspensions, we include in our Grad's solution nonlinear terms in the stress tensor in the collisional moment associated with the momentum transfer. This theoretical enhancement allows us for the detection and evaluation of the normal stress differences in the plane normal to the laminar flow. In addition, the exact solution of the kinetic model gives the explicit form of the velocity moments of the velocity distribution function. Comparison between our theoretical and numerical results shows in general a good agreement for the non-Newtonian rheological properties, the kurtosis (fourth velocity moment of the distribution function) and the velocity distribution of the kinetic model for quite strong inelasticity and not too large values of the (scaled) friction coefficient characterizing the viscous drag force. This shows the accuracy of our analytical results that allows us to describe in detail the flow dynamics of the granular suspension with zero heat flux throughout the paper.
△ Less
Submitted 10 November, 2015; v1 submitted 7 July, 2015;
originally announced July 2015.
-
Thermal diffusion segregation of an impurity in a driven granular fluid
Authors:
Francisco Vega Reyes,
Vicente Garzó
Abstract:
We study segregation of an impurity in a driven granular fluid under two types of \emph{steady} states. In the first state, the granular gas is driven by a stochastic volume force field with a Fourier-type profile while in the second state, the granular gas is sheared in such a way that inelastic cooling is balanced by viscous heating. We compare theoretical results derived from a solution of the…
▽ More
We study segregation of an impurity in a driven granular fluid under two types of \emph{steady} states. In the first state, the granular gas is driven by a stochastic volume force field with a Fourier-type profile while in the second state, the granular gas is sheared in such a way that inelastic cooling is balanced by viscous heating. We compare theoretical results derived from a solution of the (inelastic) Boltzmann equation at Navier-Stokes (NS) order with those obtained from the Direct Monte Carlo simulation (DSMC) method and molecular dynamics (MD) simulations. Good agreement is found between theory and simulation, which provides strong evidence of the reliability of NS granular hydrodynamics for these steady states (including the dynamics of the impurity), even at high inelasticities. In addition, preliminary results for thermal diffusion in granular fluids at moderate densitis are also presented. As for dilute gases \cite{VGK14}, excellent agreement is also found in this more general case.
△ Less
Submitted 23 July, 2014;
originally announced July 2014.
-
Properties of the Homogeneous Cooling State of a Gas of Inelastic Rough Particles
Authors:
Francisco Vega Reyes,
Andrés Santos,
Gilberto M. Kremer
Abstract:
In this work we address the question of whether a low-density system composed of identical rough particles may reach hydrodynamic states (also called \textit{normal} states), even if energy is not conserved in particle collisions. As a way to measure the ability of the system to present a hydrodynamic behavior, we focus on the so-called homogeneous cooling state of the granular gas and look at the…
▽ More
In this work we address the question of whether a low-density system composed of identical rough particles may reach hydrodynamic states (also called \textit{normal} states), even if energy is not conserved in particle collisions. As a way to measure the ability of the system to present a hydrodynamic behavior, we focus on the so-called homogeneous cooling state of the granular gas and look at the corresponding relaxation time as a function of inelasticity and roughness. We report computer simulation results of the sixth- and eighth-order cumulants of the particle velocity distribution function and study the influence of roughness on their relaxation times and asymptotic values. This extends the results of a previous work [Phys. Rev. E \textbf{89}, 020202(R) (2014], where lower-order cumulants were measured. Our results confirm that the relaxation times are not necessarily longer for stronger inelasticities. This implies that inelasticity by itself does not preclude hydrodynamics.
It is also observed that the cumulants associated with the angular velocity distribution may reach very high values in a certain region of (small) roughness and that these maxima coincide with small orientational correlation points.
△ Less
Submitted 7 September, 2018; v1 submitted 23 July, 2014;
originally announced July 2014.
-
Navier-Stokes transport coefficients for driven inelastic Maxwell models
Authors:
Moisés G. Chamorro,
Vicente Garzó,
Francisco Vega Reyes
Abstract:
We calculate in this work the Navier-Stokes transport coefficients from the Boltzmann equation for $d$-dimensional inelastic Maxwell models. By granular gas we mean here a low density system of identical spheres that lose a fraction of their kinetic energy after collisions. In the present work, the granular gas is fluidized by the presence of a thermostat that aides the system to reach a steady st…
▽ More
We calculate in this work the Navier-Stokes transport coefficients from the Boltzmann equation for $d$-dimensional inelastic Maxwell models. By granular gas we mean here a low density system of identical spheres that lose a fraction of their kinetic energy after collisions. In the present work, the granular gas is fluidized by the presence of a thermostat that aides the system to reach a steady state. The thermostat is composed by two terms: a random force and a drag force. The combined action of both forces, that act homogeneously on the granular gas, tries to mimic the interaction of the set of particles with a surrounding fluid. The Chapman-Enskog method is applied to solve the inelastic Boltzmann equation to first order in the deviations of the hydrodynamic fields from their values in the homogeneous steady state. Since the collisional cooling cannot be compensated locally for by the heat produced by the driving forces, the reference (zeroth-order) distribution function $f^{(0)}$ depends on time through its dependence on the granular temperature. To simplify the analysis and obtain explicit forms for the transport coefficients, the steady state conditions are considered. A comparison with previous results obtained for inelastic hard spheres is also carried out.
△ Less
Submitted 15 April, 2014; v1 submitted 19 February, 2014;
originally announced February 2014.
-
Navier-Stokes hydrodynamics in a granular gas fluidized by a non-uniform stochastic field
Authors:
Francisco Vega Reyes
Abstract:
We show a case of steady flow in a granular gas that, for small shear rates, is accurately described by Navier-Stokes hydrodynamics, even for high inelasticity. The (low density) granular gas is composed of identical inelastic spheres and is fluidized by a temperature field produced by a stochastic volume force plus thermal walls. We prove that the theoretical Navier-Stokes transport coefficients…
▽ More
We show a case of steady flow in a granular gas that, for small shear rates, is accurately described by Navier-Stokes hydrodynamics, even for high inelasticity. The (low density) granular gas is composed of identical inelastic spheres and is fluidized by a temperature field produced by a stochastic volume force plus thermal walls. We prove that the theoretical Navier-Stokes transport coefficients exhibit good agreement with numerical solutions of the corresponding the kinetic equation (direct simulation Monte Carlo method) in the steady state.
△ Less
Submitted 14 October, 2013;
originally announced October 2013.
-
Hydrodynamic granular segregation induced by boundary heating and shear
Authors:
Francisco Vega Reyes,
Vicente Garzó,
Nagi Khalil
Abstract:
Segregation induced by a thermal gradient of an impurity in a driven low-density granular gas is studied. The system is enclosed between two parallel walls from which we input thermal energy to the gas. We study here steady states occurring when the inelastic cooling is exactly balanced by some external energy input (stochastic force or viscous heating), resulting in a uniform heat flux. A segrega…
▽ More
Segregation induced by a thermal gradient of an impurity in a driven low-density granular gas is studied. The system is enclosed between two parallel walls from which we input thermal energy to the gas. We study here steady states occurring when the inelastic cooling is exactly balanced by some external energy input (stochastic force or viscous heating), resulting in a uniform heat flux. A segregation criterion based on Navier-Stokes granular hydrodynamics is written in terms of the tracer diffusion transport coefficients, whose dependence on the parameters of the system (masses, sizes and coefficients of restitution) is explicitly determined from a solution of the inelastic Boltzmann equation. The theoretical predictions are validated by means of Monte Carlo and molecular dynamics simulations, showing that Navier-Stokes hydrodynamics produces accurate segregation criteria even under strong shearing and/or inelasticity.
△ Less
Submitted 28 April, 2014; v1 submitted 2 October, 2013;
originally announced October 2013.
-
Role of roughness on the hydrodynamic homogeneous base state of inelastic spheres
Authors:
F. Vega Reyes,
A. Santos,
G. M. Kremer
Abstract:
A gas of inelastic rough spheres admits a spatially homogeneous base state which turns into a hydrodynamic state after a finite relaxation time. We show that this relaxation time is hardly dependent on the degree of inelasticity but increases dramatically with decreasing roughness. An accurate description of translational-rotational velocity correlations at all times is also provided. At a given i…
▽ More
A gas of inelastic rough spheres admits a spatially homogeneous base state which turns into a hydrodynamic state after a finite relaxation time. We show that this relaxation time is hardly dependent on the degree of inelasticity but increases dramatically with decreasing roughness. An accurate description of translational-rotational velocity correlations at all times is also provided. At a given inelasticity, the roughness parameter can be tuned to produce a huge distortion from the Maxwellian distribution function. The results are obtained from a Grad-like solution of the Boltzmann-Enskog equation complemented by Monte Carlo and molecular dynamics simulations.
△ Less
Submitted 19 February, 2014; v1 submitted 10 September, 2013;
originally announced September 2013.
-
Homogeneous steady states in a granular fluid driven by a stochastic bath with friction
Authors:
Moisés G. Chamorro,
Francisco Vega Reyes,
Vicente Garzó
Abstract:
The homogeneous state of a granular flow of smooth inelastic hard spheres or disks described by the Enskog-Boltzmann kinetic equation is analyzed. The granular gas is fluidized by the presence of a random force and a drag force. The combined action of both forces, that act homogeneously on the granular gas, tries to mimic the interaction of the set of particles with a surrounding fluid. The first…
▽ More
The homogeneous state of a granular flow of smooth inelastic hard spheres or disks described by the Enskog-Boltzmann kinetic equation is analyzed. The granular gas is fluidized by the presence of a random force and a drag force. The combined action of both forces, that act homogeneously on the granular gas, tries to mimic the interaction of the set of particles with a surrounding fluid. The first stochastic force thermalizes the system, providing for the necessary energy recovery to keep the system in its gas state at all times, whereas the second force allows us to mimic the action of the surrounding fluid viscosity. After a transient regime, the gas reaches a steady state characterized by a \emph{scaled} distribution function $\varphi$ that does not only depend on the dimensionless velocity $\mathbf{c}\equiv \mathbf{v}/v_0$ ($v_0$ being the thermal velocity) but also on the dimensionless driving force parameters. The dependence of $\varphi$ and its first relevant velocity moments $a_2$ and $a_3$ (which measure non-Gaussian properties of $\varphi$) on both the coefficient of restitution $\al$ and the driven parameters is widely investigated by means of the direct simulation Monte Carlo method. In addition, approximate forms for $a_2$ and $a_3$ are also derived from an expansion of $\varphi$ in Sonine polynomials. The theoretical expressions of the above Sonine coefficients agree well with simulation data, even for quite small values of $α$. Moreover, the third order expansion of the distribution function makes a significant accuracy improvement for larger velocities and inelasticities. Results also show that the non-Gaussian corrections to the distribution function $\varphi$ are smaller than those observed for undriven granular gases.
△ Less
Submitted 22 July, 2013; v1 submitted 23 May, 2013;
originally announced May 2013.
-
Transport properties for driven granular fluids in situations close to homogeneous steady states
Authors:
V. Garzó,
M. G. Chamorro,
F. Vega Reyes
Abstract:
The transport coefficients of a granular fluid driven by a stochastic bath with friction are obtained by solving the inelastic Enskog kinetic equation from the Chapman-Enskog method. The heat and momentum fluxes as well as the cooling rate are determined to first order in the deviations of the hydrodynamic field gradients from their values in the homogeneous steady state. Since the collisional coo…
▽ More
The transport coefficients of a granular fluid driven by a stochastic bath with friction are obtained by solving the inelastic Enskog kinetic equation from the Chapman-Enskog method. The heat and momentum fluxes as well as the cooling rate are determined to first order in the deviations of the hydrodynamic field gradients from their values in the homogeneous steady state. Since the collisional cooling cannot be compensated locally for the heat produced by the external driving force, the reference distribution $f^{(0)}$ (zeroth-order approximation) depends on time through its dependence on temperature. This fact gives rise to conceptual and practical difficulties not present in the undriven case. On the other hand, to simplify the analysis and given that we are interested in computing transport in the first order of deviations from the reference state, the steady-state conditions are considered to get explicit forms for the transport coefficients and the cooling rate. A comparison with recent molecular dynamics simulations for driven granular fluids shows an excellent agreement for the kinematic viscosity although some discrepancies are observed for the longitudinal viscosity and the thermal diffusivity at large densities. Finally, a linear stability analysis of the hydrodynamic equations with respect to the homogeneous steady state is performed. As expected, no instabilities are found thanks to the presence of the external bath.
△ Less
Submitted 4 January, 2017; v1 submitted 21 November, 2012;
originally announced November 2012.
-
Diffusion transport coefficients for granular binary mixtures at low density. Thermal diffusion segregation
Authors:
Vicente Garzó,
J. Aaron Murray,
Francisco Vega Reyes
Abstract:
The mass flux of a low-density granular binary mixture obtained previously by solving the Boltzmann equation by means of the Chapman-Enskog method is considered further. As in the elastic case, the associated transport coefficients $D$, $D_p$ and $D'$ are given in terms of the solutions of a set of coupled linear integral equations which are approximately solved by considering the first and second…
▽ More
The mass flux of a low-density granular binary mixture obtained previously by solving the Boltzmann equation by means of the Chapman-Enskog method is considered further. As in the elastic case, the associated transport coefficients $D$, $D_p$ and $D'$ are given in terms of the solutions of a set of coupled linear integral equations which are approximately solved by considering the first and second Sonine approximations. The diffusion coefficients are explicitly obtained as functions of the coefficients of restitution and the parameters of the mixture (masses, diameters and concentration) and their expressions hold for an arbitrary number of dimensions. In order to check the accuracy of the second Sonine correction for highly inelastic collisions, the Boltzmann equation is also numerically solved by means of the direct simulation Monte Carlo (DSMC) method to determine the mutual diffusion coefficient $D$ in some special situations (self-diffusion problem and tracer limit). The comparison with DSMC results reveals that the second Sonine approximation to $D$ improves the predictions made from the first Sonine approximation. We also study the granular segregation driven by a uni-directional thermal gradient. The segregation criterion is obtained from the so-called thermal diffusion factor $Λ$, which measures the amount of segregation parallel to the temperature gradient. The factor $Λ$ is determined here by considering the second-order Sonine forms of the diffusion coefficients and its dependence on the coefficients of restitution is widely analyzed across the parameter space of the system. The results obtained in this paper extend previous works carried out in the tracer limit (vanishing mole fraction of one of the species) by some of the authors of the present paper.
△ Less
Submitted 26 March, 2013; v1 submitted 8 October, 2012;
originally announced October 2012.
-
Steady base states for non-Newtonian granular hydrodynamics
Authors:
F. Vega Reyes,
A. Santos,
V. Garzó
Abstract:
We study in this work steady laminar flows in a low density granular gas modelled as a system of identical smooth hard spheres that collide inelastically. The system is excited by shear and temperature sources at the boundaries, which consist of two infinite parallel walls. Thus, the geometry of the system is the same that yields the planar Fourier and Couette flows in standard gases. We show that…
▽ More
We study in this work steady laminar flows in a low density granular gas modelled as a system of identical smooth hard spheres that collide inelastically. The system is excited by shear and temperature sources at the boundaries, which consist of two infinite parallel walls. Thus, the geometry of the system is the same that yields the planar Fourier and Couette flows in standard gases. We show that it is possible to describe the steady granular flows in this system, even at large inelasticities, by means of a (non-Newtonian) hydrodynamic approach. All five types of Couette-Fourier granular flows are systematically described, identifying the different types of hydrodynamic profiles. Excellent agreement is found between our classification of flows and simulation results. Also, we obtain the corresponding non-linear transport coefficients by following three independent and complementary methods: (1) an analytical solution obtained from Grad's 13-moment method applied to the inelastic Boltzmann equation, (2) a numerical solution of the inelastic Boltzmann equation obtained by means of the direct simulation Monte Carlo method and (3) event-driven molecular dynamics simulations. We find that, while Grad's theory does not describe quantitatively well all transport coefficients, the three procedures yield the same general classification of planar Couette-Fourier flows for the granular gas
△ Less
Submitted 14 December, 2012; v1 submitted 1 March, 2012;
originally announced March 2012.
-
Segregation of an intruder in a heated granular dense gas
Authors:
Vicente Garzó,
Fran Vega Reyes
Abstract:
A recent segregation criterion [V. Garzó, Phys. Rev. E \textbf{78}, 020301(R) (2008)] based on the thermal diffusion factor $Λ$ of an intruder in a heated granular gas described by the inelastic Enskog equation is revisited. The sign of $Λ$ provides a criterion for the transition between the Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE). The present theory incorporates two extra…
▽ More
A recent segregation criterion [V. Garzó, Phys. Rev. E \textbf{78}, 020301(R) (2008)] based on the thermal diffusion factor $Λ$ of an intruder in a heated granular gas described by the inelastic Enskog equation is revisited. The sign of $Λ$ provides a criterion for the transition between the Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE). The present theory incorporates two extra ingredients not accounted for by the previous theoretical attempt. First, the theory is based upon the second Sonine approximation to the transport coefficients of the mass flux of intruder. Second, the dependence of the temperature ratio (intruder temperature over that of the host granular gas) on the solid volume fraction is taken into account in the first and second Sonine approximations. In order to check the accuracy of the Sonine approximation considered, the Enskog equation is also numerically solved by means of the direct simulation Monte Carlo (DSMC) method to get the kinetic diffusion coefficient $D_0$. The comparison between theory and simulation shows that the second Sonine approximation to $D_0$ yields an improvement over the first Sonine approximation when the intruder is lighter than the gas particles in the range of large inelasticity. With respect to the form of the phase diagrams for the BNE/RBNE transition, the kinetic theory results for the factor $Λ$ indicate that while the form of these diagrams depends sensitively on the order of the Sonine approximation considered when gravity is absent, no significant differences between both Sonine solutions appear in the opposite limit (gravity dominates the thermal gradient). In the former case (no gravity), the first Sonine approximation overestimates both the RBNE region and the influence of dissipation on thermal diffusion segregation.
△ Less
Submitted 9 February, 2012;
originally announced February 2012.
-
Computer simulations of an impurity in a granular gas under planar Couette flow
Authors:
F Vega Reyes,
A Santos,
V Garzó
Abstract:
We present in this work results from numerical solutions, obtained by means of the direct simulation Monte Carlo (DSMC) method, of the Boltzmann and Boltzmann--Lorentz equations for an impurity immersed in a granular gas under planar Couette flow. The DSMC results are compared with the exact solution of a recent kinetic model for the same problem. The results confirm that, in steady states and ove…
▽ More
We present in this work results from numerical solutions, obtained by means of the direct simulation Monte Carlo (DSMC) method, of the Boltzmann and Boltzmann--Lorentz equations for an impurity immersed in a granular gas under planar Couette flow. The DSMC results are compared with the exact solution of a recent kinetic model for the same problem. The results confirm that, in steady states and over a wide range of parameter values, the state of the impurity is enslaved to that of the host gas: it follows the same flow velocity profile, its concentration (relative to that of the granular gas) is constant in the bulk region, and the impurity/gas temperature ratio is also constant. We determine also the rheological properties and nonlinear hydrodynamic transport coefficients for the impurity, finding a good semi-quantitative agreement between the DSMC results and the theoretical predictions.
△ Less
Submitted 13 July, 2011; v1 submitted 9 May, 2011;
originally announced May 2011.
-
Class of dilute granular Couette flows with uniform heat flux
Authors:
Francisco Vega Reyes,
Vicente Garzó,
Andrés Santos
Abstract:
In a recent paper [F. Vega Reyes et al., Phys. Rev. Lett. 104, 028001 (2010)] we presented a preliminary description of a special class of steady Couette flows in dilute granular gases. In all flows of this class the viscous heating is exactly balanced by inelastic cooling. This yields a uniform heat flux and a linear relationship between the local temperature and flow velocity. The class (referre…
▽ More
In a recent paper [F. Vega Reyes et al., Phys. Rev. Lett. 104, 028001 (2010)] we presented a preliminary description of a special class of steady Couette flows in dilute granular gases. In all flows of this class the viscous heating is exactly balanced by inelastic cooling. This yields a uniform heat flux and a linear relationship between the local temperature and flow velocity. The class (referred to as the LTu class) includes the Fourier flow of ordinary gases and the simple shear flow of granular gases as special cases. In the present paper we provide further support for this class of Couette flows by following four different routes, two of them being theoretical (Grad's moment method of the Boltzmann equation and exact solution of a kinetic model) and the other two being computational (molecular dynamics and Monte Carlo simulations of the Boltzmann equation). Comparison between theory and simulations shows a very good agreement for the non-Newtonian rheological properties, even for quite strong inelasticity, and a good agreement for the heat flux coefficients in the case of Grad's method, the agreement being only qualitative in the case of the kinetic model.
△ Less
Submitted 15 February, 2011; v1 submitted 19 October, 2010;
originally announced October 2010.
-
Segregation by thermal diffusion in granular shear flows
Authors:
Vicente Garzó,
Francisco Vega Reyes
Abstract:
Segregation by thermal diffusion of an intruder immersed in a sheared granular gas is analyzed from the (inelastic) Boltzmann equation. Segregation is induced by the presence of a temperature gradient orthogonal to the shear flow plane and parallel to gravity. We show that, like in analogous systems without shear, the segregation criterion yields a transition between upwards segregation and downwa…
▽ More
Segregation by thermal diffusion of an intruder immersed in a sheared granular gas is analyzed from the (inelastic) Boltzmann equation. Segregation is induced by the presence of a temperature gradient orthogonal to the shear flow plane and parallel to gravity. We show that, like in analogous systems without shear, the segregation criterion yields a transition between upwards segregation and downwards segregation. The form of the phase diagrams is illustrated in detail showing that they depend sensitively on the value of gravity relative to the thermal gradient. Two specific situations are considered: i) absence of gravity, and ii) homogeneous temperature. We find that both mechanisms (upwards and downwards segregation) are stronger and more clearly separated when compared with segregation criteria in systems without shear.
△ Less
Submitted 30 June, 2010; v1 submitted 13 April, 2010;
originally announced April 2010.
-
The effects of forcing and dissipation on phase transitions in thin granular layers
Authors:
Alexander E. Lobkovsky,
Francisco Vega Reyes,
J. S. Urbach
Abstract:
Recent experimental and computational studies of vibrated thin layers of identical spheres have shown transitions to ordered phases similar to those seen in equilibrium systems. Motivated by these results, we carry out simulations of hard inelastic spheres forced by homogenous white noise. We find a transition to an ordered state of the same symmetry as that seen in the experiments, but the clea…
▽ More
Recent experimental and computational studies of vibrated thin layers of identical spheres have shown transitions to ordered phases similar to those seen in equilibrium systems. Motivated by these results, we carry out simulations of hard inelastic spheres forced by homogenous white noise. We find a transition to an ordered state of the same symmetry as that seen in the experiments, but the clear phase separation observed in the vibrated system is absent. Simulations of purely elastic spheres also show no evidence for phase separation. We show that the energy injection in the vibrated system is dramatically different in the different phases, and suggest that this creates an effective surface tension not present in the equilibrium or randomly forced systems. We do find, however, that inelasticity suppresses the onset of the ordered phase with random forcing, as is observed in the vibrating system, and that the amount of the suppression is proportional to the degree of inelasticity. The suppression depends on the details of the energy injection mechanism, but is completely eliminated when inelastic collisions are replaced by uniform system-wide energy dissipation.
△ Less
Submitted 6 January, 2010;
originally announced January 2010.
-
An exact solution of the inelastic Boltzmann equation for the Couette flow with uniform heat flux
Authors:
Andrés Santos,
Vicente Garzó,
Francisco Vega Reyes
Abstract:
In the steady Couette flow of a granular gas the sign of the heat flux gradient is governed by the competition between viscous heating and inelastic cooling. We show from the Boltzmann equation for inelastic Maxwell particles that a special class of states exists where the viscous heating and the inelastic cooling exactly compensate each other at every point, resulting in a uniform heat flux. In t…
▽ More
In the steady Couette flow of a granular gas the sign of the heat flux gradient is governed by the competition between viscous heating and inelastic cooling. We show from the Boltzmann equation for inelastic Maxwell particles that a special class of states exists where the viscous heating and the inelastic cooling exactly compensate each other at every point, resulting in a uniform heat flux. In this state the (reduced) shear rate is enslaved to the coefficient of restitution $α$, so that the only free parameter is the (reduced) thermal gradient $ε$. It turns out that the reduced moments of order $k$ are polynomials of degree $k-2$ in $ε$, with coefficients that are nonlinear functions of $α$. In particular, the rheological properties ($k=2$) are independent of $ε$ and coincide exactly with those of the simple shear flow. The heat flux ($k=3$) is linear in the thermal gradient (generalized Fourier's law), but with an effective thermal conductivity differing from the Navier--Stokes one. In addition, a heat flux component parallel to the flow velocity and normal to the thermal gradient exists. The theoretical predictions are validated by comparison with direct Monte Carlo simulations for the same model.
△ Less
Submitted 31 May, 2010; v1 submitted 9 December, 2009;
originally announced December 2009.
-
Non-Newtonian granular hydrodynamics. What do the inelastic simple shear flow and the elastic Fourier flow have in common?
Authors:
Francisco Vega Reyes,
Andrés Santos,
Vicente Garzó
Abstract:
We describe a special class of steady Couette flows in dilute granular gases admitting a non-Newtonian hydrodynamic description for strong dissipation. The class occurs when viscous heating exactly balances inelastic cooling, resulting in a uniform heat flux. It includes the Fourier flow of ordinary gases and the simple or uniform shear flow (USF) of granular gases as special cases. The rheologi…
▽ More
We describe a special class of steady Couette flows in dilute granular gases admitting a non-Newtonian hydrodynamic description for strong dissipation. The class occurs when viscous heating exactly balances inelastic cooling, resulting in a uniform heat flux. It includes the Fourier flow of ordinary gases and the simple or uniform shear flow (USF) of granular gases as special cases. The rheological functions have the same values as in the USF and generalized thermal conductivity coefficients can be identified. These points are confirmed by molecular dynamics simulations, Monte Carlo simulations of the Boltzmann equation, and analytical results from Grad's 13-moment method.
△ Less
Submitted 13 January, 2010; v1 submitted 3 July, 2009;
originally announced July 2009.
-
Mass transport of impurities in a moderately dense granular gas
Authors:
Vicente Garzo,
Francisco Vega Reyes
Abstract:
Transport coefficients associated with the mass flux of impurities immersed in a moderately dense granular gas of hard disks or spheres described by the inelastic Enskog equation are obtained by means of the Chapman-Enskog expansion. The transport coefficients are determined as the solutions of a set of coupled linear integral equations recently derived for polydisperse granular mixtures [V. Gar…
▽ More
Transport coefficients associated with the mass flux of impurities immersed in a moderately dense granular gas of hard disks or spheres described by the inelastic Enskog equation are obtained by means of the Chapman-Enskog expansion. The transport coefficients are determined as the solutions of a set of coupled linear integral equations recently derived for polydisperse granular mixtures [V. Garzó, J. W. Dufty and C. M. Hrenya, Phys. Rev. E {\bf 76}, 031304 (2007)]. With the objective of obtaining theoretical expressions for the transport coefficients that are sufficiently accurate for highly inelastic collisions, we solve the above integral equations by using the second Sonine approximation. As a complementary route, we numerically solve by means of the direct simulation Monte Carlo method (DSMC) the inelastic Enskog equation to get the kinetic diffusion coefficient $D_0$ for two and three dimensions. We have observed in all our simulations that the disagreement, for arbitrarily large inelasticity, in the values of both solutions (DSMC and second Sonine approximation) is less than 4%. Moreover, we show that the second Sonine approximation to $D_0$ yields a dramatic improvement (up to 50%) over the first Sonine approximation for impurity particles lighter than the surrounding gas and in the range of large inelasticity. The results reported in this paper are of direct application in important problems in granular flows, such as segregation driven by gravity and a thermal gradient. We analyze here the segregation criteria that result from our theoretical expressions of the transport coefficients.
△ Less
Submitted 3 March, 2009; v1 submitted 17 December, 2008;
originally announced December 2008.
-
Steady base states for Navier-Stokes granular hydrodynamics with boundary heating and shear
Authors:
F. Vega Reyes,
J. S. Urbach
Abstract:
We study the Navier-Stokes steady states for a low density monodisperse hard sphere granular gas (i.e, a hard sphere ideal monatomic gas with inelastic interparticle collisions). We present a classification of the uniform steady states that can arise from shear and temperature (or energy input) applied at the boundaries (parallel walls). We consider both symmetric and asymmetric boundary conditi…
▽ More
We study the Navier-Stokes steady states for a low density monodisperse hard sphere granular gas (i.e, a hard sphere ideal monatomic gas with inelastic interparticle collisions). We present a classification of the uniform steady states that can arise from shear and temperature (or energy input) applied at the boundaries (parallel walls). We consider both symmetric and asymmetric boundary conditions and find steady states not previously reported, including sheared states with linear temperature profiles. We provide explicit expressions for the hydrodynamic profiles for all these steady states. Our results are validated by the numerical solution of the Boltzmann kinetic equation for the granular gas obtained by the direct simulation Monte Carlo method, and by molecular dynamics simulations. We discuss the physical origin of the new steady states and derive conditions for the validity of Navier-Stokes hydrodynamics.
△ Less
Submitted 25 April, 2009; v1 submitted 31 July, 2008;
originally announced July 2008.
-
Modified Sonine approximation for granular binary mixtures
Authors:
Vicente Garzó,
Francisco Vega Reyes,
José María Montanero
Abstract:
We evaluate in this work the hydrodynamic transport coefficients of a granular binary mixture in $d$ dimensions. In order to eliminate the observed disagreement (for strong dissipation) between computer simulations and previously calculated theoretical transport coefficients for a monocomponent gas, we obtain explicit expressions of the seven Navier-Stokes transport coefficients with the use of…
▽ More
We evaluate in this work the hydrodynamic transport coefficients of a granular binary mixture in $d$ dimensions. In order to eliminate the observed disagreement (for strong dissipation) between computer simulations and previously calculated theoretical transport coefficients for a monocomponent gas, we obtain explicit expressions of the seven Navier-Stokes transport coefficients with the use of a new Sonine approach in the Chapman-Enskog theory. Our new approach consists in replacing, where appropriate in the Chapman-Enskog procedure, the Maxwell-Boltzmann distribution weight function (used in the standard first Sonine approximation) by the homogeneous cooling state distribution for each species. The rationale for doing this lies in the fact that, as it is well known, the non-Maxwellian contributions to the distribution function of the granular mixture become more important in the range of strong dissipation we are interested in. The form of the transport coefficients is quite common in both standard and modified Sonine approximations, the distinction appearing in the explicit form of the different collision frequencies associated with the transport coefficients. Additionally, we numerically solve by means of the direct simulation Monte Carlo method the inelastic Boltzmann equation to get the diffusion and the shear viscosity coefficients for two and three dimensions. As in the case of a monocomponent gas, the modified Sonine approximation improves the estimates of the standard one, showing again the reliability of this method at strong values of dissipation.
△ Less
Submitted 3 March, 2009; v1 submitted 11 June, 2008;
originally announced June 2008.
-
Rheological properties of a granular impurity in the Couette flow
Authors:
Francisco Vega Reyes,
Vicente Garzo,
Andres Santos
Abstract:
We discuss in this work the validity of the theoretical solution of the nonlinear Couette flow for a granular impurity obtained in a recent work [preprint arXiv:0802.0526], in the range of large inelasticity and shear rate. We show there is a good agreement between the theoretical solution and Monte Carlo simulation data, even under these extreme conditions. We also discuss an extended theoretic…
▽ More
We discuss in this work the validity of the theoretical solution of the nonlinear Couette flow for a granular impurity obtained in a recent work [preprint arXiv:0802.0526], in the range of large inelasticity and shear rate. We show there is a good agreement between the theoretical solution and Monte Carlo simulation data, even under these extreme conditions. We also discuss an extended theoretical solution that would work for large inelasticities in ranges of shear rate $a$ not covered by our previous work (i.e., below the threshold value $a_{th}$ for which uniform shear flow may be obtained) and compare also with simulation data. Preliminary results in the simulations give useful insight in order to obtain an exact and general solution of the nonlinear Couette flow (both for $a\ge a_{th}$ and $a<a_{th}$).
△ Less
Submitted 16 April, 2008;
originally announced April 2008.
-
Effect of inelasticity on the phase transitions of a thin vibrated granular layer
Authors:
Francisco Vega Reyes,
Jeffrey S. Urbach
Abstract:
We describe an experimental and computational investigation of the ordered and disordered phases of a vibrating thin, dense granular layer composed of identical metal spheres. We compare the results from spheres with different amounts of inelasticity and show that inelasticity has a strong effect on the phase diagram. We also report the melting of an ordered phase to a homogeneous disordered liq…
▽ More
We describe an experimental and computational investigation of the ordered and disordered phases of a vibrating thin, dense granular layer composed of identical metal spheres. We compare the results from spheres with different amounts of inelasticity and show that inelasticity has a strong effect on the phase diagram. We also report the melting of an ordered phase to a homogeneous disordered liquid phase at high vibration amplitude or at large inelasticities. Our results show that dissipation has a strong effect on ordering and that in this system ordered phases are absent entirely in highly inelastic materials.
△ Less
Submitted 20 November, 2008; v1 submitted 7 March, 2008;
originally announced March 2008.
-
Impurity in a granular gas under nonlinear Couette flow
Authors:
F. Vega Reyes,
V. Garzo,
A. Santos
Abstract:
We study in this work the transport properties of an impurity immersed in a granular gas under stationary nonlinear Couette flow. The starting point is a kinetic model for low-density granular mixtures recently proposed by the authors [Vega Reyes F et al. 2007 Phys. Rev. E 75 061306]. Two routes have been considered. First, a hydrodynamic or normal solution is found by exploiting a formal mappin…
▽ More
We study in this work the transport properties of an impurity immersed in a granular gas under stationary nonlinear Couette flow. The starting point is a kinetic model for low-density granular mixtures recently proposed by the authors [Vega Reyes F et al. 2007 Phys. Rev. E 75 061306]. Two routes have been considered. First, a hydrodynamic or normal solution is found by exploiting a formal mapping between the kinetic equations for the gas particles and for the impurity. We show that the transport properties of the impurity are characterized by the ratio between the temperatures of the impurity and gas particles and by five generalized transport coefficients: three related to the momentum flux (a nonlinear shear viscosity and two normal stress differences) and two related to the heat flux (a nonlinear thermal conductivity and a cross coefficient measuring a component of the heat flux orthogonal to the thermal gradient). Second, by means of a Monte Carlo simulation method we numerically solve the kinetic equations and show that our hydrodynamic solution is valid in the bulk of the fluid when realistic boundary conditions are used. Furthermore, the hydrodynamic solution applies to arbitrarily (inside the continuum regime) large values of the shear rate, of the inelasticity, and of the rest of parameters of the system. Preliminary simulation results of the true Boltzmann description show the reliability of the nonlinear hydrodynamic solution of the kinetic model. This shows again the validity of a hydrodynamic description for granular flows, even under extreme conditions, beyond the Navier-Stokes domain.
△ Less
Submitted 5 September, 2008; v1 submitted 4 February, 2008;
originally announced February 2008.
-
Granular mixtures modeled as elastic hard spheres subject to a drag force
Authors:
Francisco Vega Reyes,
Vicente Garzo,
Andres Santos
Abstract:
Granular gaseous mixtures under rapid flow conditions are usually modeled by a multicomponent system of smooth inelastic hard spheres with constant coefficients of normal restitution. In the low density regime an adequate framework is provided by the set of coupled inelastic Boltzmann equations. Due to the intricacy of the inelastic Boltzmann collision operator, in this paper we propose a simple…
▽ More
Granular gaseous mixtures under rapid flow conditions are usually modeled by a multicomponent system of smooth inelastic hard spheres with constant coefficients of normal restitution. In the low density regime an adequate framework is provided by the set of coupled inelastic Boltzmann equations. Due to the intricacy of the inelastic Boltzmann collision operator, in this paper we propose a simpler model of elastic hard spheres subject to the action of an effective drag force, which mimics the effect of dissipation present in the original granular gas. The Navier--Stokes transport coefficients for a binary mixture are obtained from the model by application of the Chapman--Enskog method. The three coefficients associated with the mass flux are the same as those obtained from the inelastic Boltzmann equation, while the remaining four transport coefficients show a general good agreement, especially in the case of the thermal conductivity. Finally, the approximate decomposition of the inelastic Boltzmann collision operator is exploited to construct a model kinetic equation for granular mixtures as a direct extension of a known kinetic model for elastic collisions.
△ Less
Submitted 24 April, 2007; v1 submitted 23 January, 2007;
originally announced January 2007.
-
The dynamics of thin vibrated granular layers
Authors:
P. Melby,
F. Vega Reyes,
A. Prevost,
R. Robertson,
P. Kumar,
D. A. Egolf,
J. S. Urbach
Abstract:
We describe a series of experiments and computer simulations on vibrated granular media in a geometry chosen to eliminate gravitationally induced settling. The system consists of a collection of identical spherical particles on a horizontal plate vibrating vertically, with or without a confining lid. Previously reported results are reviewed, including the observation of homogeneous, disordered l…
▽ More
We describe a series of experiments and computer simulations on vibrated granular media in a geometry chosen to eliminate gravitationally induced settling. The system consists of a collection of identical spherical particles on a horizontal plate vibrating vertically, with or without a confining lid. Previously reported results are reviewed, including the observation of homogeneous, disordered liquid-like states, an instability to a `collapse' of motionless spheres on a perfect hexagonal lattice, and a fluctuating, hexagonally ordered state. In the presence of a confining lid we see a variety of solid phases at high densities and relatively high vibration amplitudes, several of which are reported for the first time in this article. The phase behavior of the system is closely related to that observed in confined hard-sphere colloidal suspensions in equilibrium, but with modifications due to the effects of the forcing and dissipation. We also review measurements of velocity distributions, which range from Maxwellian to strongly non-Maxwellian depending on the experimental parameter values. We describe measurements of spatial velocity correlations that show a clear dependence on the mechanism of energy injection. We also report new measurements of the velocity autocorrelation function in the granular layer and show that increased inelasticity leads to enhanced particle self-diffusion.
△ Less
Submitted 22 December, 2004;
originally announced December 2004.