Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Techniques”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Production of recycled pulsars in globular clusters via two-body tidal capture

The numbers of binary systems containing a neutron star that are formed in the cores of globular clusters via two-body tidal capture are investigated. The subsequent evolution of the resulting binary systems is examined. Monte Carlo techniques are employed to trace the histories of many neutron stars in the core of a model globular cluster. Various scenarios for mass stripping during capture and orbital circularization, and the effects of magnetic braking on the circularized orbits are explored. Calculations are carried out for two model clusters: 47 Tuc and Omega Cen. Distributions of orbital periods for newly formed binaries containing a neutron star in orbit with either a main-sequence star or a giant are presented. Each binary is followed until contact between the captured star and the neutron star is established. It is argued that these systems should be regarded as potential progenitors of 'recycled' pulsars.

Di Stefano, R↗
NASA NTRS1993

Numerical simulation of rarefied flow through a slit. I - Direct simulation Monte Carlo results

The pressure-driven flow of a rarefied monatomic gas through a two-dimensional slit is simulated using the direct simulation Monte Carlo technique. Of particular interest is the change in flow field structure as pressure ratio and Knudsen number are varied. Comparisons are made to quantify the limits of validity of free-molecular theory and approximate, nearly free-molecular iterative methods. Also addressed is the sensitivity of the numerical solutions to grid structure and boundary conditions. The free-molecular theory is found to predict quantitative flow field properties (e.g., centerline velocities or downstream flux profiles) reasonably well for large finite Knudsen number with the error dependent on the pressure ratio. The nearly free-molecular corrections are shown to have limited range of applicability. A previously derived parameter is found to correlate total mass flux well as a function of pressure ratio and Knudsen number over a large portion of the transitional regime.

Wadsworth, D. C.↗
NASA NTRS1990

Statistical circuit design for yield improvement in CMOS circuits

This paper addresses the statistical design of CMOS integrated circuits for improved parametric yield. The work uses the Monte Carlo technique of circuit simulation to obtain an unbiased estimation of the yield. A simple graphical analysis tool, the yield factor histogram, is presented. The yield factor histograms are generated by a new computer program called SPICENTER. Using the yield factor histograms, the most sensitive circuit parameters are noted, and their nominal values are changed to improve the yield. Two basic CMOS example circuits, one analog and one digital, are chosen and their designs are 'centered' to illustrate the use of the yield factor histograms for statistical circuit design.

Kamath, H. J.↗
NASA NTRS1992

Acceleration of positrons in supernova shocks

During this project we investigated the acceleration of leptons (electrons and positrons) in collisionless shock waves. In particular, we were interested in how leptons are accelerated in the blast waves existing in the remnants of supernova explosions. Supernova remnants (SNRs) have long been considered as the most likely source of galactic cosmic rays but no definite connection between SNRs and the cosmic rays seen at earth can be made. Only by understanding lepton acceleration in shocks can the rich SNR data base be properly used to understand cosmic ray origins. Our project was directed at the neglected aspects of lepton acceleration. We showed that the efficiency of lepton acceleration depended critically on the lepton injection energy. We showed that, even when infection effects are not important, that proton and lepton distribution functions produced by shocks are quite different in the critical energy range for producing the observed synchrotron emission. We also showed that transrelativistic effects produced proton spectra that were not in agreement with standard results from radio observations, but that the lepton spectra were, in fact, consistent with observations. We performed simulations of relativistic shocks (shocks where the flow speed is a sizable fraction of the speed of light) and discovered some interesting effects. We first demonstrated the power of the Monte Carlo technique by determining the shock jump conditions in relativistic shocks. We then proceeded to determine how relativistic shocks accelerate particles. We found that nonlinear relativistic shocks treat protons and leptons even more differently than nonrelativistic shocks. The transrelativistic effects on the shock structure from the heavy ion component reduces the lepton acceleration to a tiny fraction of the ion acceleration. This effect is dramatic even if high energy leptons (many times thermal energy) are injected, and was totally unexpected. Our results have important consequences for astrophysical environments expected to harbor relativistic flows such as extra-galactic radio sources and accretion onto compact objects.

Ellison, Donald C.↗
NASA NTRS1991

Low altitude plume impingement handbook

Plume Impingement modeling is required whenever an object immersed in a rocket exhaust plume must survive or remain undamaged within specified limits, due to thermal and pressure environments induced by the plume. At high altitudes inviscid plume models, Monte Carlo techniques along with the Plume Impingement Program can be used to predict reasonably accurate environments since there are usually no strong flowfield/body interactions or atmospheric effects. However, at low altitudes there is plume-atmospheric mixing and potential large flowfield perturbations due to plume-structure interaction. If the impinged surface is large relative to the flowfield and the flowfield is supersonic, the shock near the surface can stand off the surface several exit radii. This results in an effective total pressure that is higher than that which exists in the free plume at the surface. Additionally, in two phase plumes, there can be strong particle-gas interaction in the flowfield immediately ahead of the surface. To date there have been three levels of sophistication that have been used for low altitude plume induced environment predictions. Level 1 calculations rely on empirical characterizations of the flowfield and relatively simple impingement modeling. An example of this technique is described by Piesik. A Level 2 approach consists of characterizing the viscous plume using the SPF/2 code or RAMP2/LAMP and using the Plume Impingement Program to predict the environments. A Level 3 analysis would consist of using a Navier-Stokes code such as the FDNS code to model the flowfield and structure during a single calculation. To date, Level 1 and Level 2 type analyses have been primarily used to perform environment calculations. The recent advances in CFD modeling and computer resources allow Level 2 type analysis to be used for final design studies. Following some background on low altitude impingement, Level 1, 2, and 3 type analysis will be described.

Smith, Sheldon D.↗
NASA NTRS1993

Thermal analysis of combinatorial solid geometry models using SINDA

Algorithms have been developed using Monte Carlo techniques to determine the thermal network parameters necessary to perform a finite difference analysis on Combinatorial Solid Geometry (CSG) models. Orbital and laser fluxes as well as internal heat generation are modeled to facilitate satellite modeling. The results of the thermal calculations are used to model the infrared (IR) images of targets and assess target vulnerability. Sample analyses and validation are presented which demonstrate code products.

Gerencser, Diane↗
NASA NTRS1994

Comments on the statistical analysis of excess variance in the COBE differential microwave radiometer maps

Cosmic anisotrophy produces an excess variance sq sigma(sub sky) in the Delta maps produced by the Differential Microwave Radiometer (DMR) on cosmic background explorer (COBE) that is over and above the instrument noise. After smoothing to an effective resolution of 10 deg, this excess sigma(sub sky)(10 deg), provides an estimate for the amplitude of the primordial density perturbation power spectrum with a cosmic uncertainty of only 12%. We employ detailed Monte Carlo techniques to express the amplitude derived from this statistic in terms of the universal root mean square (rms) quadrupole amplitude, (Q sq/RMS)(exp 0.5). The effects of monopole and dipole subtraction and the non-Gaussian shape of the DMR beam cause the derived (Q sq/RMS)(exp 0.5) to be 5%-10% larger than would be derived using simplified analytic approximations. We also investigate the properties of two other map statistics: the actual quadrupole and the Boughn-Cottingham statistic. Both the sigma(sub sky)(10 deg) statistic and the Boughn-Cottingham statistic are consistent with the (Q sq/RMS)(exp 0.5) = 17 +/- 5 micro K reported by Smoot et al. (1992) and Wright et al. (1992).

Wright, E. L.↗
NASA NTRS1994

Cross-correlation of point series using a new method

Traditional methods of cross-correlation of two time series do not apply to point time series. Here, a new method, devised specifically for point series, utilizes a correlation measure that is based in the rms difference (or, alternatively, the median absolute difference) between nearest neightbors in overlapped segments of the two series. Error estimates for the observed locations of the points, as well as a systematic shift of one series with respect to the other to accommodate a constant, but unknown, lead or lag, are easily incorporated into the analysis using Monte Carlo techniques. A methodological restriction adopted here is that one series be treated as a template series against which the other, called the target series, is cross-correlated. To estimate a significance level for the correlation measure, the adopted alternative (null) hypothesis is that the target series arises from a homogeneous Poisson process. The new method is applied to cross-correlating the times of the greatest geomagnetic storms with the times of maximum in the undecennial solar activity cycle.

Strothers, Richard B.↗
NASA NTRS1994

The far-ultraviolet dust albedo in the Upper Scorpius subgroup of the Scorpius OB2 association

During NRL's Far Ultraviolet Cameras experiemnt on STS-39, four images of the giant reflection nebula encompassing the Upper Scorpius subgroup of the Sco OB2 association were obtained in two ultraviolet bandpasses with lambda(sub eff) = 1362 A and 1769 A. From these images and IUE and TD-1 stellar spectra, the ratio of nebular to stellar flux was calculated. The ratio ranged from 0.577 to 0.921 at 1362 A and 0.681 to 0.916 at 1769 A with the spread in the ratio arising mainly from uncertainties in the sky background. In order to analyze these images, a model utilizing Monte Carlo techniques to describe radiative transfer in a spherical nebula with asymmetrically distributed stars was developed by elaborating on previous work by Witt. This model was used to determine the range of albedos reproducing the observed nebular-to-stellar flux ratios while allowing the scattering phase function asymmetry to vary between 0.0 and 0.8. The resulting albedos were 0.47-0.70 at 1362 A and 0.55-0.72 at 1769 A.

Gordon, Karl D.↗
NASA NTRS1994

A Coulomb collision algorithm for weighted particle simulations

A binary Coulomb collision algorithm is developed for weighted particle simulations employing Monte Carlo techniques. Charged particles within a given spatial grid cell are pair-wise scattered, explicitly conserving momentum and implicitly conserving energy. A similar algorithm developed by Takizuka and Abe (1977) conserves momentum and energy provided the particles are unweighted (each particle representing equal fractions of the total particle density). If applied as is to simulations incorporating weighted particles, the plasma temperatures equilibrate to an incorrect temperature, as compared to theory. Using the appropriate pairing statistics, a Coulomb collision algorithm is developed for weighted particles. The algorithm conserves energy and momentum and produces the appropriate relaxation time scales as compared to theoretical predictions. Such an algorithm is necessary for future work studying self-consistent multi-species kinetic transport.

Miller, Ronald H.↗
NASA NTRS1995

Monte Carlo simulations of particle acceleration at oblique shocks: Including cross-field diffusion

The Monte Carlo technique of simulating diffusive particle acceleration at shocks has made spectral predictions that compare extremely well with particle distributions observed at the quasi-parallel region of the earth's bow shock. The current extension of this work to compare simulation predictions with particle spectra at oblique interplanetary shocks has required the inclusion of significant cross-field diffusion (strong scattering) in the simulation technique, since oblique shocks are intrinsically inefficient in the limit of weak scattering. In this paper, we present results from the method we have developed for the inclusion of cross-field diffusion in our simulations, namely model predictions of particle spectra downstream of oblique subluminal shocks. While the high-energy spectral index is independent of the shock obliquity and the strength of the scattering, the latter is observed to profoundly influence the efficiency of injection of cosmic rays into the acceleration process.

Baring, M. G.↗
NASA NTRS1997

Multi-Species Test of Ion Cyclotron Resonance Heating at High Altitudes

Observations of ion distributions and plasma waves obtained by the Dynamics Explorer 1 satellite in the high-altitude, nightside auroral zone are used to study ion energization for three ion species. A number of theoretical models have been proposed to account for the transverse heating of these ion populations. One of these, the ion cyclotron resonance heating (ICRH) mechanism, explains ion conic formation through ion cyclotron resonance with broadband electromagnetic wave turbulence in the vicinity of the characteristic ion cyclotron frequency. The cyclotron resonant heating of the ions by low- frequency electromagnetic waves is an important energy source for the transport of ions from the ionosphere to the magnetosphere. In this paper we test the applicability of the ICRH mechanism to three simultaneously heated and accelerated ion species by modelling the ion conic formation in terms of a resonant wave-particle interaction in which the ions extract energy from the portion of the broadband electromagnetic wave spectrum which includes the ion cyclotron frequency. Using a Monte Carlo technique we evaluate the ion heating produced by the electromagnetic turbulence at low frequencies and find that the wave amplitudes near the ion cyclotron frequencies are sufficient to explain the observed ion energies.

Persoon, A. M.↗
NASA NTRS1997

Multi-Species Test of Ion Cyclotron Resonance Heating at High Altitudes

Observations of ion distributions and plasma waves obtained by the Dynamics Explorer 1 satellite in the high-altitude, nightside auroral zone are used to study ion energization for three ion species. A number of theoretical models have been proposed to account for the transverse heating of these ion populations. One of these, the ion cyclotron resonance heating (ICRH) mechanism, explains ion conic formation through ion cyclotron resonance with broadband electromagnetic wave turbulence in the vicinity of the characteristic ion cyclotron frequency. The cyclotron resonant heating of the ions by low-frequency electromagnetic waves is an important energy source for the transport of ions from the ionosphere to the magnetosphere. In this paper we test the applicability of the ICRH mechanism to three simultaneously heated and accelerated ion species by modelling the ion conic formation in terms of a resonant wave-particle interaction in which the ions extract energy from the portion of the broadband electromagnetic wave spectrum which includes the ion cyclotron frequency. Using a Monte Carlo technique we evaluate the ion heating produced by the electromagnetic turbulence at low frequencies and find that the wave amplitudes near the ion cyclotron frequencies are sufficient to explain the observed ion energies.

Persoon, A. M.↗
NASA NTRS1998

Monte Carlo Simulation of Endlinking Oligomers

This report describes initial efforts to model the endlinking reaction of phenylethynyl-terminated oligomers. Several different molecular weights were simulated using the Bond Fluctuation Monte Carlo technique on a 20 x 20 x 20 unit lattice with periodic boundary conditions. After a monodisperse "melt" was equilibrated, chain ends were linked whenever they came within the allowed bond distance. Ends remained reactive throughout, so that multiple links were permitted. Even under these very liberal crosslinking assumptions, geometrical factors limited the degree of crosslinking. Average crosslink functionalities were 2.3 to 2.6; surprisingly, they did not depend strongly on the chain length. These results agreed well with the degrees of crosslinking inferred from experiment in a cured phenylethynyl-terminated polyimide oligomer.

Hinkley, Jeffrey A.↗
NASA NTRS1960

An Exploratory Statistical Analysis of a Planet Approach-Phase Guidance Scheme Using Angular Measurements with Significant Error

An exploratory analysis of vehicle guidance during the approach to a target planet is presented. The objective of the guidance maneuver is to guide the vehicle to a specific perigee distance with a high degree of accuracy and minimum corrective velocity expenditure. The guidance maneuver is simulated by considering the random sampling of real measurements with significant error and reducing this information to prescribe appropriate corrective action. The instrumentation system assumed includes optical and/or infrared devices to indicate range and a reference angle in the trajectory plane. Statistical results are obtained by Monte-Carlo techniques and are shown as the expectation of guidance accuracy and velocity-increment requirements. Results are nondimensional and applicable to any planet within limits of two-body assumptions. The problem of determining how many corrections to make and when to make them is a consequence of the conflicting requirement of accurate trajectory determination and propulsion. Optimum values were found for a vehicle approaching a planet along a parabolic trajectory with an initial perigee distance of 5 radii and a target perigee of 1.02 radii. In this example measurement errors were less than i minute of arc. Results indicate that four corrections applied in the vicinity of 50, 16, 15, and 1.5 radii, respectively, yield minimum velocity-increment requirements. Thrust devices capable of producing a large variation of velocity-increment size are required. For a vehicle approaching the earth, miss distances within 32 miles are obtained with 90-percent probability. Total velocity increments used in guidance are less than 3300 feet per second with 90-percent probability. It is noted that the above representative results are valid only for the particular guidance scheme hypothesized in this analysis. A parametric study is presented which indicates the effects of measurement error size, initial perigee, and initial energy on the guidance requirements. Measurement error size significantly affects both guidance accuracy and velocity-increment expenditure. The initial trajectory, as given by its perigee and energy, affects the velocity-increment expenditure but not final guidance accuracy.

Friedlander, Alan L.↗
NASA NTRS1998

Direct Acceleration of Pickup Ions at The Solar Wind Termination Shock: The Production of Anomalous Cosmic Rays

We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed, which in effect solves the Boltzmann equation and is not restricted to near-isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross- field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that modestly strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He(+)/H(+) or O(+)/H(+) to within a factor approx. 5. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft. The fact that ACR fluxes can be modeled with standard shock assumptions suggests that the much-discussed "injection problem" for highly oblique shocks stems from incomplete (either mathematical or computer) modeling of these shocks rather than from any actual difficulty shocks may have in injecting and accelerating thermal or quasi-thermal particles.

Ellison, Donald C.↗
NASA NTRS1998

An Analysis of an Improved Bus-Based Multiprocessor Architecture

This paper analyses the effectiveness of a hybrid multiprocessing/multicomputing architecture that is based upon a single-board-computer multiprocessor (SBCM) architecture. Based upon empirical analysis using discrete event simulations and Monte Carlo techniques, this hybrid architecture, called the enhanced single-board-computer multiprocessor (ESBCM), is shown to have improved performance and scalability characteristics over current SBCM designs.

Ricks, Kenneth G.↗
NASA NTRS1999

Simulation of Dynamics of PVT Growth: ZnSe

An unsteady model for the simulation of PVT growth process has been developed and applied to growth of ZnSe. The model is capable of tracking the unsteady aspects of the growth process caused by, for example, changes in the system temperature associated with growth, source-depletion, translation of the charge, out-gassing at high temperatures, as well as changes in the source composition during growth. The model predicts the evolution of growth interface morphology as well as the shape of source-gas interface. Thermo-fluid transport in the gas phase is handled by a Low-Mach number unsteady compressible flow formulation. Heat transfer in the crystal, source, and the quartz housing is calculated by the Monte-Carlo technique to capture the variations in the spectral transmittance and index of refraction of the participating media. Simulation results are presented for the growth interface morphology, thermal stresses in the crystal, and growth rate at various stages of growth. The influence of impurity gases on the growth process is also discussed, The simulation results are compared with the experimental observations. The capabilities of the developed tool in addressing other PVT growth processes such as SiC are discussed.

Worlikar, A.↗