Approximation of the grad-div operator in non-convex domains