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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Effect of magnetic helicity upon rectilinear propagation of charged particles in random magnetic fields

When charged particles spiral along a large constant magnetic field, their trajectories are scattered by any random field components that are superposed on the guiding field. If the random field configuration embodies helicity, the scattering is asymmetrical with respect to a plane perpendicular to the guiding field, for particles moving into the forward hemisphere are scattered at different rates from those moving into the backward hemisphere. This asymmetry gives rise to new terms in the transport equations that describe propagation of charged particles. Helicity has virtually no impact on qualitative features of the diffusive mode of propagation. However, characteristic velocities of the coherent modes that appear after a highly anisotropic injection exhibit an asymmetry related to helicity. Explicit formulas, which embody the effects of helicity, are given for the anisotropies, the coefficient diffusion, and the coherent velocities. Predictions derived from these expressions are in good agreement with Monte Carlo simulations of particle transport, but the simulations reveal certain phenomena whose explanation calls for further analytical work.

Earl, James A.↗
NASA NTRS2013

Direct Simulation of Extinction in a Slab of Spherical Particles

The exact multiple sphere superposition method is used to calculate the coherent and incoherent contributions to the ensemble-averaged electric field amplitude and Poynting vector in systems of randomly positioned nonabsorbing spherical particles. The target systems consist of cylindrical volumes, with radius several times larger than length, containing spheres with positional configurations generated by a Monte Carlo sampling method. Spatially dependent values for coherent electric field amplitude, coherent energy flux, and diffuse energy flux, are calculated by averaging of exact local field and flux values over multiple configurations and over spatially independent directions for fixed target geometry, sphere properties, and sphere volume fraction. Our results reveal exponential attenuation of the coherent field and the coherent energy flux inside the particulate layer and thereby further corroborate the general methodology of the microphysical radiative transfer theory. An effective medium model based on plane wave transmission and reflection by a plane layer is used to model the dependence of the coherent electric field on particle packing density. The effective attenuation coefficient of the random medium, computed from the direct simulations, is found to agree closely with effective medium theories and with measurements. In addition, the simulation results reveal the presence of a counter-propagating component to the coherent field, which arises due to the internal reflection of the main coherent field component by the target boundary. The characteristics of the diffuse flux are compared to, and found to be consistent with, a model based on the diffusion approximation of the radiative transfer theory.

spheres↗
NASA NTRS2000

A Validation Summary of the NCC Turbulent Reacting/non-reacting Spray Computations

This pper provides a validation summary of the spray computations performed as a part of the NCC (National Combustion Code) development activity. NCC is being developed with the aim of advancing the current prediction tools used in the design of advanced technology combustors based on the multidimensional computational methods. The solution procedure combines the novelty of the application of the scalar Monte Carlo PDF (Probability Density Function) method to the modeling of turbulent spray flames with the ability to perform the computations on unstructured grids with parallel computing. The calculation procedure was applied to predict the flow properties of three different spray cases. One is a nonswirling unconfined reacting spray, the second is a nonswirling unconfined nonreacting spray, and the third is a confined swirl-stabilized spray flame. The comparisons involving both gas-phase and droplet velocities, droplet size distributions, and gas-phase temperatures show reasonable agreement with the available experimental data. The comparisons involve both the results obtained from the use of the Monte Carlo PDF method as well as those obtained from the conventional computational fluid dynamics (CFD) solution. Detailed comparisons in the case of a reacting nonswirling spray clearly highlight the importance of chemistry/turbulence interactions in the modeling of reacting sprays. The results from the PDF and non-PDF methods were found to be markedly different and the PDF solution is closer to the reported experimental data. The PDF computations predict that most of the combustion occurs in a predominantly diffusion-flame environment. However, the non-PDF solution predicts incorrectly that the combustion occurs in a predominantly vaporization-controlled regime. The Monte Carlo temperature distribution shows that the functional form of the PDF for the temperature fluctuations varies substantially from point to point. The results also bring to the fore some of the deficiencies associated with the use of assumed-shape PDF methods in spray computations.

Raju, M. S.↗
NASA NTRS1991

A study of hydrogen diffusion flames using PDF turbulence model

The application of probability density function (pdf) turbulence models is addressed in this work. For the purpose of accurate prediction of turbulent combustion, an algorithm that combines a conventional CFD flow solver with the Monte Carlo simulation of the pdf evolution equation has been developed. The algorithm has been validated using experimental data for a heated turbulent plane jet. The study of H2-F2 diffusion flames has been carried out using this algorithm. Numerical results compared favorably with experimental data. The computuations 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.↗
NASA NTRS1974

Some new results on electron transport in the atmosphere

The penetration, diffusion and slowing down of electrons in a semiinfinite air medium has been studied by the Monte Carlo method. The results are applicable in the atmosphere at altitudes up to 300 km. Most of the results pertain to monoenergetic electron beams, with energies between 2 keV and 2 MeV, injected into the atmosphere at a height of 300 km, either vertically downwards or with a pitch-angle distribution isotropic over the downward hemisphere. Some results were also obtained for various initial pitch angles between 0 and 90 deg. Information has been generated concerning the backscattering of electrons from the atmosphere, the altitude dependence of energy deposition by electrons and by secondary bremsstrahlung, and the evolution of electron flux spectra as function of the atmospheric depth.

Berger, M. J.↗
NASA NTRS1976

Impulsive solar X-ray bursts. 4: Polarization, directivity and spectrum of the reflected and total bremsstrahlung radiation from a beam of electrons directed toward the photosphere

A Monte Carlo method is described for evaluation of the spectrum, directivity and polarization of X-rays diffusely reflected from stellar photospheres. the accuracy of the technique is evaluated through comparison with analytic results. Using the characteristics of the incident X-rays of the model for solar X-ray flares, the spectrum, directivity and polarization of the reflected and the total X-ray fluxes are evaluated. The results are compared with observations.

Langer, S. H.↗
NASA NTRS1983

Modeling of Fluctuating Mass Flux in Variable Density Flows

The approach solves for both Reynolds and Favre averaged quantities and calculates the scalar pdf. Turbulent models used to close the governing equations are formulated to account for complex mixing and variable density effects. In addition, turbulent mass diffusivities are not assumed to be in constant proportion to turbulent momentum diffusivities. The governing equations are solved by a combination of finite-difference technique and Monte-Carlo simulation. Some preliminary results on simple variable density shear flows are presented. The differences between these results and those obtained using conventional models are discussed.

So, R. M. C.↗
NASA NTRS1971

Heat transfer and forces on concave surfaces in free molecule flow.

A Monte Carlo modeling technique is described for mathematically simulating free molecular flows over a concave spherical surface and a concave cylindrical surface of finite length. The half-angle of the surfaces may vary from 0 to 90 degrees, and the incident flow may have an arbitrary speed ratio and an arbitrary angle of attack. Partial diffuse reflection and imperfect energy accommodation for molecules colliding with the surfaces are also considered. Results of heat transfer, drag and lift coefficients are presented for a variety of flow conditions. The present Monte Carlo results are shown to be in very good agreement with certain available theoretical solutions.

Fan, C.↗
NASA NTRS1977

Impulsive solar X-ray bursts. III - Polarization, directivity, and spectrum of the reflected and total bremsstrahlung radiation from a beam of electrons directed toward the photosphere

The paper presents the spectrum, directivity, and state of polarization of the bremsstrahlung radiation expected from a beam of high-energy electrons spiraling along radial magnetic field lines toward the photosphere. A Monte Carlo method is then described for evaluation of the spectrum, directivity, and polarization of X-rays diffusely reflected from stellar photospheres. The accuracy of the technique is evaluated through comparison with analytic results. The calculated characteristics of the incident X-rays are used to evaluate the spectrum, directivity, and polarization of the reflected and total X-ray fluxes. The results are compared with observations.

Langer, S. H.↗
NASA NTRS1981

A numerical model study of long-term planetary wave predictability

A numerical model was developed to study planetary wave behavior over periods up to 4 months under the influences of stochastic and diffusive forcing by synoptic scale motions and of sea surface temperature anomalies. The model is suitable for Monte Carlo experiments which provide ensemble mean and variance predictions and permit the response to heating anomalies to be partly separated from random variability. The model incorporates the linear-balanced dynamical equations for wind at two tropospheric levels and temperature at one level. The model atmosphere is forced by semiempirical January heating linearly dependent on local air temperature. Topography and the lower boundary layer are represented by vertical motion and drag functions, respectively, of the lower level wind. The linear phase speeds of Rossby waves are modified to reflect the effect of the stratosphere.

Wobus, R. L.↗
NASA NTRS1995

A case study of view-factor rectification procedures for diffuse-gray radiation enclosure computations

The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.

Taylor, Robert P.↗
NASA NTRS1972

Some new results on electron transport in the atmosphere

The penetration, diffusion and slowing down of electrons in a semi-infinite air medium has been studied by the Monte Carlo method. The results are applicable to the atmosphere at altitudes up to 300 km. Most of the results pertain to monoenergetic electron beams injected into the atmosphere at a height of 300 km, either vertically downwards or with a pitch-angle distribution isotropic over the downward hemisphere. Some results were also obtained for various initial pitch angles between 0 deg and 90 deg. Information has been generated concerning the following topics: (1) the backscattering of electrons from the atmosphere, expressed in terms of backscattering coefficients, angular distributions and energy spectra of reflected electrons, for incident energies T(o) between 2 keV and 2 MeV; (2) energy deposition by electrons as a function of the altitude, down to 80 km, for T(o) between 2 keV and 2 MeV; (3) the corresponding energy depostion by electron-produced bremsstrahlung, down to 30 km; (4) the evolution of the electron flux spectrum as function of the atmospheric depth, for T(o) between 2 keV and 20 keV. Energy deposition results are given for incident electron beams with exponential and power-exponential spectra.

Berger, M. J.↗
NASA NTRS1992

Laboratory scaled simulation of lidar cloud sounding experiments

The results of lidar measurements carried out on laboratory scale models of clouds are presented. Measurements on laboratory scale models are important since one has the knowledge of the relevant parameters of the diffusing medium, such as: scattering and absorption coefficients, phase function, homogeneity, shape, etc. Knowledge of these parameters enables one to use the results to test the reliability of theoretical and numerical investigations. To obtain a laboratory scaled model of a lidar system sounding a cloud, it is necessary to scale down all the geometrical quantities by the same factor to reduce distances of the order of kilometers to the order of meters, keeping the size and the optical depth of the diffusers unchanged. If a time resolution of the order of nanoseconds is necessary for a lidar sounding actual clouds, the corresponding time resolution for the laboratory model should be of the order of picoseconds. It is possible to obtain this resolution by using picosecond laser systems and fast electrooptical detectors like the streak camera. The results of the laboratory measurements showed that the multiple scattering effect strongly depends on the size of the diffusers, as well as on the concentration. The experimental results were compared with the numerical results of a Monte Carlo code. A generally good agreement was obtained.

Zaccanti, G.↗
NASA NTRS2012

Taking the Pulse of PyroCumulus Clouds

Forest fires can burn large areas, but can also inject smoke into the upper troposphere/lower stratosphere (UT/LS), where stakes are even higher for climate, because emissions tend to have a longer lifetime, and can produce significant regional and even global climate effects, as is the case with some volcanoes. Large forest fires are now believed to be more common in summer, especially in the boreal regions, where pyrocumulus (pyroCu), and occasionally pyrocumuionimbus (pyroCb) clouds are formed, which can transport emissions into the UT/LS. A major difficulty in developing realistic fire plume models is the lack of observational data within fire plumes that resolves structure at a few 100 m scales, which can be used to validate these models. Here, we report detailed airborne radiation measurements within strong pyroCu taken over boreal forest fires in Saskatchewan, Canada during the Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) summer field campaign in 2008. We find that the angular distribution of radiance within the pyroCu is closely related to the diffusion domain in water clouds and can be described by very similar simple cosine functions. We demonstrate with Monte Carlo simulations that radiation transport in pyroCu is inherently a 3D phenomenon and must account for particle absorption. However, the simple cosine function promises to offer an easy solution for climate models. The presence of a prominent smoke core, defined by strong extinction in the UV, VIS and NIR, suggests that the core might be an important pathway for emission transport to the upper troposphere and lower stratosphere. We speculate that this plume injection core is generated and sustained by complex processes not yet well understood, but not necessarily related directly to the intense fires that originally initiated the plume rise.

Gatebe, C. K.↗
NASA NTRS1976

The partially averaged field approach to cosmic ray diffusion

The kinetic equation for particles interacting with turbulent fluctuations is derived by a new nonlinear technique which successfully corrects the difficulties associated with quasilinear theory. In this new method the effects of the fluctuations are evaluated along particle orbits which themselves include the effects of a statistically averaged subset of the possible configurations of the turbulence. The new method is illustrated by calculating the pitch angle diffusion coefficient D sub Mu Mu for particles interacting with slab model magnetic turbulence, i.e., magnetic fluctuations linearly polarized transverse to a mean magnetic field. Results are compared with those of quasilinear theory and also with those of Monte Carlo calculations. The major effect of the nonlinear treatment in this illustration is the determination of D sub Mu Mu in the vicinity of 90 deg pitch angles where quasilinear theory breaks down. The spatial diffusion coefficient parallel to a mean magnetic field is evaluated using D sub Mu Mu as calculated by this technique. It is argued that the partially averaged field method is not limited to small amplitude fluctuating fields and is hence not a perturbation theory.

Jones, F. C.↗
NASA NTRS2016

Probing Acceleration and Turbulence at Relativistic Shocks in Blazar Jets

Diffusive shock acceleration (DSA) at relativistic shocks is widely thought to be an important acceleration mechanism in various astrophysical jet sources, including radio-loud active galactic nuclei such as blazars. Such acceleration can produce the non-thermal particles that emit the broad-band continuum radiation that is detected from extragalactic jets. An important recent development for blazar science is the ability of Fermi-Large Area Telescope spectroscopy to pin down the shape of the distribution of the underlying non-thermal particle population. This paper highlights how multiwavelength spectra spanning optical to X-ray to gamma-ray bands can be used to probe diffusive acceleration in relativistic, oblique, magnetohydrodynamic (MHD) shocks in blazar jets. Diagnostics on the MHD turbulence near such shocks are obtained using thermal and non-thermal particle distributions resulting from detailed Monte Carlo simulations of DSA. These probes are afforded by the characteristic property that the synchrotron vF(sub v), peak energy does not appear in the gamma-ray band above 100 MeV. We investigate self-consistently the radiative synchrotron and inverse Compton signatures of the simulated particle distributions. Important constraints on the diffusive mean free paths of electrons, and the level of electromagnetic field turbulence are identified for three different case study blazars, Mrk 501, BL Lacertae and AO 0235+164. The X-ray excess of AO 0235+164 in a flare state can be modelled as the signature of bulk Compton scattering of external radiation fields, thereby tightly constraining the energy-dependence of the diffusion coefficient for electrons. The concomitant interpretations that turbulence levels decline with remoteness from jet shocks, and the probable significant role for non-gyroresonant diffusion, are posited.

Baring, Matthew G.↗
NASA NTRS1976

X-ray reverberations and the giant X-ray bursts

It is shown that the shape and spectral evolution of the giant X-ray bursts from the source 3U 1820-30 can be explained by Compton scattering of a short, intense X-ray pulse in a cloud surrounding the source. Pulse shapes due to Thomson scattering of an X-ray burst in an electron cloud were calculated for the (1) optically thin case on the assumption of one scattering per photon, (2) intermediate case with optical depth of about unity, and (3) optically thick case where the process is regarded as diffusion of photons through a uniform sphere. For the intermediate case, the effects of the reverberation were determined explicitly by Monte Carlo calculation. For an optical depth of 3, square pulse duration of 2 sec, characteristic cloud radius of 70,000 km, characteristic cloud density of 4 times 10 to the 14th per cu cm, and temperature of 5-30 keV, the calculations give a reasonably accurate description of X-ray bursts from 3U 1820-30. The scattering model does not imply the existence of a supermassive, central black hole.

Canizares, C. R.↗
NASA NTRS1987

Particle acceleration at shocks - A Monte Carlo method

A Monte Carlo method is presented for the problem of acceleration of test particles at relativistic shocks. The particles are assumed to diffuse in pitch angle as a result of scattering off magnetic irregularities frozen into the fluid. Several tests are performed using the analytic results available for both relativistic and nonrelativistic shock speeds. The acceleration at relativistic shocks under the influence of radiation losses is investigated, including the effects of a momentum dependence in the diffusion coefficient. The results demonstrate the usefulness of the technique in those situations in which the diffusion approximation cannot be employed, such as when relativistic bulk motion is considered, when particles are permitted to escape at the boundaries, and when the effects of the finite length of the particle mean free path are important.

Kirk, J. G.↗