Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods:statistical”

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Constraining Present‐Day Anthropogenic Total Iron Emissions Using Model and Observations

Abstract Iron emissions from human activities, such as oil combustion and smelting, affect the Earth's climate and marine ecosystems. These emissions are difficult to quantify accurately due to a lack of observations, particularly in remote ocean regions. In this study, we used long‐term, near‐source observations in areas with a dominance of anthropogenic iron emissions in various parts of the world to better estimate the total amount of anthropogenic iron emissions. We also used a statistical source apportionment method to identify the anthropogenic components and their sub‐sources from bulk aerosol observations in the United States. We find that the estimates of anthropogenic iron emissions are within a factor of 3 in most regions compared to previous inventory estimates. Under‐ or overestimation varied by region and depended on the number of sites, interannual variability, and the statistical filter choice. Smelting‐related iron emissions are overestimated by a factor of 1.5 in East Asia compared to previous estimates. More long‐term iron observations and the consideration of the influence of dust and wildfires could help reduce the uncertainty in anthropogenic iron emissions estimates.

Meteorology & Atmospheric Sciences↗
DOE OSTI2018

A Locally Adapting Technique for Edge Detection using Image Segmentation

Rapid growth in the field of quantitative digital image analysis is paving the way for scientific researchers to make precise measurements about objects in an image. To compute quantities from an image such as the density of compressed materials or the velocity of a shockwave, object boundaries must first be determined. Images containing regions that each have a spatial trend in intensity are of particular interest here. For edge detection, we present a supervised, statistical image segmentation method that incorporates spatial information to locate boundaries between regions with overlapping intensity histograms, specifically for images where the regions are known but precise boundary locations are unknown. The segmentation of a pixel is determined by comparing its intensity to distributions from nearby pixel intensities, and a gradient of the segmented image indicates edge locations. Because of the statistical nature of the algorithm, we use maximum likelihood estimation to quantify uncertainty about each boundary. We demonstrate the success of this algorithm at locating boundaries and providing uncertainty bands on a radiograph of a multicomponent cylinder and on an optical image of a laser-induced shockwave.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Entropy-based feature selection for capturing impacts in Earth system models with abrupt forcing

This paper presents the development of a new entropy-based feature selection method for identifying and quantifying impacts. Here, impacts are defined as statistically significant differences in spatio-temporal fields when comparing datasets with and without an external forcing in an Earth system model. Temporal feature selection is performed by first computing the cross-fuzzy entropy to quantify similarity of patterns between two datasets and then applying changepoint detection to identify regions of statistically constant entropy. The method is used to capture temperate north surface cooling from a 9-member simulation ensemble of the Mt. Pinatubo volcanic eruption, which injected 10 Tg of SO 2 into the stratosphere. The results estimate a mean difference decrease in near surface air temperature of -0.560 K with a 99% confidence interval between -0.864 K and -0.257 K between April and November of 1992, one year following the eruption. A sensitivity analysis with decreasing SO 2 injection revealed that the impact is statistically significant at 5 Tg but not at 3 Tg. Using identified features, a dependency graph model based on a 9-day lag had significantly fewer nodes than a graph based on monthly means. Furthermore, this demonstrates our method’s ability to perform dimension reduction while still uncovering source-to-impact pathways.

Changepoint detection↗
DOE OSTI2018

A physical-based gas–surface interaction model for rarefied gas flow simulation

Empirical gas–surface interaction models, such as the Maxwell model and the Cercignani–Lampis model, are widely used as the boundary condition in rarefied gas flow simulations. The accuracy of these models in the prediction of macroscopic behavior of rarefied gas flows is less satisfactory in some cases especially the highly non-equilibrium ones. Molecular dynamics simulation can accurately resolve the gas–surface interaction process at atomic scale, and hence can predict accurate macroscopic behavior. They are however too computationally expensive to be applied in real problems. In this work, a statistical physical-based gas–surface interaction model, which complies with the basic relations of boundary condition, is developed based on the framework of the washboard model. In virtue of its physical basis, this new model is capable of capturing some important relations/trends for which the classic empirical models fail to model correctly. As such, the new model is much more accurate than the classic models, and in the meantime is more efficient than MD simulations. Therefore, it can serve as a more accurate and efficient boundary condition for rarefied gas flow simulations. - Highlights: • The new physical-based gas–surface interaction model satisfies the basic relations of boundary condition so that it can be used in rarefied gas flow simulations. • The new physical model possesses several advantages comparing to the empirical gas–surface interaction models, hence it's a more accurate boundary condition for rarefied gas flow study. • The new physical model is a statistical model which can be efficiently implemented by Monte Carlo method. It's much more efficient comparing to the MD simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Recent advances, and unresolved issues, in the application of computational modelling to the prediction of the biological effects of nanomaterials

Nanomaterials research is one of the fastest growing contemporary research areas. The unprecedented properties of these materials have meant that they are being incorporated into products very quickly. Regulatory agencies are concerned they cannot assess the potential hazards of these materials adequately, as data on the biological properties of nanomaterials are still relatively limited and expensive to acquire. Computational modelling methods have much to offer in helping understand the mechanisms by which toxicity may occur, and in predicting the likelihood of adverse biological impacts of materials not yet tested experimentally. This paper reviews the progress these methods, particularly those QSAR-based, have made in understanding and predicting potentially adverse biological effects of nanomaterials, and also the limitations and pitfalls of these methods. - Highlights: • Nanomaterials regulators need good information to make good decisions. • Nanomaterials and their interactions with biology are very complex. • Computational methods use existing data to predict properties of new nanomaterials. • Statistical, data driven modelling methods have been successfully applied to this task. • Much more must be learnt before robust toolkits will be widely usable by regulators.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

A method of surface marker location optimization for tumor motion estimation in lung stereotactic body radiation therapy

Purpose: Accurately localizing lung tumor localization is essential for high-precision radiation therapy techniques such as stereotactic body radiation therapy (SBRT). Since direct monitoring of tumor motion is not always achievable due to the limitation of imaging modalities for treatment guidance, placement of fiducial markers on the patient’s body surface to act as a surrogate for tumor position prediction is a practical alternative for tracking lung tumor motion during SBRT treatments. In this work, the authors propose an innovative and robust model to solve the multimarker position optimization problem. The model is able to overcome the major drawbacks of the sparse optimization approach (SOA) model. Methods: The principle-component-analysis (PCA) method was employed as the framework to build the authors’ statistical prediction model. The method can be divided into two stages. The first stage is to build the surrogate tumor matrix and calculate its eigenvalues and associated eigenvectors. The second stage is to determine the “best represented” columns of the eigenvector matrix obtained from stage one and subsequently acquire the optimal marker positions as well as numbers. Using 4-dimensional CT (4DCT) and breath hold CT imaging data, the PCA method was compared to the SOA method with respect to calculation time, average prediction accuracy, prediction stability, noise resistance, marker position consistency, and marker distribution. Results: The PCA and SOA methods which were both tested were on all 11 patients for a total of 130 cases including 4DCT and breath-hold CT scenarios. The maximum calculation time for the PCA method was less than 1 s with 64 752 surface points, whereas the average calculation time for the SOA method was over 12 min with 400 surface points. Overall, the tumor center position prediction errors were comparable between the two methods, and all were less than 1.5 mm. However, for the extreme scenarios (breath hold), the prediction errors for the PCA method were not only smaller, but were also more stable than for the SOA method. Results obtained by imposing a series of random noises to the surrogates indicated that the PCA method was much more noise resistant than the SOA method. The marker position consistency tests using various combinations of 4DCT phases to construct the surrogates suggested that the marker position predictions of the PCA method were more consistent than those of the SOA method, in spite of surrogate construction. Marker distribution tests indicated that greater than 80% of the calculated marker positions fell into the high cross correlation and high motion magnitude regions for both of the algorithms. Conclusions: The PCA model is an accurate, efficient, robust, and practical model for solving the multimarker position optimization problem to predict lung tumor motion during SBRT treatments. Due to its generality, PCA model can also be applied to other imaging guidance system whichever using surface motion as the surrogates.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Scaled Vecchia Approximation for Fast Computer-Model Emulation

Many scientific phenomena are studied using computer experiments consisting of multiple runs of a computer model while varying the input settings. Gaussian processes (GPs) are a popular tool for the analysis of computer experiments, enabling interpolation between input settings, but direct GP inference is computationally infeasible for large datasets. We adapt and extend a powerful class of GP methods from spatial statistics to enable the scalable analysis and emulation of large computer experiments. Specifically, we apply Vecchia’s ordered conditional approximation in a transformed input space, with each input scaled according to how strongly it relates to the computer-model response. The scaling is learned from the data by estimating parameters in the GP covariance function using Fisher scoring. Our methods are highly scalable, enabling estimation, joint prediction, and simulation in near-linear time in the number of model runs. In several numerical examples, our approach substantially outperformed existing methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Expanded method for the determination of burnup in nuclear fuels using multiple neodymium isotopes

Established methods for the determination of burnup in nuclear fuels commonly rely on the measurement of 148Nd in the spent fuel in concert with known fission product yields to determine the atom percent of fissions in the fuel. This isotope of neodymium is used for various reasons, including chemical and radioactive stability, ease of measurement, and low rates of formation and destruction due to neutron flux apart from fission. However, careful calculation of effective cumulative fission yields and correction factors for (n,γ) capture reactions allows for additional stable and long-lived isotopes of neodymium to be used to provide additional independent measurements of burnup, reducing statistical uncertainty. This method was developed and successfully applied to measure the burnup of compacts from the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program. The mean burnup measured using 143Nd, 145Nd, 146Nd, 148Nd, and 150Nd was statistically observed to be the same as that measured using 148Nd alone, but the statistical uncertainty in the measurement was reduced by a factor of 2, providing a tighter confidence interval in the final results.

Helmreich, Grant [ORNL] (ORCID:0000000330464394)↗
DOE OSTI2015

Bayesian inference on the sphere beyond statistical isotropy

We present a general method for Bayesian inference of the underlying covariance structure of random fields on a sphere. We employ the Bipolar Spherical Harmonic (BipoSH) representation of general covariance structure on the sphere. We illustrate the efficacy of the method as a principled approach to assess violation of statistical isotropy (SI) in the sky maps of Cosmic Microwave Background (CMB) fluctuations. SI violation in observed CMB maps arise due to known physical effects such as Doppler boost and weak lensing; yet unknown theoretical possibilities like cosmic topology and subtle violations of the cosmological principle, as well as, expected observational artefacts of scanning the sky with a non-circular beam, masking, foreground residuals, anisotropic noise, etc. We explicitly demonstrate the recovery of the input SI violation signals with their full statistics in simulated CMB maps. Our formalism easily adapts to exploring parametric physical models with non-SI covariance, as we illustrate for the inference of the parameters of a Doppler boosted sky map. Our approach promises to provide a robust quantitative evaluation of the evidence for SI violation related anomalies in the CMB sky by estimating the BipoSH spectra along with their complete posterior.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

End-Use Savings Shapes Measure Documentation: Dispatch Schedule Generation for Demand Flexibility Measures

This supplemental document describes the methodology used for determining the dispatch timing of various EUSS demand flexibility measures. Demand flexibility measures are designed to reduce/dispatch electricity demand in buildings during especially beneficial/critical times. The method used in this work utilizes predictions of building loads to generate a schedule that reflects the periods when the building's daily peak load occurs to support decision making in demand flexibility measures. The dispatch schedule generation method described in this document creates an hourly schedule that includes a load dispatch (peak) window for each day for a whole year based on load prediction, with options using different prediction methods: perfect prediction, bin-sampling method, fixed schedule, and outdoor air temperature (OAT)-based prediction method. The perfect prediction method performs a simulation to obtain the annual load profile as predicted load, representing the scenario of perfect load prediction. The bin-sampling method (1) categorizes days into representative bins by temperature characteristics, (2) performs simulations on sample days from each of those bins to create representative (or predicted) load, and (3) assigns representative loads for all days in a year based on the bin categorization. The fixed schedule method defines uniform start and end time of peak window with assumed fixed daily peak time, for all days in a season or a year. The OAT-based prediction method uses the statistics of OAT (minimum and maximum) as the indicators of peak load, with specified delay response time from building loads to temperature. Given the load prediction, daily peak periods are determined as a time window with specified length in each day that include the predicted daily peak load and with a secondary rule such as maximizing energy saving potential. The dispatch schedule generation method is not a standalone measure and is intended to be combined with other demand flexibility measures that could leverage the peak schedule and apply demand controls on specific systems or devices for demand response, such as measures described in "Measure Documentation - Thermostat Control for Load Shedding" and "Measure Documentation - Thermostat Control for Load Shifting".

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRSDate not supplied

Development of Solar Flare and Energetic Particle Prediction Portal (SEP 3 )

Solar activity is a primary factor determining the state of the Earth’s space environment, geomagnetic and ionospheric disturbances, and radiation hazards. In the current state of knowledge, machine learning (ML) methods provide essential tools for processing data, investigating relationships among various physical properties and characteristics, uncovering hidden connections, and predicting hazardous solar events. The primary difficulty in developing and applying modern machine-learning tools in heliophysics is that the essential data are scattered among over a hundred data repositories developed by instrument teams of space missions and ground-based observatories. In addition, statistical and ML methods require long time series of homogeneous measurements. To facilitate ML-ready data preparation and access, we have developed an interactive database of solar flares integrating the most essential datasets (https://solarflare.njit.edu/). The database performs an initial data processing and is automatically updated. In addition, we are developing the Solar Energetic Particle Prediction Portal (SEP3, https://sun.njit.edu/SEP3), which hosts web applications that allow users to retrieve the database records. The Portal has a search page for browsing the events from the most widely used catalogs and a dedicated space to share the most recent achievements of the team. The interactive widget can display soft X-ray and proton flux time series from GOES satellites and the flare records. The data portal has been used to evaluate the forecasts of solar proton events and investigate machine-learning approaches to SEP prediction.

SMD↗
DOE OSTI2015

Structural fluctuation of proteins induced by thermodynamic perturbation

A theory to describe structural fluctuations of protein induced by thermodynamic perturbations, pressure, temperature, and denaturant, is proposed. The theory is formulated based on the three methods in the statistical mechanics: the generalized Langevin theory, the linear response theory, and the three dimensional interaction site model (3D-RISM) theory. The theory clarifies how the change in thermodynamic conditions, or a macroscopic perturbation, induces the conformational fluctuation, which is a microscopic property. The theoretical results are applied, on the conceptual basis, to explain the experimental finding by Akasaka et al., concerning the NMR experiment which states that the conformational change induced by pressure corresponds to structural fluctuations occurring in the ambient condition. A method to evaluate the structural fluctuation induced by pressure is also suggested by means of the 3D-RISM and the site-site Kirkwood-Buff theories.

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

Robust statistical reconstruction for charged particle tomography

Systems and methods for charged particle detection including statistical reconstruction of object volume scattering density profiles from charged particle tomographic data to determine the probability distribution of charged particle scattering using a statistical multiple scattering model and determine a substantially maximum likelihood estimate of object volume scattering density using expectation maximization (ML/EM) algorithm to reconstruct the object volume scattering density. The presence of and/or type of object occupying the volume of interest can be identified from the reconstructed volume scattering density profile. The charged particle tomographic data can be cosmic ray muon tomographic data from a muon tracker for scanning packages, containers, vehicles or cargo. The method can be implemented using a computer program which is executable on a computer.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1978

Dynamic response and input identification of MDOF structures subjected to coupled random vector inputs

A method of random dynamic analysis is presented which is based on the classical approach applied to a discretized structure. The method assumes that the system identification is available in the form of natural modes and frequencies. These modes and frequencies can be found from either available solutions or approximately from finite element programs. A computer program has been developed to perform the computations required in the analysis. All computations are performed with transformed modal variables, which results in significant economy since the number of modal degrees of freedom is almost always less than the number of physical degrees of freedom. The program computes the response to force and base inputs which are statistically coupled. A method is also presented for predicting one of the inputs if the second input and the acceleration response at a point on the structure are known. Finally, results are presented for the random response of a rectangular plate subjected to a random pressure and a random base input. The inputs are considered individually and with various degrees of statistical coupling.

Ocallahan, J. C.↗
NASA NTRS2000

A statistical model of the human core-temperature circadian rhythm

We formulate a statistical model of the human core-temperature circadian rhythm in which the circadian signal is modeled as a van der Pol oscillator, the thermoregulatory response is represented as a first-order autoregressive process, and the evoked effect of activity is modeled with a function specific for each circadian protocol. The new model directly links differential equation-based simulation models and harmonic regression analysis methods and permits statistical analysis of both static and dynamical properties of the circadian pacemaker from experimental data. We estimate the model parameters by using numerically efficient maximum likelihood algorithms and analyze human core-temperature data from forced desynchrony, free-run, and constant-routine protocols. By representing explicitly the dynamical effects of ambient light input to the human circadian pacemaker, the new model can estimate with high precision the correct intrinsic period of this oscillator ( approximately 24 h) from both free-run and forced desynchrony studies. Although the van der Pol model approximates well the dynamical features of the circadian pacemaker, the optimal dynamical model of the human biological clock may have a harmonic structure different from that of the van der Pol oscillator.

NASA Discipline Regulatory Physiology↗
DOE OSTI2024

Design, Detection, and Countermeasure of Frequency Spectrum Attack and Its Impact on Long Short-Term Memory Load Forecasting and Microgrid Energy Management

This paper introduces a frequency-domain false data injection attack called Frequency Spectrum Attack (FSA) and explores its effects on load forecasting and the energy management system (EMS) in a microgrid. The FSA analyzes time-series signals in the frequency domain to identify patterns in their frequency spectrum. It learns the distribution of dominant frequencies in a dataset of healthy signals. Subsequently, it manipulates the amplitudes of dominant frequencies within this healthy distribution, ensuring a stealthy attack against statistical analysis of the signal spectrum. We evaluated the performance of FSA on LSTM, a state-of-the-art network for load forecasting. The results show that FSA can triple the Mean Absolute Error (MAE) of predictions compared to the normal case and increase it by 70% compared to noise injection attacks. Furthermore, FSA indirectly enhances battery utilization in the EMS by 45%. We then proposed a detection method that combines statistical analysis and machine-learning-based classification techniques with features. The model effectively distinguishes FSA from healthy and noisy signals, achieving an accuracy of 98.7% and an F1-score of 98.1% on a load dataset, covering healthy, FSA, and noisy load data. Finally, a countermeasure was introduced based on the statistical analysis of the frequency spectrum of healthy signals to mitigate the impact of FSA. This countermeasure successfully reduces the MAE of the attacked model from 0.135 to 0.053, validating its effectiveness in mitigating FSA.

Nazeri, Amirhossein↗
DOE OSTI2023

Measurement of the muon anomalous precession frequency in runs 2 \& 3 of the Muon $g-2$ Experiment at Fermilab

This dissertation presents a measurement of the muon anomalous precession frequency for the Runs 2 \& 3 data of the E989 Muon $g-2$ Experiment at Fermilab. The muon anomalous precession frequency is one of two key inputs, the other being the magnetic field, used to precisely determine the muon magnetic anomaly, $a_{\mu}$. In April 2021, the E989 collaboration reported its first measurement of $a_{\mu}$ to an unprecedented precision of 460 parts-per-billion (ppb). This result is in agreement with the previous measurement performed by the E821 collaboration at Brookhaven National Laboratory, and the combined experimental value is in tension with the Standard Model prediction at 4.2$\sigma$, a possible hint of new physics. The first result from E989 was based on the Run-1 data, which was collected in 2018 and comprises 6\% of the experiment’s target statistics; the Run-1 result was statistics-limited. The Runs 2 \& 3 data were collected in 2019-2020 and constitute a four-fold increase i n statistics compared to Run-1 and consequently, a factor of two reduction in the measurement’s statistical uncertainty. This reduced statistical uncertainty, as well as continued understanding of systematic effects, will result in an even more precise measurement of $a_{\mu}$ and will help clarify the observed tension between theory and experiment. The precession frequency analysis presented in this dissertation includes a number of improvements compared to Boston University's Run-1 analysis: the implementation of a more robust pileup-subtraction procedure, the implementation of a kernel ratio method, the adoption of the statistically optimal asymmetry-weighted method, and investigations that led to an improved understanding of an outstanding residual slow effect. This dissertation motivates a measurement of the muon magnetic anomaly, describes the experimental principle, gives an overview of the E989 experiment, and presents a precession frequency analysis with full systematic unc ertainty evaluation for the Runs 2 \& 3 data. The combined precession frequency measurement, using the ratio asymmetry-weighted method, has a statistical uncertainty of 201 ppb and a systematic uncertainty of 25 ppb, constituting the most precise determination of the muon anomalous precision frequency to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1989

Statistical fluctuations in Monte Carlo calculations

The time counter and modified Nanbu simulation techniques are analyzed, with emphasis placed on the convergence of the calculations to a steady macroscopic state. Such variables as translational and rotational temperature, and flow velocity, sampled at several points in the flowfield, are considered. Both macroscopic averages and molecular distribution functions are analyzed. The calculation of inelastic collisions, in which transfer of energy between translational and internal energy modes is performed, is achieved through the use of the Larsen-Borgnakke phenomenological model. It is noted that, with reference to translational temperature, the time counter method shows less statistical scatter than that found with the modified Nanbu simulation technique.

Boyd, I. D.↗