1 / 38

Computer Simulation of Gravity-Driven Granular Flow

Computer Simulation of Gravity-Driven Granular Flow. Western Canadian Conference for Young Researchers in Mathematics. University of Western Ontario Department of Applied Mathematics John Drozd and Dr. Colin Denniston. Granular Matter. Granular matter definition

isabellej
Download Presentation

Computer Simulation of Gravity-Driven Granular Flow

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. Computer Simulation of Gravity-Driven Granular Flow Western Canadian Conference for Young Researchers in Mathematics University of Western Ontario Department of Applied Mathematics John Drozd and Dr. Colin Denniston

  2. Granular Matter • Granular matter definition • Small discrete particles vs. continuum • Granular matter interest • Biology, engineering, geology, material science, physics. • Mathematics and computer science. • Granular motion • Energy input and dissipation. • Granular experiments • Vibration

  3. Small Amplitude Surface Waves  Large Amplitude Surface Waves  3-Node Arching  C. Wassgren et al. 1996

  4. Simulation

  5. Other Phenomena in Granular Materials • Shear flow • Vertical shaking • Horizontal shaking • Conical hopper • Rotating drum • Cylindrical pan

  6. Cylindrical Pan Oscillations Oleh Baran et al. 2001

  7. Vibratory Drum Grinder Raw Feed Grinding Media, Rods Reactor Springs Vibrator Motor Ground Product Isolation Spring Vibratory Drum Grinder

  8. Goal • Find optimum oscillation that results in a force between the rods which achieves the ultimate stress of a particular medium that is to be crushed between the rods. • Minimize the total energy required to grind the medium. • Mixing is also important.

  9. Event-driven Simulation Without Gravity

  10. Event-driven Simulation With Gravity

  11. Typical Simulation

  12. Harry Swinney et al. 1997

  13. Harry Swinney et al. 1997

  14. Granular Materials: A window to studying the Transition from a non-Newtonian Granular Fluid To A "Glassy" system

  15. Savage and Jeffrey J. Fluid Mech. 130, 187, 1983. Velocity q Collision rules for dry granular media as modelled by inelastic hard spheres As collisions become weaker (relative velocity vn small), they become more elastic. C. Bizon et. al., PRL 80, 57, 1997.

  16. 300 (free fall region) Donev et alPRL96 "Do Binary Hard Disks Exhibit an Ideal Glass Transition?"  polydisperse monodisperse 250 (fluid region) vy 200 (glass region) dvy/dt 150 P y Polydispersity means Normal distribution of particle radii x z y Donev et al PRE 71

  17. The density in the glassy region is a constant. In the interface between the fluid and the glass does the density approach the glass density exponentially? 0 = 0.9 0 = 0.95 0 = 0.99 Interface width seems to increase as 0  1  vy y How does  depend on (1  0) ?

  18. Density vs Height in Fluid-Glass Transition

  19. Length Scale in Transition Slope = 0.42 poly Slope = 0.46 mono "interface width diverges"

  20. Y Velocity Distribution 300 (free fall region) Poiseuille flow 250 (fluid region) Plug flow snapshot 200 (glass region) Mono-disperse (crystallized) only Mono kink fracture 150 y z x

  21. 300 (free fall region) Granular Temperature fluid 250 (fluid region) 235 (At Equilibrium Temperature) 200 (glass region) equilibrium 150 y glass x z

  22. Fluctuating and Flow Velocity Experiment by N. Menon and D. J. Durian, Science, 275, 1997. 16 x 16 32 x 32 v Simulation results In Glassy Region ! J.J. Drozd and C. Denniston vf Europhysics Letters, 76 (3), 360, 2006 "questionable" averaging over nonuniform regions gives 2/3

  23.   1 in fluid glass transition For 0 = 0.9,0.95,0.96,0.97,0.98,0.99 Subtracting of Tg and vc and not averaging over regions of different vx2 Down centre Slope  = 1.0

  24. "Particle Dynamics in Sheared Granular Matter" Physical Review Letters 85, Number 7, p. 1428 (2000) Experiment: (W. Losert, L. Bocquet, T.C. Lubensky and J.P. Gollub)

  25. Velocity Fluctuations vs. Shear Rate U Slope = 0.406  0.018 Experiment Slope = 0.4 From simulation Must Subtract Tg ! Physical Review Letters 85, Number 7, (2000)

  26. q Shear Stress y x

  27. Viscosity vs Temperature Slope ~ 2  1 1.92  0.084 Transformation from a liquid to a glass takes place in a continuous manner. Relaxation times of a liquid and its Shear Viscosity increase …can do slightly better… very rapidly as Temperature is lowered. …"anomalous" viscosity. Is a fluid with "infinite" viscosity a useful description of the interior phase?

  28. Experimental data from the book: “Sands, Powders, and Grains: An Introduction to the Physics of Granular Materials” By Jacques Duran.

  29. Related to Forces: Impulse Distribution Most frequent collisions contributing to smallest impulses Simulation  Impulse defined: Magnitude of momentum after collision minus momentum beforecollision. Experiment  (Longhi, Easwar) Quasi-1d Theory  (Coppersmith, et al)

  30. Power Laws for Collision Times Collision time = time between collisions 1) spheres in 2d 2) 2d disks 3) 3d spheres Similar power laws for 2d and 3d simulations!

  31. Comparison With Experiment : experiment 1.5 vs. simulation 2.8 Discrepancy as a result of Experimental response time and sensitivity of detector.  Experiment “Spheres in 2d”: 3d Simulation with front and back reflecting walls separated one diameter apart  Pressure Transducer  Figure from experimental paper: “Large Force Fluctuations in a Flowing Granular Medium” Phys. Rev. Lett. 89, 045501 (2002) E. Longhi, N. Easwar, N. Menon

  32. Probability Distribution for Impulses vs. Collision Times (log scale)  = 2.75  = 1.50

  33. random packing at early stage  = 2.75 Is there any difference between this glass and a crystal? Answer: Look at Monodisperse grains crystallization  at later stage  = 4.3 Disorder has a universal effect on Collision Time power law.

  34. Summary of Power Laws

  35. Conclusions • A gravity-driven hard sphere simulation was used to study the glass transition from a granular hard sphere fluid to a jammed glass. • We get the same 2/3 power law for velocity fluctuations vs. flow velocity as found in experiment, when each data point is averaged over a nonuniform region. • When we look at data points averaged from a uniform region we find a power law of 1 as expected. • We found a diverging length scale at this jamming (glass) to unjamming (granular fluid) transition. • We compared our simulation to experiment on the connection between local velocity fluctuations and shear rate and found quantitative agreement. • We resolved a discrepancy with experiment on the collision time power law which we found depends on the level of disorder (glass) or order (crystal).

  36. Momentum Conservationkik+gi = 0 Normal Stresses Along Height Weight not supported by a pressure gradient.

  37. Momentum Conservation

More Related