A local projection stabilized method for fictitious domains

In this work a local projection stabilization method is proposed to solve a fictitious domain problem. The method adds a suitable fluctuation term to the formulation thus rendering the natural space for the Lagrange multiplier stable. Stability and convergence are proved and these results are illustrated by a numerical experiment.Comment: Submitted Preprin

Stabilised hybrid discontinuous Galerkin methods for the Stokes problem with non-standard boundary conditions

In several studies it has been observed that, when using stabilised $\mathbb{P}_k^{}\times\mathbb{P}_k^{}$ elements for both velocity and pressure, the error for the pressure is smaller, or even of a higher order in some cases, than the one obtained when using inf-sup stable $\mathbb{P}_k^{}\times\mathbb{P}_{k-1}^{}$ (although no formal proof of either of these facts has been given). This increase in polynomial order requires the introduction of stabilising terms, since the finite element pairs used do not stability the inf-sup condition. With this motivation, we apply the stabilisation approach to the hybrid discontinuous Galerkin discretisation for the Stokes problem with non-standard boundary conditions

Finite element eigenvalue enclosures for the Maxwell operator

We propose employing the extension of the Lehmann-Maehly-Goerisch method developed by Zimmermann and Mertins, as a highly effective tool for the pollution-free finite element computation of the eigenfrequencies of the resonant cavity problem on a bounded region. This method gives complementary bounds for the eigenfrequencies which are adjacent to a given real parameter. We present a concrete numerical scheme which provides certified enclosures in a suitable asymptotic regime. We illustrate the applicability of this scheme by means of some numerical experiments on benchmark data using Lagrange elements and unstructured meshes.Comment: arXiv admin note: substantial text overlap with arXiv:1306.535

Stabilised finite element methods for the Oseen problem on anisotropic quadrilateral meshes

In this work we present and analyse new inf-sup stable, and stabilised, finite element methods for the Oseen equation in anisotropic quadrilateral meshes. The meshes are formed of closed parallelograms, and the analysis is restricted to two space dimensions. Starting with the lowest order QIn this work we present and analyse new inf-sup stable, and stabilised, finite element methods for the Oseen equation in anisotropic quadrilateral meshes. The meshes are formed of closed parallelograms, and the analysis is restricted to two space dimensions. Starting with the lowest order Q2 1 × P0 pair, we first identify the pressure components that make this finite element pair to be non-inf-sup stable, especially with respect to the aspect ratio. We then propose a way to penalise them, both strongly, by directly removing them from the space, and weakly, by adding a stabilisation term based on jumps of the pressure across selected edges. Concerning the velocity stabilisation, we propose an enhanced grad-div term. Stability and optimal a priori error estimates are given, and the results are confirmed numerically. Q21 × P0 pair, we first identify the pressure components that make this finite element pair to be non-inf-sup stable, especially with respect to the aspect ratio. We then propose a way to penalise them, both strongly, by directly removing them from the space, and weakly, by adding a stabilisation term based on jumps of the pressure across selected edges. Concerning the velocity stabilisation, we propose an enhanced grad-div term. Stability and optimal a priori error estimates are given, and the results are confirmed numerically

A symmetric nodal conservative finite element method for the Darcy equation

This work introduces and analyzes novel stable Petrov-Galerkin EnrichedMethods (PGEM) for the Darcy problem based on the simplest but unstable continuous P1/P0 pair. Stability is recovered inside a Petrov-Galerkin framework where element-wise dependent residual functions, named multi-scale functions, enrich both velocity and pressure trial spaces. Unlike the velocity test space that is augmented with bubble-like functions, multi-scale functions correct edge residuals as well. The multi-scale functions turn out to be the well-known lowest order Raviart-Thomas basis functions for the velocity and discontinuous quadratics polynomial functions for the pressure. The enrichment strategy suggests the way to recover the local mass conservation property for nodal-based interpolation spaces. We prove that the method and its symmetric version are well-posed and achieve optimal error estimates in natural norms. Numerical validations confirm claimed theoretical results

Analysis of algebraic flux correction schemes

A family of algebraic flux correction schemes for linear boundary value problems in any space dimension is studied. These methods’ main feature is that they limit the fluxes along each one of the edges of the triangulation, and we suppose that the limiters used are symmetric. For an abstract problem, the existence of a solution, existence and uniqueness of the solution of a linearized problem, and an a priori error estimate, are proved under rather general assumptions on the limiters. For a particular (but standard in practice) choice of the limiters, it is shown that a local discrete maximum principle holds. The theory developed for the abstract problem is applied to convection–diffusion–reaction equations, where in particular an error estimate is derived. Numerical studies show its sharpness

Some analytical results for an algebraic flux correction scheme for a steady convection-diffusion equation in 1D

Algebraic flux correction schemes are nonlinear discretizations of convection dominated problems. In this work, a scheme from this class is studied for a steady-state convection--diffusion equation in one dimension. It is proved that this scheme satisfies the discrete maximum principle. Also, as it is a nonlinear scheme, the solvability of the linear subproblems arising in a Picard iteration is studied, where positive and negative results are proved. Furthermore, the non-existence of solutions for the nonlinear scheme is proved by means of counterexamples. Therefore, a modification of the method, which ensures the existence of a solution, is proposed. A weak version of the discrete maximum principle is proved for this modified method

Stabilisation of high aspect ratio mixed finite elements for incompressible flow

Anisotropically refined mixed finite elements are beneficial for the resolution of local features such as boundary layers. Unfortunately, the stability of the resulting scheme is highly sensitive to the aspect ratio of the elements. Previous analysis revealed that the degeneration arises from a relatively small number of spurious (piecewise constant) pressure modes. The present article is concerned with resolving the problem of how to suppress the spurious pressure modes in order to restore stability yet at the same time not incur any deterioration in the approximation properties of the reduced pressure space. Two results are presented. The first gives the minimal constraints on the pressure space needed to restore stability with respect to aspect ratio and shows that the approximation properties of the constrained pressure space and the unconstrained pressure space are essentially identical. Alternatively, one can impose the constraint weakly through the use of a stabilized finite element scheme. A second result shows that the stabilized finite element scheme is robust with respect to the aspect ratio of the elements and produces an approximation that satisfies an error bound of the same type to the mixed finite element scheme using the constrained space

Analysis of algebraic flux correction schemes

A family of algebraic flux correction schemes for linear boundary value problems in any space dimension is studied. These methods' main feature is that they limit the fluxes along each one of the edges of the triangulation, and we suppose that the limiters used are symmetric. For an abstract problem, the existence of a solution, existence and uniqueness of the solution of a linearized problem, and an a priori error estimate, are proved under rather general assumptions on the limiters. For a particular (but standard in practice) choice of the limiters, it is shown that a local discrete maximum principle holds. The theory developed for the abstract problem is applied to convection-diffusion-reaction equations, where in particular an error estimate is derived. Numerical studies show its sharpness