Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fracturing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Fracture resistance of a TiB2 particle/SiC matrix composite at elevated temperature

The fracture resistance of a commercial TiB 2 particle/SiC matrix composite was evaluated at temperatures ranging from 20 to 1400 °C. A laser interferometric strain gauge (LISG) was used to continuously monitor the crack mouth opening displacement (CMOD) of the chevron-notched and straight-notched, three-point bend specimens used. Crack growth resistance curves (R-curves) were determined from the load versus displacement curves and displacement calibrations. Fracture toughness, work-of-fracture, and R-curve levels were found to decrease with increasing temperature. Microstructure, fracture surface, and oxidation coat were examined to explain the fracture behavior.

Fracture↗
DOE OSTI2018

An Overview of the EGS Collab Project: Field Validation of Coupled Process Modeling of Fracturing and Fluid Flow at the Sanford Underground Research Facility, Lead, SD

Enhanced Geothermal Systems (EGS) development will require an ability to accurately predict the flow rates and temperatures of the production wells. While simple in concept, complex heterogeneous fracture pathways can lead to channeling, short-circuiting, and premature thermal breakthrough complicating EGS. The EGS Collab project will establish a suite of intermediate-scale (~10-20 m) field test beds coupled with stimulation and interwell flow tests to provide a basis to better understand the fracture geometries and processes that control heat transfer between rock and stimulated fractures. These experiments will provide a means of testing tools, codes, and concepts that could later be employed under geothermal reservoir conditions at FORGE and in EGS. We will perform well-controlled, in situ experiments focused on rock fracture behavior and permeability enhancement. Pre- and post-test modeling of each test will allow for model prediction and validation. Comprehensive instrumentation will be used to collect high-quality, high-resolution geophysical and fracture characterization and fluid flow data. ,These data will be analyzed and compared with models and field observations to further elucidate the basic relationships between stress, induced seismicity, and permeability enhancement. To the maximum extent achievable, we will observe and quantify other key governing parameters that impact permeability, and attempt to understand how these parameters might change throughout the development and operation of an EGS project with the goal of enabling commercial viability of EGS. The Sanford Underground Research Facility (SURF) in South Dakota was selected as the EGS Collab project experimental site based on the evaluation of information from several sites. Our team has designed and constructed the first field experiment planned for this project, which is supported by the US Department of Energy’s Geothermal Technologies Office.

EGS, FORGE, field fracture stimulation and flow ex↗
DOE OSTI2025

A review of thermo-hydro-mechanical modeling of coupled processes in fractured rock: From continuum to discontinuum perspective

Coupled thermo-hydro-mechanical (THM) processes in fractured rock are playing a crucial role in geoscience and geoengineering applications. Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling. This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions. By comparing various models, one can better assess the uncertainties in predictions, particularly those related to the conceptual models. The review explores how these methodologies have significantly enhanced the fundamental understanding of how fractures respond to fluid injection and production, and improved predictive capabilities pertaining to coupled processes within fractured systems. It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience. The selection and calibration of model parameters should be based on typical ranges and applied to the specific conditions of applications. The challenges arising from inherent heterogeneity and uncertainties, nonlinear THM coupled processes, scale dependence, and computational limitations in representing field scale fractures are discussed. Realizing potential advances on computational capacity calls for methodical conceptualization, mathematical modeling, selection of numerical solution strategies, implementation, and calibration to foster simulation outcomes that intricately reflect the nuanced complexities of geological phenomena. Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical field of study.

Coupling scheme↗
DOE OSTI2016

Investigating the age distribution of fracture discharge using multiple environmental tracers, Bedrichov Tunnel, Czech Republic

We investigated the transit time distribution (TTD) of discharge collected from fractures in the Bedrichov Tunnel, Czech Republic, using lumped parameter models and multiple environmental tracers. We then utilize time series of δ 18 O, δ 2 H and 3 H along with CFC measurements from individual fractures in the Bedrichov Tunnel of the Czech Republic to investigate the TTD, and the uncertainty in estimated mean travel time in several fracture networks of varying length and discharge. We also compare several TTDs, including the dispersion distribution, the exponential distribution, and a developed TTD which includes the effects of matrix diffusion. The effect of seasonal recharge is explored by comparing several seasonal weighting functions to derive the historical recharge concentration. We identify best fit mean ages for each TTD by minimizing the error-weighted, multi-tracer χ2 residual for each seasonal weighting function. We use this methodology to test the ability of each TTD and seasonal input function to fit the observed tracer concentrations, and the effect of choosing different TTD and seasonal recharge functions on the mean age estimation. We find that the estimated mean transit time is a function of both the assumed TTD and seasonal weighting function. Best fits as measured by the χ2 value were achieved for the dispersion model using the seasonal input function developed here for two of the three modeled sites, while at the third site, equally good fits were achieved with the exponential model and the dispersion model and our seasonal input function. The average mean transit time for all TTDs and seasonal input functions converged to similar values at each location. The sensitivity of the estimated mean transit time to the seasonal weighting function was equal to that of the TTD. These results indicated that understanding seasonality of recharge is at least as important as the uncertainty in the flow path distribution in fracture networks and that unique identification of the TTD and mean transit time is difficult given the uncertainty in the recharge function. But, the mean transit time appears to be relatively robust to the structural model uncertainty. The results presented here should be applicable to other studies using environmental tracers to constrain flow and transport properties in fractured rock systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

TOUGH-UDEC: A simulator for coupled multiphase fluid flows, heat transfers and discontinuous deformations in fractured porous media

A numerical simulator entitled TOUGH-UDEC is presented for the assessment of coupled thermal-hydraulic-mechanical processes in fractured porous media. Two existing well-established codes, TOUGH2 and UDEC, are coupled to model multiphase fluid flows, heat transfers, and discontinuous deformations in fractured porous media by means of discrete fracture representation. TOUGH2 is widely used for the modeling of heat transfers and multiphase multicomponent fluid flows, and UDEC is a well-known distinct element code for rock mechanics. The two codes are solved sequentially, with coupling parameters passed to each equation at specific intervals. After solving thermal-hydraulic equations within the TOUGH2 code, pressure and temperature information is imported into the UDEC model. Upon solving the mechanical equation within the UDEC code the calculated fracture apertures are converted to the equivalent permeability and porosity values for a TOUGH2 flow analysis. The solution is calculated by iteratively following an explicit sequence for numerical efficiency. Verifications are presented to demonstrate the capabilities of the coupled TOUGH-UDEC simulator. Three application examples of (1) shear dilation due to increased pore pressure, (2) thermal stress and (3) CO 2 injection, show that the new simulator can be an effective tool for geoengineering applications involving shear activation of fractures and faults.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Effect of oxidative breakers on organic matter degradation, contaminant mobility and critical mineral release during shale-fracturing fluid interactions in the Marcellus Shale

Oil and gas production from organic-rich shale formations has become viable through advancements in multistage hydraulic fracturing. However, fluctuations in oil and gas prices, coupled with the sharp decline in gas production after the initial days of fracturing operations, drive operators to devise strategies for enhancing hydrocarbon production. The use of highly reactive fracturing fluids that include strong oxidizing agents, also known as breakers, can potentially increase well productivity. The oxidative breakers are used to reduce the viscosity of gel-based fluids after the proppant is transported into fracture zones of the target formation. These breakers can also degrade the organic matter, enhancing the release of hydrocarbons. However, chemical byproducts generated by the interaction of oxidative breakers with the shale matrix have not been extensively studied. Here, this study investigated the fluid-rock interactions between Marcellus Shale and synthetic hydraulic fracturing fluid (HFF) solutions comprising three different oxidative breakers, i.e., ammonium persulfate, sodium bromate, and sodium hypochlorite, commonly used in the Appalachian Basin, USA. Our results demonstrate that the type of oxidizing breakers used in the HFF controlled the type and amount of volatile organic compounds (VOCs) and semi-volatile organic compounds (SVOCs) produced. In all HFF reacted effluents, we observed the transformation of VOCs and the presence of organic acids in variable amounts. However, effluents from HFF containing sodium bromate and sodium hypochlorite breakers showed the presence of several halogenated organic compounds. Changes in major ions and mineralogy indicate that carbonate dissolution and barite precipitation were ubiquitous in all shale reacted effluents. Our results also demonstrate that the addition of oxidative breakers increased the concentration of several major and trace elements in the effluents. These elements fall under the critical mineral (CM) or critical element category due to their high demand in emerging technologies and susceptibility to supply chain disruption due to variety of factors. However, before oxidative breakers can be used at a larger field scale to enhance the release and recovery of CM and hydrocarbons, a better understanding is required of the potential environmental impacts associated with the generation and transformation of contaminants during breaker-fluid-rock interaction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Simulations of CO 2 injection into fractures and faults for improving their geophysical characterization at EGS sites

Here, we propose the use of CO 2 in push-pull well tests to improve geophysical identification and characterization of fractures and faults at enhanced geothermal system (EGS) sites. Using TOUGH2/ECO2N, we carried out numerical experiments of push-pull injection-production cycling of CO 2 into idealized vertical fractures and faults to produce pressure-saturation-temperature conditions that can be analyzed for their geophysical response. Our results show that there is a strong difference between injection and production mainly because of CO 2 buoyancy. While the CO 2 -plume grows laterally and upward during injection, not all CO 2 is recovered during the subsequent production phase. Even under the best conditions for recovery, at least 10% of the volume of the pores still remains filled with CO 2 . To improve EGS characterization, comparisons can be made of active seismic methods carried out before and after (time lapse mode) CO 2 injection into the fracture or fault. Here we find that across the CO 2 saturation range, C11 (the normal stiffness in the horizontal direction perpendicular to the fracture plane) varies between maximum and minimum values by about 15%. It reaches a maximum at around 6% gas saturation, decreasing exponentially to a minimum at higher saturations. Finally, our results suggest that CO 2 injection can be effectively used to infiltrate fault and fracture zones reaching about optimal saturation values in order to enhance seismic imaging at EGS sites.

58 GEOSCIENCES↗
DOE OSTI2017

Influence of injection-induced cooling on deviatoric stress and shear reactivation of preexisting fractures in Enhanced Geothermal Systems

Cold water injection into a hot, fractured, geothermal reservoir may trigger shear activation of pre-existing fractures that can help to enhance reservoir permeability, but may also result in unwanted seismicity. In this study, we investigate through numerical modeling of a hypothetical geothermal reservoir how injection-induced cooling may influence the potential for shear activation, paying special attention to the evolution of deviatoric stress under various stress regimes. In each case, we consider either a reservoir with homogeneous hydraulic properties or the presence of a more permeable fracture zone intersecting the injection well. This fracture zone is either oriented in the maximum (S Hmax ) or minimum (S hmin ) horizontal stress direction. Finally, our main finding is that depending on the configuration, injection-induced cooling stresses can favor or prevent shear reactivation of the preexisting fracture, and this effect can vary temporally and spatially.

58 GEOSCIENCES↗
DOE OSTI2019

Effects of Thermal Aging on the Fracture Toughness of Cast Stainless Steel CF8

Cast stainless steels are widely used in the primary coolant systems of nuclear power plants because of their high strength, toughness, and corrosion resistance. Mechanical property degradation in the steels, primarily due to a thermally-embrittled ferrite, is of major concern in an extended-term operation. To investigate the aging-induced loss of fracture resistance in cast stainless steels, static fracture resistance (J-R) testing was performed for three different CF8 cast stainless steels with significantly different d-ferrite contents. The materials were thermally-aged at 290-400°C for up to 10kh and fracture-tested at 25-400°C. It was found that the fracture toughness in high-ferrite materials generally decreased with thermal aging while the relationship between aging and fracture toughness in low-ferrite materials was rather complicated, showing non-monotonic variation with the degree of aging. While significant degradation was observed for the highly-aged high-ferrite materials, no material was observed having a K0.2mm less than 100 MPavm.

Collins, David A.↗
DOE OSTI2023

The effects of radiation-induced grain subdivision and dislocations on the fracture properties of uranium dioxide

Microcantilever bending tests were applied on uranium dioxide samples irradiated by 84 MeV Xe ions. Depending on the initial grain sizes, different levels of grain subdivision and dislocation development were observed near the irradiated surface. These radiation damage features were shown to degrade the fracture properties of both samples. The fundamental aspects of the radiation-microstructure-property relationships were discussed in this paper, taking a holistic consideration of the microcantilever bending data, fractography and microstructures. The separate effects of the irradiation-induced grain subdivision and dislocations on the fracture initiation and propagation of uranium dioxide were discussed. In conclusion, it was observed that radiation-induced dislocations and grain subdivision, without the presence of Xe bubbles co-located with the defects, cause reductions in the Young's modulus, fracture strain, and fracture stress, and a transition to intergranular fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2026

3D probabilistic fracture mechanics / computational fluid dynamics simulation of a reactor pressure vessel under transient conditions

Reactor pressure vessels (RPVs) are safety-critical light-water-reactor components that, under irradiation, experience long-term material degradation in the form of embrittlement. This can increase their susceptibility to fracture under thermal-shock conditions, which could occur during off-normal transients such as loss-of-coolant accidents (LOCAs). During a LOCA, the most severe conditions for the RPV occur when emergency core cooling water is injected through the cold legs into the water-and-steam-filled RPV. The rapid cooling of the downcomer and internal RPV surface causes decreased temperature and elevated thermally driven tensile stresses in the RPV wall. This, combined with long-term material embrittlement, may cause fracture initiation at pre-existing flaws, challenging the integrity of the RPV. Assessing RPV integrity during transients with a large spatial variation in the coolant temperature requires a modeling approach that considers the effects of spatially varying coolant temperature on the fracture probability of a population of flaws distributed throughout the RPV, accounting for spatially varying embrittlement. Here, the present study addresses this need by demonstrating first-of-its-kind coupling of a high-fidelity 3D computational fluid dynamics code with 3D probabilistic fracture mechanics This was accomplished using representative models of a pressurized-water reactor subjected to small- and medium-break LOCA conditions, both of which can result in large spatial temperature variations. While the observed impact of accounting for 3D effects was minimal under the small-break LOCA this study indicates a significant increase in the probability of fracture initiation under the medium-break LOCA when 3D effects are considered, relative to a spatially uniform cooling scenario.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Thicknesses of Chemically Altered Zones in Shale Matrices Resulting from Interactions with Hydraulic Fracturing Fluid

Hydraulic fracturing of unconventional shale reservoirs increases the fracture network surface area to access hydrocarbons from the low permeability rock matrix. Porosity and permeability of the matrix, through which hydrocarbons migrate to fractures, are important for determining production efficiency and can be altered by chemical interactions between shale and hydraulic fracturing fluids (HFFs). Here, we present results from an experimental study that characterizes the thickness of the alteration zone in the shale matrix after shale–HFF interactions. Experiments were conducted with whole cores submerged in HFF both with and without added barium and sulfate to promote barite scale formation. After 3 weeks of reaction at 77 bar and 80°C, the cores were characterized using X-ray microtomography, synchrotron X-ray fluorescence microprobe imaging, and synchrotron X-ray absorption spectroscopy. Our results show that the thickness of the altered zone depends on shale mineralogical composition and varies for different chemical reactions. Reactions between the low-carbonate Marcellus shale and HFF, pyrite (FeS 2 ) oxidation manifests as both a thick zone of sulfur oxidation (>0.5 cm) and a thinner zone of iron oxidation (100–150 $μ$m). Carbonate dissolution extended 100–200 μm into the matrix, with the resulting observable secondary porosity localized at the shale–fluid interface where mineral grains were removed by either dissolution or mechanical erosion. In solutions oversaturated with respect to barite, barite precipitates were observed in the reaction fluid and at the shal–HFF interface. In contrast, the carbonate dissolution zone in the high-carbonate Eagle Ford was only 30–40 μm thick, within which a uniform texture of increased porosity was observed. Pyrite oxidation in the Eagle Ford was evident from an iron oxidation zone (150–200 μm thick), while sulfur oxidation was minor and hard to observe. Barite precipitation extended 1–2 mm into the matrix when the initial HFF was oversaturated with respect to barite, filling shale microcracks down to the submicrometer length scale. Our findings provide a scientific basis to predict the extent of chemical alteration in shale reservoirs during hydraulic fracturing and its impacts on hydrocarbon production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A Critical Review of the Physicochemical Impacts of Water Chemistry on Shale in Hydraulic Fracturing Systems

Hydraulic fracturing of unconventional hydrocarbon resources involves the sequential injection of a high-pressure, particle-laden fluid with varying pH’s to make commercial production viable in low permeability rocks. This process both requires and produces extraordinary volumes of water. The water used for hydraulic fracturing is typically fresh, whereas “flowback” water is typically saline with a variety of additives which complicate safe disposal. As production operations continue to expand, there is an increasing interest in treating and reusing this high-salinity produced water for further fracturing. Here in this paper we review the relevant transport and geochemical properties of shales, and critically analyze the impact of water chemistry (including produced water) on these properties. We discuss five major geochemical mechanisms that are prominently involved in the temporal and spatial evolution of fractures during the stimulation and production phase: shale softening, mineral dissolution, mineral precipitation, fines migration, and wettability alteration. A higher salinity fluid creates both benefits and complications in controlling these mechanisms. For example, higher salinity fluid inhibits clay dispersion, but simultaneously requires more additives to achieve appropriate viscosity for proppant emplacement. In total this review highlights the nuances of enhanced hydrogeochemical shale stimulation in relation to the choice of fracturing fluid chemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Dispersion and Mixing in Three-Dimensional Discrete Fracture Networks: Nonlinear Interplay Between Structural and Hydraulic Heterogeneity: DISPERSION AND MIXING-REGULAR VERSUS DISORDERED NETWORKS

We investigate the relative impact of topological, geometric, and hydraulic heterogeneity on transport processes in three-dimensional fracture networks. Focusing on the two largest scales of heterogeneity in these systems, individual fracture and network structure, we compare transport through analogous structured and disordered three-dimensional fracture networks with varying degrees of hydraulic heterogeneity. For the moderate levels of hydraulic heterogeneity we consider, network structure is the dominant control of transport through the networks. Less dispersion, both longitudinal and transverse, is observed in structured networks than in disordered networks, due in part to the higher connectivity in the former, independent of the level of hydraulic heterogeneity. However, increases in dispersion with higher hydraulic heterogeneity are larger in the disordered networks than in the structured networks, thereby indicating that the interplay between structural and hydraulic heterogeneity is nonlinear. Here, we propose a measure of disorder in fracture networks by computing the Shannon entropy of the spectrum of the Laplacian of a weighted graph representation of the networks, where the weights are given by a combination of topological, geometric, and hydraulic properties. Finally, this metric, as a relative indicator by comparison between two networks, is a first approach to the dispersion potential and “mixing capacity” of a fracture network.

58 GEOSCIENCES↗
DOE OSTI2019

Cyclic Permeability Evolution During Repose Then Reactivation of Fractures and Faults

Abstract Cyclic growth and decay of permeability in fractures is explored during repeated reactivation and repose of saw‐cut fractures of Green River shale. These slide‐hold‐slide experiments are supported by measurements of fracture normal deformation and optical surface profilometry. Overall, we observe continuous permeability decay during repose (holds) and significant permeability enhancement during slow reactivation (slide). The permeability decay is accompanied by fault compaction. Both hydraulic aperture change (Δ b h ) and measured compaction (Δ b s ) are consistent with time‐dependent power law closure with a power exponent of ~0.2–0.4. These dual compaction magnitudes are positively correlated but Δ b h > Δ b s in late stage holds. Permeability enhancement during reactivation is typically also accompanied by fault dilation. However, we also observe some cases where hydraulic aperture change decouples from the measured deformation, conceivably driven by mobilization of wear products and influenced by the development of flow bottlenecks. Pretest and posttest surface profiles show that the surface topography of the fractures is planed down by shear removal. The shear removal is significant with initial laboratory prepared surface (~10 μm of aperture height) but less significant following consecutive reactivations (~2 μm). The flattened surfaces retain small‐scale, ~10–20 μm wavelength, roughness. Flow simulations, constrained by the surface topography and measured deformation, indicate that small‐scale roughness may control permeability at flow bottlenecks within a dominant flow channel. These results suggest cycles of permeability creation and destruction are an intrinsic component of the natural hydraulic system present in faults and fractures and provide an improved mechanistic understanding of the evolution of permeability during fault repose and reactivation.

Im, Kyungjae↗
DOE OSTI2020

Distributed Acoustic Sensing as a Distributed Hydraulic Sensor in Fractured Bedrock

Distributed acoustic sensing (DAS) was originally intended to measure oscillatory strain at frequencies of 1 Hz or more on a fiber optic cable. Recently, measurements at much lower frequencies have opened the possibility of using DAS as a dynamic strain sensor in boreholes. A fiber optic cable mechanically coupled to a geologic formation will strain in response to hydraulic stresses in pores and fractures. A DAS interrogator can measure dynamic strain in the borehole, which can be related to fluid pressure through the mechanical compliance properties of the formation. Because DAS makes distributed measurements, it is capable of both locating hydraulically active features and quantifying the fluid pressure in the formation. We present field experiments in which a fiber optic cable was mechanically coupled to two crystalline rock boreholes. The formation was stressed hydraulically at another well using alternating injection and pumping. The DAS instrument measured oscillating strain at the location of a fracture zone known to be hydraulically active. Rock displacements of less than 1 nm were measured. Laboratory experiments confirm that displacement is measured correctly. These results suggest that fiber optic cable embedded in geologic formations may be used to map hydraulic connections in three-dimensional fracture networks. A great advantage of this approach is that strain, an indirect measure of hydraulic stress, can be measured without beforehand knowledge of flowing fractures that intersect boreholes. The technology has obvious applications in water resources, geothermal energy, CO 2 sequestration, and remediation of groundwater in fractured bedrock.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Understanding Subsurface Fracture Evolution Dynamics Using Time‐Lapse Full Waveform Inversion of Continuous Active‐Source Seismic Monitoring Data

Abstract Predicting the behavior, geometry, and flow properties of subsurface fractures remains a challenging problem. Seismic models that can characterize fractures usually suffer from low spatiotemporal resolution. Here, we develop a correlative double‐difference time‐lapse full waveform inversion of continuous active source seismic monitoring data for determining high‐spatiotemporal‐resolution time‐lapse V p models of in‐situ fracture evolution at a shallow contamination site in Wyoming, USA. Assisted by rock physics modeling, we find that (a) rapidly increasing pore pressure initializes and grows the fracture, increasing the porosity slightly (from ∼13.7% to ∼14.6%) in the tight clay formation, thus decreasing V p (∼50 m/s); (b) the fluid injection continues decreasing V p , likely through the introduction of gas bubbles in the injectate; and (c) final V p reductions reach over ∼150 m/s due to a posited ∼4.5% gas saturation. Our results demonstrate that high‐resolution V p changes are indicative of mechanical and fluid changes within the fracture zone during hydrofracturing.

58 GEOSCIENCES↗
DOE OSTI2020

Vulnerable window of yield strength for swelling-driven fracture of phase-transforming battery materials

Abstract Despite numerous experimental and theoretical investigations of the mechanical behavior of high-capacity Si and Ge Li-ion battery anodes, our basic understanding of swelling-driven fracture in these materials remains limited. Existing theoretical studies have provided insights into elasto-plastic deformations caused by large volume change phase transformations, but have not modeled fracture explicitly beyond Griffith’s criterion. Here, we use a multi-physics phase-field approach to model self-consistently anisotropic phase transformation, elasto-plastic deformation, and crack initiation and propagation during lithiation of Si nanopillars. Our computational results reveal that fracture occurs within a “vulnerable window” inside the two-dimensional parameter space of yield strength and fracture energy and highlight the importance of taking into account the surface localization of plastic deformation to accurately predict the magnitude of tensile stresses at the onset of fracture. They further demonstrate how the increased robustness of hollow nanopillars can be understood as a direct effect of anode geometry on the size of this vulnerable window. Those insights provide an improved theoretical basis for designing next-generation mechanically stable phase-transforming battery materials undergoing large volume changes.

25 ENERGY STORAGE↗