Modeling of unsteady flows of multiphase viscous fluid in a pore space

Abstract

The authors have developed and implemented a numerical algorithm to model unsteady flows of a viscous multiphase isothermal fluid by finite difference method using the projection method for the numerical solution of the Navier-Stokes equation. The projection method implies splitting the initial system of equations by physical processes, in which convective transport and the effect of the pressure gradient are separately taken into account. As a result, at each step, it is necessary to solve the Poisson equation to find the pressure field. The solution of SLAE is performed by a parallel direct solver based on LU decomposition. An explicit scheme is used to solve the Cahn-Hilliard equation to update the phase field, the parameter of which is taken into account when adding surface forces to the Navier-Stokes equation. Computational experiments showing qualitative and quantitative agreement with experimental and numerical data from the literature are presented.