Implicit treatment of molar mass equations in secondary oil migration

The hydrocarbon migration can be described by a coupled set of partial differential equations describing the dynamics of the temperature, component flow, pressure and velocity. A sequential solution procedure where the component flow is solved explicitly, gives severe restrictions on the time step given by the CFL condition. In this paper an implicit solution procedure is given and results from numerical tests are presented. The results are compared with the explicit solutions. As expected the implicit algorithm allows for substantially larger time steps.