AN IMMERSED BOUNDARY METHOD WITH DIVERGENCE-FREE INTERPOLATION FOR UNSTRUCTURED POLYHEDRAL GRIDS