ABSTRACT Effective management of particle and power exhaust is essential for fusion reactor operation, requiring accurate 3D modeling of divertor physics. The EMC3 code is a well‐established Monte Carlo tool for plasma transport in stellarator edge regions, but until now lacked and diamagnetic drifts. This work presents and verifies a first implementation of these drifts in EMC3 using the Method of Manufactured Solutions (MMS) in a cylindrical test geometry. Building on earlier developments of gradient, divergence, and tracing schemes, the present work introduces drift terms employing analytical or temperature‐based proxies for the potential field. The simulations reproduce the expected drift‐induced transport, including perpendicular fluxes and poloidal asymmetries, and show the theoretical convergence with number of particles . Using the temperature as potential proxy does not lead to accumulating statistical noise, confirming numerical stability. This implementation marks a key step toward a fully self‐consistent treatment of plasma drifts in EMC3, paving the way for inclusion of potential field computation and boundary conditions.
Wolf et al. (2026) studied this question.