Gauss’ theorem may be used to convert the
integrated value of FUX over the control volume to boundary fluxes at its sides.
Then, the flux leaving one control volume will automatically be gained by the
adjacent volume and conservation during advection is guaranteed