Source:Alim, Usman, University of Calgary, Number 2015-1079-12, Calgary, AB (2015)
Vector-field interpolation is a fundamental task in flow simulation and visualization. The common practice is to interpolate the vector field in a component-wise fashion. When the vector field of interest is solenoidal (divergence free), such an approach is not conservative and gives rise to artificial divergence. In this work, we numerically compare some recently proposed scalar interpolation methods on the Cartesian and body-centered cubic lattices, and investigate their ability to conserve the solenoidal nature of the vector field. We start with a sampled version of a synthetic solenoidal vector field and use an interpolative component-wise reconstruction method to approximate the vector field and its divergence at arbitrary locations. Our results show that an improved scalar interpolation method does not necessarily lead to a more conservative vector field approximation.