Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Simulation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Monte Carlo simulation of astronaut light flash observations on Apollo and Skylab missions

Monte Carlo calculations are performed to study the relation between cosmic ray nuclei and visual flashes observed by astronauts, in space. A satisfactory overall fit to all available data excluding the Skylab observations in the South Atlantic Anomaly (SAA) has been achieved by varying three visual system parameters: effective retinal thickness, minimum projected track length inside retina and minimum energy loss rate. All light flash data from space excluding the SAA observations can be understood fully upon the basis of the primary cosmic ray charge and energy spectra together with the fitted values of the three visual system parameters. Monte Carlo simulations of the SAA observations suggest the existence of a previously unobserved multiply-charged component of the inner radiation belt.

Osborne, W. Z.↗
NASA NTRS1997

Monte Carlo Simulation of THz Multipliers

Schottky Barrier diode frequency multipliers are critical components in submillimeter and Thz space based earth observation systems. As the operating frequency of these multipliers has increased, the agreement between design predictions and experimental results has become poorer. The multiplier design is usually based on a nonlinear model using a form of harmonic balance and a model for the Schottky barrier diode. Conventional voltage dependent lumped element models do a poor job of predicting THz frequency performance. This paper will describe a large signal Monte Carlo simulation of Schottky barrier multipliers. The simulation is a time dependent particle field Monte Carlo simulation with ohmic and Schottky barrier boundary conditions included that has been combined with a fixed point solution for the nonlinear circuit interaction. The results in the paper will point out some important time constants in varactor operation and will describe the effects of current saturation and nonlinear resistances on multiplier operation.

East, J.↗
NASA NTRS2003

Direct Simulation Monte Carlo Calculations in Support of the Columbia Shuttle Orbiter Accident Investigation

The Direct Simulation Monte Carlo method was used to provide 3-D simulations of the early entry phase of the Shuttle Orbiter. Undamaged and damaged scenarios were modeled to provide calibration points for engineering "bridging function" type of analysis. Currently the simulation technology (software and hardware) are mature enough to allow realistic simulations of three dimensional vehicles.

Gallis, Michael A.↗
NASA NTRS1985

Shuttle vertical fin flowfield by the direct simulation Monte Carlo method

The flow properties in a model flowfield, simulating the shuttle vertical fin, determined using the Direct Simulation Monte Carlo method. The case analyzed corresponds to an orbit height of 225 km with the freestream velocity vector orthogonal to the fin surface. Contour plots of the flowfield distributions of density, temperature, velocity and flow angle are presented. The results also include mean molecular collision frequency (which reaches 1/60 sec near the surface), collision frequency density (approaches 7 x 10 to the 18/cu m sec at the surface) and the mean free path (19 m at the surface).

Hueser, J. E.↗
NASA NTRS1985

Monte Carlo simulation of EAS generated by 10(14) - 10(16) eV protons

Detailed Monte Carlo simulations of extensive air showers to be detected by the Homestake Surface Underground Telescope and other similar detectors located at sea level and mountain altitudes have been performed for 10 to the 14th power to 10 to the 16th power eV primary energies. The results of these Monte Carlo calculations will provide an opportunity to compare the experimental data with different models for the composition and spectra of primaries and for the development of air showers. The results obtained for extensive air showers generated by 10 to the 14th power to 10 to the 16th power eV primary protons are reported.

Fenyves, E. J.↗
NASA NTRS1993

Spacecraft contamination investigation by direct simulation Monte Carlo - Contamination on UARS/HALOE

A three-dimensional version of the direct simulation Monte Carlo method is adapted to study the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite, with the ultimate goal of estimating the cumulative level of molecular deposits on critical optical surfaces of the Halogen Occultation Experiment. Emphasis is placed on describing the development of preprocessing codes that simplify the incorporation of geometry and outgassing data into the main run code. A novel use of parallel processing software allows use of multiple processors to solve different regions of the overall flowfield. Initial results relating to contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the Halogen Occultation telescope's aperture plane are presented, along with data related to code performance.

Rault, Didier F. G.↗
NASA NTRS1979

Monte Carlo simulation of the classical two-dimensional one-component plasma

A Monte Carlo simulation, lattice dynamics in the harmonic approximation, and solution of the hypernetted chain equation are used to study the classical two-dimensional one-component plasma. A fluid phase is found for a specified dimensionless quantity, and a solid phase for higher values of this quantity. It is reported that the solid phase shows directional long range order, noting that in the solid phase positional long range order is lost as the thermodynamic limit is approached. Finally, the results of calculations of the thermodynamic functions and one- and two-particle correlation functions are presented.

Gann, R. C.↗
DOE OSTI2025

Python Library for Monte Carlo Simulations with Ab Initio and Machine-Learned Interatomic Potentials

There is a growing need in the simulation community for software that provides a transparent, reproducible, usable, and extensible (TRUE) Monte Carlo (MC) simulation framework employing energies from ab initio methods and machine-learning interatomic potentials (MLIPs). We introduce a Python library (ASE-MC) that adds Monte Carlo functionality to the Atomic Simulation Environment (ASE) package. Now, we can combine the powerful tools used to build systems and perform ab initio and MLIP in ASE with MC simulation algorithms to sample the configurational space with a concise Python script. After presenting the design philosophy, we demonstrate the flexibility of our approach using selected examples. These example simulations include liquid water described with a message-passing MLIP in the canonical and isothermal–isobaric ensembles, sampling the characteristic dihedral angle of biphenyl and comparing an MLIP to first-principles calculations, and a grand canonical Monte Carlo simulation of ammonia adsorption on Pt(111). These examples showcase the main features of the software, which include flexibility in the choice of ab initio or MLIP engine, ab initio or MLIP grand canonical MC with cavity bias insertions and deletions, the ability to add custom MC moves to the move set, and how users can condense complex MC workflows into a single Python script. Finally, this library serves as a framework for reproducible Monte Carlo simulations, facilitating easy reproduction of the work and application to new systems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2021

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗
NASA NTRS1995

Parallel Monte Carlo Simulation for control system design

The research during the 1993/94 academic year addressed the design of parallel algorithms for stochastic robustness synthesis (SRS). SRS uses Monte Carlo simulation to compute probabilities of system instability and other design-metric violations. The probabilities form a cost function which is used by a genetic algorithm (GA). The GA searches for the stochastic optimal controller. The existing sequential algorithm was analyzed and modified to execute in a distributed environment. For this, parallel approaches to Monte Carlo simulation and genetic algorithms were investigated. Initial empirical results are available for the KSR1.

Schubert, Wolfgang M.↗
NASA NTRS1983

Monte Carlo simulation of the effect of soil moisture variation on the microwave emission from soils

In this paper, results of a Monte Carlo simulation of the effect of noise on the relationship between the microwave emissivity of soil and its moisture content are presented. It is found that whenever the magnitude of the noise for the independent variable, in this case the soil moisture, is increased, both the slope of the regression and the correlation coefficient decrease. In particular, when the noise has magnitude equivalent to a coefficient of variation of 0.25, the slope and correlation coefficient are in good agreement with those obtained from the data of a 21-cm airborne microwave radiometer which was flown over a test site in hand county, South Dakota. The comparison was made using a linear relationship to determine the estimated emissivity from the ground measurements of soil moisture. The linear relationship was derived from a radiative transfer model calculation of the microwave emissivities using realistic soil-moisture profiles. The effect of surface roughness was included in the relationship, and the variability of the surface roughness was also simulated by a Monte Carlo technique.

Mo, T.↗
NASA NTRS1993

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↗
NASA NTRS1995

Monte Carlo simulations of precise timekeeping in the Milstar communication satellite system

The Milstar communications satellite system will provide secure antijam communication capabilities for DOD operations into the next century. In order to accomplish this task, the Milstar system will employ precise timekeeping on its satellites and at its ground control stations. The constellation will consist of four satellites in geosynchronous orbit, each carrying a set of four rubidium (Rb) atomic clocks. Several times a day, during normal operation, the Mission Control Element (MCE) will collect timing information from the constellation, and after several days use this information to update the time and frequency of the satellite clocks. The MCE will maintain precise time with a cesium (Cs) atomic clock, synchronized to UTC(USNO) via a GPS receiver. We have developed a Monte Carlo simulation of Milstar's space segment timekeeping. The simulation includes the effects of: uplink/downlink time transfer noise; satellite crosslink time transfer noise; satellite diurnal temperature variations; satellite and ground station atomic clock noise; and also quantization limits regarding satellite time and frequency corrections. The Monte Carlo simulation capability has proven to be an invaluable tool in assessing the performance characteristics of various timekeeping algorithms proposed for Milstar, and also in highlighting the timekeeping capabilities of the system. Here, we provide a brief overview of the basic Milstar timekeeping architecture as it is presently envisioned. We then describe the Monte Carlo simulation of space segment timekeeping, and provide examples of the simulation's efficacy in resolving timekeeping issues.

Camparo, James C.↗
NASA NTRS1989

Direct simulation Monte Carlo analysis on parallel processors

A method is presented for executing a direct simulation Monte Carlo (DSMC) analysis using parallel processing. The method is based on using domain decomposition to distribute the work load among multiple processors, and the DSMC analysis is performed completely in parallel. Message passing is used to transfer molecules between processors and to provide the synchronization necessary for the correct physical simulation. Benchmark problems are described for testing the method and results are presented which demonstrate the performance on two commercially available multicomputers. The results show that reasonable parallel speedup and efficiency can be obtained if the problem is properly sized to the number of processors. It is projected that with a massively parallel system, performance exceeding that of current supercomputers is possible.

Wilmoth, Richard G.↗
NASA NTRS1977

Monte Carlo simulation of the jet stream process

A Monte Carlo model is formulated to simulate the orbital evolution of a system of colliding particles. It is found that inelastic collision alone (even if the impact energy dissipation from collision is very large) does not lead to the formation of a narrow ring-like jet stream; instead, a flat disk structure, similar to Saturn's rings, usually results. To produce the radial focusing effect, it is argued that additional dynamical effects, which would strengthen the collisional interaction between the particles in near-circular orbits, is needed.

Ip, W.-H.↗
NASA NTRS1977

Crater saturation and equilibrium - A Monte Carlo simulation

Conditions of crater saturation and equilibrium have been simulated with a Monte Carlo computer program which utilizes a generating function of the form N = b x D to the -a. The following conclusions have been deduced from the results of the model: (1) saturation and equilibrium for large craters occur at crater densities substantially greater than those observed on the moon, Mars, and Mercury, an indication that the observed populations are essentially production populations, (2) generating functions with a greater than 2 lead to final crater populations with a = 2, (3) generating functions with a between 2 and 1.5 lead to crater populations with the parameter a approaching 1.0, and (4) changes in the slope of crater populations occur gradually and over a wide interval of crater densities and diameters as saturation is approached.

Woronow, A.↗