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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Radio to Gamma-Ray Emission from Shell-Type Supernova Remnants: Predictions from Non-Linear Shock Acceleration Models

Supernova remnants (SNRs) are widely believed to be the principal source of galactic cosmic rays, produced by diffusive shock acceleration in the environs of the remnant's expanding blast wave. Such energetic particles can produce gamma-rays and lower energy photons via interactions with the ambient plasma. The recently reported observation of TeV gamma-rays from SN1006 by the CANGAROO Collaboration, combined with the fact that several unidentified EGRET sources have been associated with known radio/optical/X-ray-emitting remnants, provides powerful motivation for studying gamma-ray emission from SNRs. In this paper, we present results from a Monte Carlo simulation of non-linear shock structure and acceleration coupled with photon emission in shell-like SNRs. These non-linearities are a by-product of the dynamical influence of the accelerated cosmic rays on the shocked plasma and result in distributions of cosmic rays which deviate from pure power-laws. Such deviations are crucial to acceleration efficiency considerations and impact photon intensities and spectral shapes at all energies, producing GeV/TeV intensity ratios that are quite different from test particle predictions.

Baring, Matthew G.↗
NASA NTRS1995

Orientation and velocity dependence of the nonequilibrium partition coefficient

Monte Carlo simulations based on a Spin-1 Ising Model for binary alloys have been used to investigate the non-equilibrium partition coefficient (k(sub neq)) as a function of solid-liquid interface velocity and orientation. In simulations of Si with a second component k(sub neq) is greater in the [111] direction than the [100] direction in agreement with experimental results reported by Azlz et al. The simulated partition coefficient scales with the square of the step velocity divided by the diffusion coefficient of the secondary component in the liquid.

Beatty, K. M.↗
NASA NTRS2013

Uncertainty Determination for Aeroheating in Uranus and Saturn Probe Entries by the Monte Carlo Method

The 2013-2022 Decaedal survey for planetary exploration has identified probe missions to Uranus and Saturn as high priorities. This work endeavors to examine the uncertainty for determining aeroheating in such entry environments. Representative entry trajectories are constructed using the TRAJ software. Flowfields at selected points on the trajectories are then computed using the Data Parallel Line Relaxation (DPLR) Computational Fluid Dynamics Code. A Monte Carlo study is performed on the DPLR input parameters to determine the uncertainty in the predicted aeroheating, and correlation coefficients are examined to identify which input parameters show the most influence on the uncertainty. A review of the present best practices for input parameters (e.g. transport coefficient and vibrational relaxation time) is also conducted. It is found that the 2(sigma) - uncertainty for heating on Uranus entry is no more than 2.1%, assuming an equilibrium catalytic wall, with the uncertainty being determined primarily by diffusion and H(sub 2) recombination rate within the boundary layer. However, if the wall is assumed to be partially or non-catalytic, this uncertainty may increase to as large as 18%. The catalytic wall model can contribute over 3x change in heat flux and a 20% variation in film coefficient. Therefore, coupled material response/fluid dynamic models are recommended for this problem. It was also found that much of this variability is artificially suppressed when a constant Schmidt number approach is implemented. Because the boundary layer is reacting, it is necessary to employ self-consistent effective binary diffusion to obtain a correct thermal transport solution. For Saturn entries, the 2(sigma) - uncertainty for convective heating was less than 3.7%. The major uncertainty driver was dependent on shock temperature/velocity, changing from boundary layer thermal conductivity to diffusivity and then to shock layer ionization rate as velocity increases. While radiative heating for Uranus entry was negligible, the nominal solution for Saturn computed up to 20% radiative heating at the highest velocity examined. The radiative heating followed a non-normal distribution, with up to a 3x variation in magnitude. This uncertainty is driven by the H(sub 2) dissociation rate, as H(sub 2) that persists in the hot non-equilibrium zone contributes significantly to radiation.

Palmer, Grant↗
NASA NTRS1991

Numerical and experimental investigations of rarefied nozzle and plume flows of nitrogen

Numerical and experimental investigations are performed for the rarefied flow of nitrogen through a small nozzle which is expanded into near-vacuum conditions. Two different numerical studies are undertaken: the first employs a continuum approach in solving the Navier-Stokes equations, and the second employs a particle approach through use of the direct simulation Monte Carlo method (DSMC). The experimental investigation concerns the measurement of pressure, using a Pitot tube, in the nozzle exit plane and near-field of the plume. Comparison of the experimental and numerical data at the nozzle exit reveals that the DSMC technique provides the more accurate description of the expanding flow. It is discovered that the DSMC solutions are quite sensitive to the model employed to simulate the interaction between the gas and the nozzle wall surface. It is concluded that the simplistic fully diffuse model is quite satisfactory for the present application.

Boyd, Iain D.↗
NASA NTRS1985

Numerical Descriptions of Cosmic-Ray Transport

The behavior of energetic particles in the solar system is described by a well known Fokker-Planck equation. Although analytic methods yield insight into the nature of its solutions, especially in the diffusion regime, calculations that go beyond diffusion are very complicated. The reliability of these calculations is of concern, because numerical methods are notorious for their errors and artifacts. The well known Milne problem of classical transport theory was analyzed with the aid of three different numerical methods. These are: (1) The method of eigenfunctions in which the distribution function is approximated by a sum of eigenfunctions of the scattering operator, (2) Numerical solutions of a finite difference aquation; and (3) Direct simulation of the scattering and streaming of individual particles with the aid of Monte Carlo methods.

Earl, J. A.↗
NASA NTRS1980

Solar wind sputtering effects in the atmospheres of Mars and Venus

It is found through an investigation, combining Monte Carlo simulations and analytical techniques, of the direct collisional interaction of an energetic particle flux with the neutral components of a planetary atmosphere, that solar wind sputtering could provide an important exospheric mass sink on both Mars and Venus under appropriate conditions. The computed rates of helium loss in the Venusian atmosphere and of carbon, nitrogen and oxygen in the Martian atmosphere imply that sputtering would have a significant effect on the noble gas budget of Venus and dominate the chemical and photochemical loss processes of Mars. Because of diffusive separation of lighter elements and isotopes, and because the gravitational binding energy is proportional to the mass, the erosion process preferentially removes the lighter components of the atmosphere. Solar wind sputtering could therefore compete with other erosion mechanisms in generating substantial fractionation effects.

Watson, C. C.↗
NASA NTRS1987

A dust scattered halo in starburst galaxy M82?

The source of the halo about M82 has been under discussion for several years. One explanation for it is the dust model of Solinger, Morrison and Markert in which they propose a diffuse cloud of dust through the M81 group, with M82 traveling through the group holding a denser cloud of dust around it. The feasibility of the dust theory is examined in the X-ray range, using the halo in the X-ray image of M82 taken by the Einstein Observatory. To this end the X-ray cross section for dust is presented, along with the single scattered image of an X-ray source surrounded by a dust cloud; multiply scattered images were simulated with a Monte Carlo program; profiles of the halo along the major and minor axes of M82 are presented. Also presented is an accounting for line spectrographs of M82 that show unusual splitting, using the dust model. The final model proposed for the X-ray image requires dust of radius 50 to 300 A, with density on the order of 10 to the -7th power cu cm, out to a distance of about 9 kpc for some regions.

Rohan, Michael↗
NASA NTRS1983

A potential-energy scaling model to simulate the initial stages of thin-film growth

A solid on solid (SOS) Monte Carlo computer simulation employing a potential energy scaling technique was used to model the initial stages of thin film growth. The model monitors variations in the vertical interaction potential that occur due to the arrival or departure of selected adatoms or impurities at all sites in the 400 sq. ft. array. Boltzmann ordered statistics are used to simulate fluctuations in vibrational energy at each site in the array, and the resulting site energy is compared with threshold levels of possible atomic events. In addition to adsorption, desorption, and surface migration, adatom incorporation and diffusion of a substrate atom to the surface are also included. The lateral interaction of nearest, second nearest, and third nearest neighbors is also considered. A series of computer experiments are conducted to illustrate the behavior of the model.

Heinbockel, J. H.↗
NASA NTRS1973

Path-length distributions of photons diffusely reflected from a semi-infinite atmosphere.

Knowledge of the probability distribution of photon path lengths for diffuse reflection in the continuum from a planetary atmosphere greatly simplifies computation of the shape and strength of absorption lines formed under conditions of multiple scattering. Graphs of this probability distribution are presented for diffuse reflection from a semiinfinite conservative atmosphere for both isotropic and forward-directed scattering, and for several angles of incidence and reflection. Graphs showing distributions for the case of absorption in the continuum are also discussed. The distributions are obtained from an analytic expression for the path-length distribution for nth-order scattering and a tabulation of nth-order intensities. It is found that te shape of the probability distribution depends markedly on the degree of forward scattering due to the phase function. The results are compared with those obtained by a Monte Carlo computation of the sort which might be used if tabulated nth-order intensities were not available.

Appleby, J. F.↗
NASA NTRS1994

The effect of wave-particle interactions on the polar wind: Preliminary results

A Monte Carlo simulation was developed in order to study the effects of wave-particle interactions (WPI) on the plasma outflow in the polar wind. The simulation also considered the other mechanisms included in the `classical' polar wind studies such as gravity, the polarization electrostatic field, and the divergence of geomagnetic field lines. Although the plasma consisted of electrons, H(+) and O(+) ions, we emphasized the behavior of H(+) in this preliminary study. The ion distribution function, as well as the profiles of its moments (density, drift velocity, temperatures, etc.) were found for different levels of WPI, that is, for different values of normalized diffusion rates in the velocity space (normalized D(sub perpendicular) H(+)). We found that as the WPI strength increases: (1) the ion drift velocity increases and its density decreases; (2) the perpendicular temperature T(sub perpendicular) (H(+)) increases; (3) the parallel temperature T(sub parallel) (H(+)) first decreases and then increases due to the balance between the parallel adiabatic cooling and the transfer of the energy from the perpendicular to the parallel direction; and (4) the temperature anisotropy ((T(sub parallel) (H(+))/T(sub perpendicular) (H(+)) is reduced and even reversed in some cases. For strong WPI (normalized D(sub perpendicular) (H(+) much greater than 1), the ion distribution function shows weak conic features at high altitudes.

Barakat, A. R.↗
NASA NTRS1998

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

VanOverbeke, Thomas J.↗
NASA NTRS2001

Transport Phenomena of Off-Axis Sputtering Deposition

Various high quality epitaxial films, especially oxides, have been synthesized using off-axis sputtering deposition. In this presentation, we report the experiment results of ZnO films grown by the off-axis sputtering deposition. Films were synthesized in temperatures ranged from room temperature to 600 C, and pressures from 5 mTorr to 150 mTorr. Film growth rate was measured by surface profilometer, ellipsometer, and wavelength dispersive spectrometry. Due to the collisions between the sputtered species and the residue gases, the kinetic energy of species was reduced and the transport of depositing species changed from a ballistic movement for low pressure to a diffusive drift for high pressure in which the transport species were almost thermalized. The measurements show an increase of growth rates along the gravity vector when the Knudsen number of transport species is less than 0.05, which suggests that gravity affected the transport characterization in off-axis sputtering deposition. Because the product of pressure (p) and travel distance (d) of sputtered species, pd, exceeds several mTorr-cm during film deposition, the classical simulations for sputtering process in high vacuum system may not be applied. Based on these experimental measurements, a transport process of the off-axis sputtering deposition is proposed. Several models including the Monte Carlo method and gravity-driven flow dynamics simulation will be discussed.

Zhu, S.↗
NASA NTRS2001

Transport Phenomena of Off-Axis Sputtering Deposition

Various high quality epitaxial films, especially oxides, have been synthesized using off-axis sputtering deposition. In this presentation, we report the experiment results of ZnO films grown by the off-axis sputtering deposition. Films were synthesized in temperatures ranged from room temperature to 600 C, and pressures from 5 mTorr to 150 mTorr. Film growth rate was measured by surface profilometer, ellipsometer, and wavelength dispersive spectrometry. Due to the collisions between the sputtered species and the residue gases, the kinetic energy of species was reduced and the transport of depositing species changed from a ballistic movement for low pressure to a diffuse drift for high pressure in which the transport species were almost thermalized. The measurements show an increase of growth rates along the gravity vector when the Knodson (Knudsen??) number of transport species is less than 0.05, which suggests that gravity affected the transport characterization in off-axis sputtering deposition. Because the product of pressure (p) and travel distance (d) of sputtered species, p exceeds several mTorr-cm during film deposition, the classical simulations for sputtering process in high vacuum system may not be applied. Based on these experimental measurements, a transport process of the off-axis sputtering deposition is proposed. Several methods including the Monte Carlo method and gravity-driven flow dynamics simulation will be discussed.

Zhu, S.↗
NASA NTRS2014

A Radiation Chemistry Code Based on the Green's Function of the Diffusion Equation

Stochastic radiation track structure codes are of great interest for space radiation studies and hadron therapy in medicine. These codes are used for a many purposes, notably for microdosimetry and DNA damage studies. In the last two decades, they were also used with the Independent Reaction Times (IRT) method in the simulation of chemical reactions, to calculate the yield of various radiolytic species produced during the radiolysis of water and in chemical dosimeters. Recently, we have developed a Green's function based code to simulate reversible chemical reactions with an intermediate state, which yielded results in excellent agreement with those obtained by using the IRT method. This code was also used to simulate and the interaction of particles with membrane receptors. We are in the process of including this program for use with the Monte‐Carlo track structure code Relativistic Ion Tracks (RITRACKS). This recent addition should greatly expand the capabilities of RITRACKS, notably to simulate DNA damage by both the direct and indirect effect.

Plante, Ianik↗
NASA NTRS1985

Large-cell Monte Carlo renormalization of irreversible growth processes

Monte Carlo sampling is applied to a recently formulated direct-cell renormalization method for irreversible, disorderly growth processes. Large-cell Monte Carlo renormalization is carried out for various nonequilibrium problems based on the formulation dealing with relative probabilities. Specifically, the method is demonstrated by application to the 'true' self-avoiding walk and the Eden model of growing animals for d = 2, 3, and 4 and to the invasion percolation problem for d = 2 and 3. The results are asymptotically in agreement with expectations; however, unexpected complications arise, suggesting the possibility of crossovers, and in any case, demonstrating the danger of using small cells alone, because of the very slow convergence as the cell size b is extrapolated to infinity. The difficulty of applying the present method to the diffusion-limited-aggregation model, is commented on.

Nakanishi, H.↗
NASA NTRS1991

The estimation of background production by cosmic rays in high-energy gamma ray telescopes

A calculational method of estimating instrumental background in high-energy gamma-ray telescopes, using the hadronic Monte Carlo code FLUKA87, is presented. The method is applied to the SAS-2 and EGRET telescope designs and is also used to explore the level of background to be expected for alternative configurations of the proposed GRITS telescope, which adapts the external fuel tank of a Space Shuttle as a gamma-ray telescope with a very large collecting area. The background produced in proton-beam tests of EGRET is much less than the predicted level. This discrepancy appears to be due to the FLUKA87 inability to transport evaporation nucleons. It is predicted that the background in EGRET will be no more than 4-10 percent of the extragalactic diffuse gamma radiation.

Edwards, H. L.↗
NASA NTRS2012

Disentangling AGN and Star Formation in Soft X-Rays

We have explored the interplay of star formation and active galactic nucleus (AGN) activity in soft X-rays (0.5-2 keV) in two samples of Seyfert 2 galaxies (Sy2s). Using a combination of low-resolution CCD spectra from Chandra and XMM-Newton, we modeled the soft emission of 34 Sy2s using power-law and thermal models. For the 11 sources with high signal-to-noise Chandra imaging of the diffuse host galaxy emission, we estimate the luminosity due to star formation by removing the AGN, fitting the residual emission. The AGN and star formation contributions to the soft X-ray luminosity (i.e., L(sub x,AGN) and L(sub x,SF)) for the remaining 24 Sy2s were estimated from the power-law and thermal luminosities derived from spectral fitting. These luminosities were scaled based on a template derived from XSINGS analysis of normal star-forming galaxies. To account for errors in the luminosities derived from spectral fitting and the spread in the scaling factor, we estimated L(sub x,AGN) and L(sub x,SF))from Monte Carlo simulations. These simulated luminosities agree with L(sub x,AGN) and L(sub x,SF) derived from Chandra imaging analysis within a 3sigma confidence level. Using the infrared [Ne ii]12.8 micron and [O iv]26 micron lines as a proxy of star formation and AGN activity, respectively, we independently disentangle the contributions of these two processes to the total soft X-ray emission. This decomposition generally agrees with L(sub x,SF) and L(sub x,AGN) at the 3 sigma level. In the absence of resolvable nuclear emission, our decomposition method provides a reasonable estimate of emission due to star formation in galaxies hosting type 2 AGNs.

Disentangling↗
NASA NTRS1986

Direct simulation of typical AOTV entry flows

The direct simulation Monte Carlo method is applied to two flows that are relevant to the projected aeroassisted orbital transfer vehicle. The first is the axially symmetric flow about a 70 degree blunted cone shield with a cylindrical afterbody, and the second is the one-dimensional flow along the stagnation streamline of the first flow. The shock standoff distance for the stagnation streamline study is obtained from the full axially symmetric calculation. Advantage is then taken of the speed of the one-dimensional calculation to study the sensitivity of the stagnation point heat transfer to a number of physical variables. These include the chemical reaction rate coefficients, the surface catalycity, the electric field associated with ambipolar diffusion, and the nonequilibrium thermal radiation from the shock front.

Bird, G. A.↗