1. A method for constructing a velocity vector field from a given set of fluxes between cells in a grid of cells, comprising: subdividing the cells into subcells; creating a set of equations for calculating the flux across faces between subcells; and using a divergence-free correction method to reduce the number of equations to be solved per cell.