Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time of arrival estimation”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Sensitivity of Pyrocumulus Convection to Tree Mortality During the 2020 Creek Fire in California

Abstract This study investigated the sensitivity of pyrocumulonimbus (PyroCb) induced by the California Creek fire of 2020 to the amount and type of surface fuels, within the WRF‐SFIRE modeling system. Satellite data were used to derive fire arrival times to constrain fire progression, and to augment the fuel characterization with better estimates of combustible vegetation accounting for tree mortality. Machine learning was employed to classify standing dead vegetation from aerial imagery, which was then added as a custom fuel class along with the standard Anderson fuel categories. Simulations using this new fuel class produced a larger and more vigorous PyroCb than the control run, however, still under‐predicted the cloud top. Additional augmentation of fuel mass to represent the accumulation of dead vegetation on the forest floor further improved the simulations, demonstrating the efficacy of representing both dead standing and fallen vegetation to produce more realistic PyroCb and smoke simulations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

Energetic particles and coronal mass ejections in the high latitude heliosphere: Ulysses-LET observations

The COSPIN Low Energy Telescope (LET) onboard the Ulysses spacecraft measures protons, alphas and heavier ions at energies of approximately 1 to 50 MeV/n. Ulysses measurements offer favorable opportunities to study the effects of solar activity in the out-of-ecliptic regions of the heliosphere. Using LET data, we have investigated the properties of transient energetic ions at high heliographic latitudes when Ulysses was permanently immersed in high speed solar wind and magnetically connected to the Sun on open magnetic field lines. Recurrent increases in the fluxes of energetic ions at high heliographic latitudes at frequencies related to the solar rotation period were found to occur in association with co-rotating interaction regions (CIRs). Here we investigate fluxes of energetic particles that showed no relationship to ClRs. From the investigation of plasma and magnetic held data it is found that all of the transient high latitude particle events were associated with the passage of a coronal mass ejection (CME) over Ulysses. Enhancements in particle fluxes several days prior to the arrival of a CME, but with a significant time delay with respect to the estimated CME-onset at the Sun, were most probably associated with interplanetary shocks driven by fast CMEs. These particle events exhibit unusually high rho/alpha-ratios and are not observed for CMEs not driving a shock. However, not all CMEs that passed Ulysses were associated with a particle event. We find evidence that at high solar latitudes, solar flare particles cannot reach Ulysses on open magnetic field lines, but can reach the spacecraft if particles are injected into magnetic flux-ropes (CMEs) at the Sun. These findings are supported by soft X-ray observations from the Japanese Yohkoh-satellite.

Bothmer, V.↗
NASA NTRS2007

Advanced Systems for Monitoring Underwater Sounds

The term "Passive Acoustic Monitoring System" (PAMS) describes a developmental sensing-and-data-acquisition system for recording underwater sounds. The sounds (more precisely, digitized and preprocessed versions from acoustic transducers) are subsequently analyzed by a combination of data processing and interpretation to identify and/or, in some cases, to locate the sources of those sounds. PAMS was originally designed to locate the sources such as fish of species that one knows or seeks to identify. The PAMS unit could also be used to locate other sources, for example, marine life, human divers, and/or vessels. The underlying principles of passive acoustic sensing and analyzing acoustic-signal data in conjunction with temperature and salinity data are not new and not unique to PAMS. Part of the uniqueness of the PAMS design is that it is the first deep-sea instrumentation design to provide a capability for studying soniferous marine animals (especially fish) over the wide depth range described below. The uniqueness of PAMS also lies partly in a synergistic combination of advanced sensing, packaging, and data-processing design features with features adapted from proven marine instrumentation systems. This combination affords a versatility that enables adaptation to a variety of undersea missions using a variety of sensors. The interpretation of acoustic data can include visual inspection of power-spectrum plots for identification of spectral signatures of known biological species or artificial sources. Alternatively or in addition, data analysis could include determination of relative times of arrival of signals at different acoustic sensors arrayed at known locations. From these times of arrival, locations of acoustic sources (and errors in those locations) can be estimated. Estimates of relative locations of sources and sensors can be refined through analysis of the attenuation of sound in the intervening water in combination with water-temperature and salinity data acquired by instrumentation systems other than PAMS. A PAMS is packaged as a battery-powered unit, mated with external sensors, that can operate in the ocean at any depth from 2 m to 1 km. A PAMS includes a pressure housing, a deep-sea battery, a hydrophone (which is one of the mating external sensors), and an external monitor and keyboard box. In addition to acoustic transducers, external sensors can include temperature probes and, potentially, underwater cameras. The pressure housing contains a computer that includes a hard drive, DC-to- DC power converters, a post-amplifier board, a sound card, and a universal serial bus (USB) 4-port hub.

Lane, Michael↗
NASA NTRS1991

Potential flow applications to complex configurations

Recent advances in CFD methods have enabled the analytic calculation of the carriage loads for stores mounted on complex aircraft. The latest results have demonstrated excellent agreement with test data for the F-15 at M = 0.98. However, in a preliminary design environment, the necessity of generating and validating a Euler grid to fit the aircraft and store arrangement may not be feasible, particularly when effects of configuration changes are considered. For that reason alternative approaches which require less time to arrive at an answer deserve consideration. The paper presents the results of a study to determine if potential flow solutions can give acceptable estimates of store carriage loads at transonic speeds in a timely manner.

Cenko, A.↗
DOE OSTI2022

Earthquake Phase Association Using a Bayesian Gaussian Mixture Model

Earthquake phase association algorithms aggregate picked seismic phases from a network of seismometers into individual seismic events and play an important role in earthquake monitoring and research. Dense seismic networks and improved phase picking methods produce massive seismic phase datasets, particularly for earthquake swarms and aftershocks occurring closely in time and space, making phase association a challenging problem. Here, we present a new association method, the Gaussian Mixture Model Association (GaMMA), that combines the Gaussian mixture model with earthquake location, origin time, and magnitude estimation. We treat earthquake phase association as an unsupervised clustering problem in a probabilistic framework, where each earthquake corresponds to a cluster of P and S phases with a hyperbolic moveout of arrival times and a decay of amplitude with distance. We use the multivariate Gaussian distribution to model the collection of phase picks of an event; and the mean of the multivariate Gaussian distribution is given by the predicted arrival time and amplitude from the causative event. We carry out the pick assignment to each earthquake and determine earthquake source parameters (i.e., earthquake location, origin time, and magnitude) under the maximum likelihood criterion using the Expectation-Maximization algorithm. The GaMMA method does not require typical association steps of other algorithms, such as grid-search or supervised training. The results for both synthetic tests and for the 2019 Ridgecrest earthquake sequence show that GaMMA effectively associates phases from a temporally and spatially dense earthquake sequence while producing useful estimates of earthquake location and magnitude.

58 GEOSCIENCES↗
DOE OSTI2023

Adaptive Time Step Control for Multirate Infinitesimal Methods

Multirate methods have been used for decades to temporally evolve initial-value problems in which different components evolve on distinct time scales, and thus use of different step sizes for these components can result in increased computational efficiency. Generally, such methods select these different step sizes based on experimentation or stability considerations. For problems that evolve on a single time scale, adaptivity approaches that strive to control local temporal error are widely used to achieve numerical results of a desired accuracy with minimal computational effort, while alleviating the need for manual experimentation with different time step sizes. However, there is a notable gap in the publication record on the development of adaptive time step controllers for multirate methods. In this paper, we extend the single-rate controller work of Gustafsson [ACM Trans. Math. Software, 20 (1994), pp. 496-517] to the multirate method setting. Specifically, we develop controllers based on polynomial approximations to the principal error functions for both the "fast" and "slow" time scales within multirate infinitesimal (MRI) methods. We additionally investigate a variety of approaches for estimating the errors arising from each time scale within MRI methods. We then numerically evaluate the proposed multirate controllers and error estimation strategies on a range of multirate test problems, comparing their performance against an estimated optimal performance. Through this work, we combine the most performant of these approaches to arrive at a set of multirate adaptive time step controllers that robustly achieve desired solution accuracy with minimal computational effort.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

An investigation of TNAV equipped aircraft in a simulated en route metering environment

This document presents the results of an effort to estimate how often a TNAV (Time Navigation) equipped aircraft could be given a TNAV clearance in the En Route Metering (ERM) system as a function of the percentage of arriving traffic which is TNAV equipped. A fast-time simulation of Denver Stapleton international arrival traffic in the Denver Air Route Traffic Control Center route structure, including en route metering operations, was used to develop data on estimated conflicts, clearance communications and fuel usage for traffic mixes of 25, 50, 75 and 100% TNAV equipped. This study supports an overall effort by NASA to assess the benefits and required technology for using TNAV-equipped aircraft in the ERM environment.

Groce, J. L.↗
NASA NTRS2019

Adaptive Antenna Pattern Notching of Interference in Synthetic Aperture Radar Data Using Digital Beamforming

Radio Frequency Interference (RFI) is a growing problem in Synthetic Aperture Radar (SAR)systems as scientific motivations push the radars to lower frequencies and as more wireless services share the frequency spectrum. New, advanced SAR instruments, such as NASA's EcoSAR, DBSAR and DLR's Tandem-L mission, employ a multichannel architecture capable of Digital Beamforming (DBF). Radars with DBF are capable of notching the antenna pattern in specific directions, which can be utilized to suppress RFI on board or in post-processing. A well-researched beamformer for this purpose is the Minimum Variance Distortionless Response (MVDR) Beamformer. However, the number of interferers that can be removed through notching is limited by the number of receive channels. It is therefore essential to adaptively change the antenna pattern notching throughout the image in time and frequency for the best results with a given number of receive channels. In this paper, we present several methods to achieve this notching by making use of the spatial SAR signal distribution in range time, range frequency, azimuth time and azimuth Doppler that is inherent to the SAR imaging geometry. Because this distribution is time-variable and yet predictable, it can be used to improve the angle of arrival estimation of the RFI and the adaptive notching. The presented methods can be applied to a Digital Beamforming (DBF) SAR signal with multiple channels in elevation and/or in azimuth. Simulations show that the proposed methods increase the ability to notch out-of-swath interference from multiple directions and lessen the impact on in-swath interference. The improvement of the interferometric coherence for a single-pass interferogram acquired by NASA's EcoSAR system (capable of DBF in elevation) is demonstrated. The removal of periodic RFI artifacts is achieved.

RFI↗
DOE OSTI2022

Heat pulse testing at monitoring wells to estimate subsurface fluid velocities in geological CO 2 storage

Monitoring the injected CO 2 during geological CO 2 storage (GCS) is essential to assure containment and identify CO 2 leakage. Here in this work, a new approach is introduced to estimate the evolution of the downhole fluid velocity at a monitoring well and identify CO 2 arrival time using in-well heat pulse/tracer test. The proposed technique involves using a downhole heater to generate a series of heat pulses and measuring their corresponding temperature response. The surface temperature of the downhole heater is controlled by the supplied electrical power and the heat loss by convection to the surroundings. Convective heat transfer is well described using Newton's law of cooling in which the temperature difference between the heater and the surrounding fluids drives the heat transfer, for which the convection heat transfer coefficient (h) controls the magnitude of heat loss. Among various factors that control h, it depends on the type of the flowing fluid and its velocity. Through analyzing the measured temperature at different heat pulses, the changes in h - due to mobilization of the in-situ brine or CO 2 arrival - can be estimated. Consequently, the velocity of the flowing fluid across the heater can be obtained. Since heat transfer by convection is sensitive to the type of the surrounding fluid, intrusion of CO 2 can be detected from the relatively higher surface temperature obtained at CO 2 arrival. Churchill and Bernstein (1977)'s correlation is adopted to estimate the change of fluid velocity in terms of the change in h. To demonstrate the validity of the proposed technique, the results are applied and validated against those of COMSOL Multiphysics simulation tool for single-phase brine (before CO 2 arrival) and single-phase CO 2 (after CO 2 arrival). The observed temperature heating is sensitive to the flowing fluid velocity and fluid type. The temperature signal observed at CO 2 arrival is large and easily detectable using temperature monitoring tool which provides reliable indication for tracking CO 2 arrival at monitoring wells compared with passive temperature monitoring. The results obtained using the proposed technique agree very well with the numerical results obtained from the simulation tool with a maximum estimation error of 7 percent.

02 PETROLEUM↗
NASA NTRS1988

Estimating Electron Content Of The Ionosphere

Method for estimation total electron content of ionosphere. Based on measurements of signals transmitted from global positioning satellites (GPS's). Ionospheric delays obtained by measuring differential arrival times of signals at two different frequencies. Since GPS observations cover certain regions of sky, location of space probe to be calibrated might not overlap GPS field of view. Certain assumptions made about behavior of ionosphere so total electron content estimated in any direction from receiver.

Lanyi, G. E.↗
NASA NTRS2014

Exploration of the Theoretical Physical Capacity of the John F. Kennedy International Airport Runway System

A design study was completed to explore the theoretical physical capacity (TPC) of the John F. Kennedy International Airport (KJFK) runway system for a northflow configuration assuming impedance-free (to throughput) air traffic control functionality. Individual runways were modeled using an agent-based, airspace simulation tool, the Airspace Concept Evaluation System (ACES), with all runways conducting both departures and arrivals on a first-come first-served (FCFS) scheduling basis. A realistic future flight schedule was expanded to 3.5 times the traffic level of a selected baseline day, September 26, 2006, to provide a steady overdemand state for KJFK runways. Rules constraining departure and arrival operations were defined to reflect physical limits beyond which safe operations could no longer be assumed. Safety buffers to account for all sources of operational variability were not included in the TPC estimate. Visual approaches were assumed for all arrivals to minimize inter-arrival spacing. Parallel runway operations were assumed to be independent based on lateral spacing distances. Resulting time intervals between successive airport operations were primarily constrained by same-runway and then by intersecting-runway spacing requirements. The resulting physical runway capacity approximates a theoretical limit that cannot be exceeded without modifying runway interaction assumptions. Comparison with current KJFK operational limits for a north-flow runway configuration indicates a substantial throughput gap of approximately 48%. This gap may be further analyzed to determine which part may be feasibly bridged through the deployment of advanced systems and procedures, and which part cannot, because it is either impossible or not cost-effective to control. Advanced systems for bridging the throughput gap may be conceptualized and simulated using this same experimental setup to estimate the level of gap closure achieved.

Neitzke, Kurt W.↗
DOE OSTI2021

Data-driven Event Identification in the U.S. Power Systems Based on 2D-OLPP and RUSBoosting Trees

Accurate event identification is an essential part of situation awareness ability for power system operators. Therefore, this work proposes an integrated event identification algorithm for power systems. First, to obtain and filter suitable inputs for event identification, an event detection trigger based on the rate of change of frequency (RoCoF) is presented. Then, the wave arrival time difference-based triangulation method considering the anisotropy of wave propagation speed is utilized to estimate the location of the detected event. Next, the two-dimensional orthogonal locality preserving projection (2D-OLPP)-based method, which is suitable for multiple types of measured data, is employed to achieve higher effectiveness in extracting the event features compared with traditional one-dimensional projection and principle component analysis (PCA). Finally, the random undersampling boosted (RUSBoosted) trees-based classifier, which can mitigate the data sample imbalance issue, is utilized to identify the type of the detected event. Furthermore, the proposed approach is demonstrated using the actual measurement data of U.S. power systems from FNET/GridEye. Comparison results show that the proposed event identification algorithm can achieve better performance than existing approaches.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2014

A Long Distance Laser Altimeter for Terrain Relative Navigation and Spacecraft Landing

A high precision laser altimeter was developed under the Autonomous Landing and Hazard Avoidance (ALHAT) project at NASA Langley Research Center. The laser altimeter provides slant-path range measurements from operational ranges exceeding 30 km that will be used to support surface-relative state estimation and navigation during planetary descent and precision landing. The altimeter uses an advanced time-of-arrival receiver, which produces multiple signal-return range measurements from tens of kilometers with 5 cm precision. The transmitter is eye-safe, simplifying operations and testing on earth. The prototype is fully autonomous, and able to withstand the thermal and mechanical stresses experienced during test flights conducted aboard helicopters, fixed-wing aircraft, and Morpheus, a terrestrial rocket-powered vehicle developed by NASA Johnson Space Center. This paper provides an overview of the sensor and presents results obtained during recent field experiments including a helicopter flight test conducted in December 2012 and Morpheus flight tests conducted during March of 2014.

Pierrottet, Diego F.↗
NASA NTRS1995

Profiles of Gamma-Ray Bursts and Their Component Pulses

One physically informative regularity of their otherwise heterogeneous ensemble, is that many Gamma-Ray Bursts consist of well defined pulses. To objectively quantify the temporal structure of BATSE bursts, we have developed an automatic modeling procedure that separates overlapping pulses and determines the energy-dependence of the pulse-shape parameters. No binning of photon arrival times is needed, so when applied to time-tagged events (TTE) the procedure captures variability information down to the shortest time scales present in the raw data. Maximizing the Bayesian likelihood function Pr(data/model) yields estimates of the model parameters, including the number of pulses present, and allows intercomparison of models of different forms. As with any nonlinear optimization, good initial guesses are crucial to avoid convergence to undesirable local minima. We find excellent initial pulse decompositions by wavelet-denoising a cumulative distribution of the raw photon arrival data; differentiation then gives a time profile mostly free of the systematic effects of degraded resolution (as in ordinary Fourier smoothing) and binning. We present statistical information on pulse rise-time, decay-time, peakedness, and amplitudes, plus their energy dependences - both within a single burst and for a large ensemble of bursts.

Scargle, Jeff D.↗
NASA NTRS2023

Real Options Analysis for Valuation of Climate Adaptation Pathways With Application to Transit Infrastructure

Climate change and sea level rise (SLR) are expected to increase the frequency and intensity of coastal flood events, posing risks to coastal communities and infrastructure. While regional climate adaptation investments can provide substantive flood protection, existing plans often neglect uncertainty in future climate conditions and adaptation performance, consequently neglecting the option value of flexibly implementing proposed projects. Addressing this gap, we develop and employ a generalizable real options analysis (ROA) valuation framework that considers how uncertainty in adaptation project costs, SLR, flood severity, and flood losses inform the full range of adaptation performance outcomes. We further propose and apply a novel, computationally efficient flood loss sampling algorithm to estimate the consequences of randomly arriving coastal flood events. We apply this ROA framework to assess the option value of flexibly timing adaptation investments over time, investigating an adaptation pathway proposed by the City of Boston from the perspective of the regional transit system manager. Our results suggest that flexible implementation can provide significant option value in the near-to mid-term(>30 years), with highest option values under low-probability, high consequence scenarios. Our results also suggest adaptation pathway performance in the latter half of the 21stcenturyis most sensitive to uncertainty in sea level rise, flood loss estimates, and flood frequency, underscoring the importance of uncertainty quantification in the long-term valuation of adaptation investments.

Michael V. Martello↗
DOE OSTI2021

Seismic Imaging of the Mw 7.1 Ridgecrest Earthquake Rupture Zone From Data Recorded by Dense Linear Arrays

Abstract We analyze seismograms recorded by four arrays (B1–B4) with 100 m station spacing and apertures of 4–8 km that cross the surface rupture of the 2019 Mw 7.1 Ridgecrest earthquake. The arrays extend from B1 in the northwest to B4 in the southeast of the surface rupture. Delay times between P wave arrivals associated with ∼1,200 local earthquakes and four teleseismic events are used to estimate local velocity variations beneath the arrays. Both teleseismic and local P waves travel faster on the northeast than the southwest side of the fault beneath arrays B1 and B4, but the velocity contrast is less reliably resolved at arrays B2 and B3. We identify several 1–2 km wide low‐velocity zones with much slower inner cores that amplify S waveforms, inferred as damage zones, beneath each array. The damage zones at arrays B2 and B4 also generate fault‐zone head and trapped waves. An automated detector, based on peak ground velocities and durations of high‐amplitude waves, identifies candidate fault‐zone trapped waves (FZTWs) in a localized zone for ∼600 earthquakes at array B4. Synthetic waveform modeling of averaged FZTWs, generated by ∼30 events with high‐quality signals, indicates that the trapping structure at array B4 has a width of ∼300 m, depth of 3–5 km, S wave velocity reduction of ∼20% with respect to the surrounding rock, Q ‐value of ∼30, and S wave velocity contrast of ∼4% across the fault (faster on the northeast side). The results show complex fault‐zone internal structures (velocity contrasts and low‐velocity zones) that vary along fault strike.

Qiu, Hongrui↗
DOE OSTI2026

Seismic Recovery of Acoustic Source Height During the OSIRIS-REx Sample Return Capsule Re-Entry

Controlled hypersonic reentries with well‐constrained trajectories are exceedingly rare, yet they provide valuable reference cases for validating geophysical source characterization methods applicable to bolides and space debris. The 24 September 2023 re‐entry of the OSIRIS‐REx Sample Return Capsule, the first interplanetary capsule to return over the contiguous United States in nearly two decades, presented such an opportunity. Here, we present seismic observations of this event recorded at 16 stations across three sites at ground distances of 7–20 km from the trajectory. Using eigenray tracing through ground‐to‐space atmospheric specifications, we invert picked acoustic‐to‐seismic travel times for effective shock‐generated acoustic source altitudes along the known trajectory. The method recovers source heights of 57.9–59.5 km with travel‐time residuals <1 s across all stations. The seismic‐derived source heights, when compared at their along‐track positions, fall within the range of previously published infrasound‐derived heights for the corresponding trajectory segment, with intermethod differences comparable with the timing‐based uncertainty. These results suggest that air‐to‐ground coupling does not introduce a resolvable systematic timing bias in the picked near‐field seismic arrivals, supporting seismic observations as a complementary tool for altitude estimation of bolides, re‐entry vehicles, and other hypersonic atmospheric sources.

Silber, Elizabeth Allaryce [Sandia National Labora↗
NASA NTRS1985

Mini and super mini arrays for the study of highest energy cosmic rays

The chief difficulty in studying the highest energy cosmic rays is the extremely low intensity, only approximately 5 particles per sq km sr century above 10 to the 20th power eV. Instead of attacking the problem by assembling all of the available resources in one place, as has been done in the past, a way that the task can be performed at much less cost per unit sensitive area is suggested, by using numerous inexpensive mini arrays operating independently of each other. In addition to the quantities usually observed, each mini array will record shower particle arrival time distributions. At 10 to the 20th power eV the saving in cost per primary particle is estimated to be a factor of 10 or better, compared to methods now in use, even for mini arrays newly built from scratch for just this purpose.

Linsley, J.↗