Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical balances”

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 OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2023

Unlocking the price

To help balance load and power generation that ensue from the electrification of transportation and the increased connection of variable power sources to the grid, ABB has developed an RL-based dynamic price model for EV charging with the required flexibility to respond to changing grid conditions.

Suryanarayana, Harish↗
DOE OSTI2025

Probing exciton diffusion dynamics in photosynthetic supercomplexes via exciton–exciton annihilation

Photosynthesis converts solar energy into chemical energy through coordinated energy transfer between light-harvesting complexes and reaction centers (RCs). Understanding exciton motion, particularly the exciton diffusion length, is essential for optimizing energy efficiency in photosystems. In this work, we combine intensity-cycling transient absorption spectroscopy with kinetic Monte Carlo (kMC) simulation to investigate exciton motion in the C2S2 photosystem II supercomplex of spinach. Using exciton–exciton annihilation, revealed in the fifth-order response, we experimentally estimate an exciton diffusion length of 10.9 nm based on a 3D normal diffusion model, suggesting the ability of excitons to traverse the supercomplex. However, kMC simulations reveal that exciton motion is sub-diffusive because of spatial constraints and the strong RC traps. An anomalous diffusion model analysis of the experimental data yields a diffusion length of 9.7 nm, while the simulated diffusion length is 7.4 nm. The variable exciton residence time across subunits, partly influenced by their connectivity to the trap, indicates inhomogeneous annihilation probability and suggests how plants balance efficient light harvesting with photoprotection. We also explore the influence of specific assumptions in the annihilation simulation, which are challenging to access in more complex environments, such as the thylakoid membrane. Our study provides a framework for studying exciton dynamics using exciton–exciton annihilation, which can be extended to understand the light-harvesting efficiencies of larger, more complex photosynthetic assemblies.

Zhang, Kunyan (ORCID:000000026830409X)↗
DOE OSTI2022

Electric Access System Enhancement (EASE): Assessment of a Distributed Energy Resource Management System for Enabling Dynamic Hosting Capacity

The Electric Access System Enhancement (EASE) project demonstrated a scalable, interoperable, and cost-effective means of integrating high penetration of distributed energy resources (DERs). The control architecture developed leveraged a Distributed Energy Resource Management System, Distribution System Operator for transacting energy, and 3rd Party DER Aggregator platform. The project identified ways to enhance the customer interconnection process to the grid and improve access to information from DERs and optimize the usage of DERs to provide energy services in a simulated day-ahead shadow market while maintain grid reliability. By integrating these capabilities into a scalable system of systems, the DSO can effectively balance DER generation and customer demand on the distribution network. This capability allows the grid to host more DER than traditionally possible on wires alone. Hosting more DER has the added benefit of supplying increased demand growth, sometimes beyond the capacity limits of the distribution network itself. This is known as a dynamic hosting capacity, which could help utilities manage the forecasted growth in electricity demand as California switch to electric vehicles and appliances. This could help to establish energy storage and PV generation as a “pseudo firm” generation resource mix if managed appropriately and provide sufficient resource adequacy for distribution capacity upgrade deferrals.

14 SOLAR ENERGY↗
DOE OSTI2021

Magnetic excitations of the hybrid multiferroic (ND 4 ) 2 FeCl 5 ·D 2 O

We report a comprehensive inelastic neutron scattering study of the hybrid molecule-based multiferroic compound ( ND 4 ) 2 FeCl 5 · D 2 O in the zero-field incommensurate cycloidal phase and the high-field quasicollinear phase. The spontaneous electric polarization changes its direction concurrently with the field-induced magnetic transition, from mostly aligned with the crystallographic a axis to the c axis. To account for such a change in polarization direction, the underlying multiferroic mechanism was proposed to switch from the spin-current model induced via the inverse Dzyaloshinskii-Moriya interaction to the p - d hybridization model. We perform a detailed analysis of the inelastic neutron data of ( ND 4 ) 2 FeCl 5 · D 2 O using linear spin-wave theory to quantify magnetic interaction strengths and investigate the possible impact of different multiferroic mechanisms on the magnetic couplings. Overall our result reveals that the spin dynamics of both multiferroic phases can be well described by a Heisenberg Hamiltonian with easy-plane anisotropy. We do not find notable differences between the optimal model parameters of the two phases. The hierarchy of exchange couplings and the balance among frustrated interactions remain the same between two phases, suggesting that magnetic interactions in ( ND 4 ) 2 FeCl 5 · D 2 O are much more robust than the electric polarization in response to delicate reorganizations of the electronic degrees of freedom in an applied magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Probe Chiral Magnetic Effect with Signed Balance Function

In this paper a pair of observables are proposed as alternative ways, by examining the fluctuation of net momentum-ordering of charged pairs, to study the charge separation induced by the Chiral Magnetic Effect (CME) in relativistic heavy ion collisions. They are, the out-of-plane to in-plane ratio of fluctuation of the difference between signed balance functions measured in pair’s rest frame, and the ratio of it to similar measurement made in the laboratory frame. In this work, both observables have been studied with simulations including flow-related backgrounds, and for the first time, backgrounds that are related to resonance's global spin alignment. The two observables have similar positive responses to signal, and opposite, limited responses to identifiable backgrounds arising from resonance flow and spin alignment. Both observables have also been tested with two realistic models, namely, a multi-phase transport (AMPT) model and the anomalous-viscous fluid dynamics (AVFD) model. These two observables, when cross examined, will provide useful insights in the study of CME-induced charge separation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Carbon Dioxide Production in Bedrock Beneath Soils Substantially Contributes to Forest Carbon Cycling

Abstract Soils are widely considered the primary terrestrial organic matter pool mediating carbon transactions with the atmosphere and groundwater. Because soils are both a host and a product of rhizosphere activity, they are thought to mark the location where photosynthetic fixation of carbon dioxide (CO 2 ) is balanced by the oxidation of organic matter. However, in many terrestrial environments, the rhizosphere extends below soils and into fractured bedrock, and it is unknown if the resulting biological and hydrologic dynamics in bedrock have a significant impact on carbon cycling. Here we show substantial production of CO 2 in weathered bedrock at 4–8 m below the thin soils (<0.5 m thick) of a Northern California forest using innovative monitoring technology for sampling gases and water in fractured rock. The deep CO 2 production supports a persistent upward flux of CO 2(g) year‐round from bedrock to soil, constituting between 2% and 29% of the average daily CO 2 efflux from the land surface. When water is rapidly transported across the fractured bedrock vadose zone, nearly 50% of the CO 2 produced in the bedrock dissolves into water, promoting water‐rock interaction and export of dissolved inorganic carbon (DIC) from the unsaturated zone to groundwater, constituting as much as 80% of the DIC exiting the hillslope. Such CO 2 production in weathered bedrock is subject to unique moisture, temperature, biological, and mineralogical conditions which are currently missing from terrestrial carbon models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Description and Demonstration of the Coupled Community Earth System Model v2 – Community Ice Sheet Model v2 (CESM2‐CISM2)

Abstract Earth system/ice‐sheet coupling is an area of recent, major Earth System Model (ESM) development. This work occurs at the intersection of glaciology and climate science and is motivated by a need for robust projections of sea‐level rise. The Community Ice Sheet Model version 2 (CISM2) is the newest component model of the Community Earth System Model version 2 (CESM2). This study describes the coupling and novel capabilities of the model, including: (1) an advanced energy‐balance‐based surface mass balance calculation in the land component with downscaling via elevation classes; (2) a closed freshwater budget from ice sheet to the ocean from surface runoff, basal melting, and ice discharge; (3) dynamic land surface types; and (4) dynamic atmospheric topography. The Earth system/ice‐sheet coupling is demonstrated in a simulation with an evolving Greenland Ice Sheet (GrIS) under an idealized high CO 2 scenario. The model simulates a large expansion of ablation areas (where surface ablation exceeds snow accumulation) and a large increase in surface runoff. This results in an elevated freshwater flux to the ocean, as well as thinning of the ice sheet and area retreat. These GrIS changes result in reduced Greenland surface albedo, changes in the sign and magnitude of sensible and latent heat fluxes, and modified surface roughness and overall ice sheet topography. Representation of these couplings between climate and ice sheets is key for the simulation of ice and climate interactions.

Muntjewerf, Laura↗
DOE OSTI2023

Warming decreases desert ecosystem functioning by altering biocrusts in drylands

Warming and precipitation fluctuations are changing desert ecosystems in global drylands. However, the effects of climate change on keystone species such as cryptogamic biocrust in drylands remain relatively under-investigated, even though biocrusts play a vital role in desert ecosystems. Here we conducted a long-term experiment (14 years) to simulate the responses of two main types of biocrusts to warming coupled with reduction in precipitation that was achieved by open-top chambers (OTCs) to simulate the predicted warming and precipitation decreasing under climate change scenario. We also conduct a structural analysis to evaluate the resulting changes in desert ecosystem functioning (carbon and nitrogen cycling). Neither warming and corresponding rainfall reduction treatments had a negative effect on lichen species richness, but both treatments reduced lichen cover and biomass. The negative effects of warming on moss-dominated crusts were much greater than those on lichen-dominated crusts. Although mosses and lichens had varying degree responses to warming, the loss of mosses and decreased lichen cover and biomass, as well as the shortening of the wet time, resulted in a reduction in carbon and nitrogen fixation, soil enzyme activity and water-holding capacity of biocrusts and topsoil. These impacts collectively change the water balance of drylands and weaken the hydrological and biogeochemical function of biocrusts. Synthesis and applications: Results from this long-term experiment suggest that the ecosystem C and N cycling and water balance of global drylands may be highly impacted by climate change, in part because of the response of biocrusts, which contribute an important implication for both dryland restoration and earth system dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modeling the small-scale deposition of snow onto structured Arctic sea ice during a MOSAiC storm using snowBedFoam 1.0.

The remoteness and extreme conditions of the Arctic make it a very difficult environment to investigate. In these polar regions covered by sea ice, the wind is relatively strong due to the absence of obstructions and redistributes a large part of the deposited snow mass, which complicates estimates for precipitation hardly distinguishable from blowing or drifting snow. Moreover, the snow mass balance in the sea ice system is still poorly understood, notably due to the complex structure of its surface. Quantitatively assessing the snow distribution on sea ice and its connection to the sea ice surface features is an important step to remove the snow mass balance uncertainties (i.e., snow transport contribution) in the Arctic environment. In this work we introduce snowBedFoam 1.0., a physics-based snow transport model implemented in the open-source fluid dynamics software OpenFOAM. We combine the numerical simulations with terrestrial laser scan observations of surface dynamics to simulate snow deposition in a MOSAiC (Multidisciplinary Drifting Observatory for the Study of Arctic Climate) sea ice domain with a complicated structure typical for pressure ridges. The results demonstrate that a large fraction of snow accumulates in their vicinity, which compares favorably against scanner measurements. However, the approximations imposed by the numerical framework, together with potential measurement errors (precipitation), give rise to quantitative inaccuracies, which should be addressed in future work. The modeling of snow distribution on sea ice should help to better constrain precipitation estimates and more generally assess and predict snow and ice dynamics in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Understanding Adsorption and Reactions at Aqueous Oxide Interfaces with Neural Network Potential Molecular Dynamics

Chemical processes at metal oxide−water interfaces are of central importance in geochemistry, biology, and energy technologies. A better understanding of these processes would allow us to make a significant step toward optimizing and controlling them, which could in turn lead to broader impacts. Computational modeling is indispensable to accomplishing this task because complexity and disorder often make it difficult to extract atomistic information from experiments. Balancing computational cost and accuracy, simulation schemes based on efficient machine learning representations of the potential energy surface (PES) predicted by ab initio calculations have become increasingly popular over the past decade. In particular, several studies have demonstrated the ability of machine learning models to accurately reproduce the complex ab initio PESs of aqueous oxide interfaces, allowing simulations of systems and processes that are not accessible with ab initio methods. In this Account, we review our recent efforts to understand adsorption processes and reactions at aqueous oxide interfaces using deep potential molecular dynamics (DPMD), a simulation scheme employing deep neural networks (DNNs), which has proven to be quite successful in accurately describing many different systems in the condensed phase. After summarizing the DPMD methodology, we first review our work on the acid−base chemistry of oxide surfaces in contact with water, a fundamental characteristic that controls proton transfer and surface charge at the interface. We focus on the aqueous interface of rutile IrO 2 , an oxide material thus far considered the best catalyst for the oxygen evolution reaction (OER). We show that this interface is characterized by a large fraction of dissociated water and a strong Brønsted acidity of the surface sites, in good agreement with the experimentally measured value of the point of zero proton charge. In our second example, we investigate how the adsorption of organic species from ambient air or water affects the structure and wettability of the aqueous interfaces of TiO 2 , a prototypical photocatalytic material. This is a question that is relevant to understanding the UV-induced hydrophilicity of TiO 2 surfaces, a property at the basis of self-cleaning windows and related applications. Specifically focusing on formic and acetic acids, the two most common atmospheric organic acids, our simulations reveal that these acids control the wettability of TiO 2 largely through acid−base chemistry at the interface rather than chemisorption on the oxide surface, a finding that could help improve the design of self-cleaning surfaces and photocatalytic devices. Finally, we review our recent study of methanol at TiO 2 −water interfaces, a system whose interest is largely motivated by the role of methanol in enhancing photocatalytic hydrogen evolution on TiO 2 . Our simulations provide mechanistic insights into the coupled roles of the organic adsorbate and water at the TiO 2 interface, with implications for how methanol enhances the activity of H 2 evolution.

adsorption↗
DOE OSTI2024

Permafrost Carbon: Progress on Understanding Stocks and Fluxes Across Northern Terrestrial Ecosystems

Significant progress in permafrost carbon science made over the past decades include the identification of vast permafrost carbon stocks, the development of new pan-Arctic permafrost maps, an increase in terrestrial measurement sites for CO 2 and methane fluxes, and important factors affecting carbon cycling, including vegetation changes, periods of soil freezing and thawing, wildfire, and other disturbance events. Process-based modeling studies now include key elements of permafrost carbon cycling and advances in statistical modeling and inverse modeling enhance understanding of permafrost region C budgets. By combining existing data syntheses and model outputs, the permafrost region is likely a wetland methane source and small terrestrial ecosystem CO 2 sink with lower net CO 2 uptake toward higher latitudes, excluding wildfire emissions. For 2002–2014, the strongest CO2 sink was located in western Canada (median: -52 g C m -2 y -1 ) and smallest sinks in Alaska, Canadian tundra, and Siberian tundra (medians: -5 to -9 g C m -2 y -1 ). Eurasian regions had the largest median wetland methane fluxes (16–18 g CH4 m -2 y -1 ). Quantifying the regional scale carbon balance remains challenging because of high spatial and temporal variability and relatively low density of observations. More accurate permafrost region carbon fluxes require: (a) the development of better maps characterizing wetlands and dynamics of vegetation and disturbances, including abrupt permafrost thaw; (b) the establishment of new year-round CO 2 and methane flux sites in underrepresented areas; and (c) improved models that better represent important permafrost carbon cycle dynamics, including non-growing season emissions and disturbance effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Comprehensive Characterization of Water Group Ion Composition and Distributions in Saturn's Magnetosphere With Cassini Plasma Spectrometer Data

Saturn's magnetosphere is continuously supplied with neutrals from the Enceladus plume and the icy rings, which undergo ionization and charge-exchange to form a complex water-group plasma environment. While the Cassini Plasma Spectrometer (CAPS) instrument has provided extensive compositional information, detailed separation of individual water-group ion species in time-of-flight (TOF) data has not previously been achieved. In this study, we perform forward modeling of CAPS-IMS energy-per-charge (E/Q) and TOF spectra obtained between 2004 and 2012 to resolve O + , OH + , H 2 O + , and H 3 O + and to characterize their plasma properties, including number density, temperature, and thermodynamic κ. Our results demonstrate that O + is the dominant thermal ion species throughout Saturn's magnetosphere, comprising up to ∼70% of the total ion population beyond ∼5 Saturn radii (R S ). In contrast, molecular ions such as OH + , H 2 O + , and H 3 O + dominate closer to Enceladus but rapidly dissociate into atomic ions between ∼5 and 10 R S . This radial region is also characterized by the steepest increase in plasma flow speed, which rises from ∼40% to ∼80% of rigid corotation. Simultaneously, ion velocity distributions approach Maxwell–Boltzmann equilibrium, as indicated by high kappa values. These findings provide new constraints on the ion–neutral chemistry that regulates the balance between molecular and atomic ions in Saturn's magnetosphere. They also emphasize the critical role of the 5–10 R S region as a transition zone for both plasma composition and dynamics. Our results refine previous CAPS-based studies and underscore the need to incorporate seasonal variability and ionospheric coupling into future global models of Saturn's plasma environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Defining a Platform Approach and Market Participation: Data Driven Business Models for Solid State Transformer-Based Synthetic Inertia and Voltage Stability Controls (CRADA Final Report, Project 1, Mod 1)

The primary objective of this project is to determine the incremental value created with the medium voltage solid-state transformer (MV SST) technology to different stakeholders in view of the updated DER grid regulations. This includes studying the benefits of the MV SST technology in a range of use cases for EV and DER penetration including (1) “corridor charging” for EVs and (2) solar plus storage (FERC 2222). The potential customers of this technology include utilities for EV charging, DER installers who must meet utility interconnection requirements, balancing authorities, and DER aggregators. The traditional transformers on the grid could be a limiting factor for the EV-grid integration as the distribution transformers were not designed to handle the dynamic and fluctuating EV charging loads. Thus, the issues such as voltage fluctuations, increased losses and reduced efficiency [1] can negatively impact the grid operation. To address these challenges, transformers with flexibility and adaptability become imperative to meet the evolving energy demands. In this regard, the concept of Medium Voltage Solid-State Transformers.

14 SOLAR ENERGY↗
DOE OSTI2026

Simulation of Divertor Performance in ST40 Under Dynamic Double-Null Plasmas

A power fraction model was implemented for the simultaneous prediction of 3-D surface temperature evolution at all four divertor targets in near-double-null (DN) tokamak configurations, which is especially important for compact high-field devices that may not have the ability to dissipate large amounts of power on the high-field side. Evaluating the power-sharing between the four divertor strike points in a disconnected DN configuration is important for understanding the overall power balance, as well as for optimizing the power exhaust performance and prolonging the survivability of the plasma facing components (PFCs). This power-sharing is typically evaluated in terms of the separation between the primary and secondary separatrices at the outboard midplane, $\textit {dR}_{\text {sep}}$. The Heat flux Engineering Analysis Toolkit (HEAT) is coupled with Brunner’s power fraction model to simulate the deposited heat flux and resultant temperature change on 3-D divertor targets in a dynamic DN (DDN) pulse operation in ST40, a high-field spherical tokamak. The simulation results showed that with DDN operation, the operation time has significantly increased compared with single-null geometry configurations.

ST40↗
DOE OSTI2026

Contrasting carbon dynamics in grazed and flood-prone grasslands on mineral and degraded peat soils

Ecosystem-scale methane (CH₄) flux measurements from grazed grasslands remain scarce, despite their importance for understanding grassland contributions to the global carbon budget. In this study, we present full annual budgets of both carbon dioxide (CO₂) and CH₄ derived from eddy covariance measurements at two contrasting grazed grasslands: a floodplain grassland at Marchegg, Austria, and an intensively grazed pasture at Sherman Barn, California. By combining continuous, year-round observations of CO₂ and CH₄, this study provides a rare, comparative assessment of greenhouse gas dynamics across distinct climatic, hydrological, and management regimes. Our results highlight how environmental conditions, grazing intensity, and hydrology jointly regulate CO₂ exchange and CH₄ emissions, underscoring the importance of including methane alongside net ecosystem exchange when evaluating grassland carbon balances. At Marchegg, characterized by seasonal flooding and moderate horse grazing on mineral soils, the ecosystem acted as a weak carbon sink in 2024, with a GWP 100 of 42.2 ± 113.4 g CO₂ eq m⁻². Net CO₂ uptake (–27.3 ± 30.2 g C m⁻² yr⁻¹) was partly offset by CH₄ emissions (1.6 ± 0.04 g C m⁻² yr⁻¹), which were strongly linked to soil moisture and inundation events. In contrast, Sherman Barn, a cattle pasture on degraded peat soils, was a consistent carbon source in 2019, with a GWP 100 of 567.6 ± 120.5 g CO₂ eq m⁻², associated with high ecosystem respiration during summer. The site released 125.3 ± 32.2 g C m⁻² yr⁻¹ of CO₂ and 3 ± 0.06 g C m⁻² yr⁻¹ of CH₄. Across both sites, net ecosystem exchange was primarily linked to photosynthetically active radiation and vegetation greenness, while CH₄ fluxes were related to soil moisture rather than grazing intensity. Flooding at Marchegg reduced CO₂ uptake but enhanced CH₄ emissions, highlighting the critical role of flood timing within the growing season. Moreover, inundation appeared to suppress the spread of invasive plant species, emphasizing the ecological value of dynamic hydrological regimes. Together, these findings reveal the complexity of grassland carbon budgets and the need for site-specific, year-round CO₂ and CH₄ monitoring to inform climate-adaptive management strategies.

Carbon sequestration↗
DOE OSTI2026

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗
DOE OSTI2021

Favorable morphology and electronic conductivity of functional sublayers for highly efficient water splitting electrodes

Low electronic conductivities and improper morphologies of anode electrodes greatly limit the reaction area, catalyst utilization and efficiency in proton exchange membrane water electrolyzers. In this study, conductive sublayers with different conductivities and morphologies were introduced into anode electrodes in membrane-based water electrolyzers. In-situ and ex-situ investigation results showed that conductive sublayers (Au mesh and carbon nanotube (CNT) film) augmented the sheet conductivity of anode electrodes by up to 4000 times (from 2000 to 0.5 ohm square -1 ), and the ohmic resistance of water electrolyzers was reduced to 1/3 when inserting conductive sublayers. In addition, CNT film provided a higher electrochemical active area than Au mesh, because of favorable morphologies (large porosity and surface area) of CNT fibers on CNT films. Therefore, the current density of water splitting was increased by 3 times (from 4.55 mA cm -2 to 14.83 mA cm -2 ) at 2.5 V compared to a conventional anode electrode. Visualizations on bubble dynamics showed improved performances with conductive sublayers; this was mainly due to greatly increased number of reaction sites, highly spread reaction area (from 50 to 1000 um), and reduced activation overpotential. Therefore, a balance between high electronic conductivity and nanoporous morphology is essential to the anode electrode for larger reaction sites and areas in highly efficient water electrolyzers.

25 ENERGY STORAGE↗