Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “unit dynamics”

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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Lithium-ion Battery Strain Gauge Monitoring and Depth of Discharge Estimation

Lithium-ion battery diagnostics and prognostics rely on measurements of electrical impedance, capacity, and voltage to infer the internal state of the battery. Mechanical changes to the cell structure represent an additional measure of the battery’s state because these changes are related to the overall battery health. As lithium-ion batteries are charged and discharged, lithium ions are inserted or removed from the anode and cathode, a process called intercalation and de-intercalation. As lithium ions intercalate and de-intercalate, they can cause changes to the lattice of the electrode particles, resulting in volumetric changes. These volumetric changes cause mechanical stresses and strains on the lithium-ion battery electrodes, and subsequently, the whole cell’s thickness varies as it is charged and discharged. This paper presents a study on the use of surface-mounted strain gauges for in-situ measurement of structural changes to lithium-ion batteries, along with a characterization of the unit-to-unit differences in strain response. A neural network modeling structure is then used to predict the battery’s depth of discharge under dynamic discharge conditions.

Lithium-ion batteries, depth of discharge, strain ↗
NASA NTRSDate not supplied

Development of 2.05 µm Fiber Lasers for CO2 DIAL Lidar Measuring Martian CO2 and Pressure

The Decadal Survey by the National Academy of Sciences in the United States and NASA’s Science Mission Directorate (SMD) Science Plan both require global CO 2 observations for the Martian atmospheric pressure, dynamics and chemistry. However, there are significant observational gaps in polar regions and during nighttime for global air pressure and CO 2 on Mars. Therefore, recently we proposed a new concept of Martian differential absorption lidar (DIAL) operating in the 2.05 µm CO 2 absorption band for global, including poles, atmospheric CO 2 and pressure observations day and night [1]. Based on the concept, we are awarded to develop 2.05 µm fiber lasers by the NASA Planetary Instrument Concepts for the Advancement of Solar System Observations (PICASSO) Program. The laser is designed to be an all-fiber master oscillator and power amplifier (MOPA) system with laser output of ~3 mJ at a repetition frequency of 2 kHz. The primary master oscillator (PMO) is locked to the center (2.0504280 µm) of the selected CO 2 absorption line. The frequency of a second MO (SMO) is locked to that of PMO. This SMO frequency is adjustable and switched between the online (2.05044156 µm) and offline (2.05050812 µm) wavelengths. The online wavelength is optimally selected so that the CO 2 absorption optical depth (AOD) is ~1.1 at 3 km and the measurement in the low Martian atmosphere has largest signal-to-noise ratio. The online wavelength can also be adjusted to a line slope location where AOD is larger to observe atmospheric pressure at higher altitudes. We will present more detail about this project and instrument development at the conference.

Zhaoyan Liu↗
DOE OSTI2025

Asymptotic-preserving dynamical low-rank method for the stiff nonlinear Boltzmann equation

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only r 2 times, where r, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

FARM supervisory capabilities for thermal energy storage

The FARM (Feasible Actuator Range Modifier) module is a component of the RAVEN-based FORCE framework for the analysis of Integrated Energy Systems (IES). FARM aids HERON in the solution of the power dispatch problem by evaluating feasible set-point signals to be issued to the control systems of the different IES unit components. Set-points need to satisfy limits on both production variables (i.e., the variables to be optimized such as the electrical power, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). To enforce all these limits, a two-stage approach is adopted. First, the power dispatcher algorithm in HERON module estimates set-points that meet the constraints on the production variables, e.g., power levels and power ramp rate limits. These constraints are called explicit constraints. Then, if necessary, FARM adjusts these set-points to ensure the respect of the limits on the process variables of interest, given the knowledge of the system dynamics acquired through machine learning algorithms. These constraints are called implicit constraints. From this standpoint, FARM constitutes a bridge between the HERON power dispatcher that adopts a simplified description of the IES unit (low-resolution physics) and the HYBRID high-fidelity models (high-resolution physics). In this report, an overview of the major capabilities of the latest release of FARM is provided, along with a summary of the tool demonstration campaign conducted at the Dynamic Energy Technology and Integration Laboratory (DETAIL) facility. These results assess the performance of the control system architecture embedding FARM both as a Validator of the HERON power dispatcher and as a real time Supervisory control scheme. Additionally, the report outlines the areas that FARM might benefit from, along with proposed solutions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

An atomistic study connecting underlying dislocation behavior with superior mechanical properties of NiCoCr medium entropy alloy

NiCoCr-based medium-entropy alloy (MEA) with a simple face-centered cubic crystal phase exhibits excellent mechanical properties, often attributed to the synergy of multiple deformation mechanisms. However, the atomistic origin of their outstanding mechanical response, including microstructural evolution and dislocation behavior under varying strain-rates and orientation, remains unclear. In this work, we employ large-scale molecular dynamics (MD) simulations to investigate the changes in deformation mechanisms along three distinct orientations ([110], [111], [100]) under varying strain rates (1 ×10 8 /sec, 1 ×10 10 /sec, 1 ×10 12 /sec) in the NiCoCr MEA. The presence of the stair-rod and the Shockley partial dislocations under uniaxial tensile strain are found to play a key role in the formation of deformation twinning and ε-martensite, which positively correlates with strain-rate dependent dislocation analysis. These findings further establish the role of the dislocations in controlling the superior mechanical response and excellent fracture toughness of the NiCoCr MEA. Systematic transmission-electron microscopy tests performed on the [111]-oriented crystals, deformed at different strain levels, at room temperature provide clear evidence of both the extended stacking-fault and the stair rods, confirming the predicted microstructural features. Finally, this study offers key insights into the complex nucleation mechanisms of deformation twinning and ε-martensite, such as twinning – and transformation–induced plasticity (TWIP-TRIP), providing valuable guidelines for studying similar material classes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Trichomonas vaginalis C2-XYPPX-repeat protein with a structured C2 domain displaying dampened flexibility upon binding calcium

C2 domains are ubiquitous membrane-binding modules of ∼130 residues in eukaryotes that are often associated with proteins involved in membrane trafficking and lipid modification. The genome of Trichomonas vaginalis, the most common, non-viral, sexually transmitted human pathogen, encodes eight genes that contain a N-terminal C2 module linked to a XYPPX-repeat domain of more than four XYPPX repeats (C2-XYPPX). While the function of the XYPPX-repeat domain remains unknown, its multiple association with C2 domains in T. vaginalis suggests it is important. Here, the C2 domain from one of these C2-XYPPX-repeat proteins, Tv-C2-1, was structurally and physically characterized using X-ray crystallography and NMR spectroscopy. The crystal structure for Tv-C2-1 shows that this domain shares a fold common to all C2 domains, a compact Greek-key motif composed of eight anti-parallel β-strands in the type-2 topology. An NMR chemical shift perturbation study with Ca 2+ showed that Tv-C2-1 bound two Ca 2+ atoms primarily via two loops (loop-1 and loop-3) on the predicted calcium binding face of the protein with K d s of 58.0 ± 0.1 μM and 232 ± 6 μM. Estimations of the overall rotational correlation time, τ c , in the apo (11.1 ns) and Ca 2+ -bound (9.2 ns) state suggests the protein becomes more compact upon Ca 2+ binding, consistent with a decrease in dynamics in loop-3 and marginally in loop-1 suggested by amide 15 N heteronuclear steady-state { 1 H}- 15 N NOEs. Showing Tv-C2-1 binds calcium and adopts a compact Greek-key motif structure, two primary features of C2 domains, suggests understanding the function of the XYPPX-repeat domain may be warranted.

NMR spectroscopy↗
DOE OSTI2026

Multibody for Everybody (M4E): A Symbolic Dynamics Modeling Tool with Applications in Simulation, Control, and Optimization

Developing the analytical model of a multibody system is often the initial step in control and optimization. The analytical model (equations of motion) describes a system’s time evolution under specified forcing conditions. Although developing these equations is easy for simple systems, this process becomes more complex for systems composed of multiple bodies. Deriving equations of motion for complex multibody systems requires specialized expertise in multibody dynamics, is time-consuming, and is susceptible to error. To address this issue, this paper presents an open-source, easy-to-use, systematic framework to derive symbolic equations of motion in both Python and MATLAB using the joint coordinate formulation. This formulation results in a set of ordinary differential equations that use the minimum set of coordinates needed to model a system. The symbolic representation provides better insight into the influence of design parameters on system performance, facilitates sensitivity analysis and parameter studies, and supports direct implementation of control and optimization routines. The tool enables numerical simulation for specified parameter sets, is modular for straightforward integration with other tools and libraries, and allows incorporation of hydrodynamics, mooring, and other external forces. The result is a reproducible, extensible pipeline for modeling, simulation, and design of complex multibody systems. The proposed tool is versatile and can be applied to domains such as robotics, control, and design. In addition, we integrated external libraries that provide capabilities for modeling offshore systems such as underwater robots and marine energy converters.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

CMIP6-based Multi-model Hydropower Projection over the Conterminous US, Version 1.1

This dataset presents a suite of hydropower projections for the conterminous United States (CONUS), derived from multiple downscaled and bias-corrected Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). The CMIP6 GCMs are downscaled using either statistical (DBCCA) or dynamical (RegCM) approaches, based on two meteorological reference datasets (Daymet and Livneh). The resulting downscaled precipitation, temperature, and wind speed data are then used to drive two calibrated hydrologic models (VIC and PRMS), enabling simulations of projected future hydrologic responses across the CONUS. Simulated total runoff is subsequently employed to drive two hydropower models (WMP, now implemented as mosartwmpy-power, and WRES) to evaluate how climate change may affect future hydropower production for both federal and non-federal hydropower fleets. This dataset was developed to support the SECURE Water Act Section 9505 Assessment for the U.S. Department of Energy (DOE) Water Power Technologies Office (WPTO). For further details, see Broman et al. (2024), Thurber et al. (2024), Kao et al. (2022), and Zhou et al. (2023).

Voisin, Nathalie [Pacific Northwest National Labor↗
DOE OSTI2025

Exploring Novel Configuration Designs Alternatives for Lazy Wave Submarine Power Cables

Floating offshore wind turbines require a complex integrated design to balance the platform response and station-keeping requirements in a highly dynamic environment. Dynamic subsea power cables are one of the critical components that have to be considered in detail to ensure both extreme, as well as repeated operational conditions do not exceed the cable design limits. This paper presents the design feasibility limits of a 66kV dynamic cable for a lazy wave configuration. The work explores an innovative staggered cable configuration that has the potential to provide additional compliance in shallow water conditions. Both configurations are evaluated using the motions of a representative floating wind platform under development for deployment in relatively shallow water depth of 150m. The assessment was conducted using the hydrodynamic software Orcaflex, incorporating a range of design inputs and design variations. Key parameters analysed include cable curvature and effective tension and bending stress. The alternative cable configuration provides a more gradual crossing of the water column, which would allow to support increased platform excursions in shallow water depths. The design uses a combined buoyancy and bend protection element, offering designers additional options for challenging operational conditions. The paper will be of interest to designers and developers of floating wind projects, providing a comparative assessment between the conventional lazy wave configuration and an innovative alternative that enables larger cable and mooring excursion envelopes in challenging conditions.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Analog and symbolic computation through the Koopman framework

We develop a Koopman operator framework for studying the computational structure of dynamical systems. Specifically, we show that the resolvent of the Koopman operator provides a natural abstraction of halting, yielding a ‘Koopman halting problem’ that is recursively enumerable in general. For symbolic systems, such as those defined on Cantor space, this operator formulation captures reachability between clopen sets, while for equicontinuous systems we prove that the Koopman halting problem is decidable. Our framework demonstrates that absorbing (halting) states in coarse-grained finite automata correspond to Koopman eigenfunctions with eigenvalue one, while cycles in the transition graph impose spectral constraints associated with periodic dynamics. These results provide a unifying perspective on computation in symbolic and analog systems, showing how computational universality is reflected in operator spectra, invariant subspaces, and algebraic structures. Beyond symbolic dynamics, this operator-theoretic lens opens pathways to analyze the computational properties of a broader class of dynamical systems, including polynomial and analog models, and suggests that computational hardness may admit dynamical signatures in terms of Koopman spectral structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Experimental and theoretical investigation into the high pressure deflagration products of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

Diamond anvil cell (DAC) laser ignition experiments and reactive ab initio molecular dynamics (AIMD) simulations were performed on the high explosive (HE) LLM-105 to investigate its high pressure (HP) deflagration chemistry. Raman and optical spectroscopy measurements reveal LLM-105 reacts into an opaque carbonaceous product at 4–25 GPa. At pressures >~ 27 GPa, the reaction product consists of an amorphous optically transparent solid and nitrogen (N 2 ) in the solid phase. While not a one-to-one comparison due to the small time and length scales, the HP AIMD simulations show that some of the product is molecular N 2 , in qualitative agreement with experiment, while above 20 GPa most of the product consists of large amorphous C x H y N z O k clusters. Clustering is enhanced with pressure and reduces with temperature. In the experiments with initial sample pressure >~ 25 GPa, the pressure within the DAC decreases with minimal change in DAC cavity area. At initial sample pressures of 43.9 GPa, when quenched to 0 K, simulations predict a product experiencing a lower pressure consistent with the experimental measurement at lower load pressures. In conclusion, the results are important for understanding the HP deflagration chemistry of LLM-105.

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

Acceleration of Power System Dynamic Simulations Using a Deep Equilibrium Layer and Neural ODE Surrogate

The dominant paradigm for power system dynamic simulation is to build system-level simulations by combining physics-based models of individual components. The sheer size of the system along with the rapid integration of inverter-based resources exacerbates the computational burden of running time domain simulations. Here, in this paper, we propose a data-driven surrogate model based on implicit machine learningspecifically deep equilibrium layers and neural ordinary differential equationsto learn a reduced order model of a portion of the full underlying system. The data-driven surrogate achieves similar accuracy and reduction in simulation time compared to a physics-based surrogate, without the constraint of requiring detailed knowledge of the underlying dynamic models. This work also establishes key requirements needed to integrate the surrogate into existing simulation workflows; the proposed surrogate is initialized to a steady state operating point that matches the power flow solution by design.

Neural ordinary differential equations↗
DOE OSTI2025

California annual grass phenology and allometry influence ecosystem dynamics and fire regime in a vegetation demography model

Grass-dominated ecosystems cover wide areas of the land surface yet have received far less attention from the Earth System Model (ESM) community. This limits model projections of ecosystem dynamics in response to global change and coupled vegetation–climate dynamics. We used the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), a dynamic vegetation demography model, to determine ecosystem sensitivity to alternate, observed grass allometries and biophysical traits, and evaluated model performance in capturing California C 3 annual grasslands structure and fire regimes. Grass allometry, leaf physiology, plant phenology, and plant mortality all drove the seasonal variation in matter and energy exchange and fire dynamics in California annual grasslands. Allometry influenced grassland structure and function mainly through canopy architecture-mediated space and light competition instead of through carbon partitioning strategy. Regional variation in grassland annual burned area was driven by variation in ecosystem productivity. Our study advances the modeling of grassy ecosystems in ESMs by establishing the importance of grass allometry and plant phenology and mortality in driving C 3 annual grassland seasonal dynamics and fire regime. The calibrated annual grass allometry and biophysical traits presented can be applied in future studies to project climate–vegetation–fire feedbacks in annual grass-dominant ecosystems under global change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Nature of molybdenum carbide surfaces for catalytic hydrogen dissociation using machine-learned potentials: an ensemble-averaged perspective

Molybdenum carbides with an electronic structure similar to noble metals have gained attention as a promising low-cost catalyst for biomass valorization and the hydrogen evolution reaction. However, our fundamental understanding of the catalyst surface and how different phases of these catalysts behave at varying reaction conditions is limited to ground state density functional theory calculations as ab initio molecular dynamics (AIMD) is computationally prohibitive at relevant length and time scales. Here, in this work, we train a multi-atomic cluster expansion (MACE) machine-learned interatomic potentials (MLIP) to study hydrogen dissociation and dynamics over Mo, δ-MoC, α-Mo 2 C, and β-Mo 2 C surfaces at varying temperatures and hydrogen partial pressures. Our simulations identify unique and different molecular and atomic hydrogen adsorption sites on different surfaces that do not depend on the temperature. At low hydrogen pressures, the surface coverage is monolayer, which transitions to two-layer adsorption at higher pressures. We find that atomic hydrogen diffusion and recombinations are preferred over molybdenum atom hollow sites, while the diffusion over carbon-terminated facets was negligible, signifying particularly strong C–H interactions. In contrast, molecular hydrogen adsorption occurs mostly atop Mo or the bridging sites. At a comparable hydrogen loading, β-Mo 2 C (001) is the most active surface for hydrogen dissociation reaction. This work provides insights into the dynamic nature of the hydrogen dissociation chemistry and the diversity of hydrogen adsorption sites on molybdenum carbides.

08 HYDROGEN↗
DOE OSTI2025

Modeling COVID-19 Impacts on Prison Population

The goal of this project was to develop a unique, high-level, concept model that will (1) help examine the consequences COVID-19 has had on prison population and flow dynamics in Illinois and (2) explore (to the extent possible) how that will unfold in the short- and mid-term future. Argonne has unique expertise in dynamic systems modeling and highly specialized simulation capabilities, which will be helpful in developing such a model. The ultimate objective of this CRADA is to develop with the Illinois Public Sentencing Policy Advisory Council (SPAC) the capability to identify and understand factors and interactions associated with prison population and flow dynamics in Illinois. This understanding will then potentially lead to improvements in the allocation of various resources and services. This approach would also be useful to other states.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Integrated Metrics for County-Level Resilience Ranking Using Entropy and TOPSIS

In the face of atypical weather events, power infrastructure failures, and limited resources for resilience investment, energy decision-makers need data-driven metrics to allocate resilience investments and maximize the reduction of power outage impacts. For state-level planning, for instance, ranking the resilience of each county is key to ensuring effective distribution of resources. In such cases, resilience for each spatial unit is multifaceted and is captured by a set of indicators (i.e., metrics) that can be combined into an overall score that reduces the complexity of power outage dynamics to a single decision metric. However, weighting of these indicators is often addressed by simplifying assumptions (i.e., equal weights) or semi-subjective methods that rely on user-defined weights that can introduce biases (e.g., weighted average score). Within the disaster risk reduction and resilience engineering community, a recurring challenge in multicriteria decision-making is the objective weighting of indicators for composite indices. To address this issue, we have leveraged a Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS) combined with an entropy-based weighting approach to calculated the integrated scores. This method objectively determines the importance of each metric, better discerns between spatial units (i.e., counties), and offers a more reliable ranking of counties according to their relative resilience attributes. By improving methods for integrating resilience indicators, our approach helps planners and decision-makers prioritize resources more effectively for more efficient resilience investments.

Bhusal, Narayan [Oak Ridge National Laboratory (OR↗
DOE OSTI2025

Data-Driven Analysis of Multipactor Dynamics via Dynamic Mode Decomposition

Multipactor effect is a performance-limiting kinetic plasma effect that can occur in high-power microwave and radio frequency (RF) devices. Multipactor effect is of special concern in vacuum or near-vacuum conditions such as those in particle accelerators and spaceborne devices. In this work, we present a data-driven reduced-order model (ROM) based on dynamic mode decomposition (DMD) for modeling of multipactor effects. We study multipactor effects and the resulting nonlinear harmonic generation by processing high-fidelity data generated from electromagnetic particle-in-cell (EMPIC) simulations using the DMD algorithm. We also investigate time-delay embedding extensions of DMD with improved generalizability and accuracy for modeling the electron plasma current density behavior. Here, the results show that DMD provides valuable insights into multipactor phenomena by extracting relevant modal spatiotemporal patterns and frequencies. In addition, DMD offers the potential to time extrapolate EMPIC simulations at a minimal cost, thereby reducing overall simulation time.

43 PARTICLE ACCELERATORS↗