Trapping of Vibrational Energy in Crumpled Sheets

We investigate the propagation of transverse elastic waves in crumpled media. We set up the wave equation for transverse waves on a generic curved, strained surface via a Langrangian formalism and use this to study the scaling behaviour of the dispersion curves near the ridges and on the flat facets. This analysis suggests that ridges act as barriers to wave propagation and that modes in a certain frequency regime could be trapped in the facets. A simulation study of the wave propagation qualitatively supported our analysis and showed interesting effects of the ridges on wave propagation.Comment: RevTex 12 pages, 7 figures, Submitted to PR

Anomalous strength of membranes with elastic ridges

We report on a simulational study of the compression and buckling of elastic ridges formed by joining the boundary of a flat sheet to itself. Such ridges store energy anomalously: their resting energy scales as the linear size of the sheet to the 1/3 power. We find that the energy required to buckle such a ridge is a fixed multiple of the resting energy. Thus thin sheets with elastic ridges such as crumpled sheets are qualitatively stronger than smoothly bent sheets.Comment: 4 pages, REVTEX, 3 figure

Properties of Ridges in Elastic Membranes

When a thin elastic sheet is confined to a region much smaller than its size the morphology of the resulting crumpled membrane is a network of straight ridges or folds that meet at sharp vertices. A virial theorem predicts the ratio of the total bending and stretching energies of a ridge. Small strains and curvatures persist far away from the ridge. We discuss several kinds of perturbations that distinguish a ridge in a crumpled sheet from an isolated ridge studied earlier (A. E. Lobkovsky, Phys. Rev. E. 53 3750 (1996)). Linear response as well as buckling properties are investigated. We find that quite generally, the energy of a ridge can change by no more than a finite fraction before it buckles.Comment: 13 pages, RevTeX, acknowledgement adde

Phase field model of premelting of grain boundaries

We present a phase field model of solidification which includes the effects of the crystalline orientation in the solid phase. This model describes grain boundaries as well as solid-liquid boundaries within a unified framework. With an appropriate choice of coupling of the phase field variable to the gradient of the crystalline orientation variable in the free energy, we find that high angle boundaries undergo a premelting transition. As the melting temperature is approached from below, low angle grain boundaries remain narrow. The width of the liquid layer at high angle grain boundaries diverges logarithmically. In addition, for some choices of model coupling, there may be a discontinuous jump in the width of the fluid layer as function of temperature.Comment: 6 pages, 9 figures, RevTeX

Crumpling a Thin Sheet

Crumpled sheets have a surprisingly large resistance to further compression. We have studied the crumpling of thin sheets of Mylar under different loading conditions. When placed under a fixed compressive force, the size of a crumpled material decreases logarithmically in time for periods up to three weeks. We also find hysteretic behavior when measuring the compression as a function of applied force. By using a pre-treating protocol, we control this hysteresis and find reproducible scaling behavior for the size of the crumpled material as a function of the applied force.Comment: revtex 4 pages, 6 eps figures submitted to Phys Rev. let

Scaling of the buckling transition of ridges in thin sheets

When a thin elastic sheet crumples, the elastic energy condenses into a network of folding lines and point vertices. These folds and vertices have elastic energy densities much greater than the surrounding areas, and most of the work required to crumple the sheet is consumed in breaking the folding lines or ``ridges''. To understand crumpling it is then necessary to understand the strength of ridges. In this work, we consider the buckling of a single ridge under the action of inward forcing applied at its ends. We demonstrate a simple scaling relation for the response of the ridge to the force prior to buckling. We also show that the buckling instability depends only on the ratio of strain along the ridge to curvature across it. Numerically, we find for a wide range of boundary conditions that ridges buckle when our forcing has increased their elastic energy by 20% over their resting state value. We also observe a correlation between neighbor interactions and the location of initial buckling. Analytic arguments and numerical simulations are employed to prove these results. Implications for the strength of ridges as structural elements are discussed.Comment: 42 pages, latex, doctoral dissertation, to be submitted to Phys Rev

Low temperature dynamics of kinks on Ising interfaces

The anisotropic motion of an interface driven by its intrinsic curvature or by an external field is investigated in the context of the kinetic Ising model in both two and three dimensions. We derive in two dimensions (2d) a continuum evolution equation for the density of kinks by a time-dependent and nonlocal mapping to the asymmetric exclusion process. Whereas kinks execute random walks biased by the external field and pile up vertically on the physical 2d lattice, then execute hard-core biased random walks on a transformed 1d lattice. Their density obeys a nonlinear diffusion equation which can be transformed into the standard expression for the interface velocity v = M[(gamma + gamma'')kappa + H]$, where M, gamma + gamma'', and kappa are the interface mobility, stiffness, and curvature, respectively. In 3d, we obtain the velocity of a curved interface near the orientation from an analysis of the self-similar evolution of 2d shrinking terraces. We show that this velocity is consistent with the one predicted from the 3d tensorial generalization of the law for anisotropic curvature-driven motion. In this generalization, both the interface stiffness tensor and the curvature tensor are singular at the orientation. However, their product, which determines the interface velocity, is smooth. In addition, we illustrate how this kink-based kinetic description provides a useful framework for studying more complex situations by modeling the effect of immobile dilute impurities.Comment: 11 pages, 10 figure

Limitations on the smooth confinement of an unstretchable manifold

We prove that an m-dimensional unit ball D^m in the Euclidean space {\mathbb R}^m cannot be isometrically embedded into a higher-dimensional Euclidean ball B_r^d \subset {\mathbb R}^d of radius r < 1/2 unless one of two conditions is met -- (1)The embedding manifold has dimension d >= 2m. (2) The embedding is not smooth. The proof uses differential geometry to show that if d<2m and the embedding is smooth and isometric, we can construct a line from the center of D^m to the boundary that is geodesic in both D^m and in the embedding manifold {\mathbb R}^d. Since such a line has length 1, the diameter of the embedding ball must exceed 1.Comment: 20 Pages, 3 Figure

Near-source observations and modeling of the Kuril Islands tsunamis of 15 November 2006 and 13 January 2007

International audienceTwo major earthquakes near the Central Kuril Islands (Mw=8.3 on 15 November 2006 and Mw=8.1 on 13 January 2007) generated trans-oceanic tsunamis recorded over the entire Pacific Ocean. The strongest oscillations, exceeding several meters, occurred near the source region of the Kuril Islands. Tide gauge records for both tsunamis have been thoroughly examined and numerical models of the events have been constructed. The models of the 2006 and 2007 events include two important advancements in the simulation of seismically generated tsunamis: (a) the use of the finite failure source models by Ji (2006, 2007) which provide more detailed information than conventional models on spatial displacements in the source areas and which avoid uncertainties in source extent; and (b) the use of the three-dimensional Laplace equation to reconstruct the initial tsunami sea surface elevation (avoiding the usual shallow-water approximation). The close agreement of our simulated results with the observed tsunami waveforms at the open-ocean DART stations support the validity of this approach. Observational and model findings reveal that energy fluxes of the tsunami waves from the source areas were mainly directed southeastward toward the Hawaiian Islands, with relatively little energy propagation into the Sea of Okhotsk. A marked feature of both tsunamis was their high-frequency content, with typical wave periods ranging from 2?3 to 15?20 min. Despite certain similarities, the two tsunamis were essentially different and had opposite polarity: the leading wave of the November 2006 trans-oceanic tsunami was positive, while that for the January 2007 trans-oceanic tsunami was negative. Numerical modeling of both tsunamis indicates that, due to differences in their seismic source properties, the 2006 tsunami was more wide-spread but less focused than the 2007 tsunami

Dynamic ductile to brittle transition in a one-dimensional model of viscoplasticity

We study two closely related, nonlinear models of a viscoplastic solid. These models capture essential features of plasticity over a wide range of strain rates and applied stresses. They exhibit inelastic strain relaxation and steady flow above a well defined yield stress. In this paper, we describe a first step in exploring the implications of these models for theories of fracture and related phenomena. We consider a one dimensional problem of decohesion from a substrate of a membrane that obeys the viscoplastic constitutive equations that we have constructed. We find that, quite generally, when the yield stress becomes smaller than some threshold value, the energy required for steady decohesion becomes a non-monotonic function of the decohesion speed. As a consequence, steady state decohesion at certain speeds becomes unstable. We believe that these results are relevant to understanding the ductile to brittle transition as well as fracture stability.Comment: 10 pages, REVTeX, 12 postscript figure