Simulation of nonlinear viscous fingering in miscible displacement

The nonlinear behavior of viscous fingering in miscible displacements is studied. A Fourier spectral method is used as the basic scheme for numerical simulation. In its simplest formulation, the problem can be reduced to two algebraic equations for flow quantities and a first‐order ordinary differential equation in time for the concentration. There are two parameters, the Peclet number (Pe) and mobility ratio (M), that determine the stability characteristics. The result shows that at short times, both the growth rate and the wavelength of fingers are in good agreement with predictions from our previous linear stability theory. However, as the time goes on, the nonlinear behavior of fingers becomes important. There are always a few dominant fingers that spread and shield the growth of other fingers. The spreading and shielding effects are caused by a spanwise secondary instability, and are aided by the transverse dispersion. It is shown that once a finger becomes large enough, the concentration gradient of...