  Strong Beltrami fields have long played a key role in fluid mechanics and
magnetohydrodynamics. In particular, they are the kind of stationary solutions
of the Euler equations where one has been able to show the existence of vortex
structures (vortex tubes and vortex lines) of arbitrarily complicated topology.
On the contrary, there are very few results about the existence of generalized
Beltrami fields, that is, divergence-free fields whose curl is the field itself
times a non-constant function. In fact, generalized Beltrami fields (which are
also stationary solutions to the Euler equations) have been recently shown to
be rare, in the sense that for "most" proportionality factors there are no
nontrivial Beltrami fields of high enough regularity (e.g., of class
$C^{6,\alpha}$), not even locally.
  We show that, nevertheless, there are "many" Beltrami fields with
non-constant factor, even realizing arbitrarily complicated vortex structures.
The core results are an "almost global" stability theorem for strong Beltrami
fields, which ensures that a global strong Beltrami field with suitable decay
at infinity can be perturbed to get "many" Beltrami fields with non-constant
factor of arbitrarily high regularity and defined in the exterior of an
arbitrarily small ball, and a "local" stability theorem for generalized
Beltrami fields, which is an analogous perturbative result that is valid for
any kind of Beltrami field (not just with a constant factor) but only applies
to small enough domains.
  The proof relies on an iterative scheme of Grad-Rubin type. For this purpose,
we study the Neumann problem for the inhomogeneous Beltrami equation in
exterior domains via a boundary integral equation method and we obtain H\"older
estimates, a sharp decay at infinity and some compactness properties for these
sequences of approximate solutions.
