High-accuracy embedded boundary grid generation using the divergence theorem

We present an algorithm to produce the necessary geometric information for finite volume calculations in the context of Cartesian grids with embedded boundaries. Given an order of accuracy for the overall calculation, we show what accuracy is required for each of the geometric quantities and we demonstrate how to calculate the moments using the divergence theorem. We demonstrate that, for a known flux, these moments can be used to create a flux divergence of the expected order.