Variational multiscale modeling with discretely divergence-free subscales. (1st December 2020)