This paper reveals that the duality in finite element exterior calculus on simplicial meshes is fundamentally Hodge duality on the sphere, providing explicit correspondences and new duality isomorphisms between finite element spaces.
Contribution
The authors demonstrate that the duality in finite element exterior calculus is essentially Hodge duality on the sphere, clarifying the geometric nature of these relationships.
Findings
01
Established explicit correspondences between finite element spaces and differential forms on the sphere.
02
Derived new pointwise duality isomorphisms for finite element spaces.
03
Illustrated duality examples on the sphere.
Abstract
Finite element exterior calculus refers to the development of finite element methods for differential forms, generalizing several earlier finite element spaces of scalar fields and vector fields to arbitrary dimension n, arbitrary polynomial degree r, and arbitrary differential form degree k. The study of finite element exterior calculus began with the PrβΞk and PrββΞk families of finite element spaces on simplicial triangulations. In their development of these spaces, Arnold, Falk, and Winther rely on a duality relationship between PrβΞk and PΛr+k+1ββΞnβk and between PrββΞk and PΛr+kβΞnβk. In this article, we show that this duality relationship is, in essence, Hodge duality of differential forms on the standard n-sphere, disguised by aβ¦
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Full text
Duality in Finite Element Exterior Calculus and Hodge Duality on the Sphere
Yakov Berchenko-Kogan
Abstract.
Finite element exterior calculus refers to the development of finite element methods for differential forms, generalizing several earlier finite element spaces of scalar fields and vector fields to arbitrary dimension n, arbitrary polynomial degree r, and arbitrary differential form degree k. The study of finite element exterior calculus began with the PrβΞk and PrββΞk families of finite element spaces on simplicial triangulations. In their development of these spaces, Arnold, Falk, and Winther rely on a duality relationship between PrβΞk and PΛr+k+1ββΞnβk and between PrββΞk and PΛr+kβΞnβk. In this article, we show that this duality relationship is, in essence, Hodge duality of differential forms on the standard n-sphere, disguised by a change of coordinates. We remove the disguise, giving explicit correspondences between the PrβΞk, PrββΞk, PΛrβΞk and PΛrββΞk spaces and spaces of differential forms on the sphere. As a direct corollary, we obtain new pointwise duality isomorphisms between PrβΞk and PΛr+k+1ββΞnβk and between PrββΞk and PΛr+kβΞnβk. These isomorphisms can be implemented via a simple computation, which we illustrate with examples.
Key words and phrases:
finite element exterior calculus, Hodge duality, differential forms, finite element method
Given a simplicial triangulation T, Arnold, Falk, and Winther [4] constructed two families of spaces of piecewise polynomial differential forms on T, namely the PrβΞk(T) and PrββΞk(T) spaces, where r is the polynomial degree and k is the differential form degree. This task amounts to finding a way to uniquely specify a differential form bβ£Tβ on each simplex T of the triangulation T, while at the same time ensuring interelement continuity, in the sense that if T1β and T2β share a face f, then the restrictions of bβ£T1ββ and bβ£T2ββ to f agree. Arnold, Falk, and Winther accomplish this task by specifying b using the values of β«fβaβ§b for each face f of the triangulation T, where a comes from an appropriately chosen space of polynomial differential forms on f. Specifically, for bβPrβΞk(T), they take aβPr+kβdimfββΞdimfβk(f), and for bβPrββΞk(T), they take aβPr+kβdimfβ1βΞdimfβk(f). With these choices, given arbitrary desired values for β«fβaβ§b for all f and a, there is a unique differential form b on T with those values satisfying the interelement continuity condition above. See also [1, Theorem 3.5] for a short proof of this result.
As discussed in their later paper [5, Theorem 4.3], one can reduce the above task of finding appropriate spaces for a to a much simpler task. Namely, given a single simplex Tn of arbitrary dimension n, we consider the spaces PΛrβΞk(Tn) and PΛrββΞk(Tn) of polynomial differential forms on Tn whose restriction to the boundary βTn vanishes. For b coming from one of those spaces, the task is to find an appropriate space for a such that, given desired values of β«Tnβaβ§b for all a, there is a unique b with those values. To accomplish this task, for bβPΛrβΞk(Tn), we can take aβPr+kβnββΞnβk(Tn), and for bβPΛrββΞk(Tn), we can take aβPr+kβnβ1βΞnβk(Tn). As the authors discuss in [5], the choices of spaces for a in [4] follow by simply setting f=Tn.
As discussed in [4], we should think of the above claim as a duality result. For each bβPΛrβΞk(Tn), we have a linear functional aβ¦β«Tnβaβ§b. Hence, each bβPΛrβΞk(Tn) determines an element of the dual space Pr+kβnββΞnβk(Tn)β. Thus, we have a map PΛrβΞk(Tn)βPr+kβnββΞnβk(Tn)β, and likewise we have a map PΛrββΞk(Tn)βPr+kβnβ1βΞnβk(Tn)β. With this perspective, the statement that we can find a unique b that will give desired values for β«Tnβaβ§b for all a is simply the statement that these maps are isomorphisms.
In this paper, it will be convenient to consider the equivalent isomorphisms PrβΞk(Tn)βPΛr+k+1ββΞnβk(Tn)β and PrββΞk(Tn)βPΛr+kβΞnβk(Tn)β, which we can get from the ones above by taking the dual map and reindexing r and k. Additionally, rather that claiming that these maps are isomorphisms, we can also take the equivalent perspective of claiming that the bilinear paring (a,b)β¦β«Tnβaβ§b is nondegenerate. That is, given a nonzero aβPrβΞk(Tn), there exists a form bβPΛr+k+1ββΞnβk(Tn) such that β«Tnβaβ§b>0, and, conversely, for every b there exists such an a, and likewise for the other duality relationship.
1.1. Main results
In this article, we show that, with a change of coordinates, we can reveal the duality relationships between PrβΞk(Tn) and PΛr+k+1ββΞnβk(Tn) and between PrββΞk(Tn) and PΛr+kβΞnβk(Tn) to be Hodge duality on the standard n-sphere combined with multiplication by the bubble function uNβ:=u1ββ―un+1β defined in Notation 2.3. Specifically, using the change of coordinates, the Hodge star on the sphere βSnβ, and the bubble function uNβ, we define an injective tensorial map from k-forms Ξk(Tn) to (nβk)-forms Ξnβk(Tn). We show that, under this map, the image of PrβΞk(Tn) is PΛr+k+1ββΞnβk(Tn), and the image of PrββΞk(Tn) is PΛr+kβΞnβk(Tn), giving us natural isomorphisms between these pairs of spaces. This result appears as Corollary 3.3. We remark here that defining a natural map from PrβΞk(Tn) to PΛr+k+1ββΞnβk(Tn) is very different from defining a natural map from PrβΞk(Tn) to the dual space PΛr+k+1ββΞnβk(Tn)β. As discussed above, the latter map follows immediately from the paring β«Tnβaβ§b.
Our map also allows us to strengthen the above statement of nondegeneracy of the pairing (a,b)β¦β«Tnβaβ§b. Given a nonzero aβPrβΞk(Tn), rather than simply saying that there exists a form bβPΛr+k+1ββΞnβk(Tn) such that β«Tnβaβ§b>0, our isomorphism map explicitly constructs this form b. Moreover, the tensoriality of this map means that b only depends on a pointwise, in the sense that the value of b at xβTn only depends on the value of a at x, not on its values at other points of Tn. We have an analogous result for the other pairing, and both of these results appear as Corollary 3.4.
Computing our map is simple and straightforward to implement, as we illustrate in Examples 3.5 and 3.6. We also give Examples 3.7 and 3.8, where we specialize to the classical case of scalar fields, computing our isomorphism in the cases k=n and k=0.
The key ingredient is the transformation xiβ=ui2β, which sends the unit n-sphere u12β+β―+un+12β=1 to the standard n-simplex x1β+β―+xn+1β=1. This transformation induces a correspondence between differential forms on the simplex Tn and differential forms on the sphere Sn. In Theorems 3.1 and 3.2, we determine the spaces of differential forms on the sphere that correspond to the PrβΞk(Tn), PrββΞk(Tn), PΛrβΞk(Tn), and PΛrββΞk(Tn) spaces of differential forms on the simplex. The previously discussed Corollaries 3.3 and 3.4 quickly follow from this characterization. We note that these results are equally valid for an arbitrary simplex by mapping it to the standard simplex using barycentric coordinates.
The isomorphism maps PrβΞk(Tn)ββPΛr+k+1ββΞnβk(Tn) and PrββΞk(Tn)ββPΛr+kβΞnβk(Tn) given in this paper are the same maps as the ones given in my earlier preprint [6], and to the best of my knowledge are the only such maps in the literature that are defined pointwise. This article can be viewed as providing a new interpretation of these maps in terms of Hodge duality on the n-sphere. For a different construction of isomorphism maps between these spaces, see Martin Lichtβs recent work [9]. The isomorphisms constructed in [9, Equations (34) and (35)] involve decomposing differential forms in PrβΞk(Tn) and PrββΞk(Tn) as linear combinations in terms of canonical spanning sets defined in that paper. In contrast, the isomorphism maps in this article can be written directly in terms of the Hodge star, without needing to decompose the differential forms as linear combinations.
1.2. Outline of this paper
We discuss our notation and definitions in Section 2. In particular, we use a new definition of PrββΞk(Tn). We show in Appendix A that our definition is equivalent to the definition given by Arnold, Falk, and Winther. We also define several spaces of differential forms on the sphere.
In Section 3, we present the main results discussed above, along with examples that illustrate them. Specifically, in Theorem 3.1, we give the spaces of differential forms on the sphere that correspond to the PrβΞk(Tn), PrββΞk(Tn), PΛrβΞk(Tn), and PΛrββΞk(Tn) spaces of differential forms on the simplex via the coordinate transformation xiβ=ui2β. Then, in Theorem 3.2, we reexpress these spaces of differential forms on the sphere in terms of the Hodge star on the sphere βSnβ and the bubble function uNβ. From Theorem 3.2, we quickly prove the aforementioned Corollaries 3.3 and 3.4.
Section 4 is devoted to the proof of Theorem 3.1. Then, in Section 5, we study the PsββΞk(Sn) spaces of differential forms on the sphere given in Definition 2.9, and characterize these spaces in terms of the Hodge star βSnβ. Likewise, in Section 6, we study the P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Sn) spaces given in Definition 2.11, and characterize these spaces in terms of the bubble function uNβ. Together, these results allow us to restate Theorem 3.1 solely in terms of the space of polynomial differential forms on the sphere PsβΞk(Sn), the Hodge star βSnβ, and the bubble function uNβ, without reference to the PsββΞk(Sn) and P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Sn) spaces, giving us Theorem 3.2.
2. Preliminaries
In this section, we discuss the concepts that we use in our main results. We begin by setting our notation for differential forms on the simplex Tn, the sphere Sn, and Euclidean space Rn+1. Next, we define the spaces of polynomial differential forms that are the subject of this paper. We discuss how to compute the Hodge star on the sphere βSnβ, and we extend this operator to differential forms on Rn+1. Finally, we define notation for the transformation xiβ=ui2β, and we define even and odd differential forms on the sphere.
2.1. Notation
Notation 2.1**.**
Let Tn denote the standard n-simplex.
[TABLE]
Let Sn denote the unit sphere.
[TABLE]
Let S>0nβ denote the part of the unit sphere with strictly positive coordinates.
[TABLE]
Define T>0nβ similarly.
Notation 2.2**.**
Let Ξk(Tn) denote the space of differential k-forms on Tn. If aβΞk(Tn) and xβTn, let axβ denote a evaluated at x. That is, axβ is an antisymmetric k-linear tensor on the tangent space TxβTn.
Likewise, for Ξ±βΞk(Sn) and uβSn we have Ξ±uβ, an antisymmetric k-linear tensor on TuβSn. Similarly, for Ξ±^βΞk(Rn+1) and uβRn+1, we have Ξ±^uβ, an antisymmetric k-linear tensor on TuβRn+1β Rn+1.
Notation 2.3**.**
For I={i1β<i2β<β―<ikβ}β{1,β¦,n+1}, let uIβ denote the product ui1βββ―uikββ and let duIβ denote the wedge product dui1βββ§β―β§duikββ. Let N={1,β¦,n+1}, so
[TABLE]
In the literature, the function uNβ is called a bubble function.
Notation 2.4**.**
Given a vector field V on Rn+1 and Ξ±^βΞk(Rn+1), we will let iVβΞ±^βΞkβ1(Rn+1) denote the interior product of Ξ±^ with V.
We will use the same notation for the interior product on Tn and Sn.
Notation 2.5**.**
Let a^βΞk(Rn+1). Pulling back a^ via the inclusion map TnβͺRn+1, we obtain a differential form aβΞk(Tn). Following standard terminology, we refer to a as the restriction of a^ to Tn. We will also say that a^ is an extension of a.
We will also use this terminology for other inclusions, such as SnβͺRn+1.
2.2. Spaces of polynomial differential forms
Definition 2.6**.**
Let PrβΞk(Rn+1) denote the space of differential k-forms on Rn+1 whose coefficients are polynomials of degree at most r. Let PrβΞk(Tn) and PrβΞk(Sn) denote the restrictions of PrβΞk(Rn+1) to Tn and Sn, respectively. By convention, PrβΞk(Rn+1)={0} if r<0.
One must be careful with the definition of PrβΞk(Sn), as illustrated by the following example.
Example 2.7**.**
Using the notation (u,v,w) to represent a point in R3, we can consider the polynomial u3+uv2+uw2. Despite having degree three, this function is in P1βΞ0(Sn) because, on Sn, we have u3+uv2+uw2=u(u2+v2+w2)=u, which has degree one. In other words, as a function on Sn, u3+uv2+uw2 is the restriction of a degree one polynomial on Rn+1 to Sn.
Definition 2.8**.**
Let X be the radial vector field in Rn+1.
[TABLE]
When necessary to avoid confusion, we will also denote the radial vector field by
[TABLE]
Note that via the transformation xiβ=ui2β, we have U=2X.
Definition 2.9**.**
Let PrββΞk(Rn+1) denote those forms that are in the image under iXβ of forms of lower degree. That is,
[TABLE]
Let PrββΞk(Tn) and PrββΞk(Sn) denote the restrictions of PrββΞk(Rn+1) to Tn and Sn, respectively.
Note that this definition of PrββΞk(Tn) differs from that of Arnold, Falk, and Winther [4]. We show the equivalence of the two definitions in Proposition A.2.
We now discuss two different notions of βvanishing traceβ.
Definition 2.10**.**
Let PΛrβΞk(Tn) and PΛrββΞk(Tn) denote those forms in PrβΞk(Tn) and PrββΞk(Tn), respectively, whose restriction to βT vanishes. That is, aβPΛrβΞk(Tn) if axβ(V1β,β¦,Vkβ) for any xββTn and any V1β,β¦,Vkβ tangent to βT.
Note that saying that Ξ±^ is in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1) is stronger than simply saying that the restriction of Ξ±^ to each Ξiβ vanishes. Saying that Ξ±^βP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1) means that Ξ±^uβ(V1β,β¦,Vkβ)=0 for any uβΞiβ and arbitrary vectors V1β,β¦,VkββTuβRn+1, not just vectors tangent to Ξiβ as with restriction.
Example 2.12**.**
Using coordinates (u,v) on R2, consider Ξ±^=udvβP1βΞ1(R2). Observe that the restriction of Ξ±^ to {v=0} vanishes because the tangent space of {v=0} is spanned by βuββ, and udv(βuββ)=0. However, Ξ±^β/P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vss1βΞ1(R2) because it does not vanish on vectors normal to {v=0}: we have udv(βvββ)=u, which is not identically zero on the line {v=0}.
Note also that if Ξ±^βPsβΞk(Rn+1) is an extension of Ξ±βPsβΞk(Sn), checking that Ξ±uβ=0 only involves checking that Ξ±uβ(V1β,β¦,Vkβ)=0 for vectors tangent to Sn, whereas checking that Ξ±^uβ=0 involves checking that Ξ±^uβ(V1β,β¦,Vkβ)=0 for all vectors.
Example 2.13**.**
Using coordinates (u,v) on R2, consider Ξ±^=vdvβP1βΞ1(R2), and let Ξ± be the restriction of Ξ±^ to S1={(u,v)β£u2+v2=1}.
We claim that, on the circle, we have that Ξ±βP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vss1βΞ1(S1). Indeed, u2+v2=1 implies that udu+vdv=0, so Ξ±=βudu. Thus Ξ± vanishes when we set u=0 or v=0.
In contrast, Ξ±^β/P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vss1βΞ1(R2), since Ξ±^ does not vanish when u=0.
2.3. The Hodge star
Notation 2.14**.**
Let βSnβ:Ξk(Sn)βΞnβk(Sn) denote the Hodge star with respect to the standard metric on Sn, and let βRn+1β:Ξk(Rn+1)βΞn+1βk(Rn+1) denote the standard Hodge star on Rn+1.
As we will see, computing βSnβΞ± is a straightforward computation.
Definition 2.15**.**
Let Ξ½ denote the outward unit conormal to the sphere, so
[TABLE]
Proposition 2.16**.**
Let Ξ±βΞk(Sn). Let Ξ±^βΞk(Rn+1) be an extension of Ξ±. Then βSnβΞ± is the restriction to Sn of
[TABLE]
Proof.
Proposition B.2 gives the corresponding result for hyperplanes of oriented inner product spaces. We apply it to the hyperplane TuβSnβTuβRn+1.
β
Example 2.17**.**
If n=2 and k=1, we can think of Ξ±βΞ1(S2) as a vector field on the sphere. With this interpretation, βS2βΞ± rotates Ξ± by 90β counterclockwise, and Proposition 2.16 states that one way to compute βS2βΞ± is to take the cross product of the normal vector with Ξ±.
Proposition 2.16 motivates extending the definition of βSnβ:Ξk(Sn)βΞnβk(Sn) to βSnβ:Ξk(Rn+1)βΞnβk(Rn+1) as follows.
Definition 2.18**.**
For Ξ±^βΞk(Rn+1), let βSnβΞ±^βΞnβk(Rn+1) denote
[TABLE]
With this definition, Proposition 2.16 states that if Ξ±βΞk(Sn) is the restriction of Ξ±^, then βSnβΞ± is the restriction of βSnβΞ±^.
2.4. The coordinate transformation between Sn and Tn
Definition 2.19**.**
Consider the transformation Ξ¦:Rn+1βRn+1 defined by
[TABLE]
Observe that xβTn if and only if uβSn. In fact, Ξ¦ is a diffeomorphism when restricted to S>0nββT>0nβ.
We will use the notation Ξ¦β to refer to the pullback map both in the context of Ξ¦:Rn+1βRn+1 and in the context of Ξ¦:SnβTn, so we have pullback maps Ξ¦β:Ξk(Rn+1)βΞk(Rn+1) and Ξ¦β:Ξk(Tn)βΞk(Sn).
Observe that if aβΞk(Tn) is the restriction of a^βΞk(Rn+1) to Tn, then Ξ¦βaβΞk(Sn) is the restriction of Ξ¦βa^βΞk(Rn+1) to Sn.
2.5. Even and odd differential forms
Because Ξ¦(u1β,β¦,uiβ,β¦,un+1β)=Ξ¦(u1β,β¦,βuiβ,β¦,un+1β), any differential form in the image of Ξ¦β is invariant under all coordinate reflections. We call such forms even in all variables, or simply even.
Definition 2.20**.**
Let the ith coordinate reflection Riβ:Rn+1βRn+1 be the map
[TABLE]
We call Ξ±^βΞk(Rn+1)even in all variables or simply even if RiββΞ±^=Ξ±^ for every i. Similarly, we call Ξ±^βΞk(Rn+1)odd in all variables or simply odd if RiββΞ±^=βΞ±^ for every i. We denote these spaces by Ξekβ(Rn+1) and Ξokβ(Rn+1), respectively.
We define Ξekβ(Sn) and Ξokβ(Sn) similarly. We will also use the notation PsβΞekβ(Sn) and PsβΞokβ(Sn) to denote the even and odd forms in PsβΞk(Sn), respectively, and likewise for the other spaces of polynomial differential forms.
As with functions in one variable, we can take even and odd parts.
Definition 2.21**.**
If Ξ±^βΞk(Rn+1), let the even partΞ±^eββΞekβ(Rn+1) of Ξ±^ denote the average of all possible reflections of Ξ±^. That is,
[TABLE]
Likewise, let the odd part of Ξ±^oββΞokβ(Rn+1) of Ξ±^ denote the signed average
[TABLE]
We use the same equations to define the even and odd parts of a differential form on the sphere, Ξ±βΞk(Sn).
If Ξ±^ is a differential form with polynomial coefficients, Ξ±^eβ simply extracts those terms that are even k-forms, and Ξ±^oβ extracts those terms that are odd k-forms. Note that Ξ±^ is not the sum of its even and odd parts; there will generally be terms that are even in some variables and odd in others.
Proposition 2.22**.**
Restriction commutes with taking even and odd parts. That is, if Ξ±βΞk(Sn) is the restriction of Ξ±^βΞk(Rn+1) to Sn, then Ξ±eβ and Ξ±oβ are the restrictions of Ξ±^eβ and Ξ±^oβ, respectively.
Proof.
If Ξ± is the restriction of Ξ±^, then RiββΞ± is the restriction of RiββΞ±^. Thus, when we restrict equations (1) and (2) to Sn, we obtain the corresponding equations for Ξ±eβ and Ξ±oβ.
β
In the following proposition, we summarize how the operations we have considered interact with even and odd parts.
Proposition 2.23**.**
The following operations preserve even and odd parts.
[TABLE]
The following operations interchange even and odd parts.
[TABLE]
Proof.
Observe that
[TABLE]
Reflections reverse orientation, so
[TABLE]
Consequently, since βSnβΞ±^=βRn+1β(Ξ½β§Ξ±^),
[TABLE]
The result follows by applying these operations to equations (1) and (2).
β
3. Results and examples
We state our main result that the coordinate transformation Ξ¦ induces correspondences between spaces of polynomial differential forms on Tn and polynomial differential forms on Sn. As a consequence of these correspondences, we easily obtain the duality relationships between the P and Pβ spaces. We work with the standard simplex, but all of our results apply equally well to an arbitrary simplex, by using barycentric coordinates.
Theorem 3.1**.**
The map Ξ¦β:Ξk(Tn)βΞk(Sn) induced by the coordinate transformation Ξ¦ gives the following correspondences between polynomial differential forms on the simplex and polynomial differential forms on the sphere.
[TABLE]
where these spaces of differential forms on Sn are defined in Definitions 2.6, 2.9, 2.11, and 2.20.
Proof.
We prove this theorem in Section 4. Specifically, the four isomorphisms are proved in Theorems 4.3, 4.5, 4.6, and 4.7, respectively.
β
The results of Sections 5 and 6 allow us to express the above spaces of differential forms on Sn solely in terms of PsβΞe/okβ(Sn), the Hodge star βSnβ, and the bubble function uNβ. We obtain the following version of the theorem.
Theorem 3.2**.**
The map Ξ¦β:Ξk(Tn)βΞk(Sn) induced by the coordinate transformation Ξ¦ gives the following correspondences between polynomial differential forms on the simplex and polynomial differential forms on the sphere.
[TABLE]
Proof.
The first isomorphism is the same as in Theorem 3.1. In light of Theorem 3.1, in order to obtain the second isomorphism, we must show that PsββΞekβ(Sn)=βSnβPsβ1βΞonβkβ(Sn). We prove this claim in Corollary 5.3. To prove the third isomorphism, we must show that P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞekβ(Sn)=uNβP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβnβ1βΞokβ(Sn) if sβ₯2+k. This claim follows from Corollary 6.11. Finally, we can prove the last isomorphism by showing that P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞekβ(Sn)=uNβPsβnβ1ββΞokβ(Sn) and Psβnβ1ββΞokβ(Sn)=βSnβPsβnβ2βΞenβkβ(Sn). We prove the first claim in Corollary 6.8, and the second claim follows from Corollary 5.3.
β
As an immediate corollary, we obtain an explicit pointwise construction of the duality isomorphisms of Arnold, Falk, and Winther.
Corollary 3.3**.**
The pointwise-defined map (Ξ¦β)β1β(uNββSnβ)βΞ¦β is an isomorphism between the following spaces.
[TABLE]
These spaces are isomorphic via pointwise-defined maps to P2r+kβΞekβ(Sn) and P2r+kβ1βΞonβkβ(Sn), respectively.
Proof.
These isomorphisms follow directly from Theorem 3.2 along with the fact that βSnβ(βSnβΞ±)=(β1)k(nβk)Ξ±, so βSnβ(βSnβP2r+kβ1βΞonβkβ(Sn))=P2r+kβ1βΞonβkβ(Sn). The operations of pullback, Hodge star, and multiplication by a scalar function are all pointwise-defined maps.
β
We show that these isomorphisms give forms that are dual to one another with respect to integration.
Corollary 3.4**.**
Consider nonzero aβPrβΞk(Tn) and bβPΛr+k+1ββΞnβk(Tn) corresponding to one another via the first isomorphism in Corollary 3.3. Then β«Tnβaβ§b>0. Moreover, a and b depend on each other only pointwise, in the sense that for xβTn, bxβ depends only on axβ and vice versa.
The same holds for aβPrββΞk(Tn) and bβPΛr+kβΞnβk(Tn) corresponding to one another via the second isomorphism in Corollary 3.3.
Proof.
Let Ξ±=Ξ¦βa and Ξ²=Ξ¦βb, so Ξ²=uNβ(βSnβΞ±). By u-substitution, we have
[TABLE]
because Ξ± is not identically zero and uNβ>0 on S>0nβ.
β
We provide an example of each of the two isomorphisms in the case n=2. To simplify notation, we use coordinates (x,y,z) and (u,v,w) on R3. We compute the Hodge star βSnβ using the formula in Proposition 2.16. In the case of one-forms on the two-sphere as in these examples, the Hodge star can be interpreted as rotating a vector field on the sphere 90 degrees counterclockwise, and the formula in Proposition 2.16 states that one can do so by taking the cross product of the unit normal with the vector field.
Example 3.5**.**
[TABLE]
Example 3.6**.**
[TABLE]
In the two examples below, we apply the isomorphism in Corollary 3.3 to the classical case of scalar fields, that is, n-forms and [math]-forms.
Example 3.7**.**
Consider a general polynomial n-form on Tn, which can be expressed as
Consider now a general polynomial [math]-form on Tn, expressed as
[TABLE]
Then we have
[TABLE]
We can verify that, on the sphere, uNβvolSnβ=u2ββ―un+1βdu2ββ§β―β§dun+1β. Indeed, using equation (3), we see that taking βSnβ of both sides gives uNβ, so these n-forms must be equal on the sphere. We conclude that
Theorem 3.1 claims four isomorphisms. We prove the isomorphism for PrβΞk(Tn) in Section 4.1, for PrββΞk(Tn) in Section 4.2, and for PΛrβΞk(Tn) and PΛrββΞk(Tn) in Section 4.3.
4.1. The correspondence for PrβΞk(Tn)
We first prove the result for forms on Rn+1, and then use it to conclude the isomorphism between forms on Tn and forms on Sn.
Proposition 4.1**.**
The map Ξ¦β is an isomorphism between the following spaces.
[TABLE]
Proof.
Assume a^βPrβΞk(Rn+1), and let Ξ±^=Ξ¦βa^. Using Notation 2.3, let
[TABLE]
where pIβ is a polynomial in x of degree at most r. Since xiβ=ui2β, we have dxiβ=2uiβduiβ, so
[TABLE]
The coefficient pIβ(u12β,β¦,un+12β) is even of degree at most 2r, and uIβduIββPkβΞekβ(Rn+1), so Ξ±^βP2r+kβΞekβ(Rn+1).
Conversely, assume Ξ±^βP2r+kβΞekβ(Rn+1). Let
[TABLE]
Observe that Riββ(duIβ)=βduIβ if iβI and Riββ(duIβ)=duIβ if iβ/I. Since RiββΞ±^=Ξ±^, we conclude that RiββqIβ=βqIβ if iβI and RiββqIβ=qIβ if iβ/I.
In other words, qIβ is odd in the variable uiβ for iβI and even in uiβ for iβ/I.
Thus, qIβ is divisible by uiβ for iβI, so we can write qIβ=2krIβuIβ for a polynomial rIβ, including a factor of 2k for convenience. We see that rIβ is even in all variables uiβ for 1β€iβ€n+1. Thus, we can write rIβ(u1β,β¦,un+1β)=pIβ(u12β,β¦,un+12β). We conclude that
[TABLE]
which is Ξ¦βa^ for a^=βIβpIβdxIβ as shown above.
β
We use this result to show that this isomorphism holds between forms on Tn and Sn. We begin with a simple lemma.
Lemma 4.2**.**
The map Ξ¦β:Ξk(Tn)βΞk(Sn) is injective.
Proof.
Let aβΞk(Tn), and assume that Ξ¦βa=0. Because Ξ¦:S>0nββT>0nβ is a diffeomorphism, we know that Ξ¦β:Ξk(T>0nβ)βΞk(S>0nβ) is a bijection, so a must be zero on T>0nβ. Because a is continuous, it must therefore be zero on all of Tn.
β
The map Ξ¦β is an isomorphism between the following spaces.
[TABLE]
Proof.
Let aβPrβΞk(Tn). By definition, there exists an extension a^βPrβΞk(Rn+1). Let Ξ±^=Ξ¦βa^, which is in P2r+kβΞekβ(Rn+1) by Proposition 4.1. Let Ξ± be the restriction of Ξ±^ to Sn, so we have Ξ±=Ξ¦βa, and Ξ±βP2r+kβΞekβ(Sn) by definition. We conclude that Ξ¦β does indeed map PrβΞk(Tn) to P2r+kβΞekβ(Sn), and we know that this map is injective by Lemma 4.2.
To show surjectivity, let Ξ±βP2r+kβΞekβ(Sn). By definition, there exists an extension Ξ±^βP2r+kβΞekβ(Rn+1). By Proposition 4.1, there exists a^βPrβΞk(Rn+1) such that Ξ±^=Ξ¦βa^. Letting a be the restriction of a^ to Tn, we see that aβPrβΞk(Tn) and Ξ±=Ξ¦βa.
β
4.2. The correspondence for PrββΞk(Tn)
As before, we begin by considering forms on Rn+1. We first prove that the change of coordinates induces an isomorphism between PrββΞk(Rn+1) and P2r+kββΞekβ(Rn+1), and then conclude that it also induces an isomorphism between PrββΞk(Tn) and P2r+kββΞekβ(Sn).
Proposition 4.4**.**
The map Ξ¦β is an isomorphism between the following spaces.
[TABLE]
Proof.
The key fact needed here is that Ξ¦ preserves the radial vector field up to a constant scalar factor. That is, one can compute that Ξ¦ββU=2X.
If a^βPrββΞk(Rn+1), then by definition a^=iXβb^ for some b^βPrβ1βΞk+1(Rn+1). Let Ξ²^β=Ξ¦βb^, which is in P2rβ2+k+1βΞek+1β(Rn+1) by Proposition 4.1. Then
[TABLE]
Thus Ξ¦βa^=21βiUβΞ²^β. Since Ξ²^ββP2r+kβ1βΞek+1β(Rn+1), we know that Ξ¦βa^βP2r+kββΞekβ(Rn+1) by Definition 2.9 and Proposition 2.23.
Conversely, if Ξ±^βP2r+kββΞekβ(Rn+1), then by definition Ξ±^=21βiUβΞ²^ββ² for some Ξ²^ββ²βP2r+kβ1βΞk+1(Rn+1). Let Ξ²^β be the even part of Ξ²^ββ². By Proposition 2.23, Ξ±^=Ξ±^eβ=21βiUβΞ²^βeβ²β=21βiUβΞ²^β. By Proposition 4.1, Ξ²^β=Ξ¦βb^ for some b^βPrβ1βΞk+1(Rn+1). Set a^=iXβb^. Then a^βPrββΞk(Rn+1) by definition, and Ξ¦βa^=Ξ±^ by equation (4).
β
We conclude the corresponding isomorphism between forms on Tn and forms on Sn.
Theorem 4.5** (Theorem 3.1, second isomorphism).**
The map Ξ¦β is an isomorphism between the following spaces.
[TABLE]
Proof.
Assume aβPrββΞk(Tn). By definition, a has an extension a^βPrββΞk(Rn+1). Then Ξ¦βa^ is in P2r+kββΞekβ(Rn+1) by Proposition 4.4. Since Ξ¦βa is the restriction of Ξ¦βa^ to Sn, we have that Ξ¦βaβP2r+kββΞekβ(Sn) by definition.
Conversely, if Ξ±βP2r+kββΞekβ(Sn), then by definition it has an extension Ξ±^βP2r+kββΞekβ(Rn+1). By Proposition 4.4, Ξ±^=Ξ¦βa^ for a^βPrββΞk(Rn+1). Letting a be the restriction of a^ to Tn, we have that aβPrββΞk(Tn) by definition, and Ξ±=Ξ¦βa.
β
4.3. The correspondence for PΛrβΞk(Tn) and PΛrββΞk(Tn)
We now show that forms on Tn with vanishing trace correspond to forms on Sn with vanishing trace. However, the notions of βvanishing traceβ are different; see Definitions 2.10 and 2.11 and the discussion and examples that follow.
By definition, aΛ=0. Since Ξ¦ maps Sinβ1β to Tinβ1β, we have that Ξ±Λ=Ξ¦βaΛ=0. In other words, we know that Ξ± vanishes on vectors tangent to Sinβ1β. It remains to to show that Ξ± vanishes on the vector βuiβββ normal to Sinβ1β, which we can do using the fact that Ξ± is even. More precisely, let iβuiββββΞ±β denote the restriction of the interior product iβuiββββΞ± to Sinβ1β. We must check iβuiββββΞ±β=0.
Applying the reflection Riβ, we have that R_{i}^{*}\bigl{(}\frac{\partial{}}{\partial{u_{i}}}\bigr{)}=-\frac{\partial{}}{\partial{u_{i}}}, so R_{i}^{*}\Bigl{(}i_{\frac{\partial{}}{\partial{u_{i}}}}\alpha\Bigr{)}=-i_{\frac{\partial{}}{\partial{u_{i}}}}\alpha because Ξ± is even. On the other hand, the reflection Riβ fixes Sinβ1β, so the restrictions of R_{i}^{*}\Bigl{(}i_{\frac{\partial{}}{\partial{u_{i}}}}\alpha\Bigr{)} and iβuiββββΞ± to Sinβ1β are equal. We conclude that iβuiββββΞ±β=βiβuiββββΞ±β, so iβuiββββΞ±β=0, as desired.
Conversely, assume that Ξ±βP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vss2r+kβΞekβ(Sn). We know that Ξ±=Ξ¦βa for aβPrβΞk(Tn), and we aim to show that aβPΛrβΞk(Tn). We have in particular that the restriction of Ξ± to Sinβ1β is zero, that is, Ξ±Λ=0 in the above notation. Thus, Ξ¦βaΛ=Ξ±Λ=0. Applying Lemma 4.2 to Ξ¦β:Ξk(Tinβ1β)βΞk(Sinβ1β), we conclude that aΛ=0. Thus, the restriction of a to Tinβ1β vanishes for all i, so we conclude that aβPΛrβΞk(Tn), as desired.
β
The map Ξ¦β is an isomorphism between the following spaces.
[TABLE]
Proof.
Theorem 4.6 tells us that Ξ¦β is an isomorphism between PΛrβΞk(Tn) and P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vss2r+kβΞekβ(Sn), and Theorem 4.5 tells us that Ξ¦β is an isomorphism between PrββΞk(Tn) and P2r+kββΞekβ(Sn). Taking the intersection, we obtain the desired result.
β
5. Characterizations of PsββΞk(Sn)
In Section 4.2, we showed that Ξ¦β is an isomorphism between PrββΞk(Tn) and P2r+kββΞekβ(Sn). To complete the proof of the correspondence in Theorem 3.2 for the PrββΞk(Tn) spaces, we show in this section that P2r+kββΞekβ(Sn)=βSnβP2r+kβ1βΞonβkβ(Sn). As before, we first begin with the corresponding result on Rn+1.
Proposition 5.1**.**
The space PsββΞk(Rn+1) is the image under βSnβ of differential forms of lower polynomial degree. That is,
[TABLE]
Proof.
By definition, PsββΞk(Rn+1)=iUβPsβ1βΞk+1(Rn+1). Thus, we aim to show that any form that can be expressed as βSnβΞ³^β for Ξ³^ββPsβ1βΞnβk(Rn+1) can be expressed as iUβΞ²^β for Ξ²^ββPsβ1βΞk+1(Rn+1), and vice versa.
Because Ξ½ is dual to U with respect to the standard metric, we use Proposition B.1 to compute that for Ξ³^ββΞnβk(Rn+1),
[TABLE]
Thus, setting Ξ²^β:=(β1)nβkβRn+1βΞ³^β, we have βSnβΞ³^β=iUβΞ²^β. If Ξ³^ββPsβ1βΞnβk(Rn+1), then Ξ²^ββPsβ1βΞk+1(Rn+1). Conversely, if Ξ²^ββPsβ1βΞk+1(Rn+1), then Ξ³^β=(β1)nβkβRn+1β1βΞ²^ββPsβ1βΞnβk(Rn+1).
β
We immediately obtain the corresponding result for forms on Sn.
Proposition 5.2**.**
The space PsββΞk(Sn) is the image under βSnβ of differential forms of lower polynomial degree. That is,
[TABLE]
Proof.
The restriction of PsββΞk(Rn+1) to the sphere is PsββΞk(Sn) by definition. The restriction of βSnβPsβ1βΞnβk(Rn+1) to the sphere is βSnβPsβ1βΞnβk(Sn) because if Ξ³^β is an extension of Ξ³, then βSnβΞ³^β is an extension of βSnβΞ³.
β
Corollary 5.3**.**
We have that
[TABLE]
Proof.
We use Propositions 2.22 and 2.23 to take the even and odd parts of both sides of equation (5).
β
We will also need the following characterization of PsββΞk(Rn+1).
By definition PsββΞk(Rn+1)=iUβPsβ1βΞk+1(Rn+1). Since iUββiUβ=0, we conclude that PsββΞk(Rn+1)βkeriUβ. Since iUβ vanishes on [math]-forms, we also have that PsββΞ0(Rn+1)+P0βΞ0(Rn+1)βPsβΞ0(Rn+1)=keriUβ.
Conversely, let Ξ±^βkeriUβ. Assuming sβ₯0, decompose Ξ±^ into homogeneous components. That is, write Ξ±^=Ξ±^sβ+β―+Ξ±^0β, where the coefficients of Ξ±^jβ are homogeneous polynomials of degree j. Observe that the coefficients of iUβΞ±^jβ are homogeneous polynomials of degree j+1. Thus, iUβΞ±^=0 implies that iUβΞ±^jβ=0 for all j.
We show that Ξ±^jββPsββΞk(Rn+1) by showing that the Lie derivative LUβΞ±^jβ is in PsββΞk(Rn+1) and that LUβΞ±^jβ is a constant multiple of Ξ±^jβ.
Using the Cartan formula, we have that
[TABLE]
We have that Ξ±^jββPjβΞk(Rn+1)βPsβΞk(Rn+1), so dΞ±^jββPsβ1βΞk+1(Rn+1), and so LUβΞ±^jββPsββΞk(Rn+1) by definition.
To show that LUβΞ±^jβ is a constant multiple of Ξ±^jβ, we start with the fact that LUβuiβ=uiβ. Because the Lie derivative commutes with the exterior derivative d, we then have LUβ(duiβ)=duiβ. The form Ξ±^jβ is a linear combination of products of j terms of the form uiβ and k terms of the form duiβ. Thus, the product rule for the Lie derivative gives us that LUβΞ±^jβ=(j+k)Ξ±^jβ.
Except for the case j=k=0, we can divide LUβΞ±^jβ by j+k to conclude that Ξ±^jββPsββΞk(Rn+1). When kβ₯1, this means that Ξ±^jββPsββΞk(Rn+1) for all j, so Ξ±^βPsββΞk(Rn+1). When k=0, we have that Ξ±^jββPsββΞ0(Rn+1) for jβ₯1 and that Ξ±^0ββP0βΞ0(Rn+1), so Ξ±^βPsββΞ0(Rn+1)+P0βΞ0(Rn+1).
β
6. Characterizations of ${\mathop{\kern 0.0pt\mathcal{P}}\limits^{\vbox to-5.5757pt{
In Section 4.3, we showed that Ξ¦β is an isomorphism between PΛrβΞk(Tn) and P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞekβ(Sn), where s=2r+k. To complete the proof of the correspondence in Theorem 3.2 for the PΛrβΞk(Tn) and PΛrββΞk(Tn) spaces, we show in this section that P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞekβ(Sn)=uNβPsβnβ1βΞokβ(Sn) and P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞekβ(Sn)=uNβPsβnβ1ββΞokβ(Sn).
This section is structured as follows. We first prove the corresponding results for forms on Rn+1, showing that forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1) and P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Rn+1) are divisible by uNβ and determining the quotient. However, unlike in the previous section, the desired results for forms on Sn do not immediately follow from these results for forms on Rn+1. The missing ingredient is showing that forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Sn) can be extended to forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1) and that forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Sn) can be extended to forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Rn+1). This claim is subtle, and in fact fails for exceptional values of s and k.
Example 6.1**.**
Using coordinates (u,v) on R2, consider vdv. This form is in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vss1βΞ1(S1) as seen in Example 2.13, but it cannot be extended to a form in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vss1βΞ1(R2). Indeed, as we will show in Proposition 6.2, a form in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞ1(R2) must be divisible by uv and thus have degree at least two.
The situation is simpler for the PsββΞk(Sn) spaces compared to the PsβΞk(Sn) spaces: we will show that forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Sn) can always be extended to forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Rn+1). We can then use the result that forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Rn+1) are divisible by uNβ to get our desired result that forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Sn) are divisible by uNβ. By taking the Hodge dual of this result, we can show that forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Sn) are also divisible by uNβ, without needing to show that these forms can be extended to forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1).
6.1. Differential forms on Rn+1
Proposition 6.2**.**
We have that
[TABLE]
Proof.
If Ξ±^βP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1), then each polynomial coefficient of Ξ±^ vanishes when we set uiβ=0, so each coefficient is divisible by uiβ. Thus, Ξ±^ is divisible by uNβ. Conversely, if Ξ±^βuNβPsβnβ1βΞk(Rn+1), then Ξ±^βP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1) because uNβ=0 whenever uiβ=0.
β
We prove the corresponding claim for P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Rn+1), but there is an exception due to the fact that uNββP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vssn+1ββΞ0(Rn+1) but 1β/P0ββΞ0(Rn+1).
By Proposition 5.4, proving the above claim is equivalent to showing that
[TABLE]
Let Ξ±^βP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Rn+1). Because Ξ±^βPsββΞk(Rn+1), we know that iUβΞ±^=0. Because Ξ±^βP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1), we know by Proposition 6.2 that Ξ±^=uNβΞ²^β for Ξ²^ββPsβnβ1βΞk(Rn+1). It remains to show that Ξ²^ββkeriUβ. We have
[TABLE]
Thus uNβ(iUβΞ²^β)=0. Dividing by uNβ, we conclude that iUβΞ²^β=0, as desired.
Conversely, let Ξ±^=uNβΞ²^β where Ξ²^ββkeriUβ. Proposition 6.2 tells us that Ξ±^βP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1). Meanwhile, equation (8) gives us that iUβΞ±^=0, so Proposition 5.4 tells us that Ξ±^βPsββΞk(Rn+1). (We can rule out the P0βΞ0(Rn+1) summand in Proposition 5.4 because Ξ±^ is divisible by uNβ.)
β
6.2. Extending forms on Sn to forms on Rn+1
We would like to use these results for Rn+1 to show the corresponding results for the sphere. It is clear that the restriction of a form in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1) to Sn must be in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Sn). It is less clear that forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Sn) can be extended to forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1); in fact, in exceptional cases, they cannot, as in Example 6.1 and more generally in Example 6.13. However, as we will prove in this subsection, for the P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Sn) spaces, this extension result always holds; there are no exceptional cases.
We begin with understanding the situation when k=0, with an additional parity hypothesis that we later remove without difficulty.
Lemma 6.4**.**
Let pβP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞ0(Sn). By definition, p has an extension p^ββ²βPsβΞ0(Rn+1). Assume that the total degree of every term of p^ββ² has the same parity, either all even or all odd. Then p has an extension p^ββP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞ0(Rn+1).
Proof.
The main idea of the proof is that we obtain p^β by homogenizing p^ββ². That is, we multiply each term of p^ββ² by an appropriate power of r2:=u12β+β―+un+12β.
More precisely, we construct p^β as follows. Without loss of generality, assume that s=degp^ββ². We decompose p^ββ² into homogeneous components, so we have either of the two cases:
[TABLE]
where p^βjβ²β is a homogeneous polynomial of degree j. We multiply each term of p^ββ² by an appropriate power of r2 in order to make the polynomial homogeneous. That is, set
We now prove this extension result for P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Sn), again with an additional parity hypothesis that we remove later.
Lemma 6.5**.**
Let Ξ±βP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Sn). By definition, Ξ± has an extension Ξ±^β²βPsββΞk(Rn+1). Assume that the total degree of every term of every polynomial coefficient of Ξ±^β² has the same parity, either all even or all odd. Then Ξ± has an extension Ξ±^βP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Rn+1).
To show that Ξ±^βPsββΞk(Rn+1), for uβSn, we have that iUβΞ±^uβ=iUβΞ±^uβ²β=0. Since Ξ±^ is homogeneous, we conclude that iUβΞ±^uβ=0 for all uβRn+1. We then conclude that Ξ±^βP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Rn+1) by Proposition 5.4. (We can rule out the P0βΞ0(Rn+1) summand of Proposition 5.4 using the fact that Ξ±^ is homogeneous.)
β
We now remove the parity hypothesis.
Proposition 6.6**.**
Any differential form in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Sn) can be extended to a differential form in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Rn+1).
Proof.
Let Ξ±βP\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Sn). Let
[TABLE]
Then Ξ² and Ξ³ satisfy the hypotheses of Lemma 6.5, because if Ξ±^β²βPsββΞk(Rn+1) is an extension of Ξ±, then Ξ²^βuβ²β=21β(Ξ±^uβ²β+Ξ±^βuβ²β) and Ξ³^βuβ²β=21β(Ξ±^uβ²ββΞ±^βuβ²β) are extensions of Ξ² and Ξ³ with the required parity property. Thus, we can extend Ξ² and Ξ³ to Ξ²^β and Ξ³^β in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Rn+1) and set Ξ±^=Ξ²^β+Ξ³^β.
β
6.3. Differential forms on Sn
Now that we have this extension result, we can use the fact that forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Rn+1) are divisible by uNβ to conclude that the same holds for forms in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Sn).
Moreover, P0βΞ0(Sn)βP2ββΞ0(Sn), so equation (9) holds even if k=0 as long as sβ₯n+3.
Proof.
Proposition 6.3 gives the corresponding result for forms on Rn+1. Restricting both sides to Sn using Proposition 6.6, we obtain the desired result.
To show that P0βΞ0(Sn)βP2ββΞ0(Sn), observe that iUβΞ½=u12β+β―+un+12β is in P2ββΞ0(Rn+1) by definition and restricts to the constant function 1 on Sn.
β
We take the even part of the above result.
Corollary 6.8**.**
For all k and s, we have
[TABLE]
Proof.
We use Propositions 2.22 and 2.23 to take the even part of the equations in Proposition 6.7.
The result holds even when k=0 because constant functions are even, so the odd part of P0βΞ0(Sn) is zero.
β
We now use Hodge duality to prove the corresponding result for P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Sn).
Lemma 6.9**.**
We have
[TABLE]
Proof.
Proposition 5.2 gives us PsββΞk(Sn)=βSnβPsβ1βΞnβk(Sn), and we have Ξ±uβ=0 if and only if βSnβΞ±uβ=0.
β
Lemma 6.9 and Proposition 5.2 characterize the above spaces as images under βSnβ, giving us
[TABLE]
Applying βSnβ to both sides and using the fact that βSnβ(βSnβΞ±)=(β1)k(nβk)Ξ±, we obtain
[TABLE]
Meanwhile, in the case k=n and sβ₯n, Proposition 6.7 instead gives us
[TABLE]
Taking the Hodge star as above, we instead obtain
[TABLE]
To see that volSnββP1βΞn(Sn), we have that volSnβ=βSnβ1, which means that volSnβ is the restriction to Sn of βRn+1βΞ½βP1βΞn(Rn+1).
β
Corollary 6.11**.**
Assume that k<n or sβ₯n+2. Then
[TABLE]
Proof.
We use Propositions 2.22 and 2.23 to take the even part of equation (10).
β
6.4. Additional remarks
The above results suffice for proving our theorems, but for completeness we present the extension result for P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Sn), analogous to the result we proved for P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssββΞk(Sn) in Proposition 6.6.
Corollary 6.12**.**
Assume that k<n or sβ₯n+2. Then any form in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Sn) can be extended to a form in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1).
Proof.
Proposition 6.10 gives us that P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Sn)=uNβPsβnβ1βΞk(Sn). By definition, any form in uNβPsβnβ1βΞk(Sn) can be extended to a form in uNβPsβnβ1βΞk(Rn+1), and this space is equal to P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞk(Rn+1) by Proposition 6.2.
β
The conditions on k and s in Corollary 6.12 are necessary, as we see from the following example.
Example 6.13**.**
Let Ξ±=uNβvolSnβ. As discussed above, volSnβ can be extended to βRn+1βΞ½βP1βΞn(Rn+1), so Ξ± has an extension that is in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vssn+2βΞn(Rn+1). One can show that the degree n+2 is optimal due to the fact that any form in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vsssβΞn(Rn+1) must be divisible by uNβ and the fact that Ξ± is even. However, Ξ± itself is actually of lower degree. Indeed, Ξ± can be extended to Ξ±^=u2ββ―un+1βdu2ββ§β―β§dun+1β, putting Ξ± in P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vssnβΞn(Sn). We can check that Ξ±^ does in fact restrict to Ξ± by computing that Ξ½β§Ξ±^=uNβvolRn+1β, and hence βSnβΞ±^=uNβ. See also Examples 6.1 and 3.8.
This example also explains the seemingly paradoxical claim that P\vboxtoβ5.5757ptto0.0pt\hssΛΛ\hss\vssnβΞn(Sn)=uNβ(Rβ volSnβ) given by Proposition 6.10 when k=n and s=n: we have a differential form of degree n that is divisible by uNβ, even though uNβ has degree n+1. See also Example 2.7.
7. Conclusion
The transformation xiβ=ui2β induces a correspondence between differential forms on the simplex x1β+β―xn+1β=1 and differential forms the sphere u12β+β―+un+12β=1. We completely characterized the spaces of differential forms on the sphere corresponding to the PrβΞk(Tn), PrββΞk(Tn), PΛrβΞk(Tn), and PΛrββΞk(Tn) families. This correspondence gives an explanation for the isomorphisms PrβΞk(Tn)β PΛr+k+1ββΞnβk(Tn) and PrββΞk(Tn)β PΛr+kβΞnβk(Tn) used by Arnold, Falk, and Winther in their development of finite element exterior calculus: after the change of coordinates, both of these isomorphisms reveal themselves to be the Hodge star βSnβ followed by multiplication by the bubble function u1ββ―un+1β. Our result thus gives new isomorphism maps PrβΞk(Tn)βPΛr+k+1ββΞnβk(Tn) and PrββΞk(Tn)βPΛr+kβΞnβk(Tn) that are defined pointwise. As seen in our examples, evaluating these maps is a quick computation that does not require expressing the differential forms in terms of a basis for PrβΞk(Tn) or PrββΞk(Tn).
Another direction in which one could continue this work would be to consider parallelotope meshes instead of simplicial meshes. Finite element exterior calculus for parallelotope meshes [2, 3] works in much the same way as for simplicial meshes. Following the ideas in our current work, one could apply a coordinate transformation to convert a differential form on an n-dimensional parallelotope into an even differential form on an n-torus. Could such a coordinate transformation result in a better understanding of finite element spaces such as the serendipity elements?
8. Acknowledgments
I would like to thank Ari Stern, Douglas Arnold, Evan Gawlik, and the anonymous referees for their feedback on this project. I would also like to acknowledge the support of the AMSβSimons Travel Grant.
Appendix A Equivalence of definitions of PrββΞk(Tn)
The definition of PrββΞk(Tn) that we gave in Definition 2.9 is different from the definition given by Arnold, Walk, and Winther in [4]. In this appendix, we show that the two definitions are equivalent.
The definition in [4, Subsection 3.2] involves choosing an arbitrary point z in the simplex Tn and defining a vector field on Tn such that the vector based at xβTn points away from z with magnitude β£xβzβ£. Setting z to be the center of the simplex (n+11β,β¦,n+11β), one can check that this vector field is given by the formula
[TABLE]
for xβTn. Using this vector field, we can define PrββΞk(Tn).
Definition A.1**.**
For rβ₯1, Arnold, Falk, and Winther [4] define PrββΞk(Tn) to be
[TABLE]
where ΞΊ denotes the interior product with the vector field VΞΊβ.
Proposition A.2**.**
The definitions of PrββΞk(Tn) given by Definitions 2.9 and A.1 are equivalent.
Proof.
Let t=x1β+β―+xn+1β. With the equation
[TABLE]
we can extend the definition of VΞΊβ to all of Rn+1. Geometrically, VΞΊβ is the projection of the radial vector field X to the simplicies t=const. We can extend the definition of ΞΊ to ΞΊ:Ξk+1(Rn+1)βΞk(Rn+1) to be the interior product with the vector field VΞΊβ, and so we have
[TABLE]
Assume now that a satisfies Definition 2.9, so a can be extended to a differential form a^=iXβb^, where b^βPrβ1βΞk+1(Rn+1). Then
[TABLE]
Thus a=ΞΊb+c, where b and c are the restrictions of b^ and n+1tβiβtβb^ to Tn, respectively. By definition, bβPrβ1βΞk+1(Tn). Meanwhile, since t=1 on Tn, we know that c is also the restriction to Tn of n+11βiβtβb^, which is in Prβ1βΞk(Rn+1) because βt has polynomial degree zero. Thus cβPrβ1βΞk(Tn). We conclude that aβΞΊPrβ1βΞk+1(Tn)+Prβ1βΞk(Tn), so a satisfies Definition A.1.
Conversely, assume that a satisfies Definition A.1. We consider the two summands as separate cases.
If aβPrβ1βΞk(Tn), then let a^β²βPrβ1βΞk(Rn+1) be an arbitrary extension of a, and let
[TABLE]
Using equation (13), since the restriction of t to Tn is 1 and the restriction of dt to Tn is zero, the restriction of a^ to Tn is the same as the restriction of a^β², namely a. Meanwhile, using equation (11), we have that a^βiXβPrβ1βΞk+1(Rn+1) because dtβP0βΞ1(Rn+1), so dtβ§a^β²βPrβ1βΞk+1(Rn+1). Thus a satisfies Definition 2.9.
If a=ΞΊb for some bβPrβ1βΞk+1(Tn), then let b^β²βPrβ1βΞk+1(Rn+1) be an arbitrary extension of b. Set
[TABLE]
As before, observe from equation (16) that the restriction of b^ to Tn is the same as the restriction of b^β², namely b. Because VΞΊβ is tangent to Tn, we can then conclude that the restriction of ΞΊb^ to Tn is ΞΊb=a. Next, observe from equation (14) that iβtβb^=0. Consequently, we can set
[TABLE]
from which we see that a^ is an extension of a. Finally, observe that b^βPrβ1βΞk+1(Rn+1) because dt and βt both have polynomial degree zero. Thus a satisfies Definition 2.9, as desired.
β
Appendix B Vector space identities
In this appendix, we provide proofs of two vector space identities involving the Hodge star. I believe that these identities are well-known, but I have not been able to find published references for them.
B.1. The interior product, the exterior product, and the Hodge star
Let XβV, and let Ξ½βVβ be the dual vector corresponding to X with respect to the inner product. It is a standard result that the adjoint of iXβ is Ξ½β§ in the sense that
[TABLE]
where Ξ±^ββkVβ and Ξ²^βββkβ1Vβ.
This adjoint relationship has the following consequence for the Hodge star on ββVβ.
Proposition B.1**.**
Let V be an oriented inner product space. If XβV and Ξ½βVβ are dual to one another with respect to the inner product, then
[TABLE]
for all Ξ±^ββkVβ.
Proof.
Let dimV=n+1. For all Ξ²^βββnβkVβ, we have
[TABLE]
B.2. The Hodge star on a hyperplane
In this subsection, we consider a hyperplane H that is orthogonal to a unit covector Ξ½βVβ. The inner product on V induces an inner product on H, and the orientation on V along with the choice of unit conormal Ξ½ induces an orientation on H. Thus, in addition to a Hodge star operator ββVββββVβ, we also have a Hodge star operator ββHββββHβ. We denote these by βVβ and βHβ, respectively.
We show that βHβ can be computed from βVβ and Ξ½ as follows.
Proposition B.2**.**
With notation as above, let Ξ±ββkHβ be the restriction to H of Ξ±^ββkVβ. Then βHβΞ± is the restriction to H of
[TABLE]
Proof.
Let n=dimH, let Ξ²^βββnβkVβ denote βVβ(Ξ½β§Ξ±^), and let Ξ²ββnβkHβ be the restriction of Ξ²^β to H. We prove that βHβΞ±=Ξ² by verifying it on a basis for Ξ±^.
Let Ξ½,e1β,β¦,enβ be an oriented orthonormal basis for Vβ, so e1β,β¦,enβ is an oriented orthonormal basis for Hβ, and we have that volHβ=e1ββ§β―enβ and volVβ=Ξ½β§volHβ. For I={i1β<i2β<β―<ikβ}β{1,β¦,n}, let eIβ denote ei1βββ§β―β§eikββ.
If Ξ±^=Ξ½β§eIβ, then Ξ±=0 because the restriction of Ξ½ to H is zero. Also, Ξ½β§Ξ±^=0, so Ξ²^β and Ξ² are zero as well, as desired.
Now let Ξ±^=eIβ. By the definition of βVβ with respect to the oriented basis Ξ½,e1β,β¦,enβ, we see that Ξ²^β=βVβ(Ξ½β§eIβ) has the form Β±eJβ, where J={1,β¦,n}βI and the sign is chosen so that (Ξ½β§eIβ)β§(Β±eJβ)=volVβ. Since Ξ½β§(eIββ§Β±eJβ)=volVβ, we conclude that eIββ§Β±eJβ=volHβ. Thus, by the definition of βHβ with respect to the oriented basis e1β,β¦,enβ, we have βHβeIβ=Β±eJβ, that is, βHβΞ±=Ξ².
β
Bibliography12
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[1] D. N. Arnold , Spaces of finite element differential forms , in Analysis and numerics of partial differential equations, vol. 4 of Springer I Nd AM Ser., Springer, Milan, 2013, pp. 117β140.
2[2] D. N. Arnold and G. Awanou , The serendipity family of finite elements , Found. Comput. Math., 11 (2011), pp. 337β344.
3[3] D. N. Arnold, D. Boffi, and F. Bonizzoni , Finite element differential forms on curvilinear cubic meshes and their approximation properties , Numer. Math., 129 (2015), pp. 1β20.
4[4] D. N. Arnold, R. S. Falk, and R. Winther , Finite element exterior calculus, homological techniques, and applications , Acta Numer., 15 (2006), pp. 1β155.
5[5] , Geometric decompositions and local bases for spaces of finite element differential forms , Comput. Methods Appl. Mech. Engrg., 198 (2009), pp. 1660β1672.
6[6] Y. Berchenko-Kogan , Duality in finite element exterior calculus , 2018. http://arxiv.org/abs/1807.01161 .
7[7] F. Brezzi, J. Douglas, Jr., and L. D. Marini , Two families of mixed finite elements for second order elliptic problems , Numer. Math., 47 (1985), pp. 217β235.
8[8] R. S. Falk and R. Winther , Local bounded cochain projections , Math. Comp., 83 (2014), pp. 2631β2656.