Mechanically driven glassy systems and complex fluids exhibit a wealth of rheological behaviors that call for theoretical understanding and predictive modeling. A distinct feature of these nonequilibrium systems is their dynamically evolving state of structural disorder, which determines their rheological responses. Here we highlight a recently developed nonequilibrium thermodynamic framework in which the structural state is characterized by an evolving effective disorder temperature that may differ from the ordinary thermal temperature. The specific properties of each physical system of interest are described by a small set of coarse-grained internal state variables and their associated energies and entropies. The dynamics of the internal variables, together with the flow of energy and entropy between the different parts of the driven system, determine continuum-level rheological constitutive laws. We conclude with brief descriptions of several successful applications of this framework. © 2013 The Royal Society of Chemistry.

Quantitative understanding of the fracture toughness of metallic glasses, including the associated ductile-to-brittle transitions, is not yet available. Here we use a simple model of plastic deformation in glasses, coupled to an advanced Eulerian level set formulation for solving complex free boundary problems, to calculate the fracture toughness of metallic glasses as a function of the degree of structural relaxation corresponding to different annealing times near the glass temperature. Our main result indicates the existence of an elasto-plastic crack tip instability for sufficiently relaxed glasses, resulting in a marked drop in the toughness, which we interpret as a ductile-to-brittle transition similar to experimental observations.

Understanding the fracture toughness (resistance) of glasses is a fundamental problem of prime theoretical and practical importance. Here we theoretically study its dependence on the loading rate, the age (history) of the glass, and the notch radius ρ. Reduced-dimensionality analysis suggests that the notch fracture toughness results from a competition between the initial, age- and history-dependent, plastic relaxation time scale τ0pl and an effective loading time scale τext(KI,ρ), where KI is the tensile stress-intensity-factor rate. The toughness is predicted to scale with ρ independently of ξ≡τext/τ0pl for ξ1, to scale as Tρlog(ξ) for ξ1 (related to thermal activation, where T is the temperature), and to feature a nonmonotonic behavior in the crossover region ξ∼O(1) (related to plastic yielding dynamics). These predictions are verified using 2D computations, providing a unified picture of the notch fracture toughness of glasses. The theory highlights the importance of time-scale competition and far-from-steady-state elasto-viscoplastic dynamics for understanding the toughness and shows that the latter varies quite significantly with the glass age (history) and applied loading rate. Experimental support for bulk metallic glasses is presented, and possible implications for applications are discussed.

A well-established numerical approach to solve the Navier-Stokes equations for incompressible fluids is Chorin's projection method [1], whereby the fluid velocity is explicitly updated, and then an elliptic problem for the pressure is solved, which is used to orthogonally project the velocity field to maintain the incompressibility constraint. In this paper, we develop a mathematical correspondence between Newtonian fluids in the incompressible limit and hypo-elastoplastic solids in the slow, quasi-static limit. Using this correspondence, we formulate a new fixed-grid, Eulerian numerical method for simulating quasi-static hypo-elastoplastic solids, whereby the stress is explicitly updated, and then an elliptic problem for the velocity is solved, which is used to orthogonally project the stress to maintain the quasi-staticity constraint. We develop a finite-difference implementation of the method and apply it to an elasto-viscoplastic model of a bulk metallic glass based on the shear transformation zone theory. We show that in a two-dimensional plane strain simple shear simulation, the method is in quantitative agreement with an explicit method. Like the fluid projection method, it is efficient and numerically robust, making it practical for a wide variety of applications. We also demonstrate that the method can be extended to simulate objects with evolving boundaries. We highlight a number of correspondences between incompressible fluid mechanics and quasi-static elastoplasticity, creating possibilities for translating other numerical methods between the two classes of physical problems.

Identifying heterogeneous structures in glasses-such as localized soft spots-and understanding structure-dynamics relations in these systems remain major scientific challenges. Here, we derive an exact expression for the local thermal energy of interacting particles (the mean local potential energy change caused by thermal fluctuations) in glassy systems by a systematic low-temperature expansion. We show that the local thermal energy can attain anomalously large values, inversely related to the degree of softness of localized structures in a glass, determined by a coupling between internal stresses-an intrinsic signature of glassy frustration-anharmonicity and low-frequency vibrational modes. These anomalously large values follow a fat-tailed distribution, with a universal exponent related to the recently observed universal [Formula: see text] density of states of quasilocalized low-frequency vibrational modes. When the spatial thermal energy field-a "softness field"-is considered, this power law tail manifests itself by highly localized spots, which are significantly softer than their surroundings. These soft spots are shown to be susceptible to plastic rearrangements under external driving forces, having predictive powers that surpass those of the normal modes-based approach. These results offer a general, system/model-independent, physical/observable-based approach to identify structural properties of quiescent glasses and relate them to glassy dynamics.