Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “field validation”

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 883 records · Page 49

Ultra-low field 13 C MRI of hyperpolarized pyruvate

Medicine is evolving beyond therapy largely predicated on anatomical information and towards incorporating patient-specific molecular biomarkers of disease for more accurate diagnosis and effective treatment. The complementary combination of hyperpolarization by spin-lock induced crossing signal amplification by reversible exchange (SLIC SABRE) and low field magnetic resonance imaging (MRI) can enable accessible metabolic imaging to advance personalized medicine. Hyperpolarized 13 C-enriched pyruvate has demonstrated promise for imaging metabolism in cancer, heart disease and neurodegenerative disorders; however, broader clinical adoption awaits validated clinical indications, and is further constrained by the cost and limited availability of current hyperpolarization technology. Parahydrogen-based polarization techniques, paired with low-cost high-performance MRI at millitesla fields, offer a means of broadening the reach of metabolic imaging. Here we show results demonstrating in situ hyperpolarization of pyruvate at 6.5 mT by SLIC SABRE, followed by immediate readout without field cycling or sample shuttling. We achieve 13 C signal enhancements several million times above thermal equilibrium at 6.5 mT, corresponding to polarization levels of approximately 3%. Leveraging this enhancement, we perform 13 C MRI and acquire NMR spectra with resolution sufficient to distinguish chemical shifts between pyruvate isotopomers. These results show a viable pathway towards accessible metabolic imaging with hyperpolarized 13 C MRI at ultra-low field.

Medical and clinical diagnostics↗

Science Validation for Dark Energy Research with Optical Imaging Surveys

The Universe has been expanding at an accelerating rate over the past several billion years, as though an unknown form of dark energy permeates all of space. The ultimate scientific goal of the proposed research is to distinguish between different physical mechanisms that could account for this observed accelerated expansion, for example, the zero-point energy of the vacuum, a dynamical form of energy that varies in time and/or space, or a modification to our theory of gravity. Wide- field optical imaging surveys of the night sky can test these competing models by measuring both the cosmic expansion history and the growth of large-scale structure. To this end, the Dark Energy Survey (DES) has cataloged several hundred million galaxies and thousands of supernovae. The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), will enlarge the census to billions of galaxies and hundreds of thousands of supernovae. A critical question for these dark energy experiments is whether systematic uncertainties can continue to be controlled at a level to keep pace with the statistical precision offered by such enormous datasets. The immediate research objectives of this project were (1) to prepare and validate input datasets that are the foundation of cosmological analyses with DES, and (2) to prepare for value-added characterization of Rubin Observatory commissioning data to inform early operations and accelerate the realization of dark energy science from LSST data products. For DES, we assembled and curated cosmology-ready data releases that include value-added components such as enhanced photometric and astrometric calibrations, alternative source extraction algorithms, maps of the survey coverage and survey conditions, object classifications, object quality selections, galaxy shapes, and photometric redshifts. We used the galaxy clustering technique to validate the photometric redshift distributions of various galaxy samples to be used as lenses in combined studies of galaxy clustering and weak gravitational lensing. The galaxy clustering redshift analysis was enhanced by use of a larger sample of reference galaxies from the eBOSS spectroscopic survey that extends to higher redshifts. For LSST, we prepared for science validation studies of commissioning data aimed at dark energy science capability that extend beyond the normative system-level tests to be done by the Rubin Observatory Construction Project. We identified a set of proposed survey strategies and candidate target fields that could be observed during the commissioning period to enhance science validation activities related to studies of dark energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comparative Analysis of HEATNETS for Geothermal Network Performance

Thermal energy networks (TENs), also known as 5th generation district energy systems, or more specifically geothermal networks when exchanging heat with geothermal boreholes, are an important technology for decarbonization. In these networks an ambient loop connects buildings and thermal sources, such as a borehole field, to exchange energy and maintain a desired loop temperature. Water-source heat pumps are used at the buildings to connect to the ambient or thermal loop to meet to the building heating and cooling loads and maintain comfort. A semi-transient, reduced-order technical model and techno-economic model, called HEATNETS, has been developed at NREL that captures the flow of energy around a TEN. In this work, a comparison of the HEATNETS technical model and a well-known coding platform used for modeling geothermal networks, TRNSYS, has been completed for a proposed geothermal network as a verification and validation process. Hourly data provided from the TRNSYS simulation included building loads, pumping power, heat pump power, temperature entering and leaving the borehole field, and mass flow rates. The hourly borehole temperatures were used to create a linear regression model utilized in HEATNETS to estimate the borehole field heat exchange. The building loads and mass flow rates were direct inputs to HEATNETS while the pumping power, heat pump power, borehole temperatures, and coefficients of performance were all simulated and calculated by HEATNETS, allowing for direct comparison of the thermal energy transfer HEATNETS considers the full process from design inputs to economic outputs and can provide modeling options for high-level initial system design and operational optimization. This study focuses on a validation of HEATNETS using results from TRNSYS. HEATNETS is not intended to replace other modeling tools, but this work demonstrates, via a comparison with an industry standard code, that HEATNETS can be a unique, high-level and rapid modeling tool for estimating the performance of a full geothermal network system.

15 GEOTHERMAL ENERGY↗

Development and Validation of Southeast Asia Solar Resource Data [Slides]

Lack of access to high-quality, publicly available, time series solar data to inform decisions that will transform energy sectors in Southeast Asia is a challenge. The solution is to level the playing field by offering free, high-quality, robust solar data to inform private sector investment and policymaking. This is done by (1) leveraging deep NREL expertise in atmospheric science, solar resource assessment, high-performance computing, and cloud-based data dissemination, (2) producing and validating high spatial and temporal resolution solar resource data, (3) making data available on the USAID-funded global Renewable Energy Data Explorer platform, (4) providing capacity building for data and applications, and (5) informing future demand-driven tool development.

14 SOLAR ENERGY↗

Improved Treatment of 1–4 Interactions in Force Fields for Molecular Dynamics Simulations

Traditional force fields commonly use a combination of bonded torsional terms and empirically scaled nonbonded interactions to capture 1-4 energies and forces of atoms separated by three bonds in a molecule. While this approach can yield accurate torsional energy barriers, it often leads to inaccurate forces and erroneous geometries and creates an interdependence between dihedral terms and nonbonded interactions, complicating parametrization and reducing transferability. Here, in this paper, we demonstrate that 1-4 interactions can be accurately modeled using only bonded coupling terms, eliminating the need for arbitrarily scaled nonbonded interactions altogether. Furthermore, by leveraging the automated parametrization capabilities of the Q-Force toolkit, we efficiently determine the necessary coupling terms without the need for manual adjustment. Our approach is first validated on a range of small molecule systems, encompassing both flexible and rigid structures, and shows a significant improvement in force field accuracy, obtaining subkcal/mol mean absolute error for every molecule tested. We further extend the bonded-only model for 1-4 interactions to Amber ff14sb, CHARMM36, and OPLS-AA force fields to reproduce ab initio gas and implicit solvent ϕ,ψ surfaces of alanine dipeptide.

Abdullah, Aalim S. [University of California, Berk↗

Technical and Economic Evaluation of the First Ever Polymer Flood Field Pilot to Enhance the Recovery of Heavy Oils on Alaska's North Slope via Machine Assisted History Matching

Polymer flooding has become globally established as a potential enhanced oil recovery method for heavy oils. To determine whether this technology may be useful in developing the substantial heavy oil resources on the Alaska North Slope, a polymer flood field pilot commenced at the Milne Point Unit in August 2018. This study seeks to evaluate the results of the field pilot on a technical and economic basis. A reservoir simulation model is constructed and calibrated to predict the oil recovery performance of the pilot through machine-assisted reservoir simulation techniques. To replicate the early water breakthrough observed during waterflooding, transmissibility contrasts are introduced into the simulation model, forcing viscous fingering effects. In the ensuing polymer flood, these transmissibility contrasts are reduced to replicate the restoration of injection conformance during polymer flooding. Transmissibility contrasts are later reinstated to replicate fracture overextension interpreted in one of the producing wells. The calibrated simulation models produced at each stage of the history matching process are used to forecast oil recovery. These forecasts are used as input for economic analysis, incremental to waterflooding expectations. The simulation forecasts indicate that polymer flooding significantly increases the heavy oil production for this field pilot compared to waterflooding alone, yielding attractive project economics. However, meaningful variations between simulation scenarios demonstrate that a simulation model is only valid for prediction if flow behavior in the reservoir remains consistent with that observed during the history matched period. Critically, this means that a simulation model calibrated for waterflooding may not fully capture the technical and economic benefits of an enhanced oil recovery process such as polymer flooding. Subsequently, the simulation model and economic model are used in conjunction to conduct a sensitivity analysis for polymer flood design parameters, from which recommendations are provided for both the continued operation of the current field pilot and future polymer flood designs. The results demonstrate that a higher polymer concentration can be injected due to the development of fractures in the reservoir. The throughput rate should remain high without exceeding operating constraints. A calculated point-forward polymer utilization parameter demonstrates the decreasing efficiency of the polymer flood at later times in the pattern life. Future projects will benefit from starting polymer injection earlier in the pattern life. A pattern with tighter horizontal well spacing will observe a greater incremental benefit from polymer flooding.

Keith, Cody↗

Spatial coherence of room-temperature monolayer WSe 2 exciton-polaritons in a trap

The emergence of spatial and temporal coherence of light emitted from solid-state systems is a fundamental phenomenon intrinsically aligned with the control of light-matter coupling. It is canonical for laser oscillation, emerges in the superradiance of collective emitters, and has been investigated in bosonic condensates of thermalized light, as well as exciton-polaritons. Our room temperature experiments show the strong light-matter coupling between microcavity photons and excitons in atomically thin WSe 2 . We evidence the density-dependent expansion of spatial and temporal coherence of the emitted light from the spatially confined system ground-state, which is accompanied by a threshold-like response of the emitted light intensity. Additionally, valley-physics is manifested in the presence of an external magnetic field, which allows us to manipulate K and K’ polaritons via the valley-Zeeman-effect. Our findings validate the potential of atomically thin crystals as versatile components of coherent light-sources, and in valleytronic applications at room temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnetic ground states of honeycomb lattice Wigner crystals

Lattice Wigner crystal states stabilized by long-range Coulomb interactions have recently been realized in two-dimensional moiré materials. We employ large-scale unrestricted Hartree-Fock techniques to unveil the magnetic phase diagrams of honeycomb lattice Wigner crystals. For the three lattice filling factors with the largest charge gaps, $n=2/3,1/2,1/3$, the magnetic phase diagrams contain multiple phases, including ones with non-collinear and non-coplanar spin arrangements. We discuss magnetization evolution with external magnetic field, which has potential as an experimental signature of exotic spin states. Our theoretical results could potentially be validated in moiré materials formed from group VI transition metal dichalcogenide twisted homobilayers.

36 MATERIALS SCIENCE↗

4D beam matrix reconstruction in particle accelerators

Transverse beam parameters in particle accelerators are commonly described using the Twiss parameters, which are experimentally accessible yet inherently limited because they neglect correlations between different transverse coordinates. Such correlations frequently arise from uncompensated cathode magnetic fields or misaligned focusing quadrupoles, affecting beam quality and accelerator performance. To address this limitation, we propose and validate a novel diagnostic method for the complete four-dimensional (4D) transverse beam matrix. Our method involves passing the beam through a beamline comprising both conventional and skew quadrupole magnets, followed by downstream measurements of the resulting two-dimensional (2D) beam profiles. These measurements represent distinct 2D projections of the underlying 4D transverse phase–space distribution. By systematically varying quadrupole strengths, multiple independent projections of the beam phase space are obtained. We reconstruct the original 4D beam matrix from these measured projections using an optimization-based least-square fit, providing fast and robust reconstruction regardless of the specific beamline configuration. Through extensive numerical simulations and realistic particle-tracking studies, we demonstrate the diagnostic’s accuracy, robustness, and capability to achieve reconstruction uncertainties smaller than measurement errors, particularly when employing sufficient numbers of quadrupole scans. This method presents a powerful and flexible approach for comprehensive beam characterization and accelerator tuning.

43 PARTICLE ACCELERATORS↗

How Generalizable is a Machine-Learning Approach for Modeling Hub-Height Turbulence Intensity?

Hub-height turbulence intensity is essential for a variety of wind energy applications. However, simulating it is a challenging task. Simple analytical models have been proposed in the literature, but they all come with significant limitations. Even state-of-the-art numerical weather prediction models, such as the Weather Research and Forecasting model, currently struggle to predict hub-height turbulence intensity. Here, we propose a machine-learning-based approach to predict hub-height turbulence intensity from other hub-height and ground-level atmospheric measurements, using observations from the Perdigao field campaign and the Southern Great Plains atmospheric observatory. We consider a random forest regression model, which we validate first at the site used for training and then under a more robust round-robin approach, and compare its performance to a multivariate linear regression. The random forest successfully outperforms the linear regression in modeling hub-height turbulence intensity, with a normalized root-mean-square error as low as 0.014 when using 30-minute average data. In order to achieve such low root-mean-square error values, the knowledge of hub-height turbulence kinetic energy (which can instead be modeled in the Weather Research and Forecasting model) is needed. Interestingly, we find that the performance of the random forest generalizes well when considering a round-robin validation (i.e., when the algorithm is trained at one site such as Perdigao or Southern Great Plains) and then applied to model hub-height turbulence intensity at the other location.

17 WIND ENERGY↗

Optimal Network Reconfiguration and Scheduling With Hardware-in-the-Loop Validation for Improved Microgrid Resilience

With the increased occurrence of various major extreme weather events, power outages and prompt power system restorations have recently drawn more attention to the resilience and recovery of power systems. From the perspective of a more resilient power delivery at the distribution grid, system restoration using network topology reconfiguration together with optimal scheduling of distributed energy resources are adopted in this paper. The proposed optimization model aims at minimizing the total load shedding cost and other operational costs, in which linearized topological constraints borrowed from graph theory and linearized DistFlow models are respectively used to maintain the radial network topology and power flow balance after system contingencies. To demonstrate the applicability of the proposed strategy, a real-world case study of a networked three-microgrid system in Adjuntas, Puerto Rico, is used with the consideration of different independent/interconnected microgrid scenarios, contingencies, and fairness settings. Furthermore, hardware-in-the-loop testing is conducted for the same three-microgrid network, where the closely matched results with the simulated ones have validated the effectiveness of the proposed restoration strategy, which is now ready to move one step forward towards field deployment. Finally, to test the proposed restoration strategy in a larger networked system, the modified IEEE-33 bus test distribution system is considered, and the results show a more resilient power delivery for critical loads under three and four line outages.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Modular Mechanism for Downhole Weight-on-Bit and Torque Reaction in Small Diameter Boreholes

Drilling systems that use downhole rotation must react torque either through the drill-string or near the motor to achieve effective drilling performance. Problems with drill-string loading such as buckling, friction, and twist become more severe as hole diameter decreases. Therefore, for small holes, reacting torque downhole without interfering with the application of weight-on-bit, is preferred. In this paper, we present a novel mechanism that enables effective and controllable downhole weight on bit transmission and torque reaction. This scalable design achieves its unique performance through four key features: (1) mechanical advantage based on geometry, (2) direction dependent behavior using rolling and sliding contact, (3) modular scalability by combining modules in series, and (4) torque reaction and weight on bit that are proportional to applied axial force. As a result, simple mechanical devices can be used to react large torques while allowing controlled force to be transmitted to the drill bit. We outline our design, provide theoretical predictions of performance, and validate the results using full-scale testing. The experimental results include laboratory studies as well as limited field testing using a percussive hammer. These results demonstrate effective torque reaction, axial force transmission, favorable scaling with multiple modules, and predictable performance that is proportional to applied force.

15 GEOTHERMAL ENERGY↗

Radar Wind Profiler (RWP) and Radio Acoustic Sounding System (RASS) Instrument Handbook

The radar wind profiler (RWP) is an active remote-sensing instrument that can routinely, and virtually unattended, observe wind and turbulence in the troposphere through scattering from clear-air irregularities of the atmospheric refractive index (Gage and Balsley 1978). The principle of Doppler radars in general is based on sending electromagnetic pulses in the vertical and several tilted directions (vertical and oblique beams) and measuring the signal that is scattered back by atmospheric turbulence at all heights and received at the antenna. The RWPs provides measurements of backscattered signal strength and wind profiles nominally between 0.1 km and 6 km. The RWP operation assumes that the wind field is homogeneous over the spatial separation of the antenna beams, which is a safe assumption under stable atmospheric conditions. However, the wind may be different between the beams due to significant changes in the wind field, resulting in erroneous wind calculations. Therefore, the radial measurements are averaged over a sufficient time period to validate the assumption of homogeneity.

54 ENVIRONMENTAL SCIENCES↗

Assess and Benchmark Magneto-Inertial Fusion (MIF) Scaling. Final Technical Report for the SNL/Rochester ALPHA Follow-on Project

This project was a follow-on to the Sandia National Laboratories (SNL) and the Laboratory for Laser Energetics (LLE) ARPA-E ALPHA project entitled “Demonstrating Fuel Magnetization and Laser Heating Tools for Low-Cost Fusion Energy”. The primary purpose of this follow-on project was to obtain additional data at the OMEGA facility to help better understand how MagLIF, a platform that has already demonstrated the scientific viability of magneto-inertial fusion, scales across a factor of 1000 in driver energy. A secondary aspect of this project was to extend simulations and analysis at SNL to cover a wider magneto-inertial fusion (MIF) parameter space and test scaling of those models across this wide range of input energies and conditions of the target. This work was successful in improving understanding of how key physics elements of MIF scales and improves confidence in setting requirements for fusion gain with larger drivers. The OMEGA experiments at the smaller scale verified the hypothesis that preheating the fuel plays a significant role in introducing wall contaminants that mix into the fuel and significantly degrade fusion performance. This contamination not only impacts target performance but the optimal input conditions for the target. However, analysis at the Z-scale showed that target performance at high preheat levels is limited by the Nernst effect, which advects magnetic flux from the hot spot, reducing magnetic insulation and consequently reduces the temperature of the fuel. The combination of MagLIF experiments at the disparate scales of OMEGA and Z along with a multiscale 3D simulation analysis has led to new insight into the physical mechanisms responsible for limiting target performance and provides important benchmarks to assess target scaling more generally for MIF schemes. Finally, in addition to the MagLIF related work, a semi-analytic model of liner driven Field Reversed Configuration (FRC) was developed that predicts the fusion gain for such systems. This model was also validated with 2D radiation magneto-hydrodynamic simulations and predicts that fusion gains of near unity could be driven by the Z machine.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Innovative Subsurface Learning and Hawaiian Exploration using Advanced Tomography (ISLAND HEAT) Phase 1 (Final Report)

The Innovative Subsurface Learning and Hawaiian Exploration using Advanced Tomography (ISLAND HEAT) project aims to leverage and enhance existing Hawai'i geothermal play fairway analysis (PFA) results and validate a conceptual model-driven, optimized, least-cost exploration and geophysical suite at a proven geothermal field with known permeability, the Puna Geothermal Venture, located in the Lower East Rift Zone. The abundance of geophysical scrutinization before, during, and after the 2018 Kilauea eruption at and around the broader Puna system offers a rare opportunity to elucidate geophysical expressions within a dynamic magmatic rift setting. Development of the methodology continues with its application at a second prospective site identified by the PFA, Mauna Kea, which serves as confirmation of the integrated approach.

15 GEOTHERMAL ENERGY↗

HFTS-1 (Final Report)

Despite the long history of hydraulic fracturing, the optimal number of fracturing stages during multi-stage fracture stimulation in horizontal wells can be a complex subject. In addition to the increased expense of multistage fracturing in horizontal wells, increasing the number of fracturing stages does not always correlate with a commensurate increase in production. The problem is the application of a uniform fracture stimulation design to all stages with no consideration for geological variations along the wellbore. The result is an inefficient use and costly waste of energy and water. Optimization of the fracturing process requires an understanding of the cause-and-effect relationship between fracturing parameters and local geological properties at a given location along the wellbore. Realizing that the generalized rock mechanics theories and hypotheses are not truly applicable to fractured and laminated shales, quantifiable impacts of a shale’s geomechanical and depositional features are a prerequisite for design and implementation of optimized hydraulic fractures. The overarching goal of this project is to understand and define the relationships of shale geology and fracture dynamics using detailed field data that includes coring of the fracture domain. Finally, resource recovery from shale formations is estimated to be less than fifteen percent. Research proposed in this project will establish the foundations for investigating enhanced recovery techniques for increased resource recovery in existing fracture treated wells. Natural gas as an EOR fluid is in the initial stages of being used broadly in the oil and gas industry and as such, many aspects of the process need to be researched and addressed prior to widespread acceptance. In addition, re-fracturing techniques, particularly using cemented liners also need to be further investigated. This project explores several of these EOR techniques in greater detail, not just at the theoretical level, but through actual field test sites. This report delves into the research program, identified tasks, several studies, observations, results, and conclusions under the Hydraulic Fracturing Test Site (HFTS) project. This test project is subdivided into three distinct phases and each phase is discussed in detail under separate sections. The original test site for Phase 1 as well as the Phase 2 site are in Midland Basin in Texas and operated by Laredo Petroleum. In the first two phases, the program targeted the Wolfcamp formation in Midland Basin. Phase 3 of the study targeted the EagleFord play in South Texas with Devon Energy as the operator. Significant conclusions from this study include: 1) lateral and vertical fracture growth mapping as well as proppant transport behavior, 2) Use of field gas and its applicability for huff n puff EOR injection in the Midland Basin, 3) Development and validation of proppant log, both at core scale as well as in-fill well implementation by analyzing drilling mud return samples, 4) Applicability of liner refracturing including best practices for enhanced unit level productivity. This report will detail all these studies as well as technical results and conclusions from said studies. These include general subsurface characterization work, completion designs as implemented, studies looking at fracture geometry, high resolution microseismic study, advanced crosswell strain mapping studies including drainage characterization, core characterization results, proppant analysis (core and observation wells), geochemistry, etc. to name a few.

02 PETROLEUM↗

Particle-in-cell Simulations of Mildly Relativistic Outflows in Kilonova Emissions

Abstract The electromagnetic emission from neutron star mergers is comprised of multiple components. Synchrotron emission from the disk-powered jet and thermal emission from the merger ejecta (powered by a variety of sources) are among the most studied sources. The low masses and high velocities of the merger ejecta quickly develop conditions where emission from collisionless shocks becomes critical and synchrotron emission from the merger ejecta constitutes a third component to the observed signal. The aim of this project is to examine shock development, magnetic field generation, and particle acceleration in the case of mildly relativistic shocks, which are expected when the tidal ejecta of neutron star mergers drive a shock into the external medium. Using LANL’s vector particle-in-cell (VPIC) code, we have run simulations of such mildly relativistic, collisionless, weakly magnetized plasmas and computed the resultant magnetic fields and particle energy spectra. We show the effects of varying plasma conditions, as well as explore the validity of using different proton-to-electron mass ratios in VPIC. Our results have implications for observing late-time electromagnetic counterparts to gravitational wave detections of neutron star mergers.

79 ASTRONOMY AND 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↗