Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Techniques”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Collisional evolution of the mass-distribution spectrum of planetesimals. II

An investigation of conditions for the early growth of planetesimals through their mutual collisions using Monte Carlo technique is presented. Rebound, erosion, catastrophic breakup, and coagulation are considered, mean random velocity is assumed to vary with time, and growth of km-sized planetesimals is shown to occur only when they have the mechanical properties similar to iron meteorites. The formation time of several hundred-km sized planetesimals is much longer than previously reported by Greenberg (1978), and an expression was derived for the temporal variation of the mean mass.

Matsui, T.↗
NASA NTRS1980

Spatial aspects of low- and medium-energy electron degradation in N2

Spatial (radial and longitudinal) yield spectra for electron energy degradation in molecular nitrogen gas for 25-eV to 10-keV incident electrons have been generated by using a Monte Carlo technique. These spatial yield spectra associated with the electron degradation process can be employed to calculate a 'yield' for any inelastic state at any position in the medium. These have been analytically represented in terms of a model containing three simple 'microplumes'. Five-dimensional yield spectra which contain the information about the polar angle of the electron have also been analytically represented within the framework of the microplume model. Aeronomical and radiological applications of our model are discussed.

Singhal, R. F.↗
NASA NTRS1980

Absorption coefficient instrument for turbid natural waters

The paper presents an instrument that directly measures multispectral absorption coefficient of turbid natural water. Attention is given to the design, which is shown to incorporate methods for the compensation of variation in the internal light source intensity, correction of the spectrally dependent nature of the optical elements, and correction for variation in the background light level. In addition, when used in conjunction with a spectrally matched total attenuation instrument, the spectrally dependent scattering coefficient can also be derived. Finally, it is reported that systematic errors associated with multiple scattering have been estimated using Monte Carlo techniques.

Friedman, E.↗
NASA NTRS1980

Neutral cometary atmospheres. I - An average random walk model for photodissociation in comets

A method for constructing physically realistic photochemical models, taking into account the isotropic ejection of dissociated molecular fragments as well as radiation pressure acceleration, has been developed using Monte Carlo techniques. The effect of the isotropic ejection, as opposed to the radial motion arbitrarily assumed in Haser's model, is adequately described by a simple average random walk model. It is shown that measured radial (Haser) scale lengths are in fact only lower limits to a range of possible true scale lengths for a given brightness profile, which explains the current discrepancies between observed scale lengths and those predicted by photochemistry.

Combi, M. R.↗
NASA NTRS1980

Two-dimensional O2 adsorbed on graphite

Properties of two-dimensional O2 adsorbed on graphite are calculated in the extremely low-coverage delta region and for monolayers, with use of pattern-recognition optimization and Monte Carlo techniques. Equilibrium configurations and orientations, orientational order-disorder, melting, and dissociation transitions are predicted at various O2 densities. Phase characteristics, including a plastic crystallite phase, are compared with experiment.

Etters, R. D.↗
NASA NTRS1980

Simulation studies of a model of high-density metallic hydrogen

Upper bounds for the ground-state energies of liquid and solid phases of metallic hydrogen and metallic deuterium have been calculated with variational methods and Monte Carlo techniques. At four densities (0.8, 1.2, 1.36, and 1.488) crystalline phases are clearly preferred in the sense that the energy difference, when compared to the liquid, is in excess of the errors inherent in the numerical procedures. At a fifth density (1.6), the energy differences between solid and liquid phases are smaller than these errors.

Mon, K. K.↗
NASA NTRS1981

Inelastic scattering in planetary atmospheres. I - The Ring effect, without aerosols

The contribution of inelastic molecular scattering (Rayleigh-Brillouin and rotational Raman scattering) to the filling-in of Fraunhofer lines in the light of the blue sky is studied. Aerosol fluorescence is shown to be negligible, and aerosol scattering is ignored. The angular and polarization dependences of the filling-in detail for single scattering are discussed. An approximate treatment of multiple scattering, using a backward Monte Carlo technique, makes it possible to investigate the effects of the ground albedo. As the molecular scatterings alone produce more line-filling than is observed, it seems likely that aerosols dilute the effect by contributing unaltered sunlight to the observed spectra.

Kattawar, G. W.↗
NASA NTRS1981

A model for X-ray pulsar spectra

The transfer of X-rays through a hot, magnetized, scattering atmosphere is studied using a Monte Carlo technique. The process considered is Compton scattering, including cyclotron resonance scattering. Unlike previous work on the subject, the effects of frequency shifts, relativistic corrections, recoil, and polarization mode changes have been included. It is found that the emergent X-ray spectra harden toward the magnetic equator, are cut off at a frequency less than the fundamental cyclotron frequency, and depend on the magnetic field strength. The relevance of the results to observed pulsar spectra is indicated.

Pravdo, S. H.↗
NASA NTRS1981

The effect of the charge exchange source on the velocity and 'temperature' distributions and their anisotropies in the earth's exosphere

The velocity distribution of atomic hydrogen in the earth's exosphere is calculated as a function of altitude and direction taking into account both the classic exobase source and the higher-altitude plasmaspheric charge exchange source. Calculations are performed on the basis of a Monte Carlo technique in which random ballistic trajectories of individual atoms are traced through a three-dimensional grid of audit zones, at which relative concentrations and momentum or energy fluxes are obtained. In the case of the classical exobase source alone, the slope of the velocity distribution is constant only for the upward radial velocity component and increases dramatically with altitude for the incoming radial and transverse velocity components, resulting in a temperature decrease. The charge exchange source, which produces the satellite hydrogen component and the hot ballistic and escape components of the exosphere, is found to enhance the wings of the velocity distributions, however this effect is not sufficient to overcome the temperature decreases at altitudes above one earth radius. The resulting global model of the hydrogen exosphere may be used as a realistic basis for radiative transfer calculations.

Hodges, R. R., Jr.↗
NASA NTRS1981

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

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

Annihilation radiation from a hot e/+/-e/-/ plasma

The radiation from electron-positron annihilations in a plasma of temperature above 10 to the 8th K is investigated as a possible source of the emission line at energies between 400 and 460 keV frequently seen in gamma-ray bursts. The annihilation rate and luminosity of an optically thin electron-positron plasma and the energy distribution of the resulting annihilation radiation are calculated by the use of a Monte Carlo technique as functions of temperature. Results indicate the annihilation spectrum to be peaked at an energy of 0.511 MeV plus a temperature-dependent blueshift, and the annihilation line to be significantly temperature-broadened. The widths of the observed burst emission lines set an upper limit of 3 x 10 to the 8th K on the temperature of any pair annihilation region in burst sources, which is considerably lower than the typical kinetic temperatures of the radiating particles. It is thus inferred that either the annihilation region is nonthermal, or spatially distinct from the burst site.

Ramaty, R.↗
NASA NTRS1981

Monte Carlo simulation of collisionless shocks showing preferential acceleration of high A/Z particles

A collisionless quasi-parallel shock is simulated by Monte Carlo techniques. The scattering of all velocity particles from thermal to high energy is assumed to occur so that the mean free path is directly proportional to velocity times the mass-to-charge-ratio, and inversely proporational to the plasma density. The shock profile and velocity spectra are obtained, showing preferential acceleration of high A/Z particles relative to protons. The inclusion of the back pressure of the scattering particles on the inflowing plasma produces a smoothing of the shock profile, which implies that the spectra are steeper than for a discontinuous shock.

Ellison, D. C.↗
NASA NTRS1980

Cosmic ray production curves below reworking zones

A method is presented for calculating cosmic ray production profiles below reworking zones. The method uses an input reworking depth determined from data such as signatures in the depth profile of ferromagnetic resonance intensity and input cosmic ray production profiles for an undisturbed surface. Reworking histories are simulated using Monte Carlo techniques, and depth profiles are used to determine cosmic ray exposure age limits with a specified probability. It is shown that the track density profiles predict cosmic ray exposure ages in lunar cores that are consistent with values determined by other methods. Results applied to neutron fluence and spallation rare gases eliminate the use of reworking depth as an adjustable parameter and give cosmic ray exposure ages that are compatible with each other.

Blanford, G. E.↗
NASA NTRS1982

Statistics of Delta v magnitude for a trajectory correction maneuver containing deterministic and random components

A number of interplanetary missions now being planned involve placing deterministic maneuvers along the flight path to alter the trajectory. Lee and Boain (1973) examined the statistics of trajectory correction maneuver (TCM) magnitude with no deterministic ('bias') component. The Delta v vector magnitude statistics were generated for several values of random Delta v standard deviations using expansions in terms of infinite hypergeometric series. The present investigation uses a different technique (Monte Carlo simulation) to generate Delta v magnitude statistics for a wider selection of random Delta v standard deviations and also extends the analysis to the case of nonzero deterministic Delta v's. These Delta v magnitude statistics are plotted parametrically. The plots are useful in assisting the analyst in quickly answering questions about the statistics of Delta v magnitude for single TCM's consisting of both a deterministic and a random component. The plots provide quick insight into the nature of the Delta v magnitude distribution for the TCM.

Bollman, W. E.↗
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 NTRS1983

Dynamical evolution of the Oort cometary cloud

The dynamical evolution of comets in the Oort cloud under the influence of stellar perturbations has been modeled using Monte Carlo techniques. It is shown that the cloud has been depleted over the history of the solar system. Comets are lost from the cloud by direct ejection due to close stellar encounters, diffusion of aphelia to distances beyond the sun's sphere of influence, or diffusion of perihelia into the planetary region where Jupiter and Saturn perturbations either eject them on hyperbolic trajectories or capture them to short-period orbits. The population of the cloud is estimated to be 1.0 - 1.5 x 10 to the 12th comets and the total mass is on the order of 1.9 earth masses. In addition to random passing stars, less frequent encounters with giant molecular clouds may play a significant role in randomizing the orbits of comets in the cloud and reducing the effective radius of the sun's sphere of influence.

Weissman, P. R.↗
NASA NTRS1984

Computing approximate random Delta v magnitude probability densities

This paper describes the development and use of an algorithm to compute approximate statistics of the magnitude of a single random trajectory correction maneuver (TCM) Delta v vector. The TCM Delta v vector is modeled as a three component Cartesian vector each of whose components is a random variable having a normal (Gaussian) distribution with zero mean and possibly unequal standard deviations. The algorithm uses these standard deviations as input to produce approximations to (1) the mean and standard deviation of the magnitude of Delta v, (2) points of the probability density function of the magnitude of Delta v, and (3) points of the cumulative and inverse cumulative distribution functions of Delta v. The approximates are based on Monte Carlo techniques developed in a previous paper by the author and extended here. The algorithm described is expected to be useful in both pre-flight planning and in-flight analysis of maneuver propellant requirements for space missions.

Chadwick, C.↗