A finite volume/discontinuous Galerkin method for the advective Cahn–Hilliard equation with degenerate mobility on porous domains stemming from micro-CT imaging

Frank F, Liu C, Alpak FO, Rivière B (2018)


Publication Language: English

Publication Type: Journal article, Original article

Publication year: 2018

Journal

Book Volume: 22

Pages Range: 543–563

Journal Issue: 2

URI: https://link.springer.com/article/10.1007/s10596-017-9709-1

DOI: 10.1007/s10596-017-9709-1

Abstract

A numerical method is formulated for the solution of the advective Cahn–Hilliard (CH) equation with constant and degenerate mobility in three-dimensional porous media with non-vanishing velocity on the exterior boundary. The CH equation describes phase separation of an immiscible binary mixture at constant temperature in the presence of a mass constraint and dissipation of free energy. Porous media/pore-scale problems specifically entail images of rocks in which the solid matrix and pore spaces are fully resolved. The interior penalty discontinuous Galerkin method is used for the spatial discretization of the CH equation in mixed form, while a semi-implicit convex–concave splitting is utilized for temporal discretization. The spatial approximation order is arbitrary, while it reduces to a finite volume scheme for the choice of elementwise constants. The resulting nonlinear systems of equations are reduced using the Schur complement and solved via inexact Newton's method. The numerical scheme is first validated using numerical convergence tests and then applied to a number of fundamental problems for validation and numerical experimentation purposes including the case of degenerate mobility. First-order physical applicability and robustness of the numerical method are shown in a breakthrough scenario on a voxel set obtained from a micro-CT scan of a real sandstone rock sample.

Authors with CRIS profile

Involved external institutions

How to cite

APA:

Frank, F., Liu, C., Alpak, F.O., & Rivière, B. (2018). A finite volume/discontinuous Galerkin method for the advective Cahn–Hilliard equation with degenerate mobility on porous domains stemming from micro-CT imaging. Computational Geosciences, 22(2), 543–563. https://dx.doi.org/10.1007/s10596-017-9709-1

MLA:

Frank, Florian, et al. "A finite volume/discontinuous Galerkin method for the advective Cahn–Hilliard equation with degenerate mobility on porous domains stemming from micro-CT imaging." Computational Geosciences 22.2 (2018): 543–563.

BibTeX: Download