An augmented velocity–vorticity–pressure formulation for the Brinkman equations. (12th May 2015)

• Record Type: Journal Article
• Authors: Anaya, Verónica; Gatica, Gabriel N.; Mora, David; Ruiz‐Baier, Ricardo
• Is Part Of: International journal for numerical methods in fluids. Volume 79:Number 3(2015)
• Page Start: 109
• Page End: 137
• Publication Date: 2015-05-12
• Subjects: Fluid dynamics -- Mathematics -- Periodicals; 532
• DOI: 10.1002/fld.4041
• Languages: English
• ISSNs: 0271-2091

Abstract

<abstract abstract-type="main" id="fld4041-abs-0001"> <title>Summary</title> <p id="fld4041-para-0001">This paper deals with the analysis of a new augmented mixed finite element method in terms of vorticity, velocity, and pressure, for the Brinkman problem with nonstandard boundary conditions. The approach is based on the introduction of Galerkin least‐squares terms arising from the constitutive equation relating the aforementioned unknowns and from the incompressibility condition. We show that the resulting augmented bilinear form is continuous and elliptic, which, thanks to the Lax–Milgram theorem, and besides proving the well‐posedness of the continuous formulation, ensures the solvability and stability of the Galerkin scheme with any finite element subspace of the continuous space. In particular, Raviart–Thomas elements of any order <inline-formula><alternatives><inline-graphic mimetype="image" xlink:href="ark:/27927/pgj2f6xxw52" xlink:type="simple" xmlns:xlink="http://www.w3.org/1999/xlink" /><mml:math display="inline" altimg="urn:x-wiley:fld:media:fld4041:fld4041-math-0001" overflow="scroll" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>k</mml:mi><mml:mo>≥</mml:mo><mml:mn>0</mml:mn></mml:math></alternatives></inline-formula> for the velocity field, and piecewise continuous polynomials of degree <italic>k</italic> + 1 for both the vorticity and the pressure, can be utilized. A priori error estimates and the corresponding rates of convergence are also given here. Next, we derive two reliable and efficient residual‐based a posteriori error estimators for this problem. The ellipticity of the bilinear form together with the local approximation properties of the Clément interpolation operator are the main tools for showing the reliability. In turn, inverse inequalities and the localization technique based on triangle‐bubble and edge‐bubble functions are utilized to show the efficiency. Finally, several numerical results illustrating the good performance of the method, confirming the properties of the estimators and showing the behavior of the associated adaptive algorithms, are reported. Copyright © 2015 John Wiley & Sons, Ltd.</p> </abstract>