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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

First principles free energy model with dynamic magnetism for δ -plutonium

We present an ab initio free energy model derived from a fully relativistic density functional theory (DFT) electronic structure with dynamic magnetism for δ -plutonium (face-centered cubic, fcc). The DFT model is extended with orbital-orbital interaction in a parameter free orbital polarization (OP) mechanism consistent with previous modeling of plutonium. Gibbs free energy is built from components associated with the temperature dependence of the electronic structure and the corresponding electronic entropy, lattice vibrations within an anharmonic lattice dynamics model, and dynamical fluctuations of the magnetization density, i.e. magnetic fluctuations. The fluctuation model consists of transverse and longitudinal modes driven by temperature induced excitations of the DFT + OP electronic structure. The ab initio model thus incorporates fluctuating states beyond the electronic ground state. Thanks to the dynamic magnetism, the theory predicts excellent thermodynamic properties and a Gibbs free energy in accord with CALPHAD and semi-empirical modeling developed from the thermodynamic observables. The magnetic fluctuations further explain anomalous behaviors of the thermal expansion in plutonium. Specifically, a thermal expansion for the δ -plutonium system turning from positive to negative at temperatures above room temperature, a tendency for gallium to reduce and remove the negative thermal expansion depending on composition, and a positive thermal expansion for the high temperature ϵ phase.

dynamic↗
DOE OSTI2025

Finite Element Analysis of the TRUST Nonlinear Dynamics Testbed

This paper builds on prior work conducted within the Los Alamos National Laboratory (LANL) Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) program. Specifically, it builds on finite element (FE) modeling efforts for the TRUST program’s Nonlinear Dynamics (ND) testbed. Historically, the FE model for the ND testbed has exclusively utilized a Lanczos eigensolver that linearly extracts the system’s natural frequencies. This paper investigates the Abaqus 2024’s explicit dynamic solver, which implements a central difference explicit solver. The central difference method used in Abaqus can capture nonlinear material responses in structural dynamic simulations making it suitable for the ND FE model.

42 ENGINEERING↗
DOE OSTI2025

Efficient sparse state preparation via quantum walks

Continuous-time quantum walks (CTQWs) on dynamic graphs, referred to as dynamic CTQWs, are a recently introduced universal model of computation that offers a new paradigm in which to envision quantum algorithms. In this work, we develop an algorithm that converts single-edge and self-loop dynamic CTQWs to the gate model of computation. We use this mapping to introduce an efficient sparse quantum state preparation framework based on dynamic CTQWs. Our approach utilizes combinatorics techniques such as minimal hitting sets, minimum spanning trees, and shortest Hamiltonian paths to reduce the number of controlled gates required to prepare sparse states. We show that our framework encompasses the current state of the art ancilla-free sparse state preparation method by reformulating this method as a CTQW. This CTQW-based framework offers an alternative to the uniformly controlled rotation method used by Qiskit by requiring fewer CX gates when the target state has a polynomial number of non-zero amplitudes.

dynamic continuous time quantum walks↗
DOE OSTI2025

Field-Driven Out-of-Equilibrium Collective Patterns for Swarm Micro-Robotics

Soft robotics has been rapidly advancing, offering significant improvements over traditional rigid robotic systems through the use of compliant materials that enhance adaptability and interaction with the environment. However, current approaches face critical challenges, including the reliance on complex “top-down” fabrication techniques and the difficulty of wireless powering and control at the microscale. Swarm robotics introduces a paradigm shift, leveraging collective dynamics to achieve cooperative and adaptable behaviors among multiple robotic units. Inspired by nature, this “bottom-up” approach enables swarm robots to execute task-specific reconfigurations, enhancing flexibility and robustness. Field-driven active colloids emerge as a promising platform for swarm microrobotics, capable of self-propulsion and self-organization into dynamic collective patterns under external field excitation and manipulation. These systems mimic biologically inspired swarm behaviors, such as flocking and vortex formation, providing a versatile foundation for designing innovative swarm microrobots. Here, this review discusses the principles of electric and magnetic field-driven collective self-organization, focusing on the particle dynamics, the emergence of collective swarm patterns, and illustrative examples of functional swarm microrobots. It concludes with future perspectives on harnessing these systems for adaptive, scalable, and multifunctional microrobotic applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Variational AutoEncoders Reveal Intensifying GPP Extremes in Continental United States based on CESM2 Simulations

Climate extremes significantly impact terrestrial carbon cycle dynamics, necessitating robust methods for detecting and analyzing anomalous behavior in plant productivity. This study presents a novel application of variational autoencoders (VAE) for identifying extreme events in gross primary productivity (GPP) from Community Earth System Model version 2 simulations across four AR6 regions in the Continental United States. We compare VAE-based anomaly detection with traditional singular spectral analysis (SSA) methods across three time periods: 1850-80, 1950-80, and 2050-80 under SSP5-8.5 scenario. The VAE architecture employs three dense layers and a latent space with input sequence length of 12 months, training on normalized GPP time series to reconstruct the GPP and identify anomalies based on reconstruction errors. Extreme events are defined using 5th percentile thresholds applied to both VAE and SSA anomalies. Results demonstrate strong regional agreement between VAE and SSA methods in spatial patterns of extreme event frequencies, despite VAE consistently producing higher threshold values (179-756 GgC for VAE vs. 100-784 GgC for SSA across regions and periods). Both methods reveal increasing magnitudes and frequencies of negative carbon cycle extremes toward 2050-80, particularly in Western and Central North America. The VAE approach shows comparable performance to established SSA techniques while offering computational advantages and enhanced capability for capturing non-linear temporal dependencies in carbon cycle variability. This research demonstrates the potential of deep learning approaches for extremes detection and provides a foundation for improved understanding of future carbon cycle risks under future conditions.

Sharma, Bharat [ORNL] (ORCID:0000000266982487)↗
DOE OSTI2026

Quantum mechanical closure of partial differential equations with symmetries

We develop a statistical framework for the dynamical closure of spatiotemporal dynamics governed by partial differential equations. Employing the mathematical framework of quantum mechanics to embed the original classical dynamics into a quantum mechanical representation, we use the space of quantum density operators to model the unresolved degrees of freedom of the original dynamics in a statistical sense, and the framework of quantum measurement to predict their contributions to the resolved dynamics. The embedded dynamics is discretized by a positivity preserving process, leading to a compressed representation that is invariant under the dynamical symmetries of the resolved dynamics. We present a data based formulation of the closure scheme and apply it to a closure problem for the shallow water equations. The numerical results demonstrate that our closure model can accurately predict the main features of the true dynamics, including for out of sample initial conditions.

Delay embedding↗
DOE OSTI2025

Frontal Maintenance in Submesoscale Flows

Classic deformation theory includes parameters—divergence, total strain, and vorticity—that are invariant to changes in the coordinate system. However, these parameters are sometimes ambiguous with respect to characterizing how fronts are formed and maintained because the presence of a front imposes a reference coordinate system. To help remedy this ambiguity, we propose a framework in frontal coordinates based on along- and cross-front velocity gradients to better characterize frontal maintenance, which can also be used to define divergence and normal strain in frontal coordinates. The framework with these four parameters (along-, cross-front velocity gradients, divergence, and normal strain in frontal coordinate) defines eight characteristic flow types at a front, providing a complete characterization of the flow that strengthens or weakens a front. This framework highlights the importance of the “strain efficiency” concept, which unambiguously defines the contribution of total strain to frontogenesis. Two examples, one based on a realistic simulation of submesoscales in the northern Gulf of Mexico and the other based on an idealized model with similar flow characteristics, are provided to demonstrate how this framework can be used to enhance our understanding of frontal dynamics in submesoscale flows.

frontogenesis↗
DOE OSTI2022

Computational Fluid Dynamics analysis of a jet-induced attrition

The flow dynamics, collisional behavior and attrition rates of a grid jet attrition unit are computationally investigated using a kinetic theory-based approach. Profiles of the flow statistical properties are presented and discussed to describe the jet dynamics. Attrition behavior and rates are analyzed from the local reconstruction of the collision energy spectra. The commonly used jet-induced attrition mechanisms which involve the entrainment of the particles into the jet followed by their acceleration and collision at the jet tip is reproduced by the CFD calculations in this work. Here, supported by the local attrition rates, it is shown that the shear-based collisions around the jets, due to counter-current particle axial motions, are an additional jet-induced mechanism generating more fines than at the jet tips.

42 ENGINEERING↗
DOE OSTI2025

A randomized sketching trust-region secant method for low-memory dynamic optimization

The numerical solution of dynamic optimization problems is often limited by the memory required to store the state trajectory, which is used to evaluate the objective function and its derivatives. Recently, [R. Muthukumar et al., SIAM Journal on Optimization 31(2), pp. 1242–1275 (2021)] introduced a trust-region method for dynamic optimization that employs randomized sketching to compress the state trajectory, resulting in inexact derivative computations. By adaptively learning the sketch rank, the trust-region algorithm achieves rigorous convergence guarantees. Here, we extend this approach to use secant Hessian approximations. Due to the randomness introduced by the sketch, the traditional secant update formulae can produce poor Hessian approximations. In particular, the difference of two gradients, computed from two different sketches, may be inconsistent. To overcome this, we employ a sketched approximation of the Hessian application, in lieu of computing the gradient difference. We numerically demonstrate the improved stability of this approach on an example from PDE-constrained optimization.

dynamic optimization↗
DOE OSTI2025

A machine-learning approach to measure 3D sample properties from 2D Transmission Electron Microscopy images

Transmission Electron Microscopy (TEM) is a powerful tool for the characterization of materials at the nanoscale; however, its inherent two-dimensional (2D) nature poses significant challenges to accurately measure three-dimensional (3D) properties. We introduce a supervised machine-learning model that predicts 3D structural information, such as sample thickness and curvature, from a series of conventional 2D TEM images. The model, a U-Net convolutional neural network, is trained on a large synthetic dataset generated from dynamical diffraction simulations that model TEM’s complex, nonlinear image formation, accounting for sample thickness and curvature. This physically realistic framework enables exploration of a broad parameter space impractical to sample experimentally. We demonstrate that the trained model has accurate predictions for experimental single-crystal silicon samples, achieving performance comparable to established measurement techniques. This work highlights the critical role of robust, simulation-based training in overcoming the limitations of real-world imaging artifacts and inconsistent sample geometries. By integrating machine learning with numerical simulations, we offer an efficient and scalable framework for quantitative TEM analysis, paving the way for more sophisticated 3D characterization of complex materials.

Dynamical diffraction↗
DOE OSTI2025

Propagating synthetic populations with dynamic Bayesian networks: a framework for long-horizon demographic forecasting

This study presents a dynamic demographic microsimulator using dynamic Bayesian networks to forecast long–term changes in household and individual life events. Leveraging longitudinal Panel Study of Income Dynamics (PSID) data, two networks for individuals and households were modeled to simulate transitions in employment, income, education, marriage, childbirth, leaving the parental home, home ownership, mortality, and household formation or dissolution. Across 1,000 simulation runs spanning 24 years, household–level outcomes remain highly accurate and individual–level predictions reasonable. Although accuracy naturally declines with projection horizon, performance remains promising at both levels. This study addresses a key limitation of existing population synthesis models, which typically generate only a single static snapshot of the population. In conclusion, by introducing a framework that propagates cross-sectional outputs into the future, the microsimulator enables the tracking of demographic evolution over time, enhances realism in population-based simulations, and supplies credible inputs to agent-based travel demand models.

Demographic modeling↗
DOE OSTI2025

Hydrogen density mapping in biomolecular crystals through dynamic nuclear polarization

Many fundamental biological processes, including those in photosynthetic reaction centers and enzyme active sites, involve charge and energy transfer, bond cleavage, protonation and hydrogen bonding. Because H atoms play such central roles in these reactions, accurately determining their positions is essential. Yet, conventional X-ray crystallography primarily resolves the heavy atoms in biological structures and provides limited insight into hydrogen, even at atomic resolution. Neutron macromolecular crystallography (NMC) overcomes this limitation by offering exceptional sensitivity to hydrogen and deuterium. Here, we present a theoretical framework for the development of dynamic nuclear polarization NMC (DNP-NMC) techniques, which exploit the alignment of neutron and proton nuclear spins to enhance and tune the hydrogen signal contribution. The DNP-NMC approach advances the resolution of H atoms within biomolecular crystals, whether bound to protein residues or present in solvent. The method establishes key relationships for the coherent structure factor of polarized neutron scattering from hydrogenous matter. It theoretically achieves full accuracy in phase reconstruction and offers a path to improve neutron structure determination, achieving accuracies exceeding ≳80% by incorporating titration states. Using a variant of the hybrid input/output phase-retrieval algorithm, it allows recovery of the hydrogen density with ≳90% phase accuracy. In conclusion, we further discuss sources of experimental uncertainty for the upcoming DNP-enabled, quasi-Laue IMAGINE-X experiment at Oak Ridge National Laboratory's High Flux Isotope Reactor.

dynamic nuclear polarization↗
DOE OSTI2026

A new cw-NMR Q-meter for dynamically polarized targets for particle physics

Polarized solid targets produced via Dynamic Nuclear Polarization rely on Continuous-Wave Nuclear Magnetism Resonance measurements to accurately determine the degree of polarization of bulk samples polarized to nearly 100%. Since the late 1970's phase sensitive detection methods have been utilized to observe the magnetization of a sample as a small change in inductance under RF excitation near the Larmor frequency of the nuclear species of interest, using a device known as a Q-meter. Liverpool Q-meters, produced in the UK in the 80's and 90's, have been the workhorse devices for these targets for decades, however their age and scarcity has meant new systems are needed. In conclusion, we describe a Q-meter system designed and built at Jefferson Lab in the Liverpool style to have comparable electronic performance with several improvements to update and adapt the devices for modern use.

Dynamic nuclear polarization↗
DOE OSTI2021

In Silico Discovery of Multistep Chemistry Initiated by a Conical Intersection: The Challenging Case of Donor–Acceptor Stenhouse Adducts

Detailed mechanistic understanding of multistep chemical reactions triggered by internal conversion via a conical intersection is a challenging task that emphasizes limitations in theoretical and experimental techniques. We present a discovery-based, hypothesis-free computational approach based on first-principles molecular dynamics to discover and refine the switching mechanism of donor–acceptor Stenhouse adducts (DASAs). We simulate the photochemical experiment in silico, following the “hot” ground state dynamics for 10 ps after photoexcitation. Using state-of-the-art graphical processing units-enabled electronic structure calculations we performed in total ~2 ns of nonadiabatic ab initio molecular dynamics discovering (a) critical intermediates that are involved in the open-to-closed transformation, (b) several competing pathways which lower the overall switching yield, and (c) key elements for future design strategies. Finally, our dynamics describe the natural evolution of both the nuclear and electronic degrees of freedom that govern the interconversion between DASA ground-state intermediates, exposing significant elements for future design strategies of molecular switches.

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

The impact of stacking and phonon environment on energy transfer in organic chromophores: computational insights

Energy transfer in organic materials is extensively studied due to many applications in optoelectronics. The electronic and vibrational relaxations within molecular assemblies can be influenced by stacking arrangements or additions of a backbone that unites them. Here, we present the computational study of the photoexcitation dynamics of a perylene diimide monomer, and face-to-face stacked dimer and trimer. By using non-Energy transfer in organic materials is extensively studied due to many applications in optoelectronics. The electronic and vibrational relaxations within molecular assemblies can be influenced by stacking arrangements or additions of a backbone that unites them. In this work, we present the computational study of the photoexcitation dynamics of a perylene diimide monomer, and face-to-face stacked dimer and trimer. By using non-adiabatic excited-state molecular dynamics simulations, we show that the non-radiative relaxation is accelerated with the number of stacked molecules. This effect is explained by differences in the energy splitting between states that impacts their corresponding nonadiabatic couplings. Additionally, our analysis of the vibronic dynamics reveals that the passage through the different conical intersections that participate in the relaxation of the stacked systems, activate a positive feedback mechanism. This effect involves a narrow set of vibrational normal modes that accelerate the process by increasing the efficiency of its vibronic dynamics. In contrast, an addition of a biologically inspired backbone slows down the relaxation rate due to its participation in the vibronic dynamics of the molecular stacking arrangements. Our results suggest the stacking arrangements and common backbones as strategies to modulate the efficiency of electronic and vibrational relaxation of diimide-based systems and other molecular aggregates.

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

CFD modeling of non-catalytic, partial-oxidation engine reformer for flare mitigation

Flaring associated natural gas is commonly employed in the oil and gas industry to reduce methane (CH 4 ) emissions but generates carbon dioxide (CO 2 ) and harmful pollutants, significantly contributing to air pollution and posing risks to public health. To mitigate this impact, M2X Energy Inc. has developed a small-scale, modular gas-to-methanol system. This system features an engine reformer that performs fuel-rich partial oxidation of wellhead gas to produce syngas—a mixture of carbon monoxide (CO) and hydrogen (H 2 )—followed by a downstream reactor for methanol synthesis. This study focused on computational fluid dynamics (CFD) modeling of the engine reformer to simulate partial oxidation chemistry, predict the rich-burn operating limit, and assess syngas quality, ultimately aiding in design and operational optimization. The CFD model, developed within a Reynolds-Averaged Navier-Stokes (RANS) turbulence framework, incorporated sub-models for turbulent combustion, a chemical mechanism with polycyclic aromatic hydrocarbon (PAH) pathways, and soot emissions to accurately capture the fuel-rich, turbulent jet ignition and combustion processes. Model validation against experimental data showed good agreement across pre- and main-chamber pressures, apparent heat release rates, and exhaust gas concentrations of key species (H 2 , CO, CO 2 , CH 4 ) for varying intake equivalence ratios. Here, the model identified a rich-burn operating limit near a fuel-air equivalence ratio of 2.35, consistent with experimental observations. Furthermore, syngas quality analysis revealed that extending the rich-burn limit through engine reformer optimization could enhance syngas production, contributing to higher methanol synthesis efficiency.

Computational Fluid Dynamics↗
DOE OSTI2025

A multi model ensemble reveals net climate benefits from regenerative practices in US Midwest croplands

Process-based cropping systems models (CSMs) are key components of measurement, monitoring, reporting, and verification frameworks of carbon markets, but model-specific differences limit their applicability across diverse pedo-climatic conditions and agronomic practices. Multi-model ensemble (MME) provides an opportunity to better estimate changes in soil organic carbon (SOC) and nitrous oxide (N 2 O) emissions from agronomic practices at scale. We used an MME across 46 million hectares of US Midwest cropland at a resolution of 4-km 2 to assess the aggregate ability of different regenerative practices to sequester SOC and N 2 O emissions compared to their counterfactual dynamic baselines. MME was validated against long-term trials and compared to its constituent CSMs, showing greater accuracy and lower uncertainty. The results show that adopting no-till combined with cover crops increased SOC stocks by 0.36 ± 0.12 Mg ha -1 yr -1 , corresponding to a net regional SOC gain of 16.4 Tg C yr -1 compared to business-as-usual baselines. These benefits are halved when each management is practiced individually, and the SOC gains are only fully realized with low initial carbon stock. By including N 2 O emissions, we can assess the overall climate mitigation potential, specifically, the extent to which carbon sequestration can offset direct N 2 O emissions. The magnitude of this potential varies depending on management practices and geographic location with net climate benefits on average ranging from 0 to 3 Mg CO 2 -eq ha -1 yr -1 . High-resolution MME results allow for robust estimates of climate mitigation, reducing barriers to carbon market participation and supporting regenerative agriculture initiatives at scale.

carbon credits↗
DOE OSTI2023

DeePMD-kit v2: A software package for deep potential models

DeePMD-kit is a powerful open-source software package that facilitates molecular dynamics simulations using machine learning potentials known as Deep Potential (DP) models. This package, which was released in 2017, has been widely used in the fields of physics, chemistry, biology, and material science for studying atomistic systems. The current version of DeePMD-kit offers numerous advanced features, such as DeepPot-SE, attention-based and hybrid descriptors, the ability to fit tensile properties, type embedding, model deviation, DP-range correction, DP long range, graphics processing unit support for customized operators, model compression, non-von Neumann molecular dynamics, and improved usability, including documentation, compiled binary packages, graphical user interfaces, and application programming interfaces. This article presents an overview of the current major version of the DeePMD-kit package, highlighting its features and technical details. Additionally, this article presents a comprehensive procedure for conducting molecular dynamics as a representative application, benchmarks the accuracy and efficiency of different models, and discusses ongoing developments.

97 MATHEMATICS AND COMPUTING↗