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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Highly interactive, steered scientific workflows on HPC systems: optimizing design solutions

Scientific workflows are becoming increasingly important in high performance computing (HPC) settings, as the feasibility and appeal of many simultaneous heterogeneous tasks increases with increasing hardware capabilities. Currently no HPC-based workflow platform supports a dynamically adaptable workflow with interactive steering and analysis at run-time. Furthermore, for most workflow programs, compute resources are fixed for a given instance, resulting in a possible waste of expensive allocation resources when tasks are spawned and killed. Here we describe the design and testing of a run-time-interactive, adaptable, steered workflow tool capable of executing thousands of parallel tasks without an MPI programming model, using a database management system to facilitate task management through multiple live connections. We find that on the Oak Ridge Leadership Computing Facility pre-exascale Summit supercomputer it is possible to launch and interactively steer workflows with thousands of simultaneous tasks with negligible latency. For the case of particle simulation and analysis tasks that run for minutes to hours, this paradigm offers the prospect of a robust and efficient means to perform simulation-space exploration with on-the-fly analysis and adaptation.

Ossyra, John R↗
DOE OSTI2016

Atomistic details of protein dynamics and the role of hydration water

The importance of protein dynamics for their biological activity is nowwell recognized. Different experimental and computational techniques have been employed to study protein dynamics, hierarchy of different processes and the coupling between protein and hydration water dynamics. But, understanding the atomistic details of protein dynamics and the role of hydration water remains rather limited. Based on overview of neutron scattering, molecular dynamic simulations, NMR and dielectric spectroscopy results we present a general picture of protein dynamics covering time scales from faster than ps to microseconds and the influence of hydration water on different relaxation processes. Internal protein dynamics spread over a wide time range fromfaster than picosecond to longer than microseconds. We suggest that the structural relaxation in hydrated proteins appears on the microsecond time scale, while faster processes present mostly motion of side groups and some domains. Hydration water plays a crucial role in protein dynamics on all time scales. It controls the coupled protein-hydration water relaxation on 10 100 ps time scale. Our process defines the friction for slower protein dynamics. Analysis suggests that changes in amount of hydration water affect not only general friction, but also influence significantly the protein's energy landscape.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

On the mechanical behavior of a Nb 3 Sn superconducting coil during a quench: Two-dimensional finite element analysis of a quench heater protected magnet

New high-field accelerator magnets based on $Nb_3Sn$ superconductors are pushing the boundaries of magnet design and quench protection towards new limits. While their large stored energies and current densities result in a very challenging scenario for magnet protection, their great electromagnetic forces create also new requirements in terms of magnet design and stress management techniques. Furthermore, the strain sensitivity of $Nb_3Sn$ cables turns the electro-mechanical limits of the conductor into a parameter of the highest importance, where conductor degradation becomes a critical aspect in magnet operation. The coupling of all the above-mentioned considerations during quench is a case of special interest that adds further complexity to the design of $Nb_3Sn$ magnets. The objective of this paper is to provide a complete two-dimensional investigation of the coil and magnet structure mechanics during a quench event. Here, the analysis is performed using a combination of finite element codes that provide the necessary input for the mechanical study. The core of the modelling strategy relies on a thermal-electric model, whose results are directly used as loads in the mechanical simulation. In doing so, the stress evolution during and after a quench is obtained. We focus for this time in the analysis of a quench heater protected magnet, where dynamic effects are of less importance in contrast to other protection systems like CLIQ.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Revisiting the late-time growth of single-mode Rayleigh–Taylor instability and the role of vorticity

Growth of the single-fluid single-mode Rayleigh-Taylor instability (RTI) is revisited in 2D and 3D using fully compressible high-resolution simulations. We conduct a systematic analysis of the effects of perturbation Reynolds number (Rep) and Atwood number (A) on RTI’s late-time growth. Contrary to the common belief that single-mode RTI reaches a terminal bubble velocity, we show that the bubble re-accelerates when Rep is sufficiently large, consistent with [Ramaparabhu et al. 2006, Wei and Livescu 2012]. However, unlike in [Ramaparabhu et al. 2006], we find that for a sufficiently high Rep, the bubble’s late-time acceleration is persistent and does not vanish. Analysis of vorticity dynamics shows a clear correlation between vortices inside the bubble and re-acceleration. Due to symmetry around the bubble and spike (vertical) axes, the self-propagation velocity of vortices points in the vertical direction. If viscosity is sufficiently small, the vortices persist long enough to enter the bubble tip and accelerate the bubble [Wei and Livescu 2012]. A similar effect has also been observed in ablative RTI [Betti and Sanz 2006]. As the spike growth increases relative to that of the bubble at higher A, vorticity production shifts downward, away from the centerline and toward the spike tip. We modify the Betti-Sanz model for bubble velocity by introducing a vorticity efficiency factor η = 0.45 to accurately account for re-acceleration caused by vorticity in the bubble tip. It had been previously suggested that vorticity generation and the associated bubble re-acceleration are suppressed at high A. However, we present evidence that if the large Rep limit is taken first, bubble re-acceleration is still possible. Our results also show that re-acceleration is much easier to occur in 3D than 2D, requiring smaller Rep thresholds.

42 ENGINEERING↗
DOE OSTI2022

Superconvergence of Online Optimization for Model Predictive Control

We develop a one-Newton-step-per-horizon, online, lag-L, model predictive control (MPC) algorithm for solving discrete-time, equality-constrained, nonlinear dynamic programs. Based on recent sensitivity analysis results for the target problems class, we prove that the approach exhibits a behavior that we call superconvergence; that is, the tracking error with respect to the full horizon solution is not only stable for successive horizon shifts, but also decreases with increasing shift order to a minimum value that decays exponentially in the length of the receding horizon. The key analytical step is the decomposition of the one-step error recursion of our algorithm into algorithmic error and perturbation error. We show that the perturbation error decays exponentially with the lag between two consecutive receding horizons, while the algorithmic error, determined by Newton’s method, achieves quadratic convergence instead. Overall this approach induces our local exponential convergence result in terms of the receding horizon length for suitable values of L. In conclusion, numerical experiments validate our theoretical findings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

PMU-Based Monitoring and Mitigation of Delayed Voltage Recovery Using Admittances

This project analyzes the delayed voltage recovery phenomenon by simplifying the Western Electricity Coordinating Council (WECC) composite load model using the load admittance and the thermal relay dynamics. From this analysis, a closed form expression approximating the recovery time is derived and the key load parameters impacting the behavior of voltage recovery are identified. A monitoring scheme based on the measured load admittance is then proposed to detect the onset of stalling, even in the presence of voltage oscillations, and estimate the duration of the delayed voltage recovery. A mitigation scheme utilizing smart thermostats and offline learning is also derived to ensure that the voltage recovers to the pre-contingency voltage within a specified time. Both the monitoring and mitigation schemes only need local measurements at a substation making them promising for online applications. Results for the monitoring and mitigation schemes are described in detail for the IEEE 162 bus system validating the various assumptions used for the analysis and establishing the connection between the delayed voltage recovery phenomenon and load admittance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Syntrophic Co-Cultures of Clostridium Organisms to Produce Higher Alcohols & Other C6-C8 Metabolites (Final Report)

The goal of this project was to advance the systems biology understanding and predictive modeling of synthetic & syntrophic Clostridium microbial consortia, focusing on elucidation of metabolic networks and environmental signals in the consortia. The project has direct applicability to lignocellulosic-biomass based production of higher alcohols as advanced biofuels and C6-C8 metabolites, that can be used as chemicals or serve as biofuel precursors. This project studied and optimized three synthetic syntrophic systems. It examined the population dynamics using flow-cytometry, time-lapse microscopy and PCR analysis. 13C-based tracer analyses was used to examine the metabolite exchange between the syntrophic cell populations and the impact of those interactions on the transcriptome of the individual populations. To enhance our analytical and predictive capabilities, genome-scale models (GSMs) for these syntrophies was developed. RNAseq data for these syntrophic coculture systems were acquired to enable a molecular level understanding of the syntrophies aiming to identify the genetic networks of each organism in the co-culture and compare those against the networks of pure cultures.

09 BIOMASS FUELS↗
DOE OSTI2016

Identification of Metabolic Pathway Systems

The estimation of parameters in even moderately large biological systems is a significant challenge. This challenge is greatly exacerbated if the mathematical formats of appropriate process descriptions are unknown. To address this challenge, the method of dynamic flux estimation (DFE) was proposed for the analysis of metabolic time series data. Under ideal conditions, the first phase of DFE yields numerical representations of all fluxes within a metabolic pathway system, either as values at each time point or as plots against their substrates and modulators. However, this numerical result does not reveal the mathematical format of each flux. Thus, the second phase of DFE selects functional formats that are consistent with the numerical trends obtained from the first phase. While greatly facilitating metabolic data analysis, DFE is only directly applicable if the pathway system contains as many dependent variables as fluxes. Because most actual systems contain more fluxes than metabolite pools, this requirement is seldom satisfied. Auxiliary methods have been proposed to alleviate this issue, but they are not general. Here we propose strategies that extend DFE toward general, slightly underdetermined pathway systems.

Genetics & Heredity↗
DOE OSTI2023

Dynamic Network Analysis of Nuclear Science Literature for Research Influence Assessment

Analyzing nuclear science literature via data-driven methods is a critical step for assessing research influence and technology advancements. Indicators of scholarly activities may be buried in large volumes of nuclear research publications and collaboration networks over time. Mining for relevant scholarly influence trends in large volumes of text can be computationally challenging; however, open-source information on research collaborations over time can offer opportunities to extract meaningful insights. While network centrality analysis of scholarly research provides topology-based insights, additional emphasis on dynamics associated with the diffusion of information through these networks is important. Here this paper represents a step in that direction through the development of a novel dynamic network analysis framework and computational engine to identify key entities and capabilities over time within global scholarly nuclear science collaboration networks. Network theoretic, stochastic simulation, and optimization methods are leveraged to address variability in scholarly interactions, influence propagation, and collaboration patterns via network connections. A topic-aware influence maximization algorithm is developed to address the goal of identifying key influential authors in diverse research topics over time. Efficient parallelized implementation of the algorithm is applied to reduce computational costs. A proof-of-concept case study using open-source Scopus data with 33,517 published nuclear research papers from 2000-2019 is presented and representative analytic insights are generated. Broad implications of these insights are discussed and future research directions are also identified.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

Safety Benefits Assessment for Accident Tolerant Fuels in Consideration of Steam Generator Tube Degradation Using Dynamic Event Tree Analysis

Accident tolerant fuel (ATF) is expected to delay or prevent core damage by providing additional coping time under accidents involving loss of core cooling. The effect of extended coping time may vary depending on the plant response to accidents. Age-related component degradation that deteriorates plant performance over time could have an impact on the actual advantages of ATF. The potential safety benefits of two near-term ATF candidates, including Cr-coated Zr cladding and FeCrAl cladding, are assessed for a 2-in. loss-of-coolant accident with failed high-pressure safety injection using the dynamic event tree (DET) approach considering possible stress corrosion cracking of steam generator (SG) tubing under aging. The DET approach allows likelihood quantification of accident sequences leading to core damage, including stochastic variation of system response and human actions during accident mitigation. The safety benefits of the selected ATF claddings in terms of additional coping time and the core damage frequency reduction rate under specified accident situations were quantitatively estimated. The results show that the deployment of the two selected ATF claddings is expected to lead to longer coping times and lower core damage frequency due to the wider safety margin to peak cladding temperature they provide. The safety advantages would be greater as SG tube degradation proceeds. Thus, the two ATF candidates would lead to less severe consequences in terms of likelihood of core damage and susceptibility to the SG tube degradation than UO 2 -Zr fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Probability density adjoint for sensitivity analysis of the Mean of Chaos

Sensitivity analysis, especially adjoint based sensitivity analysis, is a powerful tool for engineering design which allows for the efficient computation of sensitivities with respect to many parameters. However, these methods break down when used to compute sensitivities of long-time averaged quantities in chaotic dynamical systems. This paper presents a new method for sensitivity analysis of ergodic chaotic dynamical systems, the density adjoint method. The method involves solving the governing equations for the system's invariant measure and its adjoint on the system's attractor manifold rather than in phase-space. This new approach is derived for and demonstrated on one-dimensional chaotic maps and the three-dimensional Lorenz system. It is found that the density adjoint computes very finely detailed adjoint distributions and accurate sensitivities, but suffers from large computational costs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Parallel replica dynamics method for bistable stochastic reaction networks: Simulation and sensitivity analysis

Stochastic reaction networks that exhibit bistable behavior are common in systems biology, materials science, and catalysis. Sampling of stationary distributions is crucial for understanding and characterizing the long-time dynamics of bistable stochastic dynamical systems. However, simulations are often hindered by the insufficient sampling of rare transitions between the two metastable regions. In this paper, we apply the parallel replica method for a continuous time Markov chain in order to improve sampling of the stationary distribution in bistable stochastic reaction networks. The proposed method uses parallel computing to accelerate the sampling of rare transitions. Furthermore, it can be combined with the path-space information bounds for parametric sensitivity analysis. With the proposed methodology, we study three bistable biological networks: the Schlögl model, the genetic switch network, and the enzymatic futile cycle network. We demonstrate the algorithmic speedup achieved in these numerical benchmarks. More significant acceleration is expected when multi-core or graphics processing unit computer architectures and programming tools such as CUDA are employed

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Multi-mode quasi-static excitation for systems with nonlinear joints

Finite element models can be used to model and predict the hysteresis and energy dissipation exhibited by nonlinear joints in structures. As a result of the nonlinearity, the frequency and damping of a mode is dependent on excitation amplitude, and when the modes remain uncoupled, quasi-static modal analysis has been shown to efficiently predict this behavior. However, in some cases the modes have been observed to couple such that the frequency and damping of one mode is dependent on the amplitude of other modes. To model the interactions between modes, one must integrate the dynamic equations in time, which is several orders of magnitude more expensive than quasi-static analysis. This work explores an alternative where quasi-static forces are applied in the shapes of two or more modes of vibration simultaneously, and the resulting load–displacement curves are used to deduce the effect of other modes on the effective frequency and damping of the mode in question. This methodology is demonstrated on a simple 2D cantilever beam structure with a single bolted joint which exhibits micro-slip nonlinearity over a range of vibration amplitudes. The predicted frequency and damping are compared with those extracted from a few expensive dynamic simulations of the structure, showing that the quasi-static approach produces reasonable albeit highly conservative bounds on the observed dynamics. This framework is also demonstrated on a 3D structure where dynamic simulations are infeasible.

42 ENGINEERING↗
DOE OSTI2016

V-T theory for the Self-Intermediate Scattering Function in a Monatomic Liquid

In V-T theory the atomic motion is harmonic vibrations in a liquid-specific potential energy valley, plus transits, which move the system rapidly among the multitude of such valleys. Here, in its first application to the self intermediate scattering function (SISF), V-T theory produced an accurate account of molecular dynamics (MD) data at all wave numbers q and time t. Recently, analysis of the mean square displacement (MSD) resolved a crossover behavior that was not observed in the SISF study. Our purpose here is to apply the more accurate MSD calibration to the SISF, and assess the results. We derive and discuss the theoretical equations for vibrational and transit contributions to the SISF. The time evolution is divided into three successive intervals: the vibrational interval when the vibrational contribution alone accurately accounts for the MD data; the crossover when the vibrational contribution saturates and the transit contribution becomes resolved; and the diffusive interval when the transit contribution alone accurately accounts for the MD data. Finally, the resulting theoretical error is extremely small at all q and t. V-T theory is compared to mode-coupling theories for the MSD and SISF, and to recent developments in Brownian motion experiments and theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Shock induced phase transition of water: Molecular dynamics investigation

Molecular dynamics simulations were carried out using numerous force potentials to investigate the shock induced phenomenon of pure bulk liquid water. Partial phase transition was observed at single shock velocity of 4.0 km/s without requirement of any external nucleators. Change in thermodynamic variables along with radial distribution function plots and spectral analysis revealed for the first time in the literature, within the context of molecular dynamic simulations, the thermodynamic pathway leading to formation of ice VII from liquid water on shock loading. The study also revealed information for the first time in the literature about the statistical time-frame after passage of shock in which ice VII formation can be observed and variations in degree of crystallinity of the sample over the entire simulation time of 100 ns.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Estimating Material Properties Under Extreme Conditions by Using Bayesian Model Calibration with Functional Outputs

Summary Dynamic material properties experiments provide access to the most extreme temperatures and pressures attainable in a laboratory setting; the data from these experiments are often used to improve our understanding of material models at these extreme conditions. We apply Bayesian model calibration to dynamic material property applications where the experimental output is a function: velocity over time. This framework can accommodate more uncertainties and facilitate analysis of new types of experiments relative to techniques traditionally used to analyse dynamic material experiments. However, implementation of Bayesian model calibration requires more sophisticated statistical techniques, because of the functional nature of the output as well as parameter and model discrepancy identifiability. We propose a novel Bayesian model calibration process to simplify and improve the estimation of the material property calibration parameters. Specifically, we propose scaling the likelihood function by an effective sample size rather than modelling the auto-correlation function to accommodate the functional output. Additionally, we propose sensitivity analyses by using the notion of 'modularization' to assess the effect of experiment-specific nuisance input parameters on estimates of the physical parameters. The Bayesian model calibration framework proposed is applied to dynamic compression of tantalum to extreme pressures, and we conclude that the procedure results in simple, fast and valid inferences on the material properties for tantalum.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A scale-wise analysis of intermittent momentum transport in dense canopy flows

We investigate the intermittent dynamics of momentum transport and its underlying time scales in the near-wall region of the neutrally stratified atmospheric boundary layer in the presence of a vegetation canopy. This is achieved through an empirical analysis of the persistence time scales (periods between successive zero-crossings) of momentum flux events, and their connection to the ejection–sweep cycle. Using high-frequency measurements from the GoAmazon campaign, spanning multiple heights within and above a dense canopy, the analysis suggests that, when the persistence time scales ( $t_p$ ) of momentum flux events from four different quadrants are separately normalized by $\varGamma _{w}$ (integral time scale of the vertical velocity), their distributions $P(t_p/\varGamma _{w})$ remain height-invariant. This result points to a persistent memory imposed by canopy-induced coherent structures, and to their role as an efficient momentum-transporting mechanism between the canopy airspace and the region immediately above. Moreover, $P(t_p/\varGamma _{w})$ exhibits a power-law scaling at times $t_{p}<\varGamma _{w}$ , with an exponential tail appearing for $t_{p} \geq \varGamma _{w}$ . By separating the flux events based on $t_p$ , we discover that around 80 % of the momentum is transported through the long-lived events ( $t_{p} \geq \varGamma _{w}$ ) at heights immediately above the canopy, while the short-lived ones ( $t_{p} < \varGamma _{w}$ ) only contribute marginally ( $\approx 20\,\%$ ). To explain the role of instantaneous flux amplitudes in momentum transport, we compare the measurements with newly developed surrogate data and establish that the range of time scales involved with amplitude variations in the fluxes tends to increase as one transitions from within to above the canopy.

Mechanics↗
DOE OSTI2023

Three-Dimensional Electromagnetic and Structural Analysis of Disruptions in the COMPASS Upgrade Vacuum Vessel

Structural analyses of disruptions and the electromagnetic (EM) analyses to support them typically use symmetry to produce a manageable model size. Nonaxisymmetric halo loads require larger analysis models. Here, in this paper, the results of transient EM analyses of halo strikes during a vertical displacement event are presented. The plasma motions and halo characteristics are prescribed based on analyses performed by the Institute of Plasma Physics of the Czech Republic. The time transient EM solution provides loads to a full three-dimensional transient dynamic analysis of the vacuum vessel. The responses to large lateral halo loads are altered and mitigated by the mode responses of the vessel. Dynamic load factors are computed for specific locations in the vessel and supports.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗