Characteristic block-centered finite difference method for simulating incompressible wormhole propagation

In this paper, the characteristic block-centered finite difference method is introduced and analyzed to solve the incompressible wormhole propagation. Error estimates for the pressure, velocity, porosity, concentration and its flux in different discrete norms are established rigorously and carefully on non-uniform grids. Finally, some numerical experiments are presented to show that the convergence rates are in agreement with the theoretical analysis.