Search NASA⌕ Search

SEARCH · Search NASA

Results for “['Monte”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 325 records · Page 18

Application of Monte Carlo techniques to transient thermal modeling of cavity radiometers having diffuse-specular surfaces

A viable alternative to the net exchange method of radiative analysis which is equally applicable to diffuse and diffuse-specular enclosures is presented. It is particularly more advantageous to use than the net exchange method in the case of a transient thermal analysis involving conduction and storage of energy as well as radiative exchange. A new quantity, called the distribution factor is defined which replaces the angle factor and the configuration factor. Once obtained, the array of distribution factors for an ensemble of surface elements which define an enclosure permits the instantaneous net radiative heat fluxes to all of the surfaces to be computed directly in terms of the known surface temperatures at that instant. The formulation of the thermal model is described, as is the determination of distribution factors by application of a Monte Carlo analysis. The results show that when fewer than 10,000 packets are emitted, an unsatisfactory approximation for the distribution factors is obtained, but that 10,000 packets is sufficient.

Mahan, J. R.↗

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.↗

Chemical application of diffusion quantum Monte Carlo

The diffusion quantum Monte Carlo (QMC) method gives a stochastic solution to the Schroedinger equation. This approach is receiving increasing attention in chemical applications as a result of its high accuracy. However, reducing statistical uncertainty remains a priority because chemical effects are often obtained as small differences of large numbers. As an example, the single-triplet splitting of the energy of the methylene molecule CH sub 2 is given. The QMC algorithm was implemented on the CYBER 205, first as a direct transcription of the algorithm running on the VAX 11/780, and second by explicitly writing vector code for all loops longer than a crossover length C. The speed of the codes relative to one another as a function of C, and relative to the VAX, are discussed. The computational time dependence obtained versus the number of basis functions is discussed and this is compared with that obtained from traditional quantum chemistry codes and that obtained from traditional computer architectures.

Reynolds, P. J.↗

A highly optimized vectorized code for Monte Carlo simulations of SU(3) lattice gauge theories

New methods are introduced for improving the performance of the vectorized Monte Carlo SU(3) lattice gauge theory algorithm using the CDC CYBER 205. Structure, algorithm and programming considerations are discussed. The performance achieved for a 16(4) lattice on a 2-pipe system may be phrased in terms of the link update time or overall MFLOPS rates. For 32-bit arithmetic, it is 36.3 microsecond/link for 8 hits per iteration (40.9 microsecond for 10 hits) or 101.5 MFLOPS.

Barkai, D.↗

Analysis of large solid propellant rocket engine exhaust plumes using the direct simulation Monte Carlo method

A new solution procedure has been developed to analyze the flowfield properties in the vicinity of the Inertial Upper Stage/Spacecraft during the 1st stage (SRMI) burn. Continuum methods are used to compute the nozzle flow and the exhaust plume flowfield as far as the boundary where the breakdown of translational equilibrium leaves these methods invalid. The Direct Simulation Monte Carlo (DSMC) method is applied everywhere beyond this breakdown boundary. The flowfield distributions of density, velocity, temperature, relative abundance, surface flux density, and pressure are discussed for each species for 2 sets of boundary conditions: vacuum and freestream. The interaction of the exhaust plume and the freestream with the spacecraft and the 2-stream direct interaction are discussed. The results show that the low density, high velocity, counter flowing free-stream substantially modifies the flowfield properties and the flux density incident on the spacecraft. A freestream bow shock is observed in the data, located forward of the high density region of the exhaust plume into which the freestream gas does not penetrate. The total flux density incident on the spacecraft, integrated over the SRM1 burn interval is estimated to be of the order of 10 to the 22nd per sq m (about 1000 atomic layers).

Hueser, J. E.↗

A Monte Carlo comparison of the recovery of winds near upwind and downwind from the SASS-1 model function by means of the sum of squares algorithm and a maximum likelihood estimator

Backscatter measurements at upwind and crosswind are simulated for five incidence angles by means of the SASS-1 model function. The effects of communication noise and attitude errors are simulated by Monte Carlo methods, and the winds are recovered by both the Sum of Square (SOS) algorithm and a Maximum Likelihood Estimater (MLE). The SOS algorithm is shown to fail for light enough winds at all incidence angles and to fail to show areas of calm because backscatter estimates that were negative or that produced incorrect values of K sub p greater than one were discarded. The MLE performs well for all input backscatter estimates and returns calm when both are negative. The use of the SOS algorithm is shown to have introduced errors in the SASS-1 model function that, in part, cancel out the errors that result from using it, but that also cause disagreement with other data sources such as the AAFE circle flight data at light winds. Implications for future scatterometer systems are given.

Pierson, W. J., Jr.↗

Monte Carlo Modeling of Lunar Megaregolith Development

A Monte Carlo model of lunar megaregolith development is proposed. Minimum megaregolith depths are obtained in the preliminary model. The model impacts a gridded, initially flat surface area randomly, with points not constrained to grid intersections. Information is recorded at grid points on evolving topography, deepest depth impacted to, and number of hits at that location. Two types of lunar surfaces ar simulated, maria and highlands. Crater production functions are after run to densities that are multiples of observed lunar crater densities for craters of 8 km in diameter and larger. all craters are considered to have formed in a gravity scaled regime. Edge effect problems are corrected in this model by allowing craters forming outside the grid area or along its edges to effect the gridded area.

Cashore, J.↗

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.↗

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.↗

Result of Monte-Carlo simulation of electron-photon cascades in lead and layers of lead-scintillator

Results of Monte-Carlo simulation of electromagnetic cascade development in lead and lead-scintillator sandwiches are analyzed. It is demonstrated that the structure function for core approximation is not applicable in the case in which the primary energy is higher than 100 GeV. The simulation data has shown that introducing an inhomogeneous chamber structure results in subsequent reduction of secondary particles.

Wasilewski, A.↗

A three-dimensional Monte Carlo calculation of the photon initiated showers and Kiel result

The Kiel experimental results indicate an existence of the ultra high-energy gamma-rays coming from Cyg. X-3. However the result indicates that the number of the muons included in the photon initiated shower is the same as the number included in the proton initiated showers. According to our Monte Carlo calculation as shown in the graph of underpart, the number of muons included in the photon initiated showers should be less than 1/15 of the photon's. The previous simulation was made under one dimensional approximation. This time the result of three dimensional calculation is reported.

Okada, A.↗

The research program of the Liquid Scintillation Detector (LSD) in the Mont Blanc Laboratory

A massive (90 tons) liquid scintillation detector (LSD) has been running since October 1984 in the Mont Blanc Laboratory at a depth of 5,200 hg/sq cm of standard rock. The research program of the experiment covers a variety of topics in particle physics and astrophysics. The performance of the detector, the main fields of research are presented and the preliminary results are discussed.

Dadykin, V. L.↗

Regularized Monte-Carlo calculations of the optical polarization induced by scattering in X-ray binaries

The results of regularized Monte Carlo calculations of the optical polarization Stokes parameters expected as a result of scattering from gas streams or accretion wakes in close binary systems whose properties resemble those of X-ray binaries are presented. The physical processes occurring in the system are represented by a model which follows a selected sample of individual photon paths from each emitting element on the surface of the primary to the photon's detection at earth. The result represents the variations in the intensity from each individual emitting element caused by the effects of limb darkening, gravity darkening, and variation in the projected area of the emitting element as seen from the scattering element.

Dolan, J. F.↗

Monte Carlo shock-like solutions to the Boltzmann equation with collective scattering

The results of Monte Carlo simulations of steady state shocks generated by a collision operator that isotropizes the particles by means of elastic scattering in some locally defined frame of reference are presented. The simulations include both the back reaction of accelerated particles on the inflowing plasma and the free escape of high-energy particles from finite shocks. Energetic particles are found to be naturally extracted out of the background plasma by the shock process with an efficiency in good quantitative agreement with an earlier analytic approximation (Eichler, 1983 and 1984) and observations (Gosling et al., 1981) of the entire particle spectrum at a quasi-parallel interplanetary shock. The analytic approximation, which allows a self-consistent determination of the effective adiabatic index of the shocked gas, is used to calculate the overall acceleration efficiency and particle spectrum for cases where ultrarelativistic energies are obtained. It is found that shocks of the strength necessary to produce galactic cosmic rays put approximately 15 percent of the shock energy into relativistic particles.

Ellison, D. C.↗

Monte Carlo simulation of a new gamma ray telescope

A new Monte Carlo code has been written to simulate the response of the new University of California double scatter gamma ray telescope. This package of modular software routines, written in VAX FORTRAN 77 simulates the detection of 0.1 to 35 MeV gamma rays. The new telescope is flown from high altitude balloons to measure medium energy gamma radiation from astronomical sources. This paper presents (1) the basic physics methods in the code, and (2) the predicted response functions of the telescope. Gamma ray processes include Compton scattering, pair production and photoelectric absorption in plastic scintillator, NaI(Tl) and aluminum. Electron transport processes include ionization energy loss, multiple scattering, production of bremsstrahlung photons and positron annihilation.

Simone, J.↗