Dr. Caucao-Paillán, Sergio

In this paper we present and analyze a new mixed finite element method for the nonlinear problem given by the stationary convective Brinkman–Forchheimer equations with varying porosity. Our approach is based on the introduction of the pseudostress and the gradient of the porosity times the velocity, as further unknowns. As a consequence, we obtain a mixed variational formulation within a Banach spaces framework, with the velocity and the aforementioned tensors as the only unknowns. The pressure, the velocity gradient, the vorticity, and the shear stress can be computed afterwards via postprocessing formulae. A fixed-point strategy, along with monotone operators theory and the classical Banach theorem, are employed to prove the well-posedness of the continuous and discrete systems. Specific finite element subspaces satisfying the required discrete stability condition are defined, and optimal a priori error estimates are derived. Finally, several numerical examples illustrating the performance and flexibility of the method and confirming the theoretical rates of convergence, are reported.

In this paper we introduce and analyze new Banach spaces-based mixed finite element methods for the stationary nonlinear problem arising from the coupling of the convective Brinkman-Forchheimer equations with a double diffusion phenomenon. Besides the velocity and pressure variables, the symmetric stress and the skew-symmetric vorticity tensors are introduced as auxiliary unknowns of the fluid. Thus, the incompressibility condition allows to eliminate the pressure, which, along with the velocity gradient and the shear stress, can be computed afterwards via postprocessing formulae depending on the velocity and the aforementioned new tensors. Regarding the diffusive part of the coupled model, and additionally to the temperature and concentration of the solute, their gradients and pseudoheat/pseudodiffusion vectors are incorporated as further unknowns as well. The resulting mixed variational formulation, settled within a Banach spaces framework, consists of a nonlinear perturbation of, in turn, a nonlinearly perturbed saddle-point scheme, coupled with a usual saddle-point system. A fixed-point strategy, combined with classical and recent solvability results for suitable linearizations of the decoupled problems, including in particular, the Banach-Nečas-Babuška theorem and the Babuška-Brezzi theory, are employed to prove, jointly with the Banach fixed-point theorem, the well-posedness of the continuous and discrete formulations. Both PEERS and AFW elements of order l>0 for the fluid variables, and piecewise polynomials of degree

In this paper we consider a Banach spaces-based fully-mixed variational formulation that has been recently proposed for the coupling of the stationary Brinkman–Forchheimer and double-diffusion equations, and develop the first reliable and efficient residual-based a posteriori error estimator for the 2D and 3D versions of the associated mixed finite element scheme. For the reliability analysis, and due to the nonlinear nature of the problem, we employ the strong monotonicity of the operator involving the Forchheimer term, in addition to inf-sup conditions of some of the resulting bilinear forms, along with a stable Helmholtz decomposition in nonstandard Banach spaces, which, in turn, having been recently derived, constitutes another distinctive feature of the paper, and local approximation properties of the Raviart–Thomas and Clément interpolants. On the other hand, inverse inequalities, the localization technique through bubble functions, and known results from previous works, are the main tools yielding the efficiency estimate. Finally, several numerical examples confirming the theoretical properties of the estimator and illustrating the performance of the associated adaptive algorithms, are reported. In particular, the case of flow through a 2D porous media with an irregular channel networks is considered.

We propose and analyze a new mixed finite element method for the coupling of the Stokes equations with a transport problem modelled by a scalar nonlinear convection–diffusion problem. Our approach is based on the introduction of the Cauchy fluid stress and two vector unknowns involving the gradient and the total flux of the concentration. The introduction of these further unknowns lead to a mixed formulation in a Banach space framework in both Stokes and transport equations, where the aforementioned stress tensor and vector unknowns, together with the velocity and the concentration, are the main unknowns of the system. In this way, and differently from the techniques previously developed for this and related coupled problems, no augmentation procedure needs to be incorporated now into the formulation nor into the solvability analysis. The resulting non-augmented scheme is then written equivalently as a fixed-point equation, so that the well-known Banach theorem, combined with Babuška–Brezzi’s theory in Banach spaces, classical results on nonlinear monotone operators and certain regularity assumptions, are applied to prove the unique solvability of the continuous system. As for the associated Galerkin scheme, whose solvability is established similarly to the continuous case by using the Brouwer fixed-point theorem, we employ Raviart–Thomas approximations of order for the stress and total flux, and discontinuous piecewise polynomials of degree k for the velocity, concentration, and concentration gradient. With this choice of spaces, momentum is conserved in both Stokes and transport equations if the external forces belong to the piecewise constants and concentration discrete space, respectively, which constitutes one of the main features of our approach. Finally, we derive optimal a priori error estimates and provide several numerical results illustrating the good performance of the scheme and confirming the theoretical rates of convergence.

In this paper we consider a partially augmented fully-mixed variational formulation that has been recently proposed for the coupling of the stationary Brinkman–Forchheimer and double-diffusion equations, and develop an a posteriori error analysis for the 2D and 3D versions of the associated mixed finite element scheme. Indeed, we derive two reliable and efficient residual-based a posteriori error estimators for this problem on arbitrary (convex or non-convex) polygonal and polyhedral regions. The reliability of the proposed estimators draws mainly upon the uniform ellipticity and inf-sup condition of the forms involved, a suitable assumption on the data, stable Helmholtz decompositions in Hilbert and Banach frameworks, and the local approximation properties of the Clément and Raviart–Thomas operators. In turn, inverse inequalities, the localization technique based on bubble functions, and known results from previous works, are the main tools yielding the efficiency estimate. Finally, several numerical examples confirming the theoretical properties of the estimators and illustrating the performance of the associated adaptive algorithms, are reported. In particular, the case of flow through a 3D porous media with channel networks is considered.

We propose and analyze a new mixed finite element method for the nonlinear problem given by the coupling of the steady Brinkman–Forchheimer and double-diffusion equations. Besides the velocity, temperature, and concentration, our approach introduces the velocity gradient, the pseudostress tensor, and a pair of vectors involving the temperature/concentration, its gradient and the velocity, as further unknowns. As a consequence, we obtain a fully mixed variational formulation presenting a Banach spaces framework in each set of equations. In this way, and differently from the techniques previously developed for this and related coupled problems, no augmentation procedure needs to be incorporated now into the formulation nor into the solvability analysis. The resulting non-augmented scheme is then written equivalently as a fixed-point equation, so that the well-known Banach theorem, combined with classical results on nonlinear monotone operators and Babuˇ ska–Brezzi’s theory in Banach spaces, are applied to prove the unique solvability of the continuous and discrete systems. Appropriate finite element subspaces satisfying the required discrete inf-sup conditions are specified, and optimal a priori error estimates are derived. Several numerical examples confirm the theoretical rates of convergence and illustrate the performance and flexibility of the method.

We propose and analyze a new mixed finite element method for the problem of steady double-diffusive convection in a fluid-saturated porous medium. More precisely, the model is described by the coupling of the Brinkman–Forchheimer and double-diffusion equations, in which the originally sought variables are the velocity and pressure of the fluid, and the temperature and concentration of a solute. Our approach is based on the introduction of the further unknowns given by the fluid pseudostress tensor, and the pseudoheat and pseudodiffusive vectors, thus yielding a fully-mixed formulation. Furthermore, since the nonlinear term in the Brinkman–Forchheimer equation requires the velocity to live in a smaller space than usual, we partially augment the variational formulation with suitable Galerkin type terms, which forces both the temperature and concentration scalar fields to live in \(\mathrm {L}^4\). As a consequence, the aforementioned pseudoheat and pseudodiffusive vectors live in a suitable \(\mathrm {H}(\mathrm {div})\)-type Banach space. The resulting augmented scheme is written equivalently as a fixed point equation, so that the well-known Schauder and Banach theorems, combined with the Lax–Milgram and Banach–Nečas–Babuška theorems, allow to prove the unique solvability of the continuous problem. As for the associated Galerkin scheme we utilize Raviart–Thomas spaces of order \(k\ge 0\) for approximating the pseudostress tensor, as well as the pseudoheat and pseudodiffusive vectors, whereas continuous piecewise polynomials of degree \(\le k + 1\) are employed for the velocity, and piecewise polynomials of degree \(\le k\) for the temperature and concentration fields. In turn, the existence and uniqueness of the discrete solution is established similarly to its continuous counterpart, applying in this case the Brouwer and Banach fixed-point theorems, respectively. Finally, we derive optimal a priori error estimates and provide several numerical results confirming the theoretical rates of convergence and illustrating the performance and flexibility of the method.

In this work we present and analyse a mixed finite element method for the coupling of fluid flow with porous media flow. The flows are governed by the Navier–Stokes and the Darcy–Forchheimer equations, respectively, and the corresponding transmission conditions are given by mass conservation, balance of normal forces, and the Beavers–Joseph–Saffman law. We consider the standard mixed formulation in the Navier–Stokes domain and the dual-mixed one in the Darcy–Forchheimer region, which yields the introduction of the trace of the porous medium pressure as a suitable Lagrange multiplier. The well-posedness of the problem is achieved by combining a fixed-point strategy, classical results on nonlinear monotone operators and the well-known Schauder and Banach fixed-point theorems. As for the associated Galerkin scheme we employ Bernardi–Raugel and Raviart–Thomas elements for the velocities, and piecewise constant elements for the pressures and the Lagrange multiplier, whereas its existence and uniqueness of solution is established similarly to its continuous counterpart, using in this case the Brouwer and Banach fixed-point theorems, respectively. We show stability, convergence, and a priori error estimates for the associated Galerkin scheme. Finally, we report some numerical examples confirming the predicted rates of convergence, and illustrating the performance of the method.