Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic timing analysis”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

What is (quantitative) system dynamics modeling? Defining characteristics and the opportunities they create

A clear definition of system dynamics modeling can provide shared understanding and clarify the impact of the field. We introduce a set of characteristics that define quantitative system dynamics, selected to capture core philosophy, describe theoretical and practical principles, and apply to historical work but be flexible enough to remain relevant as the field progresses. The defining characteristics are: (1) models are based on causal feedback structure, (2) accumulations and delays are foundational, (3) models are equation-based, (4) concept of time is continuous, and (5) analysis focuses on feedback dynamics. We discuss the implications of these principles and use them to identify research opportunities in which the system dynamics field can advance. These research opportunities include causality, disaggregation, data science and AI, and contributing to scientific advancement. Progress in these areas has the potential to improve both the science and practice of system dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Detailed biomass fast pyrolysis kinetics integrated to computational fluid dynamic (CFD) and discrete element modeling framework: Predicting product yields at the bench-scale

Fast pyrolysis is an intricate process due to the variability and anisotropy of lignocellulosic biomass and the complicated chemistry and physics during conversion in a bubbling fluidized bed reactor (BFBR). The complexity of biomass fast pyrolysis lends itself well to computational fluid dynamics (CFD) and discrete element (DEM) analysis, which promises to reduce experimental time and its associated cost. This work investigated switchgrass fast pyrolysis simulated by computational fluid dynamics coupled with a discrete element method to track individual reacting biomass particles throughout a bench-scale BFBR reactor. We accounted for the fast pyrolysis chemistry through a comprehensive reaction scheme with secondary cracking reactions. We performed a three-step reduction for secondary cracking reactions to convert the full cracking scheme into a reduced scheme easily incorporated into our model. We assessed the impact of operational conditions on the steady-state yields of liquid bio-oil, non-condensable gases (NCG), at 550 °C over a range of fluidization numbers (2 – 6 Umf), reported as a ratio to the minimum fluidization velocity (Umf). At steady-state, the volatile bio-oil yield had a range of 49.3–50.4 wt%. Levoglucosan was the primary volatile component present with 21 wt% of the bio-oil while water was the second largest with 20 wt%. The reduction of the secondary reaction schemes did not appreciably affect the overall yields of switchgrass pyrolysis compared to the full secondary scheme.

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

The Discovery of the Long-Period, Eccentric Planet Kepler-88 d and System Characterization with Radial Velocities and Photodynamical Analysis

We present the discovery of Kepler-88 d (P{sub d} = 1403±14 days, Msini{sub d}=965±44 M{sub ⊕}=3.04±0.13 M{sub J}, e{sub d}=0.41±0.03) based on six years of radial velocity (RV) follow-up from the W. M. Keck Observatory High Resolution Echelle Spectrometer spectrograph. Kepler-88 has two previously identified planets. Kepler-88 b (KOI-142.01) transits in the NASA Kepler photometry and has very large transit timing variations (TTVs). Nesvorný et al. performed a dynamical analysis of the TTVs to uniquely identify the orbital period and mass of the perturbing planet (Kepler-88 c), which was later was confirmed with RVs from the Observatoire de Haute-Provence (OHP). To fully explore the architecture of this system, we performed photodynamical modeling on the Kepler photometry combined with the RVs from Keck and OHP and stellar parameters from spectroscopy and Gaia. Planet d is not detectable in the photometry, and long-baseline RVs are needed to ascertain its presence. A photodynamical model simultaneously optimized to fit the RVs and Kepler photometry yields the most precise planet masses and orbital properties yet for b and c: P{sub b}=10.91647±0.00014 days, M{sub b}=9.5±1.2 M{sub ⊕}, P{sub c}=22.2649±0.0007 days, and M{sub c}=214.0±5.3 M{sub ⊕}. The photodynamical solution also finds that planets b and c have low eccentricites and low mutual inclination, are apsidally anti-aligned, and have conjunctions on the same hemisphere of the star. Continued RV follow-up of systems with small planets will improve our understanding of the link between inner planetary system architectures and giant planets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Disruption event characterization and forecasting in tokamaks

Disruption prediction and avoidance is a critical need for next-step tokamaks, such as ITER. Disruption Event Characterization and Forecasting (DECAF) research fully automates analysis of tokamak data to determine chains of events that lead to disruptions and to forecast their evolution allowing sufficient time for mitigation or complete avoidance of the disruption. Disruption event chains related to local rotating or global magnetohydrodynamic (MHD) modes and vertical instability are examined with warnings issued for many off-normal physics events, including density limits, plasma dynamics, confinement transitions, and profile variations. Along with Greenwald density limit evaluation, a local radiative island power balance theory is evaluated and compared to the observation of island growth. Automated decomposition and analysis of rotating tearing modes produce physical event chains leading to disruptions. A total MHD state warning model comprised of 15 separate criteria produces a disruption forecast about 180 ms before a standard locked mode detector warning. Single DECAF event analyses have begun on KSTAR, MAST, and NSTX/-U databases with thousands of shot seconds of device operation using from 0.5 to 1 × 10 6 tested sample times per device. An initial multi-device database comparison illustrates a highly important result that plasma disruptivity does not need to increase as βN increases. Global MHD instabilities, such as resistive wall modes (RWMs), can give the briefest time period of warning before disruption compared to other physics events. In an NSTX database with unstable RWMs, the mode onset, loss of boundary and current control, and disruption event warnings are found in all cases and vertical displacement events are found in 91% of cases. An initial time-dependent reduced physics model of kinetic RWM stabilization created to forecast the disruption chain predicts instability 84% of the time for experimentally unstable cases with a relatively low false positive rate. Instances of the disruption event chain analysis illustrate dynamics including H–L back transitions for rotating MHD and global RWM triggering events. Furthermore, disruption warnings are issued with sufficient time before the disruption (on transport timescales) to potentially allow active profile control for disruption avoidance, active mode control, or mitigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dynamic mode decomposition for gyrokinetic eigenmode analysis

Dynamic mode decomposition (DMD) is a post-processing approach to decompose a complex time series into a set of modes via spectral analysis. DMD provides a new and powerful method to recover gyrokinetic drift-wave eigenfrequencies and eigenfunctions based only on the solution of the gyrokinetic-Maxwell initial value problem with almost no added cost to the initial value solver. In the present paper, DMD is applied to the CGYRO gyrokinetic code using a newly-developed CGYRO-DMD post-processor. CGYRO-DMD is numerically efficient, even on a single CPU. It does not set any restrictions on the plasma shape, beta (ratio of the plasma pressure to the magnetic field pressure), collisionality or number of species, and allows one to resolve numerous eigenmodes, even of comparable growth rates. In addition, DMD is not limited to unstable modes, but rather can capture stable and unstable branches simultaneously. In this work, we illustrate the accuracy of DMD through gyrokinetic analysis of mode transition for electromagnetic drift wave instabilities.

drift-wave eigenmodes↗
DOE OSTI2025

Engineering micromotion in Floquet prethermalization via space-time symmetries

We present a systematic framework for Floquet prethermalization under strong resonant driving, emphasizing the pivotal role of dynamical space-time symmetries. Our approach demonstrates how dynamical space-time symmetries map onto the projective static symmetry group of the prethermal Hamiltonian governing the prethermal regime. We introduce techniques for detecting dynamical symmetries through the time evolution of local observables, facilitating a detailed analysis of micromotion within each period and surpassing the limitations of conventional stroboscopic Floquet prethermal dynamics. To implement this framework, we present a prethermal protocol that preserves order-2 dynamical symmetry in a spin-ladder model, confirming the predicted relationships between the expectation values of local observables at distinct temporal points in the Floquet cycle, linked by this symmetry.

floquet systems↗
DOE OSTI2021

CSPlib

This computational singular perturbation library (CSPlib) software is useful for analysis of chemical kinetic systems. The code is used to analyze and describe the dynamical landscape underlying the time evolution of a chemical system. Given a database of state vectors, it allows the analysis of the source term and Jacobian matrix of the chemical system, the identification and decoupling of fast and slow processes, and the analysis of cause-and-effect relationships in the coupled system dynamics. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-2391 O

Najm, HabibN.↗
DOE OSTI2022

Training data selection for event classification in a highly variable environment

A problem of interest for nuclear nonproliferation is monitoring activities at nuclear facilities, where proliferation events may only take place a few times and often under variable conditions. Machine learning has revolutionized data analytics by enabling the use of measurable signatures to generate predictive models of facility operations. However, traditional methods for training these models require large, reliable data sets with labeled observations, a challenge for nonproliferation. Highly variable conditions further complicate this as events from training data may have occurred in conditions quite different from the event of interest. Our hypothesis is that when events occur in a highly variable environment, careful training data selection for each test event could outperform the standard approach of using all available training data. We developed a method to optimize training data selection for the given test event and applied it to predicting the power level of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory. In this study, the reactor startup exhibits variability between occurrences due to natural variability in environmental conditions and operational procedures. Using a combination of analysis techniques, a similitude assessment was performed on data collected from HFIR to isolate clusters that were optimal for training a predictive model. Concepts such as dynamic time warping and Jaccard similarity were used in conjunction with clustering analysis. In order to validate this approach, the model was trained on every combination of unique training events and the predictive performance was compared to the performance using a subset of the training data selected by isolated clusters found through the similitude assessment.

Iyer, A↗
DOE OSTI2026

Enabling topography-resolving structural dynamic contact simulation

Damping of structures and systems is often dominated by frictional dissipation in connections, the prediction of which remains a longstanding scientific challenge. Previous studies have shown that the actual topography of contact interfaces may have a strong effect, especially in the partial slip/liftoff regime. We recently proposed a multi-scale method, which couples finite element and boundary element modeling. The primary benefit of this approach is that it permits to analyze the effect of the actual contact topography on the dynamics of jointed structures. While this multi-scale modeling method was initially developed for quasi-static analysis, we demonstrate herein how it can be used for time step integration and Harmonic Balance analysis. We cross-verify those fully dynamic analysis methods against each other and quasi-static results, for the S4 Beam benchmark. We compare the multi-scale method against state-of-the-art full-FE analysis, in terms of numerical damping and computational performance. Some discrepancy is found to be of physical origin. Depending on the load history, it is shown that the system settles to a slightly different equilibrium. Finally, transient multi-scale simulations enable the prediction of this interesting phenomenon, for the first time, for a structure with bolted joints.

Frictional-unilateral contact↗
DOE OSTI2024

Analog Computing for Science

Conventional digital computing faces fundamental physical limits: large scale computing systems already con sume tens of Megawatts of power, Dennard scaling has ended, and data movement costs dominate application performance. Next generation experimental facilities generate data at rates that overwhelm conventional pro cessing and demand real-time analysis at the source. Analog computing, which exploits the continuous dynamics of physical systems to perform computation, promises a transformative path toward orders-of-magnitude gains in energy efficiency and time-to-solution for scientific workloads.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Nonadiabatic dynamics of photoexcited thiopyridone isomers: An interplay between El-Sayed’s conditions and energy gap law

We investigated the non-adiabatic dynamics of photoexcited thiopyridone systems across their ortho-, meta-, and para-isomeric forms. The relaxation pathways of the three isomers in both gas phase and solvent environments are mapped using surface hopping dynamics based on time-dependent density functional theory. Our analysis highlights the influence of isomeric structures on photophysical behavior, offering insights into design principles to control photochemical phenomena. The simulations suggest a systematic reduction in the rate of intersystem crossing (ISC) from ortho- to meta- to para-isomer. Comparisons with multiconfigurational wave function methods in the gas phase further demonstrate how electronic structure influences the predicted dynamical pathways. The simulated dynamics demonstrates that the spin–orbit coupling strength alone does not determine the rate of ISC, as both state energetics and underlying electronic and structural features play decisive roles. These aspects explain the much slower ISC in the para-isomer, as well as the non-negligible role of the El-Sayed forbidden pathway in the computed ISC dynamics.

Complete-active space self-consistent field↗
DOE OSTI2023

Temperature and nutrients drive eco-phenotypic dynamics in a microbial food web

Anthropogenic increases in temperature and nutrient loads will likely impact food web structure and stability. Although their independent effects have been reasonably well studied, their joint effects—particularly on coupled ecological and phenotypic dynamics—remain poorly understood. Here we experimentally manipulated temperature and nutrient levels in microbial food webs and used time-series analysis to quantify the strength of reciprocal effects between ecological and phenotypic dynamics across trophic levels. We found that (1) joint—often interactive—effects of temperature and nutrients on ecological dynamics are more common at higher trophic levels, (2) temperature and nutrients interact to shift the relative strength of top-down versus bottom-up control, and (3) rapid phenotypic change mediates observed ecological responses to changes in temperature and nutrients. Our results uncover how feedback between ecological and phenotypic dynamics mediate food web responses to environmental change. As a result, this suggests important but previously unknown ways that temperature and nutrients might jointly control the rapid eco-phenotypic feedback that determine food web dynamics in a changing world.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Data-Driven Dynamic Modeling in Power Systems: A Fresh Look on Inverter-Based Resource Modeling

Here, this article introduces ways to identify dynamic system models using measurement data. In power system analysis, a static model represents the time-invariant input-output relationship of a system, while a dynamic model describes the behavior of the system over time. For example, how will a system transit from one steady-state operation point to another?

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

A prototype scintillator real‐time beam monitor for ultra‐high dose rate radiotherapy

Background: FLASH Radiotherapy (RT) is an emergent cancer RT modality where an entire therapeutic dose is delivered at more than 1000 times higher dose rate than conventional RT. For clinical trials to be conducted safely, a precise and fast beam monitor that can generate out-of-tolerance beam interrupts is required. This paper describes the overall concept and provides results from a prototype ultra-fast, scintillator-based beam monitor for both proton and electron beam FLASH applications. Purpose: A FLASH Beam Scintillator Monitor (FBSM) is being developed that employs a novel proprietary scintillator material. The FBSM has capabilities that conventional RT detector technologies are unable to simultaneously provide: (1) large area coverage; (2) a low mass profile; (3) a linear response over a broad dynamic range; (4) radiation hardness; (5) real-time analysis to provide an IEC-compliant fast beam-interrupt signal based on true two-dimensional beam imaging, radiation dosimetry and excellent spatial resolution. Methods: The FBSM uses a proprietary low mass, less than 0.5 mm water equivalent, non-hygroscopic, radiation tolerant scintillator material (designated HM: hybrid material) that is viewed by high frame rate CMOS cameras. Folded optics using mirrors enable a thin monitor profile of ∼10 cm. A field programmable gate array (FPGA) data acquisition system generates real-time analysis on a time scale appropriate to the FLASH RT beam modality: 100–1000 Hz for pulsed electrons and 10–20 kHz for quasi-continuous scanning proton pencil beams. An ion beam monitor served as the initial development platform for this work and was tested in low energy heavy-ion beams ( 86 Kr +26 and protons). A prototype FBSM was fabricated and then tested in various radiation beams that included FLASH level dose per pulse electron beams, and a hospital RT clinic with electron beams. Results: Results presented in this report include image quality, response linearity, radiation hardness, spatial resolution, and real-time data processing. Furthermore, the HM scintillator was found to be highly radiation damage resistant. It exhibited a small 0.025%/kGy signal decrease from a 216 kGy cumulative dose resulting from continuous exposure for 15 min at a FLASH compatible dose rate of 237 Gy/s. Measurements of the signal amplitude versus beam fluence demonstrate linear response of the FBSM at FLASH compatible dose rates of >40 Gy/s. Comparison with commercial Gafchromic film indicates that the FBSM produces a high resolution 2D beam image and can reproduce a nearly identical beam profile, including primary beam tails. The spatial resolution was measured at 35–40 µm. Tests of the firmware beta version show successful operation at 20 000 Hz frame rate or 50 µs/frame, where the real-time analysis of the beam parameters is achieved in less than 1 µs. Conclusions: The FBSM is designed to provide real-time beam profile monitoring over a large active area without significantly degrading the beam quality. A prototype device has been staged in particle beams at currents of single particles up to FLASH level dose rates, using both continuous ion beams and pulsed electron beams. Using a novel scintillator, beam profiling has been demonstrated for currents extending from single particles to 10 nA currents. Radiation damage is minimal and even under FLASH conditions would require ≥50 kGy of accumulated exposure in a single spot to result in a 1% decrease in signal output. Beam imaging is comparable to radiochromic films, and provides immediate images without hours of processing. Real-time data processing, taking less than 50 µs (combined data transfer and analysis times), has been implemented in firmware for 20 kHz frame rates for continuous proton beams.

2D beam imaging↗
DOE OSTI2022

Random Walks and Sticky Surfaces: Single-Molecule Measurements of Solute Diffusion in Ethanol/Water-Filled Anodic Alumina Nanopores

Nanoporous anodic aluminum oxide (AAO) membranes are now being explored for use in advanced chemical separations, including in the dehydration of biofuels such as ethanol. Optimization of membrane performance requires an in-depth understanding of how solvent mixtures and solutes behave under nanoconfinement. In this work, the diffusion of rhodamine B (RhB) dye through 10 nm and 20 nm AAO nanopores filled with a series of ethanol/water mixtures (0 - 33% water) is explored by fluorescence correlation spectroscopy (FCS). RhB was found to diffuse through the pores by two distinct mechanisms with mean diffusion coefficients, $D_f$ and $D_s$, reflecting fast and slow diffusive motions, respectively, with values that differ by nearly two orders of magnitude. Further, both $D_f$ and $D_s$ increased with pore size and were significantly smaller than $D_b$, the RhB diffusion coefficient in bulk liquid. Mean $D_f$ values follow a composition-dependent trend that closely mimics the viscosity dependence of $D_b$. Additional slowing of fast RhB diffusion is attributed to both hydrodynamic drag and electrostatic interactions with the nanopore surface. The mean $D_s$ values exhibit a different trend with increasing water content, revealing an increase in $D_s$ and a decrease in the contributions of slow diffusion to the observed dynamics. The fluorescence time transient data used in the analysis show that the slow diffusion process is strongly hindered and likely involves frequent adsorption of RhB to the pore surfaces. These results provide new insights into the detailed molecular-level mechanisms of mass transport in nanoporous AAO membranes.

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

Non-equilibrium entropy production and information dissipation in a non-Markovian quantum dot

This study measures trajectory-level entropy production and information dissipation in a driven, non-Markovian quantum dot using time-resolved optical dynamics and machine-learning-based analysis. Although not a 2D-material system, it is relevant because it demonstrates quantitative extraction of nonequilibrium dynamics from nanoscale optical fluctuations, which is conceptually connected to the proposed studies of transient charge and spin dynamics at interfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Revisiting Power Systems Time-Domain Simulation Methods and Models

The changing nature of power systems dynamics is challenging present practices related to modeling and study of system-level dynamic behavior. While developing new techniques and models to handle the new modeling requirements, it is also critical to review some of the terminology used to describe existing simulation approaches and the embedded assumptions. This article provides a first-principles review of the simplifications and transformations commonly used in the formulation of time-domain simulation models. It introduces a taxonomy and classification of time-domain simulation models depending on their frequency bandwidth, network representation, and software availability. Furthermore, it focuses on the fundamental aspects of averaging techniques, and model reduction approaches that result in modeling choices, and discusses the associated challenges and opportunities of applying these methods in systems with large shares of Inverter Based Resources (IBRs). The article concludes with an illustrative simulation that compares the trajectories of an IBR-dominated system.

behavioral sciences↗
DOE OSTI2024

Retrieval Augmented Generation for Robust Cyber Defense

In cybersecurity, the ability to efficiently analyze and respond to vulnerabilities, weaknesses, attack patterns, and threat tactics is critical for effective defense strategies. With the increasing complexity and volume of cybersecurity data, traditional methods of querying and retrieving information are often inadequate. To address this challenge, we implemented Retrieval-Augmented Generation (RAG) systems—CyRAG and GraphCyRAG—that integrate large language models (LLMs) with both structured data from relational databases and knowledge graphs such as Neo4j. CyRAG is designed to handle structured data, focusing on CVE (Common Vulnerabilities and Exposures) and CWE (Common Weakness Enumeration) entities to generate accurate and context-rich responses. In contrast, GraphCyRAG leverages Neo4j knowledge graphs to retrieve interconnected information from CVE, CWE, CAPEC (Common Attack Pattern Enumeration and Classification), and ATT&CK (Adversarial Tactics, Techniques, and Common Knowledge) datasets. By utilizing Neo4j’s graph-based framework, GraphCyRAG enables deeper traversal of relationships between vulnerabilities and attack patterns, providing cybersecurity analysts with more comprehensive insights into potential attack vectors and mitigation strategies. Our preliminary results demonstrate that integrating knowledge graphs with RAG significantly enhances both the accuracy and depth of threat analysis, allowing for the retrieval of dynamic, real-time data and the generation of contextually aware responses. This approach helps analysts uncover hidden relationships between cyber entities, predict exploit paths, and prioritize mitigation efforts effectively. The integration of RAG with cybersecurity knowledge graphs represents a significant advancement in cybersecurity threat intelligence, enabling more informed decision-making and stronger defense strategies.

97 MATHEMATICS AND COMPUTING↗