Search

Jammed packings of granular materials display complex mechanical response. For example, the ensemble-averaged shear modulus $\left\langle G \right\rangle$ increases as a power-law in pressure $p$ for static packings of spherical particles that can rearrange during compression. We seek to design granular materials with shear moduli that can either increase {\it or} decrease with pressure without particle rearrangements even in the large-system limit. To do this, we construct {\it tessellated} granular metamaterials by joining multiple particle-filled cells together. We focus on cells that contain a small number of bidisperse disks in two dimensions. We first study the mechanical properties of individual disk-filled cells with three types of boundaries: periodic boundary conditions, fixed-length walls, and flexible walls. Hypostatic jammed packings are found for disk-filled cells with flexible walls, but not in cells with periodic boundary conditions and fixed-length walls, and they are stabilized by quartic modes of the dynamical matrix. The shear modulus of a single cell depends linearly on $p$. We find that the slope of the shear modulus with pressure, $\lambda_c < 0$ for all packings in single cells with periodic boundary conditions where the number of particles per cell $N \ge 6$. In contrast, single cells with fixed-length and flexible walls can possess $\lambda_c > 0$, as well as $\lambda_c < 0$, for $N \le 16$. We show that we can force the mechanical properties of multi-cell granular metamaterials to possess those of single cells by constraining the endpoints of the outer walls and enforcing an affine shear response. These studies demonstrate that tessellated granular metamaterials provide a novel platform for the design of soft materials with novel mechanical properties.

Most granular packings possess shear moduli ($G$) that increase with the applied external pressure, and bulk moduli ($B$) that increase or remain constant with pressure. This paper presents "tessellated" granular metamaterials for which both $G$ and the ratio $G/B$ decrease with increasing pressure. The granular metamaterials are made from flexible tessellations forming a ring of closed cells, each containing a small number of solid particles. For under-constrained tessellations, the dominant contributions to $G$ and $B$ are the particle-particle and particle-cell interactions. With specific particle configurations in the cells, we limit the number of possible particle rearrangements to achieve decreasing $G$ as we increase the pressure difference between the inside and outside of the tessellation, leading to $G/B \ll 1$ at large pressures. We further study tessellated granular metamaterials with cells containing a single particle and many particles to determine the variables that control the mechanical response of particle-filled tessellations as a function of pressure., Comment: 16 pages, 6 figures

During epithelial wound healing, cell morphology near the healed wound and the healing rate vary strongly among different developmental stages even for a single species like \textit{Drosophila}. We develop deformable particle (DP) model simulations to understand how variations in cell mechanics give rise to distinct wound closure phenotypes in the \textit{Drosophila} embryonic ectoderm and larval wing disc epithelium. We find that plastic deformation of the cell membrane can generate large changes in cell shape consistent with wound closure in the embryonic ectoderm. Our results show that the embryonic ectoderm is best described by cell membranes with an elasto-plastic response, whereas the larval wing disc is best described by cell membranes with an exclusively elastic response. By varying the mechanical response of cell membranes in DP simulations, we recapitulate the wound closure behavior of both the embryonic ectoderm and the larval wing disc.

Experimental data on the compressive strength $\sigma_{\rm max}$ versus strain rate ${\dot \varepsilon}_{\rm eng}$ for metallic glasses undergoing uniaxial compression shows significantly different behavior for different alloys. For some metallic glasses, $\sigma_{\rm max}$ decreases with increasing ${\dot \varepsilon}_{\rm eng}$, for others, $\sigma_{\rm max}$ increases with increasing ${\dot \varepsilon}_{\rm eng}$, and for others $\sigma_{\rm max}$ versus ${\dot \varepsilon}_{\rm eng}$ is nonmonotonic. Using numerical simulations of metallic glasses undergoing uniaxial compression at nonzero strain rate and temperature, we show that they obey a universal relation for the compressive strength versus temperature, which determines their mechanical response. At low ${\dot \varepsilon}_{\rm eng}$, increasing strain rate leads to increases in temperature and decreases in $\sigma^*_{\rm max}$, whereas at high ${\dot \varepsilon}_{\rm eng}$, increasing strain rate leads to decreases in temperature and increases in $\sigma^*_{\rm max}$. This non-monotonic behavior of $\sigma^*_{\rm max}$ versus temperature causes the nonmonotonic behavior of $\sigma^*_{\rm max}$ versus ${\dot \varepsilon}_{\rm eng}$. Variations in the internal dissipation change the characteristic strain rate at which the nonmonotonic behavior occurs. These results are general for a wide range of metallic glasses with different atomic interactions, damping coefficients, and chemical compositions., Comment: 5 pages, 3 figures

We employ numerical simulations to study active transistor-like switches made from two-dimensional (2D) granular crystals containing two types of grains with the same size, but different masses. We tune the mass contrast and arrangement of the grains to maximize the width of the frequency band gap in the device. The input signal is applied to a single grain on one side of the device, and the output signal is measured from another grain on the other side of the device. Changing the size of one or many grains tunes the pressure, which controls the vibrational response of the device. Switching between the on and off states is achieved using two mechanisms: 1) pressure-induced switching where the interparticle contact network is the same in the on and off states, and 2) switching through contact breaking. In general, the performance of the acoustic switch, as captured by the gain ratio and switching time between the on and off states, is better for pressure-induced switching. We show that in these acoustic switches the gain ratio between the on and off states can be larger than $10^4$ and the switching time (multiplied by the driving frequency) is comparable to that obtained recently for sonic crystals and less than that for photonic transistor-like switches. Since the self-assembly of grains with different masses into 2D granular crystals is challenging, we describe simulations of circular grains with small circular knobs placed symmetrically around the perimeter mixed with circular grains without knobs. Using umbrella sampling techniques, we show that devices with grains with $3$ knobs most efficiently form the hexagonal crystals that yield the largest band gap., Comment: 16 pages, 23 figures

We develop a numerical method to reconstruct systems of ordinary differential equations (ODEs) from time series data without {\it a priori} knowledge of the underlying ODEs using sparse basis learning and sparse function reconstruction. We show that employing sparse representations provides more accurate ODE reconstruction compared to least-squares reconstruction techniques for a given amount of time series data. We test and validate the ODE reconstruction method on known 1D, 2D, and 3D systems of ODEs. The 1D system possesses two stable fixed points; the 2D system possesses an oscillatory fixed point with closed orbits; and the 3D system displays chaotic dynamics on a strange attractor. We determine the amount of data required to achieve an error in the reconstructed functions to less than $0.1\%$. For the reconstructed 1D and 2D systems, we are able to match the trajectories from the original ODEs even at long times. For the 3D system with chaotic dynamics, as expected, the trajectories from the original and reconstructed systems do not match at long times, but the reconstructed and original models possess similar Lyapunov exponents. Now that we have validated this ODE reconstruction method on known models, it can be employed in future studies to identify new systems of ODEs using time series data from deterministic systems for which there is no currently known ODE model., Comment: 15 pages, 20 figures

We perform molecular dynamics simulations to compress binary hard spheres into jammed packings as a function of the compression rate $R$, size ratio $\alpha$, and number fraction $x_S$ of small particles to determine the connection between the glass-forming ability (GFA) and packing efficiency in bulk metallic glasses (BMGs). We define the GFA by measuring the critical compression rate $R_c$, below which jammed hard-sphere packings begin to form "random crystal" structures with defects. We find that for systems with $\alpha \gtrsim 0.8$ that do not de-mix, $R_c$ decreases strongly with $\Delta \phi_J$, as $R_c \sim \exp(-1/\Delta \phi_J^2)$, where $\Delta \phi_J$ is the difference between the average packing fraction of the amorphous packings and random crystal structures at $R_c$. Systems with $\alpha \lesssim 0.8$ partially de-mix, which promotes crystallization, but we still find a strong correlation between $R_c$ and $\Delta \phi_J$. We show that known metal-metal BMGs occur in the regions of the $\alpha$ and $x_S$ parameter space with the lowest values of $R_c$ for binary hard spheres. Our results emphasize that maximizing GFA in binary systems involves two competing effects: minimizing $\alpha$ to increase packing efficiency, while maximizing $\alpha$ to prevent de-mixing., Comment: 5 pages, 4 figures