Accelerating split-exponential track length estimator on GPU for Monte Carlo simulations

In the context of computing 3D volumetric tallies for nuclear applications, the combination of Monte Carlo methods and high-performance computing is essential to achieve accurate yet computationally feasible simulations that meet industrial time constraints. The next-event Split Exponential Track-Length Estimator ( se TLE) is particularly well suited for estimating tallies on meshes. To alleviate the computational burden associated with se TLE, such as sampling numerous outgoing pseudoparticles at each collision, estimating cross sections, performing ray tracing through complex geometries, and accumulating scores across the geometry, we leverage the parallel computing capabilities of Graphics Processing Units (GPUs). We assess the performance of our implementation using two shielding configurations and one criticality benchmark. Both photon and neutron transport simulations are considered. Scores are evaluated over Cartesian meshes, material volumes, and energy group structures. In all cases, acceleration factors greater than unity are observed in the detectors, reaching several hundred in selected regions of the phase space. In a final experiment, we demonstrate that our GPU-based implementation achieves a net energy gain (in Watt) even when compared to a conventional CPU-based TLE, despite the additional computational cost of GPU use.