May 22 – 25, 2023
Europe/London timezone

Global implicit solver for multiphase multicomponent flow in porous media with multiple gas components and general reactions

May 24, 2023, 2:30 PM
Oral Presentation (MS07) Mathematical and numerical methods for multi-scale multi-physics, nonlinear coupled processes MS07


Dr Markus Knodel (Simulation in Technology TechSim)


In order to study the efficiency of the various forms of trapping including mineral trapping scenarios for CO$_2$ storage behavior in deep layers of porous media, highly non-linear coupled diffusion-advection-reaction partial differential equations (PDEs) including kinetic and equilibrium reactions modeling the miscible multiphase multicomponent flow have to be solved. We apply the globally fully implicit PDE reduction method (PRM) developed 2007 by Kräutle and Knabner for one-phase flow, which was extended 2019 to the case of two-phase flow with a pure gas in the study of Brunner and Knabner. We extend the method to the case of an arbitrary number of gases in gaseous phase, because CO$_2$ is not the only gas that threats the climate, and usually is accompanied by other climate killing gases. The application of the PRM leads to an equation system consisting of PDEs, ordinary differential equations, and algebraic equations. The Finite Element discretized / Finite Volume stabilized equations are separated into a local and a global system but nevertheless coupled by the resolution function and evaluated with the aid of a nested Newton solver, so our solver is fully global implicit. For the phase disappearance, we use persistent variables which lead to a semismooth formulation that is solved with a semismooth Newton method. We present scenarios of the injection of a mixture of various gases into deep layers, we investigate phase change effects in the context of various gases, and study the mineral trapping effects of the storage technique. The technical framework also applies to other fields such as nuclear waste storage or oil recovery.


Global implicit solver for multiphase multicomponent flow in porous media with multiple gas components and general reactions. Comput Geosci 26, 697-724 (2022). M.M. Knodel, S. Kräutle, and P. Knabner.

Participation In-Person
Country Germany
Energy Transition Focused Abstracts This abstract is related to Energy Transition
MDPI Energies Student Poster Award No, do not submit my presenation for the student posters award.
Acceptance of the Terms & Conditions Click here to agree

Primary authors

Dr Markus Knodel (Simulation in Technology TechSim) Prof. Serge Kräutle (Universität Erlangen-Nürnberg, Department Mathematik) Prof. Peter Knabner (Universität Erlangen-Nürnberg, Department Mathematik)

Presentation materials

There are no materials yet.