Search NASA⌕ Search

SEARCH · Search NASA

Results for “Evolution operators”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

A stochastic reorganizational bath model for electronic energy transfer

Environmentally induced fluctuations of the optical gap play a crucial role in electronic energy transfer dynamics. One of the simplest approaches to incorporate such fluctuations in energy transfer dynamics is the well known Haken-Strobl-Reineker (HSR) model, in which the energy-gap fluctuation is approximated as white noise. Recently, several groups have employed molecular dynamics simulations and excited-state calculations in conjunction to account for excitation energies’ thermal fluctuations. On the other hand, since the original work of HSR, many groups have employed stochastic models to simulate the same transfer dynamics. Here, we discuss a rigorous connection between the stochastic and the atomistic bath models. If the phonon bath is treated classically, time evolution of the exciton-phonon system can be described by Ehrenfest dynamics. To establish the relationship between the stochastic and atomistic bath models, we employ a projection operator technique to derive the generalized Langevin equations for the energy-gap fluctuations. The stochastic bath model can be obtained as an approximation of the atomistic Ehrenfest equations via the generalized Langevin approach. Based on this connection, we propose a novel scheme to take account of reorganization effects within the framework of stochastic models. The proposed scheme provides a better description of the population dynamics especially in the regime of strong exciton-phonon coupling. Finally, we discuss the effect of the bath reorganization in the absorption and fluorescence spectra of ideal J-aggregates in terms of the Stokes shifts. We find a simple expression that relates the reorganization contribution to the Stokes shifts – the reorganization shift – to the ideal or non-ideal exciton delocalization in a J-aggregate. The reorganization shift can be described by three parameters: the monomer reorganization energy, the relaxation time of the optical gap, and the exciton delocalization length. This simple relationship allows one to understand the physical origin of the Stokes shifts in molecular aggregates.

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

Pathways Toward Efficient and Durable Anion Exchange Membrane Water Electrolyzers Enabled By Electro‐Active Porous Transport Layers

Abstract Green hydrogen, produced via water electrolysis using renewable electricity, will play a crucial role in decarbonizing industrial and heavy‐duty transportation sectors. Anion exchange membrane water electrolyzers (AEMWEs) can overcome many of the performance and cost limitations of incumbent technologies, however, still suffer from durability challenges due to oxidative instability of anion‐exchange ionomers. Herein, the use of an electro‐active porous transport layer as anode (PTL‐electrode) is demonstrated to enable efficient and durable AEMWEs. The stainless‐steel PTL‐electrodes are shown to have superior performance and durability compared to traditional catalyst layers containing ionomer and nanoparticle catalysts. An AEMWE cell operating at 2 A cm −2 for over 600 h exhibited a degradation rate of just 5 µV h −1 . During operation, the surface composition of the stainless steel transforms into a mixture of iron and nickel oxyhydroxides, contributing to enhanced oxygen‐evolution reaction activity. The combination of experimental work and modeling elucidates how the bulk structure of the PTL‐electrode offers an additional design dimension to further improve electrolyzer performance. Lastly, a surface modification strategy is applied to a PTL‐electrode to achieve an even higher performing AEMWE (2.3 vs 2.0 A cm −2 at 1.8 V). Overall, this work lays out pathways toward more efficient, durable, and affordable AEMWEs.

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

Cosmological moduli and non-perturbative production of axions

Cosmological moduli generically come to dominate the energy density of the early universe, and thereby trigger an early matter dominated era. Such non-standard cosmological histories are expected to have profound effects on the evolution and production of axion cold dark matter and dark radiation, as well as their prospects for detection. We consider moduli-axion couplings and investigate the early history of the coupled system, considering closely the evolution of the homogeneous modulus field, the back-reaction from the axion, and the energy densities of the two fields. A particular point of interest is the enhancement of axion production from modulus decay, due to tachyonic and parametric resonant instabilities, and the implications of such production on the cosmological moduli problem, axion dark radiation, and the available parameter space for axion dark matter. Using an effective field theory approach, WKB-based semi-analytical analysis, and detailed numerical estimates of the co-evolution of the system, we evaluate the expected decay efficiency of the modulus to axions. The effects of higher-order operators are studied and implications for UV-complete frameworks such as the Large Volume Scenarios in Type IIB string theory are considered in detail.

Axions and ALPs↗
DOE OSTI2022

Effect of fuel composition and EGR on spark-ignited engine combustion with LPG fueling: Experimental and numerical investigation

This paper presents an experimental and numerical investigation of a spark-ignited (SI) cooperative fuel research (CFR) engine fueled with different liquefied petroleum gas (LPG) fuels and exhaust gas recirculation (EGR). Here, the effects of LPG fuel composition on engine combustion characteristics are initially evaluated at two different compression ratios (CR). Results show normal combustion at CR 7 and heavy knocking combustion at CR 10 for all the tested fuels, with a more substantial impact for the LPG fuel with high proportions of n-butane species. The Livengood-Wu (LW) integral method is then used to analyze the knock occurrence risk of individual fuel based on the reactivity of the tested fuels. The introduction of EGR then demonstrates the potential of knock intensity reduction below the borderline knock limit. A zonal-based kinetic interactions study is also performed to understand the knock mitigation effectiveness of EGR over the pressure–temperature domain relevant to SI engine operation. Finally, a multidimensional, computational fluid dynamics (CFD) simulation model is shown to predict the LPG combustion characteristics and presents the evolution of in-cylinder temperature and chemical species to demonstrate the development of end-gas autoignition events without and with EGR.

02 PETROLEUM↗
DOE OSTI2025

X-ray Micro-Computed Tomography for Structural Analysis of All-Solid-State Battery at Pouch Cell Level

Characterizing the microstructure of all-solid-state batteries (ASSBs) during fabrication and operation is vital for their advancement, particularly as scaling to pouch cell levels introduces challenges in probing large-scale microstructural evolution. This work highlights the potential of synchrotron X-ray micro-computed tomography (sXCT) as a nondestructive, rapid (<30 min), and high-resolution technique for visualizing and quantifying key microstructural features, including overhang, porosity, contact loss, active surface area, and tortuosity, in all-solid-state pouch cells. The large field of view (up to millimeters) of sXCT enables detailed analysis at an industry-relevant scale, bridging the gap between laboratory research and commercial applications. Furthermore, integrating realistic sXCT-derived 3D models into multiphysics simulations could provide insights into chemo-mechanical degradation, particularly at the edges of the pouch cells, offering a pathway for designing robust, high-performance ASSBs. This perspective establishes sXCT as an indispensable tool for advancing both the understanding and the engineering of next-generation energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2015

Effects of electrolyte, catalyst, and membrane composition and operating conditions on the performance of solar-driven electrochemical reduction of carbon dioxide

Solar-driven electrochemical cells can be used to convert carbon dioxide, water, and sunlight into transportation fuels or into precursors to such fuels. The voltage efficiency of such devices depends on the (i) physical properties of its components (catalysts, electrolyte, and membrane); (ii) operating conditions (carbon dioxide flowrate and pressure, current density); and (iii) physical dimensions of the cell. The sources of energy loss in a carbon dioxide reduction (CO2R) cell are the anode and cathode overpotentials, the difference in pH between the anode and cathode, the difference in the partial pressure of carbon dioxide between the bulk electrolyte and the cathode, the ohmic loss across the electrolyte and the diffusional resistances across the boundary layers near the electrodes. We analyze the effects of these losses and propose optimal device configurations for the efficient operation of a CO2R electrochemical cell operating at a current density of 10 mA cm -2 . Cell operation at near-neutral bulk pH offers not only lower polarization losses but also better selectivity to CO2R versus hydrogen evolution. Addition of supporting electrolyte to increase its conductivity has a negative impact on cell performance because it reduces the electric field and the solubility of CO 2 . Addition of a pH buffer reduces the polarization losses but may affect catalyst selectivity. The carbon dioxide flowrate and partial pressure can have severe effects on the cell efficiency if the carbon dioxide supply rate falls below the consumption rate. The overall potential losses can be reduced by use of an anion, rather than a cation, exchange membrane. We also show that the maximum polarization losses occur for the electrochemical synthesis of CO and that such losses are lower for the synthesis of products requiring a larger number of electrons per molecule, assuming a fixed current density. We also find that the reported electrocatalytic activity of copper below -1 V vs. RHE is strongly influenced by excessive polarization of the cathode and, hence, does not represent its true activity at bulk conditions. This article provides useful guidelines for minimizing polarization losses in solar-driven CO2R electrochemical cells and a method for predicting polarization losses and obtaining kinetic overpotentials from measured partial current densities.

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

Nonuniqueness of generalized quantum master equations for a single observable

When deriving exact generalized master equations for the evolution of a reduced set of degrees of freedom, one is free to choose what quantities are relevant by specifying projection operators. However, obtaining a reduced description does not always need to be achieved through projections-one can also use conservation laws for this purpose. Such an operation should be considered as distinct from any kind of projection; that is, projection onto a single observable yields a different form of master equation compared to that resulting from a projection followed by the application of a constraint. Here, we give a simple example to show this point and give relationships that the different memory kernels must satisfy to yield the same dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

SU(2) non-Abelian gauge field theory in one dimension on digital quantum computers

A dynamical quantum simulation of SU(2) non-Abelian gauge field theory on a digital quantum computer is presented here. This was enabled on current quantum hardware by introducing a mapping of the field onto a register of qubits that utilizes local gauge symmetry while preserving local constraints on the fields, reducing the dimensionality of the calculation. Controlled plaquette operators and gauge-variant completions in the unphysical part of the Hilbert space were designed and used to implement time evolution. The new techniques developed in this work generalize to quantum simulations of higher dimensional gauge field theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Complex Structure of Molten FLiBe (2 Li F – Be F 2 ) Examined by Experimental Neutron Scattering, X-Ray Scattering, and Deep-Neural-Network Based Molecular Dynamics

The use of molten salts as coolants, fuels, and tritium breeding blankets in the next generation of fission and fusion nuclear reactors benefits from furthering the characterization of the molecular structure of molten halide salts, paving the way to predictive capability of the chemical and thermophysical properties of molten salts. Due to its neutronic, chemical, and thermochemical properties, 2 Li F - Be F 2 is a candidate molten salt for several fusion- and fission-reactor designs. We performed neutron and x-ray total-scattering measurements to determine the atomic structure of liquid 2 Li F - Be F 2 . We also performed and neural-network molecular-dynamics simulations to predict the structure obtained by neutron- and x-ray-diffraction experiments. The use of machine learning provides improvements to the efficiency in predicting the structure at a longer length scales than is achievable with simulations at significantly lower computational expense while retaining near accuracy. We found that the NNMD simulations accurately predicted the Be F 4 2 − oligomer formations seen in the experimental first-structure-factor peak. Our combination of high-resolution measurements with large-scale molecular dynamics provided an avenue to explore and experimentally verify the intermediate-range ordering beyond the first-nearest neighbor that has posed too many experimental and computational challenges in previous works. With a deeper understanding of the salt structure and ion ordering, the evolution of salt chemistry over the lifetime of a reactor can be better predicted, which is crucial to the licensing and operation of advanced fission and fusion reactors that employ molten salts. To this end, this work will serve as a reference for future studies of salt structure and macroscopic properties with and without the addition of solutes. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Computational Elucidation of Mechanical Degradation in NMC Cathodes: Impact on Cell Performance

Next-generation lithium ion batteries are expected to demonstrate superior energy and power density with longer cycle life for successful electrification of the automobile, aviation, and marine industries. Adoption of lithium metal anodes with solid electrolytes can help to achieve that goal given that the dendrite-related issues are solved eventually. Another possibility is to use Ni-rich high-capacity NMC cathode materials with liquid and/or solid electrolytes, which presently experiences rapid capacity fade while charged to higher voltages. Several mechanical and chemical degradation mechanisms are active within these NMC-based cathode particles. Recent experimental research activities attempted to correlate the mechanical damage with the capacity fade experienced by Ni-rich LiNi x Mn y Co z O 2 (x+y+z = 1) (NMC) cathodes. Here, a computational framework is developed in this study capable of quantifying the evolution of inter primary particle and cathode/electrolyte interfacial fracture experienced by the poly- and single-crystalline NMC cathodes during charge/discharge operation. Influences of mechanical degradation on the overall cell capacity, while operating with liquid and/or solid electrolytes, are successfully characterized. Decreasing the size of the cathode primary particles, or the size of the single-crystalline cathodes, can mitigate the overall mechanical degradation, and subsequent capacity fade, experienced by NMC cathodes. The developed theoretical methodology can help the engineers and scientists to better understand the mechanical degradation mechanism prevalent in Ni-rich NMC cathodes and build superior lithium ion-based energy storage devices for the application in next-generation devices.

25 ENERGY STORAGE↗
DOE OSTI2016

Organ and Tissue-specific Sucrose Transporters. Important Hubs in Gene and Metabolite Networks Regulating Carbon Use in Wood-forming Tissues of Populus

The overall project objective was to probe the relationship between sucrose transporters and plant productivity in the biomass for biofuels woody perennial, Populus. At the time the proposal was written, sucrose transporters had already been investigated in many plant model systems, primarily with respect to the export of photosynthate sucrose from source leaves, and the uptake of sucrose in storage organs and seeds. Preliminary findings by the PI found that in Populus, sucrose transporter genes (SUTs) were well expressed in wood-forming tissues that comprise the feedstock for biofuels production. Because sucrose comprises by far the predominant form in which photosynthate is delivered from source organs to sink organs like roots and wood-forming tissues, SUTs control a gate that nominally at least could impact the allocation or partitioning of sucrose for potentially competing end uses like growth (stem biomass) and storage. In addition, water use might be conditioned by the way in which sucrose is distributed throughout the plant, and/or by the way in which sucrose is partitioned intracellularly. Several dozen transgenic lines were produced in year 1 of the project to perturb the expression ratio of multiple plasma membrane (PM) SUTs (intercellular trafficking), versus the single tonoplast membrane (TM) sucrose transporter that effectively regulates intracellular trafficking of sucrose. It was possible to obtain transgenic lines with dual SUT gene knockdown using the 35S promoter, but not the wood-specific TUA1 promoter. By the end of project year 2, a decision was made to work with the 35S plants while archiving the TUA1 plants. The PhD candidate charged with producing the transgenic lines abandoned the project during its second year, substantially contributing to the decision to operate with just the 35S lines. That student’s interests ranged more toward evolutionary topics, and a report on SUT gene evolution was published (Peng et al 2014).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Analysis Report for Exascale Storage Requirements for Scientific Data.

Over the next 10 years, the Department of Energy will be transitioning from Petascale to Exascale Computing resulting in data storage, networking, and infrastructure requirements to increase by three orders of magnitude. The technologies and best practices used today are the result of a relatively slow evolution of ancestral technologies developed in the 1950s and 1960s. These include magnetic tape, magnetic disk, networking, databases, file systems, and operating systems. These technologies will continue to evolve over the next 10 to 15 years on a reasonably predictable path. Experience with the challenges involved in transitioning these fundamental technologies from Terascale to Petascale computing systems has raised questions about how these will scale another 3 or 4 orders of magnitude to meet the requirements imposed by Exascale computing systems. This report is focused on the most concerning scaling issues with data storage systems as they relate to High Performance Computing- and presents options for a path forward. Given the ability to store exponentially increasing amounts of data, far more advanced concepts and use of metadata will be critical to managing data in Exascale computing systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Cesium Ion Exchange with Crystalline Silicotitanate Literature Review

Crystalline silicotinante (CST) is planned for use as the ion exchange media to accomplish the cesium (Cs) decontamination of low-activity waste (LAW) in the current evolution of the Direct Feed Low-Activity Waste flowsheet. CST has yet to be deployed in a production-scale operating waste treatment facility on the Hanford Site. However, CST was used nearly 20 years ago to successfully treat approximately 270,000 gallons of supernatant waste in modular columns at Oak Ridge National Laboratory. As part of the technical maturation of the media for use at Hanford, an assessment of the state-of-the-art was performed by reviewing work in the technical literature. The review considered work from the U.S. Department of Energy (DOE) complex (i.e., national laboratories and government documents), the CST manufacturer, and the open technical literature nationally and internationally. The review focused on information most relevant to the caustic, high-sodium chemical environment typical of Hanford LAW. The topic areas explored in this literature review document include the following: CST structural and physical properties; Isotherm and batch contact data; Mass transport: intraparticle, interparticle, dispersion in the column; Column performance at multiple scales; Hydraulic performance; The effect of process conditions; Gas generation and retention behavior; Post-processing options for spent CST media The reviewed literature demonstrates that CST is well-suited for Cs decontamination of Hanford LAW. However, after consideration of the available literature data in these subject areas, some uncertainties regarding the application of CST for the treatment of Hanford LAW are identified and summarized. These uncertainties are areas where additional information may permit better prediction of CST performance or provide margin to facility operations and safety basis analyses.

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

Characteristics of Ultralow Aspect Ratio Toroidal Plasmas and Non-Solenoidal Startup and Edge Stability Studies at Near-Unity Aspect Ratio (Final Report)

This grant supported a long-term research and development program to explore the physics and technology of spherical tokamak (ST) plasma behavior at near-unity aspect ratio ( A ). The PEGASUS Toroidal Experiment is a university-scale ( R ) ~ 0.40 m, α ~ .35 m, I p ≤ 0.3 MA, B T ≤ 0.15 T) experiment that accesses A → 1, to provide large values of the normalized plasma current, I N = I p /α B T , or equivalently the toroidal field utilization factor I p / I TF . The research activities focused on: 1) developing the technique of Local Helicity Injection (LHI) as a viable means of initiating and growing an ST plasma without the use of a central solenoid; 2) exploring the equilibrium and stability properties of spherical tokamak plasmas at extremely low aspect ratio; and 3) obtaining new insights in the properties of the enhanced confinement H-mode regime at A ~ 1. The initial (Phase I) facility employed a unique high-stress central solenoid Ohmic for current drive and heating with simple magnet current waveforms. A subsequent major upgrade (Phase II) offered more flexible control of the plasma evolution along with a powerful LHI system for plasma initiation, heating and sustainment. The LHI technique consists of injection of high-density current streams into vacuum toroidal and vertical magnetic fields. These streams become unstable and support magnetic reconnection to allow the system to relax into a tokamak-like magnetic configuration. This program focused on: 1) development of electron current sources that provide sufficiently high injection currents in the challenging edge plasma region; 2) validation of theoretical limits to achievable I p arising from magnetic helicity and energy conservation plus absolute limits arising from relaxation to a minimal energy states; 3) tests of a power balance model to describe the evolution of the plasma current I p (t); and 4) study of magnetic fluctuations and their role in the I p evolution. This culminated in the production of toroidal currents of I p ≥ 0.2 MA with only ~8 kA of injected current. Operation with Ohmic induction alone in Phase I created plasmas with soft limits on the achievable I p and field utilization ( I p / I TF ~ 1) due to the onset of large low-order (2/1 or 3/2) internal resonant tearing modes. With the expanded capabilities in Phase II operation, methods employed to mitigate those modes and achieve I p / I TF > 1 were: 1) I TF ramp-downs; 2) non-solenoidal startup via LHI; and 3) extremely low I TF Ohmic operations assisted by LHI current injector pre-ionization. Extensive use of LHI-driven startup and drive for non-solenoidal operation produced favorable current profiles to stably operate tokamak plasmas with unprecedentedly high I N = I p / aB T ~ 15 and I p / I TF > 2. With strong reconnection-driven anomalous ion heating, this in turn provided access to tokamak plasmas with B t up to 100%, a world record 2 times higher than that previously achieved in advanced or spherical tokamaks. This high B t regime had a large fraction (~50%) of the plasma volume in a region of an absolute minimum- B configuration. The MHD stability limit in this regime was the ideal external kink mode without a stabilizing wall. These results provided the first demonstration of access to tokamak plasmas with very high B t , high i N , high elongation κ, and low internal inductance ℓ i , and as such served as a proof-ofprinciple for access to this long-sought ST regime. The naturally low B T needed for stability at A ≤ 1.2 reduces the L-H mode power threshold, P LH , so that H-mode plasmas were achieved with Ohmic heating alone. The power threshold for the transition to the H-mode regime was found to be 10× the accepted values from the ITER international scaling, confirming initial trends indicated by experiments on MAST and NSTX, and suggesting that the physics behind this transition is not yet fully understood. The density dependence of the transition and the observed insensitivity to whether the plasma is bounded by a limiter or magnetic divertor are consistent with the recent FM 3 model. Multiple unstable toroidal magnetic modes were measured during ELM (Edge Localized Mode) crashes, with generally lower toroidal mode numbers than seen at higher A , which was attributed to the larger peeling mode instability drive at low A . First-ever edge current profile measurements through an ELM crash showed current-carrying filaments were generated, as seen in ELM simulations and in electromagnetic blob transport theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Local limits of detection for anthropogenic aerosol-cloud interactions

Ship tracks are quasi-linear cloud patterns produced from the interaction of ship emissions with low boundary layer clouds. They are visible throughout the diurnal cycle in satellite images from space-borne assets like the Advanced Baseline Imagers (ABI) aboard the National Oceanic and Atmospheric Administration Geostationary Operational Environmental Satellites (GOES-R). However, complex atmospheric dynamics often make it difficult to identify and characterize the formation and evolution of tracks. Ship tracks have the potential to increase a cloud's albedo and reduce the impact of global warming. Thus, it is important to study these patterns to better understand the complex atmospheric interactions between aerosols and clouds to improve our climate models, and examine the efficacy of climate interventions, such as marine cloud brightening. Over the course of this 3-year project, we have developed novel data-driven techniques that advance our ability to assess the effects of ship emissions on marine environments and the risks of future marine cloud brightening efforts. The three main innovative technical contributions we will document here are a method to track aerosol injections using optical flow, a stochastic simulation model for track formations and an automated detection algorithm for efficient identification of ship tracks in large datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Existing Methods for Grid Strength Assessment and Role of Hydropower in Future Grids

The power system is undergoing rapid evolution with increasing penetration of inverter-based resources (IBRs), large loads, microgrids, and power electronic devices. Ensuring reliable and stable operation of the modern power grid is a multi-faceted challenge that requires detailed understanding of this complex system. Dynamic stability is a major concern in maintaining the security of power grids as the generation mixes and large loads transitions to include high shares of power electronic devices. The operation of IBRs in regions with low system strength and higher grid impedances has been found to be the main reason behind many of the power system instabilities that manifest themselves in various types of oscillations and interactions that, if not properly addressed and damped, can jeopardize the reliable operation of the power system. Weak grid conditions compound such stability problems, particularly when many IBRs operate in proximity to and connect to weak power grids. Therefore, it is important to assess grid strength for planning, integration and operation of IBRs for a given power system. In this report, we aim to understand and classify existing methodologies for grid strength assessment, along with their limitations and future needs. Further, we aim to utilize the huge untapped potential in utilizing hydro energy resources in addressing some of the pertinent challenges of grid strength and reliable operation of modern power systems. Hydropower, historically valued for its flexibility and dispatchability, now faces new constraints due to reduced share of synchronous machines in the generation mix and seasonal variability of available water resources. Yet, these same plants present untapped potential beyond energy generation - notably, as providers of critical grid services. This report explores how hydropower plants, particularly through operation as synchronous condensers, can play a pivotal role in strengthening the grid amid evolving system dynamics.

13 HYDRO ENERGY↗
DOE OSTI2025

A new metrics framework for quantifying and intercomparing atmospheric rivers in observations, reanalyses, and climate models

We present a new atmospheric river (AR) analysis and benchmarking tool, namely Atmospheric River Metrics Package (ARMP). It includes a suite of new AR metrics that are designed for quick analysis of AR characteristics via statistics in gridded climate datasets such as model output and reanalysis. This package can be used for climate model evaluation in comparison with reanalysis and observational products. Integrated metrics such as mean bias and spatial pattern correlation are efficient for diagnosing systematic AR biases in climate models. For example, the package identifies the fact that, in CMIP5 and CMIP6 (Coupled Model Intercomparison Project Phases 5 and 6) models, AR tracks in the South Atlantic are positioned farther poleward compared to ERA5 reanalysis, while in the South Pacific, tracks are generally biased towards the Equator. For the landfalling AR peak season, we find that most climate models simulate a completely opposite seasonal cycle over western Africa. This tool can also be used for identifying and characterizing structural differences among different AR detectors (ARDTs). For example, ARs detected with the Mundhenk algorithm exhibit systematically larger size, width, and length compared to the TempestExtremes (TE) method. The AR metrics developed from this work can be routinely applied for model benchmarking and during the development cycle to trace performance evolution across model versions or generations and set objective targets for the improvement of models. They can also be used by operational centers to perform near-real-time climate and extreme event impact assessments as part of their forecast cycle.

58 GEOSCIENCES↗
DOE OSTI2017

Modeling of CO 2 sequestration in coal seams: Role of CO 2 -induced coal softening on injectivity, storage efficiency and caprock deformation: Original Research Article: Modeling of CO 2 sequestration in coal seams

An effective and safe operation for sequestration of CO 2 in coal seams requires a clear understanding of injection-induced coupled hydromechanical processes such as the evolution of pore pressure, permeability, and induced caprock deformation. In this study, CO 2 injection into coal seams was studied using a coupled flow-deformation model with a new stress-dependent porosity and permeability model that considers CO 2 -induced coal softening. Based on triaxial compression tests of coal samples extracted from the site of the first series of enhanced coalbed methane field tests in China, a softening phenomenon that a substantial (one-order-of-magnitude) decrease of Young's modulus and an increase of Poisson's ratio with adsorbed CO 2 content was observed. Such softening was considered in the numerical simulation through an exponential relation between elastic properties (Young's modulus and Poisson's ratio) and CO 2 pressure considering that CO 2 content is proportional to the CO 2 pressure. Our results of the numerical simulation show that the softening of the coal strongly affects the CO 2 sequestration performance, first by impeding injectivity and stored volume (cumulative injection) during the first week of injection, and thereafter by softening mediated rebound in permeability that tends to increase injectivity and storage over the longer term. A sensitivity study shows that stronger CO 2 -induced coal softening and higher CO 2 injection pressure contribute synergistically to increase a significant increase of CO 2 injectivity and adsorption, but also result in larger caprock deformations and uplift. This study demonstrates the importance of considering the CO 2 -induced softening when analyzing the performance and environmental impact of CO 2 -sequestration operations in unminable coal seams.

54 ENVIRONMENTAL SCIENCES↗