Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ray-tracing code”

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

Advanced Simulation of ITER Core X-ray Crystal Spectroscopy

X-Ray Simulation Analysis (XRSA) is an analytical ray-tracing mixed code developed specifically for the ITER Core X-Ray Crystal Spectroscopy (XRCS-Core) diagnostic, which employs a dual-reflection configuration incorporating multiple pre-reflectors made of Highly Oriented Pyrolytic Graphite (HOPG) and spherically curved analyzing crystals. The ITER XRCS-Core is designed for high spectral resolution measurement in specific wavelength ranges, including narrow bands around 1.354 Å for W 64+ , 2.19 Å for Xe 51+ , and 2.555 Å for Xe 44+ and Xe 47+ , enabling diagnostic capability across a broad electron temperature range in the ITER plasma. XRSA facilitates efficient simulation of the spectral performance of this complex X-ray spectroscopic system. Recent updates to the XRSA code have incorporated two critical effects: auto-focusing, which specifically applies to HOPG, and polarization. These two effects are particularly important in the dual-reflection configuration used in the ITER XRCS-Core system to provide more accurate modeling results. Here, simulations conducted with the updated code demonstrate that polarization has a substantial impact on the performance of the dual-reflection system. Additionally, the combined influence of polarization and system layout introduces performance variations across channels through the same crystal.

Computer graphics↗
DOE OSTI2022

Predicted performance of a tangential viewing hard x-ray camera for the DIII-D high field side lower hybrid current drive experiment

High field side launch of lower hybrid current drive (LHCD) has improved accessibility and penetration over low field side launch on DIII-D. Simulations predict single pass absorption under a wide range of plasma conditions. Hard x-ray (HXR) measurement of LHCD generated fast electron bremsstrahlung (50–250 keV) will validate wave propagation and absorption. Emissivity profiles are recovered from one-dimensional inversion of HXR brightness to determine LH damping location, fast electron slowing down time, and some indication of the fast electron energy. The camera will be implemented by populating 32 tangential sightlines of the existing Gamma Ray Imager with Kromek SPEARTM Cadmium Zinc Telluride (CZT) detectors sensitive to 10–1000 keV photons with 10 keV energy resolution. Expected count rates allow for <0.5 ms time resolution. Pulses are processed using 50 ns shaping time Cremat CR-200 Gaussian shaping modules and are digitized by 25 MHz D-TACQ ACQ216 digitizers. The performance of the HXR camera is evaluated by comparing predicted fast electron density profiles and inverted synthetic brightnesses obtained from the ray-tracing/Fokker–Planck codes GENRAY/CQL3D. Inversions closely matched predicted fast electron profiles for a range of experimental parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Study of turbulence-induced refraction of lower hybrid waves using synthetic scrape-off layer filaments

Here, turbulence-induced refraction effects to lower hybrid (LH) wave propagation and current drive are studied using synthetic scrape-off layer (SOL) blob/filament fields. A synthetic 3D, field-following, blob turbulence model is implemented in the ray-tracing/Fokker-Planck (RTFP) codes GENRAY/CQL3D. In Alcator C-Mod, the blob field is shown to significantly affect LH ray-trajectories, leading to increased on-axis damping and smoother current profiles. This effect depends on the average blob size and amplitude. In addition, the diffusion of ray-trajectories in phase-space caused by turbulence increases the robustness of the RTFP model. A modified $N_{||}$ launch spectrum, acting as a proxy for parametric decay instability (PDI) effects, is included in simulations with the blob model. A synergy between the modified launch spectrum and turbulence-induced refraction results in synthetic hard x-ray profiles that agree with experiment. Lastly, the blob model is used to predict the effect of SOL turbulence on DIII-D high-field side (HFS) LH launch. Assuming low turbulence amplitude in the HFS SOL (~5%), turbulence-induced refraction is predicted to have little effect on current drive efficiency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Optimization of the N ∥ Upshift in the DIII-D high field side lower hybrid current drive experiment

High field side lower hybrid current drive (LHCD) is one potential candidate for efficient non-inductive current drive in tokamak power plants, and the first test of this technology will occur on the DIII-D tokamak during the 2024 campaign. Previous LFS launch experiments operated in the multi-pass regime and relied on scrape-off layer interactions to close the spectral gap. In the DIII-D experiment, single-pass damping is achievable via an upshift in the parallel refractive index N ∥ caused by mode converting twice (slow → fast → slow). This mode conversion affects the ray trajectories and can lead to enhanced N ∥ upshift depending on where mode conversion occurs. Compared to multi-pass absorption experiments, the optimization of launched N ∥ and plasma parameters can be counter-intuitive: increased density may increase efficiency and smaller N ∥,launch tend to damp closer to the separatrix. A hard x-ray camera installed to measure the bremsstrahlung (50–250 keV) radiation from LHCD-generated fast electrons is capable of verifying the trends reporting in this paper through comparison to the ray-tracing/Fokker–Planck codes GENRAY/CQL3D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design workflow of a symmetric traveling wave antenna for fast ion production on DD tokamaks

Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Parametrizations of coupled betatron motion for strongly coupled lattices

The coupling of transverse motion is a natural occurrence in particle accelerators, either in the form of a residual coupling arising from imperfections or originating by design from strong systematic coupling fields. While the first can be treated perturbatively, the latter requires a robust approach adapted to strongly coupled optics, and a parametrization of the linear optics must be performed to explore beam dynamics in such peculiar lattices. This work highlights the key physical interpretations of the main parametrization formalisms to describe linear coupled optics, along with explicit links and comparisons of these parametrizations. Concepts rarely illustrated in other works, such as forced mode flips and local coupling, are explored in detail, clarifying some anomalies that can arise in lattice functions. The analytical methods have been implemented in a reference Python package and connected with ray-tracing and integration codes to explore examples of strongly coupled lattices, which are discussed in detail to highlight the key physical interpretations of the parametrizations and characteristics of the lattices.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Coupling laser physics to radiation-hydrodynamics

In order to accurately model implosion hydrodynamics in a radiation-hydrodynamics code, it is essential to include accurate accounting for energy deposition physics. In inertial confinement fusion (ICF), where capsules are driven by lasers or laser-driven x-rays, energy deposition profiles and energy transport have a strong impact on the development and evolution of capsule dynamics and hydrodynamic instabilities. Nevertheless, accurately modeling laser beam propagation in radiation-hydrodynamics codes presents unique challenges associated with disparate resolution requirements, the potential to seed spurious noise in highly unstable systems, and computational expense. We discuss a new method for coupling laser ray-tracing physics to a radiation hydrodynamics code, developed in the process of implementing the Mazinisin laser ray-trace into the xRAGE radiation hydrodynamics code. In contrast to previous approaches, in which laser ray-tracing is performed on the radiation-hydrodynamics mesh, our method involves a mesh generation and evolution strategy that addresses the unique requirements of the laser ray-trace in a separate mesh, enabling performance enhancements and strategies to reduce noise seeded by the discretization of beams into computational rays. In addition, we have employed several methods to ensure that spurious mesh imprinting is minimized. These involved optimizing the laser and radiation-hydrodynamics meshes as well as interpolation between them and requires the use of an exact initialization method for the radiation-hydrodynamics mesh. Furthermore, these techniques have enabled efficient computation of laser-driven implosions and other experiments with minimal introduction of spurious noise.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

3D hybrid fluid-particle jet simulations and the importance of synchrotron radiative losses

Context. Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization emission. Aims. Through the use of 3D hybrid fluid-particle jet simulations (with the PLUTO code), we overcome some of the commonly made assumptions in relativistic magnetohydrodynamic (RMHD) simulations by using non-thermal particle attributes to account for the resulting synchrotron radiation. Polarized radiative transfer and ray-tracing (via the RADMC-3D code) highlight the differences in total intensity maps when (i) the jet is simulated purely with the RMHD approach, (ii) a jet tracer is considered in the RMHD approach, and (iii) a hybrid fluid-particle approach is used. The resulting emission maps were compared to the example of the radio galaxy Centaurus A. Methods. We applied the Lagrangian particle module implemented in the latest version of the PLUTO code. This new module contains a state-of-the-art algorithm for modeling diffusive shock acceleration and for accounting for radiative losses in RMHD jet simulations. The module implements the physical postulates missing in RMHD jet simulations by accounting for a cooled ambient medium and strengthening the central jet emission. Results. We find a distinction between the innermost structure of the jet and the back-flowing material by mimicking the radio emission of the Seyfert II radio galaxy Centaurus A when considering an edge-brightened jet with an underlying purely toroidal magnetic field. We demonstrate the necessity of synchrotron cooling as well as the improvements gained when directly accounting for non-thermal synchrotron radiation via non-thermal particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Cathode to Target Simulations for Scorpius: I. Simulation Codes and Models

The new Scorpius linear induction electron accelerator is under development for multi-pulse flash radiography of large, explosively-driven hydrodynamic experiments. Beam physics from the cathode to the target was examined with computer simulations, including envelope, ray-trace, and particle-in cell (PIC) codes. Beam instabilities investigated included beam breakup (BBU), image displacement, diocotron, parametric envelope, ion hose, and the resistive wall instability. Beam corkscrew motion and emittance growth from beam mismatch were also studied. The results of these simulations is documented in a series of reports. In this report the computer codes and physical models used for these simulations are described. The conclusion of this study is that Scorpius will produce and accelerate a beam with radiographic quality equivalent to the present accelerators at Los Alamos National Laboratory if the same engineering standards and construction details are upheld.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Measuring the Stress Factors for Photovoltaic (PV) Backsheet Degradation

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

14 SOLAR ENERGY↗
DOE OSTI2025

A three-dimensional laser ray-tracing methodology for radiation-hydrodynamics simulations

We report on a methodology for performing laser ray-tracing in three spatial dimensions for radiation-hydrodynamics simulation codes. Our method, which is an extension of that developed in Haines et al., Comput. Fluids 201, 104478 (2020), utilizes an automatically generated separate mesh for the laser ray-tracing from the radiation-hydrodynamics mesh. This enables the laser mesh to be tailored to minimize ray noise with significantly fewer rays than would be required when the ray-tracing is performed on the radiation-hydrodynamics mesh, primarily by allowing the use of high-aspect-ratio cells that are not suitable for hydrodynamics solvers. For a planar target, we show that our method provides a ≈ 100× reduction in computational expense to achieve a fixed level of ray noise relative to ray-tracing directly on the radiation-hydrodynamics mesh. The relatively low ray requirement also enables efficient computation of cross-beam energy transfer. Each cell in the logically cubic laser mesh is a non-convex dodecahedron with triangular sides, and numerical integration of the ray trajectories and inverse bremsstrahlung is performed by mapping each cell to the unit cube. We will describe our methodology in detail as well as its implementation in the xRAGE radiation-hydrodynamics code, discuss performance, and present the results from applying the methodology to test problems with analytic solutions for laser ray-tracing through a quadratic density gradient with an analytic solution as well as for a laser-driven heat front. In 3D radiation-hydrodynamics simulations of laser-driven experiments performed on the National Ignition Facility, laser ray-tracing with our methodology uses less than 1% of total computational time while introducing acceptably low levels of ray noise.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Initial conditions for simulations of beam physics in linear induction accelerators

Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the DualAxis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Correct Initial Conditions for Simulations of Beam Physics in Linear Induction Accelerators

Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the Dual Axis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.

43 PARTICLE ACCELERATORS↗
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↗
DOE OSTI2023

Advanced Modeling and Simulation Methods for Evaluation of Thermal Neutron Scattering Materials

With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A ray-tracing algorithm for ab-initio calculation of thermal load in undulator-based synchrotron beamlines

The OASYS suite and its powerful integration features are used to implement a ray-tracing algorithm to accurately calculate the thermal load in any component of an undulator-based synchrotron beamline. This is achieved by sampling and converting the SRW source of a given energy into a Shadow source and using the latter code to ray trace the full beamline. The accuracy of the algorithm is proved by reconstructing the full undulator radiation distribution through an aperture and comparing the result with direct calculaton of the total power using SRW. The algorithm is particularly suited to analyze cases with complex beamline layouts and optical elements, such as crystals, multilayers, and compound refractive lenses. Examples of its use to calculate the power load on elements of two of the feature beamlines at the Advanced Photon Source Upgrade Project and a comparison of the results with analytical calculations are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

New modelling capabilities in IDT

This work concerns the enhanced modelling capabilities of the discrete ordinates transport solver IDT. The novelties introduced allow for modelling unstructured geometries composed by a collection of X/Y segments and circles, and the use of reciprocity and conservation relations reduce the memory imprint as well as the computational cost of the method. IDT decomposes geometries in modular Cartesian patterns, which are the so-called Heterogeneous Cartesian Cells (HCCs), containing a chunk of the original unstructured geometries. Each HCC can be then discretized by superimposing a XY grid to refine locally the HCC. Unlike the most popular MOC, IDT performs the spatial sweeping by directional collision probabilities instead of trajectories. The sources and interface angular fluxes are expanded up to linear order. The accuracy of ray-tracing, the memory imprints together with the novel mesh refinement capabilities have been verified. A first set of preliminary results on PWR lattice problems will be presented in this paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Implementation of Triply Periodic Minimal Surfaces (TPMS) as surface objects in OpenMC

Triply Periodic Minimal Surfaces (TPMS) represent a promising geometry for future fuel designs due to their significant surface-to-volume ratio, which facilitates efficient cooling of nuclear fuel, a crucial factor for safety and efficiency. Demonstrating the remarkable capabilities of TPMS fuel requires initial modeling and simulation. This paper presents an implementation of TPMS in the Monte Carlo code OpenMC, enabling reactor physics modeling of TPMS. Here, the primary advantages over traditional methods using CAD files include reduced memory requirements for computations and high-fidelity implementation. This implementation has been tested against CAD files loaded in Serpent2, yielding promising results with low biases in the $k_{\textrm{eff}}$, comparable to biases in the material balance sheet. The implementation presented in this work will be used in future reactor physics computations related to new reactor designs involving TPMS-based fuels.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗