Fully consistent lowest-order finite element methods for generalised Stokes flows with variable viscosity

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Galarce, Felipe, Pacheco, Douglas R. Q.
Format: Preprint
Veröffentlicht: 2024
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866917828048715776
author Galarce, Felipe
Pacheco, Douglas R. Q.
author_facet Galarce, Felipe
Pacheco, Douglas R. Q.
contents Variable viscosity arises in many flow scenarios, often imposing numerical challenges. Yet, discretisation methods designed specifically for non-constant viscosity are few, and their analysis is even scarcer. In finite element methods for incompressible flows, the most popular approach to allow equal-order velocity-pressure interpolation are residual-based stabilisations. For low-order elements, however, the viscous part of that residual cannot be approximated, often compromising accuracy. Assuming slightly more regularity on the viscosity field, we can construct stabilisation methods that fully approximate the residual, regardless of the polynomial order of the finite element spaces. This work analyses two variants of this fully consistent approach, with the generalised Stokes system as a model problem. We prove unique solvability and derive expressions for the stabilisation parameter, generalising some classical results for constant viscosity. Numerical results illustrate how our method completely eliminates the spurious pressure boundary layers typically induced by low-order PSPG-like stabilisations.
format Preprint
id arxiv_https___arxiv_org_abs_2411_02100
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Fully consistent lowest-order finite element methods for generalised Stokes flows with variable viscosity
Galarce, Felipe
Pacheco, Douglas R. Q.
Numerical Analysis
Variable viscosity arises in many flow scenarios, often imposing numerical challenges. Yet, discretisation methods designed specifically for non-constant viscosity are few, and their analysis is even scarcer. In finite element methods for incompressible flows, the most popular approach to allow equal-order velocity-pressure interpolation are residual-based stabilisations. For low-order elements, however, the viscous part of that residual cannot be approximated, often compromising accuracy. Assuming slightly more regularity on the viscosity field, we can construct stabilisation methods that fully approximate the residual, regardless of the polynomial order of the finite element spaces. This work analyses two variants of this fully consistent approach, with the generalised Stokes system as a model problem. We prove unique solvability and derive expressions for the stabilisation parameter, generalising some classical results for constant viscosity. Numerical results illustrate how our method completely eliminates the spurious pressure boundary layers typically induced by low-order PSPG-like stabilisations.
title Fully consistent lowest-order finite element methods for generalised Stokes flows with variable viscosity
topic Numerical Analysis
url https://arxiv.org/abs/2411.02100