Search NASA⌕ Search

SEARCH · Search NASA

Results for “The Random Ray Method”

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.

69 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Formation and Evolution of Binary Systems Containing Collapsed Stars

This research includes theoretical studies of the formation and evolution of five types of interacting binary systems. Our main focus has been on developing a number of comprehensive population synthesis codes to study the following types of binary systems: (i) cataclysmic variables (#3, #8, #12, #15), (ii) low- and intermediate-mass X-ray binaries (#13, #20, #21), (iii) high-mass X-ray binaries (#14, #17, #22), (iv) recycled binary millisecond pulsars in globular clusters (#5, #10, #ll), and (v) planetary nebulae which form in interacting binaries (#6, #9). The numbers in parentheses refer to papers published or in preparation that are listed in this paper. These codes take a new unified approach to population synthesis studies. The first step involves a Monte Carlo selection of the primordial binaries, including the constituent masses, and orbital separations and eccentricities. Next, a variety of analytic methods are used to evolve the primary star to the point where either a dynamical episode of mass transfer to the secondary occurs (the common envelope phase), or the system evolves down an alternate path. If the residual core of the primary is greater than 2.5 solar mass, it will evolve to Fe core collapse and the production of a neutron star and a supernova explosion. In the case of systems involving neutron stars, a kick velocity is chosen randomly from an appropriate distribution and added to the orbital dynamics which determine the state of the binary system after the supernova explosion. In the third step, all binaries which commence stable mass transfer from the donor star (the original secondary in the binary system) to the compact object, are followed with a detailed binary evolution code. Finally, we include all the relevant dynamics of the binary system. For example, in the case of LMXBs, the binary system, with its recoil velocity from the supernova explosion, is followed in time through its path in the Galactic potential. For our globular cluster study, detailed three-body calculations are used to follow the interactions of the binary with passing field stars.

Rappaport, Saul↗
NASA NTRS2005

Strong-field Gravitational Effects in the Timing and Spectroscopy of X-Ray Binaries

Using the ray-tracing code described, we developed additional features to the hot spot model to explain the broadening of the quasi-periodic oscillation (QPO) peaks as well as the damping of higher frequency harmonics in the power spectrum. We derived a number of analytic results that agree closely with more detailed numerical calculations. Three primary pieces are developed: the superposition of multiple hot spots with finite lifetimes and random phases, a width in the radial distribution of geodesic orbits, and the scattering of photons from the hot spot through a corona around the black hole. The complete model is used to fit the observed power spectra of both type A and type B QPOs seen in XTE 51550-564, giving confidence limits on each of the model parameters. These results were presented in a conference paper and in a refereed publication. We also explored the possibility of using higher-order statistics, particularly the bispectrum and bicoherence, to distinguish between competing QPO models. These methods were applied to the hot spot model and a method was proposed to map out the spacetime metric around Kerr black holes in Ref. 3 below.

Bertschinger, Edmund↗
NASA NTRS1997

Csa-19, a radiation-responsive human gene, identified by an unbiased two-gel cDNA library screening method in human cancer cells

A novel polymerase chain reaction (PCR)-based method was used to identify candidate genes whose expression is altered in cancer cells by ionizing radiation. Transcriptional induction of randomly selected genes in control versus irradiated human HL60 cells was compared. Among several complementary DNA (cDNA) clones recovered by this approach, one cDNA clone (CL68-5) was downregulated in X-irradiated HL60 cells but unaffected by 12-O-tetradecanoyl phorbol-13-acetate, forskolin, or cyclosporin-A. DNA sequencing of the CL68-5 cDNA revealed 100% nucleotide sequence homology to the reported human Csa-19 gene. Northern blot analysis of RNA from control and irradiated cells revealed the expression of a single 0.7-kilobase (kb) messenger RNA (mRNA) transcript. This 0.7-kb Csa-19 mRNA transcript was also expressed in a variety of human adult and corresponding fetal normal tissues. Moreover, when the effect of X- or fission neutron-irradiation on Csa-19 mRNA was compared in cultured human cells differing in p53 gene status (p53-/- versus p53+/+), downregulation of Csa-19 by X-rays or fission neutrons was similar in p53-wild type and p53-null cell lines. Our results provide the first known example of a radiation-responsive gene in human cancer cells whose expression is not associated with p53, adenylate cyclase or protein kinase C.

NASA Discipline Radiation Health↗
NASA NTRS1999

Lateral and Vertical Heterogeneity of Thorium in the Procellarum KREEP Terrane: As Reflected in the Ejecta Deposits of Post-Imbrium Craters

The Procellarum KREEP Terrane displays the highest concentrations of Th on the Moon. However, locations of elevated Th in this region appear to be random. As observed in the 5 deg per pixel equal-area Th data, and made more evident in the preliminary 2 deg data, Th is enhanced around the craters Aristillus, Aristarchus, Kepler, Mairan, the Apennine Bench formation, and the Fra Mauro region, while noticeably and unexpectedly lower in other locations (e.g., Archimedes, Copernicus, Eratosthenes, and Plato). We have examined the composition of the materials present in these regions with the goal of understanding the patchy nature to the distribution of Th and ultimately to decipher the geologic processes that have concentrated the Th. At present time, the published resolution of the Lunar Prospector Th gamma-ray data is low (5 deg per pixel), but this will soon be superceded by significantly higher-resolution data (2 deg per pixel). Even at this improved resolution, however, it is difficult to resolve the units that are the major source of Th. In an attempt to circumvent this problem, we employ the higher-resolution Clementine multispectral data for those regions mentioned above. We use the UV-VIS-derived compositional information and the spectral properties of craters, and their ejecta as drill holes through the mare-basalt surface to investigate the thickness and composition of underlying material. With this information we attempt to piece together the stratigraphy and geologic history of the Imbrium-Procellanim region. We processed the five-band multispectral data from the Clementine Mission (415, 750,900,950, and 1000nm) using ISIS software and calibration parameters developed by the USGS, Flagstaff, Arizona. Final image mosaics are in equal-area sinusoidal projection, and have a resolution of 250 m/pixel. Using the method of we produced maps of FeO and Ti02 composition. Here we examine the Th, FeO, and Ti02 composition and spectral properties of the craters discussed above and their ejecta, with the goal of describing the materials they excavate. One interpretation for the origin of the high-Th material is that subsurface KREEPy materials have been excavated by impact craters. The material excavated may be either volcanic KREEP (e,g., Apennine Bench Formation), KREEPy impact-melt breccia formed by the Imbrium impact (e.g., Fra Mauro Formation), or other KREEP-rich crustal material. Determining which type of material is responsible for the elevated Th and its extent is important to understanding the premare and possibly the prebasin stratigraphy of the Imbrium-Procellarum Region. Merging the 5 deg. Th data with the shaded relief map, we observe that the highest Th concentrations are not related to pre-Imbrium upper crustal materials. The Apennines, Alpes, and Caucasus Mountains represent the pre-Imbrian highlands material and do not express concentrations of Th, FeO, and TiO2 as high as the most Th-fich materials exposed within the Procellarum KREEP Terrane. We observe that, in general, these massifs contain 10-14 wt% FeO and 4-7 ppm Th. Determining whether the Th signal is from KREEP basalts or KREEPy impact-melt breccias cannot be done with the Clementine data because the two rock types are compositionally and mineralogically too similar (e.g., the Th-rich, mafic impact-melt breccias in the Apollo sample collection are dominated by a KREEP-basalt like component. Mapping-the distribution and sizes of craters and whether they display elevated Th concentrations or not, should reveal the depth and thickness of the KREEP-rich materials, and whether they are ubiquitous (i.e., impact-melt breccia) or more randomly distributed; this might be taken as an indicator of localized KREEP-basalt flows. Within the southeastern region of the Imbrium basin, there are two Th hot spots. The first is associated with the crater Aristillus, and the latter with the Apennine Bench Formation. Adjacent to these two hot spots are craters with a lower Th signature: Archimedes and Autolycus. We observe in the ejecta of Aristillus, a region of significantly lower FeO (10-14 wt%) relative to the surrounding mare basalt. The crater Autolycus, 50 km to the south, did not excavate similar low-FeO material. We suggest that the lower-FeO material in the ejecta of Aristillus corresponds to Th-rich material; the FeO content observed in Aristillus ejecta is comparable to that of KREEP basalt or mafic impact melt breccia (10-12 wt% FeO). We determine that this low FeO, Th-rich material is volcanic KREEP, as opposed to Imbrium impact melt, on the basis that the low-Fe material is exposed more prominently in ejecta in the northern portion of Aristillus. Our assumption is that if the layer underlying Aristillus was continuous, a more widespread and uniform low-Fe signature wouldbe observed in the ejecta deposit. Archimedes, 110 km southwest of Aristillus, impacted the northern portion of the Apennine Bench prior to the eruption of KREEP basalt. Archimedes rim material is not as enriched in Th as the Apennine Bench, and there are differences between the two in FeO concentration and in their continuum slope. Archimedes exhibits a much steeperor "redder" continuum slope than the Apennine Bench. This steepslope suggests the presence of glassy material. The glassy material is concentrated around an unnamed crater on the southern rim of Archimedes (4.5W, 28.2N) and along the northern rim of Archimedes. We suggest two possibilities, or a combination of the two, to explain the low-Th signal from Archimedes: (1) The Apennine Bench prior to KREEP basalt eruption was lower in Th (4-7 ppm, e.g., similar to the Apennine massifs) and KREEP basalts are absent in the rim of Archimedes; or (2) the glassy (possibly pyroclastic) material layering the rim of the Archimedes, dilutes any high-Th material present with low-Th material. (Additional information is contained in the original)

Gillis, J. J.↗
NASA NTRS2004

Brazing SiC/SiC Composites to Metals

Experiments have shown that active brazing alloys (ABAs) can be used to join SiC/SiC composite materials to metals, with bond strengths sufficient for some structural applications. The SiC/SiC composite coupons used in the experiments were made from polymerbased SiC fiber preforms that were chemical-vapor-infiltrated with SiC to form SiC matrices. Some of the metal coupons used in the experiments were made from 304 stainless steel; others were made from oxygen-free, high-conductivity copper. Three ABAs were chosen for the experiments: two were chosen randomly from among a number of ABAs that were on hand at the time; the third ABA was chosen because its titanium content (1.25 percent) is less than those of the other two ABAs (1.75 and 4.5 percent, respectively) and it was desired to evaluate the effect of reducing the titanium content, as described below. The characteristics of ABAs that are considered to be beneficial for the purpose of joining SiC/SiC to metal include wettability, reactivity, and adhesion to SiC-based ceramics. Prior to further development, it was verified that the three chosen ABAs have these characteristics. For each ABA, suitable vacuum brazing process conditions were established empirically by producing a series of (SiC/SiC)/ABA wetting samples. These samples were then sectioned and subjected to scanning electron microscopy (SEM) and energy-dispersive x-ray spectrometry (EDS) for analysis of their microstructures and compositions. Specimens for destructive mechanical tests were fabricated by brazing of lap joints between SiC/SiC coupons 1/8-in. (.3.2- mm) thick and, variously, stainless steel or copper tabs. The results of destructive mechanical tests and the SEM/EDS analysis were used to guide the development of a viable method of brazing the affected materials.

Steffier, Wayne S.↗
NASA NTRS1987

Computation of asymmetric supersonic flows around cones at large incidence

The Schiff-Steger parabolized Navier-Stokes (PNS) code has been modified to allow computation of conical flowfields around cones at high incidence. The improved algorithm of Degani and Schiff has been incorporated with the PNS code. This algorithm adds the cross derivative and circumferential viscous terms to the original PNS code and modifies the algebraic eddy viscosity turbulence model to take into account regions of so called cross-flow separation. Assuming the flowfield is conical (but not necessarily symmetric) a marching stepback procedure is used: the solution is marched one step downstream using improved PNS code and the flow variables are then scaled to place the solution back to the original station. The process is repeated until no change in the flow variables is observed with further marching. The flow variables are then constant along rays of the flowfield. The experiments obtained by Bannik and Nebbeling were chosen as a test case. In these experiments a cone of 7.5 deg. half angle at Mach number 2.94 and Reynolds number 1.372 x 10(7) was tested up 34 deg. angle of attack. At high angle of attack nonconical asymmetric leeward side vortex patterns were observed. In the first set of computations, using an earlier obtained solution of the above cone for angle of attack of 22.6 deg. and at station x=0.5 as a starting solution, the angle of attack was gradually increased up to 34 deg. During this procedure the grid was carfully adjusted to capture the bow shock. A stable, converged symmetric solution was obtained. Since the numerical code converged to a symmetric solution which is not the physical one, the stability was tested by a random perturbation at each point. The possible effect of surface roughness or non perfect body shape was also investigated. It was concluded that although the assumption of conical viscous flows can be very useful for certain cases, it can not be used for the present case. Thus the second part of the investigation attempted to obtain a marching (in space) solution with the PNS method using the conical solution as initial data. Finally, the solution of the full Navier-Stokes equations was carried out.

Degani, David↗
NASA NTRS1995

An experimental study of light scattering by large, irregular particles

The intensity and polarization of light scattered by a variety of types of artificial partices large compared to the wavelength were measured as a function of phase angle. Shape, surface roughness, absorption coefficient, and internal scattering coefficient were varied systematically and their effects studied. Scattering by clear, smooth-surfaced spheres is in quantitative agreement with the predictions of the geometrical optics (ray theory) approximation to physical optics (Mie theory). The phase functions of almost all of the particles measured have both forward and backward scattering lobes. A two-parameter, double Henyey-Greenstein function generally provides reasonably good descriptions of the data, while keeping the number of free parameters to the minimum necessary. On a double Henyey- Greenstein parameter plot all of the particles fall into an L-shaped area of restricted size in which the location is characteristic of the particle type. Formalisms based on the equivalent slab model are also given for estimating the scattering efficiency of a large, irregular particle. For most dielectric particles the transmitted, forward scattered light is partially negatively polarized. It is this component that is respopnsible for the well-known maximum in the polarization curves of planetary regoliths at phase angles around 100 deg. For phase angles between about 30 deg and 70 deg the internally scattered light is found to be randomly polarized in the particles studied here, so that the only contribution to the second component of the Stokes vector is by Fresnel reflection from the particle surface. If this empirical result is general, measurement of the second Stokes vector of the light scattered from a regolith at these angles may provide a method of remotely measuring the mean refractive index.

Mcguire, Audrey F.↗
NASA NTRS1996

Defect Structure of Beta NiAl Using the BFS Method for Alloys

The semiempirical BFS method for alloys is generalized by replacing experimental input with first-principles results thus allowing for the study of complex systems. In order to examine trends and behavior of a system in the vicinity of a given point of the phase diagram a search procedure based on a sampling of selected configurations is employed. This new approach is applied to the study of the beta phase of the Ni-Al system, which exists over a range of composition from 45-60 at.% Ni. This methodology results in a straightforward and economical way of reproducing and understanding the basic features of this system. At the stoichiometric composition, NiAl should exist in a perfectly ordered B2 structure. Ni-rich alloys are characterized by antisite point defects (with Ni atoms in the Al sites) with a decrease in lattice parameters. On the Al-rich side of stoichiometry there is a steep decrease in lattice parameter and density with increasing Al content. The presence of vacancies in Ni sites would explain such behavior. Recent X-ray diffraction experiments suggest a richer structure: the evidence, while strongly favoring the presence of vacancies in Ni sites, also suggests the possibility of some vacancies in Al sites in a 3:1 ratio. Moreover, local ordering of vacant sites may be preferred over a random distribution of individual point defects.

Bozzolo, Guillermo↗
NASA NTRS2002

Large Coded Aperture Mask for Spaceflight Hard X-ray Images

The 2.6 square meter coded aperture mask is a vital part of the Burst Alert Telescope on the Swift mission. A random, but known pattern of more than 50,000 lead tiles, each 5 mm square, was bonded to a large honeycomb panel which projects a shadow on the detector array during a gamma ray burst. A two-year development process was necessary to explore ideas, apply techniques, and finalize procedures to meet the strict requirements for the coded aperture mask. Challenges included finding a honeycomb substrate with minimal gamma ray attenuation, selecting an adhesive with adequate bond strength to hold the tiles in place but soft enough to allow the tiles to expand and contract without distorting the panel under large temperature gradients, and eliminating excess adhesive from all untiled areas. The largest challenge was to find an efficient way to bond the > 50,000 lead tiles to the panel with positional tolerances measured in microns. In order to generate the desired bondline, adhesive was applied and allowed to cure to each tile. The pre-cured tiles were located in a tool to maintain positional accuracy, wet adhesive was applied to the panel, and it was lowered to the tile surface with synchronized actuators. Using this procedure, the entire tile pattern was transferred to the large honeycomb panel in a single bond. The pressure for the bond was achieved by enclosing the entire system in a vacuum bag. Thermal vacuum and acoustic tests validated this approach. This paper discusses the methods, materials, and techniques used to fabricate this very large and unique coded aperture mask for the Swift mission.

Vigneau, Danielle N.↗
NASA NTRS1998

Repair of clustered DNA damage caused by high LET radiation in human fibroblasts

It has recently been demonstrated experimentally that DNA damage induced by high LET radiation in mammalian cells is non-randomly distributed along the DNA molecule in the form of clusters of various sizes. The sizes of such clusters range from a few base-pairs to at least 200 kilobase-pairs. The high biological efficiency of high LET radiation for induction of relevant biological endpoints is probably a consequence of this clustering, although the exact mechanisms by which the clustering affects the biological outcome is not known. We discuss here results for induction and repair of base damage, single-strand breaks and double-strand breaks for low and high LET radiations. These results are discussed in the context of clustering. Of particular interest is to determine how clustering at different scales affects overall rejoining and fidelity of rejoining of DNA double-strand breaks. However, existing methods for measuring repair of DNA strand breaks are unable to resolve breaks that are close together in a cluster. This causes problems in interpretation of current results from high LET radiation and will require new methods to be developed.

NASA Discipline Radiation Health↗
NASA NTRS1994

Hot accretion disks with pairs: Effects of magnetic field and thermal cyclocsynchrotron radiation

We show the effects of thermal cyclosynchrotron radiation and magnetic viscosity on the structure of hot, two-temperature accretion disks. Magnetic field, B, is assumed to be randomly oriented and the ratio of magnetic pressure to either gas pressure, alpha = P(sub mag)/P(sub gas), or the sum of the gas and radiation pressures, alpha = (P(sub mag)/P(sub gas) + P(sub rad)), is fixed. We find those effects do not change the qualitative properties of the disks, i.e., there are still two critical accretion rates related to production of e(sup +/-) pairs, (M dot)((sup U)(sub cr)) and (M dot)((sup L)(sub cr)), that affect the number of local and global disk solutions, as recently found by Bjoernsson and Svensson for the case with B = 0. However, a critical value of the alpha-viscosity parameter above which those critical accretion rates disappear becomes smaller than alpha(sub cr) = 1 found in the case of B = 0, for P(sub mag) = alpha(P(sub gas) + P(sub rad)). If P(sub mag) = alpha P(sub gas), on the other hand, alpha(sub cr) is still about unity. Moreover, when Comptonized cyclosynchrotron radiation dominates Comptonized bremsstrahlung, radiation from the disk obeys a power law with the energy spectral index of approximately 0.5, in a qualitative agreement with X-ray observations of active galactic nuclei (AGNS) and Galactic black hole candidates. We also extend the hot disk solutions for P(sub mag) = alpha(P(sub gas) + P(sub rad)) to the effectively optically thick region, where they merge with the standard cold disk solutions. We find that the mapping method by Bjoernsson and Svensson gives a good approximation to the disk structure in the hot region and show where it breaks in the transition region. Finally, we find a region in the disk parameter space with no solutions due to the inability of Coulomb heating to supply enough energy to electrons.

Kusunose, Masaaki↗
NASA NTRS2015

Modeling the Swift BAT Trigger Algorithm with Machine Learning

To draw inferences about gamma-ray burst (GRB) source populations based on Swift observations, it is essential to understand the detection efficiency of the Swift burst alert telescope (BAT). This study considers the problem of modeling the Swift BAT triggering algorithm for long GRBs, a computationally expensive procedure, and models it using machine learning algorithms. A large sample of simulated GRBs from Lien et al. (2014) is used to train various models: random forests, boosted decision trees (with AdaBoost), support vector machines, and artificial neural networks. The best models have accuracies of approximately greater than 97% (approximately less than 3% error), which is a significant improvement on a cut in GRB flux which has an accuracy of 89:6% (10:4% error). These models are then used to measure the detection efficiency of Swift as a function of redshift z, which is used to perform Bayesian parameter estimation on the GRB rate distribution. We find a local GRB rate density of eta(sub 0) approximately 0.48(+0.41/-0.23) Gpc(exp -3) yr(exp -1) with power-law indices of eta(sub 1) approximately 1.7(+0.6/-0.5) and eta(sub 2) approximately -5.9(+5.7/-0.1) for GRBs above and below a break point of z(sub 1) approximately 6.8(+2.8/-3.2). This methodology is able to improve upon earlier studies by more accurately modeling Swift detection and using this for fully Bayesian model fitting. The code used in this is analysis is publicly available online.

gamma rays: general↗
NASA NTRS2016

Modeling the Swift Bat Trigger Algorithm with Machine Learning

To draw inferences about gamma-ray burst (GRB) source populations based on Swift observations, it is essential to understand the detection efficiency of the Swift burst alert telescope (BAT). This study considers the problem of modeling the Swift / BAT triggering algorithm for long GRBs, a computationally expensive procedure, and models it using machine learning algorithms. A large sample of simulated GRBs from Lien et al. is used to train various models: random forests, boosted decision trees (with AdaBoost), support vector machines, and artificial neural networks. The best models have accuracies of greater than or equal to 97 percent (less than or equal to 3 percent error), which is a significant improvement on a cut in GRB flux, which has an accuracy of 89.6 percent (10.4 percent error). These models are then used to measure the detection efficiency of Swift as a function of redshift z, which is used to perform Bayesian parameter estimation on the GRB rate distribution. We find a local GRB rate density of n (sub 0) approaching 0.48 (sup plus 0.41) (sub minus 0.23) per cubic gigaparsecs per year with power-law indices of n (sub 1) approaching 1.7 (sup plus 0.6) (sub minus 0.5) and n (sub 2) approaching minus 5.9 (sup plus 5.7) (sub minus 0.1) for GRBs above and below a break point of z (redshift) (sub 1) approaching 6.8 (sup plus 2.8) (sub minus 3.2). This methodology is able to improve upon earlier studies by more accurately modeling Swift detection and using this for fully Bayesian model fitting.

gamma-ray burst: general – gamma-rays: general â↗
NASA NTRS1999

Bidirectional Reflectance of a Macroscopically Flat, High-Albedo Particulate Surface: An Efficient Radiative Transfer Solution and Applications to Regoliths

Many remote sensing applications rely on accurate knowledge of the bidirectional reflection function (BRF) of surfaces composed of discrete, randomly positioned scattering particles. Theoretical computations of BRFs for plane-parallel particulate layers are usually reduced to solving the radiative transfer equation (RTE) using one of existing exact or approximate techniques. Since semi-empirical approximate approaches are notorious for their low accuracy, violation of the energy conservation law, and ability to produce unphysical results, the use of numerically exact solutions of RTE has gained justified popularity. For example, the computation of BRFs for macroscopically flat particulate surfaces in many geophysical publications is based on the adding-doubling (AD) and discrete ordinate (DO) methods. A further saving of computer resources can be achieved by using a more efficient technique to solve the plane-parallel RTE than the AD and DO methods. Since many natural particulate surfaces can be well represented by the model of an optically semi-infinite, homogeneous scattering layer, one can find the BRF directly by solving the Ambartsumian's nonlinear integral equation using a simple iterative technique. In this way, the computation of the internal radiation field is avoided and the computer code becomes highly efficient and very accurate and compact. Furthermore, the BRF thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. In this paper, we discuss numerical aspects and the computer implementation of this technique, examine the applicability of the Henyey-Greenstein phase function and the sigma-Eddington approximation in BRF and flux calculations, and describe sample applications demonstrating the potential effect of particle shape on the bidirectional reflectance of flat regolith surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web, and can be applied to a wide range of remote sensing problems. BRF computations for undulated (macroscopically rough) surfaces are more complicated and often rely on time consuming Monte Carlo procedures. This approach is especially inefficient for optically thick, weakly absorbing media (e.g., snow and desert surfaces at visible wavelengths since a photon may undergo many internal scattering events before it exists the medium or is absorbed. However, undulated surfaces can often be represented as collections of locally flat tilted facets characterized by the BRF found from the traditional plane parallel RTE. In this way the MOnte Carlo procedure could be used only to evaluate the effects of surface shadowing and multiple surface reflections, thereby bypassing the time-consuming ray tracing inside the medium and providing a great savings of CPU time.

Mishchenko, Michael I.↗
NASA NTRS2000

Sputtering Erosion in the Ion Thruster

During the first phase of this research, the sputtering yields of molybdenum by low energy (100 eV and higher) xenon ions were measured by using the methods of secondary neutral mass spectrometry (SNMS) and Rutherford backscattering spectrometry (RBS). However, the measured sputtering yields were found to be far too low to explain the sputtering erosions observed in the long-duration tests of ion thrusters. The only difference between the sputtering yield measurement experiments and the ion thruster tests was that the later are conducted at high ion fluences. Hence, a study was initiated to investigate if any linkage exists between high ion fluence and an enhanced sputtering yield. The objective of this research is to gain an understanding of the causes of the discrepancies between the sputtering rates of molybdenum grids in an ion thruster and those measured from our experiments. We are developing a molecular dynamics simulation technique for studying low-energy xenon ion interactions with molybdenum. It is difficult to determine collision sequences analytically for primary ions below the 200 eV energy range where the ion energy is too low to be able to employ a random cascade model with confidence and it is too high to have to consider only single collision at or near the surface. At these low energies, the range of primary ions is about 1 to 2 nm from the surface and it takes less than 4 collisions on the average to get an ion to degrade to such an energy that it can no longer migrate. The fine details of atomic motion during the sputtering process are revealed through computer simulation schemes. By using an appropriate interatomic potential, the positions and velocities of the incident ion together with a sufficient number of target atoms are determined in small time steps. Hence, it allows one to study the evolution of damages in the target and its effect on the sputtering yield. We are at the preliminary stages of setting up the simulation program.

Ray, Pradosh K.↗