Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

On a Critical Acceleration Scale of Dark Matter in ΛCDM and Dynamical Dark Energy

Abstract Universal acceleration a 0 emerges in various empirical laws, yet its fundamental nature remains unclear. Using Illustris and Virgo N -body simulations, we focus on the velocity and acceleration fluctuations in collisionless dark matter involving long-range gravity. For comparison, in the kinetic theory of gases, molecules undergo random elastic collisions involving short-range interactions, where only velocity fluctuations are relevant. Hierarchical structure formation proceeds through the merging of smaller halos to form larger halos, which facilitates a continuous energy cascade from small to large halos at a constant rate ε u ≈ −10 −7 m 2 s −3 . Velocity fluctuations involve a critical velocity u c ∝ (1 + z ) −3/4 . Acceleration fluctuations involve a critical acceleration a c ∝ (1 + z ) 3/4 . Two critical quantities are related by the rate of energy cascade ε u ≈ − a c u c /[2(3 π ) 2 ], where factor 3 π is from the angle of incidence during merging. With critical velocity u c on the order of 300 km s −1 at z = 0, the critical acceleration is determined to be a c 0 ≡ a c ( z = 0) ≈ 10 −10 m s −2 , suggesting a c might explain the universal acceleration a 0 ≈ 10 −10 m s −2 in the empirical Tully–Fisher relation or modified Newtonian dynamics. The redshift evolution a c ∝ (1 + z ) 3/4 is in good agreement with Magneticum and EAGLE simulations and in reasonable agreement with limited observations. This suggests a larger a 0 at a higher redshift such that galaxies of fixed mass rotate faster at a higher redshift. Note that for dark energy (DE) density ρ DE 0 ≈ a c 0 2 / G = 1 0 − 10 J m −3 , we postulate an entropic origin of the DE from acceleration fluctuations of dark matter, analogous to the gas pressure from velocity fluctuations. This leads to a dynamical DE coupled to the structure evolution involving a relatively constant DE density followed by a slow weakening phase, suggesting possible deviations from the standard ΛCDM paradigm.

N-body simulations↗
DOE OSTI2026

Complex Dynamics in Argyrodite Solid-State Ion Conductors

Argyrodites are a compositionally diverse family of materials that exhibit remarkable ion transport properties. While the average crystal structures of argyrodites have been extensively studied, ion transport in these materials is governed by a confluence of dynamic processes spanning the cation, anion, and polyanionic sublattices. This Perspective synthesizes recent advances in understanding the role of dynamics in structural behavior and ion transport properties. We examine the compositional and structural motifs that govern order−disorder transitions within the argyrodite family and further explore how ion hopping is facilitated by lattice dynamics, from long-range phonons to local rotational dynamics of polyanionic species. Through the lens of dynamics spanning multiple time and length scales, we establish guiding principles that govern transport phenomena and highlight avenues of future study for the argyrodite family of ion conductors.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Study of the interaction between $Ξ$ baryons and light mesons via femtoscopy at the LHC

Meson-baryon systems with strangeness content provide a unique laboratory for investigating the strong interaction and testing theoretical models of hadron structure and dynamics. In this work, the measured correlation functions for oppositely charged $Ξ$-K and $Ξ$ − 𝜋 pairs obtained in high-multiplicity pp collisions at $\sqrt{𝑠}$ = 13 TeV at the LHC are presented. For the first time, high-precision data on the $Ξ$-K interaction are delivered at small relative momenta. The scattering lengths, extracted via the Lednický–Lyuboshits expression of the pair wavefunction, indicate a repulsive and a shallow attractive strong interaction for the $Ξ$-K and $Ξ$ − 𝜋 systems, respectively. The $Ξ$(1620) and $Ξ$(1690) states are observed in the $Ξ$ − 𝜋 correlation function and their properties, mass and width, are determined. These measurements are in agreement with other available results. Such high-precision data can help refine the understanding of these resonant states, provide stronger constraints for chirally motivated potentials, and address the key challenge of describing the coupled-channel dynamics that may give rise to molecular configurations .

Femtoscopy↗
DOE OSTI2026

Thermal Reservoir Networks for Modularly Expandable Thermal Microgrids

The Department of Defense (DoD) faces the substantial challenge of cost-effectively retrofitting one to two installations per month, each comprising approximately 1,000 buildings, to improve resilience, reduce energy consumption, and enhance energy supply security. Achieving these objectives requires optimal system selection and effective risk mitigation during system integration. To address this need, we introduce Platform-Based Design (PBD), a structured, hierarchical methodology adapted from other industrial sectors to the domain of energy system retrofits. We demonstrate the effectiveness of PBD through a techno-economic feasibility study comparing geothermal-coupled thermal energy networks (TENs) with conventional energy systems for heating, cooling, and powering 17 buildings at Joint Base Andrews (JBA) in Maryland. Our analysis illustrates that the PBD approach enables rigorous, data-driven, sequential decision making, resulting in a family of Pareto-optimal systems, among which the TEN emerged as the most promising solution. The selected TEN design integrates geothermal borefields, heat recovery heat pumps, photovoltaic (PV) arrays, and battery storage. Compared to the baseline system – gas heating combined with air-source chillers – the proposed TEN reduces annual imported energy by 74% and peak electricity demand by 45%, achieves a levelized cost of energy of $\$0.210$/kWh, and substantially enhances resilience. Life-cycle costs increase by approximately 6%, and initial investment costs are about 2.5 times higher than the baseline. However, if central plant infrastructure, district loops, and utility-scale PV and battery systems are privately funded and operated, the initial investment would fall below the baseline system cost. Critical to achieving these significant performance improvements were detailed nonlinear dynamic simulations coupling geothermal heat transfer, energy system operation, and realistic feedback control logic. These simulations identified essential design modifications and control strategy refinements that substantially reduced energy use, peak demand, and compressor shortcycling, thereby improving durability and reliability—issues that would have been significantly more expensive to resolve during operation. Additionally, the verification step highlighted sensitivities to key design parameters that could reduce initial investment by approximately $\$2$ million and reduce annual life-cycle costs more than $\$300,000$. We recommend adopting the PBD methodology for future feasibility studies and TEN pilot projects to gain valuable operational experience. Furthermore, we recommend that DoD invest in transferring and scaling the PBD methodology to other installations. This entails developing standardized computational frameworks and component libraries as well as training industry in conducting PBD. Such investments would enable rapid, robust, reliable, and cost-effective retrofits, supporting DoD’s ambitious energy system modernization goals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Coarse-graining Hamiltonian systems using WSINDy

Abstract Weak form equation learning and surrogate modeling has proven to be computationally efficient and robust to measurement noise in a wide range of applications including ODE, PDE, and SDE discovery, as well as in coarse-graining applications, such as homogenization and mean-field descriptions of interacting particle systems. In this work we extend this coarse-graining capability to the setting of Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. A smooth $$\varepsilon$$ ε -dependent Hamiltonian vector field $$X_\varepsilon$$ X ε possesses an approximate symmetry if the limiting vector field $$X_0=\lim _{\varepsilon \rightarrow 0}X_\varepsilon$$ X 0 = lim ε → 0 X ε possesses an exact symmetry. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large perturbations imparted in the $$\varepsilon >0$$ ε > 0 regime, while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computationally efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. In this way, we argue that weak-form equation learning is particularly well-suited for Hamiltonian coarse-graining. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. While our main contribution is computational, we also provide a contribution to the literature on averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. This provides theoretical justification for our approach as WSINDy’s computations occur at the level of Hamiltonian vector fields. We illustrate the efficacy of our proposed method using physically relevant examples, including coupled oscillator dynamics, the Hénon–Heiles system for stellar motion within a galaxy, and the dynamics of charged particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Charged- and Multi-Exciton Dynamics in Colloidal Quantum Dot Molecules

Multicarrier states in quantum dots are confined to small volumes, resulting in increased nonradiative Auger recombination rates with implications for different optoelectronic applications. Recently, the fusion of two core-shell quantum dots into a dimer has provided a new physical landscape for multiexciton states, since the excitons may share a core (intradot, localized) or occupy different cores (interdot, segregated). Here, in this study, we employ transient absorption spectroscopy to investigate the multiexciton dynamics in coupled quantum dot dimers. We observe that multiexciton populations in the dimers live significantly longer in comparison to the parent monomers, in contrast to the single exciton regime. A kinetic model that accounts for the statistical differences between monomers and dimers reveals that, while intradot multiexcitons show Auger rates similar to the monomers, interdot states have reduced Auger rates. These results pave the way for the rational design of new quantum dot molecules with tailored multiexciton properties.

Auger recombination↗
DOE OSTI2026

Topology-Informed Design Rules for Deconstructable Thermoset Copolymer Networks

Existing models of thermoset deconstruction facilitated by incorporating cleavable comonomers rely on a mean-field reverse gel point paradigm, which predicts network dissolution once cleavable bonds reach a critical stoichiometric threshold, but does not account for where those bonds reside within the network architecture. Using reactive coarse-grained molecular dynamics simulations coupled with graph-theoretic analysis, we extend this stoichiometric picture to show that deconstructability is governed by the curing-imprinted network topology rather than stoichiometry alone. This topological organization is hierarchical: at the local scale, the elastic effectiveness of cross-link junctions determines which cross-links constitute the load-bearing scaffold; at the mesoscale, the cross-linking rate kinetically templates that scaffold into topologically modular communities─densely cross-linked clusters connected by sparse bridging strands that sustain network connectivity. Using betweenness centrality to identify nodes that disproportionately lie on intercommunity shortest paths, we demonstrate that effective deconstruction of the network into macromolecular fragments requires cleavable comonomers to intercept these high-centrality bridging strands. We further find that under uniform, disassortative comonomer incorporation, this topological requirement provides a mechanistic basis for extending the reverse gel point to incorporate network topology. We also show that modularity imposes a fundamental limit on fragment uniformity that persists even when the centrality requirement is met. Finally, we demonstrate that chain stiffness provides a nearly independent lever to suppress mechanically redundant cross-links and raise the glass transition temperature without significantly altering the deconstruction outcome. Together, these findings reframe the thermoset design space around network topology and provide actionable guidelines for engineering thermoset copolymers with predictable deconstructability and targeted thermomechanical performance.

coarse-grained molecular dynamics↗
DOE OSTI2025

Polarization Switching on the Open Surfaces of the Wurtzite Ferroelectric Nitrides: Ferroelectric Subsystems and Electrochemical Reactivity

Binary ferroelectric nitrides are promising materials for information technologies and power electronics. However, polarization switching in these materials is highly unusual. From the structural perspective, polarization reversal is associated with the change of the effective polarity at the surfaces and interfaces from N‐to‐M terminated, suggesting strong coupling between ferroelectric and chemical phenomena. Phenomenologically, macroscopic studies demonstrate the presence of complex time dependent phenomena including wake‐up. Here, in this study, the polarization switching using the multidimensional high‐resolution piezoresponse force microscopy (PFM) and spectroscopy is explored, detecting both the evolution of induced ferroelectric domain, electromechanical response, and surface deformation during first‐order reversal curve measurements. The presence of two weakly coupled ferroelectric subsystems are demonstrated and the bias‐induced electrochemical reactivity. The observed behaviors are very similar to the recent studies of other wurtzite system but additionally include electrochemical reactivity, suggesting the universality of these behaviors for the wurtzite binary ferroelectrics. These studies suggest potential of high‐resolution multimodal PFM spectroscopies to resolve complex coupled polarization dynamics in materials. Furthermore, these PFM based studies are fully consistent with the recent electron microscopy observations of the shark‐teeth like ferroelectric domains in nitrides. Hence, it is believed that these studies establish the universal phenomenological picture of polarization switching in binary wurtzite.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Engineering design of a kW-scale continuous reactor-heat exchanger for high temperature discharge of particle-based thermochemical energy storage

This study investigates the theoretical design parameters and thermal performance of a kW-scale continuous oxidation reactor for high temperature (~1000 °C) thermochemical energy storage (TCES) applications. The concept comprises a counter-current particle-based system that includes a reaction zone with a heat exchanger to extract the heat produced from the oxidation reaction. Both above and below the hot reactive volume are sensible heat recuperation zones to enable the feed and removal of particles and oxidizing gas near ambient temperature during steady state operation. Two operation types for the reaction zone are studied, a fluidized bed reactor (FBR) and a moving bed reactor (MBR). The results of the parametric analysis suggest that the MBR requires a smaller volume per kW of heat produced, achieving power densities in excess of 2500 kW/m3 compared to ~ 900 kW/m 3 in the FBR. Additionally, the MBR achieves between 0.71 and 0.99 oxidation conversions compared to between 0.23 and 0.38 conversions in the FBR with the same volumes and flowrates. However, the FBR has the potential to maintain a uniform reactor temperature which can produce heat transfer fluid (HTF) outlet temperatures as high as the reactor temperature, i.e., ~1000 °C, whereas the MBR produces variable reactor temperatures that can create overheating zones and low HTF outlet temperatures (< 800 °C) depending on the operating conditions selected. Future work should aim at understanding the coupled fluid dynamics, heat and mass transfer, and thermochemical reaction for any given combination of reactor volume and contacting patterns. Here, these studies should be complemented by experimental work on particle-gas TCES reactors.

25 ENERGY STORAGE↗
DOE OSTI2025

Visualizing the Three-Dimensional Arrangement of Hydrogen Atoms in Organic Molecules by Coulomb Explosion Imaging

Structure-sensitive methods based on femtosecond light or electron pulses are now making it possible to measure how molecular structures change during light-induced processes. Despite significant progress, high-fidelity imaging of nuclear positions remains a challenge even for relatively small molecular systems and, notably, regarding the positions of hydrogen atoms. As demonstrated in recent work, X-ray-induced Coulomb explosion imaging (CEI) may overcome this obstacle, as its sensitivity does not depend on the mass of the imaged atoms. The photoinduced ring opening of the heterocyclic molecule 2(5 H )-thiophenone has attracted recent interest. Here, in this work, we show that CEI offers a powerful route to imaging the peripheral H atoms in this molecule and thus, more generally, to tracking detailed nuclear motions (e.g., isomerizations) in organic molecules on ultrafast time scales. Specifically, we record momentum-space Coulomb explosion images that report on the three-dimensional positioning of all nuclei within the molecule, for instance, distinguishing H atoms in C–H bonds that lie within or are directed out of the plane defined by the heavy atoms. The prospect of imaging peripheral H atoms to probe photochemical dynamics is explored by coupling ab initio molecular dynamics with classical Coulomb explosion simulations, thereby differentiating potential photoproduct isomers, including those whose structures primarily differ in the position of the hydrogens.

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

Active species in chloroaluminate ionic liquids catalyzing low-temperature polyolefin deconstruction

Abstract Chloroaluminate ionic liquids selectively transform (waste) polyolefins into gasoline-range alkanes through tandem cracking-alkylation at temperatures below 100 °C. Further improvement of this process necessitates a deep understanding of the nature of the catalytically active species and the correlated performance in the catalyzing critical reactions for the tandem polyolefin deconstruction with isoalkanes at low temperatures. Here, we address this requirement by determining the nuclearity of the chloroaluminate ions and their interactions with reaction intermediates, combining in situ 27 Al magic-angle spinning nuclear magnetic resonance spectroscopy, in situ Raman spectroscopy, Al K-edge X-ray absorption near edge structure spectroscopy, and catalytic activity measurement. Cracking and alkylation are facilitated by carbenium ions initiated by AlCl 3 - tert -butyl chloride (TBC) adducts, which are formed by the dissociation of Al 2 Cl 7 − in the presence of TBC. The carbenium ions activate the alkane polymer strands and advance the alkylation cycle through multiple hydride transfer reactions. In situ 1 H NMR and operando infrared spectroscopy demonstrate that the cracking and alkylation processes occur synchronously; alkenes formed during cracking are rapidly incorporated into the carbenium ion-mediated alkylation cycle. The conclusions are further supported by ab initio molecular dynamics simulations coupled with an enhanced sampling method, and model experiments using n-hexadecane as a feed.

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

Nanoscale strain wave generation by a piezoelectric grating from polar vortices

Nanostructures formed by spontaneously broken symmetry have provided new ways to manipulate quantum states. Specifically, topological structures with periodic spatial ordering, such as polar vortices and skyrmions, can be ideal hosts for creating engineered responses in both spatial and frequency domains. So far, however, only a few examples of such hierarchical engineering have been reported in the literature. Here we demonstrate that the spatially modulated piezoelectric response of a polar vortex structure can create strain waves with a characteristic nanoscale wavefront. Using time-resolved pump–probe resonant X-ray scattering and diffraction measurements, coupled with dynamical phase-field simulations, we show that the piezoelectric modulation of the spontaneously formed polar vortex crystal functions as an acoustic diffraction grating. This system converts incoming laterally uniform strain waves into outgoing waves with a characteristic sub-terahertz frequency, driven by an intrinsic excitation of the polar vortex crystal. Moreover, our phase-field simulations suggest that the dynamic mechanical displacements exhibiting vortex textures are generated from both space- and time-varying piezoelectric responses. Our findings illustrate a new method for generating nanoscale strain waves with unique spatial textures by tuning the hierarchical order of polar topologies to engineer new collective modes, allowing for a wide range of control through the topological lattice.

ferroelectrics↗
DOE OSTI2025

Dynamically downscaled seasonal heat wave projections in the CONUS

Heat waves are a well-documented hazard that are projected to increase in intensity, duration, and frequency with climate change. Regions of the US experience widely varying temperatures; for example, 35 °C is extremely hot for spring in the Northeast but not for summer in the Southeast. It is important to evaluate projections within a regional context and at a high enough resolution to understand the risks to populations. We identify heat waves across the Conterminous US (CONUS) under SSP5–8.5 from 2020 to 2059 with an ensemble of dynamically downscaled Coupled Model Intercomparison Project Phase 6 (CMIP6) model outputs. We demonstrate that there are regional differences caused by seasonal and local drivers of persistent hot temperatures. Summer heat waves are increasing in intensity and duration faster than winter heat waves because of the atmospheric conditions that promote these events. Our analysis emphasizes the value of fine-resolution modeling for projecting future climate risks.

Rubin, Hannah [University of Tennessee, Knoxville ↗
DOE OSTI2026

Thermodynamically consistent incorporation of the Langmuir adsorption model into compressible fluctuating hydrodynamics

For a gas–solid interfacial system where chemical species undergo reversible adsorption, we develop a mesoscopic stochastic modeling method that simulates both gas-phase hydrodynamics and surface coverage dynamics by coupling the Langmuir adsorption model with compressible fluctuating hydrodynamics. To this end, we derive a thermodynamically consistent mass–energy update scheme that accounts for how the mass and energy variables in the gas and surface subsystems should be updated according to the changes in the number of molecules of each species in each subsystem due to adsorption and desorption events. By performing a stochastic analysis for the ideal Langmuir model and the full hydrodynamic system, we analytically confirm that our mass–energy update scheme captures thermodynamic equilibrium predicted by equilibrium statistical mechanics. We find that an internal energy correction term is needed, which is attributed to the difference in the mean kinetic energy of gas molecules colliding with the surface from that computed from the Maxwell–Boltzmann distribution. By performing an equilibrium simulation study for an ideal gas mixture of CO and Ar, with CO undergoing reversible adsorption, we validate our overall simulation method and implementation.

Adsorption↗
DOE OSTI2025

Analyzer-based X-ray phase contrast imaging using the forward-diffracted o-beam in a few-millimetre-thick Bragg-case asymmetrically cut crystal

An analyzer-based X-ray phase contrast imaging experiment employing the forward-diffracted o-beam in a few-millimetre-thick Bragg-case asymmetrically cut analyzer crystal has been implemented and tested in a geometry similar to that used in conventional radiography. The high angular stability requirements were overcome using specially designed crystals and a closed-loop feedback intensity system, monitoring the intensity of the diffracted h-beam. Analyzer-based images (ABis) taken at different angular positions on the forward-diffracted o-beam rocking curve show different contrasts, as expected. However, the ABis did not show sharp borders. This was closely investigated by X-ray topography and Δ d / d mapping and was attributed to the stresses caused by the analyzer crystal's own weight in the bulk. Further investigation of the crystal design using finite element analysis coupled with dynamical theory of X-ray diffraction is envisaged.

Honnicke, Marcelo Goncalves (ORCID:000000033276432↗
DOE OSTI2025

Widening of Wind Stress Anomalies Amplifies ENSO in a Warming Climate

Abstract Climate change simulations generally indicate the strengthening of El Niño–Southern Oscillation (ENSO) sea surface temperature (SST) variability through the twenty-first century, yet a robust physical mechanism explaining this change across different models is still lacking, and the projections for ENSO amplitude exhibit a large spread. Most commonly, changes in the background state of the tropical Pacific are invoked to explain these changes of ENSO. Here, we show that changes in the structure of wind stress anomalies associated with ENSO are potentially as important as these background state changes. Specifically, changes in the magnitude, meridional width, and zonal structure of wind stress anomalies can explain approximately 53% of the intermodel variance in the projected change of ENSO magnitude through the twenty-first century as well as 43% in ENSO periodicity changes. Among these changes in the wind structure, the meridional widening of wind anomalies plays the most important role. To demonstrate that these changes are indeed critical, we develop a hybrid model of ENSO based on the Community Earth System Model, version 2, which incorporates a dynamical ocean coupled to a simplified statistical atmosphere within the tropical Pacific. In the absence of external forcing and corresponding mean-state changes, the imposed changes in wind stress anomalies in this hybrid model result in an increase of ENSO amplitudes of nearly 10% along the equator. Our results are also theoretically supported by a recharge-oscillator model that incorporates the meridional wind structure. Thus, changes in the structure of wind stress anomalies, together with changes in the mean state, likely play a critical role in the projected strengthening of ENSO.

Stuivenvolt-Allen, Jacob [Yale University, New Hav↗
DOE OSTI2024

Development of a MOOSE-based crystal plasticity model with irradiation defect evolution for irradiation creep in 316

Irradiation creep and irradiation swelling are two of the lifetime limiting factors for structural materials in nuclear reactors. These mechanical effects are driven by irradiation defect evolution and the interaction of those defects with dislocations in the microstructure. We present here a coupled cluster dynamics – crystal plasticity approach to model irradiation swelling and creep behavior in 316 SS. The time-dependent evolution of irradiation defects is calculated with a cluster dynamics approach and passed to the crystal plasticity model to compute the dislocation evolution. We show the impact of the irradiation defect evolution on the stress state in the material, which drives inelastic deformation through dislocation motion. The inelastic deformation in the 316 SS is dependent on the dose rate, where the inelastic deformation driven by the early-stage irradiation defect evolution determines the mechanical behavior of the 316 SS.

316 Stainless Steels↗