Summary
A 2D miscible displacement computer code based on a pseudospectral method that uses Chebyshev expansions and a grid adaptive to the position of the front was developed. In this method, the effect of physical dispersion is explicitly included, while the effect of numerical dispersion is greatly suppressed. The code was shown to have exceptional convergence, accuracy, stability, and computational economy. It accurately predicts the position and sharpness of the displacement front for ID miscible and immiscible displacements. The code also predicts oil recoveries and breakthrough times for 2D miscible displacements with unfavorable mobility ratios that are in agreement with experimental data and is shown to be effective in resolving details of viscous fingering phenomena.
