Comparison of discrete Hodge star operators for surfaces

We investigate the performance of various discrete Hodge star operators for discrete exterior calculus (DEC) using circumcentric and barycentric dual meshes. The performance is evaluated through the DEC solution of Darcy and incompressible Navier-Stokes flows over surfaces. While the circumcentric Hodge operators may be favorable due to their diagonal structure, the barycentric (geometric) and the Galerkin Hodge operators have the advantage of admitting arbitrary simplicial meshes. Numerical experiments reveal that the barycentric and the Galerkin Hodge operators retain the numerical convergence order attained through the circumcentric (diagonal) Hodge operators. Furthermore, when the barycentric or the Galerkin Hodge operators are employed, a super-convergence behavior is observed for the incompressible flow solution over unstructured simplicial surface meshes generated by successive subdivision of coarser meshes. Insofar as the computational cost is concerned, the Darcy flow solutions exhibit a moderate increase in the solution time when using the barycentric or the Galerkin Hodge operators due to a modest decrease in the linear system sparsity. On the other hand, for the incompressible flow simulations, both the solution time and the linear system sparsity do not change for either the circumcentric or the barycentric and the Galerkin Hodge operators. Various definitions for discrete Hodge star operators are investigated.Darcy and incompressible Navier-Stokes flows are used for numerical experiments.Barycentric Hodge star reproduces the convergence rate of the circumcentric Hodge.The difference in the computational cost is generally minor.

[1]  J. Marsden,et al.  Structure-preserving discretization of incompressible fluids , 2009, 0912.3989.

[2]  Yiying Tong,et al.  Stable, circulation-preserving, simplicial fluids , 2006, SIGGRAPH Courses.

[3]  B. Auchmann,et al.  A geometrically defined discrete hodge operator on simplicial cells , 2006, IEEE Transactions on Magnetics.

[4]  Maks Ovsjanikov,et al.  Discrete Derivatives of Vector Fields on Surfaces -- An Operator Approach , 2015, ACM Trans. Graph..

[5]  H. Whitney.  Geometric Integration Theory , 1957 .

[6]  Keenan Crane,et al.  Energy-preserving integrators for fluid animation , 2009, SIGGRAPH 2009.

[7]  Yiying Tong,et al.  Stable, circulation-preserving, simplicial fluids , 2007, TOGS.

[8]  Yiying Tong,et al.  Model-reduced variational fluid simulation , 2015, ACM Trans. Graph..

[9]  Evan S. Gawlik,et al.  Geometric, variational discretization of continuum theories , 2010, 1010.4851.

[10]  Anil N. Hirani,et al.  Numerical Method for Darcy Flow Derived Using Discrete Exterior Calculus , 2008, ArXiv.

[11]  Chandrajit L. Bajaj,et al.  Dual formulations of mixed finite element methods with applications , 2010, Comput. Aided Des..

[12]  Anil N. Hirani,et al.  Delaunay Hodge star , 2012, Comput. Aided Des..

[13]  Maks Ovsjanikov,et al.  Functional Fluids on Surfaces , 2014 .

[14]  R. Nicolaides.  Direct discretization of planar div-curl problems , 1992 .

[15]  L. Kettunen,et al.  Geometric interpretation of discrete approaches to solving magnetostatic problems , 2004, IEEE Transactions on Magnetics.

[16]  Mathieu Desbrun,et al.  Barycentric coordinates for convex sets , 2007, Adv. Comput. Math..

[17]  D. Arnold,et al.  Finite element exterior calculus: From hodge theory to numerical stability , 2009, 0906.4325.

[18]  Anil N. Hirani,et al.  Discrete exterior calculus discretization of incompressible Navier-Stokes equations over surface simplicial meshes , 2015, J. Comput. Phys..

[19]  Anil N. Hirani,et al.  Discrete exterior calculus for variational problems in computer vision and graphics , 2003, 42nd IEEE International Conference on Decision and Control (IEEE Cat. No.03CH37475).

[20]  Anil N. Hirani,et al.  Discrete exterior calculus , 2005, math/0508341.

[21]  A. Bossavit.  Computational Electromagnetism: Variational Formulations, Complementarity, Edge Elements , 1997 .

[22]  Peter Schröder,et al.  An algorithm for the construction of intrinsic delaunay triangulations with applications to digital geometry processing , 2007, Computing.