Modeling flow in porous media with double porosity/permeability: A stabilized mixed formulation, error analysis, and numerical solutions

The flow of incompressible fluids through porous media plays a crucial role in many technological applications such as enhanced oil recovery and geological carbon-dioxide sequestration. The flow within numerous natural and synthetic porous materials that contain multiple scales of pores cannot be adequately described by the classical Darcy equations. It is for this reason that mathematical models for fluid flow in media with multiple scales of pores have been proposed in the literature. However, these models are analytically intractable for realistic problems. In this paper, a stabilized mixed four-field finite element formulation is presented to study the flow of an incompressible fluid in porous media exhibiting double porosity/permeability. The stabilization terms and the stabilization parameters are derived in a mathematically and thermodynamically consistent manner, and the computationally convenient equal-order interpolation of all the field variables is shown to be stable. A systematic error analysis is performed on the resulting stabilized weak formulation. Representative problems, patch tests and numerical convergence analyses are performed to illustrate the performance and convergence behavior of the proposed mixed formulation in the discrete setting. The accuracy of numerical solutions is assessed using the mathematical properties satisfied by the solutions of this double porosity/permeability model. Moreover, it is shown that the proposed framework can perform well under transient conditions and that it can capture well-known instabilities such as viscous fingering

Finite element approximation of the viscoelastic flow problem: a non-residual based stabilized formulation

In this paper, a three-field finite element stabilized formulation for the incompressible viscoelastic fluid flow problem is tested numerically. Starting from a residual based formulation, a non-residual based one is designed, the benefits of which are highlighted in this work. Both formulations allow one to deal with the convective nature of the problem and to use equal interpolation for the problem unknowns View the MathML sources-u-p (deviatoric stress, velocity and pressure). Additionally, some results from the numerical analysis of the formulation are stated. Numerical examples are presented to show the robustness of the method, which include the classical 4: 1 planar contraction problem and the flow over a confined cylinder case, as well as a two-fluid formulation for the planar jet buckling problem.Peer ReviewedPostprint (author's final draft

Numerical analysis of a stabilized finite element approximation for the three-field linearized viscoelastic fluid problem using arbitrary interpolations

The original publication is available at www.esaimm2an.org.In this paper we present the numerical analysis of a three-field stabilized finite element formulation recently proposed to approximate viscoelastic flows. The three-field viscoelastic fluid flow problem may suffer from two types of numerical instabilities: on the one hand we have the two inf-sup conditions related to the mixed nature problem and, on the other, the convective nature of the momentum and constitutive equations may produce global and local oscillations in the numerical approximation. Both can be overcome by resorting from the standard Galerkin method to a stabilized formulation. The one presented here is based on the subgrid scale concept, in which unresolvable scales of the continuous solution are approximately accounted for. In particular, the approach developed herein is based on the decomposition into their finite element component and a subscale, which is approximated properly to yield a stable formulation. The analyzed problem corresponds to a linearized version of the Navier-Stokes/Oldroyd-B case where the advection velocity of the momentum equation and the non-linear terms in the constitutive equation are treated using a fixed point strategy for the velocity and the velocity gradient. The proposed method permits the resolution of the problem using arbitrary interpolations for all the unknowns. We describe some important ingredients related to the design of the formulation and present the results of its numerical analysis. It is shown that the formulation is stable and optimally convergent for small Weissenberg numbers, independently of the interpolation used.Peer ReviewedPostprint (author's final draft

Dynamical and cohomological obstructions to extending group actions

We study cohomological obstructions to extending group actions on the boundary $\partial M$ of a $3$-manifold to a $C^0$-action on $M$ when $\partial M$ is diffeomorphic to a torus or a sphere. In particular, we show that for a $3$-manifold $M$ with torus boundary which is not diffeomorphic to a solid torus, the torus action on the boundary does not extend to a $C^0$-action on $M$.Comment: Minor correction to statement of Theorem 1.

A local to global argument on low dimensional manifolds

For an oriented manifold $M$ whose dimension is less than $4$, we use the contractibility of certain complexes associated to its submanifolds to cut $M$ into simpler pieces in order to do local to global arguments. In particular, in these dimensions, we give a different proof of a deep theorem of Thurston in foliation theory which says that the natural map between classifying spaces $\mathrm{B}\text{Homeo}^{\delta}(M)\to \mathrm{B}\text{Homeo}(M)$ induces a homology isomorphism where $\text{Homeo}^{\delta}(M)$ denotes the group of homeomorphisms of $M$ made discrete. Our proof shows that in low dimensions, Thurston's theorem can be proved without using foliation theory. Finally, we show that this technique gives a new perspective on the homotopy type of homeomorphism groups in low dimensions. In particular, we give a different proof of Hacher's theorem that the homeomorphism groups of Haken $3$-manifolds with boundary are homotopically discrete without using his disjunction techniques.Comment: Thoroughly revised. To appear in Transactions of the AM