Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo code”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Monte Carlo Simulation of Nonlinear Radiation Induced Plasmas

A Monte Carlo simulation model for radiation induced plasmas with nonlinear properties due to recombination was, employing a piecewise linearized predict-correct iterative technique. Several important variance reduction techniques were developed and incorporated into the model, including an antithetic variates technique. This approach is especially efficient for plasma systems with inhomogeneous media, multidimensions, and irregular boundaries. The Monte Carlo code developed has been applied to the determination of the electron energy distribution function and related parameters for a noble gas plasma created by alpha-particle irradiation. The characteristics of the radiation induced plasma involved are given.

Wang, B. S.↗
NASA NTRS1976

Statistical mechanics of light elements at high pressure. IV - A model free energy for the metallic phase

A large quantity of data on the thermodynamic properties of hydrogen-helium metallic liquids have been obtained in extended computer calculations in which a Monte Carlo code essentially identical to that described by Hubbard (1972) was used. A model free energy for metallic hydrogen with a relatively small mass fraction of helium is discussed, taking into account the definition of variables, a procedure for choosing the free energy, values for the fitting parameters, and the evaluation of the entropy constants. Possibilities concerning a use of the obtained data in studies of the interiors of the outer planets are briefly considered.

Dewitt, H. E.↗
NASA NTRS1977

Spectra and angular distributions of atmospheric gamma rays from 0.3 to 10 MeV at lambda = 40 deg

Measurements of the spectral and angular distributions of atmospheric gamma sq cm rays in the energy range 0.3-10 MeV over Palestine, Texas, at residual depths of 2.5 and 70 g/sq cm are reported. In confirmation of the general features of a model prediction, the measurements show at 2.5 g/sq cm upward moving fluxes greater than the downward moving fluxes, the effect increasing with energy, and approximate isotropy at 70 g/sq cm. Numerous characteristic gamma-ray lines were observed, most prominently at 0.511, 1.6, 2.3, 4.4, and 6.1 MeV. Their intensities were also compared with model predictions. Observations were made with an actively shielded scintillator counter with two detectors, one of aperture 50 deg FWHM and the other of 120 deg FWHM. Above 1 MeV, contributions to the counting rate from photons penetrating the shield annulus and from neutron interactions were large; they were studied by means of a Monte Carlo code and are extensively discussed.

Ling, J. C.↗
NASA NTRS1978

Consequences of narrow cyclotron emission from Hercules X-1

The implications of the recent observations of a narrow cyclotron line in the hard X-ray spectrum of Hercules X-1 are studied. A Monte Carlo code is used to simulate the X-ray transfer of an intrinsically narrow feature at approximately 56 keV through an opaque, cold magnetospheric shell. The results of this study indicate that if a narrow line can be emitted by the source region, then only about 10% of the photons remain in a narrow feature after scattering through the shell. The remaining photons are scattered into a broad feature (FWHM approximately 30 keV) that peaks near 20 keV. Thus, these calculations indicate that the intrinsic source luminosity of the cyclotron line is at least an order of magnitude greater than the observed luminosity.

Weaver, R. P.↗
NASA NTRS1975

The simulated depth history of dust grains in the lunar regolith

Trajectories giving the individual depth variations of lunar dust grains with time are randomly generated by a Monte Carlo code, where the variables are the mass and speed distribution of meteorites at the lunar surface and the geometrical shape of impact craters. A statistical analysis of a great number of such trajectories is then used to define the 'average' depth history of lunar dust grains for two grain radii: 1 and 50 microns. This yields: (1) estimates for time constants involved during the dynamic evolution of the regolith, (2) a model for the layering of the regolith, and (3) a better understanding of the basic dust-grain mechanisms responsible for the formation of the most mature lunar soil samples. The validity of various soil models proposed for the dynamic evolution of the regolith is discussed in terms of experimental constraints on the models.

Duraud, J. P.↗
NASA NTRS1983

The evolution of clusters of galaxies. I - Very rich clusters

A multiple one-body Monte Carlo code is used to investigate the evolution of galaxies in a steady cluster potential under the influence of dynamical friction, two-body relaxation, tidal stripping, and galaxy mergers. The basic assumptions, estimated time scales, method and computer program, and effects of uncertainties in galaxy encounter physics are addressed. Numerical experiments in which the mass initially carried by galaxies and the mass function are varied are reported. It is shown that the formation of a very massive cluster galaxy depends critically on the number of galaxies, the initial division of cluster material between galaxies and a smooth intracluster medium, the mass spectrum of the galaxies, and chance. By the end of all simulations, less than half of the mass in the central regions of the cluster is bound to galaxies. It appears possible to produce any of the Bautz-Morgan classes from very similar initial conditions.

Richstone, D. O.↗
NASA NTRS1985

Monte Carlo simulation of a new gamma ray telescope

A new Monte Carlo code has been written to simulate the response of the new University of California double scatter gamma ray telescope. This package of modular software routines, written in VAX FORTRAN 77 simulates the detection of 0.1 to 35 MeV gamma rays. The new telescope is flown from high altitude balloons to measure medium energy gamma radiation from astronomical sources. This paper presents (1) the basic physics methods in the code, and (2) the predicted response functions of the telescope. Gamma ray processes include Compton scattering, pair production and photoelectric absorption in plastic scintillator, NaI(Tl) and aluminum. Electron transport processes include ionization energy loss, multiple scattering, production of bremsstrahlung photons and positron annihilation.

Simone, J.↗
NASA NTRS1987

Deceleration of infalling plasma in the atmospheres of accreting neutron stars. I - Isothermal atmospheres

The stopping of infalling protons in a strong magnetic field is studied in detail using linear response theory and a Monte Carlo code. The calculations include both Coulomb collisions with atmospheric electrons and collective forces due to the resonant excitation of the atmospheric electron plasma modes by the impinging protons. Graphs are presented showing how both energy and momentum are deposited in the atmosphere by the flow. A striking dependence of the deceleration on the temperature is found which distinguishes the magnetic from the nonmagnetic process.

Miller, G. S.↗
NASA NTRS1988

On the transport of ions released in the magnetotail by the AMPTE-IRM satellite

The Ba and Li ions releasd into the magnetotail in spring 1985 by the AMPTE-IRM satellite were not observed subsequently in the inner magnetosphere with the AMPTE-CCE satellite. These results were studied by using a Monte Carlo code to compute the transport of the ions. For each release several hundred ion guiding-center trajectories were computed under simulated magnetospheric conditions, using the Tsyganenko-Usmanov (1982) magnetic-field model and the Millstone Hill convection-electric-field model (Oliver et al., 1983). The corotation and convection electric fields were mapped to altitudes above the ionosphere, assuming the magnetic-field lines to be equipotentials. The initial conditions of the ions, at the times at which the ions were picked up by the electric field, were estimated by taking into consideration the release conditions and the early-time collective effects. The results indicate that the Ba(+) ions were not observed because the CCE satellite was not along the drift paths of the ions, and the Li(+) ions were not observed because their fluxes at the satellite were too low.

Cladis, J. B.↗
NASA NTRS1990

The effect of catastrophic collisional fragmentation and diffuse medium accretion on a computational interstellar dust system

The effects of catastrophic collisional fragmentation and diffuse medium accretion on a the interstellar dust system are computed using a Monte Carlo computer model. The Monte Carlo code has as its basis an analytic solution of the bulk chemical evolution of a two-phase interstellar medium, described by Liffman and Clayton (1989). The model is subjected to numerous different interstellar processes as it transfers from one interstellar phase to another. Collisional fragmentation was found to be the dominant physical process that shapes the size spectrum of interstellar dust. It was found that, in the diffuse cloud phase, 90 percent of the refractory material is locked up in the dust grains, primarily due to accretion in the molecular medium. This result is consistent with the observed depletions of silicon. Depletions were found to be affected only slightly by diffuse cloud accretion.

Liffman, Kurt↗
NASA NTRS1991

The estimation of background production by cosmic rays in high-energy gamma ray telescopes

A calculational method of estimating instrumental background in high-energy gamma-ray telescopes, using the hadronic Monte Carlo code FLUKA87, is presented. The method is applied to the SAS-2 and EGRET telescope designs and is also used to explore the level of background to be expected for alternative configurations of the proposed GRITS telescope, which adapts the external fuel tank of a Space Shuttle as a gamma-ray telescope with a very large collecting area. The background produced in proton-beam tests of EGRET is much less than the predicted level. This discrepancy appears to be due to the FLUKA87 inability to transport evaporation nucleons. It is predicted that the background in EGRET will be no more than 4-10 percent of the extragalactic diffuse gamma radiation.

Edwards, H. L.↗
NASA NTRS1992

Theoretical studies of the outer envelopes of young stellar objects

With the Monte Carlo code developed by Whitney and Hartmann, a series of models was computed of scattering in disks around young stellar objects. The code calculates scattering by dust, including polarization, in arbitrary geometries. By computing model images, it was found that disk, by themselves, around young stellar objects would be very difficult to detect with present day imaging techniques. In comparing these images to observations of young stellar objects which show diffuse structure, little resemblance was found. A flared disk system will only give high polarization when viewed edge-on, and the position angle is always oriented perpendicular to the disk plane. This suggests that an envelope, perhaps the remnant infalling envelope, must be present to scatter more stellar light than a disk can, and obscure the star at many inclinations. A grid was computed of models of scattering in a disk+envelope system. Evidence is presented that the wind of the pre-main sequence object FU Orionis arises from the surface of the luminous prostellar accretion disk. A disk wind model calculated assuming radiative equilibrium explains the differential behavior of the observed asymmetrical absorption line profiles. The model predicts that strong lines should be asymmetric and blueshifted, while weak lines should be symmetric and doubled peaked due to disk rotation, in agreement with observations.

Hartmann, Lee↗
NASA NTRS1992

Multiple scattering effects on spaceborne lidar

A semianalytic Monte Carlo code originally developed for oceanographic calculations (Poole et al., 1981) has been modified for use in studying multiple scattering of space-based lidar. The approach is very similar to that described by Kunkel and Weinman (1976). The trajectory of each photon is followed from the transmitter through multiple scattering until the photon is either scattered backward out of the atmosphere, scattered forward into the ground and absorbed, or scattered out the sides of the cloud. The probability that the photon will return directly to the detector is computed and summed over all significant scattering events within the field of view of the detector. Multiple scattering of the lidar pulse causes an apparent increase in the transmittance of the medium. Multiple scattering effects for space-based lidar are more significant than for ground-based lidar due to the much larger beam diameter in the atmosphere. These larger diameters are due not only to the greater range between the lidar and the scattering volume, but also the need to maintain relatively large beam divergences to satisfy eye safety restrictions on the laser irradiance at the Earth's surface. The simulations presented here are for a wavelength of 1064 nm and the Deirmendjian C1 phase function, which yields an extinction coefficient of 17.259/km. We have looked at two cases: a space-based lidar at 296 km observing a C1 cloud 293 km from the lidar and, for comparison purposes, a ground-based lidar looking at a C1 cloud with a base height of either 2 km or 5 km. The C1 size distribution roughly approximates that of stratocumulus or altocumulus clouds (aufm Kampe and Weickmann, 1957).

Winker, David M.↗
NASA NTRS1992

Laboratory scaled simulation of lidar cloud sounding experiments

The results of lidar measurements carried out on laboratory scale models of clouds are presented. Measurements on laboratory scale models are important since one has the knowledge of the relevant parameters of the diffusing medium, such as: scattering and absorption coefficients, phase function, homogeneity, shape, etc. Knowledge of these parameters enables one to use the results to test the reliability of theoretical and numerical investigations. To obtain a laboratory scaled model of a lidar system sounding a cloud, it is necessary to scale down all the geometrical quantities by the same factor to reduce distances of the order of kilometers to the order of meters, keeping the size and the optical depth of the diffusers unchanged. If a time resolution of the order of nanoseconds is necessary for a lidar sounding actual clouds, the corresponding time resolution for the laboratory model should be of the order of picoseconds. It is possible to obtain this resolution by using picosecond laser systems and fast electrooptical detectors like the streak camera. The results of the laboratory measurements showed that the multiple scattering effect strongly depends on the size of the diffusers, as well as on the concentration. The experimental results were compared with the numerical results of a Monte Carlo code. A generally good agreement was obtained.

Zaccanti, G.↗
NASA NTRS1992

Investigation of gas/particle heat transfer rates in solid rocket motors

The ability of the current Nusselt number prediction technique developed by Kavanau (1955) to accurately predict alumina particle heat transfer rates in solid rocket nozzles and plumes is investigated. For the solid rocket motors SRMS) analyzed, the transitional regime is the dominant regime for the majority of particles in the flowfield. The analytical approach to determine accuracy of the Kavanau correlation utilized the G2R Direct Simulation Monte Carlo code. With this method, both sphere drag, and heat transfer rates were predicted. The sphere drag prediction were compared to the Hermsen, and Henderson drag correlations, while the heat transfer results were compared to the current theory. Results have indicated, that the predicted drag coefficient is bounded by the drag correlations considered. However, the Nusselt number varies significantly from the extrapolated profile through all flowfield regimes.

Moylan, B.↗
NASA NTRS1992

The atmospheric structures of the companion stars of eclipsing binary x ray sources

This investigation was aimed at determining structural features of the atmospheres of the massive early-type companion stars of eclipse x-ray pulsars by measurement of the attenuation of the x-ray spectrum during eclipse transitions and in deep eclipse. Several extended visits were made to ISAS in Japan by G. Clark and his graduate student, Jonathan Woo to coordinate the Ginga observations and preliminary data reduction, and to work with the Japanese host scientist, Fumiaki Nagase, in the interpretation of the data. At MIT extensive developments were made in software systems for data interpretation. In particular, a Monte Carlo code was developed for a 3-D simulation of the propagation of x-rays from the neutron star through the ionized atmosphere of the companion. With this code it was possible to determine the spectrum of Compton-scattered x-rays in deep eclipse and to subtract that component from the observed spectra, thereby isolating the software component that is attributable in large measure to x-rays that have been scattered by interstellar grains. This research has culminated in the submission of paper to the Astrophysical Journal on the determination of properties of the atmosphere of QV Nor, the BOI companion of 4U 1538-52, and the properties of interstellar dust grains along the line of sight from the source. The latter results were an unanticipated byproduct of the investigation. Data from Ginga observations of the Magellanic binaries SMC X-1 and LMC X-4 are currently under investigation as the PhD thesis project of Jonathan Woo who anticipated completion in the spring of 1993.

Clark, George W.↗
NASA NTRS1992

Effects of convection electric field on upwelling and escape of ionospheric O(+)

A Monte Carlo code is used to explore the full effects of the convection electric field on distributions of upflowing O(+) ions from the cusp/cleft ionosphere. Trajectories of individual ions/neutrals are computed as they undergo multiple charge-exchange collisions. In the ion state, the trajectories are computed in realistic models of the magnetic field and the convection, corotation, and ambipolar electric fields. The effects of ion-ion collisions are included, and the trajectories are computed with and without simultaneous stochastic heating perpendicular to the magnetic field by a realistic model of broadband, low frequency waves. In the neutral state, ballistic trajectories in the gravitational field are computed. The initial conditions of the ions, in addition to ambipolar electric field and the number densities and temperatures of O(+), H(+), and electrons as a function of height in the cusp/cleft region were obtained from the results of Gombosi and Killeen (1987), who used a hydrodynamic code to simulate the time-dependent frictional-heating effects in a magnetic tube during its motion though the convection throat. The distribution of the ion fluxes as a function of height are constructed from the case histories.

Cladis, J. B.↗
NASA NTRS1992

Particle acceleration in modified oblique non-relativistic shocks

Previous applications of the Monte Carlo technique at the quasiparallel earth bow shock has motivated the extension of this technique to oblique shock geometries typical of those found in most astrophysical shock environments. In addition, such a generalization will permit the thorough examination of theoretical predictions of rapid acceleration times at quasi-perpendicular shocks. Therefore, we have embarked on the modification of our existing Monte Carlo code and in this paper outline the major technical aspects involved in developing a simulation of cosmic-ray acceleration at modified oblique nonrelativistic shocks.

Baring, Matthew G.↗