Multidomain mixed variational analysis of transport flow through elastoviscoplastic porous media

Multidomain mixed nonlinear transport and flow phenomena through elastoviscoplastic porous media is variationally analyzed. Mixed variational formulations of the poro-mechanical system are established via composition duality methods, determining solvability results on the basis of duality principles. The conformation of the coupled physical system corresponds to constrained transport processes driven by a compressible Darcian flow, in a quasistatic elastoviscoplastic deformable subsurface porous media, modeled variationally by primal evolution mixed transport and consolidation, and dual evolution mixed flow and quasistatic deformation. For parallel computing, non-overlapping multidomain decomposition methods based on variational macro-hybridization, are presented and discussed, providing a natural multi-physics approach for the coupled transport flow and deformation system. For computational realizations, internal variational macro-hybrid mixed semi-discrete approximations are given, as well as primal and dual fully discrete semi-implicit time marching schemes. Furthermore, the corresponding coupled transport-flow-deformation system is concluded and analyzed, proposing natural resolution coupling techniques.     

Homogenization of a two‐phase incompressible fluid in crossflow filtration through a porous medium

We study the homogenization of a slow viscous two‐phase incompressible flow in a domain consisting of a free fluid domain, a porous medium, and the interface between them. We take into account the capillary forces on the fluid‐fluid interfaces. We construct boundary layers describing the flow at the interface between the free fluid and the porous medium. We derive a macroscopic model with a viscous two‐phase fluid in the free domain, a coupled Darcy law connecting two‐phase velocities in the porous medium, and boundary conditions at the permeable interface between the free fluid domain and the porous medium.     

Steady filtration problems with seawater intrusion: macro-hybrid primal mixed variational analysis

Macro-hybrid penalized, primal mixed continuous and discrete variational formulations, for steady filtration problems, with seawater intrusion, are studied. This is a mixed version of our previous paper on macro-hybrid penalized approximations (G. Alduncin, J. Esquivel-Avila, and N. Vera-Guzman, Steady filtration problems with seawater intrusion: Macro-hybrid penalized finite element approximations, Int. J. Numer. Meth. Fluids (2005) 49, pp. 935–957). Penalized pressure–velocity mixed variational formulations are introduced for non-overlapping domain decompositions, with vertical interfaces, of sections of coastal aquifers. Well-posedness and stability conditions are established at the continuous and discrete levels. Macro-hybrid primal mixed internal approximations are defined on independent subdomain grids, with transmission conditions imposed in a dual variational sense. Parallel relaxation penalty-duality algorithms are discussed from fixed-point characterizations, for iterative numerical resolution.     

Study of a numerical scheme for miscible two‐phase flow in porous media

We study the convergence of a finite volume scheme for a model of miscible two‐phase flow in porous media. In this model, one phase can dissolve into the other one. The convergence of the scheme is proved thanks to an estimate on the two pressures, which allows to prove some estimates on the discrete time derivative of some nonlinear functions of the unknowns. Monotony arguments allow to show some properties on the limits of these functions. A key point in the scheme is to use particular averaging formula for the dissolution function arising in the space term. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 723–748, 2014     

A second‐order backward difference time‐stepping scheme for penalized Navier–Stokes equations modeling filtration through porous media

We investigate the stability and convergence of a fully implicit, linearly extrapolated second‐order backward difference time‐stepping scheme for the penalized Navier–Stokes equations modeling filtration through porous media. In the penalization approach, an extended Navier–Stokes equation is used in the entire computational domain with suitable resistance terms to mimic the presence of porous medium. It is widely used as an alternative to the heterogeneous approach in which different types of partial differential equations (PDEs) are used in fluid and porous subregions along with suitable continuity conditions at the interface. However, the introduction of extra resistance terms makes the penalized Navier–Stokes equations more nonlinear. We prove that the linearly extrapolated scheme is unconditionally stable and derive optimal order error estimates without any stability condition. To show feasibility and applicability of the approach, it is used to numerically solve a passive control problem in which flow around a solid body is controlled by adding porous layers on the surface. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 681–705, 2016     

Analysis of a combined mixed finite element and discontinuous Galerkin method for incompressible two‐phase flow in porous media

We analyze a combined method consisting of the mixed finite element method for pressure equation and the discontinuous Galerkin method for saturation equation for the coupled system of incompressible two‐phase flow in porous media. The existence and uniqueness of numerical solutions are established under proper conditions by using a constructive approach. Optimal error estimates in L²(H¹) for saturation and in L ^∞ (H(div)) for velocity are derived. Copyright © 2013 John Wiley & Sons, Ltd.     

The upwind finite difference fractional steps methods for two‐phase compressible flow in porous media

The upwind finite difference fractional steps methods are put forward for the two‐phase compressible displacement problem. Some techniques, such as calculus of variations, multiplicative commutation rule of difference operators, decomposition of high‐order difference operators, and prior estimates, are adopted. Optimal order estimates in L² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of seawater intrusion and migration‐accumulation of oil resources. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 67–88, 2003