Multiphase flow in fractured porous media can be described
by discrete fracture matrix models that represent the fractures as
dimensionally reduced manifolds embedded in the bulk porous medium.
Generalizing earlier work on this approach we focus on immiscible
two-phase flow in time-dependent fracture geometries, i.e., the fracture
itself and the aperture of the fractures might evolve in time. For
dynamic fracture geometries of that kind, neglecting capillary forces,
we deduce by transversal averaging of a full dimensional description a
dimensionally reduced model that governs the geometric evolution and the
flow dynamics. The core computational contribution is a
mixed-dimensional finite-volume discretization based on a conforming
moving-mesh ansatz. This finite-volume moving-mesh (FVMM) algorithm is
tracking the fractures' motions as a family of unions of facets of the
mesh. Notably, the method permits arbitrary movement of facets of the
triangulation while keeping the mass conservation constraint. In a
series of numerical examples we investigate the modeling error of the
reduced model as it compares to the original full dimensional model.
Moreover, we show the performance of the finite-volume moving-mesh
algorithm for the complex wave pattern that is induced by the
interaction of saturation fronts and evolving fractures.
%0 Journal Article
%1 burbulla2022finitevolume
%A Burbulla, Samuel
%A Rohde, Christian
%D 2022
%J Journal of Computational Physics
%K imported from:brittalenz vorlaeufig ians
%P 111031
%R https://doi.org/10.1016/j.jcp.2022.111031
%T A finite-volume moving-mesh method for two-phase flow in
fracturing porous media
%U https://www.sciencedirect.com/science/article/pii/S0021999122000936
%X Multiphase flow in fractured porous media can be described
by discrete fracture matrix models that represent the fractures as
dimensionally reduced manifolds embedded in the bulk porous medium.
Generalizing earlier work on this approach we focus on immiscible
two-phase flow in time-dependent fracture geometries, i.e., the fracture
itself and the aperture of the fractures might evolve in time. For
dynamic fracture geometries of that kind, neglecting capillary forces,
we deduce by transversal averaging of a full dimensional description a
dimensionally reduced model that governs the geometric evolution and the
flow dynamics. The core computational contribution is a
mixed-dimensional finite-volume discretization based on a conforming
moving-mesh ansatz. This finite-volume moving-mesh (FVMM) algorithm is
tracking the fractures' motions as a family of unions of facets of the
mesh. Notably, the method permits arbitrary movement of facets of the
triangulation while keeping the mass conservation constraint. In a
series of numerical examples we investigate the modeling error of the
reduced model as it compares to the original full dimensional model.
Moreover, we show the performance of the finite-volume moving-mesh
algorithm for the complex wave pattern that is induced by the
interaction of saturation fronts and evolving fractures.
@article{burbulla2022finitevolume,
abstract = {Multiphase flow in fractured porous media can be described
by discrete fracture matrix models that represent the fractures as
dimensionally reduced manifolds embedded in the bulk porous medium.
Generalizing earlier work on this approach we focus on immiscible
two-phase flow in time-dependent fracture geometries, i.e., the fracture
itself and the aperture of the fractures might evolve in time. For
dynamic fracture geometries of that kind, neglecting capillary forces,
we deduce by transversal averaging of a full dimensional description a
dimensionally reduced model that governs the geometric evolution and the
flow dynamics. The core computational contribution is a
mixed-dimensional finite-volume discretization based on a conforming
moving-mesh ansatz. This finite-volume moving-mesh (FVMM) algorithm is
tracking the fractures' motions as a family of unions of facets of the
mesh. Notably, the method permits arbitrary movement of facets of the
triangulation while keeping the mass conservation constraint. In a
series of numerical examples we investigate the modeling error of the
reduced model as it compares to the original full dimensional model.
Moreover, we show the performance of the finite-volume moving-mesh
algorithm for the complex wave pattern that is induced by the
interaction of saturation fronts and evolving fractures.},
added-at = {2022-03-08T17:42:35.000+0100},
author = {Burbulla, Samuel and Rohde, Christian},
biburl = {https://puma.ub.uni-stuttgart.de/bibtex/2e8aad2deb8c8d600c820b08634b83cfd/mathematik},
doi = {https://doi.org/10.1016/j.jcp.2022.111031},
interhash = {1ef19dc47df248ebcee76f8655d00172},
intrahash = {e8aad2deb8c8d600c820b08634b83cfd},
issn = {0021-9991},
journal = {Journal of Computational Physics},
keywords = {imported from:brittalenz vorlaeufig ians},
pages = 111031,
timestamp = {2022-03-08T16:42:35.000+0100},
title = {A finite-volume moving-mesh method for two-phase flow in
fracturing porous media},
url = {https://www.sciencedirect.com/science/article/pii/S0021999122000936},
year = 2022
}