An Immersed Boundary Method with Divergence-Free Interpolation for Unstructured Polyhedral Grids