Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ray trace”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Measuring the Stress Factors for Photovoltaic (PV) Backsheet Degradation

Back sheet failure has resulted in power loss and large-scale recall of photovoltaic modules, resulting in billions of dollars in lost revenue. The light exposure on the backside of a photovoltaic module comes primarily from reflected light which alters the distribution of natural sunlight. Because of this, modelling the backside exposure and duplicating the exposure is much more difficult than modeling the frontside exposure. This project aims to study how various back sheets and junction box materials degrade under different conditions and to develop Python code to help model and predict degradation. The stress factors for back-sheet degradation must be quantified to extrapolate accelerated stress tests to the field. Test samples were placed in the A3, A4, and A5 conditions, as defined in IEC 62788-7-2, to assess the temperature and humidity dependence of ultraviolet (UV) induced degradation. We are utilizing a custom chamber with exposure from 0.5 UV-suns to 5 UV-suns to understand the dependence of degradation on light intensity. A group of samples put in the A3 condition had glass filters with 50% UV cut-offs of 320 nm, 335 nm, and 360 nm to assess the wavelength dependence of UV degradation. All this data is necessary to assess the impact of non-standard UV light exposure. The material evaluation tests include gloss measurements, attenuated total internal reflectance Fourier transform infrared spectroscopy (ATR-FTIR), UV-visible reflectance/transmittance utilizing a Cary Ci7000 spectrophotometer, and a nano-indenter for surface hardness and modulus measurements. Alongside the experimental work, there is a computational effort using raytracing and Python open-source tools in PVDeg , PVLib, and Bifacial_Radiance. This code will create specific exposure scenarios and enable the evaluation of chamber degradation relative to field degradation. Equation 1 is a strawman equation used to model degradation on the backside of a PV module. We will create simplified code, based on the results of ray-tracing calculations, which uses a view factor approach to provide fast calculations for the most common exposure scenarios.

14 SOLAR ENERGY↗
DOE OSTI2026

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Using Numba for GPU acceleration of Neutron Beamline Digital Twins

Digital twins of neutron instruments using Monte Carlo ray tracing have proven to be useful in neutron data analysis and verifying instrument and sample designs. However, these simulations can become quite complex and computationally demanding with tens of billions of neutrons. In this paper, we present a GPU accelerated version of MCViNE using Python and Numba to balance user extensibility with performance. Numba is an open-source just-in-time (JIT) compiler for Python using LLVM to generate efficient machine code for CPUs and GPUs with NVIDIA CUDA. The JIT nature of Numba allowed complex instrument kernels to be generated easily. Initial simulations have shown a speedup between 200-1000x over the original CPU implementation. The performance gain with Numba enables more sophisticated data analysis and impacts neutron scattering science and instrument design.

Kendrick, Coleman↗
DOE OSTI2022

Spectral Radiative Properties of Polydispersed SiO 2 Particle Beds

The focus of this work is on the measurement and analysis of the radiative properties of polycrystalline SiO 2 particle beds with various layer thicknesses. The particles are polydispersed with average diameters of 222, 150, and 40 μm . The spectral, directional–hemispherical reflectance and transmittance of the particle bed are measured at wavelengths from 0.4 to 1.8 μm using a monochromator, and the reflectance measurement is extended to 15 μm using a Fourier-transform infrared spectrometer. Particles are closely packed between two transparent windows for measuring the radiative properties. In the visible and near-infrared region up to 1.8 μm, the inverse adding–doubling method yields the effective absorption and scattering coefficients. The results suggest that short wavelength absorption needs to be included in modeling the behavior of particle beds due to multiple scattering. A discrete-scale Monte Carlo ray-tracing method is developed to model the radiative properties by assuming monodispersed spherical particles, and the simulated results compare well with measurements. The effective absorption and scattering coefficients of the particle beds obtained from the independent scattering theory are compared to those from the inverse method. As a result, the impact of dependent scattering on the packed beds is observed for smaller-sized particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Computational Optimization of A 3D Printed Collimator

This contribution describes the computational methodology behind an optimization procedure for a scattered beam collimator. The workflow includes producing a file that can be manufactured via additive methods. A conical collimator, optimized for neutron diffraction experiments in a high pressure clamp cell, is presented as an example. In such a case the scattering from the sample is much smaller than that of the pressure cell. Monte Carlo Ray tracing in MCViNE was used to model scattering from a Si powder sample and the cell. A collimator was inserted into the simulation and the number and size of channels were optimized to maximize the rejection of the parasitic signal coming from the complex sample environment. Constraints, provided by the additive manufacturing process as well as a specific neutron diffractometer, were also included in the optimization. The source code and the tutorials are available in c3dp (Islam (2019)).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

How to design a focusing guide: The large moderator case1

As continuously shaped super-mirrors are becoming available, the conceptual design of focusing guides should explore a wider range of possibilities to accomplish an efficient neutron beam extraction. Starting from a desired phase-space volume at the sample position and using an upstream ray-tracing approach, the acceptance diagram of any focusing guide can be calculated at the moderator position. To ensure high brilliance transfer and homogeneous coverage, the acceptance diagram should be fully included in the neutron source emission phase-space volume. Following this idea, the guide system can be scaled into dimensionless geometric figures that convey performance limits for a desired cross-section reduction. Moreover, if we impose a monotonic increase of the reflection angle with divergence angle at the sample position, the shape of the mirror is analytically determined. This approach was applied in the design of a focusing guide for SNAP instrument at SNS, at ORNL, USA. The results of McStas simulations are presented with different options included. Our approach facilitates finding an optimal solution for connecting multiple guide pieces to avoid excessive losses and ensure a homogeneous phase space coverage.

Frost, Matthew↗
DOE OSTI2024

Systematic Effects in Galaxy-Galaxy Lensing with DESI

The Dark Energy Spectroscopic Instrument (DESI) survey will measure spectroscopic redshifts for millions of galaxies across roughly $14,000 \, \mathrm{deg}^2$ of the sky. Cross-correlating targets in the DESI survey with complementary imaging surveys allows us to measure and analyze shear distortions caused by gravitational lensing in unprecedented detail. In this work, we analyze a series of mock catalogs with ray-traced gravitational lensing and increasing sophistication to estimate systematic effects on galaxy-galaxy lensing estimators such as the tangential shear $\gamma_{\mathrm{t}}$ and the excess surface density $\Delta\Sigma$ . We employ mock catalogs tailored to the specific imaging surveys overlapping with the DESI survey: the Dark Energy Survey (DES), the Hyper Suprime-Cam (HSC) survey, and the Kilo-Degree Survey (KiDS). Among others, we find that fiber incompleteness can have significant effects on galaxy-galaxy lensing estimators but can be corrected effectively by up-weighting DESI targets with fibers by the inverse of the fiber assignment probability. Similarly, we show that intrinsic alignment and lens magnification are expected to be statistically significant given the precision forecasted for the DESI year-1 data set. Our study informs several analysis choices for upcoming cross-correlation studies of DESI with DES, HSC, and KiDS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Lighting application efficacy: A framework for holistically measuring lighting use in buildings

Lighting consumes significant energy in buildings, but current measurements of lighting efficiency are inadequate in capturing the spatial and temporal effectiveness of light. The most widely used lighting energy efficiency measure, luminous efficacy, quantifies the amount of light generated by individual light sources. While luminous efficacy provides manufacturers and designers with a useful metric for quantifying the energy efficiency of individual lighting devices, it does not fully address the use of lighting in architectural spaces. Instead, application efficacy, the relationship between the electrical power consumed by lighting hardware and the amount of light that contributes to the visual perception of building occupants, must be measured. Here, we propose a framework for calculating lighting application efficacy (LAE) based on the primary pathway of light in architectural settings: the generation and emission of light from a luminaire, the travel of the light throughout the space and into occupants’ eyes, and the process of visual perception. The LAE framework will include the development of computational methods for estimating the proportion of emitted light that is directed to task areas and areas corresponding to occupants’ visual fields using ray-tracing lighting design software. Future modeling and validation work will quantify energy efficiency of buildings.

Durmus, Dorukalp↗
DOE OSTI2019

Application and Evaluation of Surrogate Models for Radiation Source Search

Surrogate models are increasingly required for applications in which first-principles simulation models are prohibitively expensive to employ for uncertainty analysis, design, or control. They can also be used to approximate models whose discontinuous derivatives preclude the use of gradient-based optimization or data assimilation algorithms. We consider the problem of inferring the 2D location and intensity of a radiation source in an urban environment using a ray-tracing model based on Boltzmann transport theory. Whereas the code implementing this model is relatively efficient, extension to 3D Monte Carlo transport simulations precludes subsequent Bayesian inference to infer source locations, which typically requires thousands to millions of simulations. Additionally, the resulting likelihood exhibits discontinuous derivatives due to the presence of buildings. To address these issues, we discuss the construction of surrogate models for optimization, Bayesian inference, and uncertainty propagation. Specifically, we consider surrogate models based on Legendre polynomials, multivariate adaptive regression splines, radial basis functions, Gaussian processes, and neural networks. We detail strategies for computing training points and discuss the merits and deficits of each method.

Cook, Jared A.↗
DOE OSTI2022

Influence of Superlattice Structure on V-Defect Distribution, External Quantum Efficiency and Electroluminescence for Red InGaN Based µLEDs on Silicon

Achieving high quantum efficiency in long-wavelength LEDs has posed a significant challenge to the solid-state lighting and display industries. In this article, we use V-defect engineering as a technique to achieve higher efficiencies in red InGaN LEDs on (111) Si through lateral injection. We investigate the effects of superlattice structure on the V-defect distribution, the electroluminescence properties, and the external quantum efficiency. Increasing the relative thickness of In in the InGaN/GaN superlattice and the total superlattice thickness correlate with a reduction of active region defects and increased external quantum efficiencies. The highest measured on-chip EQE was 0.15% and based on Monte-Carlo ray tracing simulations for light extraction we project this would correspond to a flip-chip EQE of ~2.5%.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Volumetric Rendering on Wavelet-Based Adaptive Grid

Numerical modeling of physical phenomena frequently involves processes across a wide range of spatial and temporal scales. In the last two decades, the advancements in wavelet-based numerical methodologies to solve partial differential equations, combined with the unique properties of wavelet analysis to resolve localized structures of the solution on dynamically adaptive computational meshes, make it feasible to perform large-scale numerical simulations of a variety of physical systems on a dynamically adaptive computational mesh that changes both in space and time. Volumetric visualization of the solution is an essential part of scientific computing, yet the existing volumetric visualization techniques do not take full advantage of multi-resolution wavelet analysis and are not fully tailored for visualization of a compressed solution on the wavelet-based adaptive computational mesh. Our objective is to explore the alternatives for the visualization of time-dependent data on space-time varying adaptive mesh using volume rendering while capitalizing on the available sparse data representation. Two alternative formulations are explored. The first one is based on volumetric ray casting of multi-scale datasets in wavelet space. Rather than working with the wavelets at the finest possible resolution, a partial inverse wavelet transform is performed as a preprocessing step to obtain scaling functions on a uniform grid at a user-prescribed resolution. As a result, a solution in physical space is represented by a superposition of scaling functions on a coarse regular grid and wavelets on an adaptive mesh. An efficient and accurate ray casting algorithm is based just on these coarse scaling functions. Additional details are added during the ray tracing by taking an appropriate number of wavelets into account based on support overlap with the interpolation point, wavelet coefficient magnitude, and other characteristics, such as opacity accumulation (front to back ordering) and deviation from frontal viewing direction. The second approach is based on complementing of wavelet-based adaptive mesh to the traditional Adaptive Mesh Refinement (AMR) mesh. Both algorithms are illustrated and compared to the existing volume visualization software for Rayleigh-Benard thermal convection and electron density data sets in terms of rendering time and visual quality for different data compression of both wavelet-based and AMR adaptive meshes.

Vezolainen, Alexei V.↗
DOE OSTI2016

Disk emission from magnetohydrodynamic simulations of spinning black holes

We present the results of a new series of global, three-dimensional, relativistic magnetohydrodynamic (MHD) simulations of thin accretion disks around spinning black holes. The disks have aspect ratios of H/R∼0.05 and spin parameters of a/M=0,0.5,0.9, and 0.99. Using the ray-tracing code Pandurata, we generate broadband thermal spectra and polarization signatures from the MHD simulations. We find that the simulated spectra can be well fit with a simple, universal emissivity profile that better reproduces the behavior of the emission from the inner disk, compared to traditional analyses carried out using a Novikov–Thorne thin disk model. Finally, we show how spectropolarization observations can be used to convincingly break the spin-inclination degeneracy well known to the continuum-fitting method of measuring black hole spin.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ODYSSEY: A PUBLIC GPU-BASED CODE FOR GENERAL RELATIVISTIC RADIATIVE TRANSFER IN KERR SPACETIME

General relativistic radiative transfer calculations coupled with the calculation of geodesics in the Kerr spacetime are an essential tool for determining the images, spectra, and light curves from matter in the vicinity of black holes. Such studies are especially important for ongoing and upcoming millimeter/submillimeter very long baseline interferometry observations of the supermassive black holes at the centers of Sgr A* and M87. To this end we introduce Odyssey, a graphics processing unit (GPU) based code for ray tracing and radiative transfer in the Kerr spacetime. On a single GPU, the performance of Odyssey can exceed 1 ns per photon, per Runge–Kutta integration step. Odyssey is publicly available, fast, accurate, and flexible enough to be modified to suit the specific needs of new users. Along with a Graphical User Interface powered by a video-accelerated display architecture, we also present an educational software tool, Odyssey-Edu, for showing in real time how null geodesics around a Kerr black hole vary as a function of black hole spin and angle of incidence onto the black hole.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

REFLECTED LIGHT CURVES, SPHERICAL AND BOND ALBEDOS OF JUPITER- AND SATURN-LIKE EXOPLANETS

Reflected light curves observed for exoplanets indicate that a few of them host bright clouds. We estimate how the light curve and total stellar heating of a planet depends on forward and backward scattering in the clouds based on Pioneer and Cassini spacecraft images of Jupiter and Saturn. We fit analytical functions to the local reflected brightnesses of Jupiter and Saturn depending on the planet’s phase. These observations cover broadbands at 0.59–0.72 and 0.39–0.5 μ m, and narrowbands at 0.938 (atmospheric window), 0.889 (CH4 absorption band), and 0.24–0.28 μ m. We simulate the images of the planets with a ray-tracing model, and disk-integrate them to produce the full-orbit light curves. For Jupiter, we also fit the modeled light curves to the observed full-disk brightness. We derive spherical albedos for Jupiter and Saturn, and for planets with Lambertian and Rayleigh-scattering atmospheres. Jupiter-like atmospheres can produce light curves that are a factor of two fainter at half-phase than the Lambertian planet, given the same geometric albedo at transit. The spherical albedo is typically lower than for a Lambertian planet by up to a factor of ∼1.5. The Lambertian assumption will underestimate the absorption of the stellar light and the equilibrium temperature of the planetary atmosphere. We also compare our light curves with the light curves of solid bodies: the moons Enceladus and Callisto. Their strong backscattering peak within a few degrees of opposition (secondary eclipse) can lead to an even stronger underestimate of the stellar heating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

The kinetic Sunyaev-Zel'dovich effect from reionization: simulated full-sky maps at arcminute resolution

The kinetic Sunyaev–Zel’dovich (kSZ) effect results from Thomson scattering by coherent flows in the reionized intergalactic medium. We present new results based on ray-tracing an 8 Gpc/h realization of reionization with resolution elements 2 Mpc/h (subtending ∼1′ at z = 6) on a side to create a full-sky kSZ map. The realization includes, self-consistently, the effects of reionization on scales corresponding to multipoles 10≲ℓ≲5000. We separate the kSZ map into Doppler (v), Ostriker–Vishniac (δv), patchy (xv), and third-order (xδv) components, and compute explicitly all the auto- and cross-correlations (e.g., 〈vv〉, 〈δvxv〉, etc.) that contribute to the total power. We find a complex and nonmonotonic dependence on the duration of reionization at ℓ∼300 and evidence for a non-negligible (10%–30%) contribution from connected four-point correlations, 〈xvxv〉{sub c}, usually neglected in analytical models. We also investigate the cross-correlation of linear matter and large-scale kSZ temperature fluctuations, focusing on (1) cross-power spectra with biased tracers of the matter density and (2) cold spots from infall onto large, rare H ii regions centered on peaks in the matter distribution at redshifts z>10 that are a generic non-Gaussian feature of patchy reionization. Finally, we show that the reionization history can be reconstructed at 5σ–10σ significance by correlating full-sky 21 cm maps stacked in bins with Δν = 10 MHz with existing cosmic microwave background (CMB) temperature maps at ℓ<500, raising the prospects for probing reionization by correlating CMB and LSS measurements. The resulting kSZ maps have been made publicly available at www.cita.utoronto.ca/∼malvarez/research/ksz-data/.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

2D RADIATION-HYDRODYNAMIC SIMULATIONS OF SUPERNOVA SHOCK BREAKOUT IN BIPOLAR EXPLOSIONS OF A BLUE SUPERGIANT PROGENITOR

A two-dimensional special relativistic radiation-hydrodynamics code is developed and applied to numerical simulations of supernova shock breakout in bipolar explosions of a blue supergiant. Our calculations successfully simulate the dynamical evolution of a blast wave in the star and its emergence from the surface. Results of the model with spherical energy deposition show a good agreement with previous simulations. Furthermore, we calculate several models with bipolar energy deposition and compare their results with the spherically symmetric model. The bolometric light curves of the shock breakout emission are calculated by a ray-tracing method. Our radiation-hydrodynamic models indicate that the early part of the shock breakout emission can be used to probe the geometry of the blast wave produced as a result of the gravitational collapse of the iron core.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

MASS MEASUREMENTS IN PROTOPLANETARY DISKS FROM HYDROGEN DEUTERIDE

The total gas mass of a protoplanetary disk is a fundamental, but poorly determined, quantity. A new technique has been demonstrated to assess directly the bulk molecular gas reservoir of molecular hydrogen using the HD J = 1–0 line at 112 μ m. In this work we present a Herschel Space Observatory {sup 10} survey of six additional T Tauri disks in the HD line. Line emission is detected at >3 σ significance in two cases: DM Tau and GM Aur. For the other four disks, we establish upper limits to the line flux. Using detailed disk structure and ray-tracing models, we calculate the temperature structure and dust mass from modeling the observed spectral energy distributions, and we include the effect of UV gas heating to determine the amount of gas required to fit the HD line. The ranges of gas masses are 1.0–4.7 × 10{sup -2} for DM Tau and 2.5–20.4 × 10{sup -2} for GM Aur. These values are larger than those found using CO for GM Aur, while the CO-derived gas mass for DM Tau is consistent with the lower end of our mass range. This suggests a CO chemical depletion from the gas phase of up to a factor of five for DM Tau and up to two orders of magnitude for GM Aur. We discuss how future analysis can narrow the mass ranges further.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Optimizing Simulation Parameters for Weak Lensing Analyses Involving Non-Gaussian Observables

We performed a series of numerical experiments to quantify the sensitivity of the predictions for weak lensing statistics obtained in ray-tracing dark matter (DM)-only simulations, to two hyper-parameters that influence the accuracy as well as the computational cost of the predictions: the thickness of the lens planes used to build past light cones and the mass resolution of the underlying DM simulation. The statistics considered are the power spectrum (PS) and a series of non-Gaussian observables, including the one-point probability density function, lensing peaks, and Minkowski functionals. Counterintuitively, we find that using thin lens planes (< 60 h {sup −1} Mpc on a 240 h {sup −1} Mpc simulation box) suppresses the PS over a broad range of scales beyond what would be acceptable for a survey comparable to the Large Synoptic Survey Telescope (LSST). A mass resolution of 7.2 × 10{sup 11} h {sup −1} M {sub ⊙} per DM particle (or 256{sup 3} particles in a (240 h {sup −1} Mpc){sup 3} box) is sufficient to extract information using the PS and non-Gaussian statistics from weak lensing data at angular scales down to 1′ with LSST-like levels of shape noise.

79 ASTRONOMY AND ASTROPHYSICS↗