SEARCH · Search NASA

Results for “monte carlo”

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.

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS1985

A three-dimensional Monte Carlo calculation of the photon initiated showers and Kiel result

The Kiel experimental results indicate an existence of the ultra high-energy gamma-rays coming from Cyg. X-3. However the result indicates that the number of the muons included in the photon initiated shower is the same as the number included in the proton initiated showers. According to our Monte Carlo calculation as shown in the graph of underpart, the number of muons included in the photon initiated showers should be less than 1/15 of the photon's. The previous simulation was made under one dimensional approximation. This time the result of three dimensional calculation is reported.

Okada, A.↗

NASA NTRS1999

Improving Pixel Level Cloud Optical Property Retrieval using Monte Carlo Simulations

The accurate pixel-by-pixel retrieval of cloud optical properties from space is influenced by radiative smoothing due to high order photon scattering and radiative roughening due to low order scattering events. Both are caused by cloud heterogeneity and the three-dimensional nature of radiative transfer and can be studied with the aid of computer simulations. We use Monte Carlo simulations on variable 1-D and 2-D model cloud fields to seek for dependencies of smoothing and roughening phenomena on single scattering albedo, solar zenith angle, and cloud characteristics. The results are discussed in the context of high resolution satellite (such as Landsat) retrieval applications. The current work extends the investigation on the inverse NIPA (Non-local Independent Pixel Approximation) as a tool for removing smoothing and improving retrievals of cloud optical depth. This is accomplished by: (1) Delineating the limits of NIPA applicability; (2) Exploring NIPA parameter dependences on cloud macrostructural features, such as mean cloud optical depth and geometrical thickness, degree of extinction and cloud top height variability. We also compare parameter values from empirical and theoretical considerations; (3) Examining the differences between applying NIPA on radiation quantities vs direct application on optical properties; (4) Studying the radiation budget importance of the NIPA corrections as a function of scale. Finally, we discuss fundamental adjustments that need to be considered for successful radiance inversion at non-conservative wavelengths and oblique Sun angles. These adjustments are necessary to remove roughening signatures which become more prominent with increasing absorption and solar zenith angle.

Oreopoulos, Lazaros↗

NASA NTRS1974

A Monte Carlo investigation of experimental data requirements for fitting polynomial functions

This report examines the extent to which sample size affects the accuracy of a low order polynomial approximation of an experimentally observed quantity and establishes a trend toward improvement in the accuracy of the approximation as a function of sample size. The task is made possible through a simulated analysis carried out by the Monte Carlo method, in which data are generated by using several transcendental or algebraic functions as models. Contaminated data of varying amounts are fitted to linear quadratic or cubic polynomials, and the behavior of the mean-squared error of the residual variance is determined as a function of sample size. Results indicate that the effect of the size of the sample is significant only for relatively small sample sizes and diminishes drastically for moderate and large amounts of experimental data.

Canavos, G. C.↗

NASA NTRS2014

A Monte Carlo Analysis of the Thrust Imbalance for the RSRMV Booster During Both the Ignition Transient and Steady State Operation

This paper presents the results of statistical analyses performed to predict the thrust imbalance between two solid rocket motor boosters to be used on the Space Launch System (SLS) vehicle. Two legacy internal ballistics codes developed for the Space Shuttle program were coupled with a Monte Carlo analysis code to determine a thrust imbalance envelope for the SLS vehicle based on the performance of 1000 motor pairs. Thirty three variables which could impact the performance of the motors during the ignition transient and thirty eight variables which could impact the performance of the motors during steady state operation of the motor were identified and treated as statistical variables for the analyses. The effects of motor to motor variation as well as variations between motors of a single pair were included in the analyses. The statistical variations of the variables were defined based on data provided by NASA's Marshall Space Flight Center for the upgraded five segment booster and from the Space Shuttle booster when appropriate. The results obtained for the statistical envelope are compared with the design specification thrust imbalance limits for the SLS launch vehicle

Foster, Winfred A., Jr.↗

NASA NTRS1993

Monte Carlo and experimental studies of nozzle flow in a low-power hydrogen arcjet

The flow in the nozzle of a small hydrogen arcjet intended for low-thrust propulsion is studied without arc ignition using numerical and experimental techniques. The flow conditions in the thruster indicate that low-density nonequilibrium effects will significantly influence the flow. Therefore, the numerical analysis is undertaken using a Monte Carlo approach. The experimental studies employ laser Raman scattering. Comparisons of the measured and computed results are made for total number density, rotational temperature, and for the number density of the first rotational level. The numerical results are found to be quite sensitive to the rotational relaxation rate, and a strong degree of thermal nonequilibrium is observed at the exit plane of the nozzle. Successful comparisons between experiment and analysis permit estimation of the rotational relaxation rate for hydrogen.

Boyd, Iain D.↗

NASA NTRS2013

Uncertainty Determination for Aeroheating in Uranus and Saturn Probe Entries by the Monte Carlo Method

The 2013-2022 Decaedal survey for planetary exploration has identified probe missions to Uranus and Saturn as high priorities. This work endeavors to examine the uncertainty for determining aeroheating in such entry environments. Representative entry trajectories are constructed using the TRAJ software. Flowfields at selected points on the trajectories are then computed using the Data Parallel Line Relaxation (DPLR) Computational Fluid Dynamics Code. A Monte Carlo study is performed on the DPLR input parameters to determine the uncertainty in the predicted aeroheating, and correlation coefficients are examined to identify which input parameters show the most influence on the uncertainty. A review of the present best practices for input parameters (e.g. transport coefficient and vibrational relaxation time) is also conducted. It is found that the 2(sigma) - uncertainty for heating on Uranus entry is no more than 2.1%, assuming an equilibrium catalytic wall, with the uncertainty being determined primarily by diffusion and H(sub 2) recombination rate within the boundary layer. However, if the wall is assumed to be partially or non-catalytic, this uncertainty may increase to as large as 18%. The catalytic wall model can contribute over 3x change in heat flux and a 20% variation in film coefficient. Therefore, coupled material response/fluid dynamic models are recommended for this problem. It was also found that much of this variability is artificially suppressed when a constant Schmidt number approach is implemented. Because the boundary layer is reacting, it is necessary to employ self-consistent effective binary diffusion to obtain a correct thermal transport solution. For Saturn entries, the 2(sigma) - uncertainty for convective heating was less than 3.7%. The major uncertainty driver was dependent on shock temperature/velocity, changing from boundary layer thermal conductivity to diffusivity and then to shock layer ionization rate as velocity increases. While radiative heating for Uranus entry was negligible, the nominal solution for Saturn computed up to 20% radiative heating at the highest velocity examined. The radiative heating followed a non-normal distribution, with up to a 3x variation in magnitude. This uncertainty is driven by the H(sub 2) dissociation rate, as H(sub 2) that persists in the hot non-equilibrium zone contributes significantly to radiation.

Palmer, Grant↗

NASA NTRS2014

A Monte Carlo Analysis of the Thrust Imbalance for the Space Launch System Booster During Both the Ignition Transient and Steady State Operation

Foster, Winfred A., Jr.↗

NASA NTRS1996

Three-Dimensional Electromagnetic Monte Carlo Particle-in-Cell Simulations of Critical Ionization Velocity Experiments in Space

Although the existence of the critical ionization velocity (CIV) is known from laboratory experiments, no agreement has been reached as to whether CIV exists in the natural space environment. In this paper we move towards more realistic models of CIV and present the first fully three-dimensional, electromagnetic particle-in-cell Monte-Carlo collision (PIC-MCC) simulations of typical space-based CIV experiments. In our model, the released neutral gas is taken to be a spherical cloud traveling across a magnetized ambient plasma. Simulations are performed for neutral clouds with various sizes and densities. The effects of the cloud parameters on ionization yield, wave energy growth, electron heating, momentum coupling, and the three-dimensional structure of the newly ionized plasma are discussed. The simulations suggest that the quantitative characteristics of momentum transfers among the ion beam, neutral cloud, and plasma waves is the key indicator of whether CIV can occur in space. The missing factors in space-based CIV experiments may be the conditions necessary for a continuous enhancement of the beam ion momentum. For a typical shaped charge release experiment, favorable CIV conditions may exist only in a very narrow, intermediate spatial region some distance from the release point due to the effects of the cloud density and size. When CIV does occur, the newly ionized plasma from the cloud forms a very complex structure due to the combined forces from the geomagnetic field, the motion induced emf, and the polarization. Hence the detection of CIV also critically depends on the sensor location.

Wang, J.↗

NASA NTRS1988

Monte Carlo modeling of large-scale ion-conic generation

Cyclotron resonance with observed electric field fluctuations is demonstrated to be responsible for production of the oxygen ion conics that are observed by the Dynamics Explorer 1 satellite in the central plasma sheet region of the earth's auroral zone. The ion velocity distribution is described by a quasi-linear diffusion equation which is solved using the Monte Carlo technique. The acceleration produced by the observed wave spectrum agrees well with the ion observations, in both form and magnitude.

Retterer, J. M.↗

NASA NTRS2014

James Webb Space Telescope Initial Mid-Course Correction Monte Carlo Implementation using Task Parallelism

The James Webb Space Telescope will be launched into a highly elliptical orbit that does not possess sufficient energy to achieve a proper Sun-Earth L2 libration point orbit. Three mid-course correction (MCC) maneuvers are planned to rectify the energy deficit: MCC-1a, MCC-1b, and MCC-2. To validate the propellant budget and trajectory design methods, a set of Monte Carlo analyses that incorporate MCC maneuver modeling and execution are employed. The first analysis focuses on the effects of launch vehicle injection errors on the magnitude of MCC-1a. The second on the spread of potential V based on the performance of the propulsion system as applied to all three MCC maneuvers. The final highlights the slight, but notable, contribution of the attitude thrusters during each MCC maneuver. Given the possible variations in these three scenarios, the trajectory design methods are determined to be robust to errors in the modeling of the flight system.

Libration Point Orbit↗

NASA NTRS2021

Lila: Optimal Dispatching in Probabilistic Temporal Networks using Monte Carlo Tree Search

Executing a Probabilistic Simple Temporal Network (PSTN) amounts at scheduling, i.e. \textit{dispatch}, a set of events under time uncertainty. This constitutes a NP-hard online optimization problem. The right execution time must be dynamically assigned to each event of the PSTN such that the temporal constraints are met, whereas activity durations are progressively observed as the execution unfolds. We propose a dispatching algorithm based on Monte Carlo Tree Search, called Lila, with the following characteristics: (i) it is an anytime algorithm, both offline and online, proven asymptotically optimal; (ii) it returns the current probability of success, either before or at any moment during operations; (iii) it handles any possible continuous or discrete, even non-parametric, probability distributions, as well as inter-dependencies between random variables, exogenous and endogenous uncertainty; and (iv) can be easily extended to handle probabilistic external events, PSTNs with resources, PSTNs with cutoff times and precondition chains, etc. Lila is universal in the sense that it can handle any dispatching protocol, simply by specifying it to the algorithm. It has the unlimited flexibility offered by the simulation paradigm, whilst it asymptotically converges to optimal decisions and/or robustness approximations.

Chien, Steve A.↗

NASA NTRS2014

James Webb Space Telescope Initial Mid-Course Correction Monte Carlo Implementation using Task Parallelism

The James Webb Space Telescope will be launched into a highly elliptical orbit that does not possess sufficient energy to achieve a proper Sun-Earth/Moon L2 libration point orbit. Three mid-course correction (MCC) maneuvers are planned to rectify the energy deficit: MCC-1a, MCC-1b, and MCC-2. To validate the propellant budget and trajectory design methods, a set of Monte Carlo analyses that incorporate MCC maneuver modeling and execution are employed. The first analysis focuses on the effects of launch vehicle injection errors on the magnitude of MCC-1a. The second on the spread of potential V based on the performance of the propulsion system as applied to all three MCC maneuvers. The final highlights the slight, but notable, contribution of the attitude thrusters during each MCC maneuver. Given the possible variations in these three scenarios, the trajectory design methods are determined to be robust to errors in the modeling of the flight system.

James Webb Space Telescope↗

NASA NTRS2017

MCRadar: A Monte Carlo Solver for Cloud and Precipitation Radar

Multiple scattering produces anomalous echoes in observed radar profiles that cannot be explained by other phenomena. These effects are most obvious for spaceborne platforms, at shorter wavelengths, and in convection, as multiple scattering is governed by antenna beam width, optical depth, and albe do; however, multiple scattering has been observed in a range of precipitating conditions. To account for the effects of multiple scattering, Monte Carlo integration is employed in a flexible framework to enable arbitrary radar configurations with finite Gaussian beams for three-dimensional atmospheric scenarios. The three-dimensional nature of the model, coupled with the finite antenna response, also allows for consideration of nonuniform beam filling effects that often coincide with multiple scattering. Examples of airborne and spaceborne radars at various wavelengths are used to illustrate the effects of multiple scattering and nonuniform beam filling. The MCRadar code is currently available in the development version of ARTS.

Adams, Ian Stuart↗

NASA NTRS2019

Assessing GEO and LEO Repeating Conjunctions Using High Fidelity Brute Force Monte Carlo Simulations

Probability of collision (P(sub c)) estimates for Earth-orbiting satellites typically assume a temporally-isolated conjunction event. However, under certain conditions two objects may experience multiple high-risk close approach events over the course of hours or days. In these repeating conjunction cases, the P(sub c) accumulates as each successive encounter occurs. The NASA Conjunction Assessment Risk Analysis team has updated its “brute force Monte Carlo” (BFMC) software to estimate such accumulating P(sub c) values for repeating conjunctions. This study describes the updated BFMC algorithm and discusses the implications for conjunction risk assessment.

Baars, Luis↗

NASA NTRS2019

Assessing Geo and Leo Repeating Conjunctions Using High Fidelity Brute Force Monte Carlo Simulations

Baars, Luis↗

NASA NTRS2018

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.↗

NASA NTRS2018

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↗

NASA NTRS2018

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↗