Search NASA⌕ Search

SEARCH · Search NASA

Results for “bulk fracture”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Petrology and Geochemistry of New Ureilites and Ureilite Genesis

Ureilites are C-bearing, basalt-depleted olivine+pyroxene achondrites from a differentiated asteroid. The group is heterogeneous, exhibiting ranges in O isotopic composition, Fe/Mg, Fe/Mn, pyx/ol, siderophile and lithophile trace element content, and C content and isotopic composition [1]. Some of these characteristics are nebular in origin; others were strongly overprinted by asteroidal igneous processes. The consensus view is that most ureilites are melt-residues, but some are partial cumulates or have interacted with a melt [1,2]. An "unroofing" event occurred while the parent asteroid was hot that froze in mineral core com-positions and resulted in FeO reduction at olivine grain margins. We have studied several new ureilites, but will focus here on two anomalous stones; LAR 04315 and NWA 1241. LAR 04315 is texturally unusual. It contains olivine with angular subdomains, and low-Ca pyroxene riddled with wormy inclusions of metal+troilite, graphite, and possibly other phases, and irregular inclusions of high-Ca pyroxene. Reduction occurred along olivine grain margins and internal fractures, but not along subdomain boundaries. Although texturally odd, LAR 04351 is a typical ureilite in mineral and bulk composition. The olivine is Fo80.8 and falls on the ureilite Fe/Mn-Fe/Mg trend. Its olivine composition falls within the range of the majority of ureilites, and it is typical of these ureilites in bulk rock lithophile and siderophile element contents.

Mittlefehldt, David W.↗
NASA NTRS2014

Identification of New Lithic Clasts in Lunar Breccia 14305 by Micro-CT and Micro-XRF Analysis

From 1969 to 1972, Apollo astronauts collected 382 kg of rocks, soils, and core samples from six locations on the surface of the Moon. The samples were initially characterized, largely by binocular examination, in a custom-built facility at Johnson Space Center (JSC), and the samples have been curated at JSC ever since. Despite over 40 years of study, demand for samples remains high (~500 subsamples per year are allocated to scientists around the world), particularly for plutonic (e.g., anorthosites, norites, etc.) and evolved (e.g., granites, KREEP basalts) lithologies. The reason for the prolonged interest is that as new scientists and new techniques examine the samples, our understanding of how the Moon, Earth, and other inner Solar System bodies formed and evolved continues to grow. Scientists continually clamor for new samples to test their emerging hypotheses. Although all of the large Apollo samples that are igneous rocks have been classified, many Apollo samples are complex polymict breccias that have previously yielded large (cm-sized) igneous clasts. In this work we present the initial efforts to use the non-destructive techniques of micro-computed tomography (micro-CT) and micro x-ray fluorescence (micro-XRF) to identify large lithic clasts in Apollo 14 polymict breccia sample 14305. The sample in this study is 14305,483, a 150 g slab of regolith breccia 14305 measuring 10x6x2 cm (Figure 1a). The sample was scanned at the University of Texas High-Resolution X-ray CT Facility on an Xradia MicroXCT scanner. Two adjacent overlapping volumes were acquired at 49.2 micrometer resolution and stitched together, resulting in 1766 slices. Each volume was acquired at 100 kV accelerating voltage and 98 mA beam current with a 1 mm CaF2 filter, with 2161 views gathered over 360deg at 3 seconds acquisition time per view. Micro-XRF analyses were done at Washington University in St. Louis, Missouri on an EDAX Orbis PC micro-XRF instrument. Multiple scans were made at 40 kV accelerating voltage, 800 mA beam current, 30 μm beam diameter, and a beam spacing of 30-120 micrometer. The micro-CT scan of 14305,483 (Figure 2) was able to identify several large lithic clasts (approx. 1 cm) within the interior of the slab. These clasts will be exposed by band-sawing or chipping of the slab, and their composition more fully characterized by subsequent micro-XRF analysis. In addition to lithic clasts, the micro-CT scans identified numerous mineral clasts, including many FeNi metal grains, as well as the prominent fractures within the slab. The micro-XRF analyses (Figure 1b,c) of the slab surfaces revealed the bulk chemical compositions (qualitative) of the different clast types observed. In particular, by looking at the ratios of major elements (e.g. Ca:Mg:Fe), differences among the many observed clast types are readily observed. Moreover, several clasts not apparent to the naked eye were revealed in the K:Al:Si ratio map. The scans are also able to identify small grains of Zr- and P-rich minerals (not shown), which could in turn yield important age dates for the samples.

Zeigler, Ryan A.↗
DOE OSTI2024

Carbon Storage Core Characterization Efforts at NETL

The multi-scale Computed Tomography (CT) and core flow facility in the Geocharacterization Laboratory at NETL, Morgantown yields porosity, permeability, and fracture properties of rock core samples obtained from the subsurface while maintaining the integrity of the sample. Additionally, geophysical bulk rock properties are analyzed with the laboratory’s GeoTEK multi-sensor core logger in a comparable fashion to downhole methods. NETL researchers collaborate with stakeholders within the carbon storage, oil and gas, and critical minerals sectors. Since 2017, over 1.88 miles of core have been analyzed within the laboratory and all data is publicly available through the Technical Report Series (TRS) on the Energy Data eXchange (EDX). Additionally, the website, RokBase, was curated to extrapolate and visualize the high-resolution data from field operations. The characterization of the Lively Grove #1 (LG#1) well provides a case study into the full capabilities of the Geocharacterization Laboratory. During the comprehensive study of LG#1 ~1-2 mm in diameter, vertical to bedding, cylindrical structures were identified throughout the St. Peter Formation. These structures are pervasive throughout the St. Peter Formation at depth and are characterized as the trace fossil, Skolithos.

Isom, Shelby L↗
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
DOE OSTI2023

Dual Recycling of Depolymerization Catalyst and Biodegradable Polyester that Markedly Outperforms Polyolefins

Abstract Chemically recyclable, circular polymers continue to attract increasing attention, but rendering both catalysts for depolymerization and high‐performance polymers recyclable is a more sustainable yet challenging goal. Here we introduce a dual catalyst/polymer recycling system in that recyclable inorganic phosphomolybdic acid catalyzes selective depolymerization of high‐ceiling‐temperature biodegradable poly(δ‐valerolactone) in bulk phase, which, upon reaching suitable molecular weight, exhibits outstanding mechanical performance with a high tensile strength of ≈66.6 MPa, fracture strain of ≈904 %, and toughness of ≈308 MJ m −3 , and thus markedly outperforms commodity polyolefins, recovering its monomer in pure state and quantitative yield at only 100 °C. In sharp contrast, the uncatalyzed depolymerization not only requires a high temperature of >310 °C but is also low yielding and non‐selective. Importantly, the recovered monomer can be repolymerized as is to reproduce the same polymer, thereby closing the circular loop, and the recycled catalyst can be reused repeatedly for depolymerization runs without loss of its catalytic activity and efficiency.

Li, Xin‐Lei↗
DOE OSTI2023

Dual Recycling of Depolymerization Catalyst and Biodegradable Polyester that Markedly Outperforms Polyolefins

Chemically recyclable, circular polymers continue to attract increasing attention, but rendering both catalysts for depolymerization and high-performance polymers recyclable is a more sustainable yet challenging goal. Here we introduce a dual catalyst/polymer recycling system in that recyclable inorganic phosphomolybdic acid catalyzes selective depolymerization of high-ceiling-temperature biodegradable poly(δ-valerolactone) in bulk phase, which, upon reaching suitable molecular weight, exhibits outstanding mechanical performance with a high tensile strength of ≈66.6 MPa, fracture strain of ≈904 %, and toughness of ≈308 MJ m –3 , and thus markedly outperforms commodity polyolefins, recovering its monomer in pure state and quantitative yield at only 100 °C. In sharp contrast, the uncatalyzed depolymerization not only requires a high temperature of >310 °C but is also low yielding and non-selective. Importantly, the recovered monomer can be repolymerized as is to reproduce the same polymer, thereby closing the circular loop, and the recycled catalyst can be reused repeatedly for depolymerization runs without loss of its catalytic activity and efficiency.

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

The deformation and burst behavior of Zircaloy-4 cladding tubes with hydride rim features subject to internal pressure loads

This study is motivated by the desire to conduct a separate effects test replicating specific aspects of an in-pile pulse reactor test on high exposure light water reactor (LWR) fuel rods. The principal purpose of the tests was to determine whether temperature and or hydrogen concentration thresholds exist in assessing the vulnerability of a stress relieved zirconium alloy cladding tubes to low strain ruptures. A secondary objective was to determine whether observations of the cladding stress during deformation could aid in the understanding of the crack formation and propagation in cladding tubes with hydride rim features. Here, the study utilized artificially hydrided Zircaloy-4 cladding tubes and an isothermal, internal pressure test that was controlled with an in-situ strain measurement. The test results show a clear transition in failure vulnerability between room temperature and 150 °C. Below 150 °C, cladding tubes rupture at very low uniform plastic strain—even at relatively modest cladding hydrogen levels. At and above 150 °C, cladding tubes only rupture at low uniform plastic strain if the cladding hydrogen content is greater than 500 ppm and usually experience at least some modest amount of plastic deformation prior to rupture. Post-test metallographs from non-ruptured cladding tubes show numerous primary cracks forming in the hydride rim and arresting at the metal/hydride interface. This observation combined with the decrease in bulk yielding of hydrided cladding tubes when compared to fresh non-hydrided tubes leads to the conclusion that the fracture stress of the hydride rims is related to and bounded by the yield stress of the zirconium metal. Some evidence is presented for a change in crack propagation mechanism from one of void formation on hydride plates and linkage at temperatures below 150 °C to more classical ductile shearing beneath the primary crack tip at temperatures above 150 °C. This apparent change in crack propagation mechanism may help explain the observed temperature and hydrogen concentration thresholds.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2021

Molecular dynamics simulations of ultrafast radiation induced melting at metal–semiconductor interfaces

Metal–semiconductor contacts in silicon carbide (SiC) diodes endure damages at the interface when exposed to harsh radiation environments. Due to the rapid rise in temperature and ultrafast cooling that follows the radiation impact, the structural properties of the materials can be altered through melting, recrystallization, and amorphization. A detailed understanding of the material failure modes at the interface is lacking, specifically at the nanoscale. We use molecular simulations to investigate the ultrafast melting at tungsten (W)–SiC interfaces following radiation damage and apply deep learning techniques to track the transient evolution of the local molecular structures. We show that W near the radiation track undergoes melting and, eventually, most of it recrystallizes with a noticeable degree of undercooling, while SiC is rendered permanently amorphous. The observation of local undercooling in W films is important as it can affect the device performance even before the bulk melting temperature of the material is reached. We also show that at high temperatures, the interface undergoes a fracture-like failure. The results presented here are significant in understating the different failure modes of SiC diode materials.

Ashwin Ravichandran↗
DOE OSTI2022

Exploring the Promise of Multifunctional “Zintl-Phase-Forming” Electrolytes for Si-Based Full Cells

Li-M-Si ternary Zintl phases have gained attention recently due to their high structural stability, which can improve the cycling stability compared to a bulk Si electrode. Adding multivalent cation salts (such as Mg 2+ and Ca 2+ ) in the electrolyte was proven to be a simple way to form Li-M-Si ternary phases in situ in Si-based Li-ion cells. To explore the promise of Zintl-phaseforming electrolytes, we systematically investigated their application in pouch cells via electrochemical and multiscale postmortem analysis. The introduction of multivalent cations, such as Mg 2+ , during charging can form Li x M y Si ternary phases. They can stabilize Si anions and reduce side reactions with electrolyte, improving the bulk stability. More importantly, Mg 2+ and Ca 2+ incorporate into interfacial side reactions and generate inorganic-rich solid- electrolyte interphase, thus enhancing the interfacial stability. Therefore, the full cells with Zintl-phase-forming electrolytes achieve higher capacity retentions at the C/3 rate after 100 cycles, compared to a baseline electrolyte. Additionally, strategies for mitigating the electrode-level fractures of Si were evaluated to make the best use of Zintl-phase-forming electrolytes. In conclusion, this work highlights the significance of synergistic impact of multifunctional additives to stabilize both bulk and interface chemistry in high-energy Si anode materials for Li-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2025

Dynamic disulfide bond networks enable self-healable and mechanically resilient intrinsically stretchable organic solar cells

A dynamic disulfide network introduced into donor/acceptor blends enables room-temperature self-healing and mechanical resilience in intrinsically stretchable organic solar cells, achieving performance recovery after high mechanical strain. The development of intrinsically stretchable organic solar cells (IS-OSCs) faces significant challenges in balancing mechanical durability and optoelectronic performance. Conventional π-conjugated polymer-based donor/acceptor blend films often exhibit limited stretchability and irreversible performance degradation under mechanical strain. To address these limitations, we propose a novel self-healable donor/acceptor blended film with a dual-network morphology, achieved by incorporating a dynamic disulfide bond-based crosslinked network into the bulk-heterojunction film. The resulting thin films demonstrate a power conversion efficiency (PCE) of 16.39% in rigid OSC devices and a fracture strain of 15.6%. Remarkably, the IS-OSCs retain 80% of their initial PCE under 30% strain and exhibit performance recovery after multiple stretch-release cycles at 40% strain through a room-temperature self-healing process. This work provides a proof-of-concept for highly stretchable and durable IS-OSCs, offering valuable insights for advancing the field of wearable energy systems, adaptive solar textiles, and sustainable electronics.

Yang, Wenyu↗
NASA NTRS1980

Tribological properties of silicon carbide in metal removal process

Material properties are considered as they relate to adhesion, friction, and wear of single crystal silicon carbide in contact with metals and alloys that are likely to be involved in a metal removal process such as grinding. Metal removal from adhesion between sliding surfaces in contact and metal removal as a result of the silicon carbide sliding against a metal, indenting into it, and plowing a series of grooves or furrows are discussed. Fracture and deformation characteristics of the silicon carbide surface are also covered. The adhesion, friction, and metal transfer to silicon carbide is related to the relative chemical activity of the metals. The more active the metal, the higher the adhesion and friction, and the greater the metal transfer to silicon carbide. Atomic size and content of alloying elements play a dominant role in controlling adhesion, friction, and abrasive wear properties of alloys. The friction and abrasive wear (metal removal) decrease linearly as the shear strength of the bulk metal increases. They decrease as the solute to solvent atomic radius ratio increases or decreases linearly from unity, and with an increase of solute content. The surface fracture of silicon carbide is due to cleavages of 0001, 10(-1)0, and/or 11(-2)0 planes.

Miyoshi, K.↗
NASA NTRSDate not supplied

Repurposing Drilling Control Diagnostics for Subsurface Edge Detection and Boundary Advisement During Planetary Drilling

Informed decision-making during lunar drilling and sampling missions will require data monitoring tools and specialized ground data systems. Accurate and updated situational awareness, with ongoing data monitoring, is critical for timely responses by to incoming science data. Traverse plans and scheduled activities may need to be flexibly changed in order to react to unexpected data or situations. Unlike (for example) Mars missions, the relative lightspeed closeness of the Moon allows for near-real-time ground processing of incoming mission and instrument data. An Apollo-class lunar regolith drill will in a sense “travel” a meter or two vertically at a given subsurface characterization site. As the drill penetrates into lunar regolith, it is likely to encounter a range of material densities, orientations, fracture toughness, and (perhaps) ice percentages. Lunar drill telemetry can provide science teams with a valuable first look into the subsurface structure, the regolith bulk properties, and constituents at each drilled site. Real-time AI-based recognition and reaction to downhole situations has been developed for automated deeper drilling on Mars and beyond. We can leverage the same knowledge bases and pattern-matching as areal-time interpreter of the subsurface, a situational awareness tool during drilling operations. We recently (Sept. 2019) demonstrated this AI drilling monitoring and analysis capability, in control of in-situ drilling and sampling operations, mounted on a KREX-2 rover in Chile’s Atacama Desert. Terrestrial automated drilling log analyses in oil exploration have used similar machine learning techniques in classifying and identifying features in drilling logs –but these typically are designed assuming a drilling fluid influencing downhole measurements and data (permeability, resistivity). Drilling models and existing AI software designed to detect and respond to drilling faults and hard materials can be repurposed, for near-real-time (ground-based) interpretation of drilling telemetry –a potentially valuable advisory tool for strata boundaries and changes in drilling parameters. On the Moon, this approach could be used to study the structure and to some extent the composition of lunar regolith vs. borehole depth, based on recognizable variations in fracture hardness, drilling energy and penetration rates while actively drilling. Since the early 2000s, a series of increasingly-capable real-time drilling telemetry interpretation and characterization software tools have been developed. These subsurface models and software tools have monitored the real-time drilling data received, and automatically identified changes in drill behavior (e.g., encountering a harder target layer, bit inclusions, drill choking due to infall downhole, and others) correlating these with subsurface structures and features. We discuss the mappings between drill borehole parameters, faults or events detected, and modeled changes in rock layer boundaries, in examples drawn from field testing at analog sites in an Arctic impact crater, Rio Tinto, and Chile’s Atacama Desert. These demonstrate how subsurface structural boundaries led to fault detections and responses by the software.

drilling advisor↗
NASA NTRS1981

Microstructure and mechanical properties of bulk yttria-partially-stabilized zirconia

A commercially available bulk 4.5 mole percent yttria-(Y2O3) partially stabilized zirconia (PSZ) was studied by light microscopy, X-ray analysis, microhardness measurement, and fracture toughness testing. The growth of the precipitates and the phase transformations were studied as a function of aging in air at 1500 C. Aging curves were constructed for both the as received and the solution annealed and quenched materials; the curves showed hardness peaks at 1397 and 1517 Kg/sq mm respectively. The rectangular plate shaped tetragonal precipitates were found to have a 110 habit plane. A total of twelve different types of tetragonal precipitates were found. Grinding of the Y2O3 PSZ into powder did not cause a significant amount of metastable tetragonal precipitates to transform into the monoclinc phase, thus indicating that transformation toughening is not a significant mechanism for the material.

Valentine, P. G.↗
NASA NTRS1991

Development of high Tc (greater than 110K) Bi, Tl and Y-based materials as superconducting circuit elements

This report is presented in two parts. Part 1 deals primarily with Bi-based materials and a small amount of work on a Y-based composition while Part 2 covers work on Tl-based materials. In Part 1, a reliable and reproducible process for producing bulk bismuth-based superconductors has been developed. It is noted however, that a percentage of the tapecast material experiences curling and fracturing after a 30 hour sintering period and is thus in need of further examination. The Bi-Sr-Ca-Cu-O (BSCCO) material has been characterized by critical temperature data, X-ray diffraction data, and surface morphology. In the case of T sub c, it is not critical to anneal the material. It appears that the BSCCO material has the possibility of producing a better grounding strap than that of the 123 material. Attempts to reproduce near room temperature superconductors in the Y-Ba-Cu-O system were unsuccessful. In Part 2, several methods of processing the high temperature superconductor Tl2Ba2Ca2Cu3O10 were investigated; i.e., different precursor compositions were sintered at various sintering times and temperatures. The highest superconductig temperature was found to be 117.8K when fired at 900 C for three hours. Higher sintering temperatures produced a melted sample which was nonsuperconducting at liquid nitrogen temperature. Also, a preliminary study found Li2O substitutions for copper appeared to increase the transition temperature and create fluxing action upon sintering. It was suggested that lower sintering temperatures might be obtained with lithium additions to produce reliable Tl2Ba2Ca2Cu3O10 processing methods.

Haertling, Gene H.↗
NASA NTRS2001

The Porosity of 433 Eros

Data from the NEAR mission show the bulk density of 433 Eros is 2.67 g/cm 3 . Given an L or LL composition, the bulk porosity of Eros is in the range of 25-29% and the macroporosity is 14-18%. This is consistent with a fractured, but coherent asteroid. Additional information is contained in the original extended abstract.

Britt, D. T.↗
DOE OSTI2021

An asymptotically compatible treatment of traction loading in linearly elastic peridynamic fracture

Meshfree discretizations of state-based peridynamic models are attractive due to their ability to naturally describe fracture of general materials. However, two factors conspire to prevent meshfree discretizations of state-based peridynamics from converging to corresponding local solutions as resolution is increased: quadrature error prevents an accurate prediction of bulk mechanics, and the lack of an explicit boundary representation presents challenges when applying traction loads. Herein, we develop a reformulation of the linear peridynamic solid (LPS) model to address these shortcomings, using improved meshfree quadrature, a reformulation of the nonlocal dilatation, and a consistent handling of the nonlocal traction condition to construct a model with rigorous accuracy guarantees. In particular, these improvements are designed to enforce discrete consistency in the presence of evolving fractures, whose a priori unknown location render consistent treatment difficult. In the absence of fracture, when a corresponding classical continuum mechanics model exists, our improvements provide asymptotically compatible convergence to corresponding local solutions, eliminating surface effects and issues with traction loading which have historically plagued peridynamic discretizations. When fracture occurs, our formulation automatically provides a sharp representation of the fracture surface by breaking bonds, avoiding the loss of mass. We provide rigorous error analysis and demonstrate convergence for a number of benchmarks, including manufactured solutions, free-surface, nonhomogeneous traction loading, and composite material problems. Finally, we validate simulations of brittle fracture against a recent experiment of dynamic crack branching in soda-lime glass, providing evidence that the scheme yields accurate predictions for practical engineering problems.

42 ENGINEERING↗
DOE OSTI2022

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗
DOE OSTI2023

A robust interface finite element formulation for modeling brittle material failure problems

Abstract Failure of many brittle materials and structures can be modeled using interface‐oriented finite elements combined with intrinsic cohesive zone models. The discontinuous Galerkin (DG) finite element method provides an innovative framework for modeling brittle crack propagation with zero‐thickness interface elements, which can accommodate extrinsic cohesive laws to avoid the artificial compliance required in intrinsic cohesive models. However, robust formulations and implementations of DG methods are critical in alleviating the well‐known convergence issues for both crack nucleation and propagation with reduced instability. This paper presents a robust interface element formulation by modifying the incomplete interior penalty Galerkin (IIPG) method, which successfully avoids the initial element interface penetration across elements that occurs prior to crack nucleation, and thereby greatly reduces the instability issue as cracks open. We further verified and validated our implementation by using a bar tension test and a beam fracturing benchmark. The robustness of our proposed interface element method was demonstrated by a micromechanics fiber/matrix debonding problem with 64 fibers embedded in a bulk matrix.

Liu, Ruijie↗