Abstract:Objective.This work aims to develop a graphics processing unit (GPU)-accelerated Monte Carlo code for the coupled transport of photon, electron/positron and neutron over a broad range of energies for medical applications.Approach.By separating the MC evolution of radiation into source, transport, and interaction kernels, the branch divergence was alleviated. The memory coalescence was achieved by vectorizing the access pattern in which the secondary particles were archived. To accelerate further particle tracking, ray-tracing hardware acceleration in the Nvidia OptiXTMframework was applied. For photon and electron/positron, the EGSnrc interaction modules were ported as a GPU-optimized configuration. For neutron, a group-wised transport based on NJOY21 preprocessed data was implemented. The developed code was validated against CPU-based FLUKA. Neutron, x-ray and electron beams incident on water and ICRP phantoms were simulated. The neutron energy group and the transport parameters of photon and electron were set to be the same in both codes. A single Nvidia RTX 4090 card was used in this code while all 20 threads of a single Intel Core i9-10900K node were used in FLUKA.Main results.The number of histories was set to ensure that statistical uncertainties lower than 2% for all voxels whose doses were larger than 20% of the maximum. In all cases, the dose differences in the voxels between the codes were within 2.5%. For photons and electrons, the developed code was 150-300 times faster than FLUKA in both geometries. For neutrons, the code was respectively 80 and 135 times faster in the water and ICRP phantoms than FLUKA.Significance.This study offers an appropriate solution for uncoalesced memory access and branch divergence commonly encountered in coupled MC transport on the GPU architecture. The formidable acceleration in computing times and accuracy shown in this study can promise a routine clinical use of MC simulations.

Parallel Computing for Doppler Broadening of Neutron Cross Sections

Research on Fast-Doppler-Broadening of Neutron Cross Sections

Research on On-the-Fly Doppler Broadening Method for Neutron Cross Section

Parallelization of Kinetic Theory Simulations

Three-dimensional Variational Nodal Method Parallelization for Pin Resolved Neutron Transport Calculations

Parallel Computation Research of Neutron Time-Spatial Kinetics Program for Fast Reactor

Parallelization of Mcnp Monte Carlo Neutron and Photon Transport Code in Parallel Virtual Machine and Message Passing Interface

Hybrid Windowed Networks for On-the-fly Doppler Broadening in RMC Code

Research on GPU-based Acceleration Method for Monte Carlo Neutron Geometry Treatment

A comparative study of history-based versus vectorized Monte Carlo methods in the GPU/CUDA environment for a simple neutron eigenvalue problem

ICONE23-1418 DEVELOPMENT AND VALIDATION OF THE FAST DOPPLER BROADENING MODULE COUPLED WITHIN RMC CODE

Research on GPU Acceleration for Monte Carlo Criticality Calculation

A GPU-accelerated Monte Carlo code, RT2for coupled transport of photon, electron/positron, and neutron

Development of A New Parallel S-N Code for Neutron-Photon Transport Calculation in 3-D Cylindrical Geometry

Study on Parallel Computing of Radiation Detector Simulation Based on Geant4

EVALUATION OF SPEEDUP OF MONTE CARLO CALCULATIONS OF TWO SIMPLE REACTOR PHYSICS PROBLEMS CODED FOR THE GPU/CUDA ENVIRONMENT

Optimization Treatment of Point-Wise Nuclear Data in Monte Carlo Criticality and Burnup Calculations

Data Parallelization Algorithms for the Direct Simulation Monte Carlo Method for Rarefied Gas Flows on the Basis of OpenMP Technology

Prompt: Probability-Conserved Cross Section Biasing Monte Carlo Particle Transport System

Accurate Simulation of Neutrons in Less Than One Minute Pt. 2: Sandman—GPU-Accelerated Adjoint Monte-Carlo Sampled Acceptance Diagrams

A scalable multilevel domain decomposition preconditioner with a subspace-based coarsening algorithm for the neutron transport calculations