Multilevel interpolation of divergence-free vector fields

Patricio Farrell, Kathryn Gillow, Holger Wendland

Research output: Contribution to journalArticlepeer-review

10 Scopus citations


We introduce a multilevel technique for interpolating scattered data of divergence-free vector fields with the help of matrix-valued compactly supported kernels. The support radius at a given level is linked to the mesh norm of the data set at that level. There are at least three advantages of this method: no grid structure is necessary for the implementation, the multilevel approach is computationally cheaper than solving a large one-shot system and the interpolant is guaranteed to be analytically divergence-free. Furthermore, though we will not pursue this here, our multilevel approach is able to represent multiple scales in the data if present. We will prove convergence of the scheme, stability estimates and give a numerical example. For the first time, we will also prove error estimates for derivatives and give approximation orders in terms of the fill distance of the finest data set.
Original languageEnglish (US)
Pages (from-to)332-353
Number of pages22
JournalIMA Journal of Numerical Analysis
Issue number1
StatePublished - May 2 2016
Externally publishedYes

Bibliographical note

KAUST Repository Item: Exported on 2022-06-03
Acknowledged KAUST grant number(s): KUK-C1-013-04
Acknowledgements: This publication was based on work supported in part by Award No KUK-C1-013-04, made by King Abdullah University of Science and Technology (KAUST).
This publication acknowledges KAUST support, but has no KAUST affiliated authors.

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics
  • General Mathematics


Dive into the research topics of 'Multilevel interpolation of divergence-free vector fields'. Together they form a unique fingerprint.

Cite this