Scale-free bursting activity in shrinkage induced cracking.

Based on computer simulations of a realistic discrete element model we demonstrate that shrinkage induced cracking of thin layers of heterogeneous materials, generating spectacular crack patterns, proceeds in bursts. These crackling pulses are characterized by scale free distributions of size and duration, however, with non-universal exponents depending on the system size and shrinking rate. On the contrary, local avalanches composed of micro-cracking events with temporal and spatial correlation are found to obey a universal power law statistics. Most notably, we demonstrate that the observed non-universality of the integrated signal is the consequence of the temporal superposition of the underlying local avalanches, which pop up in an uncorrelated way in homogeneous systems. Our results provide an explanation of recent acoustic emission measurements on drying induced shrinkage cracking and may have relevance for the acoustic monitoring of the electro-mechanical degradation of battery electrodes.

Scaling laws of failure dynamics on complex networks.

The topology of the network of load transmitting connections plays an essential role in the cascading failure dynamics of complex systems driven by the redistribution of load after local breakdown events. In particular, as the network structure is gradually tuned from regular to completely random a transition occurs from the localized to mean field behavior of failure spreading. Based on finite size scaling in the fiber bundle model of failure phenomena, here we demonstrate that outside the localized regime, the load bearing capacity and damage tolerance on the macro-scale, and the statistics of clusters of failed nodes on the micro-scale obey scaling laws with exponents which depend on the topology of the load transmission network and on the degree of disorder of the strength of nodes. Most notably, we show that the spatial structure of damage governs the emergence of the localized to mean field transition: as the network gets gradually randomized failed clusters formed on locally regular patches merge through long range links generating a percolation like transition which reduces the load concentration on the network. The results may help to design network structures with an improved robustness against cascading failure.

Effect of the loading condition on the statistics of crackling noise accompanying the failure of porous rocks.

We test the hypothesis that loading conditions affect the statistical features of crackling noise accompanying the failure of porous rocks by performing discrete element simulations of the tensile failure of model rocks and comparing the results to those of compressive simulations of the same samples. Cylindrical samples are constructed by sedimenting randomly sized spherical particles connected by beam elements representing the cementation of granules. Under a slowly increasing external tensile load, the cohesive contacts between particles break in bursts whose size fluctuates over a broad range. Close to failure breaking avalanches are found to localize on a highly stressed region where the catastrophic avalanche is triggered and the specimen breaks apart along a spanning crack. The fracture plane has a random position and orientation falling most likely close to the centre of the specimen perpendicular to the load direction. In spite of the strongly different strengths, degrees of 'brittleness' and spatial structure of damage of tensile and compressive failure of model rocks, our calculations revealed that the size, energy and duration of avalanches, and the waiting time between consecutive events all obey scale-free statistics with power law exponents which agree within their error bars in the two loading cases.

Temporal evolution of failure avalanches of the fiber bundle model on complex networks.

We investigate how the interplay of the topology of the network of load transmitting connections and the amount of disorder of the strength of the connected elements determines the temporal evolution of failure cascades driven by the redistribution of load following local failure events. We use the fiber bundle model of materials' breakdown assigning fibers to the sites of a square lattice, which is then randomly rewired using the Watts-Strogatz technique. Gradually increasing the rewiring probability, we demonstrate that the bundle undergoes a transition from the localized to the mean field universality class of breakdown phenomena. Computer simulations revealed that both the size and the duration of failure cascades are power law distributed on all network topologies with a crossover between two regimes of different exponents. The temporal evolution of cascades is described by a parabolic profile with a right handed asymmetry, which implies that cascades start slowly, then accelerate, and eventually stop suddenly. The degree of asymmetry proved to be characteristic of the network topology gradually decreasing with increasing rewiring probability. Reducing the variance of fibers' strength, the exponents of the size and the duration distribution of cascades increase in the localized regime of the failure process, while the localized to mean field transition becomes more abrupt. The consistency of the results is supported by a scaling analysis relating the characteristic exponents of the statistics and dynamics of cascades.

Curvature flows, scaling laws and the geometry of attrition under impacts.

Impact induced attrition processes are, beyond being essential models of industrial ore processing, broadly regarded as the key to decipher the provenance of sedimentary particles. Here we establish the first link between microscopic, particle-based models and the mean field theory for these processes. Based on realistic computer simulations of particle-wall collision sequences we first identify the well-known damage and fragmentation energy phases, then we show that the former is split into the abrasion phase with infinite sample lifetime (analogous to Sternberg's Law) at finite asymptotic mass and the cleavage phase with finite sample lifetime, decreasing as a power law of the impact velocity (analogous to Basquin's Law). This splitting establishes the link between mean field models (curvature-driven partial differential equations) and particle-based models: only in the abrasion phase does shape evolution emerging in the latter reproduce with startling accuracy the spatio-temporal patterns (two geometric phases) predicted by the former.

Evolution of anisotropic crack patterns in shrinking material layers.

Anisotropic crack patterns emerging in desiccating layers of pastes on a substrate can be exploited for controlled cracking with potential applications in microelectronic manufacturing. We investigate such possibilities of crack patterning in the framework of a discrete element model focusing on the temporal and spatial evolution of anisotropic crack patterns as a thin material layer gradually shrinks. In the model a homogeneous material is considered with an inherent structural disorder where anisotropy is captured by the directional dependence of the local cohesive strength. We demonstrate that there exists a threshold anisotropy below which crack initiation and propagation is determined by the disordered micro-structure, giving rise to cellular crack patterns. When the strength of anisotropy is sufficiently high, cracking is found to evolve through three distinct phases of aligned cracking which slices the sample, secondary cracking in the perpendicular direction, and finally binary fragmentation following the formation of a connected crack network. The anisotropic crack pattern results in fragments with a shape anisotropy which gradually gets reduced as binary fragmentation proceeds. The statistics of fragment masses exhibits a high degree of robustness described by a log-normal functional form at all anisotropies.

Impact-induced transition from damage to perforation.

We investigate the impact-induced damage and fracture of a bar-shaped specimen of heterogeneous materials focusing on how the system approaches perforation as the impact energy is gradually increased. A simple model is constructed which represents the bar as two rigid blocks coupled by a breakable interface with disordered local strength. The bar is clamped at the two ends, and the fracture process is initiated by an impactor hitting the bar in the middle. Our calculations revealed that depending on the imparted energy, the system has two phases: at low impact energies the bar suffers damage but keeps its integrity, while at sufficiently high energies, complete perforation occurs. We demonstrate that the transition from damage to perforation occurs analogous to continuous phase transitions. Approaching the critical point from below, the intact fraction of the interface goes to zero, while the deformation rate of the bar diverges according to power laws as function of the distance from the critical energy. As the degree of disorder increases, farther from the transition point the critical exponents agree with their zero disorder counterparts; however, close to the critical point a crossover occurs to a higher exponent.

Plato's cube and the natural geometry of fragmentation.

Plato envisioned Earth's building blocks as cubes, a shape rarely found in nature. The solar system is littered, however, with distorted polyhedra-shards of rock and ice produced by ubiquitous fragmentation. We apply the theory of convex mosaics to show that the average geometry of natural two-dimensional (2D) fragments, from mud cracks to Earth's tectonic plates, has two attractors: "Platonic" quadrangles and "Voronoi" hexagons. In three dimensions (3D), the Platonic attractor is dominant: Remarkably, the average shape of natural rock fragments is cuboid. When viewed through the lens of convex mosaics, natural fragments are indeed geometric shadows of Plato's forms. Simulations show that generic binary breakup drives all mosaics toward the Platonic attractor, explaining the ubiquity of cuboid averages. Deviations from binary fracture produce more exotic patterns that are genetically linked to the formative stress field. We compute the universal pattern generator establishing this link, for 2D and 3D fragmentation.

Record statistics of bursts signals the onset of acceleration towards failure.

Forecasting the imminent catastrophic failure has a high importance for a large variety of systems from the collapse of engineering constructions, through the emergence of landslides and earthquakes, to volcanic eruptions. Failure forecast methods predict the lifetime of the system based on the time-to-failure power law of observables describing the final acceleration towards failure. We show that the statistics of records of the event series of breaking bursts, accompanying the failure process, provides a powerful tool to detect the onset of acceleration, as an early warning of the impending catastrophe. We focus on the fracture of heterogeneous materials using a fiber bundle model, which exhibits transitions between perfectly brittle, quasi-brittle, and ductile behaviors as the amount of disorder is increased. Analyzing the lifetime of record size bursts, we demonstrate that the acceleration starts at a characteristic record rank, below which record breaking slows down due to the dominance of disorder in fracturing, while above it stress redistribution gives rise to an enhanced triggering of bursts and acceleration of the dynamics. The emergence of this signal depends on the degree of disorder making both highly brittle fracture of low disorder materials, and ductile fracture of strongly disordered ones, unpredictable.

System-size-dependent avalanche statistics in the limit of high disorder.

We investigate the effect of the amount of disorder on the statistics of breaking bursts during the quasistatic fracture of heterogeneous materials. We consider a fiber bundle model where the strength of single fibers is sampled from a power-law distribution over a finite range, so that the amount of materials' disorder can be controlled by varying the power-law exponent and the upper cutoff of fibers' strength. Analytical calculations and computer simulations, performed in the limit of equal load sharing, revealed that depending on the disorder parameters the mechanical response of the bundle is either perfectly brittle where the first fiber breaking triggers a catastrophic avalanche, or it is quasibrittle where macroscopic failure is preceded by a sequence of bursts. In the quasibrittle phase, the statistics of avalanche sizes is found to show a high degree of complexity. In particular, we demonstrate that the functional form of the size distribution of bursts depends on the system size: for large upper cutoffs of fibers' strength, in small systems the sequence of bursts has a high degree of stationarity characterized by a power-law size distribution with a universal exponent. However, for sufficiently large bundles the breaking process accelerates towards the critical point of failure, which gives rise to a crossover between two power laws. The transition between the two regimes occurs at a characteristic system size which depends on the disorder parameters.

Effect of disorder on the spatial structure of damage in slowly compressed porous rocks.

Faults and damage zone properties control a range of important phenomena, from the hydraulic properties of underground reservoirs to the physics of earthquakes on a larger scale. Here, we investigate the effect of disorder of porous rocks on the spatial structure of damage emerging under compression. Model rock samples are numerically generated by sedimenting particles where the amount of disorder is controlled by the particle size distribution. To obtain damage bands with a sufficiently large length along axis, we performed simulations of 'Brazilian'-type compression tests of cylindrical samples. As failure is approached, damage localization leads to the formation of two conjugate shear bands. The orientation angle of bands to the loading direction increases with disorder, implying a decrease in the internal coefficient of friction. The width of the damage band scales as a power law of the degree of disorder. Inside the damage band, the sample is crushed into a large number of pieces with a power law mass distribution. The shape of fragments undergoes a crossover at a disorder-dependent size from the isotropy of small pieces to the anisotropic flattened form of the large ones. The results provide important constraints in understanding the role of disorder in geological fractures.This article is part of the theme issue 'Statistical physics of fracture and earthquakes'.

Time-dependent fracture under unloading in a fiber bundle model.

We investigate the fracture of heterogeneous materials occurring under unloading from an initial load. Based on a fiber bundle model of time-dependent fracture, we show that depending on the unloading rate the system has two phases: for rapid unloading the system suffers only partial failure and it has an infinite lifetime, while at slow unloading macroscopic failure occurs in a finite time. The transition between the two phases proved to be analogous to continuous phase transitions. Computer simulations revealed that during unloading the fracture proceeds in bursts of local breakings triggered by slowly accumulating damage. In both phases the time evolution starts with a relaxation of the bursting activity characterized by a universal power-law decay of the burst rate. In the phase of finite lifetime the initial slowdown is followed by an acceleration towards macroscopic failure where the increasing rate of bursts obeys the (inverse) Omori law of earthquakes. We pointed out a strong correlation between the time where the event rate reaches a minimum value and of the lifetime of the system which allows for forecasting of the imminent catastrophic failure.

Size scaling of failure strength with fat-tailed disorder in a fiber bundle model.

We investigate the size scaling of the macroscopic fracture strength of heterogeneous materials when microscopic disorder is controlled by fat-tailed distributions. We consider a fiber bundle model where the strength of single fibers is described by a power law distribution over a finite range. Tuning the amount of disorder by varying the power law exponent and the upper cutoff of fibers' strength, in the limit of equal load sharing an astonishing size effect is revealed: For small system sizes the bundle strength increases with the number of fibers, and the usual decreasing size effect of heterogeneous materials is restored only beyond a characteristic size. We show analytically that the extreme order statistics of fibers' strength is responsible for this peculiar behavior. Analyzing the results of computer simulations we deduce a scaling form which describes the dependence of the macroscopic strength of fiber bundles on the parameters of microscopic disorder over the entire range of system sizes.

Effect of disorder on shrinkage-induced fragmentation of a thin brittle layer.

We investigate the effect of the amount of disorder on the shrinkage-induced cracking of a thin brittle layer attached to a substrate. Based on a discrete element model we study how the dynamics of cracking and the size of fragments evolve when the amount of disorder is varied. In the model a thin layer is discretized on a random lattice of Voronoi polygons attached to a substrate. Two sources of disorder are considered: structural disorder captured by the local variation of the stiffness and strength disorder represented by the random strength of cohesive elements between polygons. Increasing the amount of strength disorder, our calculations reveal a transition from a cellular crack pattern, generated by the sequential branching and merging of cracks, to a disordered ensemble of cracks where the merging of randomly nucleated microcracks dominate. In the limit of low disorder, the statistics of fragment size is described by a log-normal distribution; however, in the limit of high disorder, a power-law distribution is obtained.

Fragmentation and shear band formation by slow compression of brittle porous media.

Localized fragmentation is an important phenomenon associated with the formation of shear bands and faults in granular media. It can be studied by empirical observation, by laboratory experiment, or by numerical simulation. Here we investigate the spatial structure and statistics of fragmentation using discrete element simulations of the strain-controlled uniaxial compression of cylindrical samples of different finite size. As the system approaches failure, damage localizes in a narrow shear band or synthetic fault "gouge" containing a large number of poorly sorted noncohesive fragments on a broad bandwidth of scales, with properties similar to those of natural and experimental faults. We determine the position and orientation of the central fault plane, the width of the shear band, and the spatial and mass distribution of fragments. The relative width of the shear band decreases as a power law of the system size, and the probability distribution of the angle of the central fault plane converges to around 30 degrees, representing an internal coefficient of friction of 0.7 or so. The mass of fragments is power law distributed, with an exponent that does not depend on scale, and is near that inferred for experimental and natural fault gouges. The fragments are in general angular, with a clear self-affine geometry. The consistency of this model with experimental and field results confirms the critical roles of preexisting heterogeneity, elastic interactions, and finite system size to grain size ratio on the development of shear bands and faults in porous media.

Blending stiffness and strength disorder can stabilize fracture.

Quasibrittle behavior, where macroscopic failure is preceded by stable damaging and intensive cracking activity, is a desired feature of materials because it makes fracture predictable. Based on a fiber-bundle model with global load sharing we show that blending strength and stiffness disorder of material elements leads to the stabilization of fracture, i.e., samples that are brittle when one source of disorder is present become quasibrittle as a consequence of blending. We derive a condition of quasibrittle behavior in terms of the joint distribution of the two sources of disorder. Breaking bursts have a power-law size distribution of exponent 5/2 without any crossover to a lower exponent when the amount of disorder is gradually decreased. The results have practical relevance for the design of materials to increase the safety of constructions.

Record-breaking events during the compressive failure of porous materials.

An accurate understanding of the interplay between random and deterministic processes in generating extreme events is of critical importance in many fields, from forecasting extreme meteorological events to the catastrophic failure of materials and in the Earth. Here we investigate the statistics of record-breaking events in the time series of crackling noise generated by local rupture events during the compressive failure of porous materials. The events are generated by computer simulations of the uniaxial compression of cylindrical samples in a discrete element model of sedimentary rocks that closely resemble those of real experiments. The number of records grows initially as a decelerating power law of the number of events, followed by an acceleration immediately prior to failure. The distribution of the size and lifetime of records are power laws with relatively low exponents. We demonstrate the existence of a characteristic record rank k(*), which separates the two regimes of the time evolution. Up to this rank deceleration occurs due to the effect of random disorder. Record breaking then accelerates towards macroscopic failure, when physical interactions leading to spatial and temporal correlations dominate the location and timing of local ruptures. The size distribution of records of different ranks has a universal form independent of the record rank. Subsequences of events that occur between consecutive records are characterized by a power-law size distribution, with an exponent which decreases as failure is approached. High-rank records are preceded by smaller events of increasing size and waiting time between consecutive events and they are followed by a relaxation process. As a reference, surrogate time series are generated by reshuffling the event times. The record statistics of the uncorrelated surrogates agrees very well with the corresponding predictions of independent identically distributed random variables, which confirms that temporal and spatial correlation in the crackling noise is responsible for the observed unique behavior. In principle the results could be used to improve forecasting of catastrophic failure events, if they can be observed reliably in real time.

Kinetic Monte Carlo algorithm for thermally induced breakdown of fiber bundles.

Fiber bundle models are one of the most fundamental modeling approaches for the investigation of the fracture of heterogeneous materials being able to capture a broad spectrum of damage mechanisms, loading conditions, and types of load sharing. In the framework of the fiber bundle model we introduce a kinetic Monte Carlo algorithm to investigate the thermally induced creep rupture of materials occurring under a constant external load. We demonstrate that the method overcomes several limitations of previous techniques and provides an efficient numerical framework at any load and temperature values. We show for both equal and localized load sharing that the computational time does not depend on the temperature; it is solely determined by the external load and the system size. In the limit of low load where the lifetime of the system diverges, the computational time saturates to a constant value. The method takes into account the secondary failures induced by subsequent load redistributions after breaking events, with the additional advantage that breaking avalanches always start from a single broken fiber.

Universality of fragment shapes.

The shape of fragments generated by the breakup of solids is central to a wide variety of problems ranging from the geomorphic evolution of boulders to the accumulation of space debris orbiting Earth. Although the statistics of the mass of fragments has been found to show a universal scaling behavior, the comprehensive characterization of fragment shapes still remained a fundamental challenge. We performed a thorough experimental study of the problem fragmenting various types of materials by slowly proceeding weathering and by rapid breakup due to explosion and hammering. We demonstrate that the shape of fragments obeys an astonishing universality having the same generic evolution with the fragment size irrespective of materials details and loading conditions. There exists a cutoff size below which fragments have an isotropic shape, however, as the size increases an exponential convergence is obtained to a unique elongated form. We show that a discrete stochastic model of fragmentation reproduces both the size and shape of fragments tuning only a single parameter which strengthens the general validity of the scaling laws. The dependence of the probability of the crack plan orientation on the linear extension of fragments proved to be essential for the shape selection mechanism.

Fractal frontiers of bursts and cracks in a fiber bundle model of creep rupture.

We investigate the geometrical structure of breaking bursts generated during the creep rupture of heterogeneous materials. Using a fiber bundle model with localized load sharing we show that bursts are compact geometrical objects; however, their external frontiers have a fractal structure which reflects their growth dynamics. The perimeter fractal dimension of bursts proved to have the universal value 1.25 independent of the external load and of the amount of disorder in the system. We conjecture that according to their geometrical features, breaking bursts fall in the universality class of loop-erased self-avoiding random walks with perimeter fractal dimension 5/4 similar to the avalanches of Abelian sand pile models. The fractal dimension of the growing crack front along which bursts occur proved to increase from 1 to 1.25 as bursts gradually cover the entire front.