Premium
Numerical Simulation of the Formation of Granular Shock Wave Over Cylindrical Obstacle
Author(s) -
Huang Yrjö Jun,
Jalali Payman,
Hyppänen Timo
Publication year - 2012
Publication title -
particle and particle systems characterization
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.877
H-Index - 56
eISSN - 1521-4117
pISSN - 0934-0866
DOI - 10.1002/ppsc.200800045
Subject(s) - mechanics , volume fraction , obstacle , shock wave , granular material , shock (circulatory) , flow (mathematics) , physics , finite volume method , front (military) , geometry , materials science , classical mechanics , thermodynamics , meteorology , mathematics , medicine , quantum mechanics , political science , law
A two dimensional (2‐D) stream of granular flow with zero initial granular temperature passing over a cylindrical obstacle is simulated by means of both molecular dynamics (MD) simulation and finite volume method (FVM). In experiments, a bow‐shaped shock wave with higher area fraction forms in front of the obstacle that was reproduced in our simulations. Due to the different circumstances to which particles are subjected, the granular flow is divided in two zones. One is undisturbed where quantities, such as space fraction (volume fraction for 3‐D and area fraction for 2‐D geometries), velocity and granular temperature are uniformly distributed and the other is called the shock wave zone. In this region, the values of the space fraction increases and the velocity of particles changes. From the MD simulation, it is found that the area fraction of the shock wave depends on surface roughness, coefficient of restitution (COR) of particles, the obstacle diameter as well as velocity of the granular stream, and a triangular region forms with almost zero velocity, and granular temperature forms in front of the cylindrical obstacle. The bigger is the size of the obstacle, the more stable this region is. In FVM simulations solid phase velocity and area fraction distributions similar to the MD simulation results are obtained for proper parameters.