Search NASA⌕ Search

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.

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

At least 415 records · Page 23

Monte Carlo calculations of lunar regolith thickness distributions.

It is pointed out that none of the existing models of lunar regolith evolution take into account the relationship between regolith thickness, crater shape, and volume of debris ejected. The results of a Monte Carlo computer simulation of regolith evolution are presented. The simulation was designed to consider the full effect of the buffering regolith through calculation of the amount of debris produced by any given crater as a function of the amount of debris present at the site of the crater at the time of crater formation. The method is essentially an improved version of the Oberbeck and Quaide (1968) model.

Oberbeck, V. R.↗

A Monte Carlo analysis of the effects of instrumentation errors on aircraft parameter identification

An output error estimation algorithm was used to evaluate the effects of both static and dynamic instrumentation errors on the estimation of aircraft stability and control parameters. A Monte Carlo analysis, using simulated cruise flight data, was performed for a high performance military aircraft, a large commercial transport, and a small general-aviation aircraft. The effects of variations in the information content of the flight data, resulting from two different choices of control input maneuvers, were also determined. The results indicate that unmodeled instrumentation errors can cause inaccuracies in the estimated parameters which are comparable to their nominal values. Control input errors and angular accelerometer lags were found to be most significant of the instrumentation errors evaluated, and the perturbations they produce are much larger than those arising from the combined effects of static errors and white noise in the output response measurements.

Bryant, W. H.↗

Monte Carlo simulation of energy deposition by low-energy electrons in molecular hydrogen

A set of detailed atomic cross sections has been used to obtain the spatial deposition of energy by 1-20-eV electrons in molecular hydrogen by a Monte Carlo simulation of the actual trajectories. The energy deposition curve (energy per distance traversed) is quite peaked in the forward direction about the entry point for electrons with energies above the threshold of the electronic states, but the peak decreases and broadens noticeably as the electron energy decreases below 10 eV (threshold for the lowest excitable electronic state of H2). The curve also assumes a very symmetrical shape for energies below 10 eV, indicating the increasing importance of elastic collisions in determining the shape of the curve, although not the mode of energy deposition.

Heaps, M. G.↗

Imaginary time path integral Monte Carlo route to rate coefficients for nonadiabatic barrier crossing

Nonadiabatic transitions are central to many areas of chemical and condensed matter physics, ranging from biological electron transfer to the optical properties of one-dimensional conductors. Here, a path integral Monte Carlo method is used to simulate such transitions, based on the observation that nonadiabatic rate coefficients are often dominated by saddle point trajectories that correspond to an imaginary time. Simple analytic theories can be used to continue these imaginary time correlation functions to determine rate coefficients. The advantages and drawbacks of this approach are discussed.

Wolynes, Peter G.↗

Multi-Model Monte Carlo Estimators for Trajectory Simulation

Predicting landing radius and other quantities of interest (QoI) for entry, descent, andlanding (EDL) applications requires a viable uncertainty propagation method for quantifying the impact of uncertainties in aerodynamics, atmosphere, mass properties, etc. While standard Monte Carlo (MC) simulation is the de facto standard for producing robust and unbiasedstatistical estimators, it is often infeasible for expensive, high-fidelity models. Low-fidelity models are commonly constructed to replace the high-fidelity model in MC simulation for computational speedup, but at the expense of accuracy and unbiasedness. Emerging multi-model MC methods are bridging this gap by combining predictions from two or more modelsof varying fidelity and computational cost for efficient and unbiased uncertainty propagation.This works establishes a proof of concept for using multi-model MC to increase the speed and precision of trajectory simulation for EDL. It is shown that combining a high-fidelity EDL model with low-fidelity models (e.g., data-driven, reduced physics) in this manner has the potential to yield significant efficiency and accuracy gains for certain EDL QoIs versusa standard MC approach. Moreover, the unbiasedness of multi-model MC predictions ishighlighted by showing increased accuracy versus an approach that leverages a low-fidelity model alone.

James E Warner↗

Direct Monte Carlo Simulations of Hypersonic Low-Density Flows about an ASTV Including Wake Structure

Results of a numerical study concerning flow past a 70-deg blunted cone in hypersonic low-density flow environments are presented using the direct simulation Monte-Carlo method. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Results indicate that a stable vortex forms in the near wake at and below a freestream Knudsen number (based on cone diameter) of 0.01 and the size of the vortex increases with decreasing Knudsen number. The base region of the flow remains in thermal nonequilibrium for all cases considered herein.

Dogra, V. K.↗

Monte Carlo analysis of inaccuracies in estimated aircraft parameters caused by unmodeled flight instrumentation errors

An output error estimation algorithm was used to evaluate the effects of both static and dynamic instrumentation errors on the estimation of aircraft stability and control parameters. A Monte Carlo error analysis, using simulated cruise flight data, was performed for a high-performance military aircraft, a large commercial transport, and a small general aviation aircraft. The results indicate that unmodeled instrumentation errors can cause inaccuracies in the estimated parameters which are comparable to their nominal values. However, the corresponding perturbations to the estimated output response trajectories and characteristics equation pole locations appear to be relatively small. Control input errors and dynamic lags were found to be in the most significant of the error sources evaluated.

Hodge, W. F.↗

Monte Carlo exploration of Mikheyev-Smirnov-Wolfenstein solutions to the solar neutrino problem

The paper explores the impact of astrophysical uncertainties on the Mikheyev-Smirnov-Wolfenstein (MSW) solution by calculating the allowed MSW solutions for 1000 different solar models with a Monte Carlo selection of solar model input parameters, assuming a full three-family MSW mixing. Applications are made to the chlorine, gallium, Kamiokande, and Borexino experiments. The initial GALLEX result limits the mixing parameters to the upper diagonal and the vertical regions of the MSW triangle. The expected event rates in the Borexino experiment are also calculated, assuming the MSW solutions implied by GALLEX.

Shi, X.↗

Comparison of quick-look plume heating calculations and Monte Carlo direct simulation

A quick-look method has been developed for estimating convective heat flux due to plume impingement in the transition flow regime on a concave surface surrounding a rocket nozzle. Comparison with the flowfield and heat fluxes computed by the Monte Carlo Direct Simulation Method show the quick-look method to be conservative. Assumptions regarding the nature of the flowfield based on engineering judgment were shown to be essentially valid. Further, the assumption of free molecular impingement was shown to be non-conservative for portions of the surface. Potential refinements of the quick-look method are discussed.

Guernsey, C. S.↗

Monte Carlo studies of ocean wind vector measurements by SCATT: Objective criteria and maximum likelihood estimates for removal of aliases, and effects of cell size on accuracy of vector winds

The scatterometer on the National Oceanic Satellite System (NOSS) is studied by means of Monte Carlo techniques so as to determine the effect of two additional antennas for alias (or ambiguity) removal by means of an objective criteria technique and a normalized maximum likelihood estimator. Cells nominally 10 km by 10 km, 10 km by 50 km, and 50 km by 50 km are simulated for winds of 4, 8, 12 and 24 m/s and incidence angles of 29, 39, 47, and 53.5 deg for 15 deg changes in direction. The normalized maximum likelihood estimate (MLE) is correct a large part of the time, but the objective criterion technique is recommended as a reserve, and more quickly computed, procedure. Both methods for alias removal depend on the differences in the present model function at upwind and downwind. For 10 km by 10 km cells, it is found that the MLE method introduces a correlation between wind speed errors and aspect angle (wind direction) errors that can be as high as 0.8 or 0.9 and that the wind direction errors are unacceptably large, compared to those obtained for the SASS for similar assumptions.

Pierson, W. J.↗

Monte Carlo Tree Search Methods for the Earth-Observing Satellite Scheduling Problem

This work explores on-board planning for the single spacecraft, multiple ground station Earth-observing satellite scheduling problem through artificial neural network function approximation of state–action value estimates generated by Monte Carlo tree search (MCTS). An extensive hyperparameter search is conducted for MCTS on the basis of performance, safety, and downlink opportunity utilization to determine the best hyperparameter combination for data generation. A hyperparameter search is also conducted on neural network architectures. The learned behavior of each network is explored, and each network architecture’s robustness to orbits and epochs outside of the training distributions is investigated. Furthermore, each algorithm is compared with a genetic algorithm, which serves to provide a baseline for optimality. MCTS is shown to compute near-optimal solutions in comparison to the genetic algorithm. The state–action value networks are shown to match or exceed the performance of MCTS in six orders of magnitude less execution time, showing promise for execution on board spacecraft.

Adam P. Herrmann↗

Monte Carlo simulation of charged particles upstream of the earth's bow shock

A simple model of the earth's bow shock is presented and compared with ISEE-1 observations of diffuse ions upstream of the shock. Ion behavior is modeled in a test particle simulation using Monte Carlo techniques for a plane parallel shock at which cold solar wind ions are heated sufficiently to allow diffusion back upstream. Fitting of the model to the observations indicates that the diffusion coefficient must increase faster than linearly with velocity, and provides an upper limit to the acceleration efficiency defined as the proportion of incoming solar wind ions accelerated and observed upstream. Results suggest that astrophysical shocks are capable of accelerating ions directly from the thermal plasma without the need for a distinct preenergizing injection stage.

Ellison, D. C.↗

A Monte Carlo model of polarized thermal emission from particulate planetary surfaces

Direct emission from individual grains and multiple scattering between regolith particles are encompassed by the present model of particulate planetary surface thermal emission, whose randomly positioned spherical grains are large relative to the emission's wavelength scale. A Monte Carlo ray-tracing method is used to calculate the spectral and directional emissivity of the surface and the polarization of the emitted radiation, for 7-16 micron wavelengths. The effects of roughness at the scale of individual grains and scattering are separated to elucidate how each affects the emitted radiation. Implications of these results for planetary remote sensing are discussed.

Henderson, Bradley G.↗

Three-dimensional hypersonic rarefied flow calculations using direct simulation Monte Carlo method

A summary of three-dimensional simulations on the hypersonic rarefied flows in an effort to understand the highly nonequilibrium flows about space vehicles entering the Earth's atmosphere for a realistic estimation of the aerothermal loads is presented. Calculations are performed using the direct simulation Monte Carlo method with a five-species reacting gas model, which accounts for rotational and vibrational internal energies. Results are obtained for the external flows about various bodies in the transitional flow regime. For the cases considered, convective heating, flowfield structure and overall aerodynamic coefficients are presented and comparisons are made with the available experimental data. The agreement between the calculated and measured results are very good.

Celenligil, M. Cevdet↗

Monte Carlo Simulations of Background Spectra in Integral Imager Detectors

Predictions of the expected gamma-ray backgrounds in the ISGRI (CdTe) and PiCsIT (Csl) detectors on INTEGRAL due to cosmic-ray interactions and the diffuse gamma-ray background have been made using a coupled set of Monte Carlo radiation transport codes (HETC, FLUKA, EGS4, and MORSE) and a detailed, 3-D mass model of the spacecraft and detector assemblies. The simulations include both the prompt background component from induced hadronic and electromagnetic cascades and the delayed component due to emissions from induced radioactivity. Background spectra have been obtained with and without the use of active (BGO) shielding and charged particle rejection to evaluate the effectiveness of anticoincidence counting on background rejection.

Armstrong, T. W.↗

A Monte Carlo error analysis program for near-Mars, finite-burn, orbital transfer maneuvers

A computer program was developed which performs an error analysis of a minimum-fuel, finite-thrust, transfer maneuver between two Keplerian orbits in the vicinity of Mars. The method of analysis is the Monte Carlo approach where each off-nominal initial orbit is targeted to the desired final orbit. The errors in the initial orbit are described by two covariance matrices of state deviations and tracking errors. The function of the program is to relate these errors to the resulting errors in the final orbit. The equations of motion for the transfer trajectory are those of a spacecraft maneuvering with constant thrust and mass-flow rate in the neighborhood of a single body. The thrust vector is allowed to rotate in a plane with a constant pitch rate. The transfer trajectory is characterized by six control parameters and the final orbit is defined, or partially defined, by the desired target parameters. The program is applicable to the deboost maneuver (hyperbola to ellipse), orbital trim maneuver (ellipse to ellipse), fly-by maneuver (hyperbola to hyperbola), escape maneuvers (ellipse to hyperbola), and deorbit maneuver.

Green, R. N.↗

Evaluation of Lightning Incidence to Elements of a Complex Structure: A Monte Carlo Approach

There are complex structures for which the installation and positioning of the lightning protection system (LPS) cannot be done using the lightning protection standard guidelines. As a result, there are some "unprotected" or "exposed" areas. In an effort to quantify the lightning threat to these areas, a Monte Carlo statistical tool has been developed. This statistical tool uses two random number generators: a uniform distribution to generate origins of downward propagating leaders and a lognormal distribution to generate returns stroke peak currents. Downward leaders propagate vertically downward and their striking distances are defined by the polarity and peak current. Following the electrogeometrical concept, we assume that the leader attaches to the closest object within its striking distance. The statistical analysis is run for 10,000 years with an assumed ground flash density and peak current distributions, and the output of the program is the probability of direct attachment to objects of interest with its corresponding peak current distribution.

Mata, Carlos T.↗

Evaluation of Lightning Incidence to Elements of a Complex Structure: A Monte Carlo Approach

There are complex structures for which the installation and positioning of the lightning protection system (LPS) cannot be done using the lightning protection standard guidelines. As a result, there are some "unprotected" or "exposed" areas. In an effort to quantify the lightning threat to these areas, a Monte Carlo statistical tool has been developed. This statistical tool uses two random number generators: a uniform distribution to generate the origin of downward propagating leaders and a lognormal distribution to generate the corresponding returns stroke peak currents. Downward leaders propagate vertically downward and their striking distances are defined by the polarity and peak current. Following the electrogeometrical concept, we assume that the leader attaches to the closest object within its striking distance. The statistical analysis is run for N number of years with an assumed ground flash density and the output of the program is the probability of direct attachment to objects of interest with its corresponding peak current distribution.

Mata, Carlos T.↗