Abstract
Summary: This paper presents a finite-difference scheme for numerical modeling of the convection-dominated flows common in miscible or chemical flooding in porous media. The proposed difference operator possesses desirable phase and amplitude properties capable of resolving sharp fronts with negligible over- and undershoot. The scheme is generalized to multidimensions. Difference analogs have been derived for variable grid sizes for large-scale compositional simulations. Application to 3D miscible displacement problems and comparison with other standard industry methods have demonstrated the scheme's effectiveness in controlling numerical diffusion and grid-orientation sensitivity without any significant additional computational effort.