We present a modification of the Crouzeix-Raviart discretization of the Stokes equations in arbitrary dimension which is quasi-optimal, in the sense that the error of the discrete velocity field in a broken H1-norm is proportional to the error of the best approximation to the analytical velocity field. In particular, the velocity error is independent of the pressure error and the discrete velocity field is elementwise solenoidal. Moreover, the sum of the velocity error times the viscosity plus the pressure L2-error is proportional to the sum of the respective best errors. All proportionality constants are bounded in terms of shape regularity and do not depend on the viscosity. For simply connected two-dimensional domains, the velocity and pressure can be computed separately. The modification only affects the right-hand side, aka the load vector. The cost for building the modified load vector is proportional to the cost for building the standard load vector. Some numerical experiments illustrate our theoretical results.

A quasi-optimal Crouzeix-Raviart discretization of the Stokes equations / R. Verfurth, P. Zanotti. - In: SIAM JOURNAL ON NUMERICAL ANALYSIS. - ISSN 0036-1429. - 57:3(2019), pp. 1082-1099. [10.1137/18M1177688]

A quasi-optimal Crouzeix-Raviart discretization of the Stokes equations

P. Zanotti
2019

Abstract

We present a modification of the Crouzeix-Raviart discretization of the Stokes equations in arbitrary dimension which is quasi-optimal, in the sense that the error of the discrete velocity field in a broken H1-norm is proportional to the error of the best approximation to the analytical velocity field. In particular, the velocity error is independent of the pressure error and the discrete velocity field is elementwise solenoidal. Moreover, the sum of the velocity error times the viscosity plus the pressure L2-error is proportional to the sum of the respective best errors. All proportionality constants are bounded in terms of shape regularity and do not depend on the viscosity. For simply connected two-dimensional domains, the velocity and pressure can be computed separately. The modification only affects the right-hand side, aka the load vector. The cost for building the modified load vector is proportional to the cost for building the standard load vector. Some numerical experiments illustrate our theoretical results.
Crouzeix-Raviart element; Nonconforming finite elements; Pressure-robustness; Quasi-optimality; Stokes equations
Settore MAT/08 - Analisi Numerica
Article (author)
File in questo prodotto:
Non ci sono file associati a questo prodotto.
Pubblicazioni consigliate

Caricamento pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/2434/722269
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 6
  • ???jsp.display-item.citation.isi??? 7
social impact