Document Type


Publication Date



A radial basis function (RBF) method based on matrix-valued kernels is presented and analyzed for computing two types of vector decompositions on bounded domains: one where the normal component of the divergence-free part of the field is specified on the boundary, and one where the tangential component of the curl-free part of the field specified. These two decompositions can then be combined to obtain a full Helmholtz-Hodge decomposition of the field, i.e. the sum of divergence-free, curl-free, and harmonic fields. All decompositions are computed from samples of the field at (possibly scattered) nodes over the domain, and all boundary conditions are imposed on the vector fields, not their potentials, distinguishing this technique from many current methods. Sobolev-type error estimates for the various decompositions are provided and demonstrated with numerical examples.

Copyright Statement

This is a pre-copyedited, author-produced PDF of an article accepted for publication in IMA Journal of Numerical Analysis following peer review. The version of record:

Fuselier, E.J. and Wright, G.B.(2017). A Radial Basis Function Method for Computing Helmholtz-Hodge Decompositions. IMA Journal of Numerical Analysis, 37(2), 774-797

is available online at doi: 10.1093/imanum/drw027

Included in

Mathematics Commons