  Models of magnetohydrodynamic (MHD) equilibia that for computational
convenience assume the existence of a system of nested magnetic flux surfaces
tend to exhibit singular current sheets. These sheets are located on resonant
flux surfaces that are associated with rational values of the rotational
transform. We study the possibility of eliminating these singularities by
suitable modifications of the plasma boundary, which we prescribe in a fixed
boundary setting. We find that relatively straightforward iterative procedures
can be used to eliminate weak current sheets that are generated at resonant
flux surfaces by the nonlinear interactions of resonating wall harmonics. These
types of procedures may prove useful in the design of fusion devices with
configurations that enjoy improved stability and transport properties.
