Dra. Camaño-Valenzuela, Jessika

In this article, we analyse an augmented mixed finite element method for the steady Navier–Stokes equations. More precisely, we extend the recent results from Camaño et al.. (2017, Analysis of an augmented mixed-FEM for the Navier–Stokes problem. Math. Comput., 86, 589–615) to the case of mixed no-slip and traction boundary conditions in different parts of the boundary, and introduce and analyse a new pseudostress–velocity-augmented mixed formulation for the fluid flow problem. The well-posedness analysis is carried out by combining the classical Babuška–Brezzi theory and Banach’s fixed-point theorem. A proper adaptation of the arguments exploited in the continuous analysis allows us to state suitable hypotheses on the finite element subspaces ensuring that the associated Galerkin scheme is well defined. For instance, Raviart–Thomas elements of order k≥0 k≥0 and continuous piecewise polynomials of degree k+1 k+1 for the nonlinear pseudostress tensor and velocity, respectively, yield optimal convergence rates. In addition, we derive a reliable and efficient residual-based a posteriori error estimator for the proposed discretization. The proof of reliability hinges on the global inf–sup condition and the local approximation properties of the Clément interpolant, whereas the efficiency of the estimator follows from inverse inequalities and localization via edge–bubble functions. A set of numerical results exemplifies the performance of the augmented method with mixed boundary conditions. The tests also confirm the reliability and efficiency of the estimator, and show the performance of the associated adaptive algorithm.

A new stress-based mixed variational formulation for the stationary Navier-Stokes equations with constant density and variable viscosity depending on the magnitude of the strain tensor, is proposed and analyzed in this work. Our approach is a natural extension of a technique applied in a recent paper by some of the authors to the same boundary value problem but with a viscosity that depends nonlinearly on the gradient of velocity instead of the strain tensor. In this case, and besides remarking that the strain-dependence for the viscosity yields a more physically relevant model, we notice that to handle this nonlinearity we now need to incorporate not only the strain itself but also the vorticity as auxiliary unknowns. Furthermore, similarly as in that previous work, and aiming to deal with a suitable space for the velocity, the variational formulation is augmented with Galerkin-type terms arising from the constitutive and equilibrium equations, the relations defining the two additional unknowns, and the Dirichlet boundary condition. In this way, and as the resulting augmented scheme can be rewritten as a fixed-point operator equation, the classical Schauder and Banach theorems together with monotone operators theory are applied to derive the well-posedness of the continuous and associated discrete schemes. In particular, we show that arbitrary finite element subspaces can be utilized for the latter, and then we derive optimal a priori error estimates along with the corresponding rates of convergence. Next, a reliable and efficient residual-based a posteriori error estimator on arbitrary polygonal and polyhedral regions is proposed. The main tools used include Raviart-Thomas and Clément interpolation operators, inverse and discrete inequalities, and the localization technique based on triangle-bubble and edge-bubble functions. Finally, several numerical essays illustrating the good performance of the method, confirming the reliability and efficiency of the a posteriori error estimator, and showing the desired behavior of the adaptive algorithm, are reported.