Search NASA⌕ Search

SEARCH · Search NASA

Results for “Directed Evolution”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Dynamics of implosion, bounce, and re-implosion of a shell directly driven by a time-shaped laser pulse

Gas-filled plastic spherical shells have been imploded using the 60-beam OMEGA laser. The same laser pulse was used as in a previous experiment, which demonstrated the dynamic formation of a shell by irradiating a homogeneous plastic foam sphere. In the present experiment, the shell initially implodes, bounces, and expands; a new shell then forms and undergoes a second implosion. Analysis of x-ray self-emission-streaked images and sequences of 2D x-ray self-emission framing camera images show shell evolution that is in reasonable agreement with 1D simulations. Shell integrity and symmetry are well maintained until the formation of the new shell. The final implosion, instead, is affected by substantial asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quantifying the cooperative evolution of microphase segregation and nanostructural order in annealed polyurethanes of MDI-BDO-PTHF

Here, we quantify the relationship between microphase segregation and local structure development in thermoplastic polyurethanes. Samples were prepared with varying ratios of 4,40 -methylene diphenyl diisocyanate (MDI) and 1,4-butanediol (BDO) to polytetrahydrofuran (PTHF) and annealed at temperatures T a between RT–140 °C for 20 h each. Distinct populations of smaller and larger ordered domains are distinguished by pair distribution func tion analysis. The intermediate-range order of nanoscale hard domains in the paracrystalline state (form I) is directly extracted and characterized. The results suggest that form I of the MDI-BDO:PTHF system presents a conformational and packing order similar to the form III structure, typically obtained from stretching/annealing, but with limited spatial coherence of 3–7 nm. Heating above T g of the hard segments is necessary to achieve a substantial structural response from the annealing treatment. Increasing T a promotes purification of the soft phase via microphase segregation, coalescence of the hard blocks, and growth of paracrystalline phase, while a minority fraction of longer-range ordered domains left over from production remains relatively constant. We observe a composition-dependent decrease in the average nearest-neighbor dis tance that can be correlated with the greater number of C C, C N, and C O bonds with increased hard segment content. A non-linear deviation in the trend is observed to correlate with sample densification. The mechanical and thermal behavior of the samples is intimately tied to the segregation state.

processing-structure-properties↗
DOE OSTI2024

Nanoporous Carbon Coatings Direct Li Electrodeposition Morphology and Performance in Li Metal Anode Batteries

Li metal anodes could significantly improve battery energy density. However, Li generally electrodeposits in poorly controlled morphology, leading to safety and performance problems. One factor that controls Li anode performance and electrodeposition morphology is the nature of the electrolyte–current collector interface. Herein, we modify the Cu current collector interface by depositing precisely controlled nanoporous carbon (NPC) coatings using pulsed laser deposition to develop an understanding of how NPC coating density and thickness impact Li electrodeposition. We find that NPC density and thickness guide Li morphological evolution differently and dictate whether Li deposits at the NPC-Cu or NPC-electrolyte interface. NPC coatings generally lower overpotential for Li electrodeposition, though thicker NPC coatings limit kinetics when cycling at a high rate. Lower-density NPC enables the highest Coulombic efficiency (CE) during calendar aging tests, and higher-density NPC enables the highest CE during cycling tests.

25 ENERGY STORAGE↗
DOE OSTI2025

Recrystallization driven softening and heating rate dependencies of FeCrAl nuclear fuel cladding during accident transients

A refined understanding of FeCrAl cladding behavior during rapid transients is critical for its potential deployment in light-water reactors. Current assessments focus on transient burst testing metrics such as balloon geometry, burst temperature, and hoop stress, often used as proxies for simpler conventional tensile properties. However, directly correlating isothermal tensile and creep data with accident transient scenarios remains a challenge, although it is essential for high fidelity model development. Recent modeling based on tensile tests up to 800 °C, conducted with both immediate loading and a 10-minute soak, showed that immediate loading better predicts experimental burst temperatures, indicating a thermal softening effect. Building upon this observation, the current study connects transient performance, microstructural evolution, and high-temperature tensile properties by leveraging results from C26M claddings burst tests performed at heating rates of 1–50 °C/s and hoop stresses from 25 to 100 MPa. At 25 MPa, rupture temperatures varied by only 6 °C, but at 100 MPa, the difference reached 116 °C, with faster heating yielding higher burst temperatures. Microstructural analysis identified recrystallization as the primary cause of heating rate-dependent softening, eliminating prior cold-working. In-situ thermomechanical data linked ballooning onset to localized instabilities, similar to ultimate tensile strength behavior in conventional tensile tests. High heating rates correlated with immediate loading tensile data, while lower rates matched soaked data. Furthermore, by linking burst performance to microstructural evolution and tensile properties, this work provides a foundation for more accurate modeling of FeCrAl claddings and potentially other Fe-based materials under accident conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Evolution of Magnetohydrodynamic Turbulence in the Expanding Solar Wind: Residual Energy and Intermittency

We conduct 3D magnetohydrodynamic simulations of decaying turbulence in the context of the solar wind. To account for the spherical expansion of the solar wind, we implement the expanding box model. The initial turbulence comprises uncorrelated counterpropagating Alfvén waves and exhibits an isotropic power spectrum. Our findings reveal the consistent generation of negative residual energy whenever nonlinear interactions are present, independent of the normalized cross helicity σc and compressibility. The spherical expansion facilitates this process. The resulting residual energy is primarily distributed in the perpendicular direction, with S 2 (b) − S 2 (u) ∝ l ⊥ or equivalently –E r ∝ k$^{–2}_{⟂}$. Here S 2 (b) and S 2 (u) are second-order structure functions of magnetic field and velocity respectively. In most runs, S 2 (b) develops a scaling relation S 2 (b) ∝ $l$$^{1/2}_{⊥}$ (E b ∝ k$^{–3/2}_{⟂}$). In contrast, S 2 (u) is consistently shallower than S 2 (b), which aligns with in situ observations of the solar wind. We observe that the higher-order statistics of the turbulence, which act as a proxy for intermittency, depend on the initial σc and are strongly affected by the expansion effect. Generally, the intermittency is more pronounced when the expansion effect is present. Finally, we find that in our simulations, although the negative residual energy and intermittency grow simultaneously as the turbulence evolves, the causal relation between them seems to be weak, possibly because they are generated on different scales.

Interplanetary turbulence↗
DOE OSTI2024

Cosmology of dark energy radiation

In this work, we quantify the cosmological signatures of dark energy radiation—a novel description of dark energy, which proposes that the dynamical component of dark energy is comprised of a thermal bath of relativistic particles sourced by thermal friction from a slowly rolling scalar field. For a minimal model with particle production emerging from first principles, we find that the abundance of radiation sourced by dark energy can be as large as Ω DER = 0.03, exceeding the bounds on relic dark radiation by three orders of magnitude. Although the background and perturbative evolution of dark energy radiation are distinct from Quintessence, we find that current and near-future cosmic microwave background and supernova data will not distinguish these models of dark energy. We also find that our constraints on all models are dominated by their impact on the expansion rate of the Universe. Considering extensions that allow the dark radiation to populate neutrinos, axions, and dark photons, we evaluate the direct detection prospects of a thermal background comprised of these candidates consistent with cosmological constraints on dark energy radiation. Our study indicates that a resolution of ~6 meV is required to achieve sensitivity to relativistic neutrinos compatible with dark energy radiation in a neutrino capture experiment on tritium. We also find that dark matter axion experiments lack sensitivity to a relativistic thermal axion background, even if enhanced by dark energy radiation, and dedicated search strategies are required to probe new parameter space. We derive constraints arising from a dark photon background from oscillations into visible photons, and find that viable parameter space can be explored with the late dark energy radiation experiment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Acceleration of Thermochemistry Solves in MOOSE and Pronghorn

This work focuses on the development and implementation of strategies to accelerate thermochemical calculations within MOOSE-based multiphysics simulations, particularly for applications in MSRs. We highlight the inherent complexity of nuclear materials, which require a multiscale approach to accurately model their behavior across various physical domains, including mechanical, chemical, and thermal phenomena. Thermochemical equilibrium calculations are crucial for predicting material properties and enhancing the fidelity of these simulations. The integration of Thermochimica, a Gibbs energy minimizer, into MOOSE allows for the direct minimization of Gibbs energy at every point on the mesh. However, the computational cost of such integration is significant. To address this, we explored acceleration strategies such as multi-threading support and the use of a thermodynamic ValueCache to reduce redundant calculations. Additionally, we investigated modifications to Thermochimica to enable phase constraints and improve its coupling with phase-field models, which are essential for simulating microstructural evolution and corrosion in MSR. These efforts aim to optimize the computational efficiency and accuracy of multiphysics simulations, thereby supporting the development of reliable and efficient nuclear materials for next-generation reactor technologies.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Lewis Acid Site Engineering in Chromite Spinels Orchestrated Surface Reconstruction and Surpasses RuO 2 in Oxygen Evolution

Atomic-scale engineering of chromite spinels featuring redox-active tetrahedral A-sites and strong Cr–O covalency offers a promising route to superior platinum-group-metal-free oxygen evolution reaction (OER) catalysts. However, comprehensive studies addressing how cation substitution influences surface chemistry and governs OER activity and durability in chromite spinels remain limited. Here, in this work, a systematic investigation of the multicationic chromite series Ni x Fe y Cr 3−x−y O 4 is presented, identifying composition-dependent Lewis acidity as a descriptor of superior OER performance. It is further demonstrated that tuning surface acidity directly controls dynamic reconstruction processes and lattice-oxygen participation during spinel-based electrocatalysis. Following activation, the optimized Ni 0.8 Fe 0.3 Cr 1.9 O 4 catalyst delivers a current density of 10 mA cm −2 at an overpotential of 235 mV, surpassing RuO 2 , with excellent long-term stability. Integrating microscopic and spectroscopic analysis with operando impedance spectroscopy, it shows that activation generates an oxyhydroxide overlayer and reveals a previously unrecognized link between surface Lewis acidity and the growth kinetics and activity of these shells. Density functional theory calculations indicate that Fe incorporation at octahedral sites raises the O 2p-band center and lowers oxygen-vacancy formation energy, promoting lattice-oxygen activation and triggering reconstruction, yielding enhanced OER. This work integrates cation-driven surface-acidity modulation, acidity-governed reconstruction, and OER activity enhancement into a unified predictive framework for designing earth-abundant spinel-based catalysts.

operando impedance spectroscopy↗
DOE OSTI2024

Neutron Diffraction Residual Stress Study of a AA2219-T87 Self-Reacting Friction Stir Weld

2219-T87 is a precipitation hardenable aluminum-copper alloy which sees wide use in structural aerospace components. Thick panels of this alloy are joined via self-reacting friction welding (SRFSW); however, this thermomechanical process can result in significant loss of the alloy’s strengthening precipitates and large residual stresses which are detrimental to mechanical behavior and performance. High-resolution maps of the residual strain and stress states along the normal, transverse, and longitudinal directions of 2219-T87 SRFSW were obtained using neutron diffraction measurements. Residual stress had the highest tensile value in the heat-affected zone (HAZ) and the largest compressive stress in the base metal region of the normal and transverse directions. Line profile residual stress distributions displayed an “M”-shaped distribution in each direction with the longitudinal direction being the most pronounced while 2D residual stresses displayed an hourglass-shaped pattern. Further, average maximal longitudinal residual stress values ranged from 60.6 to 85.6% of yield. Microhardness testing across the transverse weld section produced a “V”-shaped curve with slight hardness recovery in the stir zone. Maximal and minimal microhardness values were observed in the base metal (154 HV) and thermomechanically affected zone (82 HV), respectively. Microstructural evolution was recorded using optical microscopy and showed decreasing grain size from the HAZ to the stir zone.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Probing the Kitaev honeycomb model on a neutral-atom quantum computer

Quantum simulations of many-body systems are among the most promising applications of quantum computers. In particular, models based on strongly correlated fermions are central to our understanding of quantum chemistry and materials problems, and can lead to exotic, topological phases of matter. However, owing to the non-local nature of fermions, such models are challenging to simulate with qubit devices. Here we realize a digital quantum simulation architecture for two-dimensional fermionic systems based on reconfigurable atom arrays. We utilize a fermion-to-qubit mapping based on Kitaev’s model on a honeycomb lattice, in which fermionic statistics are encoded using long-range entangled states. We prepare these states efficiently using measurement and feedforward, realize subsequent fermionic evolution through Floquet engineering with tunable entangling gates interspersed with atom rearrangement, and improve results with built-in error detection. Leveraging this fermion description of the Kitaev spin model, we efficiently prepare topological states across its complex phase diagram and verify the non-Abelian spin-liquid phase by evaluating an odd Chern number. We further explore this two-dimensional fermion system by realizing tunable dynamics and directly probing fermion exchange statistics. Finally, we simulate strong interactions and study the dynamics of the Fermi–Hubbard model on a square lattice. These results pave the way for digital quantum simulations of complex fermionic systems for materials science, chemistry and high-energy physics.

atomic and molecular physics↗
DOE OSTI2024

An introduction to Spent Nuclear Fuel decay heat for Light Water Reactors: a review from the NEA WPNCS

This paper summarized the efforts performed to understand decay heat estimation from existing spent nuclear fuel (SNF), under the auspices of the Working Party on Nuclear Criticality Safety (WPNCS) of the OECD Nuclear Energy Agency. Needs for precise estimations are related to safety, cost, and optimization of SNF handling, storage, and repository. The physical origins of decay heat (a more correct denomination would be decay power) are then introduced, to identify its main contributors (fission products and actinides) and time-dependent evolution. Due to limited absolute prediction capabilities, experimental information is crucial; measurement facilities and methods are then presented, highlighting both their relevance and our need for maintaining the unique current full-scale facility and developing new ones. The third part of this report is dedicated to the computational aspect of the decay heat estimation: calculation methods, codes, and validation. Different approaches and implementations currently exist for these three aspects, directly impacting our capabilities to predict decay heat and to inform decision-makers. Finally, recommendations from the expert community are proposed, potentially guiding future experimental and computational developments. One of the most important outcomes of this work is the consensus among participants on the need to reduce biases and uncertainties for the estimated SNF decay heat. If it is agreed that uncertainties (being one standard deviation) are on average small (less than a few percent), they still substantially impact various applications when one needs to consider up to three standard deviations, thus covering more than 95% of cases. The second main finding is the need of new decay heat measurements and validation for cases corresponding to more modern fuel characteristics: higher initial enrichment, higher average burnup, as well as shorter and longer cooling time. Similar needs exist for fuel types without public experimental data, such as MOX, VVER, or CANDU fuels. A third outcome is related to SNF assemblies for which no direct validation can be performed, representing the vast majority of cases (due to the large number of SNF assemblies currently stored, or too short or too long cooling periods of interest). A few solutions are possible, depending on the application. For the final repository, systematic measurements of quantities related to decay heat can be performed, such as neutron or gamma emission. This would provide indications of the SNF decay heat at the time of encapsulation. For other applications (short- or long-term cooling), the community would benefit from applying consistent and accepted recommendations on calculation methods, for both decay heat and uncertainties. This would improve the understanding of the results and make comparisons easier.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An integrated modeling framework with open architecture for phase field simulation of multi-component alloys

An integrated modeling framework (PanPhaseField) has been developed, which enables a direct and fast coupling between CALPHAD calculations and large-scale phase field simulations for multi-component alloys. Further, it adopts an open architecture allowing for integration of user-defined phase field models in a plug-and-play manner by taking full advantage of the user-friendly graphical interface of Pandat software. The developed modeling platform becomes an enabling tool that can be used to simulate the evolution of spatially varying microstructures of industrial complex alloys for various engineering applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Impact of Urbanization on Convection, Lightning, and Precipitation Over the Houston Metropolitan Area: Case Study Simulation From the TRACER Campaign

This study investigates the effects of urbanization, specifically land use change and anthropogenic emissions (AE), on convection, lightning, and surface precipitation for a case of summertime sea‐breeze convection observed over the Houston metropolitan area. The unique capabilities of the NASA‐Unified Weather Research and Forecasting model allows us to conduct a series of sensitivity experiments with complex configurations, in particular including multi‐year land model spin‐up simulations, treatment of aerosols and their precursors, and explicit cloud charging and lightning. The simulation results show that urban land use primarily alters the temporal evolution of convection, lightning, and surface precipitation, leading to late afternoon thunderstorm development. The decrease in latent heat flux from the land surface caused by urbanization weakens convection in the early afternoon, while a condition suitable for convection development is maintained in the late afternoon due to less stabilization of the lower troposphere by the weaker convection development and high sensible heat flux from the surface. On the other hand, anthropogenic aerosols directly enhance convection, lightning, and surface precipitation by increasing convective updrafts due to the aerosol‐induced convective invigoration. The combined effects of urban land use and AE lead to even stronger thunderstorms in the late afternoon, mostly consistent with observations. These results indicate that urbanization increases the probability of late afternoon thunderstorms over the Houston area during the summer season. Advanced weather forecasting models that incorporate these urbanization effects might support sustainable urban planning to better mitigate the impacts of urbanization on local weather and public safety.

Iguchi, T. [Univ. of Maryland, College Park, MD (U↗
DOE OSTI2025

Global heliospheric termination shock strength in the solar–interstellar interaction

A heliospheric termination shock (HTS) surrounds our Solar System at approximately 100 astronomical units from the Sun, where the expanding solar wind (SW) is compressed and heated before encountering the interstellar medium. HTS-accelerated particles govern the pressure balance with the interstellar medium, but little is known about the global properties of the HTS beyond in situ measurements from Voyager in only two directions of the sky. Here we fill this gap by extracting the HTS strength using particle-in-cell, test particle and magnetohydrodynamic simulations, constrained by Interstellar Boundary Explorer observations of energetic neutral atoms produced from HTS-accelerated particles. Our results reveal there is a higher compression near the poles during solar minimum compared with solar maximum due to the higher Mach number flow. North–south asymmetries arise from the disparate evolution of the polar coronal holes, while minimum compression near the flanks is probably due to SW slowing from mass loading over a greater distance to the HTS. The results imply a strong connection between the HTS strength and the SW and interstellar medium dynamics.

Astrophysical plasmas↗
DOE OSTI2026

From Spin to Pseudospin Symmetry: The Origin of Magic Numbers in Nuclear Structure

Magic numbers lie at the heart of nuclear structure, reflecting enhanced stability in nuclei with closed shells. While the emergence of magic numbers beyond 20 is commonly attributed to strong spin-orbit coupling, the microscopic origin of the spin-orbit potential remains elusive, owing to its dependence on the resolution scale and renormalization scheme of nuclear forces. Here, we investigate the evolution of nuclear shell structure with varying momentum resolution in nuclear interactions derived from chiral effective field theory, using the similarity renormalization group to link different scales. We uncover a novel transition from spin symmetry to pseudospin symmetry as the resolution scale decreases, during which magic numbers emerge naturally. A similar pattern is found in calculations using relativistic one-boson-exchange potentials, underscoring the robustness of the phenomenon. This establishes a direct connection between realistic nuclear forces with a high resolution scale and effective nuclear forces at coarse-grained scales, offering a first-principles explanation for the origin of magic numbers and pseudospin symmetry in nuclear shell structure and new insights into the structure of exotic nuclei far from stability

Energy levels↗
DOE OSTI2026

Final DOE-ASR Report for the Project “Using LASSO to bridge the gap between model and observations and to learn about atmospheric convection”

Atmospheric convection spans a wide range of spatial and temporal scales and involves complex interactions with the surrounding dynamic and thermodynamic environment, particularly over tropical continental regions. These processes remain a major source of uncertainty in weather and climate models, including persistent biases in the diurnal cycle of convective precipitation that directly affect estimates of climate sensitivity. Addressing these challenges requires the combined use of high-resolution observations and cloud-resolving modeling frameworks. In this context, the DOE Atmospheric Radiation Measurement (ARM) program’s Large-Eddy Simulation ARM Symbiotic Simulation and Observation (LASSO) activity provides a powerful platform that pairs comprehensive observations with numerical simulations to enable process-level understanding of atmospheric convection. Within this context, this Research and Development Partnership Pilot (RDPP) project was designed to initiate and expand DOE ARM/ASR research capacity at minority-serving institutions, while advancing scientific understanding of convective processes over the Amazon rainforest. Consistent with the RDPP mission, the project emphasized partnership development, training, and workforce capacity building alongside exploratory research activities. On the scientific side, the project produced two peer-reviewed journal articles, and one manuscript currently under review (see list in section 3.1). Together, these studies combine long-term ARM observations and cloud-resolving and convection-permitting modeling to investigate the environmental controls on the shallow-to-deep convective transition during the Amazon wet season. The results demonstrate the central role of early-day moisture preconditioning and large-scale dynamical forcing in regulating isolated deep convection, provide mechanistic insight into convective evolution, and establish physically informed modeling frameworks for future sensitivity experiments. These scientific outcomes are described in sections 2.1 to 2.3 and were disseminated in 8 conference presentations (see section 3.2) and 5 invited talks (see section 3.3), reflecting broad engagement with our community. Equally important, the project achieved its RDPP capacity-building objectives (see section 2.4). A sustained research partnership was established among the University of Maryland, Baltimore County (UMBC), Morgan State University (MSU), and Howard University (HU), and extended to include collaboration with Pacific Northwest National Laboratory (PNNL). The project organized multiple multi-day training events focused on ARM data, LASSO simulations, and quantitative analysis methods, directly engaging students, postdoctoral researchers, and faculty across institutions. These activities broadened participation in ASR research and led to independent adoption of LASSO workflows by students beyond the immediate project team. Finally, the project successfully positioned the participating institutions to pursue future DOE research. Preliminary scientific results, coupled with strengthened partnerships and technical capacity, enabled the submission of follow-on proposals to DOE ASR funding opportunities. In this way, the project fulfilled the RDPP goal of seeding durable research capacity and laying the foundation for larger-scale, sustained engagement with DOE ARM and ASR programs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Impact of cement composition, brine concentration, diffusion rate, reaction rate and boundary condition on self-sealing predictions for cement-CO 2 systems

Geological CO 2 storage (GCS) plays an important role in curbing CO 2 emissions by reducing the carbon footprint of difficult to decarbonize operations and achieve negative CO 2 emissions through activities like Bioenergy with Carbon Capture and Storage (BECCS) and Direct Air Carbon Capture and Storage (DACCS). Leakage of CO 2 through wells is an important concern when it comes to deployment of large-scale GCS. Existing wells in sites that are otherwise suitable for GCS can act as conduits for stored CO 2 to escape the reservoir. There is broad consensus that the main risk of leakage through wellbores is via fractures/damaged pathways. The results from several studies evaluating the permeability evolution of cement fractures in wells upon leakage of CO 2 agree that smaller fracture apertures, slower brine velocities and higher brine residence times promote self-sealing of fractures by mineral precipitation. Quantitatively, however, the differences in sealing conditions are significant and are typically attributed to differences in experimental conditions or model assumptions. Here we examine the sensitivity of our model, describing CO 2 leakage through wellbores, to cement composition, brine concentration, diffusion rates, and reaction rates. We also evaluate the impact of the boundary condition to allow comparisons between observations from experiments performed at constant flow rate and model predictions made at constant pressure conditions. Our results show that diffusion and reactions rates have the most impact on the self-sealing criteria for cement-CO 2 systems. In addition, conditions associated with self-sealing of fractures at constant flow rate require longer fractures, smaller fracture apertures, and slower velocities than under constant pressure.

58 GEOSCIENCES↗