Divergence boundary conditions for vector Helmholtz equations with divergence constraints
The idea of replacing a divergence constraint by a divergence boundary condition is investigated. The connections between the formulations are considered in detail. It is shown that the most common methods of using divergence boundary conditions do not always work properly. Necessary and sufficient conditions for the equivalence of the formulations are given.
By re-examining the arguments and counterexamples in I. Babuška, A. K. Aziz (1976) concerning the well-known maximum angle condition, we study the convergence behavior of the linear finite element method (FEM) on a family of distorted triangulations of the unit square originally introduced by H. Schwarz in 1880. For a Poisson problem with polynomial solution, we demonstrate arbitrarily slow convergence as well as failure of convergence if the distortion of the triangulations grows sufficiently fast....
An axisymmetric second order elliptic problem with mixed boundary conditions is considered. The shape of the domain has to be found so as to minimize a cost functional, which is given in terms of the cogradient of the solution. A new dual finite element method is used for approximate solutions. The existence of an optimal domain is proven and a convergence analysis presented.
The central objective of this paper is to develop reduced basis methods for parameter dependent transport dominated problems that are rigorously proven to exhibit rate-optimal performance when compared with the Kolmogorov n-widths of the solution sets. The central ingredient is the construction of computationally feasible “tight” surrogates which in turn are based on deriving a suitable well-conditioned variational formulation for the parameter dependent problem. The theoretical results are illustrated...
A mixed finite element method for the Navier–Stokes equations is introduced in which the stress is a primary variable. The variational formulation retains the mathematical structure of the Navier–Stokes equations and the classical theory extends naturally to this setting. Finite element spaces satisfying the associated inf–sup conditions are developed.
The Maxwell type elastic-plastic solids are characterized by decaying the absolute values of the principal components of the deviatoric part of the stress tensor during the plastic relaxation step. We propose a mathematical formulation of such a model which is compatible with the von Mises criterion of plasticity. Numerical examples show the ability of the model to deal with complex physical phenomena.
We consider the solution of second order elliptic PDEs in Rd with inhomogeneous Dirichlet data by means of an h–adaptive FEM with fixed polynomial order p ∈ N. As model example serves the Poisson equation with mixed Dirichlet–Neumann boundary conditions, where the inhomogeneous Dirichlet data are discretized by use of an H1 / 2–stable projection, for instance, the L2–projection for p = 1 or the Scott–Zhang projection for general p ≥ 1. For error estimation, we use a residual error estimator which...
In this paper we consider the Maxwell resolvent operator and its finite element approximation. In this framework it is natural the use of the edge element spaces and to impose the divergence constraint in a weak sense with the introduction of a Lagrange multiplier, following an idea by Kikuchi [14]. We shall review some of the known properties for edge element approximations and prove some new result. In particular we shall prove a uniform convergence in the norm for the sequence of discrete operators....
In this paper we consider the Maxwell resolvent operator and its finite element approximation. In this framework it is natural the use of the edge element spaces and to impose the divergence constraint in a weak sense with the introduction of a Lagrange multiplier, following an idea by Kikuchi [14]. We shall review some of the known properties for edge element approximations and prove some new result. In particular we shall prove a uniform convergence in the L2 norm for the sequence of discrete...
We consider mixed finite element discretizations of second order elliptic boundary value problems. Emphasis is on the efficient iterative solution by multilevel techniques with respect to an adaptively generated hierarchy of nonuniform triangulations. In particular, we present two multilevel solvers, the first one relying on ideas from domain decomposition and the second one resulting from mixed hybridization. Local refinement of the underlying triangulations is done by efficient and reliable a...
The electrowetting process is commonly used to handle very small amounts of liquid on a solid surface. This process can be modelled mathematically with the help of the shape optimization theory. However, solving numerically the resulting shape optimization problem is a very complex issue, even for reduced models that occur in simplified geometries. Recently, the second author obtained convincing results in the 2D axisymmetric case. In this paper, we propose and analyze a method that is suitable...
Les méthodes sans maillage emploient une interpolation associée à un ensemble de particules : aucune information concernant la connectivité ne doit être fournie. Un des atouts de ces méthodes est que la discrétisation peut être enrichie d’une façon très simple, soit en augmentant le nombre de particules (analogue à la stratégie de raffinement ), soit en augmentant l’ordre de consistance (analogue à la stratégie de raffinement ). Néanmoins, le coût du calcul des fonctions d’interpolation est très...
Les méthodes sans maillage emploient une interpolation associée à un ensemble de particules : aucune information concernant la connectivité ne doit être fournie. Un des atouts de ces méthodes est que la discrétisation peut être enrichie d'une façon très simple, soit en augmentant le nombre de particules (analogue à la stratégie de raffinement h), soit en augmentant l'ordre de consistance (analogue à la stratégie de raffinement p). Néanmoins, le coût du calcul des fonctions d'interpolation est...
For convection-diffusion problems with exponential layers, optimal error estimates for linear finite elements on Shishkin-type meshes are known. We present the first optimal convergence result in an energy norm for a Bakhvalov-type mesh.