On a 2D vector Poisson problem with apparently mutually exclusive scalar boundary conditions
This work is devoted to the study of a two-dimensional vector Poisson equation with the normal component of the unknown and the value of the divergence of the unknown prescribed simultaneously on the entire boundary. These two scalar boundary conditions appear prima facie alternative in a standard variational framework. An original variational formulation of this boundary value problem is proposed here. Furthermore, an uncoupled solution algorithm is introduced together with its finite element...
The Collatz method of twosided eigenvalue estimates was extended by K. Rektorys in his monography Variational Methods to the case of differential equations of the form with elliptic operators. This method requires to solve, successively, certain boundary value problems. In the case of partial differential equations, these problems are to be solved approximately, as a rule, and this is the source of further errors. In the work, it is shown how to estimate these additional errors, or how to avoid...
We present and analyse in this paper a novel colocated Finite Volume scheme for the solution of the Stokes problem. It has been developed following two main ideas. On one hand, the discretization of the pressure gradient term is built as the discrete transposed of the velocity divergence term, the latter being evaluated using a natural finite volume approximation; this leads to a non-standard interpolation formula for the expression of the pressure on the edges of the control volumes. On the other...
A simple superconvergent scheme for the derivatives of finite element solution is presented, when linear triangular elements are employed to solve second order elliptic systems with boundary conditions of Newton’s or Neumann’s type. For bounded plane domains with smooth boundary the local -superconvergence of the derivatives in the -norm is proved. The paper is a direct continuations of [2], where an analogous problem with Dirichlet’s boundary conditions is treated.
Second order elliptic systems with boundary conditions of Dirichlet, Neumann’s or Newton’s type are solved by means of linear finite elements on regular uniform triangulations. Error estimates of the optimal order are proved for the averaged gradient on any fixed interior subdomain, provided the problem under consideration is regular in a certain sense.
Second order elliptic systems with Dirichlet boundary conditions are solved by means of affine finite elements on regular uniform triangulations. A simple averagign scheme is proposed, which implies a superconvergence of the gradient. For domains with enough smooth boundary, a global estimate is proved in the -norm. For a class of polygonal domains the global estimate can be proven.
We consider the local projection finite element method for the discretization of a scalar convection-diffusion equation with a divergence-free convection field. We introduce a new fluctuation operator which is defined using an orthogonal projection with respect to a weighted inner product. We prove that the bilinear form corresponding to the discrete problem satisfies an inf-sup condition with respect to the SUPG norm and derive an error estimate for the discrete solution.
Discretization of second order elliptic partial differential equations by discontinuous Galerkin method often results in numerical schemes with penalties. In this paper we analyze these penalized schemes in the context of quite general triangular meshes satisfying only a semiregularity assumption. A new (modified) penalty term is presented and theoretical properties are proven together with illustrative numerical results.
This paper focuses on mathematical modeling and finite element simulation of fluid-structure interaction problems. A simplified problem of two-dimensional incompressible fluid flow interacting with a rigid structure, whose motion is described with one degree of freedom, is considered. The problem is mathematically described and numerically approximated using the finite element method. Two possibilities, namely Taylor-Hood and Scott-Vogelius elements are presented and implemented. Finally, numerical...
We propose an analogue of the maximum angle condition (commonly used in finite element analysis for triangular and tetrahedral meshes) for the case of prismatic elements. Under this condition, prisms in the meshes may degenerate in certain ways, violating the so-called inscribed ball condition presented by P. G. Ciarlet (1978), but the interpolation error remains of the order in the -norm for sufficiently smooth functions.