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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Double Compton scatter telescope calibration

Calibration techniques for a medium energy gamma ray telescope are described. Gain calibration using Compton edge spectra involves comparisons of pulse height spectra with spectra simulated by a Monte Carlo computer code which includes Compton scattering and pair production, plural scattering and variable energy resolution, and cell size. The telescope considered comprises 56 cells of liquid scintillator in four size groups, with a total liquid volume of 325 l; each cell has its own photomultiplier tube. Energy and angular resolutions and the PMT gain calibration procedure are verified with double scatter data for monoenergetic gamma rays at a known location. Detection probabilities for any cell combination in the two telescope arrays are calculated per steradian as a function of the scattering for a number of different energies with a Van de Graaff accelerator.

Dayton, B.↗
NASA NTRS1985

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

The polar ionosphere as a source of the storm time ring current

The transport of ions from the polar ionosphere to the inner magnetosphere during storm time conditions has been computed using a Monte Carlo diffusion code. The effect of the electrostatic turbulence assumed to be present during the substorm expansion phase was simulated by a process that accelerated the ions stochastically perpendicular to the magnetic field with a diffusion coefficient proportional to the rate of energization of the ions by the induced electric field. This diffusion process was continued as the ions were convected from the plasma sheet boundary layer to the double-spiral injection boundary. Inward of the injection boundary the ions were convected adiabatically. By using as input an O(+) flux of 2.8 x 10 to the 8th per sq cm per s (w greater than 10 eV) and an H(+) flux of 5.5 x 10 to the 8th per sq cm per s (w greater than .63 eV) the computed distribution functions of the ions in the ring current were found to be in good agreement, over a wide range in L (4-8), with measurements made with the ISEE 1 satellite during a storm. This O(+) flux and a large part of the H(+) flux appear to be consistent with the DE-1 and DE-2 satellite measurements of the polar ionospheric outflow during disturbed times.

Cladis, J. B.↗
NASA NTRS1985

Transport of ions in presence of induced electric field and electrostatic turbulence - Source of ions injected into ring current

The transport of ions from the polar ionosphere to the inner magnetosphere during stormtime conditions has been computed using a Monte Carlo diffusion code. The effect of the electrostatic turbulence assumed to be present during the substorm expansion phase was simulated by a process that accelerated the ions stochastically perpendicular to the magnetic field with a diffusion coefficient proportional to the energization rate of the ions by the induced electric field. This diffusion process was continued as the ions were convected from the plasma sheet boundary layer to the double-spiral injection boundary. Inward of the injection boundary, the ions were convected adiabatically. By using as input an O(+) flux of 2.8 x 10 to the 8th per sq cm per s (w greater than 10 eV) and an H(+) flux of 5.5 x 10 to the 8th per sq cm per s (w greater than 0.63 eV), the computed distribution functions of the ions in the ring current were found to be in good agreement, over a wide range in L (4 to 8), with measurements made with the ISEE-1 satellite during a storm. This O(+) flux and a large part of the H(+) flux are consistent with the DE satellite measurements of the polar ionospheric outflow during disturbed times.

Cladis, J. B.↗
NASA NTRS1992

X-ray scattering and fluorescence in the wind of a massive X-ray binary

Spectral data from the binary X-ray pulsar 4 U 0900 - 40 obtained with the Ginga satellite and covering a range of orbital phases are presented and interpreted using simulated spectra created with a Monte Carlo scattering code. It is found that scattering and fluorescence in a simple spherically symmetric distribution of matter surrounding the companion star can reproduce the eclipse spectrum, as well as the 'soft X-ray-excess' observed during egress and other orbital phases. Reasonably secure values are found for a number of the parameters that characterize the density profile of the stellar wind and atmosphere of the companion star. The Fe abundance is within a factor of about 1.3 of that in the solar neighborhood. It is shown that ionization zones are not critical to understanding the orbital-phase-resolved spectra in this source. It is also found that the contribution by scattering from interstellar dust grains to the observed spectra during eclipse is negligible, while that from diffuse emission from the 'Galactic ridge' is significant.

Lewis, W.↗
NASA NTRS1994

Rates for neutron-capture reactions on tungsten isotopes in iron meteorites

High-precision W isotopic analyses by Harper and Jacobsen indicate the W-182/W-183 ratio in the Toluca iron meteorite is shifted by -(3.0 +/- 0.9) x 10(exp -4) relative to a terrestrial standard. Possible causes of this shift are neutron-capture reactions on W during Toluca's approximately 600-Ma exposure to cosmic ray particles or radiogenic growth of W-182 from 9-Ma Hf-182 in the silicate portion of the Earth after removal of W to the Earth's core. Calculations for the rates of neutron-capture reactions on W isotopes were done to study the first possibility. The LAHET Code System (LCS) which consists of the Los Alamos High Energy Transport (LAHET) code and the Monte Carlo N-Particle(MCNP) transport code was used to numerically simulate the irradiation of the Toluca iron meteorite by galactic-cosmic-ray (GCR) particles and to calculate the rates of W(n, gamma) reactions. Toluca was modeled as a 3.9-m-radius sphere with the composition of a typical IA iron meteorite. The incident GCR protons and their interactions were modeled with LAHET, which also handled the interactions of neutrons with energies above 20 MeV. The rates for the capture of neutrons by W-182, W-183, and W-186 were calculated using the detailed library of (n, gamma) cross sections in MCNP. For this study of the possible effect of W(n, gamma) reactions on W isotope systematics, we consider the peak rates. The calculated maximum change in the normalized W-182/W-183 ratio due to neutron-capture reactions cannot account for more than 25% of the mass 182 deficit observed in Toluca W.

Masarik, J.↗
NASA NTRS1994

Numerical simulation of cosmogenic nuclide production in lunar rocks

The production rates of cosmogenic nuclides depend on the primary cosmic-ray particles, the irradiated-body's bulk composition, size, and shape, and the sample's composition and shielding depth. Although much work has been done on some of these dependencies, more detailed studies still need to be done on others. This work describes the influence of irradiation geometry on nuclide production in lunar rocks. In most cases, computer simulations of cosmogenic nuclide production were restricted to spherical objects irradiated with a 4 pi isotropic flux (meteoroids) or in lunar core samples irradiated by a 2 pi flux incident on semi-infinite layers or cylinders of huge sizes. Many lunar samples are rocks found on top of the lunar surface. For these rocks, neither of the above-mentioned models correspond to the real conditions. We present results of our simulations of cosmogenic nuclide production in models simulating the irradiation of rocks sitting on top of the lunar surface. The Galactic Cosmic Rays (GCR) production profiles in lunar rocks were calculated using the Los Alamos 3-D Monte Carlo LAHET Code System (LCS). The irradiated object was modeled as the union of a sphere with the radius of the Moon and a small hemisphere with radii varying from 10 to 100 g/sq cm simulating the lunar rock. These calculations for the production of cosmogenic nuclides in lunar rocks by GCR particle show that there are important differences between the results obtained by commonly used geometric irradiation models and the lunar-rock models presented. The steeper GCR production profiles for a rock could help to explain the poor agreement for Be-10 in rock 68815, where slab models give GCR profiles flatter than the observed profiles.

Reedy, R. C.↗
NASA NTRS1995

The detector response matrices of the burst and transient source experiment (BATSE) on the Compton Gamma Ray Observatory

The detector response matrices for the Burst And Transient Source Experiment (BATSE) on board the Compton Gamma Ray Observatory (CGRO) are described, including their creation and operation in data analysis. These response matrices are a detailed abstract representation of the gamma-ray detectors' operating characteristics that are needed for data analysis. They are constructed from an extensive set of calibration data coupled with a complex geometry electromagnetic cascade Monte Carlo simulation code. The calibration tests and simulation algorithm optimization are described. The characteristics of the BATSE detectors in the spacecraft environment are also described.

Pendleton, Geoffrey N.↗
NASA NTRS1997

Development of a Space Radiation Monte Carlo Computer Simulation

The ultimate purpose of this effort is to undertake the development of a computer simulation of the radiation environment encountered in spacecraft which is based upon the Monte Carlo technique. The current plan is to adapt and modify a Monte Carlo calculation code known as FLUKA, which is presently used in high energy and heavy ion physics, to simulate the radiation environment present in spacecraft during missions. The initial effort would be directed towards modeling the MIR and Space Shuttle environments, but the long range goal is to develop a program for the accurate prediction of the radiation environment likely to be encountered on future planned endeavors such as the Space Station, a Lunar Return Mission, or a Mars Mission. The longer the mission, especially those which will not have the shielding protection of the earth's magnetic field, the more critical the radiation threat will be. The ultimate goal of this research is to produce a code that will be useful to mission planners and engineers who need to have detailed projections of radiation exposures at specified locations within the spacecraft and for either specific times during the mission or integrated over the entire mission. In concert with the development of the simulation, it is desired to integrate it with a state-of-the-art interactive 3-D graphics-capable analysis package known as ROOT, to allow easy investigation and visualization of the results. The efforts reported on here include the initial development of the program and the demonstration of the efficacy of the technique through a model simulation of the MIR environment. This information was used to write a proposal to obtain follow-on permanent funding for this project.

Pinsky, Lawrence S.↗
NASA NTRS1998

Aerothermodynamics of the Mars Global Surveyor Spacecraft

The aerothermodynamics characteristics of the Mars Global Surveyor spacecraft are investigated and reported. These results have been used by the Mars Global Surveyor mission planners to design the aerobraking phase of the mission. Analytical and Direct Simulation Monte Carlo computer codes were used with a detailed, three dimensional model of the spacecraft to evaluate spacecraft aerobraking characteristics for flight in free molecular and transitional flow regimes. The spacecraft is found to be aerodynamically stable in aerobraking and planned contingency configurations. Aerodynamic forces, moments, and heating are found to be highly dependent on atmospheric density. Accommodation coefficient. is seen to strongly influence drag coefficient. Transitional flow effects are found to reduce overall solar panel heating. Attitude control thruster plumes are shown to interact with the freestream, diminishing the effectiveness of the attitude control system and even leading to thrust reversal. These plume-freestream interaction effects are found to be highly dependent on freestream density.

Shane, Russell W.↗
NASA NTRS1998

Estimation of 3-D Cloud Effects on TOMS Satellite Retrieval of Surface UV Irradiance

To improve surface UV irradiance retrieval from the Total Ozone Mapping Spectrometer (TOMS) we simulate errors of the TOMS cloud correction algorithm for summertime broken cloud conditions. Cloud scenes (50 km by 50 km) are modeled by a normal random (Gaussian) field with a fixed lower boundary and conservative scattering. The model relates stochastic field characteristics with the cloud amount, mean cloud diameter and aspect ratio. Clouds are embedded into Rayleigh atmosphere with standard ozone profile. Radiative transfer calculations of the radiance at the top of the atmosphere and irradiance at the surface were performed using 3-D Monte Carlo (MC) code. The results are averaged over the satellite field of view on the surface (50 km by 50 km) and compared with TOMS predicted surface irradiance for the same scene reflectance. The TOMS algorithm assumes horizontally homogeneous Cl-type cloud between 3 km and 5.5 km. The effective optical depth is determined by fitting observed (MC) radiance at 380 nm. Having the same radiance at the satellite the homogeneous and broken cloud models predict different average irradiances at the surface. This is due to the differences in Bidirectional Reflection Distribution Function (BRDF) for homogeneous and broken cloud scenes with the same hemispherical albedo. For typical TOMS observational geometry at mid-latitudes the simulated single pixels errors may be as large as +/- 20%. Qualitatively these errors are due to the dominance of the non-horizontal cloud surfaces, which are not accounted for in the homogeneous cloud model. However, due to high variability of the real cloud shapes and types it is unclear how these single pixel errors would affect TOMS time-integrated UV exposure over extended periods (weeks to months) for different regions.

Krotkov, Nickolay A.↗
NASA NTRS2000

X-Ray Astronomy

Dr. S. N. Zhang has lead a seven member group (Dr. Yuxin Feng, Mr. XuejunSun, Mr. Yongzhong Chen, Mr. Jun Lin, Mr. Yangsen Yao, and Ms. Xiaoling Zhang). This group has carried out the following activities: continued data analysis from space astrophysical missions CGRO, RXTE, ASCA and Chandra. Significant scientific results have been produced as results of their work. They discovered the three-layered accretion disk structure around black holes in X-ray binaries; their paper on this discovery is to appear in the prestigious Science magazine. They have also developed a new method for energy spectral analysis of black hole X-ray binaries; four papers on this topics were presented at the most recent Atlanta AAS meeting. They have also carried Monte-Carlo simulations of X-ray detectors, in support to the hardware development efforts at Marshall Space Flight Center (MSFC). These computation-intensive simulations have been carried out entirely on the computers at UAH. They have also carried out extensive simulations for astrophysical applications, taking advantage of the Monte-Carlo simulation codes developed previously at MSFC and further improved at UAH for detector simulations. One refereed paper and one contribution to conference proceedings have been resulted from this effort.

Wu, S. T.↗
NASA NTRS2001

Orbits of Saturn's Inner Moons and Other Observations Connected within the 1995-1996 Saturnian Ring Plane Crossing

We reduced and analyzed radio images of Saturn taken in 1990-1995 at 0.35, 2.0, 3.6 and 6.1 cm wavelength. The latitudinal brightness distribution varies substantially over time. Detailed modeling of these data shows that NH4SH condenses at the thermochemical equilibrium temperature of 235.5 K only at Northern mid-latitudes. Over most of the planet, condensation does not occur until the gas cools to 190 +/- 5 K. Supersaturation may also cause the dark equatorial region seen in 1995 at 6.1 cm. Observations of the rings show that the West (dusk) ansa is brighter than the East (dawn) ansa by factors up to 2. The magnitude of the asymmetry increases with increasing wavelength and with decreasing distance to the planet. This East-West asymmetry may be due to multiple scattering in gravitational wakes. We improved the Monte Carlo scattering code used to model the microwave emission from Saturn's rings at multiple wavelengths. The basic code outlined in Dunn, Molnar, and Fix describes the rings as a uniform, infinite layer of ice spheres with a finite vertical extent. The particle sizes follow a power law distribution and scatter light according to a linear combination of Mie and isotropic phase functions. These parameters may be adjusted to reflect the physical properties in the scattering properties of the ring particles. This is used to construct a model which can be directly compared to Very Large Array data. The code has been extended by including anisotropic particle distributions, primarily in the form of the wake structures that are likely to exist in the enhancements. The vertical extent of the wakes equals that of the ring layer. Between the wakes are regions of diminished density so that the optical depth integrated over many wakes matches that of the Voyager occultation experiments. This allows us to model the brightness temperature as a function of azimuth and look for asymmetric scattering behavior in the rings. By changing the width and relative density of the wakes, we can alter the nature of the gross scattering properties of the rings to match what is seen in our data.

Lissauer, Jack J.↗
NASA NTRS2001

Chemical and Flowfield Modeling for Enhanced Analysis of Contamination Experiments

This paper describes the application of a new Direct Simulation Monte Carlo (DSMC) code, the Molecular Beam Simulator (MBS), which is designed to analyze laboratory scale molecular beam-surface (and crossed-beam) experiments. The MBS is primarily intended to model experiments associated with spacecraft contamination effects, but it can also be used to simulate a variety of surface chemistry and reactive flow measurements. The MBS code is fully three-dimensional, includes a wide-range of chemical processes, and can model one or multiple pulsed (non-steady) sources. As an example application of the MBS code, a fast, pulsed, oxygen atom-surface experiment which examines the chemistry behind erosion of graphite by oxygen atoms is analyzed. Unsteady DSMC simulations show that experimental observations of excited molecular states after the pulse has hit the surface are consistent with two distinct chemical mechanisms: a direct one where the excited molecules are formed on the surface, and a two-step mechanism where ground state molecules formed on the surface are collisionally excited after they leave the surface by trailing oxygen atoms in the pulse. Further DSMC calculations suggest experiments which can distinguish between these mechanisms.

Braunstein, Matthew↗
NASA NTRS2003

Monte Carlo Methodology Serves Up a Software Success

Widely used for the modeling of gas flows through the computation of the motion and collisions of representative molecules, the Direct Simulation Monte Carlo method has become the gold standard for producing research and engineering predictions in the field of rarefied gas dynamics. Direct Simulation Monte Carlo was first introduced in the early 1960s by Dr. Graeme Bird, a professor at the University of Sydney, Australia. It has since proved to be a valuable tool to the aerospace and defense industries in providing design and operational support data, as well as flight data analysis. In 2002, NASA brought to the forefront a software product that maintains the same basic physics formulation of Dr. Bird's method, but provides effective modeling of complex, three-dimensional, real vehicle simulations and parallel processing capabilities to handle additional computational requirements, especially in areas where computational fluid dynamics (CFD) is not applicable. NASA's Direct Simulation Monte Carlo Analysis Code (DAC) software package is now considered the Agency s premier high-fidelity simulation tool for predicting vehicle aerodynamics and aerothermodynamic environments in rarified, or low-density, gas flows.

Source record↗
NASA NTRS2003

Modeling the Extended Neutral Atmosphere and Plasma Environment near Saturn

In the three years of this study we have published five papers in refereed journals. We have first examined satellite sources and their contribution to the observed neutral cloud. Based on the total calculated satellite sources and the spatial distribution of sputtered H20, we concluded that they cannot produce observed OH cloud. The E-ring contribution has been also studied in detail.In order to produce observed OH cloud we suggested that the E-ring might be the dominant source in inner Saturnian magnetosphere. We proposed a possible resupply mechanism which is needed to keep both E-ring and OH cloud in the present state: collisions between E-ring grains and remains of a disrupted satellite near Enceladus Lagrangian points. In this scenario a large amount of material, ranging from a few micrometers to hundred of meters, which is optically invisible at present, is likely to exist there. The fourth paper compares the magnetosheaths of the outer planets. A surprising result is that the hot proton component comprises about 40% of the total density, much larger than predicted by shock theory. Gas dynamic models of the boundaries show that the magnetospheres of Jupiter and Saturn are flattened at the poles. The last paper was published in GRL and is the first based of the model of neutrals developed as a main goal of this grant and which is now operational.This Monte Carlo collision code self- consistently determines the neutral distribution. from the rings and satellites until they are lost by ionization, by collisions with rings, moons, or Saturn, or by escape from Saturn. Our model is unique in that it includes the effects of plasma chemistry and both plasma-neutral and neutral neutral collisions to determine the dynamical evolution of the water group neutrals in Saturn's magnetosphere. The lifetimes of the neutrals against loss to photoionization, charge exchange, electron dissociation and electron impact dissociation are based on the model given by (Richardson et al. 1998) and vary with position in the magnetosphere. The dominant neutral dissociation channels H20->OH+H, H20-> O+H2, and OH->O+H2 are considered.

Richardson, John D.↗
NASA NTRS2004

The Role of Beam Geometry in Population Statistics and Pulse Profiles of Radio and Gamma-ray Pulsars

We present results of a pulsar population synthesis study that incorporates a number of recent developments and some significant improvements over our previous study. We have included the results of the Parkes multi-beam pulsar survey in our select group of nine radio surveys, doubling our sample of radio pulsars. More realistic geometries for the radio and gamma-ray beams are included in our Monte Carlo computer code that simulates the characteristics of the Galactic population of radio and gamma-ray pulsars. We adopted with some modifications the radio beam geometry of Arzoumanian, Chernoff & Cordes (2002). For the gamma-ray beam, we have assumed the slot gap geometry described in the work of Muslimov & Harding (2003). To account for the shape of the distribution of radio pulsars in the P(dot) - P diagram, we continue to find that decay of the magnetic field on a timescale of 2.8 Myr is needed. With all nine surveys, our model predicts that EGRET should have seen 7 radio-quiet (below the sensitivity of these radio surveys) and 19 radio-loud gamma-ray pulsars. AGILE (nominal sensitivity map) is expected to detect 13 radio-quiet and 37 radio-loud gamma-ray pulsars, while GLAST, with greater sensitivity is expected to detect 276 radio-quiet and 344 radio-loud gamma-ray pulsars. When the Parkes multi-beam pulsar survey is excluded, the ratio of radio-loud to radio-quiet gamma-ray pulsars decreases, especially for GLAST. The decrease for EGRET is 45%, implying that some fraction of EGRET unidentified sources are radio-loud gamma-ray pulsars. In the radio geometry adopted, short period pulsars are core dominated. Unlike the EGRET gamma-ray pulsars, our model predicts that when two gamma-ray peaks appear in the pulse profile, a dominant radio core peak appears in between the gamma-ray peaks. Our findings suggest that further improvements are required in describing both the radio and gamma-ray geometries.

Gonthier, Peter L.↗
NASA NTRS2004

Radiolytic Model for Chemical Composition of Europa's Atmosphere and Surface

The overall objective of the present effort is to produce models for major and selected minor components of Europa s neutral atmosphere in 1-D versus altitude and in 2-D versus altitude and longitude or latitude. A 3-D model versus all three coordinates (alt, long, lat) will be studied but development on this is at present limited by computing facilities available to the investigation team. In this first year we have focused on 1-D modeling with Co-I Valery Shematovich s Direct Simulation Monte Carlo (DSMC) code for water group species (H2O, O2, O, OH) and on 2-D with Co-I Mau Wong's version of a similar code for O2, O, CO, CO2, and Na. Surface source rates of H2O and O2 from sputtering and radiolysis are used in the 1-D model, while observations for CO2 at the Europa surface and Na detected in a neutral cloud ejected from Europa are used, along with the O2 sputtering rate, to constrain source rates in the 2-D version. With these separate approaches we are investigating a range of processes important to eventual implementation of a comprehensive 3-D atmospheric model which could be used to understand present observations and develop science requirements for future observations, e.g. from Earth and in Europa orbit. Within the second year we expect to merge the full water group calculations into the 2-D version of the DSMC code which can then be extended to 3-D, pending availability of computing resources. Another important goal in the second year would be the inclusion of sulk and its more volatile oxides (SO, SO2).

Cooper, John F.↗