Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied statistics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Searching for parity violation in SDSS DR16 Lyman-α forest data

The four-point correlation function is the lowest order correlation function for scalar fields that can be used to probe statistical parity invariance in an isotropic universe. There are intriguing claims of detection of parity violation in the 4-point function of BOSS galaxy clustering data. We apply the same estimator to the public SDSS Data Release 16 Lyman-α forest data. Lyman-α forest data probes a different redshift range and is sensitive to a different density regime using a completely different technique. A detection would therefore be a strong indication of new physics. We identify an accurate covariance matrix as a crucial impediment to performing this measurement accurately, consistent with existing literature on galaxy 4-point function. Here, we discuss several approaches to estimating the covariance matrix, several of which produce spurious detection. Using a robust, but very suboptimal, covariance matrix derived from subsample bootstrapping, we find no evidence for parity violation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

G-Mapper: Learning a Cover in the Mapper Construction

The Mapper algorithm is a visualization technique in topological data analysis (TDA) that outputs a graph reflecting the structure of a given dataset. However, the Mapper algorithm requires tuning several parameters in order to generate a “nice” Mapper graph. This paper focuses on selecting the cover parameter. We present an algorithm that optimizes the cover of a Mapper graph by splitting a cover repeatedly according to a statistical test for normality. Our algorithm is based on G-means clustering, which searches for the optimal number of clusters in 𝑘-means by iteratively applying the Anderson–Darling test. Our splitting procedure employs a Gaussian mixture model to carefully choose the cover according to the distribution of the given data. In conclusion, experiments for synthetic and real-world datasets demonstrate that our algorithm generates covers so that the Mapper graphs retain the essence of the datasets, while also running significantly faster than a previous iterative method.

G-means clustering↗
DOE OSTI2022

Differentiable Predictions for Large Scale Structure with SHAMNet

In simulation-based models of the galaxy-halo connection, theoretical predictions for galaxy clustering and lensing are typically made based on Monte Carlo realizations of a mock universe. In this paper, we use Subhalo Abundance Matching (SHAM) as a toy model to introduce an alternative to stochastic predictions based on mock population, demonstrating how to make simulation-based predictions for clustering and lensing that are both exact and differentiable with respect to the parameters of the model. Conventional implementations of SHAM are based on iterative algorithms such as Richardson-Lucy deconvolution; here we use the JAX library for automatic differentiation to train SHAMNet, a neural network that accurately approximates the stellar-to-halo mass relation (SMHM) defined by abundance matching. In our approach to making differentiable predictions of large-scale structure, we map parameterized PDFs onto each simulated halo, and calculate gradients of summary statistics of the galaxy distribution by using autodiff to propagate the gradients of the SMHM through the statistical estimators used to measure one- and two-point functions. Our techniques are quite general, and we conclude with an overview of how they can be applied in tandem with more complex, higher-dimensional models, creating the capability to make differentiable predictions for the multi-wavelength universe of galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Exploring Advanced Computational Tools and Techniques with Artificial Intelligence and Machine Learning in Operating Nuclear Plants

This report presents the project Idaho National Laboratory conducted for Nuclear Regulatory Commission to explore the advanced computational tools and techniques, such as artificial intelligence (AI) and machine learning (ML), for operating nuclear plants. The report reviews the nuclear data sources, with the focus on the operating experience data, that could be applied by advanced computational tools and techniques. Plant-specific and generic (national and international) data from different sources are described. The report describes the relationships between statistics and AI/ML and then introduces the most widely used AI/ML algorithms in both supervised and unsupervised learning. The report reviews the recent applications of advanced computational tools and techniques in various fields of nuclear industry, such as reactor system design and analysis, plant operation and maintenance, and nuclear safety and risk analysis. Finally, the report presents the insights from the project on the potential applicability of AI/ML techniques in improving advanced computational capabilities, how the advanced tools and techniques could contribute to the understanding of safety and risk, and what information would be needed to provide meaningful insights to decision makers. The report also documents an NRC survey on the current state of commercial nuclear power operations relative to the use of AI and ML tools as well as the role of AI/ML tools in nuclear power operations was published by the NRC as in FRN NRC-2021-0048 in April 2021. A summary of the survey including the survey questions, survey participants, survey responses, and the conclusions and insights derived from the survey is provided in the report. Finally, the report investigates potential applications of using AI/ML in operating NPPs and advanced reactors (both advanced LWRs and advanced NLWRs) to improve nuclear plant safety and efficiency. Three main application fields are defined and discussed: (1) plant safety and security assessments; (2) plant degradation modeling, fault and accident diagnosis and prognosis; and (3) plant operation and maintenance efficiency improvement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Assessment of Baseline Monitoring Data for the Remote-Handled Low-Level Waste Disposal Facility at Idaho National Laboratory

The purpose of this document is to summarize environmental monitoring data collected at the Remote-Handled Low-Level Waste (RHLLW) Disposal Facility during the first four years of facility operations (FY 2019 through FY 2022). This summary provides a “baseline” condition for the facility, which can be used to distinguish contaminant releases from the RHLLW Disposal Facility from pre existing contamination as well as potential future releases from other sources (i.e., upgradient aquifer sources), as specified in the facility monitoring plan. Sufficient data has been collected over the first four years of operation of the RHLLW Disposal Facility to establish baseline conditions for future compliance monitoring of aquifer wells and performance monitoring of vadose zone lysimeters. Except for seven elevated tritium measurements from lysimeter HFEF-South believed to have been impacted by waste disposals, all data was deemed appropriate for establishing baseline concentrations. Monitoring of the aquifer detected indicator analytes gross alpha and gross beta and target analytes tritium and C-14. There were a few C-14 detections prior to an increase in the required detection level (RDL) after the first round of sampling. I-129 and Tc-99 were not detected above RDLs in aquifer samples. Lysimeter samples detected indicator analytes gross alpha and gross beta and target analyte tritium. Target analytes C-14, I-129, and Tc-99 were not detected above RDLs in lysimeter samples. Except for tritium, baseline concentrations of indicator and target analytes measured during the first four years of RHLLW facility operations are expected to represent conditions through the expected 20 year operating period of the facility. Elevated tritium in the aquifer, a result of past discharges of tritium at upgradient facilities, showed a decline over the first four years of operations, and levels are expected to continue to decline with time as a result of dilution and decay. Current tritium levels in lysimeter samples may be elevated due to tritium in the water applied during construction and infiltration testing of the facility. As a result, tritium concentrations in lysimeters may decline with time as a result of dilution and decay. Statistical measures (i.e., mean, standard deviation, and 99% upper confidence level [UCL]) were calculated for all detected analytes using the four years of concentration data. These measures will be used to evaluate future monitoring results and to demonstrate the facility is performing, or not performing, as established in the facility performance assessment (PA). Future monitoring results will also be combined with the baseline data in this report to further establish temporal trends and natural variability in the measurements. Based on the data collected during the baseline period, it is recommended the gross alpha action level for performance monitoring increase from 10 to 20 pCi/L due to several baseline measurements exceeding the preliminary action level of 10. An action level of 20 pCi/L is slightly greater than the 99% UCL, is protective of the aquifer, and would reduce unnecessary sampling. It is also recommended that tritium be added to the lysimeter analyte list for post-baseline period monitoring. Tritium, while not a dose concern, is a good tracer that can provide valuable information on water flow in and around the RHLLW Disposal Facility. Recommendations for a tritium action level are provided if considered necessary.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

VitriEdge: Repairable & Durable Vitrimer Coatings for Wind Turbine Blade Leading Edges

The primary goal of this Level 1b incubator project was to explore the use of vitrimer coatings for repair of leading-edge erosion on end-of-life wind turbine blade surfaces, beyond coating strength of adhesion which has previously been demonstrated in the Level 1a project. Uniform vitrimer coatings (thickness: 400 µm) were applied to two end-of-life wind turbine blades for flexural, fatigue, and laminate tensile testing where the addition of the coating did not produce any statistical variation in tensile properties with minor drops in flexural strength for some laminate formulations. However, a <2% variation in storage modulus was measured for laminate structures (i.e., blade samples with vitrimer coatings) across 100,000 flexural cycles and upon laminate tensile failure, the vitrimer coatings displayed no visible signs of delamination. In addition, three methods to heal vitrimer coating damage was displayed: oven heating, addition of hot water, and a forced convection heat gun. All three heating and healing mechanisms demonstrated significant healing with scratch depths decreasing between 79-91% at healing times ranging between 1-min and 10-minutes. Finally, a water jet machine was used to simulate rain erosion for both the blade surfaces and vitrimer-coated blade surfaces where the diameter and depth of the damage was recorded as a function of exposure time, water pressure, height of exposure, and angle of exposure. Of interest, while the vitrimer coating did not significantly lessen the damage experienced during rain erosion, the addition of vitrimer composite coatings(5 wt.% mica addition) did result in a crack-resistant, durable coating capable of self-healing behavior and in all cases the angle of rain exposure was the most critical parameter explored. It is crucial to continue exploring this space where vitrimer coatings are of interest for both their self-healing properties and potential use as reversible adhesives.

17 WIND ENERGY↗
NASA NTRS1973

A search for periodic structure in solar 2 cm microwave radiation

A power spectral analysis of 285 hr of 2 cm microwave intensity data showed no statistically significant ( 96% confidence) periodicities in the frequency range 1 to 15 mHz. No correlation was found between 2 cm periodicities and solar activity in H alpha, X-ray, and several microwave frequencies. A small shift of power toward higher frequencies in the power spectrum of the 2 cm data was found to be correlated with solar H alpha and X-ray activity. Using the statistical properties of power spectra, an expression for the ratio of the minimum detectable peak-to-peak to ambient temperature at chromospheric heights may be derived. Applied to a model for oscillation bursts in quiescent supergranules, and using the most significant results of experiments to detect the microwave periodicities, this expression yields an upper limit of approximately .0015.

Sentman, D. D.↗
NASA NTRS1979

A new method of very long-term prediction of solar activity

A method of solar activity prediction based on the statistical relationship between the minimum value of geomagnetic activity in each 11 year cycle and the height of the next sunspot cycle is proposed. The method is applied to predict sunspot cycle 21.

Ohl, A. I.↗
NASA NTRS1983

An analysis of new techniqes for radiometric correction of LANDSAT-4 Thematic Mapper images

The utility of methods for generating TM RLUTS which can improve the quality of the resultant images was investigated. The TM-CCT-ADDS tape was changed to account for a different collection window for the calibration data. Several scenes of Terrebonne Bay, Louisiana and the Grand Bahamas were analyzed to evaluate the radiometric corrections operationally applied to the image data and to investigate several techniques for reducing striping in the images. Printer plots for the TM shutter data were produced and detector statistics were compiled and plotted. These statistics included various combinations of the average shutter counts for each scan before and after DC restore for forward and reverse scans. Results show that striping is caused by the detectors becoming saturated when they view a bright cloud and depress the DC restore level.

Kogut, J.↗
NASA NTRS1985

Evolution of the luminosity function of extragalactic objects

A nonparametric procedure for determination of the evolution of the luminosity function of extragalactic objects and use of this for prediction of expected redshift and luminosity distribution of objects is described. The relation between this statistical evolution of the population and their physical evolution, such as the variation with cosmological epoch of their luminosity and formation rate is presented. This procedure when applied to a sample of optically selected quasars with redshifts less than two shows that the luminosity function evolves more strongly for higher luminosities, indicating a larger quasar activity at earlier epochs and a more rapid evolution of the objects during their higher luminosity phases. It is also shown that absence of many quasars at redshifts greater than three implies slowing down of this evolution in the conventional cosmological models, perhaps indicating that this is near the epoch of the birth of the quasar (and galaxies).

Petrosian, V.↗
NASA NTRS1992

Bispectral analysis of equatorial spread F density irregularities

Bispectral analysis has been applied to density irregularities at frequencies 5-30 Hz observed with a sounding rocket launched from Peru in March 1983. Unlike the power spectrum, the bispectrum contains statistical information about the phase relations between the Fourier components which make up the waveform. In the case of spread F data from 475 km the 5-30 Hz portion of the spectrum displays overall enhanced bicoherence relative to that of the background instrumental noise and to that expected due to statistical considerations, implying that the observed f exp -2.5 power law spectrum has a significant non-Gaussian component. This is consistent with previous qualitative analyses. The bicoherence has also been calculated for simulated equatorial spread F density irregularities in approximately the same wavelength regime, and the resulting bispectrum has some features in common with that of the rocket data. The implications of this analysis for equatorial spread F are discussed, and some future investigations are suggested.

Labelle, J.↗
NASA NTRS1993

Investigation of a geodesy coexperiment to the Gravity Probe B relativity gyroscope program

Geodesy is the science of measuring the gravitational field of and positions on the Earth. Estimation of the gravitational field via gravitation gradiometry, the measurement of variations in the direction and magnitude of gravitation with respect to position, is this dissertation's focus. Gravity Probe B (GP-B) is a Stanford satellite experiment in gravitational physics. GP-B will measure the precession the rotating Earth causes on the space time around it by observing the precessions of four gyroscopes in a circular, polar, drag-free orbit at 650 km altitude. The gyroscopes are nearly perfect niobium-coated spheres of quartz, operating at 1.8 K to permit observations with extremely low thermal noise. The permissible gyroscope drift rate is miniscule, so the torques on the gyros must be tiny. A drag-free control system, by canceling accelerations caused by nongravitational forces, minimizes the support forces and hence torques. The GP-B system offers two main possibilities for geodesy. One is as a drag-free satellite to be used in trajectory-based estimates of the Earth's gravity field. We described calculations involving that approach in our previous reports, including comparison of laser only, GPS only, and combined tracking and a preliminary estimate of the possibility of estimating relativistic effects on the orbit. The second possibility is gradiometry. This technique has received a more cursory examination in previous reports, so we concentrate on it here. We explore the feasibility of using the residual suspension forces centering the GP-B gyros as gradiometer signals for geodesy. The objective of this work is a statistical prediction of the formal uncertainty in an estimate of the Earth's gravitation field using data from GP-B. We perform an instrument analysis and apply two mathematical techniques to predict uncertainty. One is an analytical approach using a flat-Earth approximation to predict geopotential information quality as a function of spatial wavelength. The second estimates the covariance matrix arising in a least-squares estimate of a spherical harmonic representation of the geopotential using GP-B gradiometer data. The results show that the GP-B data set can be used to create a consistent estimate of the geopotential up to spherical harmonic degree and order 60. The formal uncertainty of all coefficients between degrees 5 and 50 is reduced by factors of up to 30 over current satellite-only estimates and up to 7 over estimates which include surface data. The primary conclusion resulting from this study is that the gravitation gradiometer geodesy coexperiment to GP-B is both feasible and attractive.

Everitt, C. W. F.↗
NASA NTRS1994

Calibration of TOPEX/POSEIDON at Platform Harvest

We present estimates for the mean bias of the TOPEX/POSEIDON NASA altimeter (ALT) and the Centre National d'Etudes Spatiales altimeter (SSALT) using in-situ data gathered at Platform Harvest during the first 36 cycles of the mission. Data for 21 overflights of the ALT and six overflights of the SSALT have been analyzed. The analysis includes an independent assessment of in-situ measurements of sea level, the radial component of the orbit, wet tropospheric path delay, and ionospheric path delay. (The sign convention used is such that, to correct the geophysical data record values for sea level, add the bias algebraically. Unless otherwise stated, the uncertainty in a given parameter is depicted by +/- sigma(sub x), where sigma(sub x) is the sample standard deviation of x about the mean.) Tide gauges at Harvest provide estimates of sea level with an uncertainty of +/- 1.5 cm. The uncertainty in the radial component of the orbit is estimated to be +/- 1.3 cm. In-situ measurements of tropopsheric path delay at Harvest compare to within +/- 1.3 cm of the TOPEX/POSEIDON microwave radiometer, and in-situ measurements of the ionospheric path delay compare to within -0.4 +/- 0.7 cm of the dual-frequency ALT and 1.1 +/- 0.6 cm of Doppler orbitography and radiopositioning integrated by satellite. We obtain mean bias estimates of -14.5 +/- 2.9 cm for the ALT and +0.9 +/- 3.1 cm for the SSALT (where the uncertainties are based on the standard deviation of the estimated mean (sigma(sub bar x/y), which is derived from sample statistics and estimates for errors that cannot be observed). These results are consistent with independent estimates for the relative bias between the two altimeters. A linear regression applied to the complete set of data shows that there is a discernable secular trend in the time series for the ALT bias estimates. A preliminary analysis of data obtained through cycle 48 suggests that the apparent secular drift may be the result of a poorly sampled annual signal.

Christensen, E. J.↗
NASA NTRS1990

Analysis of MMU FDIR expert system

This paper describes the analysis of a rulebase for fault diagnosis, isolation, and recovery for NASA's Manned Maneuvering Unit (MMU). The MMU is used by a human astronaut to move around a spacecraft in space. In order to provide maneuverability, there are several thrusters oriented in various directions, and hand-controlled devices for useful groups of them. The rulebase describes some error detection procedures, and corrective actions that can be applied in a few cases. The approach taken in this paper is to treat rulebases as symbolic objects and compute correctness and 'reasonableness' criteria that use the statistical distribution of various syntactic structures within the rulebase. The criteria should identify awkward situations, and otherwise signal anomalies that may be errors. The rulebase analysis agorithms are derived from mathematical and computational criteria that implement certain principles developed for rulebase evaluation. The principles are Consistency, Completeness, Irredundancy, Connectivity, and finally, Distribution. Several errors were detected in the delivered rulebase. Some of these errors were easily fixed. Some errors could not be fixed with the available information. A geometric model of the thruster arrangement is needed to show how to correct certain other distribution nomalies that are in fact errors. The investigations reported here were partially supported by The Aerospace Corporation's Sponsored Research Program.

Landauer, Christopher↗
NASA NTRS1999

The Behavior of Filters and Smoothers for Strongly Nonlinear Dynamics

The Kalman filter is the optimal filter in the presence of known Gaussian error statistics and linear dynamics. Filter extension to nonlinear dynamics is non trivial in the sense of appropriately representing high order moments of the statistics. Monte Carlo, ensemble-based, methods have been advocated as the methodology for representing high order moments without any questionable closure assumptions (e.g., Miller 1994). Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz (1963) model as well as more realistic models of the oceans (Evensen and van Leeuwen 1996) and atmosphere (Houtekamer and Mitchell 1998). A few relevant issues in this context are related to the necessary number of ensemble members to properly represent the error statistics and, the necessary modifications in the usual filter equations to allow for correct update of the ensemble members (Burgers 1998). The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to quite puzzling in that results of state estimate are worse than for their filter analogue (Evensen 1997). In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use Lorenz (1963) model to test and compare the behavior of a variety implementations of ensemble filters. We also implement ensemble smoothers that are able to perform better than their filter counterparts. A discussion of feasibility of these techniques to large data assimilation problems will be given at the time of the conference.

Zhu, Yanqiu↗
NASA NTRS2003

Statistical Evaluation and Improvement of Methods for Combining Random and Harmonic Loads

Structures in many environments experience both random and harmonic excitation. A variety of closed-form techniques has been used in the aerospace industry to combine the loads resulting from the two sources. The resulting combined loads are then used to design for both yield/ultimate strength and high- cycle fatigue capability. This Technical Publication examines the cumulative distribution percentiles obtained using each method by integrating the joint probability density function of the sine and random components. A new Microsoft Excel spreadsheet macro that links with the software program Mathematica to calculate the combined value corresponding to any desired percentile is then presented along with a curve tit to this value. Another Excel macro that calculates the combination using Monte Carlo simulation is shown. Unlike the traditional techniques. these methods quantify the calculated load value with a consistent percentile. Using either of the presented methods can be extremely valuable in probabilistic design, which requires a statistical characterization of the loading. Additionally, since the CDF at high probability levels is very flat, the design value is extremely sensitive to the predetermined percentile; therefore, applying the new techniques can substantially lower the design loading without losing any of the identified structural reliability.

Brown, A. M.↗
NASA NTRS2003

A Lattice Boltzmann Method for Turbomachinery Simulations

Lattice Boltzmann (LB) Method is a relatively new method for flow simulations. The start point of LB method is statistic mechanics and Boltzmann equation. The LB method tries to set up its model at molecular scale and simulate the flow at macroscopic scale. LBM has been applied to mostly incompressible flows and simple geometry.

Hsu, A. T.↗