High-Order Method with Moving Frames to Compute the Covariant Derivatives of Vectors on General 2D Curved Surfaces

doi:10.1007/s42967-022-00225-x

Abstract

Abstract: The covariant derivative is a generalization of differentiating vectors. The Euclidean derivative is a special case of the covariant derivative in Euclidean space. The covariant derivative gathers broad attention, particularly when computing vector derivatives on curved surfaces and volumes in various applications. Covariant derivatives have been computed using the metric tensor from the analytically known curved axes. However, deriving the global axis for the domain has been mathematically and computationally challenging for an arbitrary two-dimensional (2D) surface. Consequently, computing the covariant derivative has been difficult or even impossible. A novel high-order numerical scheme is proposed for computing the covariant derivative on any 2D curved surface. A set of orthonormal vectors, known as moving frames, expand vectors to compute accurately covariant derivatives on 2D curved surfaces. The proposed scheme does not require the construction of curved axes for the metric tensor or the Christoffel symbols. The connectivity given by the Christoffel symbols is equivalently provided by the attitude matrix of orthonormal moving frames. Consequently, the proposed scheme can be extended to the general 2D curved surface. As an application, the Helmholtz-Hodge decomposition is considered for a realistic atrium and a bunny.

Key words: Covariant derivative, Moving frames, Curved surface, Helmholtz?Hodge decomposition

CLC Number:

Sehun Chun. High-Order Method with Moving Frames to Compute the Covariant Derivatives of Vectors on General 2D Curved Surfaces[J]. Communications on Applied Mathematics and Computation, 2023, 5(4): 1534-1563.

TrendMD

References

1. Adams, J.C., Swarztrauber, P.N.: Spherepack 2.0: a model development facility (No. NCAR/TN- 436+STR). https://opensky.ucar.edu/islandora/object/technotes:297
2. Arfken, G.B., Weber, H.J., Harris, F.E.: Mathematical Methods for Physicists, 7th edn. Elsevier India, India (2012)
3. Arnold, D.N., Falk, R.S., Winther, R.: Finite element exterior calculus from Hodge theory to numerical stability. Bull. Am. Math. Soc. 47(2), 281-354 (2010)
4. Batard, T., Bertalmio, M.: On covariant derivatives and their applications to image regularization. SIAM J. Imaging Sci. 7(4), 2393-2422 (2014)
5. Bernard, P.-E., Remacle, J.-F., Comblen, R., Legat, V., Hillewaert, K.: High-order discontinuous Galerkin schemes on general 2D manifolds applied to shallow water equations. J. Comput. Phys. 228(17), 6514-6535 (2009)
6. Bertalmio, M., Cheng, L., Osher, S., Sapiro, G.: Variational problems and partial differential equations on implicit surfaces. J. Comput. Phys. 174(2), 759-780 (2001)
7. Beyer, F., Daszuta, B., Frauendiener, J., Whale, B.: Numerical evolutions of fields on the 2-sphere using a spectral method based on spin-weighted spherical harmonics. Class. Quantum Gravity 31, 075019 (2014)
8. Bhatia, H., Norgard, G., Pascucci, V., Bremer, P.-T.: The Helmholtz-Hodge decomposition—a survey. IEEE. Trans. Vis. Compt. Graph. 19(8), 1386-1404 (2013)
9. Cantwell, C.D., Moxey, D., Comerford, A., Bolis, A., Rocco, G., Mengaldo, G., de Grazia, D., Yakovlev, S., Lombard, J.-E., Ekelschot, D., Jordi, B., Mohamied, Y., Eskilsson, C., Nelson, B., Vos, P., Biotto, C., Kirby, R.M., Sherwin, S.J.: Nektar++: an open-source spectral/hp element framework. Comput. Phys. Commun. 192, 205-219 (2015)
10. Cantwell, C.D., Yakovlev, S., Kirby, R.M., Peters, N.S., Sherwin, S.J.: High-order spectral/hp element discretisation for reaction-diffusion problems on surfaces: application to cardiac electrophysiology. J. Comput. Phys. 257, 813-829 (2014)
11. Cartan, É.: Geometry of Riemannian Spaces. Math. Sci. Press, Brookline (2001)
12. Cartan, É.: Riemannian Geometry in an Orthogonal Frame. World Scientific Pub. Co., Inc., Hackensack (2002)
13. Chun, S.: Method of moving frames to solve conservation laws on curved surfaces. J. Sci. Comput. 53(2), 268-294 (2012)
14. Chun, S.: Method of moving frames to solve (an)isotropic diffusion equations on curved surfaces. J. Sci. Comput. 59(3), 626-666 (2013)
15. Chun, S.: A mathematical model of the unidirectional block caused by the pulmonary veins for anatomically-induced atrial reentry. J. Biol. Phys. 40(3), 219-258 (2014)
16. Chun, S.: Method of moving frames to solve the time-dependent Maxwell’s equations on anisotropic curved surfaces: applications to invisible cloak and ELF propagation. J. Comput. Phys. 340, 85-104 (2017)
17. Chun, S., Eskilsson, C.: Method of moving frames to solve the shallow water equations on arbitrary rotating curved surfaces. J. Comput. Phys. 333, 1-23 (2017)
18. Chun, S., Marcon, J., Peiro, J., Sherwin, S.J.: Reducing errors caused by geometrical inaccuracy to solve partial differential equations with moving frames on curvilinear domain. Comput. Methodol. Appl. Mech. Eng. 398(1), 115261 (2022)
19. Chun, S., Oh, T.: Divergence/connection preservation scheme in the curvilinear domain with a small geometric approximation error. J. Sci. Comput. 92(1):15 (2022)
20. Coté, J.: A Lagrange multiplier approach for the metric terms of semi-Lagrangian models on the sphere. Q. J. R. Meteorol. Soc. 114(483), 1347-1352 (1988)
21. Cotter, C.J., Thuburn, J.: A finite element exterior calculus framework for the rotating shallow water equations. J. Comput. Phys. 257, 1506-1526 (2014)
22. Deckelnick, K., Dziuk, G., Elliott, C.M., Heine, C.-J.: An h-narrow band finite-element method for elliptic equations on implicit surfaces. IMA J. Numer. Anal. 30, 351-376 (2010)
23. Giraldo, F.X., Hesthaven, J.S., Warburton, T.: Nodal high-order discontinuous Galerkin methods for the spherical shallow water equations. J. Comput. Phys. 181(2), 499-525 (2002)
24. Heinze, T., Hense, A.: The shallow water equations on the sphere and their Lagrange-Galerkin-solution. Meteorol. Atmos. Phys. 81, 129-137 (2002)
25. Läuter, M., Giraldo, F.X., Handorf, D., Dethloff, K.: A discontinuous Galerkin method for the shallow water equations in spherical triangular coordinates. J. Comput. Phys. 227, 10226-10242 (2008)
26. Li, S., Xiao, F.: A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid. J. Comput. Phys. 229, 1774-1796 (2010)
27. Liu, B., Tong, Y., Goes, F., Desbrun, M.: Discrete connection and covariant derivative for vector field analysis and design. ACM. Trans. Graph. 35(3), 1-17 (2016)
28. Meyer, M., Desbrun, M., Schröder, P., Barr, A.H.: Discrete differential-geometry operators for triangulated 2-manifolds. VisMath ’02 Proceedings. Springer, Berlin/Heidelberg (2002)
29. Misner, C.W., Thorne, K.S., Wheeler, J.A.: Gravitation. Freeman, Boonton (2021)
30. Moore, T.A.: A General Relativity Workbook. University Science Books, Mill Valley (2013)
31. Needham, T.: Visual Differential Geometry and Forms. Princeton University Press, New York (1973)
32. Nestler, M., Nitschke, I., Voigt, A.: A finite element approach for vector- and tensor-valued surface PDEs. J. Comput. Phys. 389, 48-61 (2019)
33. Novak, I.L., Gao, F., Choi, Y.-S., Resasco, D., Schaff, J.C., Slepchenko, B.M.: Diffusion on a curved surface coupled to diffusion in the volume: application to cell biology. J. Comput. Phys. 226(2), 1271- 1290 (2007)
34. O’Neil, B.: Elementary Differential Geometry. Academic Press, Amsterdam (1997)
35. Pudykiewicz, J.A.: On numerical solution of the shallow water equations with chemical reactions on icosahedral geodesic grid. J. Comput. Phys. 230, 1956-1991 (2011)
36. Rossmanith, J.A., Bale, D.S., LeVeque, R.J.: A wave propagation algorithm for hyperbolic systems on curved manifolds. J. Comput. Phys. 199, 631-662 (2004)
37. Sherwin, S.J., Peiró, J.: Mesh generation in curvilinear domains using high-order elements. Int. J. Numer. Methods Eng. 53, 207-223 (2002)
38. Swarztrauber, P.N., Williamson, D.L., Drake, J.B.: The Cartesian method for solving partial differential equations in spherical geometry. Dyn. Atmos. Oceans 27, 679-706 (1997)
39. Taylor, M.A., Fournier, A.: A compatible and conservative spectral element method on unstructured grids. J. Comput. Phys. 229, 5879-5895 (2010)
40. Turner, M., Moxey, D., Sherwin, S.J., Peiró, J.: Automatic generation of 3D unstructured high-order curvilinear meshes. ECCOMAS Congress 2016, 428-443 (2016)
41. Williamson, D.L., Drake, J.B., Hack, J.J., Jakob, R., Swarztrauber, P.N.: A standard test set for numerical approximations to the shallow water equations in spherical geometry. J. Comput. Phys. 102(1), 211-224 (1992)
42. Xu, G.: Convergence of discrete Laplace-Beltrami operators over surfaces. Comput. Math. Appl. 48, 347-360 (2004)
43. Yang, J., Li, Y., Kim, J.: A practical finite difference scheme for the Navier-Stokes equation on curved surfaces in R³. J. Comput. Phys. 411(15), 109430 (2020)