Advanced Search
Article Contents
Article Contents

Continuation methods for principal foliations of embedded surfaces

Abstract Full Text(HTML) Figure(17) Related Papers Cited by
  • Continuation methods are a well established tool for following equilibria and periodic orbits in dynamical systems as a parameter is varied. Properly formulated, they locate and classify bifurcations of these key components of phase portraits. Principal foliations of surfaces embedded in $ \mathbb{R}^3 $ resemble phase portraits of two dimensional vector fields, but they are not orientable. In the spirit of dynamical systems theory, Gutierrez and Sotomayor investigated qualitative geometric features that characterize structurally stable principal foliations and their bifurcations in one parameter families. This paper computes return maps and applies continuation methods to obtain new insight into bifurcations of principal foliations.

    Umbilics are the singularities of principal foliations and lines of curvature connecting umbilics are analogous to homoclinic and heteroclinic bifurcations of vector fields. Here, a continuation method tracks a periodic line of curvature in a family of surfaces that deforms an ellipsoid. One of the bifurcations of these periodic lines of curvature are connections between lemon umbilics. Differences between these bifurcations and analogous saddle connections in two dimensional vector fields are emphasized. A second case study tracks umbilics in a one parameter family of surfaces with continuation methods and locates their bifurcations using Taylor expansions in "Monge coordinates." Return maps that are generalized interval exchange maps of a circle are constructed for generic surfaces with no monstar umbilics.

    Mathematics Subject Classification: Primary: 58F15, 58F17; Secondary: 53C35.


    \begin{equation} \\ \end{equation}
  • 加载中
  • Figure 1.  Splitting of an umbilic connection in $ S_\lambda $ when $ \lambda = 0.001 $. The separatrices of umbilics $ u_l $ and $ u_r $ are drawn green and blue. The line of curvature heading right from point $ w_0 $ is drawn black, and its first five crossings of cross-section $ \Sigma $ are labelled $ w_1,w_2,w_3,w_4,w_5 $. This line of curvature stays to the right of the green separatrix of $ u_l $ between $ w_0 $ and $ w_1 $, then passes left of $ u_l $ and right of the blue separatrix of $ u_r $ into the channel between the two separatrices, intersecting $ \Sigma $ from left to right at $ w_2 $. From $ w_2 $ it turns around $ u_r $ to the left of the separatrix of $ u_l $, next crossing $ \Sigma $ from right to left at $ w_3 $. From $ w_3 $ it turns around $ u_l $ with the separatrix of $ u_r $ on its right and the next turn of the separatrix of $ u_l $ on its left. It then crosses $ \Sigma $ from left to right at $ w_4 $ and finally crosses $ \Sigma $ again from right to left at $ w_5 $. Note that while the line of curvature originally headed right from $ \Sigma $ at $ w_0 $, it crosses $ \Sigma $ in the opposite direction at $ w_5 $

    Figure 2.  Graph of the return map $ \sigma $ and its fifth iterate for the surface $ S_{0.1857} $. The cross-section $ \Sigma $ is parameterized as $ [0,\sqrt{3}\cos(\theta),\sqrt{5}\sin(\theta)],\, \theta \in [-\pi,\pi] $. Values of the maps are displayed on a mesh of $ 1000 $ points

    Figure 3.  Graph of $ \sigma^5(\theta) - \theta $. Note that $ \sigma^5 $ has fixed points and breaks

    Figure 4.  The derivative of $ \sigma $, showing four jump discontinuities

    Figure 5.  Two closed lines of curvature on the surface $ S_{0.1857} $. The umbilics are drawn as red dots

    Figure 6.  Values of $ (\theta,\lambda) $ at which the return maps $ \sigma_{\lambda} $ of the surfaces $ S_{\lambda} $ have periodic orbits of period 5

    Figure 7.  The geometry of closed lines of curvature near an umbilic connection. The blue curve is an umbilic connection computed at $ \lambda = 0.185952861855609 $. The green curve $ \gamma $ is a closed line of curvature on $ S_{0.1857} $. Its intersections $ s_j $ with cross-section $ x = 0 $ (light black dotted curve) are marked by black dots. Umbilics are red dots. The point $ s_0 $ is the initial point of $ \gamma $, and its period $ 5 $ trajectory for the return map $ \sigma $ is $ s_{0},s_{2},s_{4},s_{6},s_{8},s_{10} = s_{0} $. Note that between $ s_2 $ and $ s_3 $, the line of curvature comes close to the umbilic $ u_1 $, following its separatrix in one direction before the turn and in the other direction after the turn. The $ \sigma $ trajectory of $ s_5 $ is $ s_{5},s_{7},s_{9},s_{1},s_{3} $, lying on the same line of curvature but following it with the opposite orientation

    Figure 8.  Values of $ (\theta,\lambda) $ at which the return maps $ \sigma_{\lambda} $ of the surfaces $ S_{\lambda} $ have periodic orbits of period 5 for $ \theta \in [-2.75,-2.65] $

    Figure 9.  The surface $ S_{0.24} $ defined by Equation (3) with $ b = 0.24 $. The surface is plotted in a tan color. The six umbilic points are red dots and segments of their separatrices up to intersections with the cross-section $ x = 2y $ are drawn as heavy blue curves. The cross-section is a thin black curve. A long segment of a single line of curvature is plotted as a brown curve. Note that the separatrix of the lowest umbilic in the figure passes close enough to another lemon umbilic that the gap between them cannot be seen at this scale

    Figure 10.  The $ z $ coordinate of umbilic points of the surfaces $ S_b $ plotted as a function of $ b $. Turning points on the curves occur at parameters where star and monstar umbilics merge and disappear, analogous to a saddle-node bifurcation of vector fields. Note that there are ranges of parameters with four, six and eight umbilic points

    Figure 11.  The $ z $ coordinate of umbilic points of the surfaces $ S_b $ plotted as a function of $ b $. Turning points on the curves, marked with blue crosses, occur at parameters where star and monstar umbilics merge and disappear. Large green dots mark monstar-lemon transitions

    Figure 12.  The return map $ \sigma $ of the surface $ S_{0.24} $. The cross-section $ \Sigma $ is a closed curve parameterized as $ (2r(\theta)\cos(\theta),r(\theta)\cos(\theta),r(\theta)\sin(\theta)) $ where $ r(\theta) $ is computed numerically so that $ \Sigma $ lies on $ S_{0.24} $. Lines of curvature starting at 1000 points on $ \Sigma $ are computed to their second return with $ \Sigma $. There are three discontinuities at intersections of the star umbilic separatrices with $ \Sigma $. Between each pair of discontinuities, $ \sigma $ is increasing. It has breaks (jumps in its first derivative) at intersections of the separatrices of lemon umbilics with $ \Sigma $. Apart from its discontinuities, $ \sigma $ is bijective

    Figure 13.  Successive returns of a long line of curvature on the surface $ S_{0.24} $. Over $ 42,000 $ returns are plotted as blue dots. If the line of curvature is dense in $ S_{0.24} $, its succesive returns would fill the graph of the return map $ \sigma $ densely. Here, quite large gaps remain despite the large number of returns. Nonetheless, the line of curvature does not seem to approach one that is closed. If that were the case, the successive returns would approach a discrete set of points in the intersection of the line of curvature and the cross-section. Here, a plot of only the last 1000 successive returns of this line of curvature produces a graph visually indistinguishable from this one

    Figure 14.  Principal foliation of $ S_{0.2233} $. Umbilic separatrices are heavy blue curves

    Figure 15.  Region near the star and monstar of $ S_{0.2233} $

    Figure 16.  The return map $ \sigma $ of $ S_{0.2233} $ to the cross-section $ x = 2y $. In comparison to the return map of $ S_{0.24} $ shown in figure 12, there are gaps in the domain and {image of $ \sigma $.} These come from lines of curvature that tend to the monstar umbilic and do not return to the cross-section

    Figure 17.  A (near) umbilic connection (blue) between a star umbilic (large red dot) and a lemon umbilic for $ S_{0.2447609874} $. The second and third separatrices tend to the same closed line of curvature, plotted green

  • [1] M. V. Berry and J. H. Hannay, Umbilic points on gaussian random surfaces, J. Physics A: Mathematical and General, 10 (1977), 1809-1821. 
    [2] J. W. Bruce and D. L. Fidal, On binary differential equations and umbilics, Proc. Roy. Soc. Edinburgh Sect. A, 111 (1989), 147-168.  doi: 10.1017/S0308210500025087.
    [3] F. CazalsJ.-C. FaugèreM. Pouget and F. Rouillier, The implicit structure of ridges of a smooth parametric surface, Comput. Aided Geom. Design, 23 (2006), 582-598.  doi: 10.1016/j.cagd.2006.04.002.
    [4] H. Dankowicz and F. Schilder, Recipes for Continuation, Computational Science & Engineering, 11. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2013. doi: 10.1137/1.9781611972573.
    [5] G. Darboux, Leç cons Sur La Théorie Générale des Surfaces et Les Applications Géométriques du Calcul Infinitésimal, Réimpression de la premiéreédition de 1896. Chelsea Publishing Co., Bronx, N. Y., 1972.
    [6] W. de Melo and S. van Strien, One-Dimensional Dynamics, Springer-Verlag, Berlin, 1993. doi: 10.1007/978-3-642-78043-1.
    [7] A. DhoogeW. Govaerts and Y. A. Kuznetsov, MATCONT: A MATLAB package for numerical bifurcation analysis of ODEs, ACM Trans. Math. Software, 29 (2003), 141-164.  doi: 10.1145/779359.779362.
    [8] M. P. do Carmo, Differential Geometry of Curves & Surfaces, Dover Publications, Inc., Mineola, NY, 2016, Revised & updated second edition of [MR0394451].
    [9] E. DoedelH. B. Keller and J.-P. Kernévez, Numerical analysis and control of bifurcation problems. I. Bifurcation in finite dimensions, Internat. J. Bifur. Chaos Appl. Sci. Engrg., 1 (1991), 493-520.  doi: 10.1142/S0218127491000397.
    [10] C. Dupin, Développements de Géométrie, Corcier, Paris, 1813.
    [11] W. J. F. Govaerts, Numerical Methods for Bifurcations of Dynamical Equilibria, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2000. doi: 10.1137/1.9780898719543.
    [12] J. Guckenheimer, Bifurcations of dynamical systems, In Dynamical systems (C.I.M.E. Summer School, Bressanone, 1978), Progr. Math., Birkhäuser, Boston, Mass., 8 (1980), 115–231.
    [13] J. Guckenheimer, Dense lines of curvature on convex surfaces, Proc. Amer. Math. Soc., 148 (2020), 3537-3549.  doi: 10.1090/proc/14981.
    [14] C. Gutiérrez and J. Sotomayor, Closed principal lines and bifurcation, Bol. Soc. Brasil. Mat., 17 (1986), 1-19.  doi: 10.1007/BF02585473.
    [15] C. Gutiérrez and J. Sotomayor, Periodic lines of curvature bifurcating from Darbouxian umbilical connections, In Bifurcations of Planar Vector Fields (Luminy, 1989), Lecture Notes in Math., Springer, Berlin, 1455 (1990), 196–229. doi: 10.1007/BFb0085394.
    [16] E. Hairer, S. P. Nørsett and G. Wanner, Solving Ordinary Differential Equations, 2$^nd$ edition, Springer Series in Computational Mathematics, 8. Springer-Verlag, Berlin, 1993.
    [17] M. W. HirschS. Smale and  R. L. DevaneyDifferential Equations, Dynamical Systems, and An Introduction to Chaos, 3$^{rd}$ edition, Elsevier/Academic Press, Amsterdam, 2013.  doi: 10.1016/B978-0-12-382010-5.00001-4.
    [18] Y. Katznelson and D. Ornstein, The differentiability of the conjugation of certain diffeomorphisms of the circle, Ergodic Theory Dynam. Systems, 9 (1989), 643-680.  doi: 10.1017/S0143385700005277.
    [19] B. Krauskopf, H. M. Osinga and J. Galán-Vioque (eds.), Numerical Continuation Methods for Dynamical Systems, Understanding Complex Systems, Springer, Dordrecht, 2007, Path following and boundary value problems, Dedicated to Eusebius J. Doedel for his 60th birthday. doi: 10.1007/978-1-4020-6356-5.
    [20] Y. A. Kuznetsov, Elements of Applied Bifurcation Theory, 3$^{rd}$ edition, Springer-Verlag, New York, 2004. doi: 10.1007/978-1-4757-3978-7.
    [21] S. MarmiP. Moussa and J.-C. Yoccoz, Linearization of generalized interval exchange maps, Ann. of Math. (2), 176 (2012), 1583-1646.  doi: 10.4007/annals.2012.176.3.5.
    [22] M. MartensS. van StrienW. de Melo and P. Mendes, On Cherry flows, Ergodic Theory Dynam. Systems, 10 (1990), 531-554.  doi: 10.1017/S0143385700005733.
    [23] G. Monge, Sur les lignes de courbure de la surface de l'ellipsoide, In Application de l'Analyse á la Géometrie (5th edition), Bachelier, Paris, 1796,139–160.
    [24] M. M. Peixoto, Structural stability on two-dimensional manifolds, Topology, 1 (1962), 101-120.  doi: 10.1016/0040-9383(65)90018-2.
    [25] I. R. PorteousGeometric Differentiation, For the Intelligence of Curves and Surfaces, 2$^{nd}$ edition, Cambridge University Press, Cambridge, 2001. 
    [26] Y. G. Sinaǐ and K. M. Khanin, Smoothness of conjugacies of diffeomorphisms of the circle with rotations, Russian Math. Surveys, 44 (1989), 69-99.  doi: 10.1070/RM1989v044n01ABEH002008.
    [27] J. Sotomayor and C. Gutiérrez, Structurally stable configurations of lines of principal curvature, In Bifurcation, Ergodic Theory and Applications (Dijon, 1981), Astérisque, Soc. Math. France, Paris, 98 (1982), 195–215.
    [28] J. Sotomayor and C. Gutiérrez, Configurations of lines of principal curvature and their bifurcations, In Colloquium on Dynamical Systems (Guanajuato, 1983), Aportaciones Mat., Soc. Mat. Mexicana, México, 1 (1985), 115–126.
    [29] J. Sotomayor and R. Garcia, Lines of curvature on surfaces, historical comments and recent developments, São Paulo J. Math. Sci., 2 (2008), 99-143.  doi: 10.11606/issn.2316-9028.v2i1p99-143.
    [30] J. Sotomayor and C. Gutiérrez, Structurally Stable Configurations of Lines of Curvature and Umbilic Points on Surfaces, Instituto de Matemática y Ciencias Afines, IMCA, Lima; Universidad Nacional de Ingenieria, Instituto de Matemáticas Puras y Aplicadas, Lima, 1998.
    [31] M. Viana, Ergodic theory of interval exchange maps, Rev. Mat. Complut., 19 (2006), 7-100.  doi: 10.5209/rev_REMA.2006.v19.n1.16621.
    [32] E. Zeeman, The umbilic bracelet and the double-cusp catastrophe, In Structural Stability, the Theory of Catastrophes, and Applications in the Sciences, Lecture Notes in Math., 525 (1976), 328–366.
  • 加载中



Article Metrics

HTML views(1299) PDF downloads(341) Cited by(0)

Access History

Other Articles By Authors



    DownLoad:  Full-Size Img  PowerPoint