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

Benchmarking Cylindrical Blast Wave Theory Against the OSIRIS-REx Sample Return Capsule Reentry

Weak shock theory based on cylindrical blast waves has been used to interpret meteor infrasound, but it has not been systematically benchmarked against a non-ablating hypersonic source with independently known parameters. The objective of this study is not to propose a new theoretical framework, but to evaluate the operational validity of the existing suite of blast radius formulations against a high-fidelity ground truth dataset. The OSIRIS-REx Sample Return Capsule reentry on 24 September 2023 provides such a benchmark because the capsule geometry, trajectory, and infrasound emission points are constrained from mission data and ray tracing, reducing source-side uncertainty associated with ablation. Using observations from 39 infrasound stations, this benchmarking study evaluates six published blast radius (${R}_{0}$ ) formulations and three weak-shock transition coefficients ( C ) within a stratified atmospheric propagation model to predict signal period and peak overpressure. The benchmarking identifies the Sakurai formulation as the best-performing formulation for non-ablating bodies, with the Jones/Plooster formulation performing comparably when a physically appropriate C is adopted. Sakurai and Jones/Plooster yield linear-period median absolute percentage residuals of 9% and 11%, respectively. The period predictions show only weak sensitivity to C at these propagation distances. The Mach-diameter approximation commonly used in meteor studies overestimates ${R}_{0}$ by more than a factor of 3 in the absence of ablation. Finally, these results establish a performance baseline for applying cylindrical blast wave theory to effectively non-ablating hypersonic bodies and demonstrate that the signal period is a robust observable for constraining ${R}_{0}$.

atmospheric entry↗
DOE OSTI2020

Uncertainty amplification due to density/refractive index gradients in background-oriented schlieren experiments

Here, we theoretically analyze the effect of density/refractive index gradients on the measurement precision of background-oriented schlieren (BOS) experiments by deriving the Cramer–Rao lower bound (CRLB) for the 2D centroid estimation process. A model is derived for the diffraction limited image of a dot viewed through a medium containing density gradients that includes the effect of the experimental parameters such as the magnification and f-number. It is shown using the model that nonlinearities in the density gradient field lead to blurring of the dot image. This blurring amplifies the effect of image noise on the centroid estimation process, leading to an increase in the CRLB and a decrease in the measurement precision. The ratio of position uncertainties of a dot in the reference and gradient images is shown to be a function of the ratio of the dot diameters and dot intensities. We termed this parameter the amplification ratio (A F ), and a methodology for reporting position uncertainties in tracking-based BOS measurements is proposed. The theoretical predictions of the dot position estimation variance from the CRLB are compared to ray tracing simulations, and agreement is obtained. The uncertainty amplification is also demonstrated on experimental BOS images of flow induced by a spark discharge, where it is seen that regions of high amplification ratio correspond to regions of density gradients. This analysis elucidates the dependence of the position uncertainty on density and refractive index gradient-induced distortion parameters, provides a methodology for accounting its effect on uncertainty quantification and provides a framework for optimizing experiment design.

42 ENGINEERING↗
DOE OSTI2020

Potential annual daylighting performance of a high-efficiency daylight redirecting slat system

While the primary role of window attachments is often to moderate glare and solar heat gains, they are also able to provide additional daylight to interior spaces. For this purpose, a variety of daylight-redirecting window systems have been developed over the past 150 years. Fixed reflective systems (slats/light shelves) or prismatic systems that rely on total internal reflection work well under specific solar conditions, but generally sacrifice performance over a much wider range of incident solar angles and sky conditions. Dynamic systems - typically reflective slats - are more responsive to sun angles but have not been able to achieve optimal performance for glare and daylight redirection efficiency. A previous investigation into an adjustable, reflective blind concept first conceived of in the late 1970s showed promise but was not reduced to practice due to lack of adequate simulation and analysis tools. In this paper, this concept is further developed and its energy and visual comfort performance evaluated for four mid-latitude, temperate climates using ray-tracing simulation techniques. Overall, results indicate significant potential lighting energy savings when compared with conventional automated reflective blinds (2.1–4.9 kWh/(m 2 ·a), or 14%–42%, depending on climate and orientation) or, especially, manually-operated matte white venetian blinds (1.4–7.9 kWh/(m 2 ·a), or 9%–54%, depending on climate and orientation), while maintaining acceptable or better visual comfort conditions throughout the interior space.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Self-aligned concentrating immersion-lens arrays for patterning and efficiency recovery in scaffold-reinforced perovskite solar cells

In contrast to existing silicon, CIGS, and multi-junction cells that exhibit remarkable mechanical durability, perovskite solar cells have been shown to delaminate when subjected to the mechanical loads that occur during processing and field exposures, limiting their potential as a reliable solar technology. Mechanical reinforcement is thus essential for durable and reliable perovskite technologies. One current strategy to overcome the thermomechanical fragility of perovskite solar cells is to extrinsically shield them by introducing reinforcing scaffolds which partition the cell into many distinct microcells, but improvements in mechanical stability coincide with a reduced device efficiency due to parasitic absorption by the scaffold. We address this reduced efficiency by integrating concentrating immersion-lens arrays (CILAs) into scaffold-reinforced solar cells. These scaffolds are lithographically patterned by a maskless process, in which ultraviolet light is used to pattern the scaffolds through the CILAs, ensuring self-alignment and optical contact with the microcells. Perovskite devices are deposited into the microcells, and during operation, light is concentrated into the microcells and away from the insulating scaffolds, resulting in mechanically resilient solar cells with efficiencies comparable to planar devices. The CILAs also exhibit passive tracking of incident light—verified experimentally and with ray-tracing simulations—which is critical for optimal power output as the sun moves across the sky during the day. The simplicity of the fabrication process and the efficiency of the resulting scaffolded devices show that a lens-integrated scaffold-reinforced structure is a potential pathway to high-performance, robust, commercially viable perovskite solar cells.

14 SOLAR ENERGY↗
DOE OSTI2025

Thermal Analysis of a Solid Particle Light-Trapping Planar Cavity Receiver Using Computational Fluid Dynamics

Concentrated solar power (CSP) is one of the most effective ways of harnessing solar power to create efficient, durable, and resilient energy systems. This study entails thermal modeling and analysis of a novel central tower receiver configuration. This receiver uses solid particles as the heat transfer fluid (HTF), a promising option for third-generation CSP systems. The configuration considered here is the light-trapping planar cavity receiver (LTPCR) introduced by the National Renewable Energy Laboratory. While heat transfer studies of various LTPCR subsystems have been done, system-level thermal analysis of the LTPCR receiver has not been attempted. This study also presents important sensitivity analyses of the operating parameters of the CSP system, which can help guide the design of future central tower receivers. This study employs Ansys Fluent as a computational fluid dynamics (CFD) tool to model fluid dynamics and heat transfer in the receiver, intending to quantify its thermal performance. The model seamlessly integrates Monte Carlo ray tracing data, which generates absorbed solar flux profiles from the heliostat field design, with the heat transfer characteristics of the fluidized particle bed. This unified model is designed to accurately predict the thermal behavior of the LTPCR. Analysis of preliminary results reveals that the primary loss mechanisms are radiative and natural convective losses, in that order. Based on observations from a baseline case, several strategies are suggested and numerically tested. These solutions include selective cooling of high-temperature regions and manipulation of particle bed parameters. Selective cooling of high-temperature regions reduced the peak temperature by 151 degrees C and decreased thermal losses by 0.9%. Improving the particle-wall heat transfer coefficient (P-W HTC) of the particle bed decreased the thermal losses by 1.7% and decreased the peak temperatures by 57 degrees C. Decreasing the particle inlet temperature (PIT) also reduced thermal losses by 3.5% and decreased peak temperatures by 29 degrees C. Compounding these strategies improved the thermal losses of the receiver from 13.5% in the baseline case to 7.5%. Additionally, the study explores the variation in thermal performance across different locations of the receiver, where a variation of thermal losses from 12.9% to 17.3% is found. This allows a comprehensive evaluation of potential improvements in efficiency and temperature management.

computational fluid dynamics↗
DOE OSTI2026

Optical diagnostic design for measuring the radiation front with a mid-leg pumped divertor on DIII-D

The design of an optical diagnostic system to localize the radiation front in a mid-leg pumped divertor configuration on DIII-D is presented. Divertor detachment is a key mechanism for handling power exhaust in tokamaks, and mid-leg pumping offers a promising approach to radiatively dissipate power while maintaining acceptable core performance. To predict the location of the radiation front and design a spectroscopic diagnostic to measure its position, a database of SOLPS-ITER simulations across a range of input powers and gas puffing rates representative of DIII-D operation was employed. These simulations provide self-consistent plasma backgrounds for a Cherab-Raysect synthetic diagnostic framework, which incorporates detailed tokamak geometry and physically accurate ray-tracing. Synthetic simulations of optical sightlines and viewing cone geometry were used to calculate line-integrated emission from the plasma, particularly of the C III 465 nm line, to serve as a proxy for the T e ≈ 7 − 10 eV temperature region associated with the onset and evolution of divertor detachment. Neutral deuterium emission is also evaluated for comparison. The synthetic diagnostic is used to assess the performance of existing DIII-D optical systems, including filterscopes and the Multichord Divertor Spectrometer, and to optimize line-of-sight placement within mechanical and installation constraints. The results provide quantitative guidance for diagnostic implementation on DIII-D and demonstrate the advantage of integrated synthetic diagnostics for divertor design studies and future advanced divertor concepts.

Cherab↗
DOE OSTI2023

A numerical study of the spectral radiative properties of packed bed with mixed bauxite and silica spheres

Bauxite and silica particles have gained increasing attention for applications in the field of concentrated solar power. In this work, a Monte Carlo ray-tracing simulation is performed to predict the radiative properties (absorptance, reflectance, and transmittance) of packed beds with mixed bauxite and silica spherical particles at wavelengths of 0.5, 2.6, 9.0, and 9.35 µm. These wavelengths are representative for the visible, near-infrared, and the mid-infrared regions that are important for solar and thermal radiation, respectively. A repeating unit-column approach is used to mathematically represent the particle bed. The effects of particle mixing ratio, volume fraction, and wavelength on the predicted radiative properties are examined. The obtained radiative properties are inputted to an inverse method to retrieve the effective absorption and scattering coefficients as well as the scattering albedo, which may be later used in a continuous-scale radiative heat transfer analysis. Furthermore, the independent scattering model is used to obtain the absorption coefficient and scattering albedo based on the absorption and scattering cross sections predicted by a Monte Carlo algorithm for a single particle. It is shown that the independent scattering model underpredicts the scattering coefficient for opaque particles but overpredicts the scattering coefficient for semitransparent particles for sufficiently high particle volume fractions. Here, the radiative properties calculated from the independent scattering model are compared to the full Monte Carlo simulation of the particle bed to examine the influence of particle mixing on dependent scattering.

42 ENGINEERING↗
DOE OSTI2021

A novel energy balance approach for a verifiable and accurate solution of radiation extinction in purely absorbing particle clouds

Here we consider the problem of radiation transport through purely absorbing particle clouds. The gold-standard solution of particle-resolved Monte Carlo ray-tracing method to this problem is computationally expensive and therefore solving the radiation transport equation (specifically, the Beer-Bouguer law (BB-law)) on a Eulerian mesh is often preferred. While the absorption coefficient in the real problem is infinite, the BB-law approximates it to be a finite number in the form of number density through a set of assumptions. For particle clouds that do not obey these assumptions, the BB-law predicts an incorrect exponential decay. Also, when the number density is computed using the nearest-neighbor approach, the BB-law solution diverges when the Eulerian mesh size becomes closer to or smaller than the particle size. This numerical divergence is due to the homogenization error. Although the filtering strategy for number density minimizes the homogenization error, it still converges to an incorrect exponential decay for particle clouds that break the BB-law constraints and the cost of the filtering is equivalent to that of the gold-standard solution. In this study, we develop a novel, highly accurate, verifiable and cost-effective solution to the radiation transport equation on a Eulerian domain using an energy balance approach where we derive an expression for the absorption coefficient as a function of particle and Eulerian mesh sizes. We apply our new method to Poisson and turbulent particle clouds that violate all the BB-law constraints and show that the solution is converging upon the mesh refinement, and eventually, we recover the same gold-standard solution for a much cheaper computational cost.

42 ENGINEERING↗
DOE OSTI2020

X-ray topography characterization of gallium nitride substrates for power device development

Gallium nitride substrates grown by the hydride vapor phase epitaxy (HVPE) method using a patterned growth process have been characterized by synchrotron monochromatic beam X-ray topography in the grazing incidence geometry. Images reveal a starkly heterogeneous distribution of dislocations with areas as large as 0 3 mm 2 containing threading dislocation densities below 10 3 cm -2 in between a grid of strain centers with higher threading dislocation densities (> 10 4 cm -2 ). Basal plane dislocation densities in these areas are as low as 104 cm -2 . By comparing the recorded images of dislocations with ray tracing simulations of expected dislocations in GaN, the Burgers vectors of the dislocations have been determined. The distribution of threading screw/mixed dislocations (TSDs/TMDs), threading edge dislocations (TEDs) and basal plane dislocations (BPDs) is discussed with implications for fabrication of power devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Design of a neutron microscope based on Wolter mirrors

The predominant geometry for a neutron imaging experiment is that of a pinhole camera. Here, this is primarily due to the difficulty in focusing neutrons due to the weak refractive index, which is also strongly chromatic. Proof of concept experiments demonstrated that neutron image forming lenses based on reflective Wolter mirrors can produce quantitative, high spatial resolution neutron images while also increasing the time resolution compared to the conventional pinhole camera geometry. Motivated by these results, we report the design of a neutron microscope where two Wolter mirrors replace condensing and objective lenses, in direct analogy with typical visible light microscopes. Ray tracing results indicate that this system will yield 3 µm spatial resolution images with an acquisition time of order <1 s (10 4 faster than currently possible at this spatial resolution) with a field of view of about 5 mm in diameter.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

A realistic guide misalignment model for the Second Target Station instruments at the Spallation Neutron Source

The Second Target Station (STS) at the Spallation Neutron Source, including a suite of initial instrument concepts, are being designed in Oak Ridge. The moderator size at STS is significantly smaller than that of the First Target Station, making the guide misalignment a more prominent concern. Here a realistic misalignment model based on Monte Carlo neutron ray tracing simulation was developed to evaluate the impacts of both the random misalignment and floor settlement on the performance of a toy model of 150 m straight guide system, and MENUS, an STS instrument concept with a typical guide design of two opposing Montel mirrors. The study shows that longer instruments, and instruments with complex curved mirrors/guides are sensitive to significant performance penalties due to misalignments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Multiscale modeling and cinematic visualization of photosynthetic energy conversion processes from electronic to cell scales

We report conversion of sunlight into chemical energy, namely photosynthesis, is the primary energy source of life on Earth. A visualization depicting this process, based on multiscale computational models from electronic to cell scales, is presented in the form of an excerpt from the fulldome show Birth of Planet Earth. This accessible visual narrative shows a lay audience, including children, how the energy of sunlight is captured, converted, and stored through a chain of proteins to power living cells. The visualization is the result of a multi-year collaboration among biophysicists, visualization scientists, and artists, which, in turn, is based on a decade-long experimental-computational collaboration on structural and functional modeling that produced an atomic detail description of a bacterial bioenergetic organelle, the chromatophore. Software advancements necessitated by this project have led to significant performance and feature advances, including hardware-accelerated cinematic ray tracing and instanced visualizations for efficient cell-scale modeling. The energy conversion steps depicted feature an integration of function from electronic to cell levels, spanning nearly 12 orders of magnitude in time scales. This atomic detail description uniquely enables a modern retelling of one of humanity’s earliest stories—the interplay between light and life.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

First Direct Observations of Propagation of Discrete Chorus Elements From the Equatorial Source to Higher Latitudes, Using the Van Allen Probes and Arase Satellites

Whistler mode chorus waves have recently been established as the most likely candidate for scattering relativistic electrons to produce the electron microbursts observed by low altitude satellites and balloons. These waves would have to propagate from the equatorial source region to significantly higher magnetic latitude in order to scatter electrons of these relativistic energies. This theoretically proposed propagation has never been directly observed. We present the first direct observations of the same discrete rising tone chorus elements propagating from a near equatorial (Van Allen Probes) to an off-equatorial (Arase) satellite. The chorus is observed first on the more equatorial satellite and is found to be more oblique and significantly attenuated at the off-equatorial satellite. Further, this is consistent with the prevailing theory of chorus propagation and with the idea that chorus must propagate from the equatorial source region to higher latitudes. Ray tracing of chorus at the observed frequencies confirms that these elements could be generated parallel to the field at the equator, and propagate through the medium unducted to Van Allen Probes A and then to Arase with the observed time delay, and have the observed obliquity and intensity at each satellite.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Influence of the Kuroshio Intrusion on Deep Flow Intraseasonal Variability in the Northern South China Sea

Abstract Interactions between the open ocean and marginal seas have been suggested to be critical to the redistribution and dissipation of global energy. Here, we propose a mechanism for the upper open ocean influencing the deep flow in marginal seas that hinges on the formation and propagation of topographic Rossby waves (TRWs). Observations and high‐resolution simulations suggest substantial intraseasonal variability with periods of 5–60 days associated with the deep flow over continental slopes of the northern South China Sea (NSCS). These fluctuations generally account for over 40% of the total deep kinetic energy variability and the number can reach 70% over the slopes to the west of the Luzon Strait, southwest of the Dongsha Islands, and northeast of the Xisha Islands. By examining the spatiotemporal features of the fluctuations, we demonstrate that the intraseasonal variability of the deep flow in the NSCS is closely associated with TRWs. Using a recently developed multiscale energetics analysis in combination with a ray tracing model, we show that the energy sources of TRWs can be traced back to the east of the Dongsha Islands, where the Kuroshio intrusion and related eddies energize the TRWs primarily through pressure work. These waves propagate westward across the NSCS and drive the intraseasonal variability of the deep flow over continental slopes. Our findings highlight an energy pathway from the open ocean western boundary current to the abyssal marginal sea that could modulate regional circulation as well as exchanges between major ocean basins.

Quan, Qi↗
DOE OSTI2022

Identification of Acoustic Wave Signatures in the Ionosphere From Conventional Surface Explosions Using MF/HF Doppler Sounding

Here, we present an experiment to detect one ton TNT-equivalent chemical explosions using pulsed Doppler radar observations of isodensity layers in the ionospheric E region during two campaigns. The first campaign, conducted on 15 October 2019, produced potential detections of all three shots. The detections closely resemble the temporal and spectral properties predicted using the InfraGA ray tracing and weakly nonlinear waveform propagation model. Here the model predicts that within 6.5–7.25 min of each shot a waveform peaking between 0.9 and 0.4 Hz will impact the ionosphere at 100 km. As the waves pass through this region, they will imprint their signal on an isodensity layer, which is detectable using a Doppler radar operating at the plasma frequency of the isodensity. Within the time windows of each of the three shots in the first campaign, we detect enhanced wave activity peaking near 0.5 Hz. These waves were imprinted on the Doppler signal probing an isodensity layer at 2.785 MHz near 100 km altitude. Despite these detections, the method appears to be unreliable as none of the six shots from the second campaign, conducted on 10 July 2020 were detected. The observations from this campaign were characterized by an increased acoustic noise environment in the microbarom band and persistent scintillation on the radar returns. These effects obscured any detectable signal from these shots and the baseline noise was well above the detection levels of the first campaign.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Off-harmonic optical probing of high intensity laser plasma expansion dynamics in solid density hydrogen jets

Due to the non-linear nature of relativistic laser induced plasma processes, the development of laser-plasma accelerators requires precise numerical modeling. Especially high intensity laser-solid interactions are sensitive to the temporal laser rising edge and the predictive capability of simulations suffers from incomplete information on the plasma state at the onset of the relativistic interaction. Experimental diagnostics utilizing ultra-fast optical backlighters can help to ease this challenge by providing temporally resolved inside into the plasma density evolution. We present the successful implementation of an off-harmonic optical probe laser setup to investigate the interaction of a high-intensity laser at 5.4 x 10 21 W/cm 2 peak intensity with a solid-density cylindrical cryogenic hydrogen jet target of 5μm diameter as a target test bed. The temporal synchronization of pump and probe laser, spectral filtering and spectrally resolved data of the parasitic plasma self-emission are discussed. The probing technique mitigates detector saturation by self-emission and allowed to record a temporal scan of shadowgraphy data revealing details of the target ionization and expansion dynamics that were so far not accessible for the given laser intensity. Plasma expansion speeds of up to (2.3 ± 0.4)x 10 7 m/s followed by full target transparency at 1.4ps after the high intensity laser peak are observed. A three dimensional particle-in-cell simulation initiated with the diagnosed target pre-expansion at -0.2ps and post processed by ray tracing simulations supports the experimental observations and demonstrates the capability of time resolved optical diagnostics to provide quantitative input and feedback to the numerical treatment within the time frame of the relativistic laser-plasma interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Optical analogues to the equatorial Kerr–Newman black hole

Optical analogues to black holes allow the investigation of general relativity in a laboratory setting. Previous works have considered analogues to Schwarzschild black holes in an isotropic coordinate system; the major drawback is that required material properties diverge at the horizon. We present the dielectric permittivity and permeability tensors that exactly reproduce the equatorial Kerr–Newman metric, as well as the gradient-index material that reproduces equatorial Kerr–Newman null geodesics. Importantly, the radial profile of the scalar refractive index is finite along all trajectories except at the point of rotation reversal for counter-rotating geodesics. Construction of these analogues is feasible with available ordinary materials. A finite-difference frequency-domain solver of Maxwell’s equations is used to simulate light trajectories around a variety of Kerr–Newman black holes. For reasonably sized experimental systems, ray tracing confirms that null geodesics can be well-approximated in the lab, even when allowing for imperfect construction and experimental error.

79 ASTRONOMY AND ASTROPHYSICS↗