A variational H (div) finite-element discretization approach for perfect incompressible fluids
File(s)Euler_paper.pdf (1.85 MB)
Accepted version
Author(s)
Natale, A
Cotter, CJ
Type
Journal Article
Abstract
We propose a finite-element discretization approach for the incompressible Euler equations which mimics
their geometric structure and their variational derivation. In particular, we derive a finite-element method
that arises from a nonholonomic variational principle and an appropriately defined Lagrangian, where finite-
element
H
(
div
)
vector fields are identified with advection operators; this is the first successful extension
of the structure-preserving discretization of
Pavlov
et al.
(
2009
) to the finite-element setting. The resulting
algorithm coincides with the energy-conserving scheme proposed by
Guzm
́
an
et al.
(
2016
). Through the
variational derivation, we discover that it also satisfies a discrete analogous of Kelvin’s circulation theorem.
Further, we propose an upwind-stabilized version of the scheme that dissipates enstrophy while preserving
energy conservation and the discrete Kelvin’s theorem. We prove error estimates for this version of the
scheme, and we study its behaviour through numerical tests.
their geometric structure and their variational derivation. In particular, we derive a finite-element method
that arises from a nonholonomic variational principle and an appropriately defined Lagrangian, where finite-
element
H
(
div
)
vector fields are identified with advection operators; this is the first successful extension
of the structure-preserving discretization of
Pavlov
et al.
(
2009
) to the finite-element setting. The resulting
algorithm coincides with the energy-conserving scheme proposed by
Guzm
́
an
et al.
(
2016
). Through the
variational derivation, we discover that it also satisfies a discrete analogous of Kelvin’s circulation theorem.
Further, we propose an upwind-stabilized version of the scheme that dissipates enstrophy while preserving
energy conservation and the discrete Kelvin’s theorem. We prove error estimates for this version of the
scheme, and we study its behaviour through numerical tests.
Date Issued
2017-07-01
Date Acceptance
2017-05-16
Citation
IMA Journal of Numerical Analysis, 2017, 38 (3), pp.1388-1419
ISSN
0272-4979
Publisher
Oxford University Press
Start Page
1388
End Page
1419
Journal / Book Title
IMA Journal of Numerical Analysis
Volume
38
Issue
3
Copyright Statement
© The authors 2017. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.
Sponsor
Engineering & Physical Science Research Council (EPSRC)
Natural Environment Research Council (NERC)
Natural Environment Research Council (NERC)
Natural Environment Research Council (NERC)
Natural Environment Research Council (NERC)
Natural Environment Research Council (NERC)
Grant Number
EP/L000407/1
NE/I016007/1
NE/I000747/1
NE/I02013X/1
NE/K006789/1
NE/M013634/1
Subjects
Science & Technology
Physical Sciences
Mathematics, Applied
Mathematics
Euler equations
perfect incompressible fluids
finite-element methods
structure-preserving methods
STABILIZED GALERKIN
EULER-EQUATIONS
ADVECTION
LIE
math.NA
math.NA
math-ph
math.MP
0102 Applied Mathematics
0103 Numerical and Computational Mathematics
Numerical & Computational Mathematics
Publication Status
Published
Date Publish Online
2017-06-29