No Arabic abstract
The $H^m$-conforming virtual elements of any degree $k$ on any shape of polytope in $mathbb R^n$ with $m, ngeq1$ and $kgeq m$ are recursively constructed by gluing conforming virtual elements on faces in a universal way. For the lowest degree case $k=m$, the set of degrees of freedom only involves function values and derivatives up to order $m-1$ at the vertices of the polytope. The inverse inequality and several norm equivalences for the $H^m$-conforming virtual elements are rigorously proved. The $H^m$-conforming virtual elements are then applied to discretize a polyharmonic equation with a lower order term. With the help of the interpolation error estimate and norm equivalences, the optimal error estimates are derived for the $H^m$-conforming virtual element method.
Several div-conforming and divdiv-conforming finite elements for symmetric tensors on simplexes in arbitrary dimension are constructed in this work. The shape function space is first split as the trace space and the bubble space. The later is further decomposed into the null space of the differential operator and its orthogonal complement. Instead of characterization of these subspaces of the shape function space, characterization of the duals spaces are provided. Vector div-conforming finite elements are firstly constructed as an introductory example. Then new symmetric div-conforming finite elements are constructed. The dual subspaces are then used as build blocks to construct divdiv conforming finite elements.
Edge (or Nedelec) finite elements are theoretically sound and widely used by the computational electromagnetics community. However, its implementation, specially for high order methods, is not trivial, since it involves many technicalities that are not properly described in the literature. To fill this gap, we provide a comprehensive description of a general implementation of edge elements of first kind within the scientific software project FEMPAR. We cover into detail how to implement arbitrary order (i.e., $p$-adaptive) elements on hexahedral and tetrahedral meshes. First, we set the three classical ingredients of the finite element definition by Ciarlet, both in the reference and the physical space: cell topologies, polynomial spaces and moments. With these ingredients, shape functions are automatically implemented by defining a judiciously chosen polynomial pre-basis that spans the local finite element space combined with a change of basis to automatically obtain a canonical basis with respect to the moments at hand. Next, we discuss global finite element spaces putting emphasis on the construction of global shape functions through oriented meshes, appropriate geometrical mappings, and equivalence classes of moments, in order to preserve the inter-element continuity of tangential components of the magnetic field. Finally, we extend the proposed methodology to generate global curl-conforming spaces on non-conforming hierarchically refined (i.e., $h$-adaptive) meshes with arbitrary order finite elements. Numerical results include experimental convergence rates to test the proposed implementation.
We develop a geometrically intrinsic formulation of the arbitrary-order Virtual Element Method (VEM) on polygonal cells for the numerical solution of elliptic surface partial differential equations (PDEs). The PDE is first written in covariant form using an appropriate local reference system. The knowledge of the local parametrization allows us to consider the two-dimensional VEM scheme, without any explicit approximation of the surface geometry. The theoretical properties of the classical VEM are extended to our framework by taking into consideration the highly anisotropic character of the final discretization. These properties are extensively tested on triangular and polygonal meshes using a manufactured solution. The limitations of the scheme are verified as functions of the regularity of the surface and its approximation.
We present a class of discretisation spaces and H(div)-conformal elements that can be built on any polytope. Bridging the flexibility of the Virtual Element spaces towards the elements shape with the divergence properties of the Raviart-Thomas elements on the boundaries, the designed frameworks offer a wide range of H(div)-conformal discretisations. As those elements are set up through degrees of freedom, their definitions are easily amenable to the properties the approximated quantities are wished to fulfil. Furthermore, we show that one straightforward restriction of this general setting share its properties with the classical Raviart-Thomas elements at each interface, for any order and any polytopial shape. Then, we investigate the shape of the basis functions corresponding to particular elements in the two dimensional case.
We establish a new H2 Korns inequality and its discrete analog, which greatly simplify the construction of nonconforming elements for a linear strain gradient elastic model. The Specht triangle [41] and the NZT tetrahedron [45] are analyzed as two typical representatives for robust nonconforming elements in the sense that the rate of convergence is independent of the small material parameter. We construct new regularized interpolation estimate and the enriching operator for both elements, and prove the error estimates under minimal smoothness assumption on the solution. Numerical results are consistent with the theoretical prediction.