but as noted in Section II, standard finite-difference approximations would lead to a smearing of the F function and interfaces would lose their definition. Fortunately, the fact that F is a step function with values of zero or one permits the use of a flux approximation that preserves its discontinuous nature.