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 55 records · Page 3

Signal Whisperers: Enhancing Wireless Reception Using DRL-Guided Reflector Arrays

This paper presents a multi-agent reinforcement learning (MARL) approach for controlling adjustable metallic reflector arrays to enhance wireless signal reception in non-line-of-sight (NLOS) scenarios. Unlike conventional reconfigurable intelligent surfaces (RIS) that require complex channel estimation, our system employs a centralized training with decentralized execution (CTDE) paradigm where individual agents corresponding to reflector segments autonomously optimize reflector element orientation in three-dimensional space using spatial intelligence based on user location information. Through extensive ray-tracing simulations with dynamic user mobility, the proposed multi-agent beam-focusing framework demonstrates substantial performance improvements over single-agent reinforcement learning baselines, while maintaining rapid adaptation to user movement within one simulation step. Comprehensive evaluation across varying user densities and reflector configurations validates system scalability and robustness. The results demonstrate the potential of learning-based approaches for adaptive wireless propagation control.

deep reinforcement learning↗

CHESS 2025: Orthorectified airborne RGB imagery from NEON AOP surveys

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

2018 NEON and 2025 CHESS Campaigns↗

Developing nonlinear laser-plasma instability models for high-fidelity, multi-physics simulation capability for ICF/HED (IC Report) [Slides]

Modeling nonlinear effects of laser plasma instabilities (LPI) is critically important for ICF/HED experiments. Our project has the following goals: • Develop physics-based nonlinear LPI models using PIC simulations • Couple the nonlinear LPI effects to macroscopic modeling of ICF/HED exp. through laser ray-tracing (LRT).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗

Simulations of Radiation Hydrodynamics in Strongly Coupled Plasmas (Final Report)

Radiation propagates throughout the Universe and is often the only information we have from distant astronomical objects. Additionally, radiation plays a key role in systems ranging from the interior of stars to the establishment of shocks at the edges of forming galaxies to neutron star mergers. Here radiative processes, absorption, emission, and scattering, are fundamental to transport of momentum and energy, as well as, the overall evolution of these and other astrophysical systems. In some cases, radiation-dominant systems are also sufficiently dense and ionized, such that, the ions are strongly coupled (ie electrostatic energy >> thermal energy or Γ >> 1), which will affect radiative processes and the transport of radiation. During this project, we developed a computational ray-tracing algorithm to study radiation transport through a simulated neutron star capsule implosion under conditions relevant to Omega-60 laser experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Frequency-Selectable Laser Source (FLS) for Cosmic Microwave Background Experiments

Cosmic Microwave Background (CMB) experiments measure remnant radiation from the early universe and use that data to determine fundamental properties of the universe. We can constrain key parameters, such as \textit{r}, the cosmic tensor-to-scalar ratio, and $N_\text{eff}$, the effective number of relativistic species, by analyzing the CMB power spectra. Improving our measurements of the CMB requires improving our instrument systematics, one of the most important of which is detector bandpass. Current experiments use a Fourier Transform Spectrometer (FTS) to measure bandpass. However, the FTS is systematics limited, and cannot achieve the accuracy needed to make improved CMB measurements. For this reason, we are developing a new instrument, the Frequency-Selectable Laser Source (FLS) to decrease the uncertainty in bandpass by an order of magnitude. In this paper, we describe work completed to support the version 2 upgrade to the FLS. Using ray-tracing software, we modeled the FLS optics to set physical tolerances for the new design. We also discuss the laser calibration, and future work to be completed in further development of the FLS upgrade.

Rosen-Turits, Gabriel M.↗

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↗

Characterizing the Impact of Heliostat Size, Focus Method, and Optical Error on Concentrating Solar Tower Systems

Concentrating Solar Power (CSP) tower systems use heliostat fields to direct solar energy to a central receiver. The efficiency of the heliostat field is crucial for achieving thermal generation targets and is also the largest cost component, making up 30-40% of total project costs. Developers face a tradeoff: reducing heliostat costs may decrease optical performance and increase land use, while improving performance could raise manufacturing costs (through higher precision). The optimal design depends on land, labor, and manufacturing costs. Additionally, inherent performance limitations exist - optical errors result in larger solar images as heliostats move farther from the receiver, increasing spillage and attenuation losses. Ultimately, there exists a point where receiver design power can no longer be achieved at design conditions for a given solar field area. This point depends on key heliostat field design parameters including field layout, heliostats size, focus method, and heliostat optical error. In this work, we explore and characterize the performance trade-offs between heliostat size, focus method, and optical error using Monte Carlo ray-tracing simulations using NREL's high performance computing resources. The goal of this fundamental investigation is to understand practical heliostat design performance impacts that can be used to produce more cost-effective heliostat fields while still achieving plant thermal generation targets.

14 SOLAR ENERGY↗

Systematic Effects in Galaxy-Galaxy Lensing with DESI

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

79 ASTRONOMY AND ASTROPHYSICS↗

Comparison of Polarized Radiative Transfer Codes Used by the EHT Collaboration

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT Collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I , Q , U , and V , respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Buzzard to Cardinal: Improved Mock Catalogs for Large Galaxy Surveys

Abstract We present the Cardinal mock galaxy catalogs, a new version of the Buzzard simulation that has been updated to support ongoing and future cosmological surveys, including the Dark Energy Survey (DES), DESI, and LSST. These catalogs are based on a one-quarter sky simulation populated with galaxies out to a redshift of z = 2.35 to a depth of m r = 27. Compared to the Buzzard mocks, the Cardinal mocks include an updated subhalo abundance matching model that considers orphan galaxies and includes mass-dependent scatter between galaxy luminosity and halo properties. This model can simultaneously fit galaxy clustering and group–galaxy cross-correlations measured in three different luminosity threshold samples. The Cardinal mocks also feature a new color assignment model that can simultaneously fit color-dependent galaxy clustering in three different luminosity bins. We have developed an algorithm that uses photometric data to further improve the color assignment model and have also developed a novel method to improve small-scale lensing below the ray-tracing resolution. These improvements enable the Cardinal mocks to accurately reproduce the abundance of galaxy clusters and the properties of lens galaxies in the DES data. As such, these simulations will be a valuable tool for future cosmological analyses based on large sky surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Stochastic Averaging of Radiative Transfer Coefficients for Relativistic Electrons

Abstract Synchrotron emissivities, absorptivities, and Faraday rotation and conversion coefficients are needed in modeling a variety of astrophysical sources, including Event Horizon Telescope (EHT) sources. We develop a method for estimating transfer coefficients that exploits their linear dependence on the electron distribution function, decomposing the distribution function into a sum of parts each of whose emissivity can be calculated easily. We refer to this procedure as stochastic averaging and apply it in two contexts. First, we use it to estimate the emissivity of an isotropicκdistribution function with a high-energy cutoff. The resulting coefficients can be evaluated efficiently enough to be used directly in ray-tracing calculations, and we provide an example calculation. Second, we use stochastic averaging to assess the effect of subgrid turbulence on the volume-averaged emissivity and along the way provide a prescription for a turbulent emissivity. We find that for parameters appropriate to EHT sources turbulence reduces the emissivity slightly. In the infrared, turbulence can dramatically increase the emissivity.

Astronomy & Astrophysics↗

A Scalable Gaussian Process Approach to Shear Mapping with MuyGPs

Analysis of cosmic shear is an integral part of understanding structure growth across cosmic time, which in turn provides us with information about the nature of dark energy. Conventional methods generate shear maps from which we can infer the matter distribution in the universe. Current methods (e.g., Kaiser–Squires inversion) for generating these maps, however, are tricky to implement and can introduce bias. Recent alternatives construct a spatial process prior for the lensing potential, which allows for inference of the convergence and shear parameters given lensing shear measurements. Realizing these spatial processes, however, scales cubically in the number of observations—an unacceptable expense as near-term surveys expect billions of correlated measurements. Therefore, we present a linearly scaling shear map construction alternative using a scalable Gaussian process prior called MuyGPs. MuyGPs avoids cubic scaling by conditioning interpolation on only nearest neighbors and fits hyperparameters using batched leave-one-out cross-validation. This work is the first step toward a full, scalable mass mapping method. We work in a simplified regime where we validate our method by interpolating and analyzing maps given noisy point-estimate data from all three shear fields, taken from a suite of N -body ray-tracing simulations. We also show that we can perform these operations at the scale of billions of galaxies on high-performance computing platforms.

79 ASTRONOMY AND ASTROPHYSICS↗

Changing-look Active Galactic Nuclei Powered by Disk Tearing

Changing-look active galactic nuclei (CLAGN) feature order-of-magnitude variability in both the continuum and broad line luminosities on months-to-years-long timescales, and are currently unexplained. Simulations have demonstrated that rotating black holes (BHs) sometimes tear apart tilted accretion disks. These tearing events violently restructure the disk on timescales much shorter than a viscous timescale, hinting at a connection to CLAGN. Here, we show that disk tearing can power changing-look events. We report synthetic observations of an extremely high-resolution, three-dimensional, general-relativistic magnetohydrodynamic simulation of a geometrically thin, tilted accretion disk around a rapidly rotating, 10 8 M ⊙ BH. We perform ray-tracing calculations that follow the thermal disk emission to both a line-of-sight camera and to a distribution of cameras in a prescribed torus-like broad-line region (BLR). The continuum photoionizes the BLR, and we calculate the resulting spectrum. Both the continuum and line luminosities undergo order-of-magnitude swings on months-to-years-long timescales. We find shorter, week-long variability driven by the geometric precession of the inner disk and an intraday quasiperiodic oscillation driven by radial breathing of the inner disk. When the torn disk precesses, it causes asymmetric illumination of the BLR, driving time-evolving red-to-blue asymmetries of the broad emission lines that may be a smoking gun for disk tearing. We also make predictions for future photometric observations from ULTRASAT and Vera Rubin Observatory, given a Shakura–Sunyaev–like effective temperature profile, both of which may play an important role in detecting future changing-look events. Future work should consider nonthermal emission, especially from the time-variable corona.

accretion↗

Development of the Generation 3 Particle Pilot Plant Falling Particle Receiver

A falling particle receiver (FPR) has been designed to integrate with the G3P3-USA pilot plant currently being constructed at the National Solar Thermal Test Facility (NSTTF) at Sandia National Laboratories. This receiver integrates several innovative design features including a converging tunnel (SNOUT), an optimized cavity geometry, and a multistage “catch-and-release” trough. Details about the integration of these features and the final G3P3-USA FPR design and construction are described. Ray-tracing models of the FPR utilizing the NSTTF heliostat field are developed leveraging previous modeling efforts using SolTrace. Models demonstrate that at least 1.5 MWth of incident radiative energy can be provided to the FPR on a clear day throughout a typical year in Albuquerque, NM. Spillage fluxes around the periphery of the aperture are within acceptable bounds for the majority of the year. Intercept factors are computed for each utilized heliostat at the Vernal equinox to provide guidance to heliostat operators during operation of the system.

Mills, Brantley↗

Frequency-Selectable Laser Source for Cosmic Microwave Background Experiments

Cosmic Microwave Background (CMB) experiments measure remnant radiation from the early universe and use that data to determine fundamental properties of the universe. We can constrain key parameters, such as \textit{r}, the cosmic tensor-to-scalar ratio, and $N_\text{eff}$, the effective number of relativistic species, by analyzing the CMB power spectra. Improving our measurements of the CMB requires improving our instrument systematics, one of the most important of which is detector bandpass. Current experiments use a Fourier Transform Spectrometer (FTS) to measure bandpass. However, the FTS is systematics limited, and cannot achieve the accuracy needed to make improved CMB measurements. For this reason, we are developing a new instrument, the Frequency-Selectable Laser Source (FLS) to decrease the uncertainty in bandpass by an order of magnitude. In this paper, we describe work completed to support the version 2 upgrade to the FLS. Using ray-tracing software, we modeled the FLS optics to set physical tolerances for the new design. We also discuss the laser calibration, and future work to be completed in further development of the FLS upgrade.

Rosen-Turits, Gabriel M. [Fermilab]↗

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↗

Verifying Ray Tracing Amplitude Methods for Global Magnetospheric Modeling

Abstract Ray tracing is a commonly used method for modeling the propagation of electromagnetic waves in Earth's magnetosphere. To apply ray tracing results to global models of wave‐particle interaction such as energetic electron scattering, it is useful to map the discrete rays to a volume filling mesh. However, some methods have inherent losses of energy from the wave source, or do not account for the full range of wave properties within a sample volume. We have developed and tested a 3D magnetospheric ray tracing code “MESHRAY” which resolves these issues. MESHRAY uses the conservation of Poynting flux through ray triplets with finite volume to determine the local field amplitudes. Electromagnetic wave energy density from all ray data points is mapped to a mesh and verified against the wave source power for energy conservation varying time step length, number of rays, and total time steps. We find that the method is self‐consistent and numerically robust. We further investigate whether the neglect of phase information and superposition has a significant impact on the accuracy of mapping wave intensity to a mesh. We find excellent agreement between the analytic solution for waves emitted by a line source in a plane‐stratified medium and an equivalent ray tracing solution. When phase information is excluded, ray tracing reproduces an average amplitude spread over regions of coherent constructive and destructive interference. This may be an important consideration for interpolating ray tracing results of longer wavelength waves such as magnetosonic, electromagnetic ion cyclotron, or ULF waves.

79 ASTRONOMY AND ASTROPHYSICS↗