Efficient calculation of the self magnetic field, self-force, and self-inductance for electromagnetic coils

Siena Hurwitz,Matt Landreman,Thomas M. Antonsen

2023-10-13

Abstract:The design of electromagnetic coils may require evaluation of several quantities that are challenging to compute numerically. These quantities include Lorentz forces, which may be a limiting factor due to stresses; the internal magnetic field, which is relevant for determining stress as well as a superconducting coil's proximity to its quench limit; and the inductance, which determines stored magnetic energy and dynamics. When computing the effect on one coil due to the current in another, these quantities can often be approximated quickly by treating the coils as infinitesimally thin. When computing the effect on a coil due to its own current (e.g., self-force or self-inductance), evaluation is difficult due to the presence of a singularity; coils cannot be treated as infinitesimally thin as each quantity diverges at zero conductor width. Here, we present novel and well-behaved methods for evaluating these quantities using non-singular integral formulae of reduced dimensions. These formulae are determined rigorously by dividing the domain of integration of the magnetic vector potential into two regions, exploiting appropriate approximations in each region, and expanding in high aspect ratio. Our formulae show good agreement to full finite-thickness calculations even at low aspect ratio, both analytically for a torus and numerically for a non-planar coil of a stellarator fusion device, the Helically Symmetric eXperiment (HSX). Because the integrands of these formulae develop fine structure as the minor radius becomes infinitely thin, we also develop a method of evaluating the self-force and self-inductance with even greater efficiency by integrating this sharp feature analytically. We demonstrate with this method that the self-force can be accurately computed for the HSX coil with as few as 12 grid points.

Applied Physics,Accelerator Physics,Classical Physics,Computational Physics,Plasma Physics

What problem does this paper attempt to address?

This paper aims to address several key issues encountered in the design of electromagnetic coils, including Lorentz forces, internal magnetic field, and self-inductance. Specifically, when calculating the effects produced by a coil's own current (such as self-force or self-inductance), traditional numerical methods become very difficult due to the presence of singularities. The paper proposes a new non-singular integral formula to efficiently compute these quantities and demonstrates good agreement between the new method and exact solutions under low aspect ratio conditions.

- **Lorentz forces**: Crucial for coil design as they determine the stress experienced by the coil and its supporting structure.
- **Internal magnetic field**: Used to determine the extent to which the coil is approaching the quench limit.
- **Self-inductance**: Determines the stored magnetic energy, which is particularly important for the design of fusion reactors.

The paper addresses the singularity problem in traditional methods by dividing the integration domain into "near" and "far" regions and using appropriate approximation methods to simplify high-dimensional integrals. Additionally, the authors propose an improved method that further enhances numerical computation efficiency, such as achieving high accuracy with only 12 grid points when calculating self-force.