English

A Hellan-Herrmann-Johnson-like method for the stream function formulation of the Stokes equations in two and three space dimensions

Numerical Analysis 2020-05-19 v2 Numerical Analysis

Abstract

We introduce a new discretization for the stream function formulation of the incompressible Stokes equations in two and three space dimensions. The method is strongly related to the Hellan-Herrmann-Johnson method and is based on the recently discovered mass conserving mixed stress formulation [J. Gopalakrishnan, P.L. Lederer, J. Sch\"oberl, IMA Journal of numerical Analysis, 2019] that approximates the velocity in an H(div)H(\operatorname{div})-conforming space and introduces a new stress-like variable for the approximation of the gradient of the velocity within the function space H(curldiv)H(\operatorname{curl}\operatorname{div}). The properties of the (discrete) de Rham complex allows to extend this method to a stream function formulation in two and three space dimensions. We present a detailed stability analysis in the continuous and the discrete setting where the stream function ψ\psi and its approximation ψh\psi_h are elements of H(curl)H(\operatorname{curl}) and the H(curl)H(\operatorname{curl})-conforming N\'ed\'elec finite element space, respectively. We conclude with an error analysis revealing optimal convergence rates for the error of the discrete velocity uh=curl(ψh)u_h = \operatorname{curl}(\psi_h) measured in a discrete H1H^1-norm. We present numerical examples to validate our findings and discuss structure-preserving properties such as pressure-robustness.

Keywords

Cite

@article{arxiv.2005.06506,
  title  = {A Hellan-Herrmann-Johnson-like method for the stream function formulation of the Stokes equations in two and three space dimensions},
  author = {Philip L. Lederer},
  journal= {arXiv preprint arXiv:2005.06506},
  year   = {2020}
}
R2 v1 2026-06-23T15:31:29.376Z