Search NASA⌕ Search

SEARCH · Search NASA

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Accurate and efficient predictions of keyhole dynamics in laser materials processing using machine learning-aided simulations

The keyhole phenomenon has been widely observed in laser materials processing, including laser welding, remelting, cladding, drilling, and additive manufacturing. Keyhole-induced defects, primarily pores, dramatically affect the performance of final products, impeding the broad use of these laser-based technologies. The formation of these pores is typically associated with the dynamic behavior of the keyhole. So far, the accurate characterization and prediction of keyhole features, particularly keyhole depth, as a function of time, has been a challenging task. In situ characterization of keyhole dynamic behavior using the synchrotron X-ray technique is informative but complicated and expensive. Current simulations are generally hindered by their poor accuracy and generalization abilities in predicting keyhole depths due to the lack of accurate laser absorptance data. In this study, we develop a machine learning-aided simulation method that accurately predicts keyhole dynamics, especially in keyhole depth fluctuations, over a wide range of processing parameters. In two case studies involving titanium and aluminum alloys, we achieve keyhole depth prediction with a mean absolute percentage error of 10 %, surpassing those simulated using the ray-tracing method with an error margin of 30 %, while also reducing computational time. This exceptional fidelity and efficiency empower our model to serve as a cost-effective alternative to synchrotron experiments. Our machine learning-aided simulation method is affordable and readily deployable for a large variety of materials, opening new doors to eliminate or reduce defects for a wide range of laser materials processing techniques.

Computational fluid dynamics↗
DOE OSTI2025

Diagnostics: Chapter 8 of the special issue: on the path to tokamak burning plasma operation

This chapter presents the activity conducted by the ITPA topical group (TG) on Diagnostics over about the last 15 years. Following a general introduction of the ITER Diagnostics led by their measurement roles, the document is organized in several subchapters detailing the design support, research and development activity conducted by each of the specialist working groups (WGs) of the TG. Please note that the magnetic diagnostics were supported at the TG without a specific WG. Their status is included in the general introduction. In the following some highlights of the subchapter’s contents are provided. Recent advances in ITER first wall (FW) diagnostics for the measurements of plasma-metallic wall interaction in support of the ITER research plan are reported. An InfraRed imaging Video Bolometer for ITER has been developed and tested on several tokamaks to measure the radiated power loss. A laser-induced breakdown spectroscopy (LIBS) technique which utilizes a pulsed laser beam to ablate locally by forming a crater, will measure local tritium inventory in the FW material. Real-time Residual Gas Analyzers will measure the neutral gas composition in a divertor port and an equatorial port during plasma operation. Due to the full metallic FW environment, the plasma-wall interaction in ITER will face several challenges such as the compromised radiated power and divertor heat flux measurements by reflection. Ray tracing and analysis codes have been developed to eliminate and correct the effects of reflection in the measurements. The characteristics of the reflecting surfaces depending on the roughness and angle of the incidence have been measured by dedicated experiments, and the results were applied to the reflection elimination. For the measurement of the metallic impurity radiation induced by eroded metallic atoms, a vacuum ultraviolet spectrometer has been developed and tested. An extensive thermonuclear diagnostic suite will be required to support the operation of ITER and the planned experimental program for future burning plasma experiments. Due to the harsh environmental conditions, the implementation of diagnostic systems in ITER is a major challenge. These conditions include high levels of neutron and gamma fluxes, neutron heating, particle bombardment. Therefore, the selection and design of diagnostic systems must take into account a number of phenomena previously unseen in diagnostic design. For this reason, the measurement of neutrons and confined or lost fast ions, with particular emphasis on alpha particles, is critical to ITER. The diagnostics associated with these measurements will be important for future plasma-burning experiments at ITER. The high neutron emission and very large plasma size in ITER make neutron diagnostics the main diagnostic method used to measure plasma parameters such as fusion power, fusion power density, ion temperature, energy of fast ions and their spatial distributions in the plasma core. Active spectroscopy techniques are methods where a neutral particle beam is injected into the plasma and information on plasma parameters is extracted from the measurement of line emission resulting from the beam-plasma interaction, either by plasma ions or by beam atoms. Spatial localization is achieved by crossing the beamline and multiple observation lines. The ITER plasma will be a high temperature, moderately dense, fully ionized collisional plasma. The plasma facing surfaces are principally metallic being fashioned from beryllium or tungsten but many other elements, arising from either structural or from operational needs, may enter this plasma. The energy range of the emitted photons range from meV (infra-red) to multi keV (x-rays) and originate from all areas of the plasma volume. The primary role of passive emission diagnostics is to identify what is in the plasma from spectral signatures. Extracting quantitative information from these measurements such as impurity content, ion temperature, rotation, degree of detachment and radiated power depends on calibrated instruments, a physics model of the atomic and molecular processes and plasma transport and an analysis workflow that takes into account environmental effects such as reflections. The particular needs for ITER have prompted a multi-machine, many-year effort to address all these aspects and this chapter reviews the work on diagnostic design, experiments and new analysis techniques. An overview of the laser diagnostics to be implemented on ITER is also provided in this paper. This includes descriptions of the Thomson scattering in the core, edge and divertor regions, polarimetry and interferometry diagnostics used for measuring plasma density and also measurements of helium density in the divertor using Laser Induced Flourescence. Techniques which can allow improvements on current measurements are also addressed in particular expanding poloidal polarimetry measurements to measure field fluctuations and proposed use of dispersion interferometery which has a number of advantages over existing methods. This paper identifies particular areas where further research and testing on existing tokamaks is useful even at this advanced stage to inform the design of diagnostics for ITER. Outstanding areas of concern for the implementation of laser diagnostics, in particular with a view to reliable operation are identified. An overview of the latest developments of microwave diagnostic systems and techniques is given. The primary focus is the contributions for ITER—the next step burning plasma experiment—which is supplemented by describing recent progress of techniques applicable for fusion experiments beyond ITER. The contributions are intentionally kept concise, and are being supplemented by a rich list of references for further studies. Radiation induced effects are receiving continuous and well-deserved attention of the ITER diagnostic community and they are in many cases one of the primary design drivers of the ITER diagnostic systems. The paper summarizes recent progress in this area focusing primarily on the ITER diagnostics but in some cases provides also outlook for the possible solutions for even more demanding radiation environment of fusion reactors beyond ITER. Despite advancements in the area of modeling and simulation of various radiation induced effects, experimental testing in a nuclear environment as close as possible to the target one is still seen as unavoidable for proper qualification of particular diagnostic functional elements. Recent advancement within three diagnostic areas: optical diagnostics, magnetics and bolometers is covered. Encouraging results on qualification of silica glass vacuum window assemblies are presented. In the area of magnetic sensors, progress of irradiation tests performed on ITER in-vessel LTCC inductive sensors is presented with outlook for novel technological approaches to inductive sensors utilizing thick printing and photolithography technologies being highlighted. Summary of advancements in the area of steady state magnetic field sensors based on Hall effect is given. New results of neutron irradiation test of the ITER borosilicate glass inserts for vacuum electrical feedthroughs are summarized finding negligible swelling at target level of neutron fluence. Off-line irradiation tests of fiber optic current sensors for plasma current measurement demonstrated that both for gamma doses up to 5 MGy and a total neutron fluence up to 10 15 cm −2 , radiation induced changes are still compatible with required measurement accuracy on ITER. The ITER bolometers are given as an example how considering radiation effects may influence the diagnostic design. Finally, outlook for future main R&D directions is outlined. All optical and laser-based diagnostics in ITER will be using mirrors to guide plasma radiation toward detectors, cameras and sensors. In the hostile plasma, radiation and particle environment the optical characteristics of diagnostic mirrors will degrade directly affecting the entire performance of involved diagnostic systems. An assessment of factors affecting mirror performance is provided. Among the prime adverse factors are deposition of plasma impurities, sputtering of mirror surface and steam ingress in the vicinity of mirrors. Within the International Tokamak Physics Activity with active support by ITER central team and domestic agencies, the structured research and development (R&D) program on mitigation of risks for diagnostic mirrors is underway. Within this program the mirror material development, the passive mitigation of mirror degradation by using diagnostic ducts and shutters along with an active mirror recovery program comprising the in-situ mirror cleaning and calibration is underway. Recent developments in diagnostic mirror R&D are described in this Chapter along with an example of their implementation of R&D solutions in ITER Infrared Thermography diagnostic. An assessment of still open engineering and physics questions, considerations on mirror risks during an early phase of ITER operation are given along with an overview of diagnostic mirror evolution in the late ITER operation stage toward the demonstration fusion power plant. Several crucial areas of diagnostic R&D outlined in ITER Research Plan are addressed. The basic control groups in a fusion reactor can be broken-down in five categories: (1) plasma position, magnetic configuration, and plasma current control, (2) profile control and confinement optimization, (3) MHD control and suppression, (4) edge dissipation control, radiation and plasma exhaust control and (5) break-down optimization. These categories are coupled via the physics (a control action in one domain will affect the other domains) and via shared actuators (e.g. ECRH for impurity accumulation avoidance, current density distribution control and MHD suppression). Consequently, a supervisory control system should determine the priority of the various control tasks, their couplings, and the interfaces with the safety and interlock system. For the systematic development of the various controllers taking the complexity of the plasma and the control system into account, a model-based approach is required. A short historical overview is given of the developments in systems and control theory and control engineering with special emphasis on those developments that are most relevant for Nuclear Fusion research and operation. An overview is given of the state of the field of fusion plasma control for the control categories. It will be shown how synthetic diagnostics are being developed in ITER and how they are used in diagnostic design and design validation and how they can be in model-based controller synthesis using relatively simple models. In modern control methods, multiple diagnostics are used to constrain relatively simple models. The constrained models provide an estimate for the state. This opens the route to state controllers, such as model predictive control. A major challenge in nuclear fusion research is the coherent combination of data from heterogeneous diagnostics and modeling codes for machine control and safety as well as physics studies. Measured data from different diagnostics often provide information about the same subset of physical parameters. Additionally, information provided by some diagnostics might be needed for the analysis of other diagnostics. A joint analysis of complementary and redundant data allows, e.g. to improve the reliability of parameter estimation, to increase the spatial and temporal resolution of profiles, to obtain synergistic effects, to consider diagnostics interdependencies and to find and resolve data inconsistencies. Physics-based modeling and parameter relationships provide additional information improving the treatment of ill-posed inversion problems. A coherent combination of all kind of available information within a probabilistic framework allows for improved data analysis results. The concept of integrated data analysis (IDA) in the framework of Bayesian probability theory is outlined and contrasted with conventional data analysis. Components of the probabilistic approach are summarized and specific ingredients beneficial for data analysis at fusion devices are discussed.

ITER↗
NASA NTRS2016

Normalizing a Relativistic Model of X-Ray Reflection Definition of the Reflection Fraction and Its Implementation in relxill

Aims. The only relativistic reflection model that implements a parameter relating the intensity incident on an accretion disk to the observed intensity is relxill. The parameter used in earlier versions of this model, referred to as the reflection strength, is unsatisfactory; it has been superseded by a parameter that provides insight into the accretion geometry, namely the reflection fraction. The reflection fraction is defined as the ratio of the coronal intensity illuminating the disk to the coronal intensity that reaches the observer. Methods. The relxill model combines a general relativistic ray-tracing code and a photoionization code to compute the component of radiation reflected from an accretion that is illuminated by an external source. The reflection fraction is a particularly important parameter for relativistic models with well-defined geometry, such as the lamp post model, which is a focus of this paper. Results. Relativistic spectra are compared for three inclinations and for four values of the key parameter of the lamp post model,namely the height above the black hole of the illuminating, on-axis point source. In all cases, the strongest reflection is produced for low source heights and high spin. A low-spin black hole is shown to be incapable of producing enhanced relativistic reflection. Results for the relxill model are compared to those obtained with other models and a Monte Carlo simulation. Conclusions. Fitting data by using the relxill model and the recently implemented reflection fraction, the geometry of a system can be constrained. The reflection-fraction is independent of system parameters such as inclination and black hole spin. The reflection-fraction parameter was implemented with the name reflec_frac all flavours of the relxill model, and the non-relativistic reflection model xillver, in v0.4a (18 January 2016).

accretion↗
NASA NTRS2019

A NICER view of PSR J0030+0541: Millisecond Pulsar Parameter Estimation

We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSRJ0030+0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray spectral timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar’s surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass M and equatorial radius Req are, respectively, -1.34+ M 0.16 0.15 and-12.71+ km 1.19 1.14 , while the compactness = -GM R c 0.156+ eq 2 0.010 0.008 is more tightly constrained; the credible interval bounds reported here are approximately the 16% and 84% quantiles in marginal posterior mass.

T E Riley↗
DOE OSTI2025

Implementation of compound refractive lenses for large field-of-view x-ray phase-contrast imaging during hypervelocity impact experiments

Synchrotron x-ray phase-contrast imaging (XPCI) offers time-resolved visualization of dynamic compression phenomena, but its intrinsically small field-of-view (FOV) limits the time that key features remain in frame. A novel approach to enlarge the FOV is achieved by positioning a two-dimensional parabolic compound refractive lens (CRL) upstream of the sample to deliberately defocus the white beam. Ray-tracing simulations and XPCI measurements show that this CRL configuration can expand the beam by ∼50% vertically and ∼15% horizontally based on the full width at half-maximum of the beam. Implementing the CRL, however, attenuates the photon flux and lowers signal-to-noise ratio (SNR). Task-based analysis using a calibration grid (30 μm dots) showed that both setups fail to consistently meet the Rose criterion (SNR ≥ 5) for features of this size in single-bunch imaging. Extrapolating the measured SNR Rose values suggests that the minimum consistently detectable feature lies closer to 30–40 μm for the standard XPCI setup and above 40 μm for CRL-XPCI. Despite this limitation, the CRL configuration nearly doubles the illuminated area, enabling simultaneous tracking of front and rear observations of boron carbide targets subjected to rod and sphere impacts at 1.0–2.6 km/s. Image tracking algorithms and photonic Doppler velocimetry were used to measure penetration and rear-surface velocity histories. Together, these measurements capture crack fronts, penetration, and material breakout, offering new benchmark data for validating high-strain-rate constitutive models of ceramic materials.

Ceramic materials↗
DOE OSTI2024

Artificial ground reflector size and position effects on energy yield and economics of single‐axis‐tracked bifacial photovoltaics

Abstract Artificial ground reflectors improve bifacial energy yield by increasing both front and rear‐incident irradiance. Studies have demonstrated an increase in energy yield due to the addition of artificial reflectors; however, they have not addressed the effect of varying reflector dimensions and placement on system performance and the impact of these parameters on the reflectors' financial viability. We studied the effect of high albedo (70% reflective) artificial reflectors on single‐axis‐tracked bifacial photovoltaic systems through ray‐trace modeling and field measurements. In the field, we tested a range of reflector configurations by varying reflector size and placement and demonstrated that reflectors increased daily energy yield up to 6.2% relative to natural albedo for PERC modules. To confirm the accuracy of our model, we compared modeled and measured power and found a root mean square error (RMSE) of 5.4% on an hourly basis. We modeled a typical meteorological year in Golden, Colorado, to demonstrate the effects of artificial reflectors under a wide range of operating conditions. Seventy percent reflective material can increase total incident irradiance by 1.9%–8.6% and total energy yield by 0.9%–4.5% annually after clipping is considered with a DC–AC ratio of 1.2. Clipping has a significant effect on reflector impact and must be included when assessing reflector viability because it reduces reflector energy gain. We calculated a maximum viable cost for these improvements of up to $2.50–4.60/m 2 , including both material and installation, in Golden. We expanded our analysis to cover a latitude range of 32–48°N and demonstrated that higher‐latitude installations with lower energy yield and higher diffuse irradiance content can support higher reflector costs. In both modeling and field tests, and for all locations, the ideal placement of the reflectors was found to be directly underneath the module due to the optimized rear irradiance increase.

14 SOLAR ENERGY↗
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 ↗