Search NASA⌕ Search

SEARCH · Search NASA

Results for “Probabilities Mathematics”

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 307 records · Page 17

DNA Repair Domain Modeling Can Predict Cell Death and Mutation Frequency for Wide Range Spectrum of Radiation

Exploration missions to Mars and other destinations raise many questions about the health of astronauts. The continuous exposure of astronauts to galactic cosmic rays is one of the main concerns for long-term missions. Cosmic ionizing radiations are composed of different ions of various charges and energies notably, highly charged energy (HZE) particles. The HZE particles have been shown to be more carcinogenic than low-LET radiation, suggesting the severity of chromosomal aberrations induced by HZE particles is one possible explanation. However, most mathematical models predicting cell death and mutation frequency are based on directly fitting various HZE dose response and are in essence empirical approaches. In this work, we assume a simple biological mechanism to model DNA repair and use it to simultaneously explain the low- and high-LET response using the exact same fitting parameters. Our work shows that the geometrical position of DNA repair along tracks of heavy ions are sufficient to explain why high-LET particles can induce more death and mutations. Our model is based on assuming DNA double strand breaks (DSBs) are repaired within repair domain, and that any DSBs located within the same repair domain cluster into one repair unit, facilitating chromosomal rearrangements and increasing the probability of cell death. We introduced this model in 2014 using simplified microdosimetry profiles to predict cell death. In this work, we collaborated with NASA Johnson Space Center to generate more accurate microdosimetry profiles derived by Monte Carlo techniques, taking into account track structure of HZE particles and simulating DSBs in realistic cell geometry. We simulated 224 data points (D, A, Z, E) with the BDSTRACKS model, leading to a large coverage of LET from ~10 to 2,400 keV/μm. This model was used to generate theoretical RBE for various particles and energies for both cell death and mutation frequencies. The RBE LET dependence is in agreement with experimental data known in human and murine cells. It suggests that cell shape and its orientation with respect to the HZE particle beam can modify the biological response to radiation. Such discovery will be tested experimentally and, if proven accurate, will be another strong supporting evidence for DNA repair domains and their critical role in interpreting cosmic radiation sensitivity.

Viger, Louise↗

Correlation Between Weather Alerts and Grid Component Failures for Grid Alert

Weather events cause most grid failures. Often, we even get notifications on our phones to take cover or be prepared for an imminent event. If electric grid utilities had a similar warning that also included probable scenarios and the equipment involved, they could prepare and minimize the effects. Recent research at Idaho National Laboratory into electric grid risk analysis methods resulted in a tool that allows for the development of the most likely scenarios given failure probabilities of grid components. INL has a project with the U.S. Department of Energy’s Cybersecurity, Energy Security, and Emergency Response (CESER) program to develop a Grid Alert application that receives messages from the existing emergency alert system, filters and determines components possibly affected by the emergency event, calculates probable scenarios uses MASTERRI and then notifies the utility if there is significant risk. Historical failure data of elements that comprise the U.S. electric grid have been compiled by utilities and organizations such as the international regulatory body North American Electric Reliability Corporation (NERC). Nominal failure rates are obtained from this data. To make this tool possible, estimated failure rates are needed for different component types given the alert type, severity, and location. Historic weather-related grid element failures are correlated with historic weather events from Integrated Public Alert & Warning System (IPAWS). These correlated events and failures are used along with Bayesian updates from the historical norms to provide a modified failure rate for grid elements in the alert areas and calculate probable scenarios. This discusses the Grid Alert project plan but focuses on the data gathered and process used in determining failure rates for possible grid failure scenarios.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

Preparing angular momentum eigenstates using engineered quantum walks

Coupled angular-momentum eigenstates are widely used in atomic and nuclear physics calculations and are building blocks for spin networks and the Schur transform. To combine two angular momenta J 1 and J 2 , forming eigenstates of their total angular momentum J=J 1 +J 2 , we develop a quantum-walk scheme that does not require inputting O(j 3 ) nonzero Clebsch–Gordan (CG) coefficients classically. In fact, our scheme may be regarded as a unitary method for computing CG coefficients on quantum computers with a typical complexity of O⁡(j) and a worst-case complexity of O⁡(j 3 ). Equivalently, our scheme provides decompositions of the dense CG unitary into sparser unitary operations. Our scheme prepares angular-momentum eigenstates using a sequence of Hamiltonians to move an initial state deterministically to desired final states, which are usually highly entangled states in the computational basis. In contrast with usual quantum walks, whose Hamiltonians are prescribed, we engineer the Hamiltonians in su⁡(2)×su⁡(2), which are inspired by, but different from, Hamiltonians that govern magnetic resonances and dipole interactions. To achieve a deterministic preparation of both ket and bra states, we use projection and destructive interference to double pinch the quantum walks, such that each step is a unit-probability population transfer within a two-level system. We test our state preparation scheme on classical computers, reproducing tables of CG coefficients. Finally, we also implement small test problems on current quantum hardware.

97 MATHEMATICS AND COMPUTING↗

Decomposition of Electron Ionization Mass Spectra for Space Application Using a Monte-Carlo approach

Rationale: Quadrupole mass spectrometers equipped with an electron ionization (EI) sources have been widely used in space exploration to investigate the composition of planetary surfaces and atmospheres. However, the complexity of the samples and the minimal calibration for the fragmentation of molecules in the ionization chambers have prevented the deconvolution of the majority of the mass spectra obtained at different targets, thus limiting the determination of the exact composition of the samples analyzed. We propose a Monte‐Carlo approach to solve this issue mathematically. Methods: We decomposed simulated mass spectra of mixtures acquired with unit resolving power mass spectrometers and EI sources into the sum of the single components fragmentation patterns weighted by their relative concentration using interior‐point least‐square fitting. To fit compounds with poorly known fragmentation patterns, we used a Monte‐Carlo method to vary the intensity of individual fragment ions. We then decomposed the spectrum thousands of times to obtain a statistical distribution. Results: By performing the deconvolution on a mixture of seven different molecules with interfering fragmentation patterns (H2O, O2, CH4, Ar, N2, C2H4, and C2H6) we show that this approach retrieves the mixing ratio of the individual components more accurately than regular mass spectra decomposition methods that rely on fragmentation patterns from general databases. It also provides the probability density function for each species's mixing ratio. Conclusions: By removing the solution degeneracy in the decomposition of mass spectra, the method described herein could significantly increase the scientific retrieval from archived space flight mass spectrometry data, where calibration of the ionization source is no longer an option.

Thomas Gautier↗

A Markov chain model for reliability growth and decay

A mathematical model is developed to describe a complex system undergoing a sequence of trials in which there is interaction between the internal states of the system and the outcomes of the trials. For example, the model might describe a system undergoing testing that is redesigned after each failure. The basic assumptions for the model are that the state of the system after a trial depends probabilistically only on the state before the trial and on the outcome of the trial and that the outcome of a trial depends probabilistically only on the state of the system before the trial. It is shown that under these basic assumptions, the successive states form a Markov chain and the successive states and outcomes jointly form a Markov chain. General results are obtained for the transition probabilities, steady-state distributions, etc. A special case studied in detail describes a system that has two possible state ('repaired' and 'unrepaired') undergoing trials that have three possible outcomes ('inherent failure', 'assignable-cause' 'failure' and 'success'). For this model, the reliability function is computed explicitly and an optimal repair policy is obtained.

Siegrist, K.↗

Real-Time Feature Tracking Using Homography

This software finds feature point correspondences in sequences of images. It is designed for feature matching in aerial imagery. Feature matching is a fundamental step in a number of important image processing operations: calibrating the cameras in a camera array, stabilizing images in aerial movies, geo-registration of images, and generating high-fidelity surface maps from aerial movies. The method uses a Shi-Tomasi corner detector and normalized cross-correlation. This process is likely to result in the production of some mismatches. The feature set is cleaned up using the assumption that there is a large planar patch visible in both images. At high altitude, this assumption is often reasonable. A mathematical transformation, called an homography, is developed that allows us to predict the position in image 2 of any point on the plane in image 1. Any feature pair that is inconsistent with the homography is thrown out. The output of the process is a set of feature pairs, and the homography. The algorithms in this innovation are well known, but the new implementation improves the process in several ways. It runs in real-time at 2 Hz on 64-megapixel imagery. The new Shi-Tomasi corner detector tries to produce the requested number of features by automatically adjusting the minimum distance between found features. The homography-finding code now uses an implementation of the RANSAC algorithm that adjusts the number of iterations automatically to achieve a pre-set probability of missing a set of inliers. The new interface allows the caller to pass in a set of predetermined points in one of the images. This allows the ability to track the same set of points through multiple frames.

Clouse, Daniel S.↗

Sally Ride EarthKAM - Automated Image Geo-Referencing Using Google Earth Web Plug-In

Sally Ride EarthKAM is an educational program funded by NASA that aims to provide the public the ability to picture Earth from the perspective of the International Space Station (ISS). A computer-controlled camera is mounted on the ISS in a nadir-pointing window; however, timing limitations in the system cause inaccurate positional metadata. Manually correcting images within an orbit allows the positional metadata to be improved using mathematical regressions. The manual correction process is time-consuming and thus, unfeasible for a large number of images. The standard Google Earth program allows for the importing of KML (keyhole markup language) files that previously were created. These KML file-based overlays could then be manually manipulated as image overlays, saved, and then uploaded to the project server where they are parsed and the metadata in the database is updated. The new interface eliminates the need to save, download, open, re-save, and upload the KML files. Everything is processed on the Web, and all manipulations go directly into the database. Administrators also have the control to discard any single correction that was made and validate a correction. This program streamlines a process that previously required several critical steps and was probably too complex for the average user to complete successfully. The new process is theoretically simple enough for members of the public to make use of and contribute to the success of the Sally Ride EarthKAM project. Using the Google Earth Web plug-in, EarthKAM images, and associated metadata, this software allows users to interactively manipulate an EarthKAM image overlay, and update and improve the associated metadata. The Web interface uses the Google Earth JavaScript API along with PHP-PostgreSQL to present the user the same interface capabilities without leaving the Web. The simpler graphical user interface will allow the public to participate directly and meaningfully with EarthKAM. The use of similar techniques is being investigated to place ground-based observations in a Google Mars environment, allowing the MSL (Mars Science Laboratory) Science Team a means to visualize the rover and its environment.

Andres, Paul M.↗

Novel Results Visualization for Dynamic PSA and New Modeling Features in EMRALD

The Event Modeling Risk Assessment Linked Diagram (EMRALD) tool, developed at the Idaho National Laboratory (INL), was designed to simplify the creation of dynamic models and support various research projects. One of the primary goals of EMRALD was to provide visual methods for modeling. EMRALD consists of two main components: a web-based user interface for model development and a solve engine for running model simulations. Over time, it has evolved to meet the diverse needs of its users. Initially, EMRALD's results were simple text outputs with final key state percentages and uncertainty bounds. However, because EMRALD utilizes a three-phase discrete event simulation and tracks the paths of each simulation run leading to a key state, there is significant potential to analyze large sets of path results data, including state paths, events, and timing. Visualizing this data meaningfully posed a challenge. To address this, a novel time-based Sankey diagram was developed. EMRALD exports results data in a format that can be opened in this Sankey viewer, allowing users to visualize paths, occurrences, events, and probability data for the entire simulation run in a single diagram. Moreover, when EMRALD was first created, there were limited tools capable of meeting its graphical requirements, many of which are no longer supported. In 2024, a new web-based interface was developed using modern graphing tools, enabling additional modeling features. This paper discusses the new dynamic PSA results visualization capability and the enhanced modeling tools available in EMRALD.

97 - MATHEMATICS AND COMPUTING↗

Interpretable Machine Learning Models for Autonomous Characterization of Analogue Ocean World Seawater Chemistry and Biosignature Potential Using Isotope Ratio Data

Background: Future missions to ocean worlds, such as Enceladus and Europa, will attempt to characterize the subsurface seawater chemistry and assess the potential for life. Such missions will be equipped with capabilities to precisely measure volatile isotopes in plumes, atmospheres, and exospheres. Motivation: While large isotopic fractionations can indicate a biological source, there are signatures resulting from abiotic geochemical processes that mimic isotopic biosignatures. While machine learning (ML) has the potential to disentangle competing effects and biotic mimicry, high-dimensional isotope ratio mass spectrometry (IRMS) data is likely to contain noise/irrelevant features and involve complex statistical interactions that make human inference and interpretation difficult. Further, ML predictions with as far-reaching implications as an extraterrestrial biosignature on an ocean world requires the use of interpretable models (i.e., not “black box” models) with physically and mathematically meaningful feature spaces along with false positive diagnostics. Methods: We use volatile CO2 IRMS data of analogue ocean world seawaters to validate an ML approach to provide biogeochemical context for biosignature detection. We employ a feature selection method called nearest-neighbor projected distance regression (NPDR) that detects statistical interactions and helps elucidate the mechanisms of the Random Forest classification models. Results: We train and validate predictive ML models on volatile CO2 IRMS data of analogue ocean world seawaters to predict major salt components (e.g., MgSO4, NaHCO3), pH, ionic strength, and the presence of biosignatures. Features derived from IRMS measurements are augmented with extracted time-series features. Our results show high test accuracy and interpretability, which is increased by interaction network visualization, sample-wise variable importance scores, and single-sample class probability estimates. We demonstrate an ML mission software solution that triggers autonomous data transmission and biogeochemical sample prediction.

geochemistry↗

Improved multifidelity Monte Carlo estimators based on normalizing flows and dimensionality reduction techniques

Here, we study the problem of multifidelity uncertainty propagation for computationally expensive models. In particular, we consider the general setting where the high-fidelity and low-fidelity models have a dissimilar parameterization both in terms of number of random inputs and their probability distributions, which can be either known in closed form or provided through samples. We derive novel multifidelity Monte Carlo estimators which rely on a shared subspace between the high-fidelity and low-fidelity models where the parameters follow the same probability distribution, i.e., a standard Gaussian. We build the shared space employing normalizing flows to map different probability distributions into a common one, together with linear and nonlinear dimensionality reduction techniques, active subspaces and autoencoders, respectively, which capture the subspaces where the models vary the most. We then compose the existing low-fidelity model with these transformations and construct modified models with an increased correlation with the high-fidelity model, which therefore yield multifidelity estimators with reduced variance. A series of numerical experiments illustrate the properties and advantages of our approaches.

97 MATHEMATICS AND COMPUTING↗

HELPR Version 1.1.0 User Guide

Hydrogen Extremely Low Probability of Rupture (HELPR) is a modular probabilistic fracture mechanics modeling platform developed to assess structural integrity of pipelines for transmission and distribution of hydrogen. HELPR couples fatigue and fracture engineering models with probabilistic methods to generate fast predictions and enables quantification of prediction uncertainty and sensitivity. This user manual serves as a guide through the various analysis features HELPR contains.

08 HYDROGEN↗

Integration of Condition-Based, Diagnostic, Prognostic, And Anomaly Detection Data into Reliability Models to Support a Predictive Maintenance Context

Reliability data employed in plant reliability models are an approximated integral representation of the past industrywide operational experience, and they neglect the present asset health status (available, for example, from online monitoring data and diagnostic assessments) and forecasted health projection (when available from prognostic models). Ideally, in a predictive maintenance context, system reliability models should support decision making by propagating actual health information from the asset to the system level in order to provide a quantitative snapshot of system health and identify the most critical assets. Asset health should be informed solely by that specific asset’s current and historical performance data and should not be an approximated integral representation of the past industrywide operational experience (as currently performed by system reliability models through Bayesian updating processes). This paper proposes a reliability modeling approach that relies on asset diagnostic and prognostic assessments, along with monitoring data to measure asset health. We show how state-of-the art condition-based, diagnostic, prognostic, and anomaly detection models can be linked to system reliability models not in probability terms, but in terms of margin where margin is defined as the “distance” between the present status and an undesired event (e.g., failure or unacceptable performance). Then, we show how the propagation of margin data from the asset to the system level is performed through classical reliability models such as fault trees or reliability block diagrams. The described method is in fact able to propagate heterogenous health data from the asset to the system level in order to analytically assess system health.

97 MATHEMATICS AND COMPUTING↗

A Eulerian three-fluid flow framework for the study of fuel dispersal behavior under loss-of-coolant accident conditions

As the nuclear industry explores new fuel designs to accommodate increased burnup, studying fuel behavior during loss-of-coolant accidents is essential for ensuring the safe operation of light-water reactors. Here, this study presents a computational model to simulate the intricate three-phase flow dynamics of fuel dispersal following cladding breach, where solid fuel particles and carrier fission gases interact with the surrounding liquid or gas medium. The model utilizes a Eulerian three-fluid framework to depict the phase interaction within interpenetrating continua, treating fragmented fuel as granular material. Closure relationships for bulk-flow properties of the solid phase and detailed interfacial exchange terms in high-temperature, high-speed gas-solid-liquid flow scenarios are proposed. The model was validated using adiabatic experimental results in the literature, focusing on characterizing fuel dispersal behavior during loss-of-coolant accidents conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched probability density functions reported in experimental studies. Key highlights of this study include the theoretical description of a coupled three-phase flow with phase change, interpretation of the high-pressure boundary conditions and insights into the transient behavior of gas-solid-liquid phase dynamics during a simulated high-speed dispersal event.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Electron-Donating para -Substituent (X) Enhances the Water Oxidation Activity of the Catalyst Ru(4'-X-terpyridine)(phenanthroline-SO 3 ) +

Recently, our group has developed a Ru-based water oxidation catalyst (WOC) with pendant sulfonate (1, Ru(4'-X-terpyridine)(phenanthroline-SO 3 )OTf (X = H, 1a) that shows high activity under both sacrificial oxidant (CAN, Ce(NH 4 ) 2 (NO 3 ) 6 , Ce IV ) and electrocatalytic conditions, in both acidic and neutral media. Here, we demonstrate that the functionalization of the 4'-X-terpyridine ligand with an electron-donating substituent X = OEt (1b) makes potentials of Ru II /Ru III redox catalysis more negative, whereas when X = NO 2 (1c) and CF 3 (1d), potentials are more positive. For 1b, full conversion of the sacrificial oxidant Ce IV occurred in 0.4 h (7 h for 1a), with an initial rate of 2.07 μmol O 2 s –1 and a turnover frequency of 7.6 s –1 , which is 30-fold faster than that for 1a at [cat] 0 = 20 μM. Under electrocatalytic conditions, water oxidation by 1b is three times faster than that by the parent catalyst 1a at close to the same potential. Extensive computations have identified differences in the initial PCET steps of the water oxidation by catalysts 1a, 1b, and 1d, and demonstrated the increased probability of the O 2 formation via the oxide relay pathway in the order 1b< 1a < 1d.

14 SOLAR ENERGY↗

DONKEY: A Flexible and Accurate Algorithm for Clustering

We propose an accurate clustering algorithm suitable for the varied and multidimensional data sets that correspond to temporal snapshots from on-the-fly nonadiabatic trajectory-based simulations of photoexcited dynamics. The algorithm approximates the underlying probability density function using variable kernel density estimation, with local maxima corresponding to cluster centers. Each data point is then assigned to one of the maxima by employing a maximization procedure. Finally, clusters artificially separated by minor fluctuations in the probability density are merged. The algorithm does not require parameter tuning, which ensures flexibility and reduces the risk of bias. It is tested on several synthetic data sets, where it consistently outperforms conventional clustering algorithms. As a final example, the algorithm is applied to the excited dynamics of the norbornadiene ⇌ quadricyclane (C 7 H 8 ) molecular photoswitch, demonstrating how distinct reaction pathways can be identified.

algorithms↗

Using intrusive approaches as a step towards accounting for stochasticity in wind turbine design

Current wind turbine design methods require tens of thousands of time-domain simulations and use different random seeds to account for the stochasticity of the environmental conditions. The account of stochasticity is nonintrusive because the sampling method calls a deterministic model multiple times without changing its underlying equations. In this work, we investigate and demonstrate using simple proof of concepts how intrusive approaches can be used to directly account for stochasticity in the equations representing a mechanical system. Our long term goal is to apply such methodology to the design of wind turbines without requiring an excessive number of simulations. Intrusive methods manipulate stochastic variables directly to provide the probability density functions (PDFs) of the states and outputs at any time as functions of the PDFs of the inputs. We illustrate how different methods can be used with a reduced-order model of a wind turbine with one degree of freedom and for linear and nonlinear models. We discuss how the methods can be extended and what it will take to apply them to a level of fidelity similar to current state-of-the-art wind turbine design tools.

17 WIND ENERGY↗

Space station integrated wall design and penetration damage control

The analysis code BUMPER executes a numerical solution to the problem of calculating the probability of no penetration (PNP) of a spacecraft subject to man-made orbital debris or meteoroid impact. The codes were developed on a DEC VAX 11/780 computer that uses the Virtual Memory System (VMS) operating system, which is written in FORTRAN 77 with no VAX extensions. To help illustrate the steps involved, a single sample analysis is performed. The example used is the space station reference configuration. The finite element model (FEM) of this configuration is relatively complex but demonstrates many BUMPER features. The computer tools and guidelines are described for constructing a FEM for the space station under consideration. The methods used to analyze the sensitivity of PNP to variations in design, are described. Ways are suggested for developing contour plots of the sensitivity study data. Additional BUMPER analysis examples are provided, including FEMs, command inputs, and data outputs. The mathematical theory used as the basis for the code is described, and illustrates the data flow within the analysis.

Coronado, A. R.↗

On the Effects of Cloud Water Content on Passive Microwave Snowfall Retrievals

The Bayesian passive microwave retrievals of snowfall often rely on mathematical matching of the observed vectors of brightness temperature with an a priori database of precipitation profiles and their corresponding brightness temperatures. Mathematical proximity does not necessarily lead to consistent retrievals due to limited information content of passive microwave observations. This paper defines imposter (genuine) vectors of brightness temperature as those that are mathematically close but physically inconsistent (consistent) and characterizes them through the Silhouette Coefficient (SC) analysis. The Neyman–Pearson (NP) hypothesis testing is used to separate the imposter and genuine brightness temperatures based on their associated values of cloud ice (IWP) and liquid water path (LWP), given by coincidences of CloudSat Profiling Radar (CPR) and the Global Precipitation Measurement (GPM) Microwave Imager (GMI). The study determines thresholds for IWP and LWP that allow optimal identification of imposter brightness temperatures of non-snowing and snowing clouds, which can mislead the passive microwave retrieval algorithms to falsely detect or miss the snowfall events. It is demonstrated that emission signal of supercooled liquid water in snowing clouds can lead to improved passive microwave retrieval of snowfall and conditioning the retrievals to the cloud IWP and LWP can result in marginal correction of the snowfall detection probability; however, reduce the probability of false alarm by 6%–8% over sea ice and open oceans.

Snowfall↗