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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

PDF approach for compressible turbulent reacting flows

The objective of the present work is to develop a probability density function (pdf) turbulence model for compressible reacting flows for use with a CFD flow solver. The probability density function of the species mass fraction and enthalpy are obtained by solving a pdf evolution equation using a Monte Carlo scheme. The pdf solution procedure is coupled with a compressible CFD flow solver which provides the velocity and pressure fields. A modeled pdf equation for compressible flows, capable of capturing shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed, and an averaging procedure is developed to provide smooth Monte-Carlo solutions to ensure convergence. Two supersonic diffusion flames are studied using the proposed pdf model and the results are compared with experimental data; marked improvements over CFD solutions without pdf are observed. Preliminary applications of pdf to 3D flows are also reported.

Hsu, A. T.↗
NASA NTRS1986

The diffusion of individual molecules within a gas

The Direct Simulation Monte Carlo method is used to study the positional history of the individual molecules in a gas that is homogeneous at the macroscopic level and is in Maxwellian equilibrium at the microscopic level. The behavior at small times is characterized by 'persistence of velocity' effects, and a 'random walk' type of dispersal occurs over a longer timescale. It is shown that the rate of dispersal can be directly related to the self-diffusion coefficient. In addition, the diffusion coefficients are obtained directly from one-dimensional calculations, and the local Knudsen number at which the Chapman-Enskog theory breaks down is determined. Results are presented for both simple gases and gas mixtures.

Bird, G. A.↗
NASA NTRS1990

Electron heating in quasi-perpendicular shocks - A Monte Carlo simulation

To study the problem of electron heating in quasi-perpendicular shocks, under the combined effects of 'reversible' motion, in the shock electric potential and magnetic field, and wave-particle interactions a diffusion equation was derived, in the drift (adiabatic) approximation and it was solved by using a Monte Carlo method. The results show that most of the observations can be explained within this framework. The simulation has also definitively shown that the electron parallel temperature is determined by the dc electromagnetic field and not by any wave particle induced heating. Wave-particle interactions are effective in smoothing out the large gradients in phase space produced by the 'reversible' motion of the electrons, thus producing a 'cooling' of the electrons. Some constraints on the wave-particle interaction process may be obtained from a detailed comparison between the simulation and observations. In particular, it appears that the adiabatic approximation must be violated in order to explain the observed evolution of the perpendicular temperature.

Veltri, Pierluigi↗
NASA NTRS2018

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan-1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017). Themodel reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U(o) = 0.5eV, U(w) = 0.078 – 0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

O. J. Tucker↗
NASA NTRS2018

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan‐1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017, https://doi.org/10.1126/sciadv.1701471). The model reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U0 = 0.5 eV, UW = 0.078–0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

O J Tucker↗
NASA NTRS2018

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan‐1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017, https://doi.org/10.1126/sciadv.1701471). The model reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U(sub 0) = 0.5 eV, U(sub W) = 0.078–0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

Tucker, O. J.↗
NASA NTRS1993

Simulations of particle acceleration in parallel shocks: Direct comparison between Monte Carlo and one-dimensional hybrid codes

We have made a direct comparison between two different computer simulations of a plane, parallel, collisionless shock including particle acceleration to energies typical of those of diffuse ions observed at the earth bow shock. Despite the fact that the one-dimensional hybrid and Monte Carlo techniques employ entirely different algorithms, they give surprisingly close agreement in the overall shapes of the complete distribution functions for protons as well as heavier ions. Both methods show that energetic ions emerge smoothly from the background thermal plasma with approximately the same relative injection rate and that the fraction of the incoming plasma's energy flux that is converted into downstream enthalpy flux of the accelerated population (i.e., the acceleration efficiency) is similar in the two cases. The fraction of the downstream proton distribution made up of superthermal particles is quite large, with at least 10% of the energy flux going into protons with energies above 10 keV. In addition, an upstream precursor, produced by backstreaming energetic particles, is present in both shocks, although the Monte Carlo precursor is considerably longer than that produced in the hybrid shock. These results offer convincing evidence that, at least in these ways, the two simulations are consistent in their description of parallel shock structure and particle acceleration, and they lay the groundwork for development of shock models employing a combination of both methods.

Ellison, Donald C.↗
NASA NTRS1981

Monte-Carlo model of pitch-angle scattering in solar cosmic ray events

Monte Carlo simulations of the propagation of solar cosmic-rays in interplanetary space are reported, including the effects of pitch-angle scattering and adiabatic focusing. Time intensity profiles agree well with the corresponding spatial diffusion approximation for models close to those used by Gombosi and Owens (1980). Monte Carlo simulations of the same model problem yield results that disagree with those of Gombosi and Owens, and the most likely possible source of error in the Monte Carlo simulation is considered to be in the region near the 90 deg pitch angle where the finite step size introduces certain inconsistencies.

Palmer, I. D.↗
NASA NTRS1987

Simulation of charged particle transport on the MPP

Computations of cosmic ray transport based upon finite difference methods are afflicted by instabilities, inaccuracies, and artifacts. To avoid these problems, a Monte Carlo formulation was developed which is closely related not only to the finite difference formulation, but also to the underlying physics of transport phenomena. Implementations of this approach are currently running on the Massively Parallel Processor at Goddard, whose enormous computing power overwhelms the poor statistical accuracy that usually limits the use of stochastic methods. In a Monte Carlo simulation of rectilinear transport, the coherent and diffusive effects that appeared are in good quantitative agreement with both finite difference and analytical calculations.

Earl, James A.↗
NASA NTRS1991

Effect of wave-particle interaction on the outflow of ions at high latitudes

The objective of the research is to use the low-frequency wave spectrum measured by the Plasma Wave Instrument on the DE-1 spacecraft to include the wave-particle interaction (WPI) in the different polar wind models that are available at Utah State University. A Monte Carlo technique was used to simulate the ion diffusion in the velocity space due to scattering by the waves. This enabled us to study the effect of WPI on the magnitude as well as the composition of the outflow of the ionospheric ions. In particular, in the first year the 1-D semi-kinetic code developed by Barakat and Schunk (1983) was modified in order to include the effect of the WPI.

Barakat, Abdallah R.↗
NASA NTRS1994

Probability density function approach for compressible turbulent reacting flows

The objective of the present work is to extend the probability density function (PDF) tubulence model to compressible reacting flows. The proability density function of the species mass fractions and enthalpy are obtained by solving a PDF evolution equation using a Monte Carlo scheme. The PDF solution procedure is coupled with a compression finite-volume flow solver which provides the velocity and pressure fields. A modeled PDF equation for compressible flows, capable of treating flows with shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed. Two super sonic diffusion flames are studied using the proposed PDF model and the results are compared with experimental data; marked improvements over solutions without PDF are observed.

Hsu, A. T.↗
NASA NTRS1975

Monte Carlo treatment of Lyman-alpha. II - Radiation in a spherical atmosphere

Intensity and state of polarization of solar L-alpha photons as they diffuse through an inhomogeneous, spherically symmetric, isothermal geocorona are theoretically determined. The fine structure of L-alpha and Doppler redistribution of frequencies are taken into account. The calculation use the Monte Carlo technique involving Stokes vectors. Comparison of the results with OGO-4 and OSO-4 observed intensities at an altitude of 650 km shows good agreement. Calculations of the polarization versus solar zenith angle show a residual polarization at large zenith angles which is mainly due to multiply scattered photons.

Modali, S. B.↗
NASA NTRS1987

The origin of the diffuse X-ray background

The Einstein Deep Survey fields are used to reexamine the problem of the origin of the extragalactic diffuse X-ray background. An approach based on fluctuation analysis is adopted, and a sophisticated Monte Carlo scheme is used to compare the X-ray images with a variety of models representing point and diffuse source contributions. Results are presented for two classes of models: those which determine the flux cutoff beyond which pointlike sources fail to contribute significantly if the log N-log S slope is fixed at the measured value of -1.5, and those which set limits on the slope of the log N-log S relation below its last measured point at the Deep Survey threshold. In both cases, it is found that the known population of X-ray-emitting quasars and any reasonable extrapolations therefrom do not comprise the entire X-ray background.

Hamilton, T. T.↗
NASA NTRS2015

An Innovative Concept for Spacebased Lidar Measurement of Ocean Carbon Biomass

Beam attenuation coefficient, c, provides an important optical index of plankton standing stocks, such as phytoplankton biomass and total particulate carbon concentration. Unfortunately, c has proven difficult to quantify through remote sensing. Here, we introduce an innovative approach for estimating c using lidar depolarization measurements and diffuse attenuation coefficients from ocean color products or lidar measurements of Brillouin scattering. The new approach is based on a theoretical formula established from Monte Carlo simulations that links the depolarization ratio of sea water to the ratio of diffuse attenuation Kd and beam attenuation C (i.e., a multiple scattering factor). On July 17, 2014, the CALIPSO satellite was tilted 30Âdeg off-nadir for one nighttime orbit in order to minimize ocean surface backscatter and demonstrate the lidar ocean subsurface measurement concept from space. Depolarization ratios of ocean subsurface backscatter are measured accurately. Beam attenuation coefficients computed from the depolarization ratio measurements compare well with empirical estimates from ocean color measurements. We further verify the beam attenuation coefficient retrievals using aircraft-based high spectral resolution lidar (HSRL) data that are collocated with in-water optical measurements.

Hu, Yongxiang↗
NASA NTRS1978

Optically thick X-ray transfer - The shell game

This paper investigates the radiative transfer of X-rays through a shell that is optically thick to Compton scattering, surrounding a point source of continuum X-rays. The emission and absorption of X-rays due to K-shell transitions of iron are included. The calculations are done in two entirely independent ways: by Monte Carlo simulation and by solving a Fokker-Planck diffusion equation. The emergent spectra agree very well for Thomson depths of at least about 2. The validity is confirmed of the modification to the Fokker-Planck equation of Kompaneets (1957) that is required when the photon energy is large compared with the average thermal energy of the electrons. A procedure is also developed for treating models of compact X-ray sources consisting of incomplete shells.

Langer, S. H.↗