Florin A. Radu

Manuel Borregales, Florin A. Radu, Kundan Kumar et al.

We consider a non-linear extension of Biot's model for poromechanics, wherein both the fluid flow and mechanical deformation are allowed to be non-linear. We perform an implicit discretization in time (backward Euler) and propose two iterative schemes for solving the non-linear problems appearing within each time step: a splitting algorithm extending the undrained split and fixed stress methods to non-linear problems, and a monolithic L-scheme. The convergence of both schemes is shown rigorously. Illustrative numerical examples are presented to confirm the applicability of the schemes and validate the theoretical results.

David Landa-Marbán, Florin A. Radu, Jan M. Nordbotten

The focus of this paper is the derivation of a non-standard model for microbial enhanced oil recovery (MEOR) that includes the interfacial area (IFA) between the oil and water. We consider the continuity equations for water and oil, a balance equation for the oil-water interface and advective-dispersive transport equations for bacteria, nutrients and surfactants. Surfactants lower the interfacial tension (IFT), which improves the oil recovery. Therefore, we include in the model parameterizations of the IFT reduction and residual oil saturation as a function of the surfactant concentration. We consider for the first time in context of MEOR, the role of IFA in enhanced oil recovery (EOR). The motivation to include the IFA in the model is to reduce the hysteresis in the capillary pressure relationship, include the effects of observed bacteria migration towards the oil-water interface and biological production of surfactants at the oil-water interface. An efficient and robust linearization scheme was implemented, in which we use an implicit scheme that considers a linear approximation of the capillary pressure gradient, resulting in an efficient and stable scheme. A comprehensive, 2D implementation based on two-point flux approximation (TPFA) has been achieved. Illustrative numerical simulations are presented. We give an explanation of the differences in the oil recovery profiles obtained when we consider the IFA and MEOR effects. The model can also be used to design new experiments in order to gain a better understanding and optimization of MEOR.

Markus Bause, Florin A. Radu, Uwe Köcher

Variational time discretization schemes are getting of increasing importance for the accurate numerical approximation of transient phenomena. The applicability and value of mixed finite element methods (MFEM) in space for simulating transport processes have been demonstrated in a wide class of works. We consider a family of continuous Galerkin-Petrov time discretization schemes that is combined with a mixed finite element (MFE) approximation of the spatial variables. The existence and uniqueness of the semidiscrete approximation and of the fully discrete solution are established. For this, the Banach-Nečas-Babuška theorem is applied in a non-standard way. Error estimates with explicit rates of convergence are proved for the scalar and vector-valued variable. An optimal order estimate in space and time is proved by duality techniques for the scalar variable. The convergence rates are analyzed and illustrated by numerical experiments, also on stochastically perturbed meshes.

Jakub W. Both, Kundan Kumar, Jan M. Nordbotten et al.

Mathematical models for flow and reactive transport in porous media often involve non-linear, degenerate parabolic equations. Their solutions have low regularity, and therefore lower order schemes are used for the numerical approximation. Here the backward Euler method is combined with a mixed finite element method scheme, which results in a stable and locally mass-conservative scheme. At the same time, at each time step one has to solve a non-linear algebraic system, for which linear iterations are needed. Finding robust and convergent ones is particularly challenging here, since both slow and fast diffusion cases are allowed. Commonly used schemes, like Newton and Picard iterations, are defined either for non-degenerate problems, or after regularising the problem in the case of degenerate ones. Convergence is guaranteed only if the initial guess is sufficiently close to the solution, which translates into severe restrictions on the time step. Here we discuss a linear iterative scheme which builds on the $L$-scheme, and does not employ any regularisation. We prove its rigourous convergence, which is obtained for mild restrictions on the time step. Finally, we give numerical results confirming the theoretical ones, and compare the behaviour of the scheme with other schemes.

Markus Bause, Uwe Köcher, Florin A. Radu et al.

We introduce and analyze a post-processing for a family of variational space-time approximations to wave problems. The discretization in space and time is based on continuous finite element methods. The post-processing lifts the fully discrete approximations in time from continuous to continuously differentiable ones. Further, it increases the order of convergence of the discretization in time which can be be exploited nicely, for instance, for a-posteriori error control. The convergence behavior is shown by proving error estimates of optimal order in various norms. A bound of superconvergence at the discrete times nodes is included. To show the error estimates, a special approach is developed. Firstly, error estimates for the time derivative of the post-processed solution are proved. Then, in a second step these results are used to establish the desired error estimates for the post-processed solution itself. The need for this approach comes through the structure of the wave equation providing only stability estimates that preclude us from using absorption arguments for the control of certain error quantities. A further key ingredient of this work is the construction of a new time-interpolate of the exact solution that is needed in an essential way for deriving the error estimates. Finally, a conservation of energy property is shown for the post-processed solution which is a key feature for approximation schemes to wave equations. The error estimates given in this work are confirmed by numerical experiments.

David Seus, Florin A. Radu, Christian Rohde

This article is a follow up of our submitted paper [11] in which a decomposition of the Richards equation along two soil layers was discussed. A decomposed problem was formulated and a decoupling and linearisation technique was presented to solve the problem in each time step in a fixed point type iteration. This article extends these ideas to the case of two-phase in porous media and the convergence of the proposed domain decomposition method is rigorously shown.

Laura Portero, Andrés Arrarás, Francisco J. Gaspar et al.

We study the slightly compressible Darcy-Forchheimer equations modeling gas flow in porous media, particularly in applications related to combustion processes. The equations are discretized in time using the backward Euler method and in space via a mixed finite element scheme. As a result, a nonlinear algebraic system is obtained at each time step. We propose and analyze a general iterative linearization scheme for the efficient solution of such systems and study its convergence properties at the discrete level. The performance and robustness of the scheme are assessed through a series of numerical experiments. The method is compared with standard iterative solvers, and further tested on problems with discontinuous permeability fields. The results demonstrate its reliability and competitiveness in regimes characterized by strong nonlinear effects.