Parallel CPR preconditioner for a multicomponent fluid flow filtration problem in a porous medium

The constrained pressure residual (CPR) preconditioning method is considered with regard to solution of systems with matrices appearing in discretization of PDE systems describing multicomponent fluid flow in porous media. New versions of algorithms are proposed. Numerical experiments using an actual parallel hydrodynamic simulator were performed for test and actual oil fields in Western Siberia, these experiments confirm the efficiency of the methods.     

Preconditioning techniques for a mixed Stokes/Darcy model in porous media applications

We study numerical methods for a mixed Stokes/Darcy model in porous media applications. The global model is composed of two different submodels in a fluid region and a porous media region, coupled through a set of interface conditions. The weak formulation of the coupled model is of a saddle point type. The mixed finite element discretization applied to the saddle point problem leads to a coupled, indefinite, and nonsymmetric linear system of algebraic equations. We apply the preconditioned GMRES method to solve the discrete system and are particularly interested in efficient and effective decoupled preconditioning techniques. Several decoupled preconditioners are proposed. Theoretical analysis and numerical experiments show the effectiveness and efficiency of the preconditioners. Effects of physical parameters on the convergence performance are also investigated.     

Forchheimer law derived by homogenization of gas flow in turbomachines

A general formulation of the homogenization problem of compressible fluid flow through a periodic porous material in turbomachines is presented here. This formulation is able to derive a Forchheimer law with a mean velocity dependent permeability as equivalent macroscopic behavior. To specify this permeability, additional flow problems are defined on the unit cell and solved by a mixed stabilized finite element discretization. The application of the Galerkin least-square (GLS) method requires the introduction of two stabilization terms with appropriate parameters. The mixed finite element discretization of these unit cell problems is finally outlined.     

Spectral discretization of Darcy’s equations with pressure dependent porosity

We consider the flow of a viscous incompressible fluid in a rigid homogeneous porous medium provided with boundary conditions on the pressure around a circular well. When the boundary pressure presents high variations, the permeability of the medium depends on the pressure, so that the model is nonlinear. We propose a spectral discretization of the resulting system of equations which takes into account the axisymmetry of the domain and of the flow. We prove optimal error estimates and present some numerical experiments which confirm the interest of the discretization.     

Hierarchical Preconditioning for PT AP-Systems arising from particulate flows

Particulate flows, i.e. flow of an incompressible, Newtonian carrier fluid loaded with (many) rigid bodies, play an important role in diverse technical applications. In [1] a finite element method was introduced to simulate such flows. The method relies on a splitting method and a subspace projection method to incorporate the rigid body motion of the particles in a so called one domain approach [2]. The resulting systems arising after discretization are ill conditioned in general. Thus preconditioning is mandatory for instance when using an iterative Krylov subspace methods. In this paper we present a Schur complement based preconditioner well suited for this type of application. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analysis of an Euler implicit‐mixed finite element scheme for reactive solute transport in porous media

In this article, we analyze an Euler implicit‐mixed finite element scheme for a porous media solute transport model. The transporting flux is not assumed given, but obtained by solving numerically the Richards equation, a model for subsurface fluid flow. We prove the convergence of the scheme by estimating the error in terms of the discretization parameters. In doing so we take into account the numerical error occurring in the approximation of the fluid flow. The article, is concluded by numerical experiments, which are in good agreement with the theoretical estimates. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Single‐phase flow in composite poroelastic media

The mathematical formulation and analysis of the Barenblatt–Biot model of elastic deformation and laminar flow in a heterogeneous porous medium is discussed. This describes consolidation processes in a fluid‐saturated double‐diffusion model of fractured rock. The model includes various degenerate cases, such as incompressible constituents or totally fissured components, and it is extended to include boundary conditions arising from partially exposed pores. The quasi‐static initial–boundary problem is shown to have a unique weak solution, and this solution is strong when the data are smoother. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical simulations of glacial rebound using preconditioned iterative solution methods

This paper discusses finite element discretization and preconditioning strategies for the iterative solution of nonsymmetric indefinite linear algebraic systems of equations arising in modelling of glacial rebound processes. Some numerical experiments for the purely elastic model setting are provided. Comparisons of the performance of the iterative solution method with a direct solution method are included as well.     

Newton method for reactive solute transport with equilibrium sorption in porous media

We present a mass conservative numerical scheme for reactive solute transport in porous media. The transport is modeled by a convection-diffusion-reaction equation, including equilibrium sorption. The scheme is based on the mixed finite element method (MFEM), more precisely the lowest-order Raviart-Thomas elements and one-step Euler implicit. The underlying fluid flow is described by the Richards equation, a possibly degenerate parabolic equation, which is also discretized by MFEM. This work is a continuation of Radu et al. (2008) and Radu et al. (2009) [1] and [2] where the algorithmic aspects of the scheme and the analysis of the discretization method are presented, respectively. Here we consider the Newton method for solving the fully discrete nonlinear systems arising on each time step after discretization. The convergence of the scheme is analyzed. In the case when the solute undergoes equilibrium sorption (of Freundlich type), the problem becomes degenerate and a regularization step is necessary. We derive sufficient conditions for the quadratic convergence of the Newton scheme.     

Python framework for hp-adaptive discontinuous Galerkin methods for two-phase flow in porous media

In this paper we present a framework for solving two-phase flow problems in porous media. The discretization is based on a Discontinuous Galerkin method and includes local grid adaptivity and local choice of polynomial degree. The method is implemented using the new Python frontend Dune-FemPy to the open source framework Dune. The code used for the simulations is made available as Jupyter notebook and can be used through a Docker container. We present a number of time stepping approaches ranging from a classical IMPES method to a fully coupled implicit scheme. The implementation of the discretization is very flexible allowing to test different formulations of the two-phase flow model and adaptation strategies.     

Simulation of uncompressible fluid flow through a porous media

Recently, a great interest has been focused for investigations about transport phenomena in disordered systems. One of the most treated topics is fluid flow through anisotropic materials due to the importance in many industrial processes like fluid flow in filters, membranes, walls, oil reservoirs, etc. In this work is described the formulation of a 2D mathematical model to simulate the fluid flow behavior through a porous media (PM) based on the solution of the continuity equation as a function of the Darcy’s law for a percolation system; which was reproduced using computational techniques reproduced using a random distribution of the porous media properties (porosity, permeability and saturation). The model displays the filling of a partially saturated porous media with a new injected fluid showing the non-defined advance front and dispersion of fluids phenomena.     

A computationally efficient modification of mixed finite element methods for flow problems with full transmissivity tensors

The mixed finite element method for approximately solving flow equations in porous media has received a good deal of attention in the literature. The main idea is to solve for the head/pressure and fluid velocity (Darcy velocity) simultaneously to obtain a higher order approximation of the fluid velocity. In the case of a diagonal transmissivity tensor the algebraic equations resulting from the discretization can be reduced to a system of algebraic equations for the head/pressure variable alone. This reduction results in a smaller number of unknows to be solved for in an iterative method such as preconditioned conjugate gradient method. The fluid velocity is then obtained from an algebraic relationship. In the case of full transmissivity tensor, the algebraic reduction is more difficult. This paper investigates some algorithms resulting from the modification of the mixed finite element that take advantage of the mixed finite element method for the diagonal tensor case. The resulting schemes are more efficient implementations that maintain the same order of accuracy as the original schemes. © 1993 John Wiley & Sons, Inc.     

Block ILU-preconditioners for problems of filtration of a multicomponent mixture in a porous medium

ILU class preconditioners (ILU(0), ILU(1), ILUT) employed for iterative algorithms for asymmetric linear systems with sparse matrices are considered. Test matrices used in this study originate from discretization of systems of partial differential equations describing a multicomponent fluid flow in porous media. New algorithms for block storage of matrices and block based ILU-factorization are described. This new integrated approach was tested on a wide range of matrices resulted from actual hydrodynamic simulations of oil fields in Western Siberia and had demonstrated significant reduction of computational time.     

A domain decomposition preconditioner for hermite collocation problems

We propose a preconditioning method for linear systems of equations arising from piecewise Hermite bicubic collocation applied to two‐dimensional elliptic PDEs with mixed boundary conditions. We construct an efficient, parallel preconditioner for the GMRES method. The main contribution of the article is a novel interface preconditioner derived in the framework of substructuring and employing a local Hermite collocation discretization for the interface subproblems based on a hybrid fine‐coarse mesh. Interface equations based on this mesh depend only weakly on unknowns associated with subdomains. The effectiveness of the proposed method is highlighted by numerical experiments that cover a variety of problems. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 135–151, 2003     

Preconditioning indefinite discretization matrices

Summary The finite element discretization of many elliptic boundary value problems leads to linear systems with positive definite and symmetric coefficient matrices. Many efficient preconditioners are known for these systems. We show that these preconditioning matrices can also be used for the linear systems arising from boundary value problems which are potentially indefinite due to lower order terms in the partial differential equation. Our main tool is a careful algebraic analysis of the condition numbers and the spectra of perturbed matrices which are preconditioned by the same matrices as in the unperturbed case.     

A Multilevel Method for the Solution of Time Dependent Optimal Transport

In this paper we present a new computationally efficient numerical scheme for the minimizing flow for the computation of the optimal $L_2$mass transport mapping using the fluid approach. We review the method and discuss its numerical properties. We then derive a new scaleable, efficient discretization and a solution technique for the problem and show that the problem is equivalent to a mixed form formulation of a nonlinear fluid flow in porous media. We demonstrate the effectiveness of our approach using a number of numerical experiments.     

Numerical modeling of fluid flow in anisotropic fractured porous media

A model of double porosity in the case of an anisotropic fractured porous medium is considered (Dmitriev, Maksimov; 2007). A function of fluid exchange between the fractures and porous blocks depending on flow direction is given. The flow function is based on the difference between the pressure gradients. This feature enables one to take into account anisotropic properties of filtration in a more general form. The results of numerical solving a model two-dimensional problem are presented. The computational algorithm is based on a finite-element space approximation and explicit-implicit time approximations.     

A cellular-automata method for studying porous media properties

A cellular-automata (CA) approach for investigating properties of porous media with tortuous channels and different smoothness of pore walls is proposed. This approach is aimed at combining two different CA models: the first one is intended for constructing the morphology of a porous material; the second, for simulating a fluid flow through it. The porous media morphology is obtained as a result of evolution of a cellular automaton, forming a “steady pattern.” The result is then used for simulating a fluid flow through a porous medium by applying the Lattice Gas CA model. The method has been tested on a small fragment of a porous material and implemented for investigating a carbon electrode of a hydrogen fuel cell on a multiprocessor cluster.