Ripples and Shear Bands in Plowed Granular Media

Monodisperse packings of dry, air-fluidized granular media typically exist between volume fractions from $\Phi$= 0.585 to 0.64. We demonstrate that the dynamics of granular drag are sensitive to volume fraction $\Phi$ and their exists a transition in the drag force and material deformation from smooth to oscillatory at a critical volume fraction $\Phi_{c}=0.605$. By dragging a submerged steel plate (3.81 cm width, 6.98 cm depth) through $300 \mu m$ glass beads prepared at volume fractions between 0.585 to 0.635 we find that below $\Phi_{c}$ the media deformation is smooth and non-localized while above $\Phi_{c}$ media fails along distinct shear bands. At high $\Phi$ the generation of these shear bands is periodic resulting in the ripples on the surface. Work funded by The Burroughs Wellcome Fund and the Army Research Lab MAST CT

Book Review: How Does the Constitution Secure Rights? Edited by Robert A. Goldwin and William Schambra.

Book review: How Does the Constitution Secure Rights? Edited by Robert A. Goldwin and William Schambra. Washington: American Enterprise Institute. 1985. Pp. ix, 125. Reviewed by: Charles Umbanhowar

Wavelength Scaling and Square/Stripe and Grain Mobility Transitions in Vertically Oscillated Granular Layers

Laboratory experiments are conducted to examine granular wave patterns near onset as a function of the container oscillation frequency f and amplitude A, layer depth H, and grain diameter D. The primary transition from a flat grain layer to standing waves occurs when the layer remains dilated after making contact with the container. With a flat layer and increasing dimensionless peak container acceleration G = 4 pi^2 f^2 A/g (g is the acceleration due to gravity), the wave transition occurs for G=2.6, but with decreasing G the waves persist to G=2.2. For 2.2<G<3.8, patterns are squares for f<f_ss and stripes for f>f_ss; H determines the square/stripe transition frequency f_ss=0.33(g/H)^0.5. The dispersion relations for layers with varying H collapse onto the curve L/H=1.0+1.1[f(H/g)^0.5]^(-1.32 +/- 0.03) (L is the wavelength) when the peak container velocity v exceeds a critical value v_gm of approximately 3 (Dg)^0.5. Local collision pressure measurements suggest that v_gm is associated with a transition in the horizontal grain mobility: for v>v_gm, there is a hydrodynamic-like horizontal sloshing motion, while for v<v_gm, the grains are essentially immobile and the stripe pattern apparently arises from a bending of the granular layer. For f at v_gm less than f_ss and v<v_gm, patterns are tenuous and disordered.Comment: 21 pages, 15 figures, submitted to Physica

The iTEBD algorithm beyond unitary evolution

The infinite time-evolving block decimation (iTEBD) algorithm [Phys. Rev. Lett. 98, 070201 (2007)] allows to simulate unitary evolution and to compute the ground state of one-dimensional quantum lattice systems in the thermodynamic limit. Here we extend the algorithm to tackle a much broader class of problems, namely the simulation of arbitrary one-dimensional evolution operators that can be expressed as a (translationally invariant) tensor network. Relatedly, we also address the problem of finding the dominant eigenvalue and eigenvector of a one-dimensional transfer matrix that can be expressed in the same way. New applications include the simulation, in the thermodynamic limit, of open (i.e. master equation) dynamics and thermal states in 1D quantum systems, as well as calculations with partition functions in 2D classical systems, on which we elaborate. The present extension of the algorithm also plays a prominent role in the infinite projected entangled-pair states (iPEPS) approach to infinite 2D quantum lattice systems.Comment: 11 pages, 16 figures, 1 appendix with algorithms for specific types of evolution. A typo in the appendix figures has been corrected. Accepted in PR

Book Review: How Does the Constitution Secure Rights? Edited by Robert A. Goldwin and William Schambra.

Book review: How Does the Constitution Secure Rights? Edited by Robert A. Goldwin and William Schambra. Washington: American Enterprise Institute. 1985. Pp. ix, 125. Reviewed by: Charles Umbanhowar

Coarsening Dynamics of Granular Heaplets in Tapped Granular Layers

A semi-continuum model is introduced to study the dynamics of the formation of granular heaplets in tapped granular layers. By taking into account the energy dissipation of collisions and screening effects due to avalanches, this model is able to reproduce qualitatively the pattern of these heaplets. Our simulations show that the granular heaplets are characterised by an effective surface tension which depends on the magnitude of the tapping intensity. Also, we observe that there is a coarsening effect in that the average size of the heaplets, V grows as the number of taps k increases. The growth law at intermediate times can be fitted by a scaling function V ~ k^z but the range of validity of the power law is limited by size effects. The growth exponent z appears to diverge as the tapping intensity is increased.Comment: 4 pages, 4 figure