Effect of chemical composition on the fracture toughness of bulk metallic glasses
Not Available
SEARCH · Search NASA
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.
Not Available
Not Available
Human activities involving subsurface reservoirs—resource extraction, carbon and nuclear waste storage—alter thermal, mechanical, and chemical steady-state conditions in these systems. Because these systems exist at lithostatic pressures, even minor chemical changes can cause chemically assisted deformation. Therefore, understanding how chemical effects control geomechanical properties is critical to optimizing engineering activities. The grand challenge in predicting the effect of chemical processes on mechanical properties lays in the fact that these phenomena take place at molecular scales, while they manifest all the way to reservoir scales. To address this fundamental challenge, we investigated chemical effects on deformation in model and real systems spanning molecular- to centimeter scales. We used theory, experiment, molecular dynamics simulation, and statistical analysis to (1) identify the effect of simple reactions, such as hydrolysis, on molecular structures in interfacial regions of stressed geomaterials; (2) quantify chemical effects on the bulk mechanical properties, fracture and displacement for granular rocks and single crystals; (3) develop initial understanding of universal scaling for individual displacement events in layered geomaterials; and (4) develop analytic approximations for the single-chain mechanical response utilizing asymptotically correct statistical thermodynamic theory. Taken together, these findings advance the challenging field of chemo-mechanics.
In many high energy physics neutrino beamlines, targets are made up of isotropic POCO graphite grade for production of neutrinos for high energy particle physics research and are bombarded with highly energetic pulsed proton beam. The pulsed proton beam with a small beam size creates thermal stress waves as well as radiation damage such as displacement damage, void formation, swelling, and gas formation amongst others. As a result, after a few years of operation one of the longest operated targets appears to have undergone bulk swelling and fractured along the beam centerline. A complex interaction of dynamic loading due to beam, material degradation due to irradiation are thought to have caused such a failure. Here, we present a numerical simulation to capture the combined effects of swelling and the dynamic effect of thermal stress waves due to beam heating, to explain the failure of target. An empirical formula has been developed to take into account swelling as a function of temperature and proton fluence and has been implemented in a commercial finite element code. Extensive X-ray diffraction (XRD) were performed on graphite samples from the failed target to understand the lattice parameter changes due to irradiation and build a valid empirical model. Simulation results show that swelling plays a major role in elevating the stress state in the material exceeding the failure strength while the pulsed beam heating also introduces fatigue loading that would explain brittle fracture emanating from inside of the material and propagating outward. Such empirical formulations will help in identifying the threshold values of critical parameters to extend the service life of future multi-megawatt neutrino targets and targets in other accelerator environments.
Rock fracture surfaces in the crust are essential habitat for microorganisms. Fracture‐groundwater interfaces provide physical substrates for biofilm growth and are sources of carbon, nutrients, and electron donors and acceptors. To better understand geochemical processes impacting fracture surfaces and the subsurface microbiome, we identified fractures in archived rock cores from the Soudan formation, which is known to host saline groundwaters and isolated microbial communities dependent on rock‐water interactions. Cores with open fractures were thin sectioned and studied via electron microprobe and synchrotron X‐ray fluorescence microprobe. Most fracture surfaces had mineralogy distinct from that of the bulk rock. Chlorite minerals were abundant on fracture surfaces and had elemental compositions suggesting deposition during late‐stage hydrothermal alteration. Fracture‐lining chlorites likely limit access to iron oxide and sulfide minerals that are active in subsurface biogeochemical cycles. Calcium‐rich rinds were also observed along fracture edges. These rinds were too thin and poorly ordered to be identified via light microscopy or X‐ray diffraction; however, Ca K‐edge micro‐X‐ray absorption near‐edge structure spectroscopy identified them as carbonates, minerals not observed in the bulk rock. Thermodynamic modeling shows that carbonate precipitation is largely unfavorable in Soudan groundwaters, indicating that fracture edge conditions differed from those in modern water samples. Because of the low carbon concentrations in Soudan groundwaters, carbonate rinds likely play an important role in subsurface carbon cycling and may mark fracture surfaces that once hosted biofilms. Overall, this study suggests that fracture alteration can both play an active role in and suppress rock‐water interactions essential to subsurface life.
Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).
The presence of fractures in caprocks can pose increased risks in subsurface energy systems and processes like CO 2 sequestration by introducing high-permeability leakage paths. Fracture apertures and permeability can be altered through mineral dissolution and precipitation reactions, but the reactive evolution of fractures is not well understood. In fractures, minerals that are otherwise inaccessible to reactive fluids can become exposed, resulting in mineral reactions unpredicted by bulk formation data. This work seeks to understand the relationship between mineralogy and fracture formation to enhance our understanding of reactive fracture evolution and CO 2 leakage potential. Here, the mineral compositions of mechanically induced fracture surfaces in samples of the Mancos and Marcellus shales have been quantified and compared to those of the near-fracture matrices using imaging and bulk X-ray diffraction (XRD) data. In the Mancos shale, the concentrations of clay minerals are enhanced along fracture surfaces with respect to the bulk, and the fracture is most likely to form at kaolinite–kaolinite interfaces. Further, evaluation of the mineralogical spatial variability through cross-correlation analysis of the surrounding matrix in images of samples cut perpendicular to the fracture shows that clay is 16.7 times more likely to be present than carbonate minerals near the fracture surface. The high correlation persists roughly 200 μm into the surrounding matrix for the Mancos sample and implies that the fracture formed within a defined clay-rich lithofacies.
Abstract Models exploring electrochemistry-mechanics coupling in liquid electrolyte lithium-ion battery anodes have traditionally incorporated stress impact on thermodynamics, bulk diffusive transport, and fracture, while stress-kinetics coupling is more explored in the context of all solid-state batteries. Here, we showcase the existence of strong link between active particle surface pressure and reaction kinetics affecting performance even in liquid electrolyte systems. Traction-free and immobile particle surface mechanical boundary conditions are used to delineate the varying pressure magnitudes in graphite host during cycling. Both tensile and compressive stresses are generated in traction-free case, while a fixed surface subjects the entire particle to a compression state. Pressure magnitudes are nearly two to three orders of magnitude higher for the latter resulting in significant depression of open circuit potential and improvement of exchange current densities compared to stress-free state. The results demonstrate the need for incorporating stress-kinetics linkage in models and provide a rationale for putting battery electrodes under compression to improve kinetics.
Fracture and crack growth is one of the main degradation mechanisms in solid oxide cells (SOCs). However, the modeling of crack growth in SOCs is challenging due to their complex microstructures and possible plasticity development within the Ni particles in Ni-based SOC electrodes. In this study, a phase-field fracture model is developed, which incorporates the SOC microstructures and phase-dependent material properties, including yield strength, fracture toughness in the bulk and at the interphase boundaries. The model is employed to study crack initiation and growth under thermal and redox cycling on the hydrogen electrode side of SOCs. The simulation results demonstrate that under thermal cycling, work-zone cracking dominates in electrolyte-supported SOCs with cracks initiated at the triple-phase boundaries, while only minor mechanical degradation occurs in hydrogen-electrode-supported SOCs after hundreds of thermal cycles. Under redox cycling, through-cracking of yttria-stabilized zirconia (YSZ) in the hydrogen electrode and electrolyte layers dominates. Finally, the simulation results suggest several crack-mitigation strategies, including decreasing the porosity in the hydrogen electrode support layer and synchronizing thermal strain to balance oxidation strain.
This work investigates the ability of graph neural networks (GNNs) to homogenize 2D fiber composite microstructures. We use different inhomogeneity and anisotropy indices to motivate and show that the Volume Elements (VEs) used in ML methods should ideally be far from their Representative Volume Element (RVE) size limit and, consequently, are notably anisotropic. Hence, training only the isotropic limit properties may not be acceptable. Another aspect is the need to normalize elastic stiffness values for ML, especially when high elastic contrast ratios are encountered between composite phases or in the material set. We introduce a normalization technique based on the mean-field method (MFM) to handle such high contrast ratios and train for the entire stiffness tensor. We show that the proposed GNN approaches exhibit high accuracy and efficiency compared to traditional methods and convolutional neural networks, utilizing unstructured graphs constructed from microstructure topology. Our model successfully predicts the stiffness tensor, peak strength under bulk damage, and brittle fracture initiation strength across diverse microstructure configurations while maintaining high accuracy even for extreme material contrasts and volume fractions. We also present a method to improve prediction accuracy for small dataset sizes using Voronoi partitioning.
A “screening analysis” is implemented at Pahute Mesa (PM), Nevada National Security Site (NNSS), to determine subsets of the 43 radionuclides listed in the radionuclide inventory of Finnegan et al. (2016) that are of relevance or potential relevance to the hydrologic source term (HST) for assessment of radionuclide transport in groundwater. Consideration is added to how levels of contamination are defined and whether any other radionuclides not in the inventory are relevant to the HST. A model is developed to estimate a range of possible source concentrations in groundwater that account for uncertainty in partitioning into melt glass and sorption into surrounding rock of the exchange zone for the 82 PM underground nuclear tests detonated in vertical shafts. The model also accounts for the varied hydrogeochemical settings. Transport parameters needed for the screening model calculations are developed from databases for hydrogeologic units, chemistry, mineralogy, fracture spacing, fracture aperture, fracture openness, matrix porosity, bulk density, saturation, and alteration. Of key importance is consideration of diagenetic zonation of the mineralogy related to hydrothermal alteration. The screening analysis compares model results with available groundwater radiochemistry data. Consistency between screening model results and available data helps validate the model for application to all 43 radionuclides, most of which have no measurements or only non-detect measurements of concentration in groundwater. The screening analysis determines that ten radionuclides – tritium, Sr-90, I-129, Cs-137, U-232, U-233, U-234, U-238, and the total of Pu-239 and Pu-240 are relevant to the HST. Determination of relevance is primarily based on data and/or model results indicating source concentration exceeding a maximum contaminant level (MCL) out to 100 years since the time of detonation. U-238 is relevant to assessment of contamination relative to MCLs for uranium and gross alpha particle activity. Five other radionuclides – C-14, Cl-36, Tc-99, Np-237, and Pu-238 – are determined potentially relevant to the HST based on a 0.1 MCL threshold. The screening analysis recommends additional attention to the natural daughter radionuclides of the uranium series (U-238) and thorium series (Th-232) decay chains, particularly Ra-226 and Ra-228, that were not included the inventory.
Internationally, it has been agreed that geologic repositories for spent fuel and radioactive waste are considered the internationally agreed upon solution for intermediate and long-term disposal. In countries where traditional nuclear waste repository host rocks (e.g., clay, salt, granite) are not available, other low permeability lithologies must be studied. Here, chalk is considered to determine its viability for disposal. Despite chalk's low bulk permeability, it may contain fracture networks that can facilitate radionuclide transport. In arid areas, groundwater salinity may change seasonally due to the mixing between brackish groundwater and fresh meteoric water. Such salinity changes may impact the radionuclides' mobility. In this study, radioactive U(VI) and radionuclide simulant tracers (Sr, Ce and Re) were injected into a naturally fractured chalk core. The mobility of tracers was investigated under abrupt salinity variations. Two solutions were used: a low ionic strength (IS) artificial rainwater (ARW; IS ~0.002) and a high IS artificial groundwater (AGW; IS ~0.2). During the experiments, the tracers were added to ARW, then the carrier was changed to AGW, and vice versa. Ce was mobile only in colloidal form, while Re was transported as a conservative tracer. Both Re and Ce demonstrated no change in mobility due to salinity changes. In contrast, U and Sr showed increased mobility when AGW was introduced and decreased mobility when ARW was introduced into the core. These experimental results, supported by reactive transport modeling, suggest that saline groundwater solutions promote U and Sr release via ion-exchange and enhance their migration in fractured chalk. Here the study emphasizes the impact of salinity variations near spent fuel repositories and their possible impact on radionuclide mobility.
Serpentine minerals have received a lot of attention because of their unique crystal structures, their wide occurrence in orogenic belts and their potential role in contributing seismic anisotropy in subducting slabs. Several studies have investigated crystal preferred orientation (CPO) in high temperature antigorite serpentinites from Japan, the Alps, Spain, Cuba and Tibet, documenting significant crystal alignment. However, only a limited number of lower grade serpentines have been explored to date. Mainly because of submicroscopic microstructural heterogeneities CPO cannot be measured with conventional methods such as optical microscopy and EBSD. In this study 15 serpentinites from different tectonic settings in California, the Central Alps and Northern Spain have been investigated, mainly with high energy synchrotron X-ray diffraction, to quantify bulk crystal alignment. We find that CPO is strong on sheared surfaces of fractured blocks and secondary veins but the bulk of most serpentinite samples, except high-grade recrystallized antigorite serpentinite, show only weak crystal alignment. Correspondingly calculated seismic anisotropy based on CPO is not very significant. This is supported by very heterogeneous microstructures as documented with SEM and TEM analyses.
Stress relaxation cracking (SRC) is considered one of the major failure mechanisms for 347H stainless steel welds at elevated service temperatures or during post weld heat treatment (PWHT), especially within the heat-affected zone (HAZ). This work focuses on the characterization of SRC susceptibility within 347H physically simulated arc welded HAZ at elevated temperatures. A four-step SRC thermomechanical test in combination with finite element modeling (FEM) of the welding and testing processes is developed to establish a susceptibility map for HAZ. The test first runs a thermal cycle with three different peak temperatures (1335, 1275, and 1150 °C) to duplicate representative HAZ subzone microstructures, followed by time-to-failure examination under a variety of pre-stress (260–600 MPa) and pre-strain conditions (0.03–0.19) as a function of reheat temperatures between 750 and 1050 °C. With the aid of FEM, SRC susceptibility maps are generated to identify the threshold stress, plastic strain, and creep strain as a function of test temperature. It was found out that HAZ subzone with a lower peak temperature (1150 °C) appears to be slightly less susceptible to SRC than the other two subzones that experienced higher peak temperatures. Generally, time-to-fracture reduces with increasing initially applied stress and strain for all test temperatures. The pre-stress thresholds decrease from about 500 to 330 MPa as the testing temperature increases from 800 to 1050 °C, while the corresponding initial plastic strain thresholds reduces from 0.15 to 0.06. The SRC susceptibility was also evaluated through the Larson–Miller Parameter (LMP) analysis as a function of plastic strain, initial stress and starting stress upon reaching the testing temperature, respectively. The 1050 °C test with a high pre-applied strain (0.1) exhibits an extremely short time to failure (t = 3 s) that lies outside the general trend in LMP analysis. Additionally, it was identified that a plastic strain above 0.07 is identified to significantly reduce the bulk creep strain tolerance to fracture and therefore increases SRC susceptibility. Hardness measurement and fractography analysis indicated that the strain aging of niobium carbonitrides and other potential phases in conjunction with intergranular precipitates contributes to an increase in microhardness and increased intergranular cracking susceptibility.
Active nematic systems consist of rod-like internally driven subunits that interact with one another to form large-scale coherent flows. They are important examples of far-from-equilibrium fluids, which exhibit a wealth of nonlinear behavior. This includes active turbulence, in which topological defects in the nematic order braid around one another in a chaotic fashion. One of the most studied examples of active nematics consists of a dense two-dimensional layer of microtubules, crosslinked by kinesin molecular motors that inject extensile deformations into the fluid. Though numerous theoretical studies have modeled microtubule-based active nematics, no consensus has emerged on how to fully and quantitatively capture the features of the experimental system. Here, to better understand the theoretical foundations for modeling this system, we propose a fundamental principle we call the nematic locking principle—individual microtubules cannot rotate without all neighboring microtubules also rotating. Physically, this is justified by the high density of the microtubules, their elongated nature, and their corresponding steric interactions. We assert that the nematic locking principle holds throughout the majority of the material, but breaks down in the neighborhood of topological defects and other regions of low density. We derive the most general nematic transport equation consistent with this principle and also derive the most general term that violates it, introducing fracturing into the material. We then examine the standard Beris–Edwards approach, commonly used to model this system, and show that it violates the nematic locking principle throughout the majority of the material due to fracturing. We then propose a modification to the Beris–Edwards model that enforces nematic locking nearly everywhere. This modification shuts off fracturing except in regions where the order parameter (a proxy for density) is reduced. In these regions fracturing is turned on. The resulting simulations in turn show strong nematic locking throughout the bulk of the material, with narrow bands of fracturing, consistent with experimental observation. One additional advantage of enforcing nematic locking is that nontrivial stationary state solutions, common in Beris–Edwards simulations but not seen in experiments, are eliminated.
Herein, the mechanical behaviors of Li 10 GeP 2 S 12 (LGPS) solid electrolytes during electrochemical cycling using operando X-ray tomography are investigated. It is demonstrated that the bulk mechanical decomposition of LGPS when cycled against lithium is a direct result of electrochemical reduction of the solid electrolyte at the LGPS/Li 0 interface. The reductive decomposition of LGPS during lithium plating results in the formation of low-density domains at the electrode/electrolyte interface, which impose sufficient mechanical stress on the underlying LGPS to crack the SE pellet. The critical stress developed prior to pellet fracture is significantly lower than the bulk shear modulus of LGPS, suggesting that the electrochemical instability of LGPS dramatically worsens the mechanical stability of the material near the LGPS/Li 0 interface. It is also shown that the application of a highly concentrated liquid electrolyte to the LGPS surface suppresses the reductive decomposition of LGPS, improving both the electrochemical performance and mechanical stability of the bulk LGPS solid electrolyte.
Fractured rocks play a crucial role in myriad natural and engineered systems, including Earth's critical zone, oil/gas/geothermal reservoirs, and geological CO 2 /H 2 /waste storage systems. While complex electrical conductivity is extensively used to estimate the pore and grain sizes of conventional porous rocks and soils, it is rarely used to predict the aperture size of fractured rocks and this remains poorly understood. Here, integrating theory, simulations, and experiments, we show that under external fields, fractured rocks follow the fracture-and-matrix coupling to make the bulk complex conductivity non-linear with respect to water conductivity. We find that the relaxation time and quadrature conductivity for porous media do not apply to fractured rocks, but, instead, reasonably accurate predictions of aperture size can be made based on the true formation factor. This study unravels the fundamental mechanism governing conduction and polarization of fractured rocks and paves the way for the non-invasive investigation of global fractured rocks.
Moiré superlattices formed by twisted stacking in van der Waals materials have emerged as a new platform for exploring the physics of strongly correlated materials and other emergent phenomena. However, there remains a lack of research on the mechanical properties of twisted-layer van der Waals materials, owing to a lack of suitable strategies for making three-dimensional bulk materials. Here we report the successful synthesis of a polycrystalline boron nitride bulk ceramic with high room-temperature deformability and strength. This ceramic, synthesized from an onion-like boron nitride nanoprecursor with conventional spark plasma sintering and hot-pressing sintering, consists of interlocked laminated nanoplates in which parallel laminae are stacked with varying twist angles. The compressive strain of this bulk ceramic can reach 14% before fracture, about one order of magnitude higher compared with traditional ceramics (less than 1% in general), whereas the compressive strength is about six times that of ordinary hexagonal boron nitride layered ceramics. The exceptional mechanical properties are due to a combination of the elevated intrinsic deformability of the twisted layering in the nanoplates and the three-dimensional interlocked architecture that restricts deformation from propagating across individual nanoplates. The advent of this twisted-layer boron nitride bulk ceramic opens a gate to the fabrication of highly deformable bulk ceramics.