47 research outputs found

    Comparison results for the Stokes equations

    Full text link
    This paper enfolds a medius analysis for the Stokes equations and compares different finite element methods (FEMs). A first result is a best approximation result for a P1 non-conforming FEM. The main comparison result is that the error of the P2-P0-FEM is a lower bound to the error of the Bernardi-Raugel (or reduced P2-P0) FEM, which is a lower bound to the error of the P1 non-conforming FEM, and this is a lower bound to the error of the MINI-FEM. The paper discusses the converse direction, as well as other methods such as the discontinuous Galerkin and pseudostress FEMs. Furthermore this paper provides counterexamples for equivalent convergence when different pressure approximations are considered. The mathematical arguments are various conforming companions as well as the discrete inf-sup condition

    Augmented mixed finite element method for the Oseen problem: A priori and a posteriori error analyses

    Get PDF
    We propose a new augmented dual-mixed method for the Oseen problem based on the pseudostress–velocity formulation. The stabilized formulation is obtained by adding to the dual-mixed approach suitable least squares terms that arise from the constitutive and equilibrium equations. We prove that for appropriate values of the stabilization parameters, the new variational formulation and the corresponding Galerkin scheme are well-posed, and a Céa estimate holds for any finite element subspaces. We also provide the rate of convergence when each row of the pseudostress is approximated by Raviart–Thomas or Brezzi–Douglas–Marini elements and the velocity is approximated by continuous piecewise polynomials. Moreover, we derive a simple a posteriori error estimator of residual type that consists of two residual terms and prove that it is reliable and locally efficient. Finally, we include several numerical experiments that support the theoretical results.Dirección de Investigación of the Universidad Católica de la Santísima Concepción (Chile) y CONICYT-Chile FONDECYT; Ministerio de Ciencia e Innovación del Gobierno de España

    The enriched Crouzeix-Raviart elements are equivalent to the Raviart-Thomas elements

    Full text link
    For both the Poisson model problem and the Stokes problem in any dimension, this paper proves that the enriched Crouzeix-Raviart elements are actually identical to the first order Raviart-Thomas elements in the sense that they produce the same discrete stresses. This result improves the previous result in literature which, for two dimensions, states that the piecewise constant projection of the stress by the first order Raviart-Thomas element is equal to that by the Crouzeix-Raviart element. For the eigenvalue problem of Laplace operator, this paper proves that the error of the enriched Crouzeix-Raviart element is equivalent to that of the Raviart-Thomas element up to higher order terms

    An augmented mixed finite element method for the Navier-Stokes equations with variable viscosity

    Get PDF
    A new mixed variational formulation for the Navier–Stokes equations with constant density and variable viscosity depending nonlinearly on the gradient of velocity, is proposed and analyzed here. Our approach employs a technique previously applied to the stationary Boussinesq problem and to the Navier-Stokes equations with constant viscosity, which consists firstly of the introduction of a modified pseudostress tensor involving the diffusive and convective terms, and the pressure. Next, by using an equivalent statement suggested by the incompressibility condition, the pressure is eliminated, and in order to handle the nonlinear viscosity, the gradient of velocity is incorporated as an auxiliary unknown. Furthermore, since the convective term forces the velocity to live in a smaller space than usual, we overcome this difficulty by augmenting the variational formulation with suitable Galerkin-type terms arising from the constitutive and equilibrium equations, the aforementioned relation defining the additional unknown, and the Dirichlet boundary condition. The resulting augmented scheme is then written equivalently as a fixed point equation, and hence the well-known Schauder and Banach theorems, combined with classical results on bijective monotone operators, are applied to prove the unique solvability of the continuous and discrete systems. No discrete inf-sup conditions are required for the well-posedness of the Galerkin scheme, and hence arbitrary finite element subspaces of the respective continuous spaces can be utilized. In particular, given an integer k ≥ 0, piecewise polynomials of degree ≤ k for the gradient of velocity, Raviart-Thomas spaces of order k for the pseudostress, and continuous piecewise polynomials of degree ≤ k + 1 for the velocity, constitute feasible choices. Finally, optimal a priori error estimates are derived, and several numerical results illustrating the good performance of the augmented mixed finite element method and confirming the theoretical rates of convergence are reported.Comisión Nacional de Investigación Científica y Tecnológica (Chile)Universidad del Bío-BíoMinistry of Education, Youth and Sports of the Czech Republi

    Guaranteed error control for the pseudostress approximation of the Stokes equations

    Get PDF
    The pseudostress approximation of the Stokes equations rewrites the stationary Stokes equations with pure (but possibly inhomogeneous) Dirichlet boundary conditions as another (equivalent) mixed scheme based on a stress in H (div) and the velocity in L2. Any standard mixed finite element function space can be utilized for this mixed formulation, e.g. the Raviart-Thomas discretization which is related to the Crouzeix-Raviart nonconforming finite element scheme in the lowest-order case. The effective and guaranteed a posteriori error control for this nonconforming velocity-oriented discretization can be generalized to the error control of some piecewise quadratic velocity approximation that is related to the discrete pseudostress. The analysis allows for local inf-sup constants which can be chosen in a global partition to improve the estimation. Numerical examples provide strong evidence for an effective and guaranteed error control with very small overestimation factors even for domains with large anisotropy

    Analysis of an augmented mixed-FEM for the Navier-Stokes problem

    Get PDF
    In this paper we propose and analyze a new augmented mixed finite element method for the Navier-Stokes problem. Our approach is based on the introduction of a “nonlinearpseudostress” tensor linking the pseudostress tensor with the convective term, which leads to a mixed formulation with the nonlinear-pseudostress tensor and the velocity as the main unknowns of the system. Further variables of interest, such as the fluid pressure, the fluid vorticity and the fluid velocity gradient, can be easily approximated as a simple postprocess of the finite element solutions with the same rate of convergence. The resulting mixed formulation is augmented by introducing Galerkin least-squares type terms arising from the constitutive and equilibrium equations of the Navier-Stokes equations and from the Dirichlet boundary condition, which are multiplied by stabilization parameters that are chosen in such a way that the resulting continuous formulation becomes well-posed. Then, the classical Banach’s fixed point Theorem and Lax-Milgram’s Lemma are applied to prove well-posedness of the continuous problem. Similarly, we establish well-posedness and the corresponding Cea’s estimate of the associated Galerkin scheme considering any conforming finite element subspace for each unknown. In particular, the associated Galerkin scheme can be defined by employing Raviart-Thomas elements of degree k for the nonlinear-pseudostress tensor, and continuous piecewise polynomial elements of degree k + 1 for the velocity, which leads to an optimal convergent scheme. In addition, we provide two iterative methods to solve the corresponding nonlinear system of equations and analyze their convergence. Finally, several numerical results illustrating the good performance of the method are provided.Comisión Nacional de Investigación Científica y TecnológicaMinistry of Education, Youth and Sports of the Czech Republi

    Guaranteed error control for the pseudostress approximation of the Stokes equations

    Get PDF
    The pseudostress approximation of the Stokes equations rewrites the stationary Stokes equations with pure (but possibly inhomogeneous) Dirichlet boundary conditions as another (equivalent) mixed scheme based on a stress in H(div) and the velocity in L2L^2. Any standard mixed finite element function space can be utilized for this mixed formulation, e.g. the Raviart-Thomas discretization which is related to the Crouzeix-Raviart nonconforming finite element scheme in the lowest-order case. The effective and guaranteed a posteriori error control for this nonconforming velocity-oriented discretization can be generalized to the error control of some piecewise quadratic velocity approximation that is related to the discrete pseudostress. The analysis allows for local inf-sup constants which can be chosen in a global partition to improve the estimation. Numerical examples provide strong evidence for an effective and guaranteed error control with very small overestimation factors even for domains with large anisotropy

    A posteriori error analysis of an augmented mixed method for the Navier-Stokes equations with nonlinear viscosity

    Get PDF
    In this work we develop the a posteriori error analysis of an augmented mixed finite element method for the 2D and 3D versions of the Navier-Stokes equations when the viscosity depends nonlinearly on the module of the velocity gradient. Two different reliable and efficient residual-based a posteriori error estimators for this problem on arbitrary (convex or non-convex) polygonal and polyhedral regions are derived. Our analysis of reliability of the proposed estimators draws mainly upon the global inf-sup condition satisfied by a suitable linearization of the continuous formulation, an application of Helmholtz decomposition, and the local approximation properties of the Raviart-Thomas and Clément interpolation operators. In addition, differently from previous approaches for augmented mixed formulations, the boundedness of the Clément operator plays now an interesting role in the reliability estimate. On the other hand, inverse and discrete inequalities, and the localization technique based on triangle-bubble and edge-bubble functions are utilized to show their efficiency. Finally, several numerical results are provided to illustrate the good performance of the augmented mixed method, to confirm the aforementioned properties of the a posteriori error estimators, and to show the behaviour of the associated adaptive algorithm.Centre for Mathematical Modeling (Universidad de Chile)Centro de Investigación en Ingeniería Matemática (Universidad de Concepción)Comisión Nacional de Investigación Científica y Tecnológica (Chile)Swiss National Science FoundationElsevier Mathematical Sciences Sponsorship FundMinistry of Education, Youth and Sports of the Czech Republi
    corecore