Once the flux has been computed by the above method, it is multiplied by the flux
boundary area to get the amount of fluid to be subtracted from the donor cell and
added to the acceptor cell. When the process is repeated for all cell boundaries in the
mesh, the resulting F values correspond to the time-advanced values satisfying
Eq. (10) and still sharply define all interfaces.