Predetermined Time-Step Solver for Rapid Quasi-Static Time Series (QSTS) of Distribution Systems.
Abstract not provided.
SEARCH · Search NASA
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.
Abstract not provided.
An expression was derived for the time transformation t - tau, where t is coordinate time in the solar system barycentric space-time frame of reference and tau is proper time obtained from a fixed atomic clock on earth. This transformation is suitable for use in the computation of high-precision earth-based range and Doppler observables of a spacecraft or celestial body located anywhere in the solar system; it can also be used in obtaining computed values of very long baseline interferometry data types. The formulation for computing range and Doppler observables, which is an explicit function of the transformation t - tau, is described briefly.
The long term performance of the eight LORAN-C chains is presented in terms of the Coordinated Universal Time (UTC) of the U.S. Naval Observatory (USNO); and the use of the LORAN-C navigation system for maintaining the user's clock to a UTC scale is described. The atomic time scale and the UTC of several national laboratories and observatories relative to the international atomic time are reported. Typical performance of several NASA tracking station clocks, relative to the USNO master clock, is also presented.
Estimating time and time-lag in time-of-flight velocimeters is investigated. Statistics of a filtered Poisson point process is given. A Maximum Likelihood Estimator is compared with suboptimum estimators in terms of robustness. For a dominating background combined spatial and temporal processing can improve the robustness compared with purely temporal processing. Schemes for the spatial filters are given.
In this paper, we present a method for estimating complex impedances using reflectometry and a modified steepest descent inversion algorithm. We simulate spread spectrum time domain reflectometry (SSTDR), which can measure complex impedances on energized systems for an experimental setup with resistive and capacitive loads. A parametric function, which includes both a misfit function and stabilizer function, is created. The misfit function is a least squares estimate of how close the model data matches observed data. The stabilizer function prevents the steepest descent algorithm from becoming unstable and diverging. Steepest descent iteratively identifies the model parameters that minimize the parametric function. We validate the algorithm by correctly identifying the model parameters (capacitance and resistance) associated with simulated SSTDR data, with added 3 dB white Gaussian noise. With the stabilizer function, the steepest descent algorithm estimates of the model parameters are bounded within a specified range. Furthermore, the errors for capacitance (220pF to 820pF) and resistance (50 Ω to 270 Ω) are < 10%, corresponding to a complex impedance magnitude |R +1/jωC| of 53 Ω to 510 Ω.
As part of their duties, Human Reliability Analysis must often evaluate if crews in nuclear power plants (NPPs) can complete tasks associated with a human-failure event within time limits. For example, the time required in NPP scenarios is determined by systematic and structured walkthroughs, feasibility studies, recorded times from training exercises, and interviews with experienced operators and experts. Typically, a point estimate is derived for the estimate (mean, maximum, or 95th percentile of time required). Using point-estimate values can mask the risk associated with variability among crews, plant conditions and set-up, environmental conditions, and other impact factors under which these actions are executed. While point estimates for time required and time available have served the industry well, without considering the uncertainty they could lead to biased understanding about the risk. The Integrated Human Event Analysis System - General Methodology (IDHEAS-G) model (developed by the US Nuclear Regulatory Commission, NRC) for human error probability calculates human error probability by summing two probabilities: insufficient time and cognitive error. As such, the model takes a more holistic approach by considering the full distributions for time required and time available to calculate the human error probability because the time available to complete the task is insufficient. In this study, we expand on the work of the NRC and discuss methods for estimating these time considerations. For example, for the time required, the impact of Performance Influencing Factors (PIFs) on the distribution was divided into impacts that are aleatory in nature, such as crew-to-crew variability, and those that are epistemic (i.e., the PIFs). Starting with the factors that introduce aleatory uncertainty, a first-order distribution was developed from a large set of time required (i.e., NPP task completion times) data for the range of operator actions that occur in the NPP control room under simulated accident conditions. The first-order distribution can then be adjusted to account for epistemic uncertainty using research associated with the impact of applicable PIFs on the time required. We also develop guidance for analysts to address the probability distributions for the time available. The guidance we developed on how to estimate time required and time available distributions is based on the identification of pertinent research and data, data analyses, and expert knowledge elicitation.
The baseline stability margins for NASA's Space Launch System (SLS) launch vehicle were generated via the classical approach of linearizing the system equations of motion and determining the gain and phase margins from the resulting frequency domain model. To improve the fidelity of the classical methods, the linear frequency domain approach can be extended by replacing static, memoryless nonlinearities with describing functions. This technique, however, does not address the time varying nature of the dynamics of a launch vehicle in flight. An alternative technique for the evaluation of the stability of the nonlinear launch vehicle dynamics along its trajectory is to incrementally adjust the gain and/or time delay in the time domain simulation until the system exhibits unstable behavior. This technique has the added benefit of providing a direct comparison between the time domain and frequency domain tools in support of simulation validation. This technique was implemented by using the Stability Aerospace Vehicle Analysis Tool (SAVANT) computer simulation to evaluate the stability of the SLS system with the Adaptive Augmenting Control (AAC) active and inactive along its ascent trajectory. The gains for which the vehicle maintains apparent time-domain stability defines the gain margins, and the time delay similarly defines the phase margin. This method of extracting the control stability margins from the time-domain simulation is relatively straightforward and the resultant margins can be compared to the linearized system results. The sections herein describe the techniques employed to extract the time-domain margins, compare the results between these nonlinear and the linear methods, and provide explanations for observed discrepancies. The SLS ascent trajectory was simulated with SAVANT and the classical linear stability margins were evaluated at one second intervals. The linear analysis was performed with the AAC algorithm disabled to attain baseline stability margins. At each time point, the system was linearized about the current operating point using Simulink's built-in solver. Each linearized system in time was evaluated for its rigid-body gain margin (high frequency gain margin), rigid-body phase margin, and aero gain margin (low frequency gain margin) for each control axis. Using the stability margins derived from the baseline linearization approach, the time domain derived stability margins were determined by executing time domain simulations in which axis-specific incremental gain and phase adjustments were made to the nominal system about the expected neutral stability point at specific flight times. The baseline stability margin time histories were used to shift the system gain to various values around the zero margin point such that a precise amount of expected gain margin was maintained throughout flight. When assessing the gain margins, the gain was applied starting at the time point under consideration, thereafter following the variation in the margin found in the linear analysis. When assessing the rigid-body phase margin, a constant time delay was applied to the system starting at the time point under consideration. If the baseline stability margins were correctly determined via the linear analysis, the time domain simulation results should contain unstable behavior at certain gain and phase values. Examples will be shown from repeated simulations with variable added gain and phase lag. Faithfulness of margins calculated from the linear analysis to the nonlinear system will be demonstrated.
Abstract The PICOSEC Micromegas precise timing detector is based on a Cherenkov radiator coupled to a photocathode operating in a semi-transparent mode and a Micromegas amplification structure. The first proof of concept single-channel prototype was able to achieve a time resolution below 25 ps. One of the crucial aspects in the development of precise timing gaseous detectors applicable in high-energy physics experiments is a modular design that enables large area coverage. The first 19-channel multi-pad prototype with an active area of approximately 10 cm 2 suffered from degraded timing resolution due to the non-uniformity of the preamplification gap thickness. A new 100 cm 2 detector module with 100 channels based on a rigid hybrid ceramic/FR4 Micromegas board for improved drift gap uniformity was developed. Initial measurements with 80 GeV/c muons showed improvements in timing response over the measured pads and a time resolution below 25 ps. More recent measurements with a thinner drift gap detector module and newly developed RF pulse amplifiers show that the pad centre resolution can be enhanced to the level of 17 ps. This work will present the development of the detector from structural simulations, design, and beam test commissioning with a focus on the timing performance of a thinner drift gap detector module in combination with new electronics using an automated timing scan method.
Phosphonic acid (PA) self-assembled monolayers (SAMs) are utilized at critical interfaces between transparent conductive oxides (TCO) and organic active layers in organic photovoltaic devices (OPVs). The effects of PA deposition method and time on the formation of close-packed, high-quality monolayers is investigated here for SAMs fabricated by solution deposition, micro-contact printing, and spray coating. The solution deposition isotherm for pentafluorinated benzylphosphonic acid (F5BnPA) on indium-doped zinc oxide (IZO) is studied using polarization modulation-infrared reflection-absorption spectroscopy (PM-IRRAS) at room temperature as a model PA/IZO system. Fast surface adsorption occurs within the first min; however, well-oriented high-quality SAMs are reached only after -48 h, presumably through a continual process of molecular adsorption/desorption and monolayer filling accompanied by molecular reorientation. Two other rapid, soak-free deposition techniques, micro-contact printing and spray coating, are also explored. SAM quality is compared for deposition of phenyl phosphonic acid (PPA), F13-octylphosphonic acid (F13OPA), and pentafluorinated benzyl phosphonic acid (F5BnPA) by solution deposition, micro-contact printing and spray coating using PM-IRRAS. In contrast to micro-contact printing and spray coating techniques, 48-168 h solution deposition at both room temperature and 70 degrees C result in contamination- and surface etch-free close-packed monolayers with good reproducibility. SAMs fabricated by micro-contact printing and spray coating are much less well ordered.
The Global Positioning System is an outstanding tool for the dissemination of time. Using mono-channel C/A-code GPS time receivers, the restitution of GPS time through the satellite constellation presents a peak-to-peak discrepancy of several tens of nanoseconds without SA but may be as high as several hundreds of nanoseconds with SA. As a consequence, civil users are more and more interested in implementing hardware and software methods for efficient restitution of GPS time, especially in the framework of the project of a real-time prediction of UTC (UTCp) which could be available in the form of time differences (UTCp - GPS time). Previous work, for improving the real-time restitution of GPS time with SA, to the level obtained without SA, focused on the implementation of a Kalman filter based on past data and updated at each new observation. An alternative solution relies upon the statistical features of the noise brought about by SA; it has already been shown that the SA noise is efficiently reduced by averaging data from numerous satellites observed simultaneously over a sufficiently long time. This method was successfully applied to data from a GPS time receiver, model AOA TTR-4P, connected to the cesium clock kept at the BIPM. This device, a multi-channel, dual frequency, P-code GPS time receiver, is one of the first TTR-4P units in operation in a civil laboratory. Preliminary comparative studies of this new equipment with conventional GPS time receivers are described in this paper. The results of an experimental restitution of GPS time, obtained in June 1993, are also detailed: 3 to 6 satellites were observed simultaneously with a sample interval of 15 s, an efficient smoothing of SA noise was realized by averaging data on all observed satellites over more than 1 hour. When the GPS system is complete in 1994, 8 satellites will be observable continuously from anywhere in the world and the same level of uncertainty will be obtained using a shorter averaging time.
The PICOSEC Micromegas (MM) detector is a precise timing gaseous detector consisting of a Cherenkov radiator combined with a photocathode and a MM amplifying structure. A 100-channel PICOSEC MM prototype with 10 × 10 cm 2 active area equipped with a Cesium Iodide (CsI) photocathode demonstrated a time resolution below σ = 18 ps. The objective of this work is to improve the PICOSEC MM detector robustness aspects, i.e. integration of resistive MM and carbon-based photocathodes, while maintaining good time resolution. The PICOSEC MM prototypes have been tested in laboratory conditions and successfully characterised with 150 GeV/c muon beams at the CERN SPS H4 beam line. The excellent timing performance below σ = 20 ps for an individual pad obtained with the 10 × 10 cm 2 area resistive PICOSEC MM of 20 MΩ/$\square$ showed no significant time resolution degradation as a result of adding a resistive layer. A single-pad prototype equipped with a 12 nm thick Boron Carbide (B 4 C) photocathode presented a time resolution below σ = 35 ps, opening up new possibilities for detectors with robust photocathodes. The results made the concept more suitable for the experiments in need of robust detectors with good time resolution.
State and parameter estimation of power transmission networks is important for monitoring power grid operating conditions and analyzing transient stability. Wind power generation depends on fluctuating input power levels, which are correlated in time and contribute to uncertainty in turbine dynamical models. The ensemble Kalman filter (EnKF), a standard state estimation technique, uses a deterministic forecast and does not explicitly model time-correlated noise in parameters such as mechanical input power. However, this uncertainty affects the probability of fault-induced transient instability and increased prediction bias. Here a novel approach is to model input power noise with time-correlated stochastic fluctuations, and integrate them with the network dynamics during the forecast. While the EnKF has been used to calibrate constant parameters in turbine dynamical models, the calibration of a statistical model for a time-correlated parameter has not been investigated. In this study, twin experiments on a standard transmission network test case are used to validate our time-correlated noise model framework for state estimation of unsteady operating conditions and transient stability analysis, and a methodology is proposed for the inference of the mechanical input power time-correlation length parameter using time-series data from PMUs monitoring power dynamics at generator buses.
SUMMARY Rapid development of time-lapse seismic monitoring instrumentations has made it possible to collect dense time-lapse data for tomographically retrieving time-lapse (even continuous) images of subsurface changes. While traditional time-lapse full waveform inversion (TLFWI) algorithms are designed for sparse time-lapse surveys, they lack of effective temporal constraint on time-lapse data, and, more importantly, lack of the uncertainty estimation of the TLFWI results that is critical for further interpretation. Here, we propose a new data assimilation TLFWI method, using hierarchical matrix powered extended Kalman filter (HiEKF) to quantify the image uncertainty. Compared to existing Kalman filter algorithms, HiEKF allows to store and update a data-sparse representation of the cross-covariance matrices and propagate model errors without expensive operations involving covariance matrices. Hence, HiEKF is computationally efficient and applicable to 3-D TLFWI problems. Then, we reformulate TLFWI in the framework of HiEKF (termed hereafter as TLFWI-HiEKF) to predict time-lapse images of subsurface spatiotemporal velocity changes and simultaneously quantify the uncertainty of the inverted velocity changes over time. We demonstrate the validity and applicability of TLFWI–HiEKF with two realistic CO2 monitoring models derived from Frio-II and Cranfield CO2 injection sites, respectively. In both 2-D and 3-D examples, the inverted high-resolution time-lapse velocity results clearly reveal a continuous velocity reduction due to the injection of CO2. Moreover, the accuracy of the model is increasing over time by assimilating more time-lapse data while the standard deviation is decreasing over lapsed time. We expect TLFWI-HiEKF to be equipped with real-time seismic monitoring systems for continuously imaging the distribution of subsurface gas and fluids in the future large-scale CO2 sequestration experiments and reservoir management.
Power system simulations with long-term data typically have large time steps varying from one second to a few minutes. However, for PV inverter semiconductors, the minimum thermal stress cycle occurs over the fundamental grid frequency (50 or 60 Hz). This requires the time step of the fatigue simulation to be around 100 µs. This small time step requires long computation times to process yearly power production profiles. This paper proposes a fast fatigue simulation for inverter semiconductors using the quasi-static time series (QSTS) simulation concept. The proposed simulation calculates the steady state of the semiconductor junction temperature by using a Fast Fourier Transform (FFT). The small thermal cycling during a switching period and even over the fundamental waveform is disregarded to further accelerate the simulation speed. The resulting time step of the fatigue simulation is 15 minutes, which is consistent with the solar dataset. The error of the proposed simulation is 0.16% compared to the fatigue simulation results using the complete thermal stress profile. A PV inverter that responds to a Transactive Energy System (TES) is simulated to demonstrate the use of the proposed fatigue simulation. The proposed simulation has the potential to co-simulate with system level simulation tools that also adopt the QSTS concept.
An improved coordinated universal time (UTC) system has been adopted by the International Radio Consultative Committee. It was implemented internationally by the standard-frequency and time-broadcast stations on 1 Jan. 1972. The new UTC system eliminates the frequency offset of 300 parts in 10 to the 10th power between the old UTC and atomic time, thus making the broadcast time interval (the UTC second) constant and defined by the resonant frequency of cesium atoms. The new time scale is kept in synchronism with the rotation of the Earth within plus or minus 0.7 s by step-time adjustments of exactly 1 s, when needed. A time code has been added to the disseminated time signals to permit universal time to be obtained from the broadcasts to the nearest 0.1 s for users requiring such precision. The texts of the International Radio Consultative Committee recommendation and report to implement the new UTC system are given. The coding formats used by various standard time broadcast services to transmit the difference between the universal time (UT1) and the UTC are also given. For users' convenience, worldwide primary VLF and HF transmissions stations, frequencies, and schedules of time emissions are also included. Actual time-step adjustments made by various stations on 1 Jan. 1972, are provided for future reference.
The Vision Zero-Safe Systems Approach prioritizes fast access to Emergency Medical Services (EMS) to improve the survivability of road users in transportation crashes, especially concerning the recent increase in pedestrian-involved crashes. Pedestrian crashes resulting in immediate or early death are considerably more severe than those taking longer. The time gap between injury and fatality is known as survival time, and it heavily relies on EMS response time. The characteristics of the crash location may be associated with EMS response and survival time. A US Department of Transportation initiative identifies communities often facing challenges. Six disadvantaged community (DAC) indicators, including economy, environment, equity, health, resilience, and transportation access, enable an analysis of how survival and EMS response times vary across DACs and non-DACs. To this end, this study created a unique and comprehensive database by linking DACs data with 2017–2021 pedestrian-involved fatal crashes. This study utilizes two-stage residual inclusion models with segmentation for DACs and non-DACs accounting for the endogenous relationship between EMS response and pedestrian survival time. The results indicate that EMS response time is higher and pedestrian survival time is lower in DACs than in non-DACs. A delayed EMS response time is associated with a greater reduction in survival time in DACs compared to non-DACs. Factors, e.g., nighttime and interstate crashes, contribute to higher EMS response time, while pedestrian drugs, driver speeding, and hit-and-run behaviors are associated with a greater reduction in survival time in DACs than non-DACs. Finally, the implications of the findings are discussed in the paper.
We explore a novel way to numerically resolve the scaling behavior of finite-time singularities in solutions of nonlinear parabolic PDEs. The Runge–Kutta–Legendre (RKL) and Runge–Kutta–Gegenbauer (RKG) super-time-stepping methods were originally developed for nonlinear complex physics problems with diffusion. These are multistage single step second-order, forward-in-time methods with no implicit solves. The advantage is that the time-step size for stability scales with stage number 𝑠 as $\mathcal{O}$(𝑠 2 ). Many interesting nonlinear PDEs have finite-time singularities, and the presence of diffusion often limits one to using implicit or semi-implicit time-step methods for stability constraints. Finite-time singularities are particularly challenging due to the large range of scales that one desires to resolve, often with adaptive spatial grids and adaptive time steps. Here, in this study, we show two examples of nonlinear PDEs for which the self-similar singularity structure has time and space scales that are resolvable using the RKL and RKG methods, without forcing even smaller time steps. Compared to commonly used implicit numerical methods, we achieve a significantly smaller run time while maintaining comparable accuracy. We also prove numerical monotonicity for both the RKL and RKG methods under their linear stability conditions for the constant coefficient heat equation, in the case of infinite domain and periodic boundary condition, leading to a theoretical guarantee of the superiority of the RKL and RKG methods over traditional super-time-stepping methods, such as the Runge-Kutta-Chebyshev and the orthogonal Runge-Kutta-Chebyshev methods. Code can be found at https://github.com/ZT220501/SRK-Singularity.
In this study, a time-periodic MGRIT algorithm is proposed as a means to reduce the time-to-solution of numerical algorithms by exploiting the time periodicity inherent to many applications in science and engineering. The time-periodic MGRIT algorithm is applied to a variety of linear and nonlinear single- and multiphysics problems that are periodic-in-time. It is demonstrated that the proposed parallel-in-time algorithm can obtain the same time-periodic steady-state solution as sequential time-stepping. It is shown that the required number of MGRIT iterations can be estimated a priori and that the new MGRIT variant can significantly and consistently reduce the time-to-solution compared to sequential time-stepping, irrespective of the number of dimensions, linear or nonlinear PDE models, single-physics or coupled problems and the employed computing resources. The numerical experiments demonstrate that the time-periodic MGRIT algorithm enables a greater level of parallelism yielding faster turnaround, and thus, facilitating more complex and more realistic problems to be solved.