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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
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.↗
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↗