Assuming that the grid is orthogonal, i.e.  that any centroid of volumes is such that the segment joining the centroids of two adjacent volumes is orthogonal to the face they across, except, maybe, for the boundary volumes, the above equation can be discretised as: