Numerical simulation of partial saturation in sands induced by flow and chemical reactivity