
Neutron wall loading in stellarator optimization
Enrique Miralles-Dolz, R. Michael Churchill, Dario Giovanni Panici, Jacob Schwartz, Tim Bohm, Connor Moreno, Paul P H Wilson
Stellarator neutronics is inherently three-dimensional, since non-axisymmetric equilibria produce spatially varying neutron wall loading (NWL). Localized NWL peaks affect first-wall lifetime, shielding requirements, magnet protection, maintenance planning, and ultimately plant availability, but resolving these peaks with high-fidelity Monte Carlo transport remains too expensive for direct use in stellarator optimization. This work integrates a deterministic line-of-sight NWL model into the DESC stellarator optimization package so that NWL can be evaluated as a fast, differentiable engineering metric during stage-I equilibrium optimization. The model treats the deuterium-tritium plasma as a volumetric neutron source, evaluates the Bosch-Hale fusion reactivity from prescribed flux-surface profiles, and maps source contributions to first-wall points through a geometry-dependent kernel. The resulting local NWL map can be used directly in optimization objectives or constraints. The DESC implementation is validated against ParaStell/OpenMC Monte Carlo calculations for circular axisymmetric, ITER-like axisymmetric, and non-axisymmetric configurations; the surface-averaged NWL agrees within about 1.5% across these cases and the deterministic maps reproduce the dominant spatial loading structures while tending to conservatively overestimate the Monte Carlo wall loading. In proof-of-concept optimizations, current-profile variation at fixed plasma boundary reduced peak NWL and the peaking factor by 13.84% and 13.29%, respectively, while quasi-symmetry optimization through plasma boundary shaping reduced peak and mean NWL by 17.22% and 12.97%. These results demonstrate that differentiable NWL objectives and constraints can reveal useful engineering design directions unavailable to conventional plasma physics objectives.


