Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo codes”

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 NTRS2008

Benchmark Analysis of Pion Contribution from Galactic Cosmic Rays

Shielding strategies for extended stays in space must include a comprehensive resolution of the secondary radiation environment inside the spacecraft induced by the primary, external radiation. The distribution of absorbed dose and dose equivalent is a function of the type, energy and population of these secondary products. A systematic verification and validation effort is underway for HZETRN, which is a space radiation transport code currently used by NASA. It performs neutron, proton and heavy ion transport explicitly, but it does not take into account the production and transport of mesons, photons and leptons. The question naturally arises as to what is the contribution of these particles to space radiation. The pion has a production kinetic energy threshold of about 280 MeV. The Galactic cosmic ray (GCR) spectra, coincidentally, reaches flux maxima in the hundreds of MeV range, corresponding to the pion production threshold. We present results from the Monte Carlo code MCNPX, showing the effect of lepton and meson physics when produced and transported explicitly in a GCR environment.

Aghara, Sukesh K.↗
NASA NTRS1989

Computer simulation of dust grain evolution

The latest results are reported from a Monte Carlo code that is being developed at NASA Ames. The goal of this program, is to derive from the observed and presumed properties of the interstellar medium (ISM) the following information: (1) the size spectrum of interstellar dust; (2) the chemical structure of interstellar dust; (3) interstellar abundances; and (4) the lifetime of a dust grain in the ISM. Presently this study is restricted to refractory interstellar material, i.e., the formation and destruction of ices are not included in the program. The program is embedded in an analytic solution for the bulk chemical evolution of a two-phase interstellar medium in which stars are born in molecular clouds, but new nucleosynthesis products and stellar return are entered into a complementary intercloud medium. The well-mixed matter of each interstellar phase is repeatedly cycled stochastically through the complementary phase and back. Refractory dust is created by thermal condensation as stellar matter flows away from sites of nucleosynthesis such as novae and supernovae and/or from the matter returned from evolved intermediate stars. The history of each particle is traced by standard Monte Carlo techniques as it is sputtered and fragmented by supernova shock waves in the intercloud medium. It also accretes an amorphous mantle of gaseous refractory atoms when its local medium joins with the molecular cloud medium. Finally it encounters the possibility of astration (destruction by star formation) within the molecular clouds.

Liffman, K.↗
NASA NTRS1962

Monte Carlo Calculations of Neutron Number Spectra and Buildup Factors in Infinite Conical Configurations

A Monte Carlo code simulating neutron transport in infinite cones of water and water-equivalent hydrogen was prepared for an IBM 704 computer. The code was essentially a modification of the point-source, infinite-medium code used in NASA TN D-850. Studies were made of differential neutron number spectra and associated buildup factors for infinite cones having apex half-angles of 15 degrees, 30 degrees, 45 degrees, and 60 degrees. The buildup factors obtained were compared with those for the appropriate infinite medium, which allowed an examination of the effect of solid angle subtended by material on the transport of 6-Mev source neutrons emanating from the cone apex. The variation of number buildup factor with distance for the various cones shows that neutron scattering out of the cones is predominant in the first 30 t o 40 centimeters of material, and that transport beyond this distance is of a similar nature in all the cones.

CONE↗
NASA NTRS1972

MORSE Monte Carlo shielding calculations for the zirconium hydride reference reactor

Verification of DOT-SPACETRAN transport calculations of a lithium hydride and tungsten shield for a SNAP reactor was performed using the MORSE (Monte Carlo) code. Transport of both neutrons and gamma rays was considered. Importance sampling was utilized in the MORSE calculations. Several quantities internal to the shield, as well as dose at several points outside of the configuration, were in satisfactory agreement with the DOT calculations of the same.

Burgart, C. E.↗
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↗