A MIXED FINITE ELEMENT METHOD FOR DARCY'S EQUATIONS WITH PRESSURE DEPENDENT POROSITY

Kód výsledku v IS VaVaI

<a href="https://www.isvavai.cz/riv?ss=detail&h=RIV%2F00216208%3A11320%2F16%3A10334325" target="_blank" >RIV/00216208:11320/16:10334325 - isvavai.cz</a>

Výsledek na webu

<a href="http://dx.doi.org/10.1090/mcom/2980" target="_blank" >http://dx.doi.org/10.1090/mcom/2980</a>

DOI - Digital Object Identifier

<a href="http://dx.doi.org/10.1090/mcom/2980" target="_blank" >10.1090/mcom/2980</a>

Jazyk výsledku

angličtina

Název v původním jazyce

A MIXED FINITE ELEMENT METHOD FOR DARCY'S EQUATIONS WITH PRESSURE DEPENDENT POROSITY

Popis výsledku v původním jazyce

In this work we develop the a priori and a posteriori error analyses of a mixed finite element method for Darcy's equations with porosity depending exponentially on the pressure. A simple change of variable for this unknown allows us to transform the original nonlinear problem into a linear one whose dual-mixed variational formulation falls into the frameworks of the generalized linear saddle point problems and the fixed point equations satisfied by an affine mapping. According to the latter, we are able to show the well-posedness of both the continuous and discrete schemes, as well as the associated Cea estimate, by simply applying a suitable combination of the classical Babuska-Brezzi theory and the Banach fixed point theorem. In particular, given any integer k >= 0, the stability of the Galerkin scheme is guaranteed by employing Raviart-Thomas elements of order k for the velocity, piecewise polynomials of degree k for the pressure, and continuous piecewise polynomials of degree k + 1 for an additional Lagrange multiplier given by the trace of the pressure on the Neumann boundary. Note that the two ways of writing the continuous formulation suggest accordingly two different methods for solving the discrete schemes. Next, we derive a reliable and efficient residual-based a posteriori error estimator for this problem. The global inf-sup condition satisfied by the continuous formulation, Helmholtz decompositions, and the local approximation properties of the Raviart-Thomas and Clement interpolation operators are the main tools for proving the reliability. In turn, inverse and discrete 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 both methods, confirming the aforementioned properties of the estimator, and showing the behaviour of the associated adaptive algorithm, are reported.

Název v anglickém jazyce

A MIXED FINITE ELEMENT METHOD FOR DARCY'S EQUATIONS WITH PRESSURE DEPENDENT POROSITY

Popis výsledku anglicky

In this work we develop the a priori and a posteriori error analyses of a mixed finite element method for Darcy's equations with porosity depending exponentially on the pressure. A simple change of variable for this unknown allows us to transform the original nonlinear problem into a linear one whose dual-mixed variational formulation falls into the frameworks of the generalized linear saddle point problems and the fixed point equations satisfied by an affine mapping. According to the latter, we are able to show the well-posedness of both the continuous and discrete schemes, as well as the associated Cea estimate, by simply applying a suitable combination of the classical Babuska-Brezzi theory and the Banach fixed point theorem. In particular, given any integer k >= 0, the stability of the Galerkin scheme is guaranteed by employing Raviart-Thomas elements of order k for the velocity, piecewise polynomials of degree k for the pressure, and continuous piecewise polynomials of degree k + 1 for an additional Lagrange multiplier given by the trace of the pressure on the Neumann boundary. Note that the two ways of writing the continuous formulation suggest accordingly two different methods for solving the discrete schemes. Next, we derive a reliable and efficient residual-based a posteriori error estimator for this problem. The global inf-sup condition satisfied by the continuous formulation, Helmholtz decompositions, and the local approximation properties of the Raviart-Thomas and Clement interpolation operators are the main tools for proving the reliability. In turn, inverse and discrete 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 both methods, confirming the aforementioned properties of the estimator, and showing the behaviour of the associated adaptive algorithm, are reported.

Druh

J_x - Nezařazeno - Článek v odborném periodiku (Jimp, Jsc a Jost)

CEP obor

BA - Obecná matematika

OECD FORD obor

—

Projekt

<a href="/cs/project/LL1202" target="_blank" >LL1202: Materiály s implicitními konstitutivními vztahy: Od teorie přes redukci modelů k efektivním numerickým metodám</a><br>

Návaznosti

P - Projekt vyzkumu a vyvoje financovany z verejnych zdroju (s odkazem do CEP)<br>I - Institucionalni podpora na dlouhodoby koncepcni rozvoj vyzkumne organizace

Rok uplatnění

2016

Kód důvěrnosti údajů

S - Úplné a pravdivé údaje o projektu nepodléhají ochraně podle zvláštních právních předpisů

Název periodika

Mathematics of Computation

ISSN

0025-5718

e-ISSN

—

Svazek periodika

85

Číslo periodika v rámci svazku

297

Stát vydavatele periodika

US - Spojené státy americké

Počet stran výsledku

33

Strana od-do

1-33

Kód UT WoS článku

000362848100001

EID výsledku v databázi Scopus

2-s2.0-85000347802