Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte”

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 145 records · Page 8

Uncertainty propagation in a stratospheric model. I - Development of a concise stratospheric model. II - Monte Carlo analysis of imprecisions due to reaction rates

The paper discusses the development of a concise stratospheric model which uses iteration to obtain coupling between interacting species. The one-dimensional, steady-state, diurnally-averaged model generates diffusion equations with appropriate sources and sinks for species odd oxygen, H2O, H2, CO, N2O, odd nitrogen, CH4, CH3Cl, CCl4, CF2Cl2, CFCl3, and odd chlorine. The model evaluates steady-state perturbations caused by injections of chlorine and NO(x) and may be used to predict ozone depletion. The model is used in a Monte Carlo study of the propagation of reaction-rate imprecisions by calculating an ozone perturbation caused by the addition of chlorine. Since the model is sensitive to only 10 of the more than 50 reaction rates considered, only about 1000 Monte Carlo cases are required to span the space of possible results.

Rundel, R. D.↗

Monte Carlo modeling of exospheric bodies - Mercury

In order to study the interaction with the surface, a Monte Carlo program is developed to determine the distribution with altitude as well as the global distribution of density at the surface in a single operation. The analysis presented shows that the appropriate source distribution should be Maxwell-Boltzmann flux if the particles in the distribution are to be treated as components of flux. Monte Carlo calculations with a Maxwell-Boltzmann flux source are compared with Mariner 10 UV spectrometer data. Results indicate that the presently operating models are not capable of fitting the observed Mercury exosphere. It is suggested that an atmosphere calculated with a barometric source distribution is suitable for more realistic future exospheric models.

Smith, G. R.↗

A comparative study of Conroy and Monte Carlo methods applied to multiple quadratures and multiple scattering

An efficient numerical method of multiple quadratures, the Conroy method, is applied to the problem of computing multiple scattering contributions in the radiative transfer through realistic planetary atmospheres. A brief error analysis of the method is given and comparisons are drawn with the more familiar Monte Carlo method. Both methods are stochastic problem-solving models of a physical or mathematical process and utilize the sampling scheme for points distributed over a definite region. In the Monte Carlo scheme the sample points are distributed randomly over the integration region. In the Conroy method, the sample points are distributed systematically, such that the point distribution forms a unique, closed, symmetrical pattern which effectively fills the region of the multidimensional integration. The methods are illustrated by two simple examples: one, of multidimensional integration involving two independent variables, and the other, of computing the second order scattering contribution to the sky radiance.

Deepak, A.↗

Monte Carlo analysis of satellite debris footprint dispersion

A comprehensive study is performed to investigate satellite debris impact point dispersion using a combination of Monte Carlo statistical analysis and parametric methods. The Monte Carlo technique accounts for nonlinearities in the entry point dispersion, which is represented by a covariance matrix of position and velocity errors. Because downrange distance of impact is a monotonic function of debris ballistic coefficient, a parametric method is useful for determining dispersion boundaries. The scheme is applied in the present analysis to estimate the Skylab footprint dispersions for a controlled reentry. A significant increase in the footprint dispersion is noticed for satellite breakup above a 200,000-ft altitude. A general discussion of the method used for analysis is presented together with some typical results obtained for the Skylab deboost mission, which was designed before NASA abandoned plans for a Skylab controlled reentry.

Rao, P. P.↗

Direct simulation Monte Carlo technique for modeling of the environment in the vicinity of the Space Shuttle Orbiter

A program for analyzing the flowfield parameters in the neighborhood of the Space Shuttle Orbiter has been developed. The program uses the direct simulation Monte Carlo method, which is a completely probabilistic Monte Carlo technique capable of analyzing 3-dimensional steady or unsteady flow with prescribed internal and external boundary conditions. The freestream flux densities incident on the external flowfield boundaries are calculated from the drifting Maxwellian gas properties of the freestream. The flux entering the flowfield from the internal boundary is calculated from the outgassing flux density distribution over the Shuttle external surface and the prescribed discrete source fluxes. This technique produces a numerical flowfield solution which is the probabilistic equivalent of a complete solution of the time-dependent, 3-dimensional Boltzmann equation. Flowfield results are presented for the following configurations: (1) Shuttle angles of attack of 0 and 90 deg, (2) freestream density 10 to the 9th - 10 to the 11th per cu cm, (3) Shuttle outgassing flux density, and (4) operation of the aft downfiring vernier Reaction Control System (RCS) engine. Results are presented for column density distribution of outgassed and engine species and for the flux density of outgassed and engine species incident on the Shuttle bay and the windshield.

Melfi, L. T., Jr.↗

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

Monte Carlo calibration of the SMM gamma ray spectrometer for high energy gamma rays and neutrons

The Gamma Ray Spectrometer (GRS) on the Solar Maximum Mission spacecraft was primarily designed and calibrated for nuclear gamma ray line measurements, but also has a high energy mode which allows the detection of gamma rays at energies above 10 MeV and solar neutrons above 20 MeV. The GRS response has been extrapolated until now for high energy gamma rays from an early design study employing Monte Carlo calculations. The response to 50 to 600 MeV solar neutrons was estimated from a simple model which did not consider secondary charged particles escaping into the veto shields. In view of numerous detections by the GRS of solar flares emitting high energy gamma rays, including at least two emitting directly detectable neutrons, the calibration of the high energy mode in the flight model has been recalculated by the use of more sophisticated Monte Carlo computer codes. New results presented show that the GRS response to gamma rays above 20 MeV and to neutrons above 100 MeV is significantly lower than the earlier estimates.

Cooper, J. F.↗

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

A Monte Carlo reflectance model for soil surfaces with three-dimensional structure

A Monte Carlo soil reflectance model has been developed to study the effect of macroscopic surface irregularities larger than the wavelength of incident flux. The model treats incoherent multiple scattering from Lambertian facets distributed on a periodic surface. Resulting bidirectional reflectance distribution functions are non-Lambertian and compare well with experimental trends reported in the literature. Examples showing the coupling of the Monte Carlo soil model to an adding bidirectional canopy of reflectance model are also given.

Cooper, K. D.↗

Large-cell Monte Carlo renormalization of irreversible growth processes

Monte Carlo sampling is applied to a recently formulated direct-cell renormalization method for irreversible, disorderly growth processes. Large-cell Monte Carlo renormalization is carried out for various nonequilibrium problems based on the formulation dealing with relative probabilities. Specifically, the method is demonstrated by application to the 'true' self-avoiding walk and the Eden model of growing animals for d = 2, 3, and 4 and to the invasion percolation problem for d = 2 and 3. The results are asymptotically in agreement with expectations; however, unexpected complications arise, suggesting the possibility of crossovers, and in any case, demonstrating the danger of using small cells alone, because of the very slow convergence as the cell size b is extrapolated to infinity. The difficulty of applying the present method to the diffusion-limited-aggregation model, is commented on.

Nakanishi, H.↗

Monte Carlo study of a model of diffusion-controlled reactions

The Monte Carlo method is used to perform averages over sink configurations in the present study of diffusion-controlled reactions occurring between solute particles and immobile spherical sinks. In order to determine the average steady state solute concentration profile in a locally perturbed solution for sink volume fractions phi of less than 0.3, the diffusion equation in the monopolar plus dipolar approximation of diffusive couplings between the sinks is solved numerically. The Monte Carlo method is shown to be the most accurate and efficient in the phi = 0.001-0.1 region, where a system of only 25 sinks suffices and the monopolar approximation alone is sufficiently accurate.

Beenakker, C. W. J.↗

A microstructural lattice model for strain oriented problems: A combined Monte Carlo finite element technique

A specialized, microstructural lattice model, termed MCFET for combined Monte Carlo Finite Element Technique, was developed which simulates microstructural evolution in material systems where modulated phases occur and the directionality of the modulation is influenced by internal and external stresses. In this approach, the microstructure is discretized onto a fine lattice. Each element in the lattice is labelled in accordance with its microstructural identity. Diffusion of material at elevated temperatures is simulated by allowing exchanges of neighboring elements if the exchange lowers the total energy of the system. A Monte Carlo approach is used to select the exchange site while the change in energy associated with stress fields is computed using a finite element technique. The MCFET analysis was validated by comparing this approach with a closed form, analytical method for stress assisted, shape changes of a single particle in an infinite matrix. Sample MCFET analytical for multiparticle problems were also run and in general the resulting microstructural changes associated with the application of an external stress are similar to that observed in Ni-Al-Cr alloys at elevated temperature.

Gayda, J.↗

Monte Carlo simulation of an expanding gas

By application of simple computer graphics techniques, the statistical performance of two Monte Carlo methods used in the simulation of rarefied gas flows are assessed. Specifically, two direct simulation Monte Carlo (DSMC) methods developed by Bird and Nanbu are considered. The graphics techniques are found to be of great benefit in the reduction and interpretation of the large volume of data generated, thus enabling important conclusions to be drawn about the simulation results. Hence, it is discovered that the method of Nanbu suffers from increased statistical fluctuations, thereby prohibiting its use in the solution of practical problems.

Boyd, Iain D.↗

Comparison of dose estimates using the buildup-factor method and a Baryon transport code (BRYNTRN) with Monte Carlo results

Continuing efforts toward validating the buildup factor method and the BRYNTRN code, which use the deterministic approach in solving radiation transport problems and are the candidate engineering tools in space radiation shielding analyses, are presented. A simplified theory of proton buildup factors assuming no neutron coupling is derived to verify a previously chosen form for parameterizing the dose conversion factor that includes the secondary particle buildup effect. Estimates of dose in tissue made by the two deterministic approaches and the Monte Carlo method are intercompared for cases with various thicknesses of shields and various types of proton spectra. The results are found to be in reasonable agreement but with some overestimation by the buildup factor method when the effect of neutron production in the shield is significant. Future improvement to include neutron coupling in the buildup factor theory is suggested to alleviate this shortcoming. Impressive agreement for individual components of doses, such as those from the secondaries and heavy particle recoils, are obtained between BRYNTRN and Monte Carlo results.

Shinn, Judy L.↗

Three-dimensional Monte Carlo simulation of gamma-ray scattering and production in the atmosphere

Results are reported from Monte Carlo numerical simulations of atmospheric gamma-ray scattering and production. The basic physical principles involved in the construction of the models are reviewed, and results are presented in extensive graphs for low-energy gamma rays with the spectra of gamma-ray bursts, solar flares, the Crab pulsar, and 511-keV line radiation. It is shown that the model accurately reproduces the characteristics of atmospheric albedo radiation, including details of the angular distribution. The potential applicability of the Monte Carlo technique to studies of the near-earth radiation environment is indicated.

Morris, Daniel J.↗

Surface entropy of liquids via a direct Monte Carlo approach - Application to liquid Si

Two methods are presented for a direct Monte Carlo evaluation of the surface entropy S(s) of a liquid interacting by specified, volume-independent potentials. The first method is based on an application of the approach of Ferrenberg and Swendsen (1988, 1989) to Monte Carlo simulations at two different temperatures; it gives much more reliable results for S(s) in liquid Si than previous calculations based on numerical differentiation. The second method expresses the surface entropy directly as a canonical average at fixed temperature.

Wang, Z. Q.↗

Comparison of Burnett, super-Burnett and Monte Carlo solutions for hypersonic shock structure

The continuum Navier-Stokes, Burnett, and super-Burnett equations are solved for one-dimensional shock structures in various monoatomic gases. Solutions for a hard sphere gas, argon, and a Maxwellian gas from Mach 1.3 to Mach 50 are obtained. A numerical method utilizing the complete time-dependent continuum equations and obtaining the steady-state shock structure by allowing the system to relax from arbitrary initial conditions is employed. Shock density, velocity, temperature, and entropy profiles are also obtained using the direct simulation Monte Carlo method, and these results are used as bases for comparison for continuum solution profiles. It is shown that the Burnett equations yield shock structure solutions in much closer agreement to both Monte Carlo and experimental results than do the Navier-Stokes equations. Solutions to the super-Burnett equations with coefficients as presently derived, however, are considered to be inferior to those of the Burnett equations.

Fiscko, Kurt A.↗

A study of microkinetic adjustments required to match shock wave experiments and Monte Carlo Direct Simulation for a wide Mach number range

Shock wave thickness predictions from Monte Carlo Direct Simulations, using differential scattering and the Maitland-Smith-Aziz interatomic potential, underpredict experiments as shock Mach numbers increase above about 4. Examination of several sources of data has indicated that at relatively high energies the repulsive portion of accepted potentials such as the Maitland-Smith-Aziz may be too steep. An Exponential-6 potential due to Ross, based on high energy molecular beam scattering data and shock velocity measurements in liquid argon, has been combined with the lower energy portion of the Maitland-Smith-Aziz potential. When this hybrid potential is used in Monte Carlo Direct Simulations, agreement with experiments is improved over the previous predictions using the pure Maitland-Smith-Aziz form.

Pham-Van-diep, Gerald C.↗