Search NASA⌕ Search

SEARCH · Search NASA

Results for “The Random Ray Method”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Negative fluxes and cell-miss errors in the random ray method

The random ray method is a recently developed stochastic method for solving neutral particle transport problems based on the method of characteristics. Perhaps surprisingly for a characteristics-based method using flat sources, we note that the random ray method can produce negative fluxes which may be numerically troublesome in several situations. These occur most severely in fixed source problems where the source is in a region with a small cross section. Additionally, we briefly discuss another source of bias which can occur in similar situations, namely a ray missing a mesh with a strong source and small cross section, resulting in the entirety of the source being unphysically deposited locally. This paper describes the mechanism by which negative fluxes may occur and several different methods to mitigate their effects. These fixes are tested on an eigenvalue problem, a ‘fusion-like’ shielding problem, and a shielding problem featuring an adjoint calculation. Even when extremely coarse random ray quadratures are used such that 20%–30% of cells are missed during a given iteration, use of the preferred fix technique ensures local flux tally errors remain trivial (below 1%). The preferred fix is now the default option in SCONE and OpenMC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

The Random Ray Method Versus Multigroup Monte Carlo: The Method of Characteristics in OpenMC and SCONE

The Random Ray Method (TRRM) is a recently developed approach to solving neutral particle transport problems based on the Method of Characteristics. While the method previously has been implemented only in closed-source or limited-functionality codes, this work describes its implementation in two open-source Monte Carlo codes: OpenMC and SCONE. The random ray implementations required small modifications to the existing Multigroup Monte Carlo (MGMC) solvers, offering a rare venue for redundant, fine-grained, "apples-to-apples" speed and accuracy comparisons between transport methods. To this end, TRRM and MGMC solvers are evaluated against each other using each code's native capabilities on reactor eigenvalue problems with different degrees of energy discretization. On the C5G7 benchmark (featuring only seven energy groups), TRRM achieves a maximum pin power error comparable to or lower than that of MGMC for a given run time. On a problem with 69 energy groups, MGMC is found to scale more efficiently, obtaining a lower pin power error for a given run time. However, the defining difference between the two transport methods is found to be their vastly different uncertainty distributions. Specifically, TRRM is found to maintain similar levels of accuracy and uncertainty throughout the simulation domain whereas MGMC can exhibit orders-of-magnitude greater errors in areas of the problem that feature low neutron flux. For instance, TRRM provided an up to 373 times speed advantage compared with MGMC for computing the flux in low-flux regions in the moderator surrounding the C5G7 core.

42 ENGINEERING↗
NASA NTRS1981

Ray methods in random media

We review ray-optical methods of analyzing short-wavelength propagation in random media. The advantages and limitations of ray methods are discussed, and results of the statistical theory of ray segment fluctuations pertinent to ray tracing are summarized. The standard method of Monte Carlo ray tracing is compared to a new method which takes into account recent results on the statistics of ray segment fluctuations.

Hornstein, J.↗
DOE OSTI2021

Immortal rays: Rethinking random ray neutron transport on GPU architectures

The Random Ray Method (TRRM) is a recently developed adaptation of the Method of Characteristics for neutral particle transport simulations. TRRM has demonstrated excellent performance on 3D nuclear reactor benchmark problems using CPU-based compute systems. When porting to GPU-based systems, however, new performance challenges arise that are unique to processors targeting massive fine-grained parallelism. For smaller problems, or for large problems that are domain decomposed across many computational nodes, the problem size per node has insufficient parallelism to saturate GPU node resources, thus greatly limiting speedup. In this study, we report on a newly developed “immortal ray” variant of TRRM. Here, the immortal ray technique exposes significantly more fine-grained parallelism by fundamentally reformulating the numerical details of ray discretization, resulting in performance tradeoffs with significant overall benefit on GPUs. For very small 2D simulation problems we found the new immortal ray variant allowed for up to a 4.4x speedup when run on a single GPU. For larger 3D simulation problems we found the new variant improved strong scaling by 3x when run on the Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Numerical and experimental analyses of the radiant heat flux produced by quartz heating systems

A method is developed for predicting the radiant heat flux distribution produced by tungsten filament, tubular fused-quartz envelope heating systems with reflectors. The method is an application of Monte Carlo simulation, which takes the form of a random walk or ray tracing scheme. The method is applied to four systems of increasing complexity, including a single lamp without a reflector, a single lamp with a Hat reflector, a single lamp with a parabolic reflector, and up to six lamps in a six-lamp contoured-reflector heating unit. The application of the Monte Carlo method to the simulation of the thermal radiation generated by these systems is discussed. The procedures for numerical implementation are also presented. Experiments were conducted to study these quartz heating systems and to acquire measurements of the corresponding empirical heat flux distributions for correlation with analysis. The experiments were conducted such that several complicating factors could be isolated and studied sequentially. Comparisons of the experimental results with analysis are presented and discussed. Good agreement between the experimental and simulated results was obtained in all cases. This study shows that this method can be used to analyze very complicated quartz heating systems and can account for factors such as spectral properties, specular reflection from curved surfaces, source enhancement due to reflectors and/or adjacent sources, and interaction with a participating medium in a straightforward manner.

Turner, Travis L.↗
NASA NTRS1973

A numerical study of rays in random media

Statistics of electromagnetic rays in a random medium are studied numerically by the Monte Carlo method. Two dimensional random surfaces with prescribed correlation functions are used to simulate the random media. Rays are then traced in these sample media. Statistics of the ray properties such as the ray positions and directions are computed. Histograms showing the distributions of the ray positions and directions at different points along the ray path as well as at given points in space are given. The numerical experiment is repeated for different cases corresponding to weakly and strongly random media with isotropic and anisotropic irregularities. Results are compared with those derived from theoretical investigations whenever possible.

Youakim, M. Y.↗
DOE OSTI2022

Parameter estimation for X-ray scattering analysis with Hamiltonian Markov Chain Monte Carlo

Bayesian-inference-based approaches, in particular the random-walk Markov Chain Monte Carlo (MCMC) method, have received much attention recently for X-ray scattering analysis. Hamiltonian MCMC, a state-of-the-art development in the field of MCMC, has become popular in recent years. It utilizes Hamiltonian dynamics for indirect but much more efficient drawings of the model parameters. We described the principle of the Hamiltonian MCMC for inversion problems in X-ray scattering analysis by estimating high-dimensional models for several motivating scenarios in small-angle X-ray scattering, reflectivity, and X-ray fluorescence holography. Hamiltonian MCMC with appropriate preconditioning can deliver superior performance over the random-walk MCMC, and thus can be used as an efficient tool for the statistical analysis of the parameter distributions, as well as model predictions and confidence analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1981

Computing diffuse reflection from particulate planetary surface with a new function

An equation is derived to compute the amount of diffuse light reflected by a particulate surface, such as on Mars or an asteroid. The method traces the paths of rays within an ensemble of randomly shaped grains and finds the eventual probability of emission. The amount of diffuse, unpolarized emitted light is obtained in terms of the real index of refraction, the imaginary index, and the average diameter of particles making up the surface. The equation is used to compute the empirical rule for obtaining the planetary albedo from the slope of its polarization curve. Accuracy of the equation, estimated at + or - 4%, seems justified because of quantitative agreement with experimental measures of the empirical rule. It is also shown that the equation can be applied to bubble-enclosing surfaces such as volcanic foams. Results for the indices of the moon, Mars, Io, and Europa are obtained and compared with other data.

Wolff, M.↗
DOE OSTI2022

Determining Infrared Optical Constants for Dolomite using Single Angle Reflectance and Spectroscopic Ellipsometry

The optical constants, namely the real (n) and imaginary (k) parts of the complex refractive index, are of particular interest to generate the infrared (IR) spectra of liquid and solid materials in different morphologies. To obtain n/k, however, most materials are typically not found in the monolithic forms necessary to easily measure n/k, and thus require the use of other methods such as pressing powders to form planar, specularly-reflective pellets. In this work dolomite crystals are measured using fixed-angle IR reflectance spectroscopy in both 1) a monolithic form with the crystal fixed in epoxy and polished, but data are also recorded from 2) pressed-pellet forms of powders of the same mineral using the same spectroscopic protocols. First results for the two methods are compared. It was found that the measured reflectance can vary by as much as a factor of two between the dolomite crystals and the pressed powder forms. For the two-sample preparation approaches the preliminary spectra are compared and the implications and limitations of each method for determining the optical constants of given materials are discussed. Comparison to literature data suggest that polarization effects likely account for the differing amplitude results for reflectance (and hence the k -vectors): Dolomite is known to be biaxial with significantly differing ordinary and extraordinary rays; the pressed pellet method measures an ensemble of microcrystals in randomly oriented positions, whereas the single crystal maintains just one orientation relative to the optical axis.

ellipsometry, dolomite, complex refractive index, ↗
NASA NTRS1973

Investigation of resonance integrals occurring in cosmic ray diffusion theory

Critical assessment of two versions of a procedure for calculating the pitch angle diffusion coefficient for cosmic rays in a static random magnetic field using the 'resonance integral' method of Hasselmann and Wiberenz (1968) and Jokipii (1972). One of these versions is shown to represent the physics of the situation more accurately than the other.

Jones, F. C.↗
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 NTRS2008

Effect of the Inhomogeneity of Ice Crystals on Retrieving Ice Cloud Optical Thickness and Effective Particle Size

Spherical or spheroidal air bubbles are generally trapped in the formation of rapidly growing ice crystals. In this study the single-scattering properties of inhomogeneous ice crystals containing air bubbles are investigated. Specifically, a computational model based on an improved geometric-optics method (IGOM) has been developed to simulate the scattering of light by randomly oriented hexagonal ice crystals containing spherical or spheroidal air bubbles. A combination of the ray-tracing technique and the Monte Carlo method is used. The effect of the air bubbles within ice crystals is to smooth the phase functions, diminish the 22deg and 46deg halo peaks, and substantially reduce the backscatter relative to bubble-free particles. These features vary with the number, sizes, locations and shapes of the air bubbles within ice crystals. Moreover, the asymmetry factors of inhomogeneous ice crystals decrease as the volume of air bubbles increases. Cloud reflectance lookup tables were generated at wavelengths 0.65 m and 2.13 m with different air-bubble conditions to examine the impact of the bubbles on retrieving ice cloud optical thickness and effective particle size. The reflectances simulated for inhomogeneous ice crystals are slightly larger than those computed for homogenous ice crystals at a wavelength of 0.65 microns. Thus, the retrieved cloud optical thicknesses are reduced by employing inhomogeneous ice cloud models. At a wavelength of 2.13 microns, including air bubbles in ice cloud models may also increase the reflectance. This effect implies that the retrieved effective particle sizes for inhomogeneous ice crystals are larger than those retrieved for homogeneous ice crystals, particularly, in the case of large air bubbles.

Xie, Yu↗
DOE OSTI2020

Complex Disorder in Type-I Clathrates: Synthesis and Structural Characterization of A8GaxSn46−x (A = Rb, Cs; 6.9 < x < 7.5)

Exploratory studies in the systems Rb–Ga–Sn and Cs–Ga–Sn yielded the cubic type-I clathrates with refined compositions Rb8GaxSn46−x and Cs8GaxSn46−x (6.9 < x < 7.5). Nearly single-phase materials with good crystallinity were obtained from stoichiometric reactions of the elements. The structures were characterized by means of single-crystal X-ray diffraction methods. Both Rb8GaxSn46−x and Cs8GaxSn46−x represents cases, where a Group 13 element randomly substitutes a Group 14 element in the structure. The extent of Ga/Sn mixing is apparently governed by the drive of the system to achieve an optimal valence electron count, and hence, Rb8GaxSn46−x and Cs8GaxSn46−x (x ≈ 8) can be regarded as Zintl phases. This notion is supported by structure refinements on a multitude of single-crystal X-ray diffraction data, which also confirm that both types of cages in the cubic type-I structure are fully occupied by Rb and Cs atoms. The open-framework, comprised of 46 nodes per formula unit, adapts to the incorporation of nearly eight Ga atoms within the matrix of Sn, whereby small, short-range distortions result. The exact nature of these effects is still unclear, as so far, the structural variations could only be modeled as both positional and occupational disorder at one of three framework sites. Since vacancies in the structures of the binary type-I clathrates A8Sn46−x☐x (A = Rb, Cs; ☐ = missing Sn atom) are also known to cause local distortions, the latter were also synthesized with the same protocols used for the synthesis of A8GaxSn46−x and structurally re-analyzed. The results from the latter studies confirm that homogeneity issues abound, and that the final structures/compositions are an intricate function of the experimental conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Signatures of multiple charge excitations in resonant inelastic x-ray scattering spectra of metals

In this work, we study how multiple charge excitations appear in the resonant inelastic x-ray scattering (RIXS) spectra of metals. Typically, the focus is upon single excitations in the problem, the plasmons and electron-hole pairs, while the multiexcitation processes are neglected. However, at small momentum transfer the multiexcitation contributions may dominate the signal and one needs to understand how to interpret the data in such cases. In particular, we demonstrate how to “decode” the total multiexcitation intensity and extract the plasmon dispersion. Furthermore, while our calculations are based on the random phase approximation for jellium, which does not allow us to obtain quantitatively precise results in the entire region of parameters, we expect them to capture semiqualitatively all features expected for charged Fermi-liquid states, including universal and singular properties of the RIXS spectra.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1984

SMC X-1 variability observed from HEAO 1

Studies are reported of the slow variability of SMC X-1 and its spectrum. An analysis of red-noise random variability is based on a method discussed by Deeter and Boynton (1982). The 0.7 s X-ray pulsar SMC X-1 is in a 3.89 day eclipsing binary system with a B0 I supergiant companion. Observations of the pulsar were conducted with the aid of the UCSD/MIT instrument on HEAO 1 from 1977 August through 1979 January. A light curve was constructed for the period 1977 September to 1978 December. The apparent tendency of SMC X-1 to be in one of two states, high or low, suggests the acquisition of average spectra separately for each state. The total (steady plus pulsed) emission from SMC X-1 displays a continuum spectrum with a dominant exponential form which implies a temperature of 17 keV for thin thermal bremsstrahlung emission or 5 keV if the other limit of a Wien spectrum is assumed.

Gruber, D. E.↗
NASA NTRS1992

The effects of small ice crystals on cirrus infrared radiative properties

Using a model that combines single-scattering properties for spheroidal and hexagonal ice crystals, the thermal infrared radiative properties of cirrus clouds have been investigated. Infrared scattering and absorption properties for randomly oriented spheroids and hexagons are parameterized based on the anomalous diffraction theory and a geometric ray-tracing method, respectively. Using observed ice crystal size distributions, upwelling radiances at the top of cirrus cloudy atmospheres have been computed. Results show that the presence of small ice crystals can produce significant brightness temperature differences between two infrared wavelengths in the 10-micron window. Theoretical results have been compared with observed brightness temperature differences between 8.35 and 11.16 microns and between 11.16 and 12 microns. The observed values were obtained from the high-spectral resolution interferometer sounder. It is shown that the use of the present nonspherical model for ice crystals in radiative transfer calculations leads to a significantly better interpretation of the observed data than does the use of the spherical model.

Takano, Y.↗