Stabilization of low-order cross-grid PkQl mixed finite elements
In this paper we analyze a low-order family of mixed finite element methods for the numerical solution of the Stokes problem and a second order elliptic problem, in two space dimensions. In these schemes, the pressure is interpolated on a mesh of rectangular elements, while the velocity is approxima...
Guardado en:
Publicado: |
2018
|
---|---|
Materias: | |
Acceso en línea: | https://bibliotecadigital.exactas.uba.ar/collection/paper/document/paper_03770427_v330_n_p340_Armentano http://hdl.handle.net/20.500.12110/paper_03770427_v330_n_p340_Armentano |
Aporte de: |
id |
paper:paper_03770427_v330_n_p340_Armentano |
---|---|
record_format |
dspace |
spelling |
paper:paper_03770427_v330_n_p340_Armentano2023-06-08T15:38:58Z Stabilization of low-order cross-grid PkQl mixed finite elements Cross-grid elements Elliptic problems Mixed finite elements Stability analysis Stokes problem Computational mechanics Mesh generation Navier Stokes equations Numerical methods Stabilization Elliptic problem Grid elements Mixed finite elements Stability analysis Stokes problem Finite element method In this paper we analyze a low-order family of mixed finite element methods for the numerical solution of the Stokes problem and a second order elliptic problem, in two space dimensions. In these schemes, the pressure is interpolated on a mesh of rectangular elements, while the velocity is approximated on a triangular mesh obtained by dividing each rectangle into four triangles by its diagonals. For the lowest order P1Q0, a global spurious pressure mode is shown to exist and so this element, as P1Q1 case analyzed in Armentano and Blasco (2010), is unstable. However, following the ideas given in Bochev et al. (2006), a simple stabilization procedure can be applied, when we approximate the solution of the Stokes problem, such that the new P1Q0 and P1Q1 methods are unconditionally stable, and achieve optimal accuracy with respect to solution regularity with simple and straightforward implementations. Moreover, we analyze the application of our P1Q1 element to the mixed formulation of the elliptic problem. In this case, by introducing the modified mixed weak form proposed in Brezzi et al. (1993), optimal order of accuracy can be obtained with our stabilized P1Q1 elements. Numerical results are also presented, which confirm the existence of the spurious pressure mode for the P1Q0 element and the excellent stability and accuracy of the new stabilized methods. © 2017 Elsevier B.V. 2018 https://bibliotecadigital.exactas.uba.ar/collection/paper/document/paper_03770427_v330_n_p340_Armentano http://hdl.handle.net/20.500.12110/paper_03770427_v330_n_p340_Armentano |
institution |
Universidad de Buenos Aires |
institution_str |
I-28 |
repository_str |
R-134 |
collection |
Biblioteca Digital - Facultad de Ciencias Exactas y Naturales (UBA) |
topic |
Cross-grid elements Elliptic problems Mixed finite elements Stability analysis Stokes problem Computational mechanics Mesh generation Navier Stokes equations Numerical methods Stabilization Elliptic problem Grid elements Mixed finite elements Stability analysis Stokes problem Finite element method |
spellingShingle |
Cross-grid elements Elliptic problems Mixed finite elements Stability analysis Stokes problem Computational mechanics Mesh generation Navier Stokes equations Numerical methods Stabilization Elliptic problem Grid elements Mixed finite elements Stability analysis Stokes problem Finite element method Stabilization of low-order cross-grid PkQl mixed finite elements |
topic_facet |
Cross-grid elements Elliptic problems Mixed finite elements Stability analysis Stokes problem Computational mechanics Mesh generation Navier Stokes equations Numerical methods Stabilization Elliptic problem Grid elements Mixed finite elements Stability analysis Stokes problem Finite element method |
description |
In this paper we analyze a low-order family of mixed finite element methods for the numerical solution of the Stokes problem and a second order elliptic problem, in two space dimensions. In these schemes, the pressure is interpolated on a mesh of rectangular elements, while the velocity is approximated on a triangular mesh obtained by dividing each rectangle into four triangles by its diagonals. For the lowest order P1Q0, a global spurious pressure mode is shown to exist and so this element, as P1Q1 case analyzed in Armentano and Blasco (2010), is unstable. However, following the ideas given in Bochev et al. (2006), a simple stabilization procedure can be applied, when we approximate the solution of the Stokes problem, such that the new P1Q0 and P1Q1 methods are unconditionally stable, and achieve optimal accuracy with respect to solution regularity with simple and straightforward implementations. Moreover, we analyze the application of our P1Q1 element to the mixed formulation of the elliptic problem. In this case, by introducing the modified mixed weak form proposed in Brezzi et al. (1993), optimal order of accuracy can be obtained with our stabilized P1Q1 elements. Numerical results are also presented, which confirm the existence of the spurious pressure mode for the P1Q0 element and the excellent stability and accuracy of the new stabilized methods. © 2017 Elsevier B.V. |
title |
Stabilization of low-order cross-grid PkQl mixed finite elements |
title_short |
Stabilization of low-order cross-grid PkQl mixed finite elements |
title_full |
Stabilization of low-order cross-grid PkQl mixed finite elements |
title_fullStr |
Stabilization of low-order cross-grid PkQl mixed finite elements |
title_full_unstemmed |
Stabilization of low-order cross-grid PkQl mixed finite elements |
title_sort |
stabilization of low-order cross-grid pkql mixed finite elements |
publishDate |
2018 |
url |
https://bibliotecadigital.exactas.uba.ar/collection/paper/document/paper_03770427_v330_n_p340_Armentano http://hdl.handle.net/20.500.12110/paper_03770427_v330_n_p340_Armentano |
_version_ |
1768542890770825216 |