Investigating the n=1 and n=2 error fields in W7-X using the newly accelerated FIELDLINES code

No fusion device can be created without any uncertainty; there is always a slight deviation from the geometric specification. These deviations can add up create a deviation of the magnetic field. This deviation is known as the (magnetic) error field. Correcting these error fields is desired as they cause asymmetries in the divertor loads and can thus cause damage to the device if they grow too large. These error fields can be defined by their toroidal (n) and poloidal number (m). The correction of the n = 1 and n = 2 fields in Wendelstein 7-X (W7-X) is investigated in this work. This investigation focuses on field line diffusion to the divertor, a proxy for divertor heat flux. Such work leverages the 25x speedup obtained through the implementation of a new particle-wall collision model. The n = 1 and n = 2 error fields of the as-built coils model of W7-X are corrected by scanning phase and amplitude of the trim and control coils. Reductions in the divertor load asymmetry by factors of four are demonstrated using error field correction. It is found that the as-built coils model has a significantly lower m⁄n = 1⁄1 error field than found in experiments.

Numerical Experimental for an Inverse Source Problem Method

This paper examines extensions of an iterative method for inverse evaluation of the source function for two elliptic systems. The method begins with a starting value for the undetermined source. Next, a background field and equations for the error field are obtained. 2-D domains are considered. This method is suitable for Helmholtz and Poisson operators. In the presence of finite-difference grid resolution, a varying amount of boundary data, and methods of filtering the noise in the boundary data and the noise intensity of the boundary data, the performance, accuracy, and iteration count of the algorithm are investigated. Keywords: Source, Inverse Problems, Poisson, Noise, Ill-Posedness, Well-Posed

Goal oriented error estimation in multi-scale shell element finite element problems

A major challenge with modern aircraft design is the occurrence of structural features of varied length scales. Structural stiffness can be accurately represented using homogenisation, however aspects such as the onset of failure may require information on more refined length scale for both metallic and composite components. This work considers the errors encountered in the coarse global models due to the mesh size and how these are propagated into detailed local sub-models. The error is calculated by a goal oriented error estimator, formulated by solving dual problems and Zienkiewicz-Zhu smooth field recovery. Specifically, the novel concept of this work is applying the goal oriented error estimator to shell elements and propagating this error field into the continuum sub-model. This methodology is tested on a simplified aluminium beam section with four different local feature designs, thereby illustrating the sensitivity to various local features with a common global setting. The simulations show that when the feature models only contained holes on the flange section, there was little sensitivity of the von Mises stress to the design modifications. However, when holes were added to the webbing section, there were large stress concentrations that predicted yielding. Despite this increase in nominal stress, the maximum error does not significantly change. However, the error field does change near the holes. A Monte Carlo simulation utilising marginal distributions is performed to show the robustness of the multi-scale analysis to uncertainty in the global error estimation as would be expected in experimental measurements. This shows a trade-off between Saint-Venant’s principle of the applied loading and stress concentrations on the feature model when investigating the response variance.

Figures of merit for stellarators near the magnetic axis

A new paradigm for rapid stellarator configuration design has been recently demonstrated, in which the shapes of quasisymmetric or omnigenous flux surfaces are computed directly using an expansion in small distance from the magnetic axis. To further develop this approach, here we derive several other quantities of interest that can be rapidly computed from this near-axis expansion. First, the $\boldsymbol {\nabla }\boldsymbol {B}$ and $\boldsymbol {\nabla }\boldsymbol {\nabla }\boldsymbol {B}$ tensors are computed, which can be used for direct derivative-based optimization of electromagnetic coil shapes to achieve the desired magnetic configuration. Moreover, if the norm of these tensors is large compared with the field strength for a given magnetic field, the field must have a short length scale, suggesting it may be hard to produce with coils that are suitably far away. Second, we evaluate the minor radius at which the flux surface shapes would become singular, providing a lower bound on the achievable aspect ratio. This bound is also shown to be related to an equilibrium beta limit. Finally, for configurations that are constructed to achieve a desired magnetic field strength to first order in the expansion, we compute the error field that arises due to second-order terms.