An official website of the United States government
Finite element spaces on a tetrahedron are constructed for div div -conforming symmetric tensors in three dimensions. The key tools of the con- struction are the decomposition of polynomial tensor spaces and the charac- terization of the trace operators. First, the div div Hilbert complex and its corresponding polynomial complexes are presented. Several decompositions of polynomial vector and tensor spaces are derived from the polynomial com- plexes. Second, traces for the divdiv operator are characterized through a Green’s identity. Besides the normal-normal component, another trace involving combination of first order derivatives of the tensor is continuous across the face. Due to the smoothness of polynomials, the symmetric tensor element is also continuous at vertices, and on the plane orthogonal to each edge. Besides, a finite element for sym curl-conforming trace-free tensors is constructed following the same approach. Putting all together, a finite element div div complex, as well as the bubble functions complex, in three dimensions is established.
more »
« less
H(div)-conforming finite element tensors with constraints
A unified construction of H(div)-conforming finite element tensors, including vector element, symmetric matrix element, traceless matrix element, and, in general, tensors with linear constraints, is developed in this work. It is based on the geometric decomposition of Lagrange elements into bubble functions on each sub-simplex. Each tensor at a sub-simplex is further decomposed into tangential and normal components. The tangential component forms the bubble function space, while the normal component characterizes the trace. Some degrees of freedom can be redistributed to (n-1)-dimensional faces. The developed finite element spaces are H(div)-conforming and satisfy the discrete inf-sup condition. Intrinsic bases of the constraint tensor space are also established.
more »
« less
- PAR ID:
- 10618537
- Publisher / Repository:
- Elsevier B.V. (Netherlands)
- Date Published:
- Journal Name:
- Results in Applied Mathematics
- Volume:
- 23
- Issue:
- C
- ISSN:
- 2590-0374
- Page Range / eLocation ID:
- 100494
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
A finite element elasticity complex on tetrahedral meshes and the corresponding commutative diagram are devised. The H 1 H^1 conforming finite element is the finite element developed by Neilan for the velocity field in a discrete Stokes complex. The symmetric div-conforming finite element is the Hu-Zhang element for stress tensors. The construction of an H ( inc ) H(\operatorname {inc}) -conforming finite element of minimum polynomial degree 6 6 for symmetric tensors is the focus of this paper. Our construction appears to be the first H ( inc ) H(\operatorname {inc}) -conforming finite elements on tetrahedral meshes without further splitting. The key tools of the construction are the decomposition of polynomial tensor spaces and the characterization of the trace of the inc \operatorname {inc} operator. The polynomial elasticity complex and Koszul elasticity complex are created to derive the decomposition. The trace of the inc \operatorname {inc} operator is induced from a Green’s identity. Trace complexes and bubble complexes are also derived to facilitate the construction. Two-dimensional smooth finite element Hessian complex and div div \operatorname {div}\operatorname {div} complex are constructed.more » « less
-
An H(div)-conforming finite element method for the Biot’s consolidation mo- del is developed, with displacements and fluid velocity approximated by elements from BDM_k space. The use of H(div)-conforming elements for flow variables ensures the local mass conservation. In the H(div)-conforming approximation of displacement, the tan- gential components are discretised in the interior penalty discontinuous Galerkin frame- work,and the normal components across the element interfaces are continuous. Having introduced a spatial discretisation, we develop a semi-discrete scheme and a fully dis- crete scheme,prove their unique solvability and establish optimal error estimates for each variable.more » « less
-
Finite element de Rham complexes and finite element Stokes complexes with varying degrees of smoothness in three dimensions are systematically constructed in this paper. Smooth scalar finite elements in three dimensions are derived through a non-overlapping decomposition of the simplicial lattice. H(div)-conforming finite elements and H(curl)-conforming finite elements with varying degrees of smoothness are devised based on these smooth scalar finite elements. The finite element de Rham complexes with corresponding smoothness and commutative diagrams are induced by these elements. The div stability of the H(div)-conforming finite elements is established, and the exactness of these finite element complexes is proven.more » « less
-
We present an implementation of the trimmed serendipity finite element family, using the open-source finite element package Firedrake. The new elements can be used seamlessly within the software suite for problems requiring H 1 , H (curl), or H (div)-conforming elements on meshes of squares or cubes. To test how well trimmed serendipity elements perform in comparison to traditional tensor product elements, we perform a sequence of numerical experiments including the primal Poisson, mixed Poisson, and Maxwell cavity eigenvalue problems. Overall, we find that the trimmed serendipity elements converge, as expected, at the same rate as the respective tensor product elements, while being able to offer significant savings in the time or memory required to solve certain problems.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1016/j.rinam.2024.100494
