A method is presented to establish regions of phase space for 3D vector fields through which pass no co-oriented invariant 2D submanifolds transverse to a given oriented 1D foliation. Refinements are given for the cases of volume-preserving or Cartan–Arnol’d Hamiltonian flows and for boundaryless submanifolds.

The stepped-pressure equilibrium code (SPEC) (Hudson et al 2012 Phys. Plasmas 19, 112 502) is extended to enable free-boundary multi-region relaxed magnetohydrodynamic (MRxMHD) equilibrium calculations. The vacuum field surrounding the plasma inside an arbitrary 'computational boundary', ${\mathcal{D}}$, is computed, and the virtual-casing…