CAVEAT - A revised version of the general purpose Monte Carlo program, COHORT. Volume 1 - Theory and techniques
Revised version of general purpose Monte Carlo computer code for radiation transport calculations - theory and techniques
Byrn, N. R.↗
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.
A Monte Carlo technique for determining sphere drag in near free molecular flow
First collision Monte Carlo method for calculating aerodynamic drag on hard sphere and Maxwell molecular models in near-free molecular flow
Robertson, S. J.↗
Shield weight optimization using Monte Carlo transport calculations
Outlines are given of the theory used in FASTER-3 Monte Carlo computer program for the transport of neutrons and gamma rays in complex geometries. The code has the additional capability of calculating the minimum weight layered unit shield configuration which will meet a specified dose rate constraint. It includes the treatment of geometric regions bounded by quadratic and quardric surfaces with multiple radiation sources which have a specified space, angle, and energy dependence. The program calculates, using importance sampling, the resulting number and energy fluxes at specified point, surface, and volume detectors. Results are presented for sample problems involving primary neutron and both primary and secondary photon transport in a spherical reactor shield configuration. These results include the optimization of the shield configuration.
Jordan, T. M.↗
A Monte-Carlo simulation of galactic cosmic ray effects in the lunar regolith
A Monte-Carlo soil mixing computer code was used to evaluate isotopic effects induced by galactic cosmic rays in samples on the top surface of the lunar regolith and in lunar core tubes. For surface samples, good agreement was found between theoretical predictions of Ne-21 concentrations, neutron fluences, and Xe-131/Xe-126 ratios, and the corresponding observations.
Langevin, Y.↗
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.↗
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.↗
Monte Carlo simulations of atmospheric spreading functions for space-borne optical sensors
A Monte Carlo radiative transfer model is used to simulate the atmospheric spreading effects. The spreading functions for several vertical aerosol profiles are obtained. The dependence of atmospheric conditions and aerosol properties are investigated, and the importance of the effect on MSS and TM measurements are assessed.
Kiang, R. K.↗
A Monte Carlo simulation of the Infrared Astronomical Satellite (IRAS) mission
The development of a Monte Carlo simulation of numerous aspects of the IRAS mission is described. The formulations used in modelling the infrared sky, the radiation environment, and the telescope are discussed. The results of running the simulation with varying mission parameters are presented.
Wolf, D. M.↗
Monte Carlo spectra from pair-equilibrium, weakly magnetized thermal plasmas
The results of Monte Carlo simulations of weakly magnetized, relativistic thermal plasmas are reported, including such physical processes as thermal synchrotron emission and absorption, bremsstrahlung, and Compton scattering, as well as pair-annihilation and pair-production. The rates of the latter two processes are kept approximately equal to each other, i.e., the plasmas are in pair equilibrium. The resulting spectra are of power-law form, and can be well described by a set of simple analytical formulae. Neither a pair annihilation feature nor a Wien hump is seen in the spectra. The possible relevance of the results to active galactic nuclei and gamma-bursts is discussed.
Zdziarski, A. A.↗
Adaptive domain decomposition for Monte Carlo simulations on parallel processors
A method is described for performing direct simulation Monte Carlo (DSMC) calculations on parallel processors using adaptive domain decomposition to distribute the computational work load. The method has been implemented on a commercially available hypercube and benchmark results are presented which show the performance of the method relative to current supercomputers. The problems studied were simulations of equilibrium conditions in a closed, stationary box, a two-dimensional vortex flow, and the hypersonic, rarefield flow in a two-dimensional channel. For these problems, the parallel DSMC method ran 5 to 13 times faster than on a single processor of a Cray-2. The adaptive decomposition method worked well in uniformly distributing the computational work over an arbitrary number of processors and reduced the average computational time by over a factor of two in certain cases.
Wilmoth, Richard G.↗
Adaptive domain decomposition for Monte Carlo simulations on parallel processors
A method is described for performing direct simulation Monte Carlo (DSMC) calculations on parallel processors using adaptive domain decomposition to distribute the computational work load. The method has been implemented on a commercially available hypercube and benchmark results are presented which show the performance of the method relative to current supercomputers. The problems studied were simulations of equilibrium conditions in a closed, stationary box, a two-dimensional vortex flow, and the hypersonic, rarefied flow in a two-dimensional channel. For these problems, the parallel DSMC method ran 5 to 13 times faster than on a single processor of a Cray-2. The adaptive decomposition method worked well in uniformly distributing the computational work over an arbitrary number of processors and reduced the average computational time by over a factor of two in certain cases.
Wilmoth, Richard G.↗
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.↗
A Monte Carlo approach to the evaluation of conditional expectation parameter estimates for nonlinear dynamic systems.
Monte Carlo methods for conditional expectation for nonlinear dynamic systems with acceptable computing time expenditure
Mcghee, R. B.↗
Monte Carlo treatment of Lyman-alpha radiation in a plane-parallel atmosphere.
A Monte Carlo technique involving Stokes vectors is used to obtain the state of polarization and intensity of solar Lyman-alpha photons as they diffuse through a plane-parallel homogeneous model of earth's hydrogen envelope. Fine structure of Lyman-alpha and Doppler redistribution of frequencies are taken into account. Comparison of the results with Heath's observed upper limit for polarization of 1.5 per cent implies an optical thickness tau greater than 7 and intensities of 8-10 kilorayleighs for a solar Lyman-alpha flux of 5.8 ergs per sq cm per sec.
Modali, S. B.↗
Monte Carlo approach to the spatial deposition of energy by electrons in molecular hydrogen
The Monte Carlo (MC) and continuous slowdown approximation (CSDA) approaches to the spatial deposition of energy by electrons are compared using the same detailed atomic cross section (DACS). It is found that the CSDA method overestimates the amount of energy that is deposited near the end of the path for electrons above a few hundred electron volts. The MC results are in approximate agreement with experimental data in such a way as to be relatively independent of the actual gas used. Our MC results are extended to obtain the three-dimensional deposition of energy by sub-keV electrons in molecular hydrogen.
Heaps, M. G.↗
Monte Carlo calculations of diatomic molecule gas flows including rotational mode excitation
The direct simulation Monte Carlo method was used to solve the Boltzmann equation for flows of an internally excited nonequilibrium gas, namely, of rotationally excited homonuclear diatomic nitrogen. The semi-classical transition probability model of Itikawa was investigated for its ability to simulate flow fields far from equilibrium. The behavior of diatomic nitrogen was examined for several different nonequilibrium initial states that are subjected to uniform mean flow without boundary interactions. A sample of 1000 model molecules was observed as the gas relaxed to a steady state starting from three specified initial states. The initial states considered are: (1) complete equilibrium, (2) nonequilibrium, equipartition (all rotational energy states are assigned the mean energy level obtained at equilibrium with a Boltzmann distribution at the translational temperature), and (3) nonequipartition (the mean rotational energy is different from the equilibrium mean value with respect to the translational energy states). In all cases investigated the present model satisfactorily simulated the principal features of the relaxation effects in nonequilibrium flow of diatomic molecules.
Yoshikawa, K. K.↗
STS-1 operational flight profile. Volume 5: Descent cycle 3. Appendix D: GRTLS six degree of freedom Monte Carlo dispersion analysis
The results of a six degree of freedom (6-DOF) nonlinear Monte Carlo dispersion analysis for the latest glide return to landing site (GRTLS) abort trajectory for the Space Transportation System 1 Flight are presented. For this GRTLS, the number two main engine fails at 262.5 seconds ground elapsed time. Fifty randomly selected simulations, initialized at external tank separation, are analyzed. The initial covariance matrix is a 20 x 20 matrix and includes navigation errors and dispersions in position and velocity, time, accelerometer bias, and inertial platform misalinements. In all 50 samples, speedbrake, rudder, elevon, and body flap hinge moments are acceptable. Transitions to autoland begin before 9,000 feet and there are no tailscrapes. Navigation derived dynamic pressure accuracies exceed the flight control system constraints above Mach 2.5. Three out of 50 landings exceeded tire specification limit speed of 222 knots. Pilot manual landings are expected to reduce landing speed by landing farther downrange.
Montez, M. N.↗
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.↗