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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS2024

Spectral Correlation in MODIS Water-Leaving Reflectance Retrieval Uncertainty

Spectral remote sensing reflectance, Rrs(λ) (sr−1), is the fundamental quantity used to derive a host of bio-optical and biogeochemical properties of the water column from satellite ocean color measurements. Estimation of uncertainty in those derived geophysical products is therefore dependent on knowledge of the uncertainty in satellite-retrieved R rs . Furthermore, since the associated algorithms require R rs at multiple spectral bands, the spectral (i.e., band-to-band)error covariance in R rs is needed to accurately estimate the uncertainty in those derived properties. This study establishes a derivative-based approach for propagating instrument random noise, instrument systematic uncertainty, and forward model uncertainty into R rs as retrieved using NASA’s multiple-scattering epsilon (MSEPS) atmospheric correction algorithm, to generate pixel-level error covariance in R rs . The approach is applied to measurements from Moderate Resolution Imaging Spectroradiometer (MODIS) on the Aqua satellite and verified using Monte Carlo (MC) analysis. We also make use of this full spectral error covariance in R rs to calculate uncertainty in phytoplankton pigment chlorophyll-a concentration (chl a , mg/m 3 ) and diffuse attenuation coefficient of downwelling irradiance at 490 nm (K d (490), m -1 ). Accounting for the error covariance in R rs generally reduces the estimated relative uncertainty in chl a by ∼1-2% (absolute value) in waters with chl a < 0.25 mg/m 3 where the color index (CI) algorithm is used. The reduction is ∼5-10% in waters with chl a > 0.35 mg/m 3 where the blue-green ratio (OCX) algorithm is used. Such reduction can be higher than 30% in some regions. For K d (490), the reduction by error covariance is generally ∼2%, but can be higher than 20% in some regions. The error covariance in R rs is further verified through forward-calculating chl a from MODIS-retrieved and in situ R rs and comparing estimated uncertainty with observed differences. An 8-day global composite of propagated uncertainty shows that the goal of 35% uncertainty in chl a can be achieved over deep ocean waters (chl a ≤ 0.1 mg/m3). While the derivative-based approach generates reasonable error covariance in R rs some assumptions should be updated as our knowledge improves. These include the inter-band error correlation in top-of-atmosphere reflectance, and uncertainties in the calibration of MODIS 869 nm band, in ancillary data, and in the in situ data used for system vicarious calibration.

Ocean color↗

NASA NTRS1983

Dynamical evolution of the Oort cometary cloud

The dynamical evolution of comets in the Oort cloud under the influence of stellar perturbations has been modeled using Monte Carlo techniques. It is shown that the cloud has been depleted over the history of the solar system. Comets are lost from the cloud by direct ejection due to close stellar encounters, diffusion of aphelia to distances beyond the sun's sphere of influence, or diffusion of perihelia into the planetary region where Jupiter and Saturn perturbations either eject them on hyperbolic trajectories or capture them to short-period orbits. The population of the cloud is estimated to be 1.0 - 1.5 x 10 to the 12th comets and the total mass is on the order of 1.9 earth masses. In addition to random passing stars, less frequent encounters with giant molecular clouds may play a significant role in randomizing the orbits of comets in the cloud and reducing the effective radius of the sun's sphere of influence.

Weissman, P. R.↗

NASA NTRS1978

Radiative transfer calculated from a Markov chain formalism

The theory of Markov chains is used to formulate the radiative transport problem in a general way by modeling the successive interactions of a photon as a stochastic process. Under the minimal requirement that the stochastic process is a Markov chain, the determination of the diffuse reflection or transmission from a scattering atmosphere is equivalent to the solution of a system of linear equations. This treatment is mathematically equivalent to, and thus has many of the advantages of, Monte Carlo methods, but can be considerably more rapid than Monte Carlo algorithms for numerical calculations in particular applications. We have verified the speed and accuracy of this formalism for the standard problem of finding the intensity of scattered light from a homogeneous plane-parallel atmosphere with an arbitrary phase function for scattering. Accurate results over a wide range of parameters were obtained with computation times comparable to those of a standard 'doubling' routine. The generality of this formalism thus allows fast, direct solutions to problems that were previously soluble only by Monte Carlo methods. Some comparisons are made with respect to integral equation methods.

Esposito, L. W.↗

NASA NTRS1979

DUMAND Summer Workshop, University of California, La Jolla, Calif., July 24-September 2, 1978, Proceedings. Volume 2 - UHE interactions, neutrino astronomy

The volume covers categories on inelastic neutrino scattering and the W-boson, and other ultra-high-energy processes, on pulsars, quasars and galactic nuclei, as well as other point sources and constants from gamma ray astronomy. Individual subjects include weak intermediate vector bosons and DUMAND, the Monte Carlo simulation of inelastic neutrino scattering in DUMAND, and Higgs boson production by very high-energy neutrinos. The observability of the neutrino flux from the inner region of the galactic disk, the diffuse fluxes of high-energy neutrinos, as well as the significance of gamma ray observations for neutrino astronomy are also among the topics covered.

Roberts, A.↗

NASA NTRS1995

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

NASA NTRS1992

DSMC analysis of species separation in rarefied nozzle flows

The direct-simulation Monte Carlo method has been used to investigate the behavior of a small amount of a harmful species in the plume and the backflow region of nuclear thermal propulsion rockets. Species separation due to pressure diffusion and nonequilibrium effects due to rapid expansion into a surrounding low-density environment are the most important factors in this type of flow. It is shown that a relatively large amount of the lighter species is scattered into the backflow region and the heavier species becomes negligible in this region due to the extreme separation between species. It is also shown that the type of molecular interaction between the species can have a substantial effect on separation of the species.

Chung, Chan-Hong↗

NASA NTRS2014

Solar Wind Implantation into Lunar Regolith: Hydrogen Retention in a Surface with Defects

Solar wind protons are implanted directly into the top 100 nm of the lunar near-surface region, but can either quickly diffuse out of the surface or be retained, depending upon surface temperature and the activation energy, U, associated with the implantation site. In this work, we explore the distribution of activation energies upon implantation and the associated hydrogen-retention times; this for comparison with recent observation of OH on the lunar surface. We apply a Monte Carlo approach: for simulated solar wind protons at a given local time, we assume a distribution of U values with a central peak, U(sub c) and width, U(sub w), and derive the fraction retained for long periods in the near-surface. We find that surfaces characterized by a distribution with predominantly large values of U (greater than 1 eV) like that expected at defect sites will retain implanted H (to likely form OH). Surfaces with the distribution predominantly at small values of U (less than 0.2 eV) will quickly diffuse away implanted H. However, surfaces with a large portion of activation energies between 0.3 eV less than U less than 0.9 eV will tend to be H-retentive in cool conditions but transform into H-emissive surfaces when warmed (as when the surface rotates into local noon). These mid-range activation energies give rise to a diurnal effect with diffusive loss of H at noontime.

regolith↗

NASA NTRS1973

A new approach to cosmic ray diffusion theory

An approach is presented for deriving a diffusion equation for charged particles in a static, random magnetic field. The approach differs from the usual, quasi-linear one, in that particle orbits in the average field are replaced by particle orbits in a partially averaged field. In this way the fluctuating component of the field significantly modifies the particle orbits in those cases where the orbits in the average field are unrealistic. The method permits the calculation of a finite value for the pitch angle diffusion coefficient for particles with a pitch angle of 90 rather than the divergent or ambiguous results obtained by quasi-linear theories. Results of the approach are compared with results of computer simulations using Monte Carlo techniques.

Jones, F. C.↗

NASA NTRS2013

A Green's Function Approach to Simulate DNA Damage by the Indirect Effect

The DNA damage is of fundamental importance in the understanding of the effects of ionizing radiation. DNA is damaged by the direct effect of radiation (e.g. direct ionization) and by indirect effect (e.g. damage by.OH radicals created by the radiolysis of water). Despite years of research, many questions on the DNA damage by ionizing radiation remains. In the recent years, the Green's functions of the diffusion equation (GFDE) have been used extensively in biochemistry [1], notably to simulate biochemical networks in time and space [2]. In our future work on DNA damage, we wish to use an approach based on the GFDE to refine existing models on the indirect effect of ionizing radiation on DNA. To do so, we will use the code RITRACKS [3] developed at the NASA Johnson Space Center to simulate the radiation track structure and calculate the position of radiolytic species after irradiation. We have also recently developed an efficient Monte‐Carlo sampling algorithm for the GFDE of reversible reactions with an intermediate state [4], which can be modified and adapted to simulate DNA damage by free radicals. To do so, we will use the known reaction rate constants between radicals (OH, eaq, H,...) and the DNA bases, sugars and phosphates and use the sampling algorithms to simulate the diffusion of free radicals and chemical reactions with DNA. These techniques should help the understanding of the contribution of the indirect effect in the formation of DNA damage and double‐strand breaks.

Plante, Ianik↗

NASA NTRS1989

Stochastic histories of dust grains in the interstellar medium

The purpose is to study an evolving system of refractory dust grains within the Interstellar Medium (ISM). This is done via a combination of Monte Carlo processes and a system of partial differential equations, where refractory dust grains formed within supernova remnants and ejecta from high mass loss stars are subjected to the processes of sputtering and collisional fragmentation in the diffuse media and accretion within the cold molecular clouds. In order to record chemical detail, the authors take each new particle to consist of a superrefractory core plus a more massive refractory mantle. The particles are allowed to transfer to and fro between the different phases of the interstellar medium (ISM) - on a time scale of 10(exp 8) years - until either the particles are destroyed or the program finishes at a Galaxy time of 6x10(exp 9) years. The resulting chemical and size spectrum(s) are then applied to various astrophysical problems with the following results. For an ISM which has no collisional fragmentation of the dust grains, roughly 10 percent by mass of the most refractory material survives the rigors of the ISM intact, which leaves open the possibility that fossilized isotopically anomalous material may have been present within the primordial solar nebula. Stuctured or layered refractory dust grains within the model cannot explain the observed interstellar depletions of refractory material. Fragmentation due to grain-grain collisions in the diffuse phase plus the accretion of material in the molecular cloud phase can under certain circumstances cause a bimodal distribution in grain size.

Liffman, Kurt↗

NASA NTRS1976

Neural pulse frequency modulation of an exponentially correlated Gaussian process

The effect of NPFM (Neural Pulse Frequency Modulation) on a stationary Gaussian input, namely an exponentially correlated Gaussian input, is investigated with special emphasis on the determination of the average number of pulses in unit time, known also as the average frequency of pulse occurrence. For some classes of stationary input processes where the formulation of the appropriate multidimensional Markov diffusion model of the input-plus-NPFM system is possible, the average impulse frequency may be obtained by a generalization of the approach adopted. The results are approximate and numerical, but are in close agreement with Monte Carlo computer simulation results.

Hutchinson, C. E.↗

NASA NTRS1977

Determination of three-dimensional radiative exchange factors for the Space Shuttle by Monte Carlo

A program was developed to predict radiative exchange factors in three-dimensional configurations with curved, specularly-reflecting surfaces using the Monte Carlo method. It was shown that this can be done efficiently by using the normal Monte Carlo method for closed configurations and/or high surface emissivities, while Energy Partitioning should be used for open configurations and/or low surface emissivities. Comparison with experiments on the Space Shuttle radiator panels showed that the simple model of gray, diffuse emission is sufficient for the prediction of exchange factors between actual (silver-coated teflon) surfaces. For solar-irradiation exchange factors on the Space Shuttle, on the other hand, spectral and directional variations of surface properties must be taken into account.

Modest, M. F.↗

NASA NTRS1959

A Numerical Method for Obtaining Monoenergetic Neutron Flux Distributions and Transmissions in Multiple-Region Slabs

This method is investigated for semi-infinite multiple-slab configurations of arbitrary width, composition, and source distribution. Isotropic scattering in the laboratory system is assumed. Isotropic scattering implies that the fraction of neutrons scattered in the i(sup th) volume element or subregion that will make their next collision in the j(sup th) volume element or subregion is the same for all collisions. These so-called "transfer probabilities" between subregions are calculated and used to obtain successive-collision densities from which the flux and transmission probabilities directly follow. For a thick slab with little or no absorption, a successive-collisions technique proves impractical because an unreasonably large number of collisions must be followed in order to obtain the flux. Here the appropriate integral equation is converted into a set of linear simultaneous algebraic equations that are solved for the average total flux in each subregion. When ordinary diffusion theory applies with satisfactory precision in a portion of the multiple-slab configuration, the problem is solved by ordinary diffusion theory, but the flux is plotted only in the region of validity. The angular distribution of neutrons entering the remaining portion is determined from the known diffusion flux and the remaining region is solved by higher order theory. Several procedures for applying the numerical method are presented and discussed. To illustrate the calculational procedure, a symmetrical slab ia vacuum is worked by the numerical, Monte Carlo, and P(sub 3) spherical harmonics methods. In addition, an unsymmetrical double-slab problem is solved by the numerical and Monte Carlo methods. The numerical approach proved faster and more accurate in these examples. Adaptation of the method to anisotropic scattering in slabs is indicated, although no example is included in this paper.

Schneider, Harold↗

NASA NTRS1999

Diffuse Reflection of Laser Light From Clouds

Laser light reflected from an aqueous suspension of particles or "cloud" with known thickness and particle size distribution defines the "cloud radiative Green's function", G. G is sensitive to cloud thickness, allowing retrieval of that important quantity. We describe a laboratory simulation of G, useful in design of an offbeam Lidar instrument for remote sensing of cloud thickness. Clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The size distribution extends from 0.5 microns to 25 microns, roughly lognormal, similar to real clouds. Density of suspended spheres is adjusted so photon mean-free-path is about 10 cm, 1000 times smaller than in real clouds. The light source is a Nd:YAG laser at 530 nm. Detectors are flux and photon-counting PMTs, with a glass probe for precise positioning. A Labview 5 VI controls position and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider,and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns. Step size is selectable. Far from the beam, the rate of exponential increase in the beam direction scales as expected from diffusion theory, linearly with cloud thickness, and inversely as the square root of the reduced optical thickness, independent of particle size. Nearer the beam the signal increases faster than exponential and depends on particle size. Results verify 3D Monte Carlo simulations that demonstrate detectability of remotely sensed offbeam returns, without filters at night, with narrow bandpass filter in day.

Cahalan, R. F.↗

NASA NTRS1977

Monte Carlo analysis of uncertainty propagation in a stratospheric model. 1: Development of a concise stratospheric model

A concise model has been developed to analyze uncertainties in stratospheric perturbations, yet uses a minimum of computer time and is complete enough to represent the results of more complex models. The steady state model applies iteration to achieve coupling between interacting species. The species are determined from diffusion equations with appropriate sources and sinks. Diurnal effects due to chlorine nitrate formation are accounted for by analytic approximation. The model has been used to evaluate steady state perturbations due to injections of chlorine and NO(X).

Rundel, R. D.↗

NASA NTRS1979

Geometry dependence of optical pulse broadening in multiple scattering media

A Monte Carlo simulation method has been used to evaluate the temporal response at a variety of receivers to an instantaneous laser pulse input to a turbid medium. The results are compared to the description of the ensemble of all transmitted photons for which a diffusion theory was developed. It was shown that the arrival times for all photons transmitted through the cloud are exponentially distributed with a mean time that varies as the product of the cloud thickness and the momentum transfer optical depth. For large optical depths the mean time delay associated with all reflected photons was found to depend only on the cloud thickness. The impulse response function for receivers on the beam axis varied with field of view and for small fields of view it was dominated by single scattering even for optically thick clouds. In these cases, the mean time delay could be up to 10 to the 6th times shorter than that associated with all transmitted photons with a related increase in the maximum modulation frequency.

Thomas, R. W. L.↗

NASA NTRS1981

Verification of reflectance models in turbid waters

Inherent optical parameters of very turbid waters were used to evaluate existing water reflectance models. Measured upwelling radiance spectra and Monte Carlo simulations of the radiative transfer equations were compared with results from models based upon two flow, quasi-single scattering, augmented isotropic scattering, and power series approximation. Each model was evaluated for three separate components of upwelling radiance: (1) direct sunlight; (2) diffuse skylight; and (3) internally reflected light. Limitations of existing water reflectance models as applied to turbid waters and possible applications to the extraction of water constituent information are discussed.

Tanis, F. J.↗

NASA NTRS1990

Numerical simulation of rarefied gas flow through a slit

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas from one reservoir to another through a two-dimensional slit. The cases considered are for hard vacuum downstream pressure, finite pressure ratios, and isobaric pressure with thermal diffusion, which are not well established in spite of the simplicity of the flow field. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, three kinds of collision sampling techniques, the time counter (TC) method, the null collision (NC) method, and the no time counter (NTC) method, are used.

Keith, Theo G., Jr.↗