Abstract:In this paper we employ a Banach spaces-based framework to introduce and analyze new mixed finite element methods for the numerical solution of the coupled Stokes and Poisson–Nernst–Planck equations, which is a nonlinear model describing the dynamics of electrically charged incompressible fluids. The pressure of the fluid is eliminated from the system (though computed afterwards via a postprocessing formula) thanks to the incompressibility condition and the incorporation of the fluid pseudostress as an auxiliary unknown. In turn, besides the electrostatic potential and the concentration of ionized particles, we use the electric field (rescaled gradient of the potential) and total ionic fluxes as new unknowns. The resulting fully mixed variational formulation in Banach spaces can be written as a coupled system consisting of two saddle-point problems, each one with nonlinear source terms depending on the remaining unknowns, and a perturbed saddlepoint problem with linear source terms, which is in turn additionally perturbed by a bilinear form. The well-posedness of the continuous formulation is a consequence of a fixed-point strategy in combination with the Banach theorem, the Babuˇska–Brezzi theory, the solvability of abstract perturbed saddle-point problems, and the Banach–Neˇcas–Babuˇska theorem. For this we also employ smallness assumptions on the data. An analogous approach, but using now both the Brouwer and Banach theorems, and invoking suitable stability conditions on arbitrary finite element subspaces, is employed to conclude the existence and uniqueness of solution for the associated Galerkin scheme. A priori error estimates are derived, and examples of discrete spaces that fit the theory, include, e.g., Raviart–Thomas elements of order k along with piecewise polynomials of degree ďk. Finally, rates of convergence are specified and several numerical experiments confirm the theoretical error bounds. These tests also illustrate the balance-preserving properties and applicability of the proposed family of methods.

A Mixed Finite Element Method for Nearly Incompressible Multiple-Network Poroelasticity

Mixed finite element method for saturated poroelastoplastic media at large strains

Unconditionally stable second order convergent partitioned methods for multiple-network poroelasticity

Mixed and multipoint finite element methods for rotation-based poroelasticity

A Generalized Multiscale Finite Element Method for Poroelasticity Problems I: Linear Problems

A new multiphysics finite element method for a quasi-static poroelasticity model

A priori and a posteriori error bounds for the fully mixed FEM formulation of poroelasticity with stress-dependent permeability

A new mixed finite element method for arbitrary element pair for a quasi-static nonlinear permeability thermo-poroelasticity model

A mixed nonlocal finite element model for thermo‐poro‐elasto‐plastic simulation of porous media with multiphase fluid flow

A mixed stabilized finite element formulation for finite deformation of a poroelastic solid saturated with a compressible fluid

Error Estimates of a Fully Discrete Multiphysics Finite Element Method for a Nonlinear Poroelasticity Model

A topology-motivated mixed finite element method for dynamic response of porous media

A stabilized total pressure-formulation of the Biot's poroelasticity equations in frequency domain: numerical analysis and applications

A finite element implementation of finite deformation surface and bulk poroelasticity

Multiphysics mixed finite element method with Nitsche's technique for Stokes poroelasticity problem

Performance of mixed formulations for the particle finite element method in soil mechanics problems

Multi-compartment poroelastic models of perfused biological soft tissues: implementation in FEniCSx

Network-inspired versus Kozeny-Carman based permeability-porosity relations applied to Biot's poroelasticity model

New mixed finite element methods for the coupled Stokes and Poisson–Nernst–Planck equations in Banach spaces

Analysis of multiphysics finite element method for quasi-static thermo-poroelasticity with a nonlinear convective transport term

A Lagrange Multiplier-based method for Stokes-linearized poro-hyperelastic interface problems