Numerical method for 3D two-component isothermal compressible flows with application to digital rock physics
A numerical algorithm is proposed for modeling two-component viscous compressible isothermal flows with surface effects in 3D domains of complex shape with a voxel representation of geometry. The regularized system of Navier-Stokes-Kahn-Hillard equations is used as a basic mathematical model. Modeling of the drop spreading on a flat substrate and the displacement of one liquid by another in the pore space of a real rock sample are carried out. The computational results demonstrate the applicability and good operability of the exploited system of equations, the corresponding explicit finite-difference scheme and algorithms for its implementation for the numerical solution of the class of problems under consideration.