Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamical System Scaling”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Developing a microwave-driven reactor for ammonia synthesis: insights into the unique challenges of microwave catalysis

Rapid development of new ammonia (NH 3 ) synthesis techniques that enable modular, intermittent production is essential to actualizing NH 3 's potential as a clean energy carrier, since contemporary methods are configured to centralized, continuous production methods with high emissions and are incompatible with renewable sources. Here, in this mission, microwave-driven catalysis is promising for its ability to enhance reaction kinetics and apply targeted heating for efficient energy use. However, owing to an incomplete understanding of the interaction between microwave fields and catalyst beds, the development of such microwave-catalysis systems remains underexplored and challenging. This paper investigates the 10× scale-up of a microwave-based NH 3 synthesis reactor using numerical and experimental approaches, achieving the largest reported microwave-driven NH 3 reactor to date. Results elucidate phenomena unique to microwave processes that are only predictable through numerical modeling, including how a catalyst's dielectric properties influence microwave field distribution by affecting penetration depth and how energy utilization can be poor even with sufficient attenuation. These dynamics change with scale, constrain reactor geometry, and potentially hamper performance. Nonetheless, we demonstrate a production rate of 56.6 g NH$_3$ per day, the highest reported NH 3 synthesis rate for laboratory-scale alternative techniques; correspondingly, the benchmark energy efficiency achieved in this paper (45.6 g NH$_3$ kW h -1 ) is the highest reported for such reactors of sufficient scale. Even with this exemplary energy efficiency, energy losses were found in excess of 50%, an issue resolvable through scale and reactor design. The efficiencies imparted by microwaves were key in these achievements, warranting further investigation toward development of microwave-driven NH 3 systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Synthetic Digital Environments for Training Robots on Earth & Beyond [Poster]

Creating physical replicas of real-world environments to train robots for challenging outdoor tasks, whether constructing energy infrastructure like solar farms on Earth or on the Moon and Mars, is prohibitively expensive. This project will prototype a high-fidelity digital twin framework using NVIDIA IsaacSim to create realistic digital representations of robotic systems and their operating conditions, including varied terrains and environmental factors, allowing robots to learn and adapt in a faster, safer, and more affordable way to tackle unpredictable challenges in terrestrial and extraterrestrial applications. The Robotic Space Exploration (RoSE) Lab at Colorado School of Mines focused on the development and testing of synthetic digital twins to explore multi-physics interactions between robots and unstructured environments, with emphasis on lunar conditions such as deformable regolith, reduced gravity, and terrain-robot contact dynamics. The Industrialized Construction Innovation (ICI) team at National laboratory of the Rockies (NLR), simulated robotic apparatus and construction workflows using synthetic digital twins to inform real-world deployment, targeting application-driven use cases such as robotic construction of a scaled prototype of a photovoltaic energy infrastructure. Joint efforts between RoSE and ICI are continuing to explore how environment-scale multi-physics modeling and application-level robotic system simulation could be integrated to support robotic construction of energy infrastructure in highly unstructured environments, including scenarios relevant to the Lunar South Pole.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Gauge-fixing quantum density operators at scale

We provide a theory, algorithms, and simulations of nonequilibrium quantum systems using a one-dimensional (1D) completely positive (CP), matrix-product (MP) density-operator (𝜌) representation. By generalizing the matrix product state's orthogonality center, to additionally store positive classical mixture correlations, the MP⁢𝜌 factorization naturally emerges. In this setting, we analytically and numerically examine the virtual gauge freedoms associated with the representation of quantum density operators. Based on this perspective, we simplify algorithms in certain limits to speed up the integration of the canonical-form master-equation dynamics. This enables us to quickly evolve under the dynamics of two-body quantum channels without resorting to optimization-based methods. In addition to this technical advance, we also scale up numerical examples and discuss implications for accurately modeling hardware architectures and predicting their performance in the near term. This includes an example of the quantum to classical transition of informationally leaky, i.e., decohering, qubits. In this setting, because of loss from environmental interactions, nonlocal complex coherence correlations are converted into global incoherent classical statistical mixture correlations. Lastly, the representation of both global and local correlations is discussed. We expect this work to have applications in additional nonequilibrium settings, beyond qubit engineering.

Gangapuram, Amit Jamadagni [Oak Ridge National Lab↗
DOE OSTI2022

A scalable planning framework of energy storage systems under frequency dynamics constraints

As the penetration of renewables increases in power systems, the declining system inertia can cause frequency stability issues. Battery energy storage systems (BESSs) respond fast and therefore can relieve the low inertia difficulty but need to be appropriately sized considering the associated cost. This paper presents a novel stochastic optimization model for economically planning BESS capacity while considering the spatial–temporal correlation of wind generation and generator outages under frequency stability constraints, which include the rate-of-change of frequency (RoCoF), frequency nadir (FN), and quasi-steady-state (QSS) frequency. A set of new FN constraints that can be easily linearized is developed. To account for renewable uncertainties, a realistic uncertainty modeling approach, Random Field, is adopted to generate wind generation scenarios by considering both spatial and temporal evolutions of wind speed profiles. The ESS sizing is formulated as a mixed-integer linear programming problem and solved by using a scalable decomposition-and-coordination approach, Surrogate Absolute Value Lagrangian Relaxation (SAVLR). To further improve the scalability and reduce computational burdens, a rolling-horizon-based update is developed and incorporated into SAVLR for providing a practical solution to the long-term planning of very large-scale power systems. Finally, a modified IEEE 118-bus system and the Polish system are used to validate the effectiveness and scalability of the model and solution methodology.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Real-Space Constrained Density Functional Theory Investigation of Site-Specific, Interfacial Charge Recombination Dynamics Across the Au Nanoparticle/TiO 2 Heterojunction

Au nanoparticle (NP)/TiO 2 heterojunction is a representative system to study interfacial charge transfer in photocatalysis and photovoltaics, where suppressing recombination from TiO 2 to Au can enhance hot carrier extraction. We apply real-space constrained density functional theory (CDFT) with Marcus theory to quantify charge recombination time scales across Au/TiO 2 . This approach enables direct control and visualization of charge-separated states, aligning with site-specific probes like time-resolved X-ray photoelectron spectroscopy (trXPS). We find that the charge-separated state features a bipolaron, with recombination dominated by TiO 2 LUMO to Au HOMO transitions, primarily at interfacial Au sites. Marcus rate predictions are benchmarked with surface hopping methods, quantifying differences in time scales and computational efficiency. Lastly, we examine how the Au cluster size affects the free energy change (ΔG) and reorganization energy (λ), explaining trends in closed-shell systems and highlighting challenges for open-shell extrapolations. Overall, CDFT + Marcus theory provides efficient, mechanistically transparent interfacial charge transfer modeling, and we clearly defined its applicability and limitation.

Glenna, Drew M. [Univ. of Idaho, Idaho Falls, ID (↗
DOE OSTI2020

Quantum Versus Classical Switching Dynamics of Driven Dissipative Kerr Resonators

Here, we report a first-principles study of the driven dissipative dynamics for Kerr oscillators in the mesoscopic regime. This regime is characterized by large Kerr nonlinearity, realized here using the nonlinear kinetic inductance of a large array of Josephson junctions. The experimentally measured nonlinear resonance lineshapes of the junction array modes show significant deviations from steady-state numerical predictions, and necessitate time-dependent numerical simulations indicative of strong measurement-induced dephasing in the system arising from the large cross-Kerr effect between array modes. Analytical and numerical calculations of switching rate corroborate this by showing the emergence of a slow time scale, which is much longer than the linear decay rate and is set by fluctuation-induced switching times in the bistable regime. Furthermore, our analysis shows that the usual quantum-activated escape treatment is inadequate for prediction of the switching rates at large frequency shifts caused by strong nonlinearities, necessitating a quantum treatment that utilizes the full system Liouvillian. Based on our analysis, we identify a universal crossover parameter that delineates the regimes of validity of semiclassical and quantum descriptions, respectively. Our work shows how dynamical switching effects in strongly nonlinear systems provide a platform to study quantum-to-classical transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Origins of Lithium/Sodium Reverse Permeability Selectivity in 12-Crown-4-Functionalized Polymer Membranes

Direct lithium extraction via membrane separations has been fundamentally limited by lack of monovalent ion selectivity exhibited by conventional polymeric membranes, particularly between sodium and lithium ions. Recently, a 12- Crown-4-functionalized polynorbornene membrane was shown to have the largest lithium/sodium permeability selectivity observed in a fully aqueous system to date. Using atomistic molecular dynamics simulations, we reveal that this selectivity is due to strong interactions between sodium ions and 12-Crown-4 moieties, which reduce sodium ion diffusivity while leaving lithium ion mobility relatively unaffected. Furthermore, the ion diffusivities in the membrane, when scaled by their respective solution diffusivities and free ion fractions, can be collapsed to an almost universal relationship depending on solvent volume fraction.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Comparisons Between Array Derived Dynamic Strain Rate (ADDS) and Fiber‐Optic Distributed Acoustic Sensing (DAS) Strain Rate

Abstract Distributed acoustic sensing (DAS) strain rate and particle velocity can be compared through approximate scaling with medium velocity. We instead performed a direct comparison between array derived dynamic strain (ADDS) rate and DAS strain rate for six frequency bands. The PoroTomo project at Brady's Hot Springs, Nevada, deployed a 240‐geophone 3C array co‐located with fiber‐optic DAS system and 8.7 km of buried cable. We selected subsets of the geophone array to create four smaller arrays and computed ADDS. The horizontal components of the ADDS were rotated into the direction of the fiber‐optic cable and then compared with the observed DAS strain rates. From three example regional earthquakes of local magnitudes 2.9, 4.1, and 4.3, the ADDS are found to be coherent with DAS for frequencies ≤1 Hz. For frequencies >1‐Hz, this correlation decays quickly. Small differences between linear and areal dynamic strains at 1‐Hz suggest poor signal‐to‐noise or localized strain that is perturbed by shallow heterogeneities compare to the average strain propagating across the geophone array. The implication is that around 1‐Hz, straight fiber DAS is measuring axial strain along the fiber and can provide good approximations to translational particle motions. However, above 1‐Hz, DAS becomes more sensitive to shallow velocity gradients that can be beneficial for geophysical imaging yet becomes a limitation for traditional seismic analysis methods depending on absolute amplitude and phase from translational particle motions.

58 GEOSCIENCES↗
DOE OSTI2024

Integrated Effects of Site Hydrology and Vegetation on Exchange Fluxes and Nutrient Cycling at a Coastal Terrestrial‐Aquatic Interface

Abstract The complex interactions among soil, vegetation, and site hydrologic conditions driven by precipitation and tidal cycles control the biogeochemical transformations and bi‐directional exchange of carbon and nutrients across the terrestrial–aquatic interfaces (TAIs) in coastal regions. This study uses a highly mechanistic model, Advanced Terrestrial Simulator (ATS)‐PFLOTRAN, to explore how these interactions affect exchanges of materials and carbon and nitrogen cycling. We used a transect in the Chesapeake Bay region that spans zones of open water, coastal wetland, transition, and upland forest. We designed several simulation scenarios to parse the effects of the individual controlling factors and the sensitivity of carbon cycling to reaction rate parameters derived from laboratory experiments. Our simulations reveal an active zone for carbon cycling under the transition zones between the wetland and the upland. Evapotranspiration is found to enhance the exchange fluxes between the surface and subsurface domains, resulting in a higher dissolved oxygen concentration in the TAIs. The transport of organic carbon derived from plant leaves and roots provide an additional source of organic carbon needed for the aerobic respiration and denitrification processes in the TAIs. The variability in reaction rate parameters associated with microbial activities is also found to play a dominant role in controlling the heterogeneity and dynamics of the simulated redox conditions. This modeling‐focused exploratory study enabled us to better understand the complex interactions among soil, water and microbes that govern the hydro‐biogeochemical processes at the TAIs, which is an important step toward representing coastal ecosystems in larger‐scale Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

AI-NERD: Elucidation of relaxation dynamics beyond equilibrium through AI-informed X-ray photon correlation spectroscopy

Abstract Understanding and interpreting dynamics of functional materials in situ is a grand challenge in physics and materials science due to the difficulty of experimentally probing materials at varied length and time scales. X-ray photon correlation spectroscopy (XPCS) is uniquely well-suited for characterizing materials dynamics over wide-ranging time scales. However, spatial and temporal heterogeneity in material behavior can make interpretation of experimental XPCS data difficult. In this work, we have developed an unsupervised deep learning (DL) framework for automated classification of relaxation dynamics from experimental data without requiring any prior physical knowledge of the system. We demonstrate how this method can be used to accelerate exploration of large datasets to identify samples of interest, and we apply this approach to directly correlate microscopic dynamics with macroscopic properties of a model system. Importantly, this DL framework is material and process agnostic, marking a concrete step towards autonomous materials discovery.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Light-driven nanoscale vectorial currents

Controlled charge flows are fundamental to many areas of science and technology, serving as carriers of energy and information, as probes of material properties and dynamics and as a means of revealing or even inducing broken symmetries. Emerging methods for light-based current control offer particularly promising routes beyond the speed and adaptability limitations of conventional voltage-driven systems. However, optical generation and manipulation of currents at nanometre spatial scales remains a basic challenge and a crucial step towards scalable optoelectronic systems for microelectronics and information science. Here we introduce vectorial optoelectronic metasurfaces in which ultrafast light pulses induce local directional charge flows around symmetry-broken plasmonic nanostructures, with tunable responses and arbitrary patterning down to subdiffractive nanometre scales. Local symmetries and vectorial currents are revealed by polarization-dependent and wavelength-sensitive electrical readout and terahertz (THz) emission, whereas spatially tailored global currents are demonstrated in the direct generation of elusive broadband THz vector beams 17 . We show that, in graphene, a detailed interplay between electrodynamic, thermodynamic and hydrodynamic degrees of freedom gives rise to rapidly evolving nanoscale driving forces and charge flows under the extremely spatially and temporally localized excitation. These results set the stage for versatile patterning and optical control over nanoscale currents in materials diagnostics, THz spectroscopies, nanomagnetism and ultrafast information processing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accuracy, transferability, and computational efficiency of interatomic potentials for simulations of carbon under extreme conditions

Large-scale atomistic molecular dynamics (MD) simulations provide an exceptional opportunity to advance the fundamental understanding of carbon under extreme conditions of high pressures and temperatures. However, the fidelity of these simulations depends heavily on the accuracy of classical interatomic potentials governing the dynamics of many-atom systems. Here, this study critically assesses several popular empirical potentials for carbon, as well as machine learning interatomic potentials (MLIPs), in their ability to simulate a range of physical properties at high pressures and temperatures, including the diamond equation of state, its melting line, shock Hugoniot, uniaxial compressions, and the structure of liquid carbon. Empirical potentials fail to accurately predict the behavior of carbon under high pressure–temperature conditions. In contrast, MLIPs demonstrate quantum accuracy, with Spectral Neighbor Analysis Potential (SNAP) and atomic cluster expansion (ACE) being the most accurate in reproducing the density functional theory results. ACE displays remarkable transferability despite not being specifically trained for extreme conditions. Furthermore, ACE and SNAP exhibit superior computational performance on graphics processing unit-based systems in billion atom MD simulations, with SNAP emerging as the fastest. In addition to offering practical guidance in selecting an interatomic potential with a fine balance of accuracy, transferability, and computational efficiency, this work also highlights transformative opportunities for groundbreaking scientific discoveries facilitated by quantum-accurate MD simulations with MLIPs on emerging exascale supercomputers.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solving the Hele–Shaw flow using the Harrow–Hassidim–Lloyd algorithm on superconducting devices: A study of efficiency and challenges

The development of quantum processors for practical fluid flow problems is a promising yet distant goal. Recent advances in quantum linear solvers have highlighted their potential for classical fluid dynamics. In this study, we evaluate the Harrow–Hassidim–Lloyd (HHL) quantum linear systems algorithm (QLSA) for solving the idealized Hele–Shaw flow. Our focus is on the accuracy and computational cost of the HHL solver, which we find to be sensitive to the condition number, scaling exponentially with problem size. This emphasizes the need for preconditioning to enhance the practical use of QLSAs in fluid flow applications. Moreover, we perform shots-based simulations on quantum simulators and test the HHL solver on superconducting quantum devices, where noise, large circuit depths, and gate errors limit performance. Error suppression and mitigation techniques improve accuracy, suggesting that such fluid flow problems can benchmark noise mitigation efforts. Finally, our findings provide a foundation for future, more complex application of QLSAs in fluid flow simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.

Cage effect↗
DOE OSTI2021

Sterile neutrinos, black hole vacuum and holographic principle

We construct an effective field theory (EFT) model that describes matter field interactions with Schwarzschild mini-black-holes (SBH’s), treated as a scalar field, B 0 (x). Fermion interactions with SBH’s require a complex spurion field, θ ij , which we interpret as the EFT description of “holographic information,” which is correlated with the SBH as a composite system. We consider Hawking’s virtual black hole vacuum (VBH) as a Higgs phase, $\langle$B 0 $\rangle$=V. Integrating sterile neutrino loops, the information field θij is promoted to a dynamical field, necessarily developing a tachyonic instability and acquiring a VEV of order the Planck scale. For N sterile neutrinos this breaks the vacuum to SU(N) × U(1)/SO(N) with N degenerate Majorana masses, and 12N(N+1) Nambu-Goldstone neutrino-Majorons. The model suggests many scalars fields, corresponding to all fermion bilinears, may exist bound nonperturbatively by gravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

ddcMD.os

ddcMD is a general purpose molecular dynamics (MD) code that supports MPI parallelism. MD codes are used for simulation of particles systems and capture all the many-body effects of the underlying particle potential that defines the physical system. Though MD can be used to model systems from the subatomic to astrological length scales ddcMD is mainly focused on the atomic scale length scale, In this release of ddcMD the support will be mainly for systems using the coarse-grain Martini potential, a particle potential for biological systems.

Glosli, JamesN↗
DOE OSTI2019

A two-phase substrate model for enzymatic hydrolysis of lignocellulose: application to batch and continuous reactors

Abstract Background Enzymatic hydrolysis continues to have a significant projected production cost for the biological conversion of biomass to fuels and chemicals, motivating research into improved enzyme and reactor technologies in order to reduce enzyme usage and equipment costs. However, technology development is stymied by a lack of accurate and computationally accessible enzymatic-hydrolysis reaction models. Enzymatic deconstruction of cellulosic materials is an exceedingly complex physico-chemical process. Models which elucidate specific mechanisms of deconstruction are often too computationally intensive to be accessible in process or multi-physics simulations, and empirical models are often too inflexible to be effectively applied outside of their batch contexts. In this paper, we employ a phenomenological modeling approach to represent rate slowdown due to substrate structure (implemented as two substrate phases) and feedback inhibition, and apply the model to a continuous reactor system. Results A phenomenological model was developed in order to predict glucose and solids concentrations in batch and continuous enzymatic-hydrolysis reactors from which liquor is continuously removed by ultrafiltration. A series of batch experiments were performed, varying initial conditions (solids, enzyme, and sugar concentrations), and best-fit model parameters were determined using constrained nonlinear least-squares methods. The model achieved a good fit for overall sugar yield and insoluble solids concentration, as well as for the reduced rate of sugar production over time. Additionally, without refitting model coefficients, good quantitative agreement was observed between results from continuous enzymatic-hydrolysis experiments and model predictions. Finally, the sensitivity of the model to its parameters is explored and discussed. Conclusions Although the phenomena represented by the model correspond to behaviors that emerge from clusters of mechanisms, and hence a set of model coefficients are unique to the substrate and the enzyme system, the model is efficient to solve and may be applied to novel reactor schema and implemented in computational fluid dynamics (CFD) simulations. Hence, this modeling approach finds the right balance between model complexity and computational efficiency. These capabilities have broad application to reactor design, scale-up, and process optimization.

09 BIOMASS FUELS↗