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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
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.↗
NASA NTRS1993

Cosmological gamma-ray burst model

A cosmological gamma-ray burst model that reproduces the observed gamma-ray spectra was developed. This model, which is an outgrowth of work on synchrotron emission from cosmological sources, creates the observed spectra from a power law spectrum through Compton attenuation in dense molecular clouds. It restricts the burst source to the centers of galaxies, and it is easily tested through comparisons with time dependent burst spectra. Dr. Brainerd continued development of a Monte Carlo code that tracks the random walk of a gamma-ray through a plasma in a strong magnetic field and used this code to show that a two temperature plasma in a strong magnetic field can suppress soft x-ray emission.

Brainerd, J. J.↗
NASA NTRS1994

Diffuse interstellar bands in reflection nebulae

A Monte Carlo code for radiation transport calculations is used to compare the profiles of the lambda lambda 5780 and 6613 Angstrom diffuse interstellar bands in the transmitted and the reflected light of a star embedded within an optically thin dust cloud. In addition, the behavior of polarization across the bands were calculated. The wavelength dependent complex indices of refraction across the bands were derived from the embedded cavity model. In view of the existence of different families of diffuse interstellar bands the question of other parameters of influence is addressed in short.

Fischer, O.↗
NASA NTRS1995

Photochemical modeling of Titan's atmosphere

We have developed a new photochemical model of Titan's atmosphere which includes all the important compounds and reactions in spherical geometry from the surface to 1240 km. Compared to the previous model of Yung et al. (1984), the most significant recent change in the reactions used is the updated methane dissociation scheme (Mordaunt et al. 1993). Moreover, the transfer of the solar radiation in the atmosphere and the photolysis rates have been calculated by using a Monte Carlo code. Finally, the eddy diffusion coefficient profile is adjusted in order to fit the mean vertical distribution of HCN retrieved from millimeter groundbased observations of Tanguy et al. (1990); using new values for the boundary flux of atomic nitrogen (Strobel et al. 1992). We have run the model in both steady-state and diurnal modes, with 62 species involved in 249 reactions. There is little difference between diurnal and steady-state results. Overall our results are in a closer agreement with the abundances inferred from the Voyager infrared measurements at the equator than the Yung et al. results. We find that the catalytic scheme for H recombination invoked by Yung et al. only slightly improves the model results and we conclude that this scheme is not essential to fit observations.

Toublanc, D.↗
NASA NTRS1996

Dust Around Herbig Ae Stars: Additional Constraints from their Photometric and Polarimetric Variability

For the Herbig Ae stars with Algol-like minima (UX Ori, WW Vul, etc), the effects of circumstellar dust include: excess infrared emission, anomalous ultraviolet extinction, the 'blueing' of the stars in minima accompanying by an increase of intrinsic polarization. Using a Monte-Carlo code for polarized radiation transfer we have simulated these effects and compared the results obtained for different models with the observational data available. We found that the photometric and polarimetric behavior of the stars provided essential additional constraints on the circumstellar dust models. The models with spheroidal shell geometry and compact (non-fluffy) dust grains do not appear to be able to explain all the data.

Krivova, N. A.↗
NASA NTRS1997

Properties of accretion disk coronae

The properties of accretion disk corona in a parameter regime suitable for Galactic black hole candidates are considered and the results of an analysis of these properties using a self-consistent Monte Carlo code are presented. Examples of the coronal temperature structure, the shape and angular dependency of the spectrum and the maximum temperature allowed for each optical depth of the corona are presented. It is shown that the observed spectrum of the Galactic black hole candidate Cygnus X-1 cannot be explained by accreting disk corona models with a slab geometry, where the accretion disk is sandwiched by the comptonizing medium.

Wilms, J.↗
NASA NTRS1997

Assay of the Martian Regolith with Neutrons

The purpose of the research is to combine experiments and Monte Carlo transport of neutrons through volume of soil in an attempt to model neutron leakage from planetary surfaces. Emphasis is given to the change of neutron spectra as a function of water content and location. During the first stage of effort, two experiments were conducted in which leakage of neutrons from a Pu-Be source through about 30 g/cm(exp 2) of soil were measured with several counters. A Monte Carlo code, MCNP, has been used to model many of the 100 individual runs of the experiment. Hydrogen is the element that has the most dramatic effect on the neutron spectrum and its effect on the neutron spectrum is almost the same whether it is in the form of water or polyethylene. In order to simulate various water configurations, sheets of polyethylene have been used between layers of soil as well as water in several concentrations up to 18%. Comparison of experimental results to theoretical predictions made with the MCNP code were disappointing for low concentrations of water. We have made extensive calculations to see if room return could be the cause of the discrepancies. Water concentrations of the 'dry' soil were measured by two different laboratories and differed only by 0.5%. We have made calculations to optimize the next experiment and are investigating other methods of determining the water content of 'dry' soil.

Drake, Darrell M.↗
NASA NTRS1997

Self-Consistent Thermal Accretion Disk Corona Models for Compact Objects. I: Properties of the Corona and the Spectrum of Escaping Radiation

We present the properties of accretion disk corona (ADC) models in which the radiation field, the temperature, and the total opacity of the corona are determined self-consistently. We use a nonlinear Monte Carlo code to perform the calculations. As an example, we discuss models in which the corona is situated above and below a cold accretion disk with a plane-parallel (slab) geometry, similar to the model of Haardt & Maraschi. By Comptonizing the soft radiation emitted by the accretion disk, the corona is responsible for producing the high-energy component of the escaping radiation. Our models include the reprocessing of radiation in the accretion disk. Here the photons either are Compton-reflected or photoabsorbed, giving rise to fluorescent line emission and thermal emission. The self- consistent coronal temperature is determined by balancing heating (due to viscous energy dissipation) with Compton cooling, determined using the fully relativistic, angle-dependent cross sections. The total opacity is found by balancing pair productions with annihilations. We find that, for a disk temperature kT(sub BB) approx. less than 200 eV, these coronae are unable to have a self-consistent temperature higher than approx. 140 keV if the total optical depth is approx. less than 0.2, regardless of the compactness parameter of the corona and the seed opacity. This limitation corresponds to the angle-averaged spectrum of escaping radiation having a photon index approx. greater than 1.8 within the 5-30 keV band. Finally, all models that have reprocessing features also predict a large thermal excess at lower energies. These constraints make explaining the X-ray spectra of persistent black hole candidates with ADC models very problematic.

Dove, James B.↗
NASA NTRS1998

Mars Pathfider Rarefied Aerodynamics: Computations and Measurements

On July 4, 1997, after traveling close to 500 million km, the Pathfinder spacecraft successfully completed entry, descent, and landing at Mars. In the present paper, the focus is on the hypersonic rarefied portion of Pathfinder's atmospheric entry where the synergy of flight measurements, aerodynamic calculations, and atmospheric modeling tools are used to extract Pathfinder's attitude and the freestream density. Accuracy of the capsule aerodynamics directly impacts the inferred atmospheric properties extracted from deceleration measurements made by on-board accelerometers. The range of rarefaction considered in this study extends from the free molecular to continuum conditions and angles of attack from O to 30 deg. The aerodynamic computations are made with free-molecular and direct simulation Monte Carlo codes. The calculations show that Pathfinder is statically unstable for much of the transitional rarefied regime. Due to the relatively modest forces and the gyroscopic motion of the spacecraft, the angle of attack excursions were less than 5 deg as inferred from force measurements for the rarefied portion of entry and approached a nominal zero degree trim angle near hypersonic continuum conditions.

Moss, James N.↗
NASA NTRS1996

The Response of the BATSE LADs to Radiation From the Crab Nebula and Plans for Radioactivity Studies on Space Station

The Burst And Transient Source Experiment (BATSE) onboard the Compton Gamma-Ray Observatory (CGRO) was designed to measure X-rays and gamma rays with energies from about 50 keV to above 2 MeV. As with many scientific investigations, the success of the original experiment lead to additional areas of research interest. In the case of BATSE the ability to observe the radiation from sources down to about 20 keV became readily apparent. This lead to a continuing program of measuring the spectrum of radiation from stellar objects at these lower energies. One of these, the Crab Nebula, has a very steady radiation flux and, thus, has become a "standard candle" for such measurements. The Large Area Detectors (LADS) on BATSE contain a 1.27-cm thick, 25.4-cm radius NaI(Tl) detector behind a 6.35-mm thick polystyrene Charged Particle Detector (CPD) used to "veto" charged particles signals. The detectors have been calibrated with a series of gamma and X-ray sources and the results carefully simulated with a Monte Carlo code. In the calibration process the computer simulation accounts for scattering from material in the counting room as well as the BATSE structure. For an orbiting detector, scattering from the entire spacecraft must be modeled as well as for all covering material over the detectors. Five years after CGRO was launched on April 5, 1991, a large body of observational data has been taken of the Crab Nebula. The technique used for these observations, and for many other X-ray sources, is Earth occultation. From the perspective of the spacecraft, the Earth occults most stellar objects once in orbit, i.e., the signal is lost as the source sets and is regained as the source rises. A careful analysis of the continuing signals from all sources measured allows for an accurate measurement of the spectrum of a given source. An analysis of this data from the Crab has indicated that the LADs are very responsive at energies as low as 20 keV-at energies below the range of calibration. While the model accounts for many of the interactions of the photons with the detectors, the observation of nonstatistical deviations at low energy and at small angles has suggested a need to recalibrate at energies where the attenuation effects are increasing exponentially.

Laird, C. E.↗
NASA NTRS1998

Assay of the Martian Regolith with Neutrons

Different aspects of assaying Martian regolith using neutrons have been investigated. The epithermal portion of moderated neutrons spectra is dramatically effected by the presence of hydrogen (usually in the form of water). A simple analytic formula has been derived to describe the amplitude of this portion of the neutron spectrum as a function of water concentration. Several demonstration experiments have been performed and modeled with a Monte Carlo code. Results of these experiments generally agreed with the calculations to within 20%. In addition to He-3 detectors, lithium-glass scintillators and U-238 fission ion chambers were investigated to determine their applicability to space experiments.

Drake, Darrell M.↗
NASA NTRS2000

An Analysis of Water Line Profiles in Star Formation Regions Observed by SWAS

We present spectral line profiles for the 557 GHz 1(sub 1,0) yields 1(sub 0,1) ground-state rotational transition of ortho-H2(16)O for 18 galactic star formation regions observed by SWAS. 2 Water is unambiguously detected in every source. The line profiles exhibit a wide variety of shapes, including single-peaked spectra and self-reversed profiles. We interpret these profiles using a Monte Carlo code to model the radiative transport. The observed variations in the line profiles can be explained by variations in the relative strengths of the bulk flow and small-scale turbulent motions within the clouds. Bulk flow (infall, outflow) must be present in some cloud cores, and in certain cases this bulk flow dominates the turbulent motions.

Ashby, Matthew L. N.↗
NASA NTRS2002

The Evolution of Disks and Winds in Dwarf Nova Outbursts - FUSE

This project was a project to study the FUV spectra of two proto-typical dwarf novae, U Gem and SS Cygni, through an outburst cycle. The luminosity of the boundary layer in the two systems, as evidenced by earlier EUVE observations, is different in the two systems. Our intensive study of the two systems was intended to (1) probe the ionization and kinematic structure of the wind as a function of system brightness, (2) isolate the contributions of the disk to the FUV spectra, and (3) examine physical conditions and abundances of material just being accreted onto the disk from the secondary. The U Gem and SS Cyg observations took place in March and October 2000, respectively. The data obtained with FUSE was of excellent quality. Analysis of the both observations is now essentially complete, although some modeling of the SS Cyg spectra is ongoing, as we complete an ApJ manuscript on this object. Our main results for U Gem are as follows: The plateau spectra have continuum shapes and fluxes that are approximated by steady state accretion disk model spectra with an accretion rate 7x10(exp 9) Msolar/yr. The spectra also show numerous absorption lines of H I, He II, and 2-5 times ionized transitions of C, N, O, P, S, and Si. There are no emission features in the spectra, with the possible exception of a weak feature on the red wing of the 0 VI doublet. The absorption lines are narrow (FWHM approx. 50 km/s), too narrow to arise from the disk photosphere, and at low velocities (less than or equal to 700 km/s). The S VI and O VI doublets are optically thick. The absorption lines in the plateau spectra show orbital variability: in spectra obtained at orbital phases between 0.53 and 0.79, low-ionization absorption lines appear and the central depths of the preexisting lines increase. The increase in line absorption occurs at the same orbital phases as previously observed EUV and X-ray light-curve dips. If the absorbing material is in (near-) Keplerian rotation around the disk, it must be located at large disk radii. The final observation occurred when U Gem was about 2 mag from optical quiescence. The spectra are dominated by emission from an approx. 43,000 K, metal-enriched white dwarf (WD). The inferred radius of the WD is 4.95x10(exp 8) cm, close to that observed in quiescence. Allowing for a hot heated region on the surface of the WD improves the fit to the spectrum at short (less than 960 A) wavelengths. Our main results for SS Cyg are as follows: The first two of four observations of SS Cyg show disk dominated spectra with accretion rates of order 10(exp -8) Msolar/yr. Except for narrow interstellar features (atomic and molecular H), the lines are all broad consistent with a disk or wind origin. The O VI line in the spectra is mostly of wind origin as detailed modeling with our Monte Carlo code (developed in part using funds from this project) show. The continua from spectra in observations 3 and 4, observed during the decline phase, are not well fit with steady-state disks, and show considerable resemblance to quiescent spectra obtained with HUT. The most probable interpretation for the emission features seen in the spectrum in the last two observations is that they arise from a photo-illuminated choronosphere above the disk, rather than a wind.

Long, Knox↗
NASA NTRS2004

The Energy Spectra of Heavy Nuclei Measured by the ATIC Experiment

ATIC (Advanced Thin Ionization Calorimeter) is a balloon-borne experiment to measure the spectra and composition of primary cosmic rays in the region of total energy from 100 GeV to near 100 TeV for Z from 1 to 26. ATIC consists of a pixelated silicon matrix detector to measure charge plus a fully active BGO calorimeter, to measure energy, located below a carbon target interleaved with three layers of scintillator hodoscope. The ATIC instrument had a second (scientific) flight from McMurdo, Antarctica from 12/29/02 to 1/18/03, yielding 20 days of good data. The GEANT 3.21 Monte Carlo code with the QGSM event generator and the FLUKA code with the DPMJET-II event generator were used to convert energy deposition measurements to primary energy. We present the preliminary energy spectra for the abundant elements C, O, Ne, Mg, Si and Fe and compare them with the results of the first (test) flight of ATIC in 2000-01 and with results from the HEAO-3 and CRN experiments.

Panov, A. D.↗
NASA NTRS2004

MGGPOD: a Monte Carlo Suite for Modeling Instrumental Line and Continuum Backgrounds in Gamma-Ray Astronomy

Intense and complex instrumental backgrounds, against which the much smaller signals from celestial sources have to be discerned, are a notorious problem for low and intermediate energy gamma-ray astronomy (approximately 50 keV - 10 MeV). Therefore a detailed qualitative and quantitative understanding of instrumental line and continuum backgrounds is crucial for most stages of gamma-ray astronomy missions, ranging from the design and development of new instrumentation through performance prediction to data reduction. We have developed MGGPOD, a user-friendly suite of Monte Carlo codes built around the widely used GEANT (Version 3.21) package, to simulate ab initio the physical processes relevant for the production of instrumental backgrounds. These include the build-up and delayed decay of radioactive isotopes as well as the prompt de-excitation of excited nuclei, both of which give rise to a plethora of instrumental gamma-ray background lines in addition t o continuum backgrounds. The MGGPOD package and documentation are publicly available for download. We demonstrate the capabilities of the MGGPOD suite by modeling high resolution gamma-ray spectra recorded by the Transient Gamma-Ray Spectrometer (TGRS) on board Wind during 1995. The TGRS is a Ge spectrometer operating in the 40 keV to 8 MeV range. Due to its fine energy resolution, these spectra reveal the complex instrumental background in formidable detail, particularly the many prompt and delayed gamma-ray lines. We evaluate the successes and failures of the MGGPOD package in reproducing TGRS data, and provide identifications for the numerous instrumental lines.

Weidenspointner, G.↗
NASA NTRS1996

Orbiter Aerodynamic Acceleration Flight Measurements in the Rarefied-Flow Transition Regime

Acceleration data taken from the Orbital Acceleration Research Experiment (OARE) during reentry on STS-62 have been analyzed using calibration factors taken on orbit. This is the first Orbiter mission which collected OARE data during the Orbiter reentry phase. The data examined include the flight regime from orbital altitudes down to about 90 km which covers the free-molecule-flow regime and the upper altitude fringes of the rarefied-flow transition into the hypersonic continuum. Ancillary flight data on Orbiter position, orientation, velocity, and rotation rates have been used in models to transform the measured accelerations to the Orbiter center-of-gravity, from which aerodynamic accelerations along the Orbiter body axes have been calculated. Residual offsets introduced in the measurements by unmodeled Orbiter forces are identified and discussed. Direct comparisons are made between the OARE flight data and an independent micro-gravity accelerometer experiment, the High Resolution Accelerometer Package (HiRAP), which also obtained flight data on reentry during the mission down to about 95 km. The resulting OARE aerodynamic acceleration measurements along the Orbiter's body axis, aid the normal to axial acceleration ratio in the free-molecule-flow and transition-flow regimes are presented and compared with numerical simulations from three direct simulation Monte Carlo codes.

Blanchard, Robert C.↗
NASA NTRS2004

Density, Velocity and Ionization Structure in Accretion-Disc Winds

This was a project to exploit the unique capabilities of FUSE to monitor variations in the wind- formed spectral lines of the luminous, low-inclination, cataclysmic variables(CV) -- RW Sex. (The original proposal contained two additional objects but these were not approved.) These observations were intended to allow us to determine the relative roles of density and ionization state changes in the outflow and to search for spectroscopic signatures of stochastic small-scale structure and shocked gas. By monitoring the temporal behavior of blue-ward extended absorption lines with a wide range of ionization potentials and excitation energies, we proposed to track the changing physical conditions in the outflow. We planned to use a new Monte Carlo code to calculate the ionization structure of and radiative transfer through the CV wind. The analysis therefore was intended to establish the wind geometry, kinematics and ionization state, both in a time-averaged sense and as a function of time.

Sonneborn, George↗