1 / 72

Liouville equation for granular gases

Liouville equation for granular gases. Hisao Hayakawa ( YITP, Kyoto Univ. ) at 2008/10/17 & Michio Otsuki ( YITP, Kyoto Univ., Dept. of Physics, Aoyama-Gakuin Univ. ). Aim of this talk. This talk is very different from others.

colemanj
Download Presentation

Liouville equation for granular gases

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Liouville equation for granular gases Hisao Hayakawa (YITP, Kyoto Univ.) at 2008/10/17 & Michio Otsuki (YITP, Kyoto Univ., Dept. of Physics, Aoyama-Gakuin Univ.)

  2. Aim of this talk • This talk is very different from others. • The purpose of this talk is what happens if local collision processes loose time-reversal symmetry.

  3. Contents • Introduction • I. What is granular materials? • II. Characteristics of sheared glassy or granular systems • Liouville equation and MCT for sheared granular gases • III. Liouville equation for sheared granular gases • IV. Generalized Langevin equation • V. MCT equation for sheared granular fluids • Spatial correlation in sheared isothermal liquids • VI. Spatial correlations in granular liquids • VII. Linearized generalized fluctuating hydrodynamics • VIII. Comparison between theory and simulation

  4. I. What is granular materials? • sand grains: grain diameter is ranged in 0.01mm-1mm. • Macroscopic particles • Energy dissipation • Repulsive systems • Granular materials • Many-body systems of dissipative particles http://science.nasa.gov/headlines/y2002/06dec_dunes.htm

  5. Granular shear flow • Coexistence of “solid” region and “fluid” region • There is creep motion in “solid” region. From H. M. Jeager, S. N. Nagel and R. P. Behringer, Rev. Mod. Phys. Vol. 68, 1259 (1996)

  6. Granular Gases (What happens if molecules are dissipative?) • Granular gases= A model of dusts • Uniform state is unstable. • (3) It is not easy to perform experiments for gases. I. Goldhirsch and G. Zanetti, Phys.Rev.Lett. 70 , 1619-1622 (1993).

  7. Simulation of a freely cooling gas The restitution 0.99118 Area fraction 0.25 # of particles640,000 Initial: equilibrium Time is scaled by the collision number The correlation grows with time. By M. Isobe(NITECH)

  8. A simple model of granular gas The shear mode for the perturbation to a uniform state is always unstable because aligned motion of particles is survived. => string-like structure

  9. Characteristics of inelastic collisions • Energy is not conserved in each collision. • Inelasticity is characterized by the restitution coefficient e<1. • There is no time reversal symmetry in each collision. • The phase volume is contracted at the instance of a collision.

  10. Characteristics of granular hydrodynamics • Theories remain in phenomenological level. • Many theories are based on eigenvalue analysis of hydrodynamic equations. • There is no sound wave in freely cooling case once inelasticity is introduced (HH and M.Otsuki,PRE2007) • There are sound waves in sheared gases.

  11. Contents • I. What is granular materials? • II. Characteristics of sheared glassy or granular systems • III. Liouville equation for sheared granular gases • IV. Generalized Langevin equation • V. MCT equation for sheared granular fluids • VI. Spatial correlations in granular liquids • VII. Linearized generalized fluctuating hydrodynamics • VIII. Comparison between theory and simulation

  12. II. Characteristics of sheared glassy or granular systems • Long time correlations: • No-decay of correlations and freezing • Correlated motion • Dynamical heterogeneity A correlated motion of a granular system (left) and a colloidal system.

  13. Similarity between jamming transition and glass transition • Granular materials exhibit “glass transition” as a jamming. • MCT can be used for sheared glass. Liu and Nagel, Nature (1998)

  14. Jamming transition • Jamming transition shows beautiful scalings (see right figs. by Otsuki and Hayakawa). • What are the properties of dense but fluidized granular liquids?

  15. Experimental relevancy of sheared systems • Recently, there are some relevant experiments of sheared granular flows.

  16. Simulation • Shear can be added with or without gravity. • For theoretical point of view, simple shear without gravity is the idealistic.

  17. Similarity between sheared granular fluids and sheared isothermal fluids • At least, the behaviors of velocity autocorrelation function, and the equal-time correlation function are common. (see M.Otsuki & HH, arXiv:0711.1421)

  18. Bagnold’s law for uniform sheared granular fluids Time scale The change of momentum Shear stress This is the relation between the temperature and the shear rate.

  19. MCT for sheared granular fluids • MCT equation can be derived for granular fluids starting from Liouville equation. • This approach ensures formal universality in granular systems and conventional glassy systems. • See HH and M. Otsuki, PTP 119, 381 (2008).

  20. Affine transformation in sheared fluids • Wave number is transferred.

  21. Contents • I. What is granular materials? • II. Characteristics of sheared glassy or granular systems • III. Liouville equation for sheared granular gases • IV. Generalized Langevin equation • V. MCT equation for sheared granular fluids • VI. Spatial correlations in granular liquids • VII. Linearized generalized fluctuating hydrodynamics • VIII. Comparison between theory and simulation

  22. III. Liouville equation for granular gases Collision operator Shear term in Liouvillian

  23. Collision operator Here, b represents the change from a collision

  24. Liouville equation

  25. Properties of Liouville operator

  26. Contents • I. What is granular materials? • II. Characteristics of sheared glassy or granular systems • III. Liouville equation for sheared granular gases • IV. Generalized Langevin equation • V. MCT equation for sheared granular fluids • VI. Spatial correlations in granular liquids • VII. Linearized generalized fluctuating hydrodynamics • VIII. Comparison between theory and simulation

  27. IV. Generalized Langevin equation

  28. Langevin equation in the steady state

  29. Some functions in generalized Langevin equation

  30. Remarks on steady state • We should note that the steady ρ(Γ)is highly nontrivial. • The steady state is determined by the balance between the external force and the inelastic collision. • Thus, the eigenvalue problem cannot be solved exactly. • In this sense, we adopt the formal argument. • I will demonstrate how to solve linearized hydrodynamics as an eigenvalue problem, later.

  31. Some formulae for hydrodynamic variables

  32. Some formulae in shear flow

  33. Generalized Langevin equation for sheared granular fluids (1) The density correlation function

  34. Generalized Langevin equation for sheared granular fluids (2)

  35. Equations for time-correlation

  36. Some formulae

  37. Contents • I. What is granular materials? • II. Characteristics of sheared glassy or granular systems • III. Liouville equation for sheared granular gases • IV. Generalized Langevin equation • V. MCT equation for sheared granular fluids • VI. Spatial correlations in granular liquids • VII. Linearized generalized fluctuating hydrodynamics • VIII. Comparison between theory and simulation

  38. V. MCT equation for sheared granular fluids Hard-core=> all terms are balanced under Bagnold’s scaling MCT approximation

  39. Preliminary simulation • We have checked the relevancy of MCT equation for sheared dense granular liquids. • MCT predicts the existence of a two-step relaxation. • Parameters: 1000 LJ particles in 3D. The system contains binary particles, and has weak shear and weak dissipation.

  40. Results of simulationfor weak shear and weak dissipation The existence of the quasi-arrested state as MCT predicts.

  41. Discussion of MCT equation • Can MCT describe the jamming transition? • The answer of the current MCT is NO. • How can we determine S(q)? • So far there is no theory to determine S(q), but it does not depend on F(q,t). No yield stress

  42. Conclusion of MCT equation for sheared granular fluids • MCT equation may be useful for very dense granular liquids. • Our model starts from hard-core liquids <=The defect of this approach • Nevertheless, our approach suggests that an unifying concept of sheared particles is useful.

  43. Contents • I. What is granular materials? • II. Characteristics of sheared glassy or granular systems • III. Liouville equation for sheared granular gases • IV. Generalized Langevin equation • V. MCT equation for sheared granular fluids • VI. Spatial correlations in granular liquids • VII. Linearized generalized fluctuating hydrodynamics • VIII. Comparison between theory and simulation

  44. VI. Spatial correlations in granular liquids • The determination of the spatial correlations in granular liquids is important in MCT. • It is known that there is a long-range velocity correlation r^{-d} (1997 Ernst, van Noije et al) for freely-cooling granular gases. • It is also known that there is long-range correlation obeying a power law in sheared isothermal liquids of elastic particles. • Lutsko and Dufty (1985,2002), Wada and Sasa (2003)

  45. Spatial correlations in sheared isothermal liquids • Let us explain how to determine the spatial correlations in terms of eigenvalue problems of linearized hydrodynamic equations. • The result is based on M. Otsuki and HH, arXiv:0809.4799.

  46. Motivation: to solve a confused situation • Lutsko (2002) obtained the structure factor of sheared molecular liquids, but his result is not consistent with the long-range correlation obtained by himself. • Many people believe that there is no contribution of the shear rate in the vicinity of glass transition. Is that true? • The spatial correlation should be determined in MCT. • Thus, we have to construct a theory to be valid for both particle scale and hydrodynamic scale.

  47. Quantities we consider

  48. Generalized fluctuating hydrodynamics (GFH) • GFH was proposed by Kirkpatrick(1985). The basic equations consists of mass and momentum conservations. We analyze an isothermal situation obtained by the balance between the heating and inelastic collisions.

  49. Properties of GFH • The effective pressure • The nonlocal viscous stress • The stress has the thermal fluctuation. The direct correlation function strain rate

  50. Characteristics of GFH • GFH includes the structure of liquids. • Generalized viscosities are represented by obtained by the eigenvalue problem of Enskog operator

More Related