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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport

In a computational fluid model of the atmosphere, the advective transport of trace species, or tracers, can be computationally expensive. For efficiency, models often use semi-Lagrangian advection methods. High-order interpolation semi-Lagrangian (ISL) methods, in particular, can be extremely efficient, if the problem of property preservation specific to them can be addressed. Atmosphere models often use geometrically and logically nonuniform grids for efficiency and, as a result, element-based discretizations. Such grids and discretizations make stability a particular problem for ISL methods. Generally, high-order, element-based ISL methods that use the natural polynomial interpolant associated with a nodal finite-element discretization are unstable. Here, we derive new bases having order of accuracy up to nine, with positive nodal weights, that stabilize the element-based ISL method. We use these bases to construct the linear advection operator in the property-preserving Interpolation Semi-Lagrangian Element-based Transport (Islet) method. Then we discuss key software implementation details. Finally, we show performance results for the Energy Exascale Earth System Model's atmosphere dynamical core, comparing the original and new transport methods. These simulations used up to 27,600 Graphical Processing Units (GPU) on the Oak Ridge Leadership Computing Facility's Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Numerical simulation of vortex-induced vibration response of a single IEA 10-MW wind turbine blade

Three-dimensional simulation of vortex-induced vibration (VIV) of a single International Energy Agency (IEA) 10-MW reference wind turbine blade with a length of 97.325 m is performed using the ExaWind stack, an open-source suite of codes. This study aims to illustrate the spanwise VIV response characteristics and cross-validate the results with an existing commercial framework. Five near-body meshes and three time steps are selected for the convergence study. To improve computational efficiency, several VIV triggering methods are also compared to shorten the VIV development period. The ExaWind-based VIV simulation strategy for a single IEA 10-MW blade is determined. First, the modal shape is validated against published results. Then, spanwise VIV responses of four blade configurations under a fixed and varied incoming flow velocity are analyzed. Results show that the VIV response is dominated by the first edgewise (second overall) mode. Little first-mode contributions appear near the second-mode node, producing a pi phase jump, and a higher harmonics response occurs near the blade root. Rotational degrees of freedom are minor compared with translational motion. The response versus reduced velocity is analyzed, showing a two-branch behavior similar to that of VIV for a bluff cylinder. Across all tested cases, the dominant frequency remains locked to the natural frequency of the second mode with no observed desynchronization. A mild deviation is observed for the case of 90-degree pitch and 310-degree azimuth rotation near a reduced velocity of 6, which will be examined with additional cases in future work. These findings indicate that severe VIV responses can arise under specific configurations and flow conditions, thereby increasing the potential for VIV fatigue damage and requiring greater attention during operation.

17 WIND ENERGY↗
DOE OSTI2025

Single particle impact explorations on the role of powder heat treatments in cold spray: Effects of strength and oxide structure

In cold spray, the optimization of the powder feedstock for improving the final properties of the coating and for improving the process itself are two different objectives. This paper explores the latter, studying the relationship between powder pre-heat treatment and the kinetic bonding between the individual particles and the substrate. Using a laser-induced particle impact test (LIPIT), we investigate the effect of different heat treatments on the critical adhesion velocity of Al 6061. Although heat treatments change the strength of the powders as well as the thickness of the native oxide layer each by as much as ∼50 %, the critical adhesion velocity remains constant to within only ∼10 %. A deeper analysis reveals that these two factors change in somewhat complementary ways for the heat treatments explored here, leading to counteracting effects. By quantitatively evaluating these two microstructural factors in tandem, we find that they are roughly equally important in determining the critical adhesion velocity. These observations thus provide insight for the development of new heat treatment regimens that might lower the adhesion velocity.

Al 6061↗
DOE OSTI2026

Regional specialization in prefrontal cortex manifests in the reliability of task progression codes

The brain has the remarkable ability to guide the performance of complex tasks. Distinct prefrontal cortical areas make specific contributions to this ability, with the orbitofrontal cortex (OFC) critical for processing information related to trial outcomes and the dorsomedial prefrontal cortex (dmPFC) critical for sustained effort and selecting the right action at the right time. Yet, in both areas, neural activity represents both outcome- and action-related quantities. How similar neural representations support different functions remains unclear. Here, we compared OFC and dmPFC activity in rats performing a spatial alternation task. We show that, in contrast to other task-related variables, task progression is represented in both areas, but with distinct patterns of across-trial reliability that match each area’s previously documented functional specialization. Our results indicate that the engagement of reliable, task-phase-specific activity patterns differs across prefrontal regions in a manner well suited to engage different computations at different times.

Biological and medical sciences↗
DOE OSTI2025

The impact of $γN$ and $γ^∗N$ interactions on our understanding of nucleon excitations

We review recent progress in our understanding of the nucleon excitation spectrum. Thanks to dedicated efforts at facilities such as ELSA, MAMI and Jefferson Lab, several new nucleon resonances have been discovered, and evidence for previously elusive states has been significantly improved. Numerous decay channels have been observed for the first time, and resonance properties are being extracted from these data by several groups through coupled-channel analyses of varying complexity. Electroproduction experiments have provided further insights into the internal structure of light baryon resonances – for example, the long-debated Roper resonance N (1440) is observed as a three-quark state with a significant meson-cloud component. While the non-relativistic quark model remains a valuable tool for organizing the spectrum of nucleon and Δ resonances, a variety of theoretical frameworks have emerged to offer deeper understanding, including phenomenological quark models, holographic QCD, functional methods, effective field theories, and lattice QCD. We examine the interplay between these approaches, highlight their respective strengths and explore how they complement each other in shaping our knowledge of light baryon resonances. We address several open questions in baryon spectroscopy, including the nature of the enigmatic Λ (1405), ongoing searches for exotic states such as hybrid baryons and pentaquarks, and the dichotomy between microscopic descriptions of baryons in terms of quarks and gluons versus effective hadronic descriptions based on meson-baryon dynamics.

Baryon resonances↗
DOE OSTI2025

Improving Bond Dissociations of Reactive Machine Learning Potentials through Physics-Constrained Data Augmentation

In the field of computational chemistry, predicting bond dissociation energies (BDEs) presents well-known challenges, particularly due to the multireference character of reactive systems. Many chemical reactions involve configurations where single-reference methods fall short, as the electronic structure can significantly change during bond breaking. As generating training data for partially broken bonds is a challenging task, even state-of-the-art reactive machine learning interatomic potentials (MLIPs) often fail to predict reliable BDEs and smooth dissociation curves. By contrast, simple and inexpensive physics-based models, such as the well-established Morse potential, do not suffer from any such limitations. This work leverages the Morse potential to improve reactive MLIPs by augmenting the training data set with inexpensive Morse data along the dissociation pathways. Further, this physics-constrained data augmentation (PCDA) approach results in MLIPs with smooth bond dissociation curves as well as near coupled-cluster level BDEs, all without requiring any expensive multireference quantum mechanical calculations. A case study for methane combustion demonstrates how the PCDA approach can improve an existing reactive MLIP, namely, ANI-1xnr. In conclusion, not only are the BDEs and bond dissociation curves for all radicals and molecules significantly improved compared to ANI-1xnr but the PCDA-trained MLIP retains the reliability of ANI-1xnr when performing reactive molecular dynamics simulations.

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

Sliceable, Moldable, and Highly Conductive Electrolytes for All-Solid-State Batteries

All-solid-state batteries (ASSBs) require solid electrolytes with high ionic conductivity, stability, and deformability for optimal energy and power density. Here, we developed lithium-deficient lithium yttrium bromide (LYB) solid electrolytes, Li 3–x YBr 6–x (0 ≤ x ≤ 0.50), using a comelting method with controlled lithium deficiency. These electrolytes exhibit favorable mechanical properties such as high moldability and sliceability. The Li 2.65 YBr 5.65 composition has an ionic conductivity of 4.49 mS cm –1 at 25 °C and an activation energy of 0.28 eV. Compared to Li 3 YBr 6 , Li 2.65 YBr 5.65 demonstrates improved rate performance and cycling stability in ASSBs. High-resolution X-ray diffraction confirms the formation of the LYB phase with a C2/m space group. Structural analysis reveals increased cation disorder and larger polyhedral volumes for x > 0 in Li 3–x YBr 6–x , contributing to reduced Li + migration energy barriers. Bond valence site energy calculations and molecular dynamics simulations reveal enhanced 3D lithium-ion transport. NMR spectroscopy further highlights increased Li + dynamics and impurity elimination.

Poudel, Tej P. [Florida State Univ., Tallahassee, ↗
DOE OSTI2025

Measuring error rates of mid-circuit measurements

High-fidelity mid-circuit measurements, which read out the state of specific qubits in a multiqubit processor without destroying them or disrupting their neighbors, are a critical component for useful quantum computing. They enable fault-tolerant quantum error correction, dynamic circuits, and other paths to solving classically intractable problems. But there are few methods to assess their performance comprehensively. In this work, we address this gap by introducing the first randomized benchmarking protocol that measures the rate at which mid-circuit measurements induce errors in many-qubit circuits. Using this protocol, we detect and eliminate previously undetected measurement-induced crosstalk in a 20-qubit trapped-ion quantum computer. Then, we use the same protocol to measure the rate of measurement-induced crosstalk error on a 27-qubit IBM Q processor, and quantify how much of that error is eliminated by dynamical decoupling.

Hothem, Daniel [Sandia National Laboratories (SNL-↗
DOE OSTI2026

Neural network potentials with effective charge separation for non-equilibrium dynamics of ionic solids: a ZnO case study

Developing neural network potentials (NNPs) accurate under non-equilibrium dynamics is challenging, as such systems require extensive sampling beyond equilibrium phases. Here we construct high-fidelity NNPs for zinc oxide (ZnO), a polymorphic ionic solid, using density functional theory (DFT) reference data. To efficiently capture transitional configurations, we combine enhanced-sampling molecular dynamics with empirical potentials, data distillation, and pretraining on short-range atomic energies (A-Train), followed by transfer learning with DFT-relabeled datasets. This hierarchical approach improves transferability across polymorphs and stress states. We further introduce effective charge separation, treating long-range Coulombic terms analytically while short-range residual interactions are learned by the NNP. The optimal effective charges fall in the range 0.5–1.0 q e , consistent with dielectric-screened values derived from formal charges but distinct from Bader estimates. Motivated by this observation, we propose a simple data-driven protocol in which effective charges are optimized by comparing DFT reference energies with explicit Coulomb calculations, without additional NNP training. This strategy improves accuracy and transferability in DFT-level predictions of energies, forces, and stress. Together, these results provide a practical charge-selection framework for robust NNP development in ionic solids, enabling reliable simulation of polymorphic phase transformations and non-equilibrium dynamics.

Chemistry↗
DOE OSTI2025

Hyperfine-to-rotational energy transfer in ultracold atom–molecule collisions of Rb and KRb

Energy transfer between different mechanical degrees of freedom in atom–molecule collisions has been studied and largely understood. However, systems involving spins remain less explored. Here, in this study, we directly observed energy transfer from atomic hyperfine to molecular rotation in the 87 Rb ($| {F}_{a},{M}_{{F}_{a}}\rangle =| 2,2\rangle$) + 40 K 87 Rb (X 1 Σ + , rotational state N = 0) ⟶ Rb ($| 1,1\rangle$) + KRb ( N = 0, 1, 2) collision with state-to-state precision. We also performed quantum scattering calculations that rigorously included the coupling between spin and rotational degrees of freedom at short range under the assumption of rigid-rotor KRb monomers moving along a single potential energy surface. The calculated product rotational state distribution deviates from the observations even after extensive tuning of the atom–molecule potential energy surface. In addition, our ab initio calculations indicate that spin–rotation coupling is enhanced close to a conical intersection that is energetically accessible at short range. This, together with the deviation, suggests that vibrational degrees of freedom and conical intersections play an important part in the coupling. Our observations confirm that spin is coupled to mechanical rotation at short range and establish a benchmark for future theoretical studies.

chemical physics↗
DOE OSTI2025

Complex water networks visualized by cryogenic electron microscopy of RNA

The stability and function of biomolecules are directly influenced by their myriad interactions with water. Here we investigated water through cryogenic electron microscopy (cryo-EM) on a highly solvated molecule: the Tetrahymena ribozyme. By using segmentation-guided water and ion modelling (SWIM), an approach combining resolvability and chemical parameters, we automatically modelled and cross-validated water molecules and Mg 2+ ions in the ribozyme core, revealing the extensive involvement of water in mediating RNA non-canonical interactions. Unexpectedly, in regions where SWIM does not model ordered water, we observed highly similar densities in both cryo-EM maps. In many of these regions, the cryo-EM densities superimpose with complex water networks predicted by molecular dynamics, supporting their assignment as water and suggesting a biophysical explanation for their elusiveness to conventional atomic coordinate modelling. Our study demonstrates an approach to unveil both rigid and flexible waters that surround biomolecules through cryo-EM map densities, statistical and chemical metrics, and molecular dynamics simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Jahn–Teller distortion controls electron transfer in photoexcited Cu( I ) donor–acceptor systems

The Jahn–Teller distortion (JTD) is a defining structural response to electronic excitation in Cu( I )-based transition metal complexes, yet its role in photoinduced electron transfer (PET) remains largely unexplored. Here, we demonstrate that the JTD governs charge separation (CS) through a vibronically controlled conical intersection in heteroleptic Cu( I ) bisphenanthroline–naphthalene diimide (CuHETPHEN-NDI) donor–acceptor dyads. Using 20-fs broadband transient absorption spectroscopy combined with coherent vibrational wavepacket (CVWP) analysis and quantum chemical calculations, we directly track nuclear motions that steer the system from the metal-to-ligand charge-transfer ( 1 MLCT) state to the CS state. Steric bulkiness of the pendant groups at the 2,9-positions of the phenanthroline ligand systematically slows both JTD and CS. Short-time Fourier transformation and Fourier filtering analyses identify two key vibrational signatures: a low-frequency breathing mode (∼100 cm −1 ) via a bond distance change between Cu and ligated Ns (Cu–N) that modulates the NDI anion absorption and acts as a vibronic coupling coordinate, and a higher-frequency mode (∼313 cm −1 ) that evolves along the PET trajectory. Normal mode analysis and potential energy surface calculations show that the JTD brings the 1 MLCT and CS states into degeneracy, while the Cu–N breathing motion dynamically modulates donor–acceptor electronic coupling to enable ultrafast nonadiabatic electron transfer. Steric hindrance exerted by the groups at the 2,9 positions of the phenanthroline ligands suppresses this vibronic coupling, leading to faster CVWP decoherence for the 313 cm −1 mode and slower CS. These findings unravel JTD-controlled vibronic coupling at conical intersection as a governing factor for CS and provide insight into designing Cu-based photosensitizers by harnessing structural dynamics to control PET.

Kim, Pyosang [Argonne National Laboratory (ANL), A↗
DOE OSTI2025

Turbulent-like flows in quasi two-dimensional dense suspensions of motile colloids

Dense bacterial suspensions exhibit turbulent-like flows at low Reynolds numbers, driven by the activity of the microswimmers. In this study, we develop a model system to examine these dynamics using motile colloids that mimic bacterial locomotion. The colloids are powered by the Quincke instability, which causes them to spontaneously roll in a random-walk pattern when exposed to a square-wave electric field. We experimentally investigate the flow dynamics in dense suspensions of these Quincke random walkers under quasi two-dimensional conditions, where the particle size is comparable to the gap between the electrodes. The results suggest a scaling regime in the energy spectrum ∼k −4 at high wavenumbers, observed consistently across activity levels and particle concentrations. We observe that velocity time correlations decay within a single period of the square-wave field, yet an anti-correlation appears between successive field applications, indicative of a dynamic structural memory of the ensemble.

Luo, Rui [Northwestern Univ., Evanston, IL (United↗
DOE OSTI2025

Genetic control of morphological transitions in a coacervating protein template

Nature routinely exploits liquid–liquid phase separation (LLPS) of proteins to control the assembly and mineralization of hybrid materials. Here, we show that fusion of the Car9 silica-binding peptide to an elastin-like polypeptide (ELP) yields temperature- and sequence-programmable soft matter templates for the synthesis of silicified architectures ranging in size from nanometers to micrometers. Specifically, we demonstrate unprecedented control over the diameter of silica nanoparticles (SiNP) in the 30–60 nm range with 4 nm precision, show that a single arginine residue (R4) in the Car9 sequence underpins the transition from micelles to proteinosomes, and find that substitutions in other basic residues modulate electrostatic repulsion and solvation to enable access to kinetically trapped species. These structures, which include interconnected micelles, small (∼200 nm) and large (>5 µm) vesicles, are readily visualized by SEM imaging following silicification. Molecular dynamics (MD) simulations and AlphaFold predictions reveal that mutations in positively charged residues alter interfacial packing, hydration, and conformational freedom of the silica-binding segments. Overall, our results establish sequence and thermal energy as synergistic levers for morphological control across length scales using solid-binding ELPs and establish mineralization as a powerful tool to visualize the structure of dynamic soft matter assemblies.

hierarchy↗
DOE OSTI2025

Kelvin–Helmholtz instability under stabilizing parallel magnetic field in nonhomogeneous compressible MHD flows

We study the Kelvin–Helmholtz instability (KHI) for the general case of a compressible, nonhomogeneous, magnetized plasma flow. The study is limited to a vortex sheet interface with an imposed parallel magnetic field. We introduce a new formalism based on a convective Mach number M c , a convective Alfvénic Mach number M Ac , and a total convective Mach number that combines the two. We derive an analytic expression of the KHI growth rate for a homogeneous flow (i.e., zero Atwood number, A=0) that converges toward both the expression for unmagnetized compressible flow and Chandrasekhar's expression for magnetized incompressible flow. Otherwise, the dispersion relation is solved numerically and allows deriving general stability diagrams of magnetized KHI for the triplet (A, M c , β −plasma) parameters. We show these parameters uniquely define all configurations for a parallel magnetic field. We also construct diagrams with respect to the convective Alfvénic Mach number, the β − plasma parameter, or the magnetic field showing which magnetic field strength is required for stabilizing a given shear flow. The theoretical growth rates are compared with 18 simulations made with the GAMERA code, currently used for 3D magnetospheric simulations. Finally, we apply our results to the analysis of a past KHI experiment performed at the OMEGA laser facility, showing linear theory succeeds to provide accurate estimates of the growth rate at early times. We further discuss how our results can inform future experiments in the high-Mach magnetized regime at the National Ignition Facility. Possible limitations of the study due to resistive, mixing, or turbulence effects are discussed.

compressible flows↗
DOE OSTI2025

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory↗
DOE OSTI2025

Phase-based velocity extraction method for photonic Doppler velocimetry with potential higher time resolution

We present an extension of the [Takeda et al., J. Opt. Soc. Am. 72, 156 (1982)] phase extraction method to heterodyne photonic Doppler velocimetry applications. The method yields results equivalent to those obtained by the short-time Fourier transform (STFT), while offering potential improvements in time resolution. Unlike STFT, which relies on window functions, such as the Hamming window, that emphasize central data points and diminish the influence of edges, the extended Takeda method utilizes all data uniformly. This uniform treatment allows for the derivation of empirical equations that directly relate velocity error to the actual time resolution rather than to the local analysis duration. The established equation provides a useful metric for both optimizing hardware configuration and guiding data analysis. Simulation and experimental results confirm that, for a given dataset, specifying a target time resolution yields consistent velocity errors for both methods. These findings underscore the Takeda method’s advantages, particularly its potential higher time resolution and reduced computational burden, making it a valuable tool for high-throughput applications such as laser dynamic compression experiments.

Computer simulation↗
DOE OSTI2026

Rapid organic carbon spiraling in a headwater stream linked with streamflow, biogeochemistry, and canopy phenology

Headwater streams are abundant worldwide and important to global biogeochemical cycles, serving as critical processors and transporters of C. C spiraling is a useful way to understand the retention and mineralization of organic C (OC) in streams. However, analyses of seasonal and interannual variability in OC spiraling are currently limited. In this study, we aimed to understand the temporal patterns and driving mechanisms of OC spiraling, which will inform our understanding of future OC changes under climate change. We used 7 y of daily data in a small headwater stream (Walker Branch, Tennessee, USA) to assess seasonal and interannual variability in OC spiraling length (S OC ) and mineralization velocity (v fOC ), as well as their potential related variables. On average, S OC in Walker Branch was ~10× shorter than in previously studied small streams, indicating strong connections between the water column and the benthic environment where OC mineralization mostly takes place. OC spiraling was faster during the more biologically active periods of spring and autumn compared with more elongated OC spiraling in summer and winter, when OC retention was lower and downstream transport was higher. Gross primary production (GPP) was most strongly related to S OC and v fOC . Photosynthetically active radiation (PAR) and NO 3 − were also positively and negatively related to v fOC , respectively. Trends toward earlier and longer canopy cover and reduced GPP and PAR may result in longer S OC and slower v fOC , reducing localized instream processing of OC and potentially shunting more OC downstream. However, long-term observations indicate reduced NO 3 − at Walker Branch, suggesting opposing effects to those of GPP and PAR, leading to faster v fOC and greater OC retention. Time-series analyses of OC spiraling in streams can enhance our understanding of current and future responses of OC processing and downstream transport to climate change, as well as implications for downstream OC dynamics.

biological activity↗