SEARCH · Engineering Papers

Results for “Strain calculation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Exceptional fracture toughness of CrCoNi-based medium- and high-entropy alloys at 20 kelvin

Medium- and high-entropy alloys based on the CrCoNi-system have been shown to display outstanding strength, tensile ductility and fracture toughness (damage-tolerance properties), especially at cryogenic temperatures. Here we examine the J Ic and (back-calculated) K JIc fracture toughness values of the face-centered cubic, equiatomic CrCoNi and CrMnFeCoNi alloys at 20 K. At flow stress values of ~1.5 GPa, crack-initiation K JIc toughnesses were found to be exceptionally high, respectively 235 and 415 MPa√m for CrMnFeCoNi and CrCoNi, with the latter displaying a crack-growth toughness K ss exceeding 540 MPa√m after 2.25 mm of stable cracking, which to our knowledge is the highest such value ever reported. Characterization of the crack-tip regions in CrCoNi by scanning electron and transmission electron microscopy reveal deformation structures at 20 K that are quite distinct from those at higher temperatures and involve heterogeneous nucleation, but restricted growth, of stacking faults and fine nanotwins, together with transformation to the hexagonal closed-packed phase. The coherent interfaces of these features can promote both the arrest and transmission of dislocations to generate respectively strength and ductility which strongly contributes to sustained strain hardening. Indeed, we believe that these nominally single-phase, concentrated solid-solution alloys develop their fracture resistance through a progressive synergy of deformation mechanisms, including dislocation glide, stacking-fault formation, nano-twinning and eventually in situ phase transformation, all of which serve to extend continuous strain hardening which simultaneously elevates strength and ductility (by delaying plastic instability), leading to truly exceptional resistance to fracture.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Bayesian Inference of State-Level COVID-19 Basic Reproduction Numbers across the United States

Although many persons in the United States have acquired immunity to COVID-19, either through vaccination or infection with SARS-CoV-2, COVID-19 will pose an ongoing threat to non-immune persons so long as disease transmission continues. We can estimate when sustained disease transmission will end in a population by calculating the population-specific basic reproduction number ℛ0, the expected number of secondary cases generated by an infected person in the absence of any interventions. The value of ℛ0 relates to a herd immunity threshold (HIT), which is given by 1−1/ℛ0. When the immune fraction of a population exceeds this threshold, sustained disease transmission becomes exponentially unlikely (barring mutations allowing SARS-CoV-2 to escape immunity). Here, we report state-level ℛ0 estimates obtained using Bayesian inference. Maximum a posteriori estimates range from 7.1 for New Jersey to 2.3 for Wyoming, indicating that disease transmission varies considerably across states and that reaching herd immunity will be more difficult in some states than others. ℛ0 estimates were obtained from compartmental models via the next-generation matrix approach after each model was parameterized using regional daily confirmed case reports of COVID-19 from 21 January 2020 to 21 June 2020. Our ℛ0 estimates characterize the infectiousness of ancestral strains, but they can be used to determine HITs for a distinct, currently dominant circulating strain, such as SARS-CoV-2 variant Delta (lineage B.1.617.2), if the relative infectiousness of the strain can be ascertained. On the basis of Delta-adjusted HITs, vaccination data, and seroprevalence survey data, we found that no state had achieved herd immunity as of 20 September 2021.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Constructing Regulatory Networks to Compare Axenic and Interspecies Microbial Gene Transcription

In this preliminary study, we constructed gene regulatory networks (GRNs) from transcriptional expression data of axenic and interspecies microbial cultures with the goal of predicting how cocultivation affected greenhouse gas respiration by these species. The specific strains of Methylotuvimicrobium alkaliphilum 20Z, a methylotroph, and Cyanobacterium stanieri HL-69, a phototroph, were chosen for their viability in industrial bioprocessing. We ranked directed interactions between gene pairs based on the ability of the input gene to predict the expression of a target gene relative to their transcriptomes. While we were able to identify topological differences between conditions, our initial findings require validation through experimental analysis and further modeling. We aimed to develop a systematic thresholding approach to optimize the accuracy of our networks. We filtered out trial networks separately from top gene interactions of the scored rankings. Parameters of unfiltered and filtered networks were used to test and develop thresholding approaches. Knee point detection of edge weight distributions was explored as an approach for separating significant interactions from insignificant interactions in unfiltered networks. While knee detection failed to produce analogous networks for broad cross-condition comparisons, the results informed us about the proportions of significant edges present in unfiltered networks. We also calculated the average mean degree for nodes in a selection of trial networks to find a thresholding value characteristic to all groups. While we did not reach a definitive conclusion, we gained insight into the coregulatory structures of our groups and made critical evaluations of systematic methods for filtering networks. We recommend an iterative process for the inference of GRNs, where the most significant results from preliminary explorations are used to improve the efficiency with which regulatory motifs are chosen for experimental characterization. Experimental results can then inform the framework of adjusted models to improve broad interpretations of GRNs.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Rapid identification of enteric bacteria from whole genome sequences using average nucleotide identity metrics

Identification of enteric bacteria species by whole genome sequence (WGS) analysis requires a rapid and an easily standardized approach. We leveraged the principles of average nucleotide identity using MUMmer (ANIm) software, which calculates the percent bases aligned between two bacterial genomes and their corresponding ANI values, to set threshold values for determining species consistent with the conventional identification methods of known species. The performance of species identification was evaluated using two datasets: the Reference Genome Dataset v2 (RGDv2), consisting of 43 enteric genome assemblies representing 32 species, and the Test Genome Dataset (TGDv1), comprising 454 genome assemblies which is designed to represent all species needed to query for identification, as well as rare and closely related species. The RGDv2 contains six Campylobacter spp., three Escherichia/Shigella spp., one Grimontia hollisae, six Listeria spp., one Photobacterium damselae, two Salmonella spp., and thirteen Vibrio spp., while the TGDv1 contains 454 enteric bacterial genomes representing 42 different species. The analysis showed that, when a standard minimum of 70% genome bases alignment existed, the ANI threshold values determined for these species were ≥95 for Escherichia/Shigella and Vibrio species, ≥93% for Salmonella species, and ≥92% for Campylobacter and Listeria species. Using these metrics, the RGDv2 accurately classified all validation strains in TGDv1 at the species level, which is consistent with the classification based on previous gold standard methods.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Slip localization behavior at triple junctions in nickel-base superalloys

Incipient slip localization in the vicinity of hundreds of grain boundary triple junctions (TJs) in a lightly deformed nickel-base superalloy IN718 is studied using a combination of three-dimensional (3D) crystal plasticity finite element (CPFE) modeling, high resolution digital image correlation (HR-DIC) and 3D electron back-scatter diffraction tomography (3D EBSD). A 3D reconstruction method enables identification of thousands of TJs and correspondence of any observed slip bands with their originating TJ lines below the specimen surface. Here, we present a large-scale CPFE model of the experimental 3D microstructure composed of high-fidelity representation of the TJ lines and the boundaries and interiors of the parent grains and use it to calculate the local micromechanical response and slip activity of all TJs at the onset of macroscopic yielding. Statistical analysis of the calculated quantities reveal TJs develop larger stress concentration and grain-average re-orientation than grain interiors and grain boundaries, however no substantial differences in cumulative slip were found among these microstructural regions. We find that TJs with observed slip bands generate lower grain-average re-orientation, fewer active slip systems, and more localized slip on a single system than those without. The distinctions in the reorientation and slip activity are stronger in TJs that experience more intense slip.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Nanometer-Scale Fullerene-Type Conjugated Covalent Cages Based on Triazine: Design, Doping with Li + , and H 2 /CO 2 Adsorption

Spherical and hollow molecular cages based on planar triazine (C 3 N 3 ) hubs and aromatic phenylene connectors have been developed. The cages exhibit topologies akin to C 20 , C 60 , and C 70 fullerenes with diameters that range from 2.3 to 4.9 nm. Apertures into the cage interiors are tuned by varying the aromatic connectors situated between the C 3 N 3 -units. The stabilities of the C 3 N 3 cages increase with size owing to reduced bending strain of planar nodes and connectors that make up the spherical aromatic networks. Doping of the cages with Li + reveals the capacity of the cages for significant adsorption of gaseous H 2 and CO 2 . The design of graphene-like spherical cages is also discussed.

DFT calculations↗

DOE OSTI2020

Three-dimensional modeling and simulations of single-crystal and bi-crystal titanium for high-strain-rate loading conditions

High purity single crystal titanium (Ti) under shock wave loading is modeled under both one-dimensional and three-dimensional cylindrical conditions. Cylinder sizes of 10 μm and 20 μm radius are both considered in order to assess influence of boundary conditions. A thermodynamically consistent single crystal model for application to shock conditions is presented. The model accounts for the coupled non-linear elastic, dislocation slip, deformation twinning, and structural phase transformation response of the titanium material. Plate impact experiment results using a copper flyer are used to compare against the simulations for crystals oriented in [0001] and [10$\overline1$1] crystallographic directions. The one-dimensional and threedimensional simulations of the two differently oriented single crystals indicate differences between the one-dimensional and three-dimensional representation, especially for the [10$\overline1$1] oriented single crystal. This orientation breaks the relative orientation symmetry between the crystal and cylinder which otherwise exists for the [0001] oriented single crystal. A significant amount of heterogeneity in the field response of the[10$\overline1$1] oriented simulation was demonstrated due to the highly coupled nature of the deformation. A bi-crystal model composed of both the [0001] and [10$\overline1$1] orientations with the boundary between the two along the axis of the cylinder is also considered for a cylinder model size of 10 μm. The results indicate a strong interaction between the two grains that affects the ω phase volume fraction achieved relative to the single crystal calculations.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Experimental and Computational Studies of Stress Corrosion Cracking of Alloys 308/309 and 82/182 Weldments in Corrosive and Radiation Environment

The goal of the project was to determine factors that influence SCC and IASCC in weldments found in LWR nuclear power plants. We focused on a SA508-304L SS weldment fabricated by EPRI using gas tungsten arc welding and used an aggressive BWR normal water chemistry (NWC) immersion environment. This weldment used 309L butter and 308L groove filler material. The microstructure of the 309L butter was non-uniform, exhibiting a 20–30μm thick martensitic layer closest to the SA508 interface, a 1–4 mm thick single γ austenite phase dilution zone, and a γ–δ duplex region extending to the 308L groove filler. The 308L groove filler had an entirely γ–δ duplex microstructure. Two approximately 1-inch thick 304L and SA508 plates approximately 12 by 6 square inches were joined using standard nuclear grade welding techniques. This included a post weld heat treatment of the 309L butter after application, a 0.32 cm fit-up root opening, and 57 bead lines of 308L groove filler applied in 18 layers. A 60 degree weld bevel angle was used and the 309L butter was 1.5 cm thick. Displacement cascade damage was induced using proton irradiation at the Michigan Ion Beam Laboratory. The incident proton energy was 2 MeV, the sample temperature was 360 ºC, and the calculated dpa value at 10 μm (60% of the Bragg peak depth of ~18 μm) was 5 dpa using the quick Kinchin-Pease model. Proton irradiation to this damage level required approximately 125 hours of beam time. Two types of samples were irradiated, tensile specimens and TEM bars. Samples were selected from all regions of the weldment, including the SA508-309L butter interface, the 309L-308L interface, and 308L-304L interface. The stainless steel alloys (304L, 308L, and 309L) within the heat affected zone are characterized by a duplex skeletal morphology of δ-ferrite and γ-austenite resulting from the recrystallization associated with weld fabrication. Approximately 7 to 8 mm of length along the specimen was irradiated. The gauge volume surfaces were mechanically polished and then electro-polished to remove mechanical damage from the mechanical polishing step prior to irradiation. Immersion tests were performed in a recirculating autoclave under BWR NWC conditions (2000 ppb wt. dissolved oxygen, neutral pH, 288 ºC, 10 MPa, and inlet water conductivity <100 nS/cm) to accelerate corrosion. Constant strain rate tests were performed either to failure or to approximately 5% strain. Strain rates of 10 -7 to 10 -6 mm/mm/s were used and typical immersion testing required four to six weeks to achieve failure or strains near 5%. Analysis primarily used advanced electron microscopy techniques of FIB lift out specimens.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗

DOE OSTI2025

Achieving High-Yield Conversion of Janus Transition Metal Dichalcogenides on Diverse Substrates

Janus transition metal dichalcogenides (TMDCs) with intrinsic broken mirror symmetry and vertical dipole moment provide an additional degree of freedom to manipulate material symmetry down to atomic-layer thickness. However, despite advances in synthesis strategies, fundamental understanding of this atomic substitution process remains limited, which has impeded their implementation in advanced devices. Here, by using a room-temperature atomic-layer substitution (RT-ALS) strategy, we systematically investigate the critical factors facilitating the high-yield conversion of Janus TMDCs on diverse substrates. Combining Raman spectroscopy probes, X-ray photoelectron spectroscopy (XPS) measurements, and density functional theory (DFT) calculations, we demonstrate that substrates with enhanced electron doping or larger surface polarity substantially benefit the conversion of Janus TMDCs reaching a near-unity yield. Intriguingly, the strong affinity between Janus TMDCs and substrates (e.g., Au) brings about abnormal Raman spectroscopic phenomena. These findings highlight the significance of substrates in achieving the reliable synthesis of Janus two-dimensional materials with improved homogeneity on various substrates. In addition, this takes us one step closer to utilizing Janus TMDCs as a versatile platform in next-generation optoelectronic devices, sensors, and quantum technologies.

Doping↗

DOE OSTI2020

Transport of coexisting Ni-Cu sulfide liquid and silicate melt in partially molten peridotite

Transport of coexisting sulfide and silicate melts in partially molten peridotite contributes to the redistribution of chalcophile elements within the upper mantle as well as the genesis of magmatic Ni-Cu sulfide deposits, but has not been investigated systematically. Using laboratory experiments, theoretical calculations, and X-ray synchrotron microtomography, this study documents the topology and considers controls on the extraction of two immiscible liquids during partial melting of mantle peridotite. Under hydrostatic conditions, the measured dihedral angle at silicate melt-mineral-mineral junctions is 13.7–21.3°. Silicate melt is distributed along grain edges forming incompletely interconnected melt channels that disconnected by some dead ends at melt fractions ~7–9 vol%. Application of theoretically predicted permeability (k ~10 -14 -10 -16 m 2 ) permits estimation of the extraction velocity of silicate melt of 0.7-11.1 μm/day within a single interconnected melt channel. In the absence of silicate melt, isolated sulfide droplets (3.77 vol%) show a sulfide-olivine-olivine dihedral angle of 91.5–101.3°. However, in the presence of silicate melt, sulfide droplets (average size ~2.53±2.14μm, 1σ) are partially surrounded by silicate melt and stranded in triple junctions or melt pockets due to the limitation of the smallest dimension (0.3μm) of melt channels. Thus, the extraction of sulfide liquid is highly restricted by these dead ends and the smallest dimension of melt channels during porous flow of silicate melt. In contrast, during large-strain shear deformation (strain ~1.6–2.5), initially stranded sulfide droplets were elongated and extracted with silicate melt into liquid-rich sheets with a length of several hundred microns, constantly oriented at 14.3±4.5°to the shear plane and antithetic to the shear direction. The angle is lower than that (18-30°) of those sheets containing sulfide liquid only, indicating that silicate melt dominates liquid-rich sheets. Driven by stress, silicate melt-dominated liquid-rich sheets open the appropriately oriented grain boundaries between silicate minerals, thus providing an efficient pathway for the extraction of sulfide liquid during deformation. When such sheet-like channels remain open, sulfide droplets (> millimeter-scale) can be potentially mobile through high strain domains of the upper mantle, contributing to the addition of chalcophile elements and the fertilization of the lithospheric mantle.

high P-T experiment↗

DOE OSTI2020

Single-Layer Dititanium Oxide Ti 2 O MOene: Multifunctional Promises for Electride, Anode Materials, and Superconductor

Using the first-principles calculations, we report the existence of the single-layer (SL) dititanium oxide Ti 2 O (labeled as MOene) that constructs a novel family of MXene based on transition-metal oxides. This MOene material strongly contrasts the conventional ones consisting of transition-metal carbides and/or nitrides. SL Ti 2 O has high thermal and dynamical stabilities because of the strong Ti–O ionic bonding interactions. Moreover, this material is an intrinsic electride and exhibits extremely low diffusion barriers of ~12.0 and 6.3 meV for Li and Na diffusion, respectively. When applied as anode materials in lithium-ion batteries and sodium-ion batteries, it possesses a high energy storage capacity (960.23 mAhg –1 ), surpassing the traditional MXenes-based anodes. The superb electrochemical performance stems from the existing anionic electron on Ti 2 O surface. Astonishingly, SL Ti 2 O is also determined to be a superconductor with a superconducting transition temperature (T c ) of ~9.8 K, which originates from the soft-mode of the first acoustic phonon branch and enhanced electron–phonon coupling in the low-frequency region. Furthermore, this soft-mode behaves much softer upon applying a compressive strain of 2%, leading to a higher T c of 11.9 K. Furthermore, our finding broadens the family of MXenes and could facilitate more experimental efforts toward future nanodevices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Experimental evaluation of shear modulus scaling of dynamic strength at extreme pressures

Recent progress in the development of dynamic strength experimental platforms is allowing for unprecedented insight into the assumptions used to construct constitutive models operating in extreme conditions. In this work, we make a quantitative assessment of how tantalum strength scales with its shear modulus to pressures of hundreds of gigapascals through a cross-platform examination of three dynamic strength experiments. Specifically, we make use of Split–Hopkinson pressure bar and Richtmyer–Meshkov instability experiments to assess the low-pressure strain and strain rate dependence. Concurrent examination of magnetically driven ramp-release experiments up to pressures of 350GPa allows us to examine the pressure dependence. In this work, using a modern description of the shear modulus, validated against both ab initio theory and experimental measurements, we then assess how the experimentally measured pressure dependence scales with shear modulus. Furthermore, we find that the common assumption of scaling strength linearly with the shear modulus is too soft at high pressures and offer discussion as to how descriptions of slip mediated plasticity could result in an alternative scaling that is consistent with the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Characterization of Deformational Isomerization Potential and Interconversion Dynamics with Ultrafast X-ray Solution Scattering

Dimeric complexes composed of d 8 square planar metal centers and rigid bridging ligands provide model systems to understand the interplay between attractive dispersion forces and steric strain in order to assist the development of reliable methods to model metal dimer complexes more broadly. [Ir 2 (dimen) 4 ] 2+ (dimen = para-diisocyanomenthane) presents a unique case study for such phenomena, as distortions of the optimal structure of a ligand with limited conformational flexibility counteract the attractive dispersive forces from the metal and ligand to yield a complex with two ground state deformational isomers. Here, in this work, we use ultrafast X-ray solution scattering (XSS) and optical transient absorption spectroscopy (OTAS) to reveal the nature of the equilibrium distribution and the exchange rate between the deformational isomers. The two ground state isomers have spectrally distinct electronic excitations that enable the selective excitation of one isomer or the other using a femtosecond duration pulse of visible light. We then track the dynamics of the nonequilibrium depletion of the electronic ground state population—often termed the ground state hole—with ultrafast XSS and OTAS, revealing a restoration of the ground state equilibrium in 2.3 ps. This combined experimental and theoretical study provides a critical test of various density functional approximations in the description of bridged d 8 –d 8 metal complexes. The results show that density functional theory calculations can reproduce the primary experimental observations if dispersion interactions are added, and a hybrid functional, which includes exact exchange, is used.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Computational design and analysis of modular cells for large libraries of exchangeable product synthesis modules

Microbial metabolism can be harnessed to produce a large library of useful chemicals from renewable resources such as plant biomass. However, it is laborious and expensive to create microbial biocatalysts to produce each new product. To tackle this challenge, we have recently developed modular cell (ModCell) design principles that enable rapid generation of production strains by assembling a modular (chassis) cell with exchangeable production modules to achieve overproduction of target molecules. Previous computational ModCell design methods are limited to analyze small libraries of around 20 products. In this study, we developed a new computational method, named ModCell-HPC, that can design modular cells for large libraries with hundreds of products with a highly-parallel and multi-objective evolutionary algorithm and enable us to elucidate modular design properties. We demonstrated ModCell-HPC to design Escherichia coli modular cells towards a library of 161 endogenous production modules. From these simulations, we identified E. coli modular cells with few genetic manipulations that can produce dozens of molecules in a growth-coupled manner with different types of fermentable sugars. These designs revealed key genetic manipulations at the chassis and module levels to accomplish versatile modular cells, involving not only in the removal of major by-products but also modification of branch points in the central metabolism. We further found that the effect of various sugar degradation on redox metabolism results in lower compatibility between a modular cell and production modules for growth on pentoses than hexoses. To better characterize the degree of compatibility, we developed a method to calculate the minimal set cover, identifying that only three modular cells are all needed to couple with up 85 compatible production modules. By determining the unknown compatibility contribution metric, we further elucidated the design features that allow an existing modular cell to be re-purposed towards production of new molecules. Altogether, ModCell-HPC is a useful tool for understanding modularity of biological systems and guiding more efficient and generalizable design of modular cells that help reduce research and development cost in biocatalysis.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Factors Affecting the Electron Conductivity in Single Crystal Li 7 La 3 Zr 2 O 12 and Li 7 P 3 S 11

One of the serious challenges in all solid-state Li ion batteries is neutral Li intrusion into the solid-state electrolyte that can ultimately cause catastrophic failure. One possibility for this is due to n-type electron conductivity that induces the reaction Li + + e – → Li 0 at sites where the potential is less than the Li + /Li potential. This paper reports hybrid density functional theory calculations of the electronic conductivity in two prototype single crystalline solid-state electrolytes, cubic Li 7 La 3 Zr 2 O 12 (c-LLZO) and Li 7 P 3 S 11 (LPS). The formation energies of important point defects that can affect electron conductivity are determined, and we find that the mechanism of n-type electron conductivity for both solid-state electrolytes is via “small” electron polaron hopping, where the quotes signify that substantial Li ion rearrangement is associated with the polaron formation and its migration. In both electrolytes, the formation energies for the small polarons at the Fermi energy are too high to generate measurable electron conductivity at room temperature. For c-LLZO, the concentration of electron polarons necessary to ensure charge neutrality from positively charged oxygen vacancies formed in synthesis can be significantly higher. Hence, the electron conductivity could be significant when measured with ion-blocking metal electrodes, and we discuss how the synthesis conditions could affect this magnitude. However, in the solid-state battery, these polarons are replaced by negatively charged Li vacancies so that the electron conductivity should remain minimal. For LPS single crystals, the inherent minimal electron conductivity is independent of synthesis conditions. Here we also show that the cost of forming Li 0 in bulk c-LLZO is enormous due to strain effects so that it could only potentially form at voids, grain boundaries, or around vacancy defects which relax the lattice strain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

High Pressure Synthesis of Rubidium Superhydrides

Through laser-heated diamond anvil cell experiments, we synthesize a series of rubidium superhydrides and explore their properties with synchrotron x-ray powder diffraction and Raman spectroscopy measurements, combined with density functional theory calculations. Upon heating rubidium monohydride embedded in H 2 at a pressure of 18 GPa, we form RbH 9 − I , which is stable upon decompression down to 8.7 GPa, the lowest stability pressure of any known superhydride. At 22 GPa, another polymorph, RbH 9 − II is synthesised at high temperature. Unique to the Rb-H system among binary metal hydrides is that further compression does not promote the formation of polyhydrides with higher hydrogen content. Instead, heating above 87 GPa yields RbH 5 , which exhibits two polymorphs ( RbH 5 − I and RbH 5 − II ). All of the crystal structures comprise a complex network of quasimolecular H 2 units and H − anions, with RbH 5 providing the first experimental evidence of linear H 3 − anions. Published by the American Physical Society 2025

Kuzovnikov, Mikhail A.↗