A discontinuous Galerkin method for a diffuse-interface model of immiscible two-phase flows with soluble surfactant.

A numerical method using discontinuous polynomial approximations is formulated for solving a phase-field model of two immiscible fluids with a soluble surfactant. The proposed scheme is shown to decay the total free Helmholtz energy at the discrete level, which is consistent with the continuous model dynamics. The scheme recovers the Langmuir adsorption isotherms at equilibrium. Simulations of spinodal decomposition, flow through a cylinder and flow through a sequence of pore throats show the dynamics of the flow with and without surfactant. Finally the numerical method is used to simulate fluid flows in the pore space of Berea sandstone obtained by micro-CT imaging. [ABSTRACT FROM AUTHOR]