Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo estimators”

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 235 records · Page 13

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan‐1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017, https://doi.org/10.1126/sciadv.1701471). The model reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U(sub 0) = 0.5 eV, U(sub W) = 0.078–0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

Tucker, O. J.↗

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan-1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017). Themodel reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U(o) = 0.5eV, U(w) = 0.078 – 0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

O. J. Tucker↗

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan‐1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017, https://doi.org/10.1126/sciadv.1701471). The model reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U0 = 0.5 eV, UW = 0.078–0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

O J Tucker↗

Robust Maneuvering Envelope Estimation Based on Reachability Analysis in an Optimal Control Formulation

This paper discusses an algorithm for estimating the safe maneuvering envelope of damaged aircraft. The algorithm performs a robust reachability analysis through an optimal control formulation while making use of time scale separation and taking into account uncertainties in the aerodynamic derivatives. Starting with an optimal control formulation, the optimization problem can be rewritten as a Hamilton- Jacobi-Bellman equation. This equation can be solved by level set methods. This approach has been applied on an aircraft example involving structural airframe damage. Monte Carlo validation tests have confirmed that this approach is successful in estimating the safe maneuvering envelope for damaged aircraft.

optimal control↗

Disentangling AGN and Star Formation in Soft X-Rays

We have explored the interplay of star formation and active galactic nucleus (AGN) activity in soft X-rays (0.5-2 keV) in two samples of Seyfert 2 galaxies (Sy2s). Using a combination of low-resolution CCD spectra from Chandra and XMM-Newton, we modeled the soft emission of 34 Sy2s using power-law and thermal models. For the 11 sources with high signal-to-noise Chandra imaging of the diffuse host galaxy emission, we estimate the luminosity due to star formation by removing the AGN, fitting the residual emission. The AGN and star formation contributions to the soft X-ray luminosity (i.e., L(sub x,AGN) and L(sub x,SF)) for the remaining 24 Sy2s were estimated from the power-law and thermal luminosities derived from spectral fitting. These luminosities were scaled based on a template derived from XSINGS analysis of normal star-forming galaxies. To account for errors in the luminosities derived from spectral fitting and the spread in the scaling factor, we estimated L(sub x,AGN) and L(sub x,SF))from Monte Carlo simulations. These simulated luminosities agree with L(sub x,AGN) and L(sub x,SF) derived from Chandra imaging analysis within a 3sigma confidence level. Using the infrared [Ne ii]12.8 micron and [O iv]26 micron lines as a proxy of star formation and AGN activity, respectively, we independently disentangle the contributions of these two processes to the total soft X-ray emission. This decomposition generally agrees with L(sub x,SF) and L(sub x,AGN) at the 3 sigma level. In the absence of resolvable nuclear emission, our decomposition method provides a reasonable estimate of emission due to star formation in galaxies hosting type 2 AGNs.

Disentangling↗

Gyro and accelerometer failure detection and identification in redundant sensor systems

Algorithms for failure detection and identification for redundant noncolinear arrays of single degree of freedom gyros and accelerometers are described. These algorithms are optimum in the sense that detection occurs as soon as it is no longer possible to account for the instrument outputs as the outputs of good instruments operating within their noise tolerances, and identification occurs as soon as it is true that only a particular instrument failure could account for the actual instrument outputs within the noise tolerance of good instruments. An estimation algorithm is described which minimizes the maximum possible estimation error magnitude for the given set of instrument outputs. Monte Carlo simulation results are presented for the application of the algorithms to an inertial reference unit consisting of six gyros and six accelerometers in two alternate configurations.

Potter, J. E.↗

Minimax failure detection and identification in redundant gyro and accelerometer systems.

Minimax algorithms for failure detection and identification for redundant noncolinear arrays of single-degree-of-freedom gyros and accelerometers are described. These algorithms are optimum in the sense that detection occurs as soon as it is no longer possible to account for the instrument outputs as the outputs of good instruments operating within their noise tolerances, and identification occurs as soon as it is true that only a particular instrument failure could account for the actual instrument outputs within the noise tolerance of good instruments. An estimation algorithm is described which minimizes the maximum possible estimation error magnitude for the given set of instrument outputs. Monte Carlo simulation results are presented for the application of the algorithms to an inertial reference unit consisting of six gyros and six accelerometers in two alternate configurations.

Potter, J. E.↗

Algorithms for adaptive stochastic control for a class of linear systems

Control of linear, discrete time, stochastic systems with unknown control gain parameters is discussed. Two suboptimal adaptive control schemes are derived: one is based on underestimating future control and the other is based on overestimating future control. Both schemes require little on-line computation and incorporate in their control laws some information on estimation errors. The performance of these laws is studied by Monte Carlo simulations on a computer. Two single input, third order systems are considered, one stable and the other unstable, and the performance of the two adaptive control schemes is compared with that of the scheme based on enforced certainty equivalence and the scheme where the control gain parameters are known.

Toda, M.↗

Proton-tissue dose buildup factors

When an object is exposed to external radiation, the dose field within the object is a complicated function of the character of the external radiation, the shape of the object (including orientation), and the object's material composition. In this note, the dose conversion factors for protons in tissue are represented using buildup factors. A parametric form for the buildup factors is obtained. The values of the parameters are derived from Monte Carlo calculations of various authors. All the necessary information to estimate nuclear reaction effects in proton irration of convex objects of arbitrary shape is included.

Wilson, J. W.↗

Heating of the moon by heterogeneous accretion

Monte Carlo models of lunar accretion are developed by randomly sampling mass and velocity distributions of objects in the solar nebula in circumsolar and circumterrestrial orbits. The thermal effect of an impact is expressed as the volume of melt produced; the melt volumes produced by each Monte Carlo case are sorted into regions to produce percentage melt estimates as a function of radius. One-hundred cases are calculated to obtain a range of plausible thermal histories for lunar accretion. The average melt percentage estimates show substantial deep melting, but the spread of melt percentage values about the mean from the random sampling procedure is wide enough to include 0% melting for radii below 100 km and 100% melting for radii below 500 km.

Ransford, G. A.↗

A minicomputer based software system for the selection of optimal subsets of Thematic Mapper channels

A software system has been developed and implemented on a minicomputer for feature selection based on two inter-dependent methods. The first is an enhancement of the traditional approach based on optimizing interclass average separabilities. The second is based on a Monte Carlo simulation of multispectral data and machine classification with subsequent estimation of classification accuracy as a function of channel subset. The two methods are mutually supportive - the first allows rapid screening whereas the second is based on the more solid theoretical foundation of maximizing classification accuracy.

Card, D. H.↗

Planetary accretion rates - Analytical derivation

New analytical formulae that are valid even in cases where Keplerian shear determines the approach velocity in planetary accretion rates are presently shown to accurately reproduce Monte Carlo study results. The formulae yield accretion rates over a great variety of planet formation-conducive conditions. Since the absolute values of impact rates yielded by the formulae were expected to have factor-of-2-to-3 uncertainties, the good agreement with Monte Carlo experiments obtained suggests the reasonableness of many of the estimates adopted, or the mutual cancelling-out of various uncertainties.

Greenberg, Richard↗

Effects of location uncertainties on the observed distribution of gamma-ray bursts detected by BATSE

The Burst and Transient Source Experiment (BATSE) on NASA's Compton Gamma-Ray Observatory has shown that the sources of gamma-ray bursts have an angular distribution consistent with isotropy, yet the bursts are spatially inhomogeneous. Detection of a statistically significant deviation from isotropy in the burst distribution would provide a clue to the distance of the burst sources. BATSE's ability to detect slight deviations of the angular distribution from isotropy is, in part, dependent on the accuracy of the computed burst locations. Using bursts whose locations have been determined accurately and independently, we estimate BATSE's location uncertainties and incorporate them into 3D Monte Carlo simulations of the burst distributions to show that, despite individual burst location uncertainties of about 5 deg, the ability of BATSE to place stringent limits on the anisotropy present in the angular distribution of the gamma-ray bursts is not significantly degraded. Observational effects of possible galactic halo distributions of the burst sources are noted.

Horack, J. M.↗

Impact and volcanism - A momentum scaling law for erosion

The controversy that raged through the 1950s to 1970s over similarities and differences between meteorite impact and volcanic processes is revisited. We propose that there are quantitative similarities in erosion caused by high-speed ejecta produced by either impacts or volcanic processes. Field and petrographic data from the Manicouagan impact crater, Canada, are used to demonstrate that, during the emplacement of the impact melt sheet, erosion occurred at a rate of 2562 kg/sq m per s. Field data for the Mount St. Helens lateral blast of May 18, 1980, suggest an erosion rate of 21 kg/sq m per s, and field data for a small pyroclastic flow on August 7, 1980, suggest an erosion rate of about 14 kg /sq m per s. It is proposed that these three rates were determined dominantly by the momentum of the ejecta, and a quantitative formulation based on lofting theory is given. A new application of the Monte Carlo approach to analysis provides minimum, most likely, and maximum estimates for both the field and the theoretical analyses. The substantial erosion that occurs by fast moving flows results in mixing of stratigraphic components over large distances and to very fine scales.

Simonds, Charles H.↗

Crack Growth Properties of Sealing Glasses

The crack growth properties of several sealing glasses were measured using constant stress rate testing in 2% and 95% RH (relative humidity). Crack growth parameters measured in high humidity are systematically smaller (n and B) than those measured in low humidity, and velocities for dry environments are approx. 100x lower than for wet environments. The crack velocity is very sensitivity to small changes in RH at low RH. Confidence intervals on parameters that were estimated from propagation of errors were comparable to those from Monte Carlo simulation.

Salem, Jonathan A.↗

Probabilistic Analysis of a Composite Crew Module

An approach for conducting reliability-based analysis (RBA) of a Composite Crew Module (CCM) is presented. The goal is to identify and quantify the benefits of probabilistic design methods for the CCM and future space vehicles. The coarse finite element model from a previous NASA Engineering and Safety Center (NESC) project is used as the baseline deterministic analysis model to evaluate the performance of the CCM using a strength-based failure index. The first step in the probabilistic analysis process is the determination of the uncertainty distributions for key parameters in the model. Analytical data from water landing simulations are used to develop an uncertainty distribution, but such data were unavailable for other load cases. The uncertainty distributions for the other load scale factors and the strength allowables are generated based on assumed coefficients of variation. Probability of first-ply failure is estimated using three methods: the first order reliability method (FORM), Monte Carlo simulation, and conditional sampling. Results for the three methods were consistent. The reliability is shown to be driven by first ply failure in one region of the CCM at the high altitude abort load set. The final predicted probability of failure is on the order of 10-11 due to the conservative nature of the factors of safety on the deterministic loads.

FROM↗

Three-Dimensional Finite Element Ablative Thermal Response and Thermostructural Design of Thermal Protection Systems

A finite element ablation and thermal response program is presented for simulation of three-dimensional transient thermostructural analysis. The three-dimensional governing differential equations and finite element formulation are summarized. A novel probabilistic design methodology for thermal protection systems is presented. The design methodology is an eight step process beginning with a parameter sensitivity study and is followed by a deterministic analysis whereby an optimum design can determined. The design process concludes with a Monte Carlo simulation where the probabilities of exceeding design specifications are estimated. The design methodology is demonstrated by applying the methodology to the carbon phenolic compression pads of the Crew Exploration Vehicle. The maximum allowed values of bondline temperature and tensile stress are used as the design specifications in this study.

Dec, John A.↗

A Quantitative Risk-Benefit Analysis of Prophylactic Surgery Prior to Extended-Duration Spaceflight

INTRODUCTION: Among otherwise healthy astronauts undertaking deep space missions, the risks for acute appendicitis (AA) and cholecystitis (AC) are not zero. If these conditions were to occur during spaceflight they may require surgery for definitive care. The proposed study quantifies and compares the risks of developing de novo AA and AC in-flight to the surgical risks of prophylactic laparoscopic appendectomy (LA) and cholecystectomy (LC) using NASA's Integrated Medical Model (IMM). METHODS: The IMM is a Monte Carlo simulation that forecasts medical events during spaceflight missions and estimates the impact of these medical events on crew health. In this study, four Design Reference Missions (DRMs) were created to assess the probability of an astronaut developing in-flight small-bowel obstruction (SBO) following prophylactic 1) LA, 2) LC, 3) LA and LC, or 4) neither surgery (SR# S-20160407-351). Model inputs were drawn from a large, population-based 2011 Swedish study that examined the incidence and risks of post-operative SBO over a 5-year follow-up period. The study group included 1,152 patients who underwent LA, and 16,371 who underwent LC. RESULTS: Preliminary results indicate that prophylactic LA may yield higher mission risks than the control DRM. Complete analyses are pending and will be subsequently available. DISCUSSION: The risk versus benefits of prophylactic surgery in astronauts to decrease the probability of acute surgical events during spaceflight has only been qualitatively examined in prior studies. Within the assumptions and limitations of the IMM, this work provides the first quantitative guidance that has previously been lacking to this important question for future deep space exploration missions.

Carroll, Danielle↗