Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion 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.

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 NTRS1994

Enhanced cosmic-ray acceleration rates in highly inclined astrophysical shocks

The theory of cosmic-ray acceleration by subluminal magnetohydrodynamics (MHD) astrophysical shocks according to the test particle approximation is extended to highly oblique shocks where flow speeds which appear nonrelativistic in the shock rest frame appear on transformation to the de Hoffmann-Teller or E identically equal to 0 frame to have upstream flow speeds approaching c. Monte Carlo simulation shows that relative to the predictions of diffusion theory, as the upstream E identically equal to 0 frame velocity approaches c, flatter spectra and faster acceleration rates occur. These spectral and acceleration time changes are similar to those found for relativistic, parallel MHD shocks and may affect all nonthermal active galactic nuclei emission from relativistic electrons. Also there is an additional means of increasing the upper limit of the cosmic-ray spectrum expected from active galactic nuclei, although the approximation used may not accurately reproduce the spectral shape.

Lieu, R.↗
NASA NTRS1973

Radiation configuration factors for annular rings and hemispherical sectors.

In the application of Monte Carlo computer techniques to the study of problems in the free molecular flow regime, the analogy between radiation heat transfer with diffuse emissions and the free molecular flow field makes it possible to use the radiation configuration factors as a test on the operation of the Monte Carlo program. A source plane which was an annular ring of arbitrary radius around a sphere with the center of the disk coincident with the center of the sphere was examined. Since various indicators seemed to verify the correct operation of the program, the program was used to calculate the radiation configuration factors. The study adds a new geometry to the general literature on radiation configuration factors.

Ballance, J. O.↗
NASA NTRS1990

The effect of catastrophic collisional fragmentation and diffuse medium accretion on a computational interstellar dust system

The effects of catastrophic collisional fragmentation and diffuse medium accretion on a the interstellar dust system are computed using a Monte Carlo computer model. The Monte Carlo code has as its basis an analytic solution of the bulk chemical evolution of a two-phase interstellar medium, described by Liffman and Clayton (1989). The model is subjected to numerous different interstellar processes as it transfers from one interstellar phase to another. Collisional fragmentation was found to be the dominant physical process that shapes the size spectrum of interstellar dust. It was found that, in the diffuse cloud phase, 90 percent of the refractory material is locked up in the dust grains, primarily due to accretion in the molecular medium. This result is consistent with the observed depletions of silicon. Depletions were found to be affected only slightly by diffuse cloud accretion.

Liffman, Kurt↗
NASA NTRS1979

Comptonization of X-rays by low-temperature electrons

A method is described for calculating the spectrum that results from the Compton scattering of a monochromatic source of X-rays by low-temperature electrons, both for initial-value relaxation problems and for steady-state spatial diffusion problems. The method gives an exact solution of the inital-value problem for evolution of the spectrum in an infinite homogeneous medium if Klein-Nishina corrections to the Thomson cross section are neglected. This, together with approximate solutions for problems in which Klein-Nishina corrections are significant and/or spatial diffusion occurs, shows spectral structure near the original photon wavelength that may be used to infer physical conditions in cosmic X-ray sources. Explicit results, shown for examples of time relaxation in an infinite medium and spatial diffusion through a uniform sphere, are compared with results obtained by Monte Carlo calculations and by solving the appropriate Fokker-Planck equation.

Illarionov, A.↗
NASA NTRS2011

Computing Temperatures in Optically Thick Protoplanetary Disks

We worked with a Monte Carlo radiative transfer code to simulate the transfer of energy through protoplanetary disks, where planet formation occurs. The code tracks photons from the star into the disk, through scattering, absorption and re-emission, until they escape to infinity. High optical depths in the disk interior dominate the computation time because it takes the photon packet many interactions to get out of the region. High optical depths also receive few photons and therefore do not have well-estimated temperatures. We applied a modified random walk (MRW) approximation for treating high optical depths and to speed up the Monte Carlo calculations. The MRW is implemented by calculating the average number of interactions the photon packet will undergo in diffusing within a single cell of the spatial grid and then updating the packet position, packet frequencies, and local radiation absorption rate appropriately. The MRW approximation was then tested for accuracy and speed compared to the original code. We determined that MRW provides accurate answers to Monte Carlo Radiative transfer simulations. The speed gained from using MRW is shown to be proportional to the disk mass.

radiative transfer↗
NASA NTRS1983

Monte Carlo simulation of steady state shock structure including cosmic ray mediation and particle escape

Both hydrodynamic calculations (Drury and Volk, 1981, and Axford et al., 1982) and kinetic simulations imply the existence of thermal subshocks in high-Mach-number cosmic-ray-mediated shocks. The injection efficiency of particles from the thermal background into the diffusive shock-acceleration process is determined in part by the sharpness and compression ratio of these subshocks. Results are reported for a Monte Carlo simulation that includes both the back reaction of accelerated particles on the inflowing plasma, producing a smoothing of the shock transition, and the free escape of particles allowing arbitrarily large overall compression ratios in high-Mach-number steady-state shocks. Energy spectra and estimates of the proportion of thermal ions accelerated to high energy are obtained.

Ellison, D. C.↗
NASA NTRS1990

Growth morphology with anisotropic surface kinetics

The morphological evolution of crystals growing from an incongruent vapor phase is studied using a Monte Carlo model, and the full range of growth morphologies is recovered. The diffusion in the bulk nutrient and the anisotropy in the interface kinetics are morphologically destabilizing and stabilizing, respectively. For a given set of simulation parameters and lattice symmetries there is a critical size, which scales linearly with the mean free path in the vapor, beyond which a crystal cannot retain its stable, macroscopically faceted growth shape. Surface diffusion stabilizes faceted growth on the shorter scale of the mean surface diffusion length. In simulations with a uniform drift superimposed on the random walk nutrient transport, crystal faces oriented toward the drift show enhanced morphological stability compared to the purely diffusive situation. Rotational drifts with periodic reversal of direction are morphologically stabilizing for all crystal facets.

Xiao, Rong-Fu↗
NASA NTRS1988

Muon motion in titanium hydride

Motional narrowing of the transverse-field muon spin rotation signal was observed in gamma-TiH(x) for x = 1.83, 1.97, and 1.99. An analysis of the data for TiH1.99 near room temperature indicates that the mechanism responsible for the motion of the muon out of the octahedral site is thermally activated diffusion with an attempt frequency comparable to the optical vibrations of the lattice. Monte Carlo calculations to simulate the effect of muon and proton motion upon the muon field-correlation time were used to interpret the motional narrowing in TiH1.97 near 500 K. The interpretation is dependent upon whether the Bloembergen, Purcell, and Pound (BPP) theory or an independent spin-pair relaxation model is used to obtain the vacancy jump rate from proton NMR T1 measurements. Use of BPP theory shows that the field-correction time can be obtained if the rate of motion of the muon with respect to the rate of the motion for the protons is decreased. An independent spin-pair relaxation model indicates that the field-correlation time can be obtained if the rate of motion for the nearest-neighbor protons is decreased.

Kempton, J. R.↗
NASA NTRS1991

Comptonization of gamma rays by cold electrons

An analytic method is developed for calculating the emergent spectrum of gamma-rays and X-rays scattered in a homogeneous medium with low-temperature electrons. The Klein-Nishina corrections of the scattering cross section and absorption processes are taken in account. The wavelength relaxation and the spatial diffusion problems are solved separately, and the emergent spectrum is calculated by convolving the evolution function of the spectrum in an infinite medium with the photon luminosity resulting from the spatial diffusion in a finite sphere. The analytic results are compared with that of Monte Carlo calculations and it is concluded that the analytic result is quite accurate.

Xu, Yueming↗
NASA NTRS2012

Multi-Scale Modeling of Liquid Phase Sintering Affected by Gravity: Preliminary Analysis

A multi-scale simulation concept taking into account impact of gravity on liquid phase sintering is described. The gravity influence can be included at both the micro- and macro-scales. At the micro-scale, the diffusion mass-transport is directionally modified in the framework of kinetic Monte-Carlo simulations to include the impact of gravity. The micro-scale simulations can provide the values of the constitutive parameters for macroscopic sintering simulations. At the macro-scale, we are attempting to embed a continuum model of sintering into a finite-element framework that includes the gravity forces and substrate friction. If successful, the finite elements analysis will enable predictions relevant to space-based processing, including size and shape and property predictions. Model experiments are underway to support the models via extraction of viscosity moduli versus composition, particle size, heating rate, temperature and time.

Olevsky, Eugene↗
NASA NTRS1978

Uncertainty propagation in a stratospheric model. I - Development of a concise stratospheric model. II - Monte Carlo analysis of imprecisions due to reaction rates

The paper discusses the development of a concise stratospheric model which uses iteration to obtain coupling between interacting species. The one-dimensional, steady-state, diurnally-averaged model generates diffusion equations with appropriate sources and sinks for species odd oxygen, H2O, H2, CO, N2O, odd nitrogen, CH4, CH3Cl, CCl4, CF2Cl2, CFCl3, and odd chlorine. The model evaluates steady-state perturbations caused by injections of chlorine and NO(x) and may be used to predict ozone depletion. The model is used in a Monte Carlo study of the propagation of reaction-rate imprecisions by calculating an ozone perturbation caused by the addition of chlorine. Since the model is sensitive to only 10 of the more than 50 reaction rates considered, only about 1000 Monte Carlo cases are required to span the space of possible results.

Rundel, R. D.↗
NASA NTRS1973

New approach to cosmic-ray diffusion theory.

We have investigated a new approach to deriving a diffusion equation for charged particles in a static, random magnetic field. Our method incorporates essential effects of the magnetic fluctuations in the lowest order particle orbits. Significant corrections to the usual quasilinear diffusion coefficient for cosmic rays with pitch angles near 90 deg are a consequence. Monte Carlo results bear out the validity of our theory.

Jones, F. C.↗
NASA NTRS1995

A deterministic particle method for one-dimensional reaction-diffusion equations

We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution.

Mascagni, Michael↗
NASA NTRS1987

A microstructural lattice model for strain oriented problems: A combined Monte Carlo finite element technique

A specialized, microstructural lattice model, termed MCFET for combined Monte Carlo Finite Element Technique, was developed which simulates microstructural evolution in material systems where modulated phases occur and the directionality of the modulation is influenced by internal and external stresses. In this approach, the microstructure is discretized onto a fine lattice. Each element in the lattice is labelled in accordance with its microstructural identity. Diffusion of material at elevated temperatures is simulated by allowing exchanges of neighboring elements if the exchange lowers the total energy of the system. A Monte Carlo approach is used to select the exchange site while the change in energy associated with stress fields is computed using a finite element technique. The MCFET analysis was validated by comparing this approach with a closed form, analytical method for stress assisted, shape changes of a single particle in an infinite matrix. Sample MCFET analytical for multiparticle problems were also run and in general the resulting microstructural changes associated with the application of an external stress are similar to that observed in Ni-Al-Cr alloys at elevated temperature.

Gayda, J.↗
NASA NTRS1981

A radiometric model of an earth radiation budget radiometer optical system with diffuse-specular surfaces

The Earth Radiation Budget Experiment (ERBE) is to fly on NASA's Earth Radiation Budget Satellite (ERBS) and on NOAA F and NOAA G. Large spatial scale earth energy budget data will be derived primarily from measurements made by the ERBE nonscanning instrument (ERBE-NS). A description is given of a mathematical model capable of simulating the radiometric response of any of the ERBE-NS earth viewing channels. The model uses a Monte Carlo method to accurately account for directional distributions of emission and reflection from optical surfaces which are neither strictly diffuse nor strictly specular. The model computes radiation exchange factors among optical system components, and determines the distribution in the optical system of energy from an outside source. Attention is also given to an approach for implementing the model and results obtained from the implementation.

Luther, M. R.↗
NASA NTRS1991

A study of hydrogen diffusion flames using PDF turbulence model

The application of probability density function (pdf) turbulence models is addressed. For the purpose of accurate prediction of turbulent combustion, an algorithm that combines a conventional computational fluid dynamic (CFD) flow solver with the Monte Carlo simulation of the pdf evolution equation was developed. The algorithm was validated using experimental data for a heated turbulent plane jet. The study of H2-F2 diffusion flames was carried out using this algorithm. Numerical results compared favorably with experimental data. The computations show that the flame center shifts as the equivalence ratio changes, and that for the same equivalence ratio, similarity solutions for flames exist.

Hsu, Andrew T.↗