The Galerkin finite element method usually produces spurious oscillations
when solving hyperbolic differential equations. A number of stabilized
methods have been proposed and the area has actively been researched in past
decades to remove these oscillations. Another issue, which arises in the
numerical simulations, is the computational cost.
Resolving phenomena such as shock, discontinuity and rarefaction waves
requires adaptivity of the mesh.
In this work we present an adaptive stabilized finite element method for
compressible Euler equations. The entropy viscosity: a new class of
high-order finite element method is used where the nonlinear viscosity is
based on a local size of an entropy production. We prove a posteriori error
estimation of a quantity of interest in terms of a dual problem for the
linearized Euler equations. The implementation in 2D and 3D as well as
different boundary conditions are discussed.
|