Mixed finite element projection methods for the unsteady Stokes equations

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Aricò, Costanza, Helmig, Rainer, Yotov, Ivan
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866914978190065664
author Aricò, Costanza
Helmig, Rainer
Yotov, Ivan
author_facet Aricò, Costanza
Helmig, Rainer
Yotov, Ivan
contents We develop $H$(div)-conforming mixed finite element methods for the unsteady Stokes equations modeling single-phase incompressible fluid flow. A projection method in the framework of the incremental pressure correction methodology is applied, where a predictor and a corrector problems are sequentially solved, accounting for the viscous effects and incompressibility, respectively. The predictor problem is based on a stress-velocity mixed formulation, while the corrector projection problem uses a velocity-pressure mixed formulation. The scheme results in pointwise divergence-free velocity computed at the end of each time step. We establish unconditional stability and first order in time accuracy. In the implementation we focus on generally unstructured triangular grids. We employ a second order multipoint flux mixed finite element method based on the next-to-the-lowest order Raviart-Thomas space $RT_1$ and a suitable quadrature rule. In the predictor problem this approach allows for a local stress elimination, resulting in element-based systems for each velocity component with three degrees of freedom per element. Similarly, in the corrector problem, the velocity is locally eliminated and an element-based system for the pressure is solved. At the end of each time step we obtain a second order accurate $H$(div)-conforming piecewise linear velocity, which is pointwise divergence free. We present a series of numerical tests to illustrate the performance of the method.
format Preprint
id arxiv_https___arxiv_org_abs_2410_14266
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Mixed finite element projection methods for the unsteady Stokes equations
Aricò, Costanza
Helmig, Rainer
Yotov, Ivan
Numerical Analysis
We develop $H$(div)-conforming mixed finite element methods for the unsteady Stokes equations modeling single-phase incompressible fluid flow. A projection method in the framework of the incremental pressure correction methodology is applied, where a predictor and a corrector problems are sequentially solved, accounting for the viscous effects and incompressibility, respectively. The predictor problem is based on a stress-velocity mixed formulation, while the corrector projection problem uses a velocity-pressure mixed formulation. The scheme results in pointwise divergence-free velocity computed at the end of each time step. We establish unconditional stability and first order in time accuracy. In the implementation we focus on generally unstructured triangular grids. We employ a second order multipoint flux mixed finite element method based on the next-to-the-lowest order Raviart-Thomas space $RT_1$ and a suitable quadrature rule. In the predictor problem this approach allows for a local stress elimination, resulting in element-based systems for each velocity component with three degrees of freedom per element. Similarly, in the corrector problem, the velocity is locally eliminated and an element-based system for the pressure is solved. At the end of each time step we obtain a second order accurate $H$(div)-conforming piecewise linear velocity, which is pointwise divergence free. We present a series of numerical tests to illustrate the performance of the method.
title Mixed finite element projection methods for the unsteady Stokes equations
topic Numerical Analysis
url https://arxiv.org/abs/2410.14266