High-order simulation of foam enhanced oil recovery