To obtain engineering-feasible designs of stellarators with permanent magnets and simplified coils, a new algorithm has been developed based on Fourier decomposition and surface magnetic charges method. The strong toroidal fields in a quasi-axisymmetric stellarator are still generated by coils. The permanent magnets are designed to compensate the normal magnetic field Bn on the plasma surface ∂P created by the coils and plasma. The normal magnetic fields created by the permanent magnets Bpmn are calculated as the difference between the magnetic fields created by the surface magnetic charges on the inner surface ∂D and the outer surface ∂Dh of the magnets. The Fourier coefficients of the magnet thickness function are computed through matrix division operation based on the least square principle with dominant Fourier components selected through 2D Fourier transformation of Bpmn. The residual uncompensated Bn is minimized through iteration to progressively optimize the thickness function. This new algorithm has been successfully applied to design the permanent magnets of an l = 2 quasi-axisymmetric stellarator with background magnetic field created by 12 identical circular planar coils for demonstrations. High accuracy has been achieved, allowing for a flux-surface-averaged residual Bn relative to the total field and a maximum residual Bn of less than 2 Gs for ∼1 T total field. This new design has some advantages in engineering implementations: all permanent magnet pieces have the same remanence Br; only one single layer of magnets are mounted perpendicular to the winding surface; the magnets can be easily inserted into the cells of a gridded frame attached to the winding surface and fixed with springs from the back, which greatly simplifies the manufacture, assembly and maintenance of the magnets, and thus facilitates precision control and cost reduction.
A novel scheme for error field correction in permanent magnet stellarators