Search NASA⌕ Search

SEARCH · Search NASA

Results for “Welch method”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Welch Method and Bootstrapping Applied to Subcritical Gamma Noise

We measured the prompt neutron decay constant 𝛼 of the CROCUS zero-power reactor at the Swiss Federal Institute of Technology Lausanne using cross-power spectral density (CPSD) analysis of gamma-gamma correlations from two trans-stilbene organic scintillators positioned near the reactor core. We measured critical and subcritical states, with water levels ranging from 960 mm (critical) to 800 mm (𝜌=−1.4 $ subcritical). Our analysis used the Welch method, dividing signal segments for fast Fourier transform (FFT) frequency analysis and applying bootstrapping uncertainty quantification that uses Welch-defined segments. Results demonstrated a clear increase in the measured 𝛼 as reactor reactivity decreased, distinguishing critical from subcritical conditions. At the 960-mm critical level, 𝛼 was estimated at 155.9 ± 0.7 s −1 , and for the 800-mm subcritical level, 𝛼 increased significantly to 367.3 ± 6.9 s –1 . A linear regression of subcritical states yielded a critical estimate of 154.0 ± 3.1 s –1 , aligning with the static 𝛼 estimate at critical. The bootstrapping method produced normally distributed 𝛼 estimates, confirming data consistency. The gamma CPSD 𝛼 estimates clearly distinguish reactor states and improve monitoring of zero-power reactors. The future deployment of modular and microreactors as potential candidates for noise analysis is demonstrated in CROCUS, particularly zero-power mock-ups of new designs. The improvement of noise analysis in the subcritical domain from this work will support experimental data for reactor deployment and procedure.

CROCUS↗
NASA NTRS1999

Multiple Point Dynamic Gas Density Measurements Using Molecular Rayleigh Scattering

A nonintrusive technique for measuring dynamic gas density properties is described. Molecular Rayleigh scattering is used to measure the time-history of gas density simultaneously at eight spatial locations at a 50 kHz sampling rate. The data are analyzed using the Welch method of modified periodograms to reduce measurement uncertainty. Cross-correlations, power spectral density functions, cross-spectral density functions, and coherence functions may be obtained from the data. The technique is demonstrated using low speed co-flowing jets with a heated inner jet.

Seasholtz, Richard↗
NASA NTRS1998

Transfer function analysis of dynamic cerebral autoregulation in humans

To test the hypothesis that spontaneous changes in cerebral blood flow are primarily induced by changes in arterial pressure and that cerebral autoregulation is a frequency-dependent phenomenon, we measured mean arterial pressure in the finger and mean blood flow velocity in the middle cerebral artery (VMCA) during supine rest and acute hypotension induced by thigh cuff deflation in 10 healthy subjects. Transfer function gain, phase, and coherence function between changes in arterial pressure and VMCA were estimated using the Welch method. The impulse response function, calculated as the inverse Fourier transform of this transfer function, enabled the calculation of transient changes in VMCA during acute hypotension, which was compared with the directly measured change in VMCA during thigh cuff deflation. Beat-to-beat changes in VMCA occurred simultaneously with changes in arterial pressure, and the autospectrum of VMCA showed characteristics similar to arterial pressure. Transfer gain increased substantially with increasing frequency from 0.07 to 0.20 Hz in association with a gradual decrease in phase. The coherence function was > 0.5 in the frequency range of 0.07-0.30 Hz and < 0.5 at < 0.07 Hz. Furthermore, the predicted change in VMCA was similar to the measured VMCA during thigh cuff deflation. These data suggest that spontaneous changes in VMCA that occur at the frequency range of 0.07-0.30 Hz are related strongly to changes in arterial pressure and, furthermore, that short-term regulation of cerebral blood flow in response to changes in arterial pressure can be modeled by a transfer function with the quality of a high-pass filter in the frequency range of 0.07-0.30 Hz.

NASA Discipline Cardiopulmonary↗
DOE OSTI2022

Magnitude-squared Coherence: A Powerful Tool for Disentangling Doppler Planet Discoveries from Stellar Activity

Abstract If Doppler searches for Earth-mass, habitable planets are to succeed, observers must be able to identify and model out stellar activity signals. Here we demonstrate how to diagnose activity signals by calculating the magnitude-squared coherence C ˆ xy 2 ( f ) between an activity-indicator time series x t and the radial-velocity (RV) time series y t . Since planets only cause modulation in RV, not in activity indicators, a high value of C ˆ xy 2 ( f ) indicates that the signal at frequency f has a stellar origin. We use Welch’s method to measure coherence between activity indicators and RVs in archival observations of GJ 581, α Cen B, and GJ 3998. High RV-H α coherence at the frequency of GJ 3998 b and high RV-S index coherence at the frequency of GJ 3998c, indicate that the planets may actually be stellar signals. We also replicate previous results showing that GJ 581 d and g are rotation harmonics and demonstrate that α Cen B has activity signals that are not associated with rotation. Welch’s power spectrum estimates have cleaner spectral windows than Lomb–Scargle periodograms, improving our ability to estimate rotation periods. We find that the rotation period of GJ 581 is 132 days, with no evidence of differential rotation. Welch’s method may yield unacceptably large bias for data sets with N < 75 observations and works best on data sets with N > 100. Tapering the time-domain data can reduce the bias of the Welch’s power spectrum estimator, but observers should not apply tapers to data sets with extremely uneven observing cadence. A software package for calculating magnitude-squared coherence and Welch’s power spectrum estimates is available on github.

Dodson-Robinson, Sarah E. (ORCID:0000000287964974)↗
DOE OSTI2025

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)↗
NASA NTRS2012

Analysis of Dynamic Data from Supersonic Retropropulsion Experiments in NASA Langley's Unitary Plan Wind Tunnel

Recent experimental supersonic retropropulsion tests were conducted at the NASA Langley Research Center Unitary Plan Wind Tunnel Test Section 2 for a range of Mach numbers from 2.4 to 4.6. A 5-inch 70-degree sphere-cone forebody model with a 10-inch cylindrical aftbody experimental model was used which is capable of multiple retrorocket configurations. These configurations include a single central nozzle on the center point of the forebody, three nozzles at the forebody half-radius, and a combination of the first two configurations with no jets being plugged. A series of measurements were achieved through various instrumentation including forebody and aftbody pressure, internal pressures and temperatures, and high speed Schlieren visualization. Specifically, several high speed pressure transducers on the forebody and in the plenum were implemented to look at unsteady flow effects. The following work focuses on analyzing frequency traits due to the unsteady flow for a range of thrust coefficients for single, tri, and quad-nozzle test cases at freestream Mach 4.6 and angle of attack ranging from -8 degrees to +20 degrees. This analysis uses Matlab s fast Fourier transform, Welch's method (modified average of a periodogram), to create a power spectral density and analyze any high speed pressure transducer frequency traits due to the unsteady flow.

Codoni, Joshua R.↗
NASA NTRS1977

Three-dimensional low Reynolds number flows with a free surface

The two-dimensional leveling problem (Degani, Gutfinger, 1976) is extended to three dimensions in the case where the flow Re number is very low and attention is paid to the free surface boundary condition with surface tension effects included. The no-slip boundary condition on the wall is observed. The numerical solution falls back on the Marker and Cell (MAC) method (Harlow and Welch, 1965) with the computation region divided into a finite number of stationary rectangular cells (or boxes in the 3-D case) and fluid flow traverses the cells (or boxes).

Degani, D.↗
DOE OSTI2021

Analytical solution for the cumulative wake of wind turbines in wind farms

This paper solves an approximate form of conservation of mass and momentum for a turbine in a wind farm array. The solution is a fairly simple explicit relationship that predicts the streamwise velocity distribution within a wind farm with an arbitrary layout. As this model is obtained by solving flow-governing equations directly for a turbine that is subject to upwind turbine wakes, no ad hoc superposition technique is needed to predict wind farm flows. A suite of large-eddy simulations (LES) of wind farm arrays is used to examine self-similarity as well as validity of the so-called conservation of momentum deficit for turbine wakes in wind farms. The simulations are performed with and without the presence of some specific turbines in the wind farm. This allows us to systematically study some of the assumptions made to develop the analytical model. A modified version of the conservation of momentum deficit is also proposed to provide slightly better results at short downwind distances, as well as in the far wake of turbines deep inside a wind farm. Model predictions are validated against the LES data for turbines in both full-wake and partial-wake conditions. While our results highlight the limitation in capturing the flow speed-up between adjacent turbine columns, the model is overall able to acceptably predict flow distributions for a moderately sized wind farm. Finally, the paper employs the new model to provide insights on the accuracy of common wake superposition methods.

17 WIND ENERGY↗
NASA NTRS1976

The fast decoding of Reed-Solomon codes using fermat theoretic transforms and continued fractions

It is shown that Reed-Solomon (RS) codes can be decoded by using a fast Fourier transform algorithm over finite fields GF(F sub n) where F sub n is a Fermat prime, and continued fractions. This new transform decoding method is simpler than the standard method for RS codes. The computing time of this new decoding algorithm in software can be faster than the standard decoding method for RS codes.

Welch, L. R.↗
NASA NTRS1989

Cloud field classification based upon high spatial resolution textural features. II - Simplified vector approaches

This paper compares the results of cloud-field classification derived from two simplified vector approaches, the Sum and Difference Histogram (SADH) and the Gray Level Difference Vector (GLDV), with the results produced by the Gray Level Cooccurrence Matrix (GLCM) approach described by Welch et al. (1988). It is shown that the SADH method produces accuracies equivalent to those obtained using the GLCM method, while the GLDV method fails to resolve error clusters. Compared to the GLCM method, the SADH method leads to a 31 percent saving in run time and a 50 percent saving in storage requirements, while the GLVD approach leads to a 40 percent saving in run time and an 87 percent saving in storage requirements.

Chen, D. W.↗
DOE OSTI2023

Automatisierte Generierung von Mikrokinetiken für heterogen katalysierte Reaktionen unter Berücksichtigung korrelierter Unsicherheiten**

Abstract Die Studie stellt eine ab initio basierte Methode für die automatische Konstruktion von detaillierten Mikrokinetiken vor, welche korrelierte Unsicherheiten in allen energetischen Parametern und Schätzroutinen berücksichtigt. Hierfür wurden innerhalb des Unsicherheitsraums des BEEF‐vdW Funktionals 2000 Mikrokinetiken für die Oxidationsreaktionen repräsentativer Abgasemissionen von stöchiometrischen Verbrennungsmotoren auf Pt(111) generiert und durch Multiskalenmodellierung mit Experimenten verglichen. Das Ensemble der Simulationsergebnisse unterstreicht die Bedeutung der Berücksichtigung von Unsicherheiten. Innerhalb dieses Ensembles ist es möglich, eine Mikrokinetik zu identifizieren, welche mit den experimentellen Daten übereinstimmt. Diese Mikrokinetik lässt sich auf ein einziges Austausch‐Korrelations‐Funktional zurückführen und legt nahe, dass Pt(111) das aktive Zentrum für die Oxidation leichter Kohlenwasserstoffe sein könnte. Die Studie bietet einen universellen Rahmen für die automatisierte Konstruktion von Reaktionsmechanismen unter Berücksichtigung korrelierter Unsicherheiten und ermöglicht eine DFT‐basierte Optimierung der Mikrokinetik auch für andere heterogen katalysierte Systeme.

Kreitz, Bjarne↗
NASA NTRS2023

Dynamically Rendering Rough Terrain with Minimal Memory Overhead

Rendering highly detailed terrain is a process with the potential to consume a great deal of a computer’s random access memory (RAM). In a browser-based application, this resource is limited even further, leading to the necessity to use alternative methods of rendering the large amount of data needed for high detail. This report describes one such method that places the onus of rendering on the speed of the graphics processing unit (GPU) rather than on the computer’s memory. By removing attribute buffers, which contribute greatly to memory costs, from the rendering pipeline and generating the requisite attributes on the fly using a heightmap texture instead, it is estimated that memory usage can be cut down to one-sixth that of the previous method.

Visualization↗
NASA NTRS1987

Hubble Space Telescope servicing - Experience base for a new era

An account is given of the conceptual evolution of the Hubble Space Telescope's (HST) maintenance methodology, and the belated recognition of its long-term support requirements, as well as measures conceived for the extension of this orbiting observatory's service life, and the prospective benefits derivable by comparable programs from HST experience. Attention is given to HST cost and capability tradeoffs, the maintenance concept's definition, the integrated logistic support planning instituted by the program and its failure-prediction methods, as well as the institutional support base employed.

Welch, James C.↗
NASA NTRS1978

The fast decoding of Reed-Solomon codes using Fermat theoretic transforms and continued fractions

It is shown that Reed-Solomon (RS) codes can be decoded by using a fast Fourier transform (FFT) algorithm over finite fields GF(F sub n), where F sub n is a Fermat prime, and continued fractions. This new transform decoding method is simpler than the standard method for RS codes. The computing time of this new decoding algorithm in software can be faster than the standard decoding method for RS codes.

Reed, I. S.↗
NASA NTRS1975

Plan for the uniform mapping of earth resources and environmental complexes from Skylab imagery. Assessment of natural vegetation, environmental, and crop analogs

The author has identified the following significant results. For interpreting a wide range of natural vegetation analogs, S-190A color infrared and the ERTS-1 color composite were consistently more useful than were conventional color or black and white photos. Color infrared was superior for five vegetation analogs while color was superior for only three. The errors in identification appeared to associate more with black and white single band images than with multiband color. For rice crop analogs, spectral and spatial discriminations both contribute to the usefulness of images for data collection. Tests and subjective analyses conducted in this study indicated that the spectral bands exploited in color infrared film were the most useful for agricultural crop analysis. Accuracy of crop identification on any single date of Skylab images was less than that of multidate analysis due to differences in crop calendar, cultural practices used, rice variety, planting date, planting method, water use, fertilization, disease, or mechanical problems.

Poulton, C. E.↗
DOE OSTI2026

Comparative Evaluation of Spectral Methods for Robust Reactor Noise Estimation

Reactor noise analysis provides a noninvasive means to determine neutron kinetic parameters from stochastic fluctuations in detector signals. However, standard cross-power spectral density (CPSD) analyses can be sensitive to numerical processing choices, which may introduce processing-dependent systematic shifts in estimates of the prompt neutron decay constant (α) and limit reproducibility. This study uses a hybrid multitaper–Welch spectral estimator to analyze subcritical noise measurements from a fast-spectrum critical assembly. The decay constant α was extracted using three frequency-domain methods: the CPSD, the magnitude-squared coherence (MSC), and the generalized magnitude-squared coherence (GMSC). These coherence-based estimators normalize detector auto-spectral structure and are expected to reduce the sensitivity of fitted α values to processing parameters. A Sobol global sensitivity analysis identified which numerical inputs most strongly influence the fitted values of α. All estimators produced a linear dependence of α on inverse count rate, with delayed-critical extrapolations near 1.7 × 10 4 s −1 , in agreement within 8% of MCNP6.3 KOPTS benchmark calculations. Sensitivity results show that while the CPSD depends on both time-bin width and taper selection, the MSC and GMSC are dominated by time-bin width alone, indicating reduced parameter coupling and greater robustness to processing variability. These findings demonstrate the feasibility and practical value of coherence-based spectral estimators for extracting α from reactor noise and support their broader application to multi-detector and irregular datasets in subcritical system characterization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗