First Principles Classical Calculations of WHAM2 and BEAM Mirror Neutron Rates using the GENRAY/ CQL3D-M codes

The CQL3D-M neutral beam/RF/nonlinear-collisional bounce-average Fokker-Planck and GENRAY ray tracing codes provide a first-principles classical calculation of e/D+ distribution functions in axisymmetric mirror geometry. For WHAM2 (under construction at UW, Madison) and BEAM (Q=1 design phase), we obtain DD neutron rates in a manner similar to [1]. A self-consistent parallel electric field keeps electron and ion densities approximately equal. The collision operator conserves particles and energy. In WHAM2 [2] with 100 kW NB power at midplane 45 deg, we obtain neutron rate 1.3e12 n/s. Adding 50 kW of RF fast wave power gives 1.4e13 n/s, a 25X increase. Electron density is 2.1e13/cc; Te is 3.3 keV. In BEAM [3], under design, we obtain equivalent Q=1 with 10MW/100 keV NB power.

References:

[1] R. W. Harvey, Y. V. Petrov, and C. B. Forest, AIP Conf. Proc.1771, 040002 (2016). https://doi.org/10.1063/1.4964187

[2] D. Endrizzi, J.K. Anderson, M. Brown, et al., [WHAM] Accepted by J. Plasma Phys. (2023).

[3] C.B. Forest, J.K. Anderson, et al., et al., [BEAM] Under consideration by J. Plasma Phys. (2023).