SEARCH · Search NASA

Results for “monte carlo”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1970

FASTER 3: A generalized-geometry Monte Carlo computer program for the transport of neutrons and gamma rays. Volume 2: Users manual

A description of the FASTER-III program for Monte Carlo Carlo calculation of photon and neutron transport in complex geometries is presented. Major revisions include the capability of calculating minimum weight shield configurations for primary and secondary radiation and optimal importance sampling parameters. The program description includes a users manual describing the preparation of input data cards, the printout from a sample problem including the data card images, definitions of Fortran variables, the program logic, and the control cards required to run on the IBM 7094, IBM 360, UNIVAC 1108 and CDC 6600 computers.

Jordan, T. M.↗

NASA NTRS1985

Monte Carlo simulations in support of the Shuttle upper atmospheric mass spectrometer experiment

This paper presents the results of calculations obtained with a Direct Simulation Monte Carlo (DSMC) method that describes both the external flow about the nose region of the Shuttle Orbiter and the internal flow in an inlet tube that provides the initial path between the shock-processed gases and a mass spectrometer mounted inside the Orbiter. A dedicated, three-dimensional version of the DSMC was developed during this study for the internal flow simulation along with the procedure for interfacing the external and internal flows. The calculations span the 140 to 95 km altitude range, which includes most of the hypersonic transitional flow encountered during reentry. All calculations were for a multicomponent gas mixture consisting of five chemical species while simulating the effects of transitional, rotational, vibrational, and chemical nonequilibrium. The results of the calculations show that within the entry region of the inlet tube where the gas has equilibrated with the sidewall tube temperature, the pressure is substantially less than the pressure at the external surface. This pressure correction for the entry region is significant for all conditions investigated and increases with altitude. The results highlight the structure of both the external and internal flow. Results of parametric studies show the sensitivity of the equilibrated inlet sidewall pressure to mass flow rates, gas-surface reflection model, tube side-wall temperature variations, and surface recombinations.

Moss, J. N.↗

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 NTRS2004

Monte Carlo Techniques for Calculations of Charge Deposition and Displacement Damage from Protons in Visible and Infrared Sensor Arrays

This slide presentation examines motivation for Monte Carlo methods, charge deposition in sensor arrays, displacement damage calculations, and future work. The discussion of charge deposition sensor arrays includes Si active pixel sensor APS arrays and LWIR HgCdTe FPAs. The discussion of displacement damage calculations includes nonionizing energy loss (NIEL), HgCdTe NIEL calculation results including variance, and implications for damage in HgCdTe detector arrays.

Marshall, Paul↗

NASA NTRS2002

Status of the Space Radiation Monte Carlos Simulation Based on FLUKA and ROOT

The NASA-funded project reported on at the first IWSSRR in Arona to develop a Monte-Carlo simulation program for use in simulating the space radiation environment based on the FLUKA and ROOT codes is well into its second year of development, and considerable progress has been made. The general tasks required to achieve the final goals include the addition of heavy-ion interactions into the FLUKA code and the provision of a ROOT-based interface to FLUKA. The most significant progress to date includes the incorporation of the DPMJET event generator code within FLUKA to handle heavy-ion interactions for incident projectile energies greater than 3GeV/A. The ongoing effort intends to extend the treatment of these interactions down to 10 MeV, and at present two alternative approaches are being explored. The ROOT interface is being pursued in conjunction with the CERN LHC ALICE software team through an adaptation of their existing AliROOT software. As a check on the validity of the code, a simulation of the recent data taken by the ATIC experiment is underway.

Andersen, Victor↗

NASA NTRS1969

Monte Carlo direct view factor and generalized radiative heat transfer programs

Computer programs find the direct view factor from one surface segment to another using the Monte carlo technique, and the radioactive-transfer coefficients between surface segments. An advantage of the programs is the great generality of problems treatable and rapidity of solution from problem conception to receipt of results.

Mc Williams, J. L.↗

NASA NTRS1974

Nonlinear panel response by a Monte Carlo approach

The vibration of clamped and simply supported elastic panels due to subsonic and supersonic turbulent boundary-layer flows is investigated by a Monte Carlo technique. The resulting generalized random forces are simulated numerically from boundary-layer turbulence spectra, and the response analysis is performed in the time domain. The mutual interaction between panel motion and external and/or internal airflow is included. Response studies are performed with respect to rms response, probability structure, peak distribution, threshold crossing and spectral density. The effect on the response statistics of in-plane loading, static pressure differential, and cavity pressure is investigated.

Vaicaitis, R.↗

NASA NTRS1982

An analysis of the OI 1304 A dayglow using a Monte Carlo resonant scattering model with partial frequency redistribution

The transport of resonance radiation under optically thick conditions is shown to be accurately described by a Monte Carlo model of the atomic oxygen 1304 A airglow triplet in which partial frequency redistribution, temperature gradients, pure absorption and multilevel scattering are accounted for. All features of the data can be explained by photoelectron impact excitation and the resonant scattering of sunlight, where the latter source dominates below 100 and above 500 km and is stronger at intermediate altitudes than previously thought. It is concluded that the OI 1304 A emission can be used in studies of excitation processes and atomic oxygen densities in planetary atmospheres.

Meier, R. R.↗

NASA NTRS1989

Monte Carlo computations of F-region incoherent radar spectra at high latitudes and the use of a simple method for non-Maxwellian spectral calculations

Monte Carlo simulations of ion velocity distributions in the high-latitude F region have been performed in order to improve the calculation of incoherent radar spectra in the auroral ionosphere. The results confirm that when the ion temperature becomes large due to frictional heating in the presence of collisions with the neutral background constituent, F region spectra evolve from a normal double hump, to a triple hump, to a spectrum with a single maximum. An empirical approach is developed to overcome the inadequacy of the Maxwellian assumption for the case of radar aspect angles of between 30 and 70 deg.

Kikuchi, K.↗

NASA NTRS1995

Theoretical study of hole initiated impact ionization in bulk silicon and GaAs using a wave-vector-dependent numerical transition rate formulation within an ensemble Monte Carlo calculation

In this paper, calculations of the hole initiated interband impact ionization rate in bulk silicon and GaAs are presented based on an ensemble Monte Carlo simulation with the inclusion of a wave-vector-dependent numerical transition rate formulation. The ionization transition rate is determined for each of the three valence bands, heavy, light, and split-off, using Fermi's golden rule with a two-body, screened Coulomb interaction. The dielectric function used within the calculation is assumed to be wave-vector-dependent. Calculations of the field-dependent impact ionization rate as well as the quantum yield are presented. It is found from both the quantum yield results and examination of the hole distribution function that the effective threshold energy for hole initiated impact ionization is relatively soft, similar to that predicted for the corresponding electron initiated ionization events occur more frequently than either heavy or split-offf initiated ionization events in bulk silicon over the applied electric field strengths examined here, 250-500 kV/cm. Conversely,in GaAs, the vast majority of hole initated ionization events originate from holes within the split-off band.

Oguzman, Ismail H.↗

NASA NTRS1983

Lagged average forecasting, an alternative to Monte Carlo forecasting

A 'lagged average forecast' (LAF) model is developed for stochastic dynamic weather forecasting and used for predictions in comparison with the results of a Monte Carlo forecast (MCF). The technique involves the calculation of sample statistics from an ensemble of forecasts, with each ensemble member being an ordinary dynamical forecast (ODF). Initial conditions at a time lagging the start of the forecast period are used, with varying amounts of time for the lags. Forcing by asymmetric Newtonian heating of the lower layer is used in a two-layer, f-plane, highly truncated spectral model in a test forecasting run. Both the LAF and MCF are found to be more accurate than the ODF due to ensemble averaging with the MCF and the LAF. When a regression filter is introduced, all models become more accurate, with the LAF model giving the best results. The possibility of generating monthly or seasonal forecasts with the LAF is discussed.

Hoffman, R. N.↗

NASA NTRS1985

Monte Carlo calculations of relativistic solar proton propagation in interplanetary space

Particle fluxes and pitch angle distributions of relativistic solar protons at 1 AU were determined by Monte Carlo calculations. The analysis covers two hours after the release of the particles from the Sun and total of eight 100000 particle trajectories were simulated. The pitch angle scattering was assumed to be isotropic ad the scattering mean free path was varied from 0.1 to 4 AU. As an application, the solar injection time and interplanetary scattering mean free path of particles that gave rise to the GLE on May, 1978 were determined. Assuming exponential form, the injection decay time was found to be about 11 minutes. The m.f.p. of pitch angle scattering during the event was about 1 AU.

Lumme, M.↗

NASA NTRS1985

EAS accompanied by gamma-families at Mt. Norikura and comparison with Monte Carlo simulation

The experimental data of extensive air showers (EAS) accompanied by gamma-families, with total energy greater than 10 TeV, were compared with a Monte Carlo simulation with a rising cross section proportional to E sup 0.04 for the p-air inelastic cross section. It is found that the absolute intensity of size spectrum of such EAS is strongly affected by the primary protons intensity at 10 to the 15th power approx. 10 to the 16th power eV region and the experimental size spectrum agrees with the simulated spectra for the p-poor primary composition better than the p-rich one.

Shima, M.↗

NASA NTRS1987

Acoustic fatigue - A Monte Carlo approach

The present paper reviews the basic concepts of nonlinear structural response to surface flow and acoustic pressure inputs. A time domain solution and the Monte Carlo method are utilized to estimate the required statistics of the nonlinear response. Simulation procedures of multi-dimensional and multi-variate random processes are discussed. Application of these procedures to acoustic fatigue of flight structures are considered. Numerical examples are included to illustrate the applicability of this method for the solution of nonlinear problems.

Vaicaitis, R.↗

NASA NTRS1989

Monte Carlo calculation of dynamical properties of the two-dimensional Hubbard model

A new method is introduced for analytically continuing imaginary-time data from quantum Monte Carlo calculations to the real-frequency axis. The method is based on a least-squares-fitting procedure with constraints of positivity and smoothness on the real-frequency quantities. Results are shown for the single-particle spectral-weight function and density of states for the half-filled, two-dimensional Hubbard model.

White, S. R.↗

NASA NTRS1993

Comparison of generalized transport and Monte-Carlo models of the escape of a minor species

The steady-state diffusion of a minor species through a static background species is studied using a Monte Carlo model and a generalized 16-moment transport model. The two models are in excellent agreement in the collision-dominated region and in the 'transition region'. In the 'collisionless' region the 16-moment solution contains two singularities, and physical meaning cannot be assigned to the solution in their vicinity. In all regions, agreement between the models is best for the distribution function and for the lower-order moments and is less good for higher-order moments. Moments of order higher than the heat flow and hence beyond the level of description provided by the transport model have a noticeable effect on the shape of distribution functions in the collisionless region.

Demars, H. G.↗

NASA NTRS1976

Multistream and Monte Carlo calculations of the sun's aureole

The scattered UV intensities in the sun's aureole are calculated using both a multi-stream scattering method and a Monte Carlo approach. Angular distributions for both Rayleigh and aerosol scatterings are obtained with a realistic atmospheric model. Moderate- and sharply-forward-peaked phase functions of aerosol scattering, corresponding to realistic analytic size distributions, are incorporated. Results from the two independent calculations are in reasonable agreement for a realistic atmospheric model. The results indicate that the scattered UV intensity in the sun's aureole is about four orders smaller than the transmitted intensity while the scattered intensity for pure Rayleigh scattering is only 1-10 millionths times that of the transmitted sun. Based on these intensity ratios, we estimate that the integrated scattered contributions from the aureole to a well-collimated sun photometer of acceptance aperture 2-3 deg are below about 1% of the direct contribution.

Furman, D. R.↗

NASA NTRS1983

A potential energy scaling Monte Carlo simulation of thin film nucleation and growth

The initial growth of thin Ge fims on the (100) surface of an Fe substrate is investigated theoretically by means of Monte Carlo simulations based on a potential-energy-scaling technique. The substrate is modeled as a 20 x 20-square array with periodic boundary conditions, as described by Heinbockel et al. (1983), and the movement of surface atoms under the influence of the substrate interaction potential and the lateral interaction of neighboring atoms is explored via continuous updating (on the time scale of single events) of the potential energy at each site in the array. Results for the clustering of nine dispersed atoms over 1.0 s at 600 K and for deposition at 5 x 10 to the -14th/sq cm s over 2.0 s at 500 K are presented graphically.

Outlaw, R. A.↗