A divergence-free low-order stabilized finite element method for a generalized steady state Boussinesq problem

Alejandro Allendes, Gabriel R. Barrenechea, Cesar Narranjo

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

In this work we propose and analyze a new stabilized finite element method for the coupled Navier-Stokes/temperature (or Boussinesq) equations. The method is built using low order conforming elements for velocity and temperature, and piecewise constant elements for pressure. With the help of the lowest order Raviart-Thomas space, a lifting of the jumps of the discrete pressure is built in such a way that when this lifting is added to the conforming velocity field, the resulting velocity is solenoidal (at the price of being non-conforming). This field is then fed to the momentum and temperature equations, guaranteeing that the convective terms in these equations are antisymmetric, without the need of altering them, thus simplifying the analysis of the resulting method. Existence of solutions, discrete stability, and optimal convergence are proved for both the conforming velocity field, and its corresponding divergence-free non-conforming counterpart. Numerical results confirm the theoretical findings, as well as the gain provided by the solenoidal discrete velocity field over the conforming one.
LanguageEnglish
Pages90-120
Number of pages31
JournalComputer Methods in Applied Mechanics and Engineering
Volume340
Early online date31 May 2018
DOIs
Publication statusPublished - 1 Oct 2018

Fingerprint

divergence
finite element method
velocity distribution
Finite element method
temperature
momentum
Temperature
Momentum

Keywords

  • Boussinesq problem
  • stabilized finite element method
  • divergence-free discrete velocity
  • Navier–Stokes/temperature equations

Cite this

@article{db3dbca880b843a6ab69bc143695903c,
title = "A divergence-free low-order stabilized finite element method for a generalized steady state Boussinesq problem",
abstract = "In this work we propose and analyze a new stabilized finite element method for the coupled Navier-Stokes/temperature (or Boussinesq) equations. The method is built using low order conforming elements for velocity and temperature, and piecewise constant elements for pressure. With the help of the lowest order Raviart-Thomas space, a lifting of the jumps of the discrete pressure is built in such a way that when this lifting is added to the conforming velocity field, the resulting velocity is solenoidal (at the price of being non-conforming). This field is then fed to the momentum and temperature equations, guaranteeing that the convective terms in these equations are antisymmetric, without the need of altering them, thus simplifying the analysis of the resulting method. Existence of solutions, discrete stability, and optimal convergence are proved for both the conforming velocity field, and its corresponding divergence-free non-conforming counterpart. Numerical results confirm the theoretical findings, as well as the gain provided by the solenoidal discrete velocity field over the conforming one.",
keywords = "Boussinesq problem, stabilized finite element method, divergence-free discrete velocity, Navier–Stokes/temperature equations",
author = "Alejandro Allendes and Barrenechea, {Gabriel R.} and Cesar Narranjo",
year = "2018",
month = "10",
day = "1",
doi = "10.1016/j.cma.2018.05.020",
language = "English",
volume = "340",
pages = "90--120",
journal = "Computer Methods in Applied Mechanics end Engineering",
issn = "0045-7825",

}

A divergence-free low-order stabilized finite element method for a generalized steady state Boussinesq problem. / Allendes, Alejandro; Barrenechea, Gabriel R.; Narranjo, Cesar.

In: Computer Methods in Applied Mechanics and Engineering, Vol. 340, 01.10.2018, p. 90-120.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A divergence-free low-order stabilized finite element method for a generalized steady state Boussinesq problem

AU - Allendes, Alejandro

AU - Barrenechea, Gabriel R.

AU - Narranjo, Cesar

PY - 2018/10/1

Y1 - 2018/10/1

N2 - In this work we propose and analyze a new stabilized finite element method for the coupled Navier-Stokes/temperature (or Boussinesq) equations. The method is built using low order conforming elements for velocity and temperature, and piecewise constant elements for pressure. With the help of the lowest order Raviart-Thomas space, a lifting of the jumps of the discrete pressure is built in such a way that when this lifting is added to the conforming velocity field, the resulting velocity is solenoidal (at the price of being non-conforming). This field is then fed to the momentum and temperature equations, guaranteeing that the convective terms in these equations are antisymmetric, without the need of altering them, thus simplifying the analysis of the resulting method. Existence of solutions, discrete stability, and optimal convergence are proved for both the conforming velocity field, and its corresponding divergence-free non-conforming counterpart. Numerical results confirm the theoretical findings, as well as the gain provided by the solenoidal discrete velocity field over the conforming one.

AB - In this work we propose and analyze a new stabilized finite element method for the coupled Navier-Stokes/temperature (or Boussinesq) equations. The method is built using low order conforming elements for velocity and temperature, and piecewise constant elements for pressure. With the help of the lowest order Raviart-Thomas space, a lifting of the jumps of the discrete pressure is built in such a way that when this lifting is added to the conforming velocity field, the resulting velocity is solenoidal (at the price of being non-conforming). This field is then fed to the momentum and temperature equations, guaranteeing that the convective terms in these equations are antisymmetric, without the need of altering them, thus simplifying the analysis of the resulting method. Existence of solutions, discrete stability, and optimal convergence are proved for both the conforming velocity field, and its corresponding divergence-free non-conforming counterpart. Numerical results confirm the theoretical findings, as well as the gain provided by the solenoidal discrete velocity field over the conforming one.

KW - Boussinesq problem

KW - stabilized finite element method

KW - divergence-free discrete velocity

KW - Navier–Stokes/temperature equations

UR - https://www.journals.elsevier.com/computer-methods-in-applied-mechanics-and-engineering

U2 - 10.1016/j.cma.2018.05.020

DO - 10.1016/j.cma.2018.05.020

M3 - Article

VL - 340

SP - 90

EP - 120

JO - Computer Methods in Applied Mechanics end Engineering

T2 - Computer Methods in Applied Mechanics end Engineering

JF - Computer Methods in Applied Mechanics end Engineering

SN - 0045-7825

ER -