Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time”

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,729 records · Page 96

Convergence of the multigrid reduction in time algorithm for the linear elasticity equations

This paper presents some recent advances for parallel-in-time methods applied to linear elasticity. With recent computer architecture changes leading to stagnant clock speeds, but ever increasing numbers of cores, future speedups will be available through increased concurrency. Thus, sequential algorithms, such as time stepping, will suffer a bottleneck. This paper explores multigrid reduction in time (MGRIT) for an important application area, linear elasticity. Previously, efforts at parallel-in-time for elasticity have experienced difficulties, for example, the beating phenomenon. As a result, practical parallel-in-time algorithms for this application area currently do not exist. This paper proposes some solutions made possible by MGRIT (e.g., slow temporal coarsening and FCF-relaxation) and, more importantly, a different formulation of the problem that is more amenable to parallel-in-time methods. Using a recently developed convergence theory for MGRIT and Parareal, we show that the changed formulation of the problem avoids the instability issues and allows the reduction of the error using two temporal grids. We then extend our approach to the multilevel case, where we demonstrate how slow temporal coarsening improves convergence. Finally, the paper ends with supporting numerical results showing a practical algorithm enjoying speedup benefits over the sequential algorithm.

97 MATHEMATICS AND COMPUTING↗

An assessment of the utility of multirate time integration for the modeling of laser powder bed fusion

Finite element simulation of the laser powder bed fusion process is made challenging by the inherently multiscale nature of the process. When using the typical global time stepping techniques, slowly-evolving regions of the domain receive the same numerical treatment as the regions with the highest temperature rates. The current work details the implementation and evaluation of an implicit multirate method which is able to advance different regions of the domain with distinct time steps, depending on their current solution rate. Previous work indicates that this representation of the temporal scales of the problem can lead to significant reductions in wall clock run time, and it is shown herein that multirate time integration, when used with a uniform process-scale mesh, can result in speedups of approximately 19, 42, and 87 times, for domains with edge lengths of approximately 1 mm, 2 mm, and 5 mm, respectively. When used in conjunction with $\ h$-refinement (limited to two levels of refinement), the resulting speedups (taken relative to the uniform mesh) are around 28, 85, and 665 times. The method is demonstrated to converge as indicated by the literature, and its use with an AM-Bench domain is demonstrated.

42 ENGINEERING↗

A quantum relaxation-time approximation for finite fermion systems

We propose a relaxation time approximation for the description of the dynamics of strongly excited fermion systems. Our approach is based on time-dependent density functional theory at the level of the local density approximation. This mean-field picture is augmented by collisional correlations handled in relaxation time approximation which is inspired from the corresponding semi-classical picture. The method involves the estimate of microscopic relaxation rates/times which is presently taken from the well established semi-classical experience. The relaxation time approximation implies evaluation of the instantaneous equilibrium state towards which the dynamical state is progressively driven at the pace of the microscopic relaxation time. As test case, we consider Na clusters of various sizes excited either by a swift ion projectile or by a short and intense laser pulse, driven in various dynamical regimes ranging from linear to strongly non-linear reactions. We observe a strong effect of dissipation on sensitive observables such as net ionization and angular distributions of emitted electrons. The effect is especially large for moderate excitations where typical relaxation/dissipation time scales efficiently compete with ionization for dissipating the available excitation energy. Technical details on the actual procedure to implement a working recipe of such a quantum relaxation approximation are given in appendices for completeness.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Time-interval distributions in a digital gamma multi-channel analyzer at extreme input rates

Here, this study explores the relationship between experimentally observed inter-event time-interval distributions (TIDs) and measured dead times for pulses generated by a HPGe gamma detector and processed by a digital spectrometer. The system utilizes a fast and slow channel, pile-up rejector, trapezoid filtering, and flash analog-to-digital converter. The experimentally derived TIDs were compared with theory for validation. The results demonstrate that the theoretical model reliably describes measured TIDs up to 40% dead time. However, significant distortion effects become increasingly pronounced at higher input rates. It appears that the deviation between the measured and calculated time-interval distributions take the shape of higher-order convolutions of the TID, which represent the time difference between counts for more than two successive events. In this work, theoretical functions for the TID are expanded to reproduce measurements up to 90%. This refinement in the interpretation and treatment of the measured TIDs provides improved accuracy and precision in the prediction of the true event rate and measured dead time in the counter. Even though it reflects occasional count loss that is specific to the set-up, it is expected that a similar adjustment may be applicable to other detection systems as well.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Time-Resolved Electrical Scanning Probe Microscopy of Layered Perovskites Reveals Spatial Variations in Photoinduced Ionic and Electronic Carrier Motion

We study light-induced dynamics in thin films comprising Ruddlesden–Popper phases of the layered 2D perovskite (C 4 H 9 NH 3 ) 2 PbI 4 . We probe ionic and electronic carrier dynamics using two complementary scanning probe methods, time-resolved G-mode Kelvin probe force microscopy and fast free time-resolved electrostatic force microscopy, as a function of position, time, and illumination. Here, we show that the average surface photovoltage sign is dominated by the band bending at the buried perovskite–substrate interface. However, the film exhibits substantial variations in the spatial and temporal response of the photovoltage. Under illumination, the photovoltage equilibrates over hundreds of microseconds, a time scale associated with ionic motion and trapped electronic carriers. Surprisingly, we observe that the surface photovoltage of the 2D grain centers evolves more rapidly in time than at the grain boundaries. We propose that the slower evolution at grain boundaries is due to a combination of ion migration occurring between PbI 4 planes, as well as electronic carriers traversing grain boundary traps, thereby changing the time-dependent band unbending at grain boundaries. These results provide a model for the photoinduced dynamics in 2D perovskites and are a useful basis for interpreting photovoltage dynamics on hybrid 2D/3D structures.

14 SOLAR ENERGY↗

River Dynamics Control Transit Time Distributions and Biogeochemical Reactions in a Dam-Regulated River Corridor

Residence Time Distributions (RTDs) exert an important control on biogeochemical translation in the hydrological systems. Previous tracer studies have revealed that RTDs often followed time-invariant exponential, lognormal, or heavy-tailed RTDs that have power-law behavior for long tails in headwater or low-order streams. However, there is a recognition that RTDs can be more complicated and time-variable in response to dynamic hydrological forcing. Here, we use particle tracking to estimate RTDs along the Hanford Reach of the Columbia River and to quantify the influences of river stage fluctuations. Particle tracking is conducted using the velocity field from high-resolution 3D groundwater flow simulations. The effects of hydrological forcing on the residence time distribution were evaluated by varying river flow boundary conditions and releasing particles in different time windows. Our results revealed that dynamic stage fluctuations created rapidly changing losing-gaining conditions in the river and led to highly transient RTDs, which contributes to multiple modes of RTDs. Dam-induced high-frequency (sub-daily) flow variation contributes to the short-time (sub-daily) component of the RTDs. Deviation of the reactant consumption under the single-mode assumption compared to the multimodal RTDs is relatively small (~5%) and appears when the Damköhler number was close to one.More specifically, high-frequency (daily) dam-induced stage variations have higher impacts on biogeochemical reactions with faster reaction rate, and can potentially increase the reactant consumption by 26.68%. Based on these findings, we suggest that current river basin models could be improved by including bank storage and more complex RTDs influenced by both short and long term river stage fluctuations.

54 ENVIRONMENTAL SCIENCES↗

Multitaper Magnitude‐Squared Coherence for Time Series With Missing Data: Understanding Oscillatory Processes Traced by Multiple Observables

To explore the hypothesis of a common source of variability in two time series, observers may estimate the magnitude-squared coherence (MSC), which is a frequency-domain view of the cross correlation. For time series that do not have uniform observing cadence, MSC can be estimated using Welch's overlapping segment averaging. However, multitaper has superior statistical properties to Welch's method in terms of the tradeoff between bias, variance, and bandwidth. The classical multitaper technique has recently been extended to accommodate time series with underlying uniform observing cadence from which some observations are missing. This situation is common for solar and geomagnetic data sets, which may have gaps due to breaks in satellite coverage, instrument downtime, or poor observing conditions. We demonstrate the scientific use of missing-data multitaper magnitude-squared coherence by detecting known solar mid-term oscillations in simultaneous, missing-data time series of solar Lyman α flux and geomagnetic Disturbance Storm Time index. Due to their superior statistical properties, we recommend that multitaper methods be used for all heliospheric time series with underlying uniform observing cadence.

Astro-statistics techniques (1886)↗

Cross-correlation based time delay estimation for turbulent flow velocity measurements: Statistical considerations

Time delay estimation methods (TDE) are well-known techniques to investigate poloidal flows in hot magnetized plasmas through the propagation properties of turbulent structures in the medium. One of these methods is based on the estimation of the time lag at which the cross-correlation function (CCF) estimation reaches its maximum value. The uncertainty of the peak location refers to the smallest determinable flow velocity modulation, and therefore the standard deviation of the time delay imposes important limitation to the measurements. In this article, the relative standard deviation of the CCF estimation and the standard deviation of its peak location are calculated analytically using a simple model of turbulent signals. This model assumes independent (non interacting) overlapping events (coherent structures) with randomly distributed spatio-temporal origins moving with background flow. The result of our calculations is the derivation of a general formula for the CCF variance, which is valid not exclusively in the high event density limit, but also for arbitrary event densities. Our formula reproduces the well known expression for high event densities previously published in the literature. In this paper we also present a derivation of the variance of time delay estimation that turns out to be inversely proportional to the applied time window. The derived formulas were tested in real plasma measurements. The calculations are an extension of the earlier work of Bencze and Zoletnik [Phys. Plasmas 12, 052323 (2005)] where the autocorrelation-width technique was developed. Additionally, we show that velocities calculated by a TDE method possess a broadband noise which originates from this variance, its power spectral density cannot be decreased by worsening the time resolution and can be coherent with noises of other velocity measurements where the same turbulent structures are used. This noise should not be confused with the impact of zero mean frequency zonal flow modulations and can be the explanation for the TEXTOR velocity spectra measured by beam emission spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data acquisition system with pulse height capability for the TOFED time-of-flight neutron spectrometer

A new time-of-flight neutron spectrometer TOFED has been constructed for installation at Experimental Advanced Superconducting Tokamak. A data acquisition system combining measurements of flight time and energy from the interaction of neutrons with the TOFED scintillators has been developed. The data acquisition system can provide a digitizing resolution better than 1.5% (to be compared with the >10% resolution of the recoil particle energy in the plastic scintillators) and a time resolution <1 ns. At the same time, it is compatible with high count rate event recording, which is an essential feature to investigate phenomena occurring on time scales faster than the slowing down time (≈100 ms) of the beam ions in the plasma. Implications of these results on the TOFED capability to resolve fast ion signatures in the neutron spectrum from EAST plasmas are discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Detecting hidden transient events in noisy nonlinear time-series

The information impulse function (IIF), running Variance, and local Hölder Exponent are three conceptually different time-series evaluation techniques. These techniques examine time-series for local changes in information content, statistical variation, and point-wise smoothness, respectively. Using simulated data emulating a randomly excited nonlinear dynamical system, this study interrogates the utility of each method to correctly differentiate a transient event from the background while simultaneously locating it in time. Computational experiments are designed and conducted to evaluate the efficacy of each technique by varying pulse size, time location, and noise level in time-series. Our findings reveal that, in most cases, the first instance of a transient event is more easily observed with the information-based approach of IIF than with the Variance and local Hölder Exponent methods. While our study highlights the unique strengths of each technique, the results suggest that very robust and reliable event detection for nonlinear systems producing noisy time-series data can be obtained by incorporating the IIF into the analysis.

42 ENGINEERING↗

Time delay between photoemission from the 2p and 2s subshells of neon

The R-matrix incorporating time (RMT) method is a method developed recently for solving the time-dependent Schroedinger equation for multielectron atomic systems exposed to intense short-pulse laser light. We have employed the RMT method to investigate the time delay in the photoemission of an electron liberated from a 2p orbital in a neon atom with respect to one released from a 2s orbital following absorption of an attosecond xuv pulse. Time delays due to xuv pulses in the range 76-105 eV are presented. For an xuv pulse at the experimentally relevant energy of 105.2 eV, we calculate the time delay to be 10.2{+-}1.3 attoseconds (as), somewhat larger than estimated by other theoretical calculations, but still a factor of 2 smaller than experiment. We repeated the calculation for a photon energy of 89.8 eV with a larger basis set capable of modeling correlated-electron dynamics within the neon atom and the residual Ne{sup +} ion. A time delay of 14.5{+-}1.5 as was observed, compared to a 16.7{+-}1.5 as result using a single-configuration representation of the residual Ne{sup +} ion.

74 ATOMIC AND MOLECULAR PHYSICS↗

Optical delay encoding for fast timing and detector signal multiplexing in PET

Purpose: The large number of detector channels in modern positron emission tomography (PET) scanners poses a challenge in terms of readout electronics complexity. Multiplexing schemes are typically implemented to reduce the number of physical readout channels, but often result in performance degradation. Novel methods of multiplexing in PET must be developed to avoid this data degradation. The preservation of fast timing information is especially important for time-of-flight PET. Methods: A new multiplexing scheme based on encoding detector interaction events with a series of extremely fast overlapping optical pulses with precise delays is demonstrated in this work. Encoding events in this way potentially allows many detector channels to be simultaneously encoded onto a single optical fiber that is then read out by a single digitizer. A two channel silicon photomultiplier-based prototype utilizing this optical delay encoding technique along with dual threshold time-over-threshold is demonstrated. Results: The optical encoding and multiplexing prototype achieves a coincidence time resolution of 160 ps full width at half maximum (FWHM) and an energy resolution of 13.1% FWHM at 511 keV with 3 × 3 × 5 mm{sup 3} LYSO crystals. All interaction information for both detectors, including timing, energy, and channel identification, is encoded onto a single optical fiber with little degradation. Conclusions: Optical delay encoding and multiplexing technology could lead to time-of-flight PET scanners with fewer readout channels and simplified data acquisition systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Multigrid reduction in time with Richardson extrapolation

The advent of exascale computing will leave many users with access to more computational resources than they can simultaneously use, e.g., billion-way parallelism. In particular, this is true for time-dependent simulations that limit parallelism to the spatial domain. One method to add parallelism in time to existing simulation codes and thus take advantage of ever larger compute resources is Multigrid Reduction in Time (MGRIT). The goal is to achieve a smaller time-to-solution through parallelism in time. In this paper, MGRIT is enhanced with Richardson extrapolation in a cost-efficient way to produce a parallel-in-time method with improved accuracy. Overall, this leads to a large improvement in the accuracy per computational cost of MGRIT.

97 MATHEMATICS AND COMPUTING↗

Localized Exponential Time Differencing Method for Shallow Water Equations: Algorithms and Numerical Study

Here, we explore the performance of the exponential time differencing (ETD) method applied to the rotating shallow water equations. Comparing with explicit time stepping of the same order accuracy in time, the ETD algorithms could reduce the computational time in many cases by allowing the use of large time step sizes while still maintaining numerical stability. To accelerate the ETD simulations, we propose a localized approach that synthesizes the ETD method and overlapping domain decomposition. By dividing the original problem into many subdomain problems of smaller sizes and solving them locally, the proposed approach could speed up the calculation of matrix exponential vector products. Several standard test cases for shallow water equations of one or multiple layers are considered. The results show great potential of the localized ETD method for high-performance computing because each subdomain problem can be naturally solved in parallel at every time step.

58 GEOSCIENCES↗

Short‐term time step convergence in a climate model

Abstract This paper evaluates the numerical convergence of very short (1 h) simulations carried out with a spectral‐element (SE) configuration of the Community Atmosphere Model version 5 (CAM5). While the horizontal grid spacing is fixed at approximately 110 km, the process‐coupling time step is varied between 1800 and 1 s to reveal the convergence rate with respect to the temporal resolution. Special attention is paid to the behavior of the parameterized subgrid‐scale physics. First, a dynamical core test with reduced dynamics time steps is presented. The results demonstrate that the experimental setup is able to correctly assess the convergence rate of the discrete solutions to the adiabatic equations of atmospheric motion. Second, results from full‐physics CAM5 simulations with reduced physics and dynamics time steps are discussed. It is shown that the convergence rate is 0.4—considerably slower than the expected rate of 1.0. Sensitivity experiments indicate that, among the various subgrid‐scale physical parameterizations, the stratiform cloud schemes are associated with the largest time‐stepping errors, and are the primary cause of slow time step convergence. While the details of our findings are model specific, the general test procedure is applicable to any atmospheric general circulation model. The need for more accurate numerical treatments of physical parameterizations, especially the representation of stratiform clouds, is likely common in many models. The suggested test technique can help quantify the time‐stepping errors and identify the related model sensitivities.

54 ENVIRONMENTAL SCIENCES↗

Real-time estimation of the electron temperature profile in DIII-D by leveraging neural-network surrogate models

Control of both the magnitude and the shape of tokamak profiles will be necessary to achieve stable, high-performance plasmas. In order to reject disturbances in real time, feedback-control algorithms rely on accurate real-time knowledge of the plasma state. When diagnostics alone are insufficient, either because they are limited in number or their measurements are too noisy, observers can be used to combine diagnostic data with a response model to provide a better estimation of different plasma properties. An observer has been developed to estimate the electron temperature profile in real time using both diagnostic data from the Thomson scattering system and a model based on the electron heat transport equation describing the evolution of the electron temperature profile. Here, neural network surrogate models are leveraged to help improve the overall model prediction while staying within computation time constraints for real-time use. The observer algorithm is shown in offline tests to produce smooth profiles that are consistent with both the diagnostic data and the electron heat transport equation. When implemented into the real-time plasma control system, this observer will provide valuable information on the electron temperature profile to many potential feedback-control applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient Multigrid Reduction-in-Time for Method-of-Lines Discretizations of Linear Advection

Parallel-in-time methods for partial differential equations (PDEs) have been the subject of intense development over recent decades, particularly for diffusion-dominated problems. It has been widely reported in the literature, however, that many of these methods perform quite poorly for advection-dominated problems. In this report we analyze the particular iterative parallel-in-time algorithm of multigrid reduction-in-time (MGRIT) for discretizations of constant-wave-speed linear advection problems. We focus on common method-of-lines discretizations that employ upwind finite differences in space and Runge-Kutta methods in time. Using a convergence framework we developed in previous work, we prove for a subclass of these discretizations that, if using the standard approach of rediscretizing the fine-grid problem on the coarse grid, robust MGRIT convergence with respect to CFL number and coarsening factor is not possible. This poor convergence and non-robustness is caused, at least in part, by an inadequate coarse-grid correction for smooth Fourier modes in space-time known as characteristic components. We propose an alternative coarse-grid operator that provides a better correction of these modes. This coarse-grid operator is related to previous work and uses a semi-Lagrangian discretization combined with an implicitly treated truncation error correction. Theory and numerical experiments show the proposed coarse-grid operator yields fast MGRIT convergence for many of the method-of-lines discretizations considered, including for both implicit and explicit discretizations of high order. Parallel results demonstrate speed-up over sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗

Multifractal detrended fluctuation analysis of soil radon ( 222 Rn) and thoron ( 220 Rn) time series

In this paper, we present results of Multifractal Detrended Fluctuation Analysis (MF-DFA) of soil Radon ( 222 Rn) and Thoron ( 220 Rn) time series to examine the scaling and multifractal features. Data consists of 15,692 measurements taken each over forty-minute interval. RTM 1688-2, 222 Rn/ 220 Rn electronic device has been used for measurement purposes. Seasonal periodicities have been removed from the data and MF-DFA was employed on de-seasonalized data. Original 222 Rn and 220 Rn time series have been converted into surrogate and shuffled time series. MFDFA has also been applied to the surrogate and shuffled series to examine the multifractality nature. Results of the study show that $^{222}_{Orig}Rn$, $^{222}_{Surr}Rn$ and $^{220}_{Orig}Rn$, $^{220}_{Surr}Rn$ time series are longrange positively correlated and an increase or decrease in either of radionuclide concentrations will be followed by another increase or decrease of concentrations in future measurements. The $^{220}_{Shuff}Rn$ time series has an independent or short-range dependent structure called white noise and the shuffled time series exhibit a very weak multifractality.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗