A MIXED FINITE ELEMENT METHOD FOR DARCY'S EQUATIONS WITH PRESSURE DEPENDENT POROSITY
The result's identifiers
Result code in 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>
Result on the web
<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>
Alternative languages
Result language
angličtina
Original language name
A MIXED FINITE ELEMENT METHOD FOR DARCY'S EQUATIONS WITH PRESSURE DEPENDENT POROSITY
Original language description
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.
Czech name
—
Czech description
—
Classification
Type
J<sub>x</sub> - Unclassified - Peer-reviewed scientific article (Jimp, Jsc and Jost)
CEP classification
BA - General mathematics
OECD FORD branch
—
Result continuities
Project
<a href="/en/project/LL1202" target="_blank" >LL1202: Implicitly constituted material models: from theory through model reduction to efficient numerical methods</a><br>
Continuities
P - Projekt vyzkumu a vyvoje financovany z verejnych zdroju (s odkazem do CEP)<br>I - Institucionalni podpora na dlouhodoby koncepcni rozvoj vyzkumne organizace
Others
Publication year
2016
Confidentiality
S - Úplné a pravdivé údaje o projektu nepodléhají ochraně podle zvláštních právních předpisů
Data specific for result type
Name of the periodical
Mathematics of Computation
ISSN
0025-5718
e-ISSN
—
Volume of the periodical
85
Issue of the periodical within the volume
297
Country of publishing house
US - UNITED STATES
Number of pages
33
Pages from-to
1-33
UT code for WoS article
000362848100001
EID of the result in the Scopus database
2-s2.0-85000347802