Search NASA⌕ Search

SEARCH · Search NASA

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 1,387 records · Page 77

Analysis of the elliptic integrable non-linear system in IOTA using tracking of a single electron

Integrable nonlinear lattices that can be realized in practical accelerators are of great interest, as they offer the potential to support high-intensity beams via Landau damping of collective instabilities. One such system, based on an elliptic potential, has been extensively studied at the IOTA storage ring at Fermilab. The analysis of strongly nonlinear dynamics with multi-particle bunches is challenging due to the rapid decoherence of kicked beams. IOTA has the capability to track single electrons using linear multi-anode photomultiplier tubes for simultaneously measuring transverse coordinates and arrival times of synchrotron-radiation pulses. This technology enables the full reconstruction of turn-by-turn positions and momenta in all three planes for a single particle. Using this apparatus, we measured the dependence of small-amplitude tunes on the strength of the nonlinear magnet, as well as tunes dependence on oscillations amplitudes.

Romanov, A. [Fermilab]↗

Offline Power Systems Applications Enabled by Phasor Measurement Units: Technical Assistance to the Power Sectors of Southeast Asia

This report provides a brief overview of several offline (non-real-time) applications facilitated by high-resolution time-synchronized measurements recorded by phasor measurement units (PMUs). The high reporting rate and time-synchronization of PMU records provide a detailed view of power system dynamics, enabling electric utilities to obtain a better understanding of their systems. In this report, the following applications have been reviewed: Power plant model validation, System model validation, Ringdown oscillation analysis, Frequency response analysis, Postmortem analysis of disturbance events Along with a brief technical background of the applications above, applicable North American Electric Reliability Corporation (NERC) standards have been discussed, and examples of implementation in North American organizations have been provided. Implementing several of the discussed applications may need a preliminary stage of data gathering from multiple entities, and several frameworks and process flows have been formulated by organizations around the world for this purpose. However, the data-gathering stage has not been considered in the scope of the present report.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Journey to Time-Variable Moment Tensors through Inversion of Acoustic and Seismoacoustic Data

We explore the capability of acoustic and seismoacoustic datasets to directly resolve a complex, time-variable source consisting of a buried mechanism, represented as a moment tensor, and a spall mechanism, represented as a vertical force at the surface. Traditionally, each component of a resolved moment tensor assumes one underlying source time function, which likely fails to capture the full evolution of a dynamic source, such as an explosion followed by slip on near-source joints or development of spallation. Specifically, we expand previous work to resolve a time-variable moment tensor using single-modality and joint-modality inversion frameworks through analysis of infrasound and seismoacoustic data recorded as part of the Source Physics Experiment Phase II: Dry Alluvium Geology (DAG). We investigate the impact of including signals from seismic-to-air coupling that are local to each infrasound sensor in comparison to mainly atmosphere-propagating acoustic signals, which occur from coupling of the wavefield from the subsurface to the atmosphere directly above the source. Additionally, we assess the ability of our inversion algorithm to fit observed infrasound data using a variety of time-variable source mechanisms. First, we consider the buried moment tensor source alone, which assumes that the determined Green’s functions incorporate effects from spallation or that the impact from spallation is minimal. Second, we examine the estimated buried moment tensor and vertical surface spallation as terms that must both be resolved in the inversion. Third, we assess the ability for an estimated vertical surface spallation source to fit the acoustic data on its own. Finally, we compare results from the joint inversion of both seismic geophone and infrasound acoustic data for the buried-only source compared to buried and spallation sources. Our results are a preliminary investigation into the applications of the inversion technique to recorded datasets and show the technique has limited capabilities using acoustic data alone. Instead, this method shows promise for seismic and seismoacoustic datasets to resolve the time-variable mechanisms of a buried source.

47 OTHER INSTRUMENTATION↗

Analysis of the Elliptic Integrable Non-Linear System in IOTA Using Tracking of a Single Electron

Integrable nonlinear lattices that can be realized in practical accelerators are of great interest, as they offer the potential to support high-intensity beams via Landau damping of collective instabilities. One such system, based on an elliptic potential, has been extensively studied at the IOTA storage ring at Fermilab. The analysis of strongly nonlinear dynamics with multi-particle bunches is challenging due to the rapid decoherence of kicked beams. IOTA has the capability to track single electrons using linear multi-anode photomultiplier tubes for simultaneously measuring transverse coordinates and arrival times of synchrotron-radiation pulses. This technology enables the full reconstruction of turn-by-turn positions and momenta in all three planes for a single particle. Using this apparatus, we measured the dependence of small-amplitude tunes on the strength of the nonlinear magnet, as well as tunes dependence on oscillations amplitudes.

Romanov, Aleksandr Leonidovich [Fermilab] (ORCID:0↗

Measurement of Greenhouse Gas Flux from Agricultural Soils Using Static Chambers

Measurement of greenhouse gas (GHG) fluxes between the soil and the atmosphere, in both managed and unmanaged ecosystems, is critical to understanding the biogeochemical drivers of climate change and to the development and evaluation of GHG mitigation strategies based on modulation of landscape management practices. The static chamber-based method described here is based on trapping gases emitted from the soil surface within a chamber and collecting samples from the chamber headspace at regular intervals for analysis by gas chromatography. Change in gas concentration over time is used to calculate flux. This method can be utilized to measure landscape-based flux of carbon dioxide, nitrous oxide, and methane, and to estimate differences between treatments or explore system dynamics over seasons or years. Infrastructure requirements are modest, but a comprehensive experimental design is essential. This method is easily deployed in the field, conforms to established guidelines, and produces data suitable to large-scale GHG emissions studies.

54 ENVIRONMENTAL SCIENCES↗

diffReplication - An Energy-Aware Fault Tolerance Model for Silent Error Detection and Mitigation in Heterogeneous Extreme-scale Computing Environment

At extreme scale, the frequency of silent errors – a class of errors that remain undetected by low-level error detection mechanisms – increases significantly with the computational complexity of the application and the scale of the computing infrastructure. As hardware and software advances are made to usher in the next scientific era of computing, developing new approaches to mitigate the impact of silent errors remains a challenging problem. In this work, we propose an energy-aware fault-tolerance model, referred to diffReplication to overcome silent errors. In the proposed model, the main process is associated with one replica that executes at the same rate as the main process, and one diffReplica that is executed at a fraction of the main process' execution rate. If the main and its replica reach consensus at the end of a computation phase, the state of the diffReplica is updated and computation is resumed. If the synchronization attempt results in a disagreement, however, the diffReplica increases its execution speed to complete the computation and quickly reach the synchronization barrier. Assuming a single error over any given synchronization interval, a majority voting is used to reach consensus and tolerate silent errors. To further enhance its performance, diffReplication is augmented with speculative execution, whereby the main or its fast replica is selected to continue execution without waiting for the diffReplica. The selection process is based on the previous behaviour of the main and its replica. A performance analysis study is carried out to assess the performance of diffReplication, in terms of the energy saving and time-to-completion reduction achieved by the diffReplication scheme. The experiment shows that speculative execution reduces the time to completion with additional energy, and dynamic decision-making balances the energy consumption and time to completion.

97 MATHEMATICS AND COMPUTING↗

Transient dynamics of terrestrial carbon storage: Mathematical foundation and numeric examples

Terrestrial ecosystems absorb roughly 30% of anthropogenic CO 2 emissions since preindustrial era, but it is unclear whether this carbon (C) sink will endure into the future. Despite extensive modeling, experimental, and observational studies, what fundamentally determines transient dynamics of terrestrial C storage under climate change is still not very clear. Here we develop a new framework for understanding transient dynamics of terrestrial C storage through mathematical analysis and numerical experiments. Our analysis indicates that the ultimate force driving ecosystem C storage change is the C storage capacity, which is jointly determined by ecosystem C input (e.g., net primary production, NPP) and residence time. Since both C input and residence time vary with time, the C storage capacity is time-dependent and acts as a moving attractor that actual C storage chases. The rate of change in C storage is proportional to the C storage potential, the difference between the current storage and the storage capacity. The C storage capacity represents instantaneous responses of the land C cycle to external forcing, whereas the C storage potential represents the internal capability of the land C cycle to influence the C change trajectory in the next time step. The influence happens through redistribution of net C pool changes in a network of pools with different residence times. Furthermore, this and our other studies have demonstrated that one matrix equation can exactly replicate simulations of most land C cycle models (i.e., physical emulators). As a result, simulation outputs of those models can be placed into a three-dimensional (3D) parameter space to measure their differences. The latter can be decomposed into traceable components to track the origins of model uncertainty. Moreover, the emulators make data assimilation computationally feasible so that both C flux- and pool-related datasets can be used to better constrain model predictions of land C sequestration. We also propose that the C storage potential be the targeted variable for research, market trading, and government negotiation for C credits.

54 ENVIRONMENTAL SCIENCES↗

Unstructured mesh-based multi-physics reactor transient simulation with iMC code

In this study, we present a time-dependent Monte Carlo simulation result about a complex fuel-loaded light water reactor system with thermal-hydraulic feedback using the KAIST iMC code. The intra-pin temperature distribution is calculated based on Monte Carlo-tallied detailed power distribution and the finite element heat transfer method. The temperature effect of fuel, absorber, and coolant are adequately reflected into the dynamic response by the on-the-fly Doppler broadened cross-section correction. For the time-dependent Monte Carlo transport simulation, the predictor-corrector quasi-static Monte Carlo method is used for this study. We tested the coupled analysis framework with a 17-by-17 fuel assembly loaded with the centrally-shielded burnable absorber (CSBA) fuel element. The temperature feedback effect on the material cross-section and the moderator density incurred negative reactivity in response to the external positive reactivity insertion, showing the inherent stability of the reactor system. The pin-wise power and temperature change during the system transient are presented and discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Hydrogenation of the nanopowders that form in a carbon-helium plasma stream during the introduction of Ni and Mg

Composite nanoparticles consisting of magnesium, nickel, and carbon atoms are studied both theoretically and experimentally. The calculations performed in terms of the density functional theory show that the jump frequency of hydrogen atoms in nickel-containing magnesium hydride increases substantially near impurity nickel atoms; as a result, the rate of hydrogen absorption by magnesium also increases. Nickel on the magnesium surface is shown to be absorbed via an island growth mechanism. Composite Mg-C, Ni-C, and Mg-Ni-C powders are produced by plasmachemical synthesis in a carbon-helium plasma stream. Hydrogen is introduced into a chamber during synthesis. It is found by X-ray photoelectron spectroscopy and thermogravimetric analysis that, among these three composites, only Mg-Ni-C contains magnesium fixed in the MgH{sub 2} compound. The process of such 'ultrarapid' hydrogenation of magnesium, which occurs in the time of formation of composite nanoparticles, can be explained by the catalytic action of nickel, which is enhanced by a high temperature. Scanning electron microscopy micrographs demonstrate the dynamics of the dehydrogenation of Mg-Ni-C composite nanoparticles in heating by an electron beam.

36 MATERIALS SCIENCE↗

Data-driven material models for atomistic simulation

The central approximation made in classical molecular dynamics simulation of materials is the interatomic potential used to calculate the forces on the atoms. Great effort and ingenuity is required to construct viable functional forms and find accurate parametrizations for potentials using traditional approaches. Machine learning has emerged as an effective alternative approach to develop accurate and robust interatomic potentials. Beginning with a very general model form, the potential is learned directly from a database of electronic structure calculations and therefore can be viewed as a multiscale link between quantum and classical atomistic simulations. Risk of inaccurate extrapolation exists outside the narrow range of time and length scales where the two methods can be directly compared. Here, we use the spectral neighbor analysis potential (SNAP) and show how a fit can be produced with minimal interpolation errors which is also robust in extrapolating beyond training. To demonstrate the approach, we have developed a tungsten-beryllium potential suitable for the full range of binary compositions. Subsequently, large-scale molecular dynamics simulations were performed of high energy Be atom implantation onto the (001) surface of solid tungsten. The machine learned W-Be potential generates a population of implantation structures consistent with quantum calculations of defect formation energies. A very shallow ( < 2 nm ) average Be implantation depth is predicted which may explain ITER diverter degradation in the presence of beryllium.

36 MATERIALS SCIENCE↗

TH-C-17A-03: Dynamic Visualization and Dosimetry of IMRT and VMAT Treatment Plans by Video-Rate Imaging of Cherenkov Radiation in Pure Water

Purpose: A novel optical dosimetry technique for the QA and verification of intensity-modulated radiation therapy (IMRT) and volumetric-modulated arc therapy (VMAT) radiotherapy plans was investigated for the first time by capturing images of the induced Cherenkov radiation in water. Methods: An intensified CCD camera (ICCD) was used to acquire a two-dimensional (2D) projection image of the Cherenkov radiation induced by IMRT and VMAT plans, based on the Task Group 119 C-Shape geometry. Plans were generated using the Varian Eclipse treatment planning system (TPS) and delivered using 6 MV x-rays from a Varian TrueBeam Linear Accelerator (Linac) incident on a water tank. The ICCD acquisition was gated to the Linac, operated for single pulse imaging, and binned to a resolution of 512×512 pixels. The resulting videos were analyzed temporally for regions of interest (ROI) covering the planning target volume (PTV) and organ at risk (OAR) and summed to obtain an overall light distribution, which was compared to the expected dose distribution from the TPS using a gammaindex analysis. Results: The chosen camera settings resulted in data at 23.5 frames per second. Temporal intensity plots of the PTV and OAR ROIs confirmed the preferential delivery of dose to the PTV versus the OAR, and the gamma analysis yielded 95.2% and 95.6% agreement between the light distribution and expected TPS dose distribution based upon a 3% / 3 mm dose difference and distance-to-agreement criterion for the IMRT and VMAT plans respectively. Conclusion: The results from this initial study demonstrate the first documented use of Cherenkov radiation for optical dosimetry of dynamic IMRT and VMAT treatment plans. The proposed modality has several potential advantages over alternative methods including the real-time nature of the acquisition, and upon future refinement may prove to be a robust and novel dosimetry method with both research and clinical applications. NIH R01CA109558 and R21EB017559.

07 ISOTOPE AND RADIATION SOURCES↗

Gravitational waves from asymmetric oscillon dynamics?

It has been recently suggested that oscillons produced in the early universe from certain asymmetric potentials continue to emit gravitational waves for a number of e-folds of expansion after their formation, leading to potentially detectable gravitational wave signals. We revisit this claim by conducting a convergence study using graphics processing unit (GPU)-accelerated lattice simulations and show that numerical errors accumulated with time are significant in low-resolution scenarios, or in scenarios where the run-time causes the resolution to drop below the relevant scales in the problem. Our study determines that the dominant, growing high frequency peak of the gravitational wave signals in the fiducial "hill-top model" in [arXiv:1607.01314] is a numerical artifact. This new finding prompts the need for a more careful analysis of the numerical validity of other similar results related to gravitational waves from oscillon dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗

Post-processing interstitialcy diffusion from molecular dynamics simulations

An algorithm to rigorously trace the interstitialcy diffusion trajectory in crystals is developed. The algorithm incorporates unsupervised learning and graph optimization which obviate the need to input extra domain specific information depending on crystal or temperature of the simulation. The algorithm is implemented in a flexible framework as a post-processor to molecular dynamics (MD) simulations. We describe in detail the reduction of interstitialcy diffusion into known computational problems of unsupervised clustering and graph optimization. We also discuss the steps, computational efficiency and key components of the algorithm. Using the algorithm, thermal interstitialcy diffusion from low to near-melting point temperatures is studied. We encapsulate the algorithms in a modular framework with functionality to calculate diffusion coefficients, migration energies and other trajectory properties. The study validates the algorithm by establishing the conformity of output parameters with experimental values and provides detailed insights for the interstitialcy diffusion mechanism. The algorithm along with the help of supporting visualizations and analysis gives convincing details and a new approach to quantifying diffusion jumps, jump-lengths, time between jumps and to identify interstitials from lattice atoms. -- Graphical abstract:.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Laser-induced non-equilibrium plasma kernel dynamics

A non-equilibrium model for laser generated plasmas is validated and used to investigate plasma kernel dynamics. Breakdowns in quiescent air at different ambient conditions (e.g. pressure) for different laser parameters (e.g. power density and beam wavelength) are considered. The hydrodynamics is described based on the chemically reactive Navier–Stokes equations where non-equilibrium effects are accounted for using a three-temperature model. Laser-plasma interactions over nano-second time-scales are modeled with a radiative transfer equation (RTE), including both multiphoton ionization (MPI) and inverse Bremsstrahlung (IB). Validation is conducted by comparing absorbed energy and plasma emission with the corresponding experimental measurements, and by observations of the structure of the breakdown. The laser generated plasma exhibits a two-lobe structure, developing both in the forward and backward directions. An analysis of the dynamics suggests that the kernel evolution results from a radiation-driven wave that is: (i) triggered by MPI, (ii) sustained by energy deposition in IB interactions, which compensates for the energy lost by free-electrons in ionizing collisions, and (iii) guided by both MPI and ionization by electron impact (IE).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation of High Frequency Response in a Dynamically Loaded, Nonlinear Soil Column

Detailed guidance on linear seismic analysis of soil columns is provided in “Seismic Analysis of Safety-Related Nuclear Structures and Commentary (ASCE 4, 1998),” which is currently under revision. A new Appendix in ASCE 4-2014 (draft) is being added to provide guidance for nonlinear time domain analysis which includes evaluation of soil columns. When performing linear analysis, a given soil column is typically evaluated with a linear, viscous damped constitutive model. When submitted to a sine wave motion, this constitutive model produces a smooth hysteresis loop. For nonlinear analysis, the soil column can be modelled with an appropriate nonlinear hysteretic soil model. For the model in this paper, the stiffness and energy absorption result from a defined post yielding shear stress versus shear strain curve. This curve is input with tabular data points. When submitted to a sine wave motion, this constitutive model produces a hysteresis loop that looks similar in shape to the input tabular data points on the sides with discontinuous, pointed ends. This paper compares linear and nonlinear soil column results. The results show that the nonlinear analysis produces additional high frequency response. The paper provides additional study to establish what portion of the high frequency response is due to numerical noise associated with the tabular input curve and what portion is accurately caused by the pointed ends of the hysteresis loop. Finally, the paper shows how the results are changed when a significant structural mass is added to the top of the soil column.

nonlinear↗

Heat load measurements for the PIP-II pHB650 cryomodule

This study presents a brief overview of the 1st and 2nd phases and an in-depth analysis of the 3rd phase heat load testing performed on the pHB650 (prototype High Beta 650 MHz) cryomodule at PIP2IT (PIP-II Injector Test Facility), with a focus on both the results and the methodological advancements that have improved testing efficiency and accuracy. A key challenge identified in the testing campaign is the higher-than-expected heat loads observed in the first PIP-II (Proton Improvement Plan II) prototype cryomodules (pSSR1 and pHB650) tested at PIP2IT. Elevated heat loads are concerning given the fixed capacity of the PIP-II cryoplant that is currently being installed at Fermilab. However, understanding the sources of these elevated heat loads offers a critical opportunity to implement effective heat load mitigations on upcoming PIP-II cryomodules to stay within the available capacity of the PIP-II cryoplant. The study includes a summary of test results, descriptions of measurement procedures, and key observations on parameters directly and indirectly related to heat load measurements. Direct observations include measured heat loads and the effectiveness of JT heat exchanger under varying conditions, while indirect observation analyze factors such as the temperature distribution on the two-phase pipe and relief piping under varying conditions. Thermal acoustic oscillations (TAO) were identified during testing, which was mitigated by replacing the original G10 stem with a stainless steel stem equipped with wipers for the cryomodule cooldown valve. A major innovation during pHB650 Phase 3 testing was the development of an automated Python script to streamline data acquisition, analysis, and reporting of heat load results. This script automatically retrieved data from ACNET (Accelerator Control Network), performed heat load calculations, and generated detailed reports featuring plots and tables. This advancement significantly reduced manual labor and enhanced the thoroughness of data analysis compared to earlier campaigns. The heat load test reports were promptly uploaded to the electronic logbook shortly after each test, enabling rapid feedback and collaboration between the SRF and cryogenic teams. The heat load measurements included various components: HTTS (high-temperature thermal shield), LTTS (low-temperature thermal shield), 2K isothermal and non-isothermal heat loads. Results were recorded both within the cryomodule and between the bayonet can supply and return. Measurements were conducted under different operating conditions such as "standard", "linac", and "simulated dynamic". Additionally, HTTS and LTTS heat loads were calculated in real time, allowing for the tracking of thermal stability and identification of changes during testing, both in steady-state and transient conditions. The results of this testing campaign not only provide valuable insights into the performance of the pHB650 cryomodule but also highlight best practices and lessons learned that will inform future cryomodule testing at PIP2IT. These include adopting automated tools for data analysis, refining real-time measurement capabilities, and emphasizing detailed pre-test planning. The framework established in this campaign aims to set an improved standard for cryomodule testing and heat load reporting in future cryomodule test campaigns.

Porwisiak, D. [Fermilab; Wroclaw Tech. U.]↗

Feedback control design for non-inductively sustained scenarios in NSTX-U using TRANSP

This paper examines a method for real-time control of non-inductively sustained scenarios in NSTX-U by using TRANSP,a time-dependent integrated modeling code for prediction and interpretive analysis of tokamak experimental data, as asimulator. The actuators considered for control in this work are the six neutral beam sources and the plasma boundaryshape. To understand the response of the plasma current, stored energy, and central safety factor to these actuatorsand to enable systematic design of control algorithms, simulations were run in which the actuators were modulated anda linearized dynamic response model was generated. A multi-variable model-based control scheme that accounts for thecoupling and slow dynamics of the system while mitigating the effect of actuator limitations was designed andsimulated. Simulations show that modest changes in the outer gap and heating power can improve the response time ofthe system, reject perturbations, and track target values of the controlled values.

Predictive simulations↗

Dynamics of iodine geminate recombination in supercritical xenon solvent: Caging effect

Understanding the dynamics of chemical reactions in solutions is vital, as their rates and kinetics are significantly affected by the solvent environment. Supercritical solvents offer extensive applications in chemical reactions by enabling the manipulation of the solution environment. Here, in this study, we investigate the geminate recombination of iodine in a supercritical xenon solvent by using ReaxFF-based molecular dynamics simulations. Our findings reveal that the highest iodine recombination rate occurs near supercritical conditions, while lower-pressure conditions lead to reduced collision rates and unstable recombination, and higher-pressure conditions hinder iodine diffusion, resulting in a lower recombination rate. Our analysis shows that the xenon local density at the time of recombination is at least 2.5 times higher than the global density, confirming the presence of xenon clusters surrounding the Iodine atoms. This observation is further supported by coordination number analysis, which confirms an elevated xenon local density during recombination. In addition, the correlation between the total energy of xenon atoms within a cluster and recombined iodine atoms underscores the kinetic energy transfer process, validating the occurrence of geminate recombination. The excess kinetic energy from the recombining iodine atoms is transferred to the surrounding xenon atoms. Our examination of geminate recombination demonstrates that iodine atoms confined within xenon clusters—whether through manual insertion of atoms or the fast dissociation of an iodine molecule within xenon clusters—are more likely to recombine as primary geminate recombination. However, extending the iodine molecule dissociation time allows iodine atoms to diffuse out of the cluster, and the recombination to shift toward secondary geminate recombination.

Cage effect↗