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 901 records · Page 50

Bifurcation and Stability Analysis of the Equilibrium States in Thermodynamic Systems in a Small Vicinity of the Equilibrium Values of Parameters

The dynamic behavior of thermodynamic system, described by one order parameter and one control parameter, in a small neighborhood of ordinary and bifurcation equilibrium values of the system parameters is studied. Using the general methods of investigating the branching (bifurcations) of solutions for nonlinear equations, we performed an exhaustive analysis of the order parameter dependences on the control parameter in a small vicinity of the equilibrium values of parameters, including the stability analysis of the equilibrium states, and the asymptotic behavior of the order parameter dependences on the control parameter (bifurcation diagrams). The peculiarities of the transition to an unstable state of the system are discussed, and the estimates of the transition time to the unstable state in the neighborhood of ordinary and bifurcation equilibrium values of parameters are given. The influence of an external field on the dynamic behavior of thermodynamic system is analyzed, and the peculiarities of the system dynamic behavior are discussed near the ordinary and bifurcation equilibrium values of parameters in the presence of external field. The dynamic process of magnetization of a ferromagnet is discussed by using the general methods of bifurcation and stability analysis presented in the paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Coccidioidomycosis Dynamics in Relation to Climate in the Southwestern United States: Cocci. and climate in the SW U.S.

Valley fever is endemic to the southwestern United States. Humans contract this fungal disease by inhaling spores of Coccidioides spp. Changes in the environment can influence the abundance and dispersal of Coccidioides spp., causing fluctuations in valley fever incidence. We combined county-level case records from state health agencies to create a regional valley fever database for the southwestern United States, including Arizona, California, Nevada, New Mexico, and Utah. We used this data set to explore how environmental factors influenced the spatial pattern and temporal dynamics of valley fever incidence during 2000–2015. We compiled climate and environmental geospatial data sets from multiple sources to compare with valley fever incidence. These variables included air temperature, precipitation, soil moisture, surface dust concentration, normalized difference vegetation index, and cropland area. We found that valley fever incidence was greater in areas with warmer air temperatures and drier soils. The mean annual cycle of incidence varied throughout the southwestern United States and peaked following periods of low precipitation and soil moisture. From year-to-year, however, autumn incidence was higher following cooler, wetter, and productive springs in the San Joaquin Valley of California. In southcentral Arizona, incidence increased significantly through time. By 2015, incidence in this region was more than double the rate in the San Joaquin Valley. Our analysis provides a framework for interpreting the influence of climate change on valley fever incidence dynamics. Our results may allow the U.S. Centers for Disease Control and Prevention to improve their estimates of the spatial pattern and intensity of valley fever endemicity.

54 ENVIRONMENTAL SCIENCES↗

Evaluating a Historical Airborne Release Test with Modern Modeling Methods

Safety analysts throughout the U.S. Department of Energy (DOE) complex rely on the data provided in the DOE Handbook, DOE-HDBK-3010, to determine source terms that may be incorporated into the document safety analyses. Most often, analysts simply take the bounding values due to time constraints or simply to bound calculations. This is a safe approach that helps avoid regulatory critique; however, it may not provide results that are meaningful or relevant to the conditions being evaluated. The derivation of the data, such as airborne release fractions (ARFs) and respirable fractions (RFs) in the Handbook often depend on very limited table-top and bench/laboratory experiments, as well as engineering judgment which may not be well substantiated, and may not be representative of the actual situation. One historical dataset included a test where a solid contaminant was sprinkled over a liquid gasoline pool and set on fire (also referred to as Mishima and Schwendiman's test later in this paper). The airborne release was measured over the duration of the burn, but then collected separately for a time after the burning was completed. It can be assumed that the collected contaminant in the latter period corresponds to re-suspended particles, whereas the earlier period collected contaminant includes a combination of directly released as well as early released re-suspended particles. This scenario was recently modeled using SIERRA/Fuego as part of a parametric sensitivity study, however there were some shortcomings of the analysis methods. First, the code at the time lacked a multi-component particle evaporation model, which is necessary to capture the particle size effect of liquid/solid particles that can dynamically evaporate the liquid component. Second, a particle sticking and resuspension model did not exist at the time of the study. Recent model implementation now permits the study be repeated with the improved physical models. In the earlier work, four potential mechanisms were considered to potentially contribute to the release of the particles. These were: 1. Evaporation Induced Entrainment (EIE); 2. Surface Agitation by Wind; 3. Surface Agitation by Boiling; 4. Residue Entrainment (Resuspension). The sensitivity study evaluated mechanism 1 and 3 in the above list. Parameters varied were selected to be a reflection of a practical range of test conditions. The varied parameters included the boiling duration, fuel height, injected mass, particle size, turbulence intensity, and initial particle height. The boiling duration was the most significant parameter. The most significant mechanism was no.3, Surface Agitation by Boiling. This paper reports the implementation of the new multi-component evaporation model and its effect on the findings from the previous effort. The Mishima and Schwendiman test was re-modeled with the new evaporation model, and the results of the re-evaluation are presented. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

An analysis of model proton-coupled electron transfer reactions via the mixed quantum-classical Liouville approach

The nonadiabatic dynamics of model proton-coupled electron transfer (PCET) reactions is investigated for the first time using a surface-hopping algorithm based on the solution of the mixed quantum-classical Liouville equation (QCLE). This method provides a rigorous treatment of quantum coherence/decoherence effects in the dynamics of mixed quantum-classical systems, which is lacking in the molecular dynamics with quantum transitions surface-hopping approach commonly used for simulating PCET reactions. Within this approach, the protonic and electronic coordinates are treated quantum mechanically and the solvent coordinate evolves classically on both single adiabatic surfaces and on coherently coupled pairs of adiabatic surfaces. Both concerted and sequential PCET reactions are studied in detail under various subsystem-bath coupling conditions and insights into the dynamical principles underlying PCET reactions are gained. Notably, an examination of the trajectories reveals that the system spends the majority of its time on the average of two coherently coupled adiabatic surfaces, during which a phase enters into the calculation of an observable. In general, the results of this paper demonstrate the applicability of QCLE-based surface-hopping dynamics to the study of PCET and emphasize the importance of mean surface evolution and decoherence effects in the calculation of PCET rate constants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ergodic theory and visualization. II. Fourier mesochronic plots visualize (quasi)periodic sets

We present an application and analysis of a visualization method for measure-preserving dynamical systems introduced by I. Mezić and A. Banaszuk [Physica D 197, 101 (2004)], based on frequency analysis and Koopman operator theory. This extends our earlier work on visualization of ergodic partition [Z. Levnajić and I. Mezić, Chaos 20, 033114 (2010)]. Our method employs the concept of Fourier time average [I. Mezić and A. Banaszuk, Physica D 197, 101 (2004)], and is realized as a computational algorithms for visualization of periodic and quasi-periodic sets in the phase space. The complement of periodic phase space partition contains chaotic zone, and we show how to identify it. The range of method's applicability is illustrated using well-known Chirikov standard map, while its potential in illuminating higher-dimensional dynamics is presented by studying the Froeschlé map and the Extended Standard Map.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Time-dependent behavioral data from zebrafish reveals novel signatures of chemical toxicity using point of departure analysis

High-content imaging of embryonic zebrafish behavior can be used as a screening approach to rapidly evaluate the relative potential for chemicals to cause toxicity. However, most statistical methods applied to these data so far transform movement values to incidence-based “hits” and calculate lowest effect levels (LELs), which loses individual fish resolution of behavior and defies hazard ranking due to reliance on applied dose levels. Here, we have developed and parallelized a workflow to calculate benchmark dose (BMD) values from dynamic, high-content zebrafish behavior data that scales for high-throughput chemical screening. To appropriately capture the dynamics of zebrafish response from light to dark stimulus, we summarized time-dependent data using both area under the curve and immediate change at the transition point into two novel metrics that characterized abnormal behavior as a function of chemical concentration. The BMD workflow was applied to calculate BMD 10 values of 1,060 ToxCast chemicals for 24 zebrafish endpoints, including behavior, mortality and morphology. The BMD 10 values provided better precision and separation than LELs for clustering chemicals since they were derived from models that best-fit their concentration-response curves. Clustering analysis revealed behavioral signatures as the most sensitive endpoints. High concordance in chemical activity was observed between ToxCast assays and zebrafish behavioral endpoints, however ToxPi analysis indicated that rankings based on in vitro data was not a reliable predictor of in vivo rankings for lower potency chemicals. The methods developed in this work will enable the use of high-content zebrafish behavioral screening data for BMD analysis in toxicological hazard assessment.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling Framework for Bulk Electric Grid Impacts from HEMP E1 and E3 Effects (Tasks 3.1 Final Report)

This report presents a framework to evaluate the impact of a high-altitude electromagnetic pulse (HEMP) event on a bulk electric power grid. This report limits itself to modeling the impact of EMP E1 and E3 components. The co-simulation of E1 and E3 is presented in detail, and the focus of the paper is on the framework rather than actual results. This approach is highly conservative as E1 and E3 are not maximized with the same event characteristics and may only slightly overlap. The actual results shown in this report are based on a synthetic grid with synthetic data and a limited exemplary EMP model. The framework presented can be leveraged and used to analyze the impact of other threat scenarios, both manmade and natural disasters. This report d escribes a Monte-Carlo based methodology to probabilistically quantify the transient response of the power grid to a HEMP event. The approach uses multiple fundamental steps to characterize the system response to HEMP events, focused on the E1 and E3 components of the event. 1) Obtain component failure data related to HEMP events testing of components and creating component failure models. Use the component failure model to create component failure conditional probability density function (PDF) that is a function of the HEMP induced terminal voltage. 2) Model HEMP scenarios and calculate the E1 coupled voltage profiles seen by all system components. Model the same HEMP scenarios and calculate the transformer reactive power consumption profiles due to E3. 3) Sample each component failure PDF to determine which grid components will fail, due to the E1 voltage spike, for each scenario. 4) Perform dynamic simulations that incorporate the predicted component failures from E1 and reactive power consumption at each transformer affected by E3. These simulations allow for secondary transients to affect the relays/protection remaining in service which can lead to cascading outages. 5) Identify the locations and amount of load lost for each scenario through grid dynamic simulation. This can be an indication of the immediate grid impacts from a HEMP event. In addition, perform more detailed analysis to determine critical nodes and system trends. 6) To help realize the longer-term impacts, a security constrained alternating current optimal power flow (ACOPF) is run to maximize critical load served. This report describes a modeling framework to assess the systemic grid impacts due to a HEMP event. This stochastic simulation framework generates a large amount of data for each Monte Carlo replication, including HEMP location and characteristics, relay and component failures, E3 GIC profiles, cascading dynamics including voltage and frequency over time, and final system state. This data can then be analyzed to identify trends, e.g., unique system behavior modes or critical components whose failure is more likely to cause serious systemic effects. The proposed analysis process is demonstrated on a representative system. In order to draw realistic conclusions of the impact of a HEMP event on the grid, a significant amount of work remains with respect to modeling the impact on various grid components.

24 POWER TRANSMISSION AND DISTRIBUTION↗

High Performance Interactive System Dynamics Visualization

This brochure describes a system dynamics simulation (SD) framework that supports an end-to-end analysis workflow that is optimized for deployment on ESIF facilities(Peregrine and the Insight Center). It includes (I) parallel and distributed simulation of SD models, (ii) real-time 3D visualization of running simulations, and (iii) comprehensive database-oriented persistence of simulation metadata, inputs, and outputs.

97 MATHEMATICS AND COMPUTING↗

High Performance Interactive System Dynamics Visualization

This presentation describes a system dynamics simulation (SD) framework that supports an end-to-end analysis workflow that is optimized for deployment on ESIF facilities(Peregrine and the Insight Center). It includes (I) parallel and distributed simulation of SD models, (ii) real-time 3D visualization of running simulations, and (iii) comprehensive database-oriented persistence of simulation metadata, inputs, and outputs.

97 MATHEMATICS AND COMPUTING↗

Comparison of damage and ablation dynamics of multilayer dielectric films initiated by few-cycle pulses versus longer femtosecond pulses

The importance of high intensity few- to single-cycle laser pulses for applications such as intense isolated attosecond pulse generation is constantly growing, and with the breakdown of the monochromatic approximation in field ionization models, the few-cycle pulse (FCP) interaction with solids near the damage threshold has ushered a new paradigm of nonperturbative light–matter interaction. In this Letter, we systematically study and contrast how femtosecond laser-induced damage and ablation behaviors of SiO 2 /HfO 2 -based reflective multilayer dielectric thin film systems vary between FCP and 110 fs pulses. With time-resolved surface microscopy and ex situ analysis, we show that there are distinct differences in the interaction depending on the pulse duration, specifically in the “blister” morphology formation at lower fluences (damage) as well as in the dynamics of debris formation at higher fluences (ablation).

42 ENGINEERING↗

An Overview of the Monte Carlo Application ToolKit (MCATK)

MCATK is a C++ component-based Monte Carlo neutron-gamma transport software library designed to build specialized applications and designed to provide new functionality in existing general-purpose Monte Carlo codes like MCNP; it was developed with Agile software engineering methodologies under the motivation to reduce costs. The characteristics of MCATK can be summarized as follows: MCATK physics – continuous energy neutron-gamma transport with multi-temperature treatment, static eigenvalue (k and α) algorithms, time-dependent algorithm, fission chain algorithms; MCATK geometry – mesh geometries, solid body geometries. MCATK provides verified, unit-tested Monte Carlo components, flexibility in Monte Carlo applications development, and numerous tools such as geometry and cross section plotters. Recent work has involved deterministic and Monte Carlo analysis of stochastic systems. Static and dynamic analysis is discussed, and the results of a dynamic test problem are given.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Anisotropic structural dynamics of monolayer crystals revealed by femtosecond surface X-ray scattering

Ultrafast X-ray scattering is one of the primary tools to track intrinsic crystallographic evolution with atomic accuracy in real time. However, its application to study nonequilibrium structural properties at the two-dimensional limit remains a long-standing challenge due to a significant reduction of diffraction volume and complexity of data analysis. Here, we report femtosecond surface X-ray diffraction in combination with crystallographic model-refinement calculations to quantify the ultrafast structural dynamics of monolayer WSe 2 crystals supported on a substrate. We found the absorbed optical photon energy is preferably coupled to the in-plane lattice vibrations within one picosecond whereas the out-of-plane lattice vibration amplitude remains unchanged during the first ten picoseconds. The model-assisted fitting suggests an asymmetric intralayer spacing change upon excitation. The observed nonequilibrium anisotropic structural dynamics agrees with first-principles modelling in both real and momentum space, marking the distinct structural dynamics of monolayer crystals from their bulk counterparts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

XFEL SASE pulses can enhance time-dependent observables

Abstract X-ray free electron lasers (XFELs) have emerged as powerful sources of short and intense x-ray pulses. We propose a simple and robust procedure which takes advantage of the inherent stochasticity of self-amplified stimulated emission (SASE) pulses to enhance the time-resolution and signal strength of the recorded data. Notably, the proposed method is able to enhance the average signal without knowledge of the signal strength of individual shots. Simple metrics for the probe pulses are introduced, such as an effective pulse duration applicable to SASE pulses characterised in the time domain using e.g. an X-band transverse cavity. The approach is evaluated using simulated and real pulse data in the context of ultrafast electron dynamics in a molecule. Utilising H 2 as a model system, we demonstrate the efficacy of the method theoretically, successfully enhancing the predicted nonresonant ultrafast x-ray scattering signal associated with electron dynamics. The method presented is broadly applicable and offers a general strategy for enhancing time-dependent observables at XFELs.

data analysis↗

Kinetic theory-based numerical modeling and analysis of bi-disperse segregated mixture fluidized bed

We discuss a series of continuum Euler-Euler simulations of an initially mixed bi-disperse fluidized bed which segregates under certain operating conditions. The simulations use the multi-phase kinetic theory-based description of the momentum and energy exchanges between the phases by Simonin’s Group [see e.g. Gourdel, Simonin and Brunier (1999). Proceedings of 6th International Conference on Circulating Fluidized Beds, Germany, pp. 205-210]. The discussion and analysis of the results focus on the fluid-particle momentum exchange (i.e. drag). Simulations using mono- and poly-disperse fluid-particle drag correlations are analyzed for the Geldart D-type size bi-disperse gas-solid experiments performed by Goldschmidt et al. [Powder Tech., pp. 135-159 (2003)]. The poly-disperse gas-particle drag correlations account for the local particle size distribution by using an effective mixture diameter when calculating the Reynolds number and then correcting the resulting force coefficient. Simulation results show very good predictions of the segregation index for bidisperse beds with the mono-disperse drag correlations contrary to the poly-disperse drag correlations for which the segregation rate is systematically under-predicted. The statistical analysis of the results shows a clear separation in the distribution of the gas-particle mean relaxation times of the small and large particles with simulations using the mono-disperse drag. In contrast, the poly-disperse drag simulations have a significant overlap and also a smaller difference in the mean particle relaxation times. This results in the small and large particles in the bed to respond to the gas similarly without enough relative time lag. The results suggest that the difference in the particle response time induce flow dynamics favorable to a force imbalance which results in the segregation.

97 MATHEMATICS AND COMPUTING↗

Molecular dynamics simulation of gold cluster growth during sputter deposition

We present a molecular dynamics simulation scheme that we apply to study the time evolution of the self-organized growth process of metal cluster assemblies formed by sputter-deposited gold atoms on a planar surface. The simulation model incorporates the characteristics of the plasma-assisted deposition process and allows for an investigation over a wide range of deposition parameters. It is used to obtain data for the cluster properties which can directly be compared with recently published experimental data for gold on polystyrene [M. Schwartzkopf et al., ACS Appl. Mater. Interfaces 7, 13547 (2015)]. While good agreement is found between the two, the simulations additionally provide valuable time-dependent real-space data of the surface morphology, some of whose details are hidden in the reciprocal-space scattering images that were used for the experimental analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Translation-rotation decoupling and nonexponentiality in room temperature ionic liquids

Using a combination of light scattering techniques and broadband dielectric spectroscopy, we have measured the temperature dependence of structural relaxation time and self diffusion in three imidazolium-based room temperature ionic liquids: [bmim][NTf2], [bmim][PF6], and [bmim][TFA]. A detailed analysis of the results demonstrates that self diffusion decouples from structural relaxation in these systems as the temperature is decreased toward Tg. The degree to which the dynamics are decoupled, however, is shown to be surprisingly weak when compared to other supercooled liquids of similar fragility. In addition to the weak decoupling, we demonstrate that the temperature dependence of the structural relaxation time in all three liquids can be well described by a single Vogel-Fulcher-Tamann function over 13 decades in time from 10 11 s up to 102 s. Furthermore, the stretching of the structural relaxation is shown to be temperature independent over the same range of time scales, i.e., time temperature superposition is valid for these ionic liquids from far above the melting point down to the glass transition temperature.We suggest that these phenomena are interconnected and all result from the same underlying mechanism strong and directional intermolecular interactions.

Griffin, Phillip↗

Cooperative macromolecular device revealed by meta-analysis of static and time-resolved structures

Here we present a meta-analysis of a large collection of static structures of a protein in the Protein Data Bank in order to extract the progression of structural events during protein function. We apply this strategy to the homodimeric hemoglobin HbI from Scapharca inaequivalvis. We derive a simple dynamic model describing how binding of the first ligand in one of the two chemically identical subunits facilitates a second binding event in the other partner subunit. The results of our ultrafast time-resolved crystallographic studies support this model. We demonstrate that HbI functions like a homodimeric mechanical device, such as pliers or scissors. Ligand-induced motion originating in one subunit is transmitted to the other via conserved pivot points, where the E and F' helices from two partner subunits are 'bolted' together to form a stable dimer interface permitting slight relative rotation but preventing sliding.

59 BASIC BIOLOGICAL SCIENCES↗

Detonation mode and frequency analysis under high loss conditions for stoichiometric propane-oxygen

In this paper, the propagation characteristics of galloping detonations were quantified with a high-time-resolution velocity diagnostic. Combustion waves were initiated in 30-m lengths of 4.1-mm inner diameter transparent tubing filled with stoichiometric propane–oxygen mixtures. Chemiluminescence from the resulting waves was imaged to determine the luminous wave front position and velocity every 83.3 μ. As the mixture initial pressure was decreased from 20 to 7 kPa, the wave was observed to become increasingly unsteady and transition from steady detonation to a galloping detonation. While wave velocities averaged over the full tube length smoothly decreased with initial pressure down to half of the Chapman–Jouguet detonation velocity (D CJ ) at the quenching limit, the actual propagation mechanism was seen to be a galloping wave with a cycle period of approximately 1.0 ms, corresponding to a cycle length of 1.3–2.0 m or 317–488 tube diameters depending on the average wave speed. The long test section length of 7300 tube diameters allowed observation of up to 20 galloping cycles, allowing for statistical analysis of the wave dynamics. In the galloping regime, a bimodal velocity distribution was observed with peaks centered near 0.4 D CJ and 0.95 D CJ . Decreasing initial pressure increasingly favored the low velocity mode. Galloping frequencies ranged from 0.8 to 1.0 kHz and were insensitive to initial mixture pressure. Wave deflagration-to-detonation transition and detonation failure trajectories were found to be repeatable in a given test and also across different initial mixture pressures. The temporal duration of wave dwell at the low and high velocity modes during galloping was also quantified. It was found that the mean wave dwell duration in the low velocity mode was a weak function of initial mixture pressure, while the mean dwell time in the high velocity mode depended exponentially on initial mixture pressure. Analysis of the velocity histories using dynamical systems ideas demonstrated trajectories that varied from stable to limit cycles to aperiodic motion with decreasing initial pressure. Finally, the results indicate that galloping detonation is a persistent phenomenon at long tube lengths.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗