Solid-oxide iron-air batteries are an emerging technology for large-scale energy storage, but mechanical degradation of Fe-based storage materials limits battery lifetime. Experimental studies have revealed cycling degradation due to large volume changes during oxidation/reduction (via H2O/H2at 800 °C), but degradation has not yet been correlated with the microstructural stress and strain evolution. Here, we implement a finite element model for oxidation of a Fe lamella to FeO (74% volumetric expansion), in a lamellar Fe foam designed for battery applications. Growth of FeO at the Fe/gas interface is coupled, via an oxidation reaction and solid-state diffusion, with the shrinkage rate of the Fe lamellar core. Using isotropic linear elasticity and plastic hardening, the model simulates deformation of a continuously growing FeO layer by dynamically switching “gas” elements into new “FeO” elements along a sharp FeO/gas interface. As oxidation progresses, the effective plastic strain and von Mises stress increase in FeO. Distribution of tensile and compressive stresses along the Fe/FeO interface are validated by oxidation theory and explain interface delamination, as observed during in operando X-ray tomography experiments. The model explains the superior stability of lamellar vs dendritic foam architectures and the improved redox lifetime of Fe-Ni foams.
- Award ID(s):
- 1914780
- NSF-PAR ID:
- 10336145
- Date Published:
- Journal Name:
- Acta Geotechnica
- ISSN:
- 1861-1125
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
more » « less
Abstract This article introduces a computational design framework for obtaining three‐dimensional (3D) periodic elastoplastic architected materials with enhanced performance, subject to uniaxial or shear strain. A nonlinear finite element model accounting for plastic deformation is developed, where a Lagrange multiplier approach is utilized to impose periodicity constraints. The analysis assumes that the material obeys a von Mises plasticity model with linear isotropic hardening. The finite element model is combined with a corresponding path‐dependent adjoint sensitivity formulation, which is derived analytically. The optimization problem is parametrized using the solid isotropic material penalization method. Designs are optimized for either end compliance or toughness for a given prescribed displacement. Such a framework results in producing materials with enhanced performance through much better utilization of an elastoplastic material. Several 3D examples are used to demonstrate the effectiveness of the mathematical framework.
-
more » « less
Summary This study presents a thermo‐hydro‐mechanical (THM) model of unsaturated soils using isogeometric analysis (IGA). The framework employs Bézier extraction to connect IGA to the conventional finite element analysis (FEA), featuring the current study as one of the first attempts to develop an IGA‐FEA framework for solving THM problems in unsaturated soils. IGA offers higher levels of interelement continuity making it an attractive method for solving highly nonlinear problems. The governing equations of linear momentum, mass, and energy balance are coupled based on the averaging procedure within the hybrid mixture theory. The Drucker‐Prager yield surface is used to limit the modified effective stress where the model follows small strain, quasi‐static loading conditions. Temperature dependency of the surface tension is implemented in the soil‐water retention curve. Nonuniform rational B‐splines (NURBS) basis functions are used in the standard Galerkin method and weak formulations of the balance equations. Displacement, capillary pressure, gas pressure, and temperature are four independent quantities that are approximated by NURBS in spatial discretization. The framework is used to simulate strain localization in an undrained dense sand subjected to plane strain biaxial compression under different temperatures and displacement velocities. Results show that an increase in the displacement rate leads to reduction in the equivalent plastic strain while an increase in the temperature leads to an increase in the equivalent plastic strain. The findings suggest that the proposed IGA‐based framework offers a viable alternative for solving THM problems in unsaturated soils.
-
null (Ed.)SUMMARY The yield surfaces of rocks keep evolving beyond the initial yield stress owing to the damage accumulation and porosity change during brittle deformation. Using a poroelastic damage rheology model, we demonstrate that the measure of coupling between the yield surface change and accumulated damage is correlated with strain localization and the Kaiser effect. Constant or minor yield surface change is associated with strong strain localization, as seen in low-porosity crystalline rocks. In contrast, strong coupling between damage growth and the yield surface leads to distributed deformation, as seen in high-porosity rocks. Assuming that during brittle deformation damage occurs primarily in the form of microcracks, we propose that the measured acoustic emission (AE) in rock samples correlates with the damage accumulation. This allows quantifying the Kaiser effect under cyclic loading by matching between the onset of AE and the onset of damage growth. The ratio of the stress at the onset of AE to the peak stress of the previous loading cycle, or Felicity Ratio (FR), is calculated for different model parameters. The results of the simulation show that FR gradually decreases in the case of weak coupling between yield surface and damage growth. For a strong damage-related coupling promoting significant yield surface change, the FR remains close to one and decreases only towards the failure. The model predicts that a steep decrease in FR is associated with a transition between distributed and localized modes of failure. By linking the evolving yield surface to strain localization patterns and the Kaiser effect, the poroelastic damage rheology model provides a new quantitative tool to study failure modes of brittle rocks.more » « less
-
Abstract Many low-angle normal faults (dip ≤30°) accommodate tens of kilometers of crustal extension, but their mechanics remain contentious. Most models for low-angle normal fault slip assume vertical maximum principal stress σ1, leading many authors to conclude that low-angle normal faults are poorly oriented in the stress field (≥60° from σ1) and weak (low friction). In contrast, models for low-angle normal fault formation in isotropic rocks typically assume Coulomb failure and require inclined σ1 (no misorientation). Here, a data-based, mechanical-tectonic model is presented for formation of the Whipple detachment fault, southeastern California. The model honors local and regional geologic and tectonic history and laboratory friction measurements. The Whipple detachment fault formed progressively in the brittle-plastic transition by linking of “minidetachments,” which are small-scale analogs (meters to kilometers in length) in the upper footwall. Minidetachments followed mylonitic anisotropy along planes of maximum shear stress (45° from the maximum principal stress), not Coulomb fractures. They evolved from mylonitic flow to cataclasis and frictional slip at 300–400 °C and ∼9.5 km depth, while fluid pressure fell from lithostatic to hydrostatic levels. Minidetachment friction was presumably high (0.6–0.85), based upon formation of quartzofeldspathic cataclasite and pseudotachylyte. Similar mechanics are inferred for both the minidetachments and the Whipple detachment fault, driven by high differential stress (∼150–160 MPa). A Mohr construction is presented with the fault dip as the main free parameter. Using “Byerlee friction” (0.6–0.85) on the minidetachments and the Whipple detachment fault, and internal friction (1.0–1.7) on newly formed Reidel shears, the initial fault dips are calculated at 16°–26°, with σ1 plunging ∼61°–71° northeast. Linked minidetachments probably were not well aligned, and slip on the evolving Whipple detachment fault probably contributed to fault smoothing, by off-fault fracturing and cataclasis, and to formation of the fault core and fractured damage zone. Stress rotation may have occurred only within the mylonitic shear zone, but asymmetric tectonic forces applied to the brittle crust probably caused gradual rotation of σ1 above it as a result of: (1) the upward force applied to the base of marginal North America by buoyant asthenosphere upwelling into an opening slab-free window and/or (2) basal, top-to-the-NE shear traction due to midcrustal mylonitic flow during tectonic exhumation of the Orocopia Schist. The mechanical-tectonic model probably applies directly to low-angle normal faults of the lower Colorado River extensional corridor, and aspects of the model (e.g., significance of anisotropy, stress rotation) likely apply to formation of other strong low-angle normal faults.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s11440-022-01605-6
