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↗
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)↗
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↗
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↗
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 ↗
DOE OSTI2020

Continuous polymer films deposited on top of porous substrates using plasma-enhanced atomic layer deposition and molecular layer deposition

Molecular layer deposition (MLD) is able to produce ultrathin polymer films with control over thickness, cross-linking, and chemical composition. With these capabilities, MLD should be useful in the fabrication of novel polymer membranes on porous supports. However, confining a continuous MLD film to the surface of porous substrates is difficult because of MLD film growth in the pores. The deposition in the pores lowers the conductance of the porous support. This paper presents a method to deposit continuous polymer films on top of porous substrates. In this method, Al2O3 plasma-enhanced atomic layer deposition (PE-ALD) using trimethylaluminum and oxygen plasma as the reactants was first used to cap the pores of the substrate. Subsequently, a polyamide MLD film was deposited on the Al2O3 PE-ALD capping layer using m-phenylenediamine and trimesoyl chloride as the reactants. The Al2O3 pore caps were then removed from the porous substrate by etching from the backside using a timed exposure to a dilute sodium hydroxide solution. This method was demonstrated using anodic aluminum oxide (AAO) and polyethersulfone (PES) porous substrates. Al2O3 PE-ALD film growth was limited to the top of the porous substrate, resulting in rapid surface recombination or high sticking coefficients for the reactive plasma species within the pores. Gas permeance measurements confirmed the pore capping of the AAO substrates. The reopening of the pores by dissolving the Al2O3 pore caps with a sodium hydroxide solution was monitored using gas permeance versus etch time. The removal of the Al2O3 pore caps from the PES substrates could also dissolve the Al2O3 layer underneath the MLD film. The loss of this Al2O3 layer led to the detachment of the MLD film from the PES substrate. However, the MLD film could be anchored to the PES support at fractures located in the Al2O3 film prior to the MLD. The Al2O3 film fracture allowed the MLD film to anchor firmly to the PES substrate by MLD in the pores of the PES porous substrate. The distance between the anchor points was a function of fracture density. This distance could be controlled by applying a tensile stress to the Al2O3 PE-ALD film to fracture the film through sample bending. This method produced firmly anchored polymer MLD films on top of the PES porous substrates.

Welch, Brian C. (ORCID:000000029747045X)↗
DOE OSTI2023

Colloidal Nafion Particles: Are Cylinders Ubiquitous?

Colloidal Nafion morphology plays a critical role in determining the performance of fuel cells and electrolyzers. While small-angle neutron scattering (SANS) studies previously described Nafion in liquid media as dispersed cylinders, the analysis remains nonunique with multiple possible morphological descriptions of the data. Here, using SANS and all-atomistic molecular dynamics, we confirm that Nafion morphology in liquid media differs substantially depending on dispersing agent and dispersion method. H + Nafion dispersed in N-methyl pyrrolidone forms swollen cluster particles with physically cross-linked ionic groups. Scattering profiles from dispersed Nafion membrane have a large structure factor feature not observed for redispersed Nafion D-521. H + Nafion dispersed in water has a highly elongated cylindrical morphology (radius = 10 ± 1.5 Å, height = 358 ± 4.7 Å) with fully dissociated and solvated sulfonic acid groups on the particle wall. In conclusion, these results highlight an important discrepancy between the methods of preparing Nafion dispersions and the use of simplified analysis techniques to describe Nafion morphology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Microgrids in Emerging Markets - Private Sector Perspectives

This quick read assesses the barriers and opportunities for private sector entry into microgrid development. Information presented in this paper is based on qualitative and quantitative data collection and analysis methods to provide an empirical understanding of barriers to private sector clean energy investment in emerging markets.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Electrically enhanced Haber-Bosch (EEHB) anhydrous ammonia synthesis

The present invention is directed to a method and system for enhancing the production of ammonia from gaseous hydrogen and nitrogen. Advantageously, the method and system does not emit carbon gases during production. The method and system enhances the production of ammonia compared to traditional Haber-Bosch reactions.

Beach, Joseph↗
DOE OSTI2021

Method and apparatus to obtain limited angle tomographic images from stationary gamma cameras

A nuclear imaging system and method for performing three-dimensional imaging of anatomical structures. The system and method includes two or more gamma ray detectors each used in combination with a variable-slant hole collimator. The detectors are positioned in close proximity to, or in contact with, the structure being imaged. The detectors remain in a stationary position during the data collection process. An imaging or reconstruction method is then used to reconstruct a three-dimensional image from the data derived from the detectors.

Kross, Brian J.↗
DOE OSTI2025

Modeled sensitivity of multi-MA accelerator performance to electrode contaminant inventory

Significant particle-in-cell code development has enabled simulations of power flow in multi-MA accelerators to include the desorption of surface contaminants, their ionization into surface plasmas, and the impact of these plasmas on efficiency. The simulations base desorption on an Arrhenius equation, whose most significant unknown is the surface contaminant inventory. The sensitivity of power-flow simulations to this inventory is studied here using Sandia National Laboratories' Z accelerator with a 7-nH MagLIF load [Phys. Plasmas 17, 056303 (2010)]. Simulations are conducted in 3D cylindrical coordinates for the current-adder, or “convolute,” region of Z and in 2D for the final feed only. Simulated contaminant inventories are varied from 1 to 32 monolayers (MLs) in 2D, and 2 to 4 ML in 3D. The results reveal sensitivities to the local ratio of E/B⁠. The high B-field, low E-field region near the short-circuit load is insensitive to the contaminant inventory, where assumed values of 4–32 ML change the load current by ≤ 2%, and agree with experiment to within 2% at peak current. A 1-ML value is the outlier, increasing the load current by 5%, but still within measurement uncertainty. In contrast, the relatively higher E-field, lower B-field convolute region has slower contaminant desorption and higher-magnitude E-field penetration of the surface plasmas. The current loss in the convolute region does increase with contaminant inventory. The loss assuming 4 ML is 12% larger than for 2 ML, with 4 ML being the better match to experiment.

Arrhenius equation↗
DOE OSTI2022

Exploring the Great Pyramid: Detector Technical Design Report with Stand-Alone Monte Carlo Simulations

Cosmic-ray muon imaging has been used to non-destructively examine the Pyramids of Khufu and Khafre on the Giza Plateau; the EGP project will continue this line of research by undertaking a full tomographic scan of the former and by doing so will increase the sensitivity of the technique by upwards of two orders of magnitude. For this purpose, a muon telescope using triangular (vernier) detector cells far outperforms one with rectangular cells, providing an angular and positional RMS improvement by a factor of 4 for the same cost per unit area. A refinement algorithm was developed to handle tracks that produce secondaries. The triangular detector yields a positional resolution for a back-projection to the center of the pyramid of less than 20 cm, ensuring that a feature large enough to be of significance will still be seen by the telescope. This method is shown to be able to accurately and precisely reconstruct tracks of muons that pass through the King s Chamber, Queen s Chamber, and Grand Gallery.

43 PARTICLE ACCELERATORS↗