Dr. Caucao-Paillán, Sergio

In this paper we consider a mixed variational formulation that have been recently proposed for the coupling of the Navier–Stokes and Darcy–Forchheimer equations, and derive, though in a non-standard sense, a reliable and efficient residual-baseda posteriorierror estimator suitable for an adaptive mesh-refinement method. For the reliability estimate, which holds with respect to the square root of the error estimator, we make use of the inf-sup condition and the strict monotonicity of the operators involved, a suitable Helmholtz decomposition in non-standard Banach spaces in the porous medium, local approximation properties of the Clément interpolant and Raviart–Thomas operator, and a smallness assumption on the data. In turn, inverse inequalities, the localization technique based on triangle-bubble and edge-bubble functions in localLpspaces, are the main tools for developing the efficiency analysis, which is valid for the error estimator itself up to a suitable additional error term. Finally, several numerical results confirming the properties of the estimator and illustrating the performance of the associated adaptive algorithm are reported.

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.