Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ray tracing”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modeling thin layers of analytes on substrates for spectral analysis: use of solid/liquid n and k values to model reflectance spectra

Since solids are only sometimes seen en masse in a pure bulk form, and for liquids other than water almost never, a capability to model reflectance spectra from analytes deposited on various substrates would be highly advantageous. If available, the real, n ( ν ) , and imaginary, k ( ν ) , components of the complex refractive index, n n + i k , can be used to simulate infrared spectra, accounting for reflection, refraction, and absorption phenomena as a function of wavelength. We focus on using the Pacific Northwest National Laboratory (PNNL) derived n / k vectors for solid and liquid analytes deposited as thin layers on different types of substrates including conductors, such as aluminum, and inorganic dielectrics, such as glass. The model is an adaptation of the Monte Carlo ray trace modeling program, TracePro, extended through the use of its macrolanguage. The model is tested using thin films of organic liquids including silicone oil and no. 2 diesel fuel, as well as organic solids such as caffeine and acetaminophen on aluminum and glass. The predicted spectra for the solid films were compared to experimental hemispherical reflectance data measured using a Fourier transform spectrometer with an integrating sphere. The thickness of the calculated layer is a parameter for predicting the (transflectance) spectra and is obtained using the areal density measured from gravimetric methods to generate the thin-layer samples. Comparison of the calculated spectra with experimental hemispherical reflectance data shows excellent agreement, indicating promise for the use of measured n / k data to synthesize reference spectral data. In particular, accounting for the inhomogeneity of the deposits greatly improved the match with experimental data. Finally, the theoretical modeling shows that for thicker layers (ca. 20 to 100 μm) of typical organics possessing moderately strong k values, the longwave infrared features are often saturated and better spectral contrast is obtained from the overtone/combination bands in the shortwave infrared.

42 ENGINEERING↗
DOE OSTI2020

Analytical expressions of the surface shape of ‘diaboloid’ mirrors

Modern deterministic polishing processes allow fabrication of x-ray optics with almost any arbitrary aspherical surface shape. Among these optics, the so called "diaboloid"mirror is of special interest. The diaboloid mirror that converts a cylindrical wave to a spherical wave would improve focusing in x-ray beamlines implementing a diffraction element between a parabolic cylinder and a toroidal mirror. The replacement of the toroidal mirror in existing beamlines by the diaboloid mirror would mitigate aberrations. The shape of the diaboloid mirror is usually calculated numerically based on a truncated polynomial solution of the optical path problem. Here, we present an exact analytical solution for the shape of a diaboloid mirror as a function of the conjugate parameters of the mirror placed in a beamline. The derived analytical expressions for the diaboloid mirror in both the canonical and mirror-based coordinate systems are implemented in ray-tracing simulations to verify the beamline performances.

Yashchuk, Valeriy V.↗
DOE OSTI2022

Near-field acoustical holography and acoustic power analysis of a simulated, highly heated supersonic jet

Although near-field acoustical holography (NAH) and acoustic intensity analysis have previously been used to investigate the apparent jet noise sources produced by military aircraft, explicit connections to supersonic jet characteristics cannot be made due to a lack of information about the exhaust plume. Here, to begin to bridge this gap and better understand the source information yielded by NAH, the current study instead applies NAH to a virtual measurement of the near-field pressures of a highly heated laboratory-scale supersonic jet generated by large-eddy simulation (LES). The holographic reconstructions of the pressure, particle velocity, and acoustic intensity are found to match the LES-generated acoustic field well and are used to calculate the acoustic power of the jet. The jet's calculated overall acoustic power is compared to the free-stream mechanical power, resulting in an acoustic efficiency of 1.5%. Ray-tracing of the acoustic intensity to the jet centerline generates an axial distribution of the acoustic power origin, showing that almost all the power originates from the supersonic portion of the flow and with the distribution peak upstream of the potential core tip. Holographic reconstruction of the pressures along the nozzle lipline captures the general spectral shape of the LES-generated pressures, though it underestimates the amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

The ArborX Library: Version 2.0

This article provides an overview of the 2.0 release of the ArborX library, a performance portable geometric search library based on Kokkos. We describe the major changes in ArborX 2.0 including a new interface for the library to support a wider range of user problems, new search data structures (brute force and distributed), support for user functions to be executed on the results (callbacks), and an expanded set of the supported algorithms (ray tracing and clustering).

GPU↗
DOE OSTI2025

ArborX 2.0

ArborX library tackles a problem of efficiently finding geometric objects that are close in space. Variations of this problem, such as finding the nearest neighbors of a point, or finding all objects within a certain distance, are inherent components of applications in many fields. The data may be large so that solving the problem efficiently may require significant computational resources, such as multiple processors or accelerators such as general purpose GPUs. ArborX' main advantage in its ability to solve large problems efficiently utilizing a combination of distributed and on-node parallelism. ArborX can be run efficiently on a wide variety of hardware, including GPUs from different vendors, which distinguishes it from other available libraries which typically choose only few of these. The other advantage is that it supports both types of user problems: spatial problems (useful for intersections and finding objects within certain distance), and nearest neighbor problems. ArborX also supports flexible interface in its interaction with a user. Particularly, it allows a user to call user's own function on a positive match, a functionality not rarely available in other libraries. ArborX implements construction and traversal algorithms using efficient tree structures, such as bounding volume hierarchy (BVH). At its core, ArborX uses linear BVH for its low construction cost and sufficient quality. ArborX implements both spatial and nearest-neighbor traversal algorithms. ArborX also provides several clustering algorithms (minimum spanning tree, DBSCAN, HDBSCAN*), interpolation using minimum least squares and ray tracing. ArborX is written using C++, and is parallelized using the message passing interface (MPI) for the distributed communication, and the Kokkos library for on-node parallelism. This approach allows ArborX to be run on a wide variety of hardware, from common laptops and desktops to supercomputers while using the same codebase.

Prokopenko, Andrey [Oak Ridge National Laboratory ↗
DOE OSTI2021

Gravitation and Mesh Adaption

The Gravitation and Mesh Adaptation (GaMA) toolbox consists of a set of Matlab classes for modeling the environment around asteroids and comets. A variety of gravitational models and supporting algorithms from are consolidated alongside a custom meshing utility tailored for the application. This allows the user to work within a single streamlined environment to import and manipulate surface definitions, probe the dynamical environment, integrate trajectories, and post-process results. For other applications, modified surface meshes can be exported. Gravity Models: 1) Werner's analytic polyhedron model 2) Mascon model 3) Gottlieb's spherical harmonic model 4) Approximate polyhedron models 5) Curvilinear surface models 6) Custom composite models Meshing Features: 1) Array-based half-edge data structure 2) Supports curvilinear surface definitions up to degree 4 3) Ray-tracing 4) Coarsening 5) Feature-based refinement 6) Projection 7) Smoothing 8) Mesh quality validity tests 9) Mesh repair Additional Features: 1) Solar radiation pressure model 2) Distant 3rd body model 3) Collision detection 4) Trajectory integration and post-processing 5) Visualization of surface fields

Pearl, Jason↗
DOE OSTI2022

ODD (Version 0.X)

The Orthogonal Diffusion Discretizations (ODD) library is a collection of low order (Diffusion, P1, Ray Tracing, and Variable Eddington Factor (VEF)) thermal radiation solvers for orthogonal structured grids driven by an application programing interface (API).

Cleveland, Mathew↗
DOE OSTI2023

Macaw

MaCaw is a MOOSE-based application that enables domain-decomposed neutral particle transport calculations in MOOSE. It leverages the ray tracing MOOSE module for unstructured mesh tracking and OpenMC for collision physics, handling material definitions, and tallying quantities.

Harbour, LoganH.↗
DOE OSTI2025

Direct numerical simulations for hybrid rocket boundary layers: Performance modeling and scaling

This paper presents a comprehensive performance and scaling analysis of direct numerical simulations for reacting boundary layers, focusing on slab burner configurations. Using a PETSc-based finite volume CFD framework, the study evaluates the scalability and computational cost of flow, chemistry, and radiation evaluations across 2D and 3D simulations. Polymethyl methacrylate (PMMA) is the fuel with pure O 2 as the oxidizer, modeled using a detailed chemical kinetics mechanism with 113 species and 660 reactions. A ray-tracing-based radiation solver, designed for distributed memory applications, is implemented to model radiation heat transfer. Parallel scalability is analyzed for the coupled flow, chemistry, and radiation heat transfer processes. Weak and strong scaling studies are conducted on up to 15,000 computational ranks, revealing robust performance when flow cells exceed 200 per rank. Chemistry evaluations dominate the computational cost in large 3D simulations, accounting for approximately 40% of the total runtime, while flow processes contribute around 35%, and radiation solver contributions remain below 10% due to reduced evaluation frequencies. GPU accelerated chemistry evaluation, implemented with Zero-RK, demonstrates significant promise, achieving up to a 4x speedup for workloads exceeding 30,000 cells per GPU. However, diminishing returns are observed for smaller workloads due to CPU-GPU communication overhead. This study identifies key challenges, including memory bottlenecks and the effects of domain partitioning on flow scalability, while highlighting the potential of GPU-accelerated chemistry to reduce computational costs. In conclusion, these findings provide realizable run configurations for 2D, 3D, and GPU-accelerated cases, offering insights for optimizing reactive flow solvers.

CFD Scalability↗
DOE OSTI2023

Sidelobe modeling and mitigation for a three mirror anastigmat cosmic microwave background telescope

Telescopes measuring cosmic microwave background (CMB) polarization on large angular scales require exquisite control of systematic errors to ensure the fidelity of the cosmological results. In particular, far-sidelobe contamination from wide angle scattering is a potentially prominent source of systematic error for large aperture microwave telescopes. Here we describe and demonstrate a ray-tracing-based modeling technique to predict far sidelobes for a three mirror anastigmat telescope designed to observe the CMB from the South Pole. Those sidelobes are produced by light scattered in the receiver optics subsequently interacting with the walls of the surrounding telescope enclosure. After comparing simulated sidelobe maps and angular power spectra for different enclosure wall treatments, we propose a highly scattering surface that would provide more than an order of magnitude reduction in the degree-scale far-sidelobe contrast compared to a typical reflective surface. Finally, we conclude by discussing the fabrication of a prototype scattering wall panel and presenting measurements of its angular scattering profile.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Design of a broadband, high-power, off-axis spiral polarization reflector

Spirally polarized beams are known to produce a strong longitudinal field under high numerical aperture focusing. Here, we report on the design of a reflective optical element with an off-axis geometry that is capable of transforming a uniformly polarized beam into a predominantly spirally polarized beam. By utilizing the retardance present in Fresnel reflection, our device can operate over a broad range of wavelengths. The output polarization and intensity structure of a beam passing through the reflector are calculated using polarization and scalar ray tracing. Vector diffraction is used to calculate the polarization and intensity structure under $f/0.656$ focusing. The off-axis reflective design is well-suited for high-power laser applications and could enable current and future laser-based particle accelerators to achieve strong longitudinal fields at focus.

Black, A. Nicholas [Univ. of Rochester, NY (United↗
DOE OSTI2022

Altering apparent optical properties with an array of semitransparent mesoscale structures

The ability to control and optimize interactions between light and matter has much utility in engineering design. A well-researched way to achieve optical property modulation is via the use of optical metamaterials, which feature sub-wavelength scale surface structures. In this work, an alternative approach for modulating optical properties is presented using a composite surface modified with a periodic array of semitransparent hemispherical shell mesoscale structures which are larger than the incident light wavelength. A ray-tracing simulation approach is used to predict the optical behavior for an arrayed surface. At oblique angles of incidence, significant increases and decreases in apparent absorptance are achieved via the use of optically thick and thin shells, respectively. Additionally, a potential application to solar cells is described with optimal spectral behavior achieved via the use of semitransparent external structures.

Kucuktas, Onur A.↗
DOE OSTI2023

Effects of a nuclear-disturbed environment on electromagnetic wave propagation through the atmosphere

This paper investigates the effects of a nuclear-disturbed environment on the transmission of electromagnetic (EM) waves through the atmosphere. An atmospheric nuclear detonation can produce heightened free electron densities in the surrounding atmosphere that can disrupt EM waves that propagate through the disturbed region. Radiation transport models simulated the ionization and free electron densities created in the atmosphere from a 1 MT detonation at heights of burst of 5 km, 25 km, and 75 km. Recombination rates for the free electrons in the atmosphere were applied, from previous work in the literature, to determine the nuclear-induced electron densities as a function of time and space after the detonation. A ray-tracing algorithm was applied to determine the refraction and reflection of waves propagating in the different nuclear-disturbed environments. The simulation results show that the free electron plasma created from an atmospheric nuclear detonation depend on the height of burst of the weapon, the weapon yield, and the time after detonation. Detonations at higher altitudes produce higher free electron densities for greater durations and over larger ranges. The larger the free electron densities, the greater the impact on EM wavelengths in regards to refraction, reflection, and absorption in the atmosphere. An analysis of modern infrastructure and the effects of nuclear-disturbed atmospheres on different signal wavelengths and systems is discussed.

42 ENGINEERING↗
DOE OSTI2024

Manufacturability-based optical design optimization for advanced Kirkpatrick–Baez X-ray focusing mirrors

The advanced Kirkpatrick–Baez (AKB) mirror setup is an effective and compelling solution to provide stable X-ray nano-focusing for synchrotron radiation or free-electron laser beamlines. We propose an AKB mirror design optimization approach to mitigate the difficulties associated with mirror fabrication by minimizing the total slope ranges of the four curved mirrors while achieving the expected focusing performance. In the optimization, we have considered geometry constraints to ensure the beam acceptance with the required clear aperture, the diffraction-limited focal size with the adequate numerical aperture, and the desired mirror gaps for adjustment and the necessary working distance for the sample stage. Additionally, practical constraints linked to mirror metrology and fabrication, such as mirror length limits and curvature uncertainty in measurement, are taken into account. Furthermore, progressive objective optimization eliminates the need for any initial guess, fully automating the AKB optimization process. This approach facilitates the development of an elegant Wolter-I or Wolter-III type AKB design solution that satisfies these multiple constraints. In cases where constraints cannot be simultaneously satisfied, the optimization results provide valuable insights into areas where trade-offs need to be considered. Simulations with ray tracing and wavefront propagation validate the optimized AKB design showing high tolerance to the beam incident angle.

36 MATERIALS SCIENCE↗
DOE OSTI2026

CHESS 2025: Orthorectified airborne RGB imagery from NEON AOP surveys

This dataset provides Level 1 (L1) and Level 3 (L3) orthorectified Red-Green-Blue (RGB) imagery collected for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). This high-resolution imagery is a photographic record of red, green, and blue visible light from sunlight reflected off of the Earth’s surface. The data comprise full-color images of the ground surface and are primarily intended to provide context to imaging spectroscopy and light detection and ranging (LiDAR) data. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. RGB images were acquired using the PhaseOne IXM-RS150F high-resolution digital camera onboard the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP). The package data include both an L1 product comprising one camera frame per file and an L3 mosaic aligned to the Universal Transverse Mercator (UTM) Zone 13N grid and the World Geodetic System (WGS) 84 projection. Both products are provided in geotif (.tif) format at 0.1 m ground resolution. The bulk of the imagery was collected during the main CHESS field campaign from June 13 to July 15, 2025. Additional images of a portion of the Upper Taylor (UPTA) domain were collected on September 18, 2025, to fill gaps in imagery identified after the main campaign was complete. RGB camera imagery is not radiometrically calibrated, and therefore pixel values should not be exploited for scientific analysis. Pixel values have undergone a manual adjustment to enhance feature identification. The imagery is rigorously geolocated which does allow for reliable geometric information to be retrieved. To generate the orthorectified imagery, the NEON AOP camera captured visible spectrum in red, green, and blue bands. The raw images were then processed using NEON’s camera orthorectification workflow. A boresight calibration flight was made to build a complete camera, distortion, and alignment model. Color balance/white balance and exposure correction were applied to the raw RGB images. The corrected images were orthorectified by ray-tracing image pixels to a lidar-derived digital surface model (DSM) mesh using the refined camera model, outputting orthorectified raster pixels on a regular grid. Flightline-level data were mosaicked by selecting per-pixel contributions from overlapping orthorectified images using line-of-sight (LOS) zenith angle minimization to reduce edge distortions. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2023

SKS Polarization Anomalies Due to the Coriolis Force

The Earth’s Coriolis force has been well–known to impact surface waves and normal modes, which is essential to accurately interpret these waves. However, the Coriolis force on body waves has been assumed to be negligible and mostly ignored. It has been previously shown that the Coriolis force impacts polarizations of shear waves, whereas the wavefronts remain unaffected. We expand on the potential influences of Earth’s Coriolis force on shear–wave polarization measurements by conducting 3D numerical simulations for elastic waves generated by earthquake and explosive sources in a radially symmetric, and 3D mantle and crustal models. The Coriolis force can produce polarization anomalies of mantle shear waves up to 7° and core phases, such as SKS and SKKS, up to 4°. Uncorrected shear–wave polarizations due to the Coriolis force can cause an additional source of error (5°–10° in fast direction, and 0.2–0.3 s delay time depending on the method and seismic phase), inaccurate interpretation of station misalignments, and imprecise estimates of the core–mantle boundary topography. We show how to correct for the Coriolis force on teleseismic shear waves using 1D ray tracing for well–isolated phases. Here, we recommend the use of full waveform simulations to accurately account for earthquake sources parameters, poorly isolated phases that could include interfering phase arrivals within the measurement time window, and the effect of the Coriolis force on the polarizations of shear waves.

58 GEOSCIENCES↗
DOE OSTI2020

A Machine Learning Approach to the Detection of Ghosting Artifacts in Dark Energy Survey Images

Unwanted artifacts that often plague astronomical images arise from a number of sources that include imperfect optics, faulty image sensors, cosmic ray hits, and even airplanes and satellites. One major source that is dicult to avoid is known as ghosting, and is caused by the scattering and internal reflections of light o of the telescope’s mechanical and optical components. Detecting ghosting artifacts eciently in the large cosmological surveys that will acquire petabytes of data can be a daunting task. In this paper, we use data from the Dark Energy Survey to develop, train, and validate a machine learning model to detect ghosts based on convolutional neural networks. The model architecture and training procedure is discussed in detail, and the performance on the training and validation set presented. Testing is performed on data and results are compared with those from a ray-tracing algorithm. A proof-of-principle is demonstrated that shows promise for Stage IV surveys and beyond.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Gamma-Ray Mirrors. Quarterly Progress Report Q3 FY20 April-June 2020

During the third quarter of FY20, the optic design, or Phase II Task 2 of the project, continued. Our effort was focused on verifying 3D ray-tracing simulations representative of the IFEL set-up with cross-code comparison. Simulations for a single reflection Kirkpatrick-Baez optic showed the feasibility of a design that provides 1D focusing at the detector. Following an increase in scope of work in June, we have been actively engaging in preparations to conduct an additional experiment at IFEL before the end of the fiscal year.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗