Simple 3-D Visualization of Some Common Mathematical Minimal Surfaces using MATLAB

Size: px

Start display at page:

Download "Simple 3-D Visualization of Some Common Mathematical Minimal Surfaces using MATLAB"

Harold Miller
5 years ago
Views:

1 Simple 3-D Visualization of Some Common Mathematical Minimal Surfaces using MATLAB Nasiha Muna 1 and Albert E. Patterson 2 1 Department of Physics, Chemistry, and Mathematics, Alabama Agricultural and Mechanical (A&M) University, Normal, Alabama, 35762, USA. nasiha.muna@aamu.edu 2 Department of Industrial and Enterprise Systems Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois, 61801, USA. pttrsnv2@illinois.edu The MATLAB code for this work is available at: Surfaces_Visualization_Code. This report and its code are published under a CC-BY 4.0 International license, so it can be used and modified with proper attribution to the authors. Abstract This report presents a simple approach for visualizing some common mathematical minimal surfaces using MATLAB tools. The studied minimal surfaces were the gyroid, lidinoid, Schwarz P surface, Schwarz D surface, Scherk tower, helicoid surface, catenoid surface, and Mobius strip. For each, the mathematic definition is given along with a simple matlab code for it and an example 3-D surface plot generated from the given code. 1 Introduction The selected minimal surfaces to be discussed in this report were the gyroid, lidinoid, Schwarz P surface, Schwarz D surface, Scherk tower, helicoid surface, catenoid surface, and Mobius strip. These were selected, as they represent a suitable subset of the available minimal surfaces and serve as good cases to demonstrate the simple approach proposed in this report. Some very complex surfaces, such as Costa surfaces, are not covered here and will be discussed in future works. Most minimal surfaces are sometimes relatively difficult to plot in three dimensions, as they are all typically defined by implicit or parametric functions (sometimes in the complex domain). Fortunately, appropriate built-in MATLAB functions can be used to plot each of the surfaces relative to its defining equation. The gyroid and lidinoid functions are triply periodic minimal surfaces of genus 3 embedded in R 3 that possess no planar symmetry curves nor straight lines. They are the distinctive embedded surfaces of the associate families of the Schwarz P and H surfaces [1 3]. The gyroid function is typically defined as f gy = sin(x) cos(y) + sin(y) cos(z) + sin(z) cos(x) = 0 (1) while the lidinoid function is most often expressed as: f lid = 1 (sin(2x) cos(y) sin(z) + sin(2y) cos(z) sin(x) + sin(2z) cos(x) sin(y)) 2 1 (2) 2 (cos(2x) cos(2y) + cos(2y) cos(2z) + cos(2z)cos(2x)) = 0 The Schwarz minimal surfaces span many different kinds of forms, but the two most commonlyseen are the Schwarz Primitive ( Schwarz P ) and Schwarz Diamond ( Schwarz D ) surfaces [3 5]. The Schwarz P surface is the triply periodic minimal surface built from a skew hexagon inscribed in the edges of a regular octahedron. The Schwarz P surface encompasses sixfold junctions and is typically expressed as: 1

2 f sp = cos(x) + cos(y) + cos(z) (3) The Schwarz D surface is similar to the Schwarz P except that it contains six skew quadrilaterals from a skew hexagon inscribed in the edges of a cube. It has two intertwined corresponding mazes, each holding the figure of a version of the diamond arrangement; unlike Schwarz P, Schwarz D had fourfold junctions. Schwarz D is most often expressed as: f sd = sin(x) sin(y) sin(z) + sin(x) cos(y) cos(z)+ cos(x) sin(y) cos(z) + cos(x) cos(y) sin(z) = 0 (4) The Scherk surface is a 1-parameter singly-periodic, complete, properly-embedded surface in R 3 and can be modified in shape by defining new handles for it [6, 7]. Several forms exist, but the most common is the Scherk tower, which can be defined as: f sc = sinh(x) sinh(y) sin(z) (5) The helicoid surface is simply a non-algebraic surface that coils around an axis and is one of the parametric minimal surfaces discussed in this report. When used as a minimal surface, the edge of the surface is bounded by a circular helix and is therefore the only ruled minimal surface besides the plane. The helicoid is very commonly encountered in nature, with examples including seashells, horns, and polymer structures [8 10]. A good parametric equation to describe the surface is: x = u cos(t) y = u sin(t) where u and t are the parameters specified by the user. z = 2 3 t (6) The catenoid surface is another parametric minimal surface, which is an infinite, boundless, and locally isometric to the helicoid surface previously discussed. The equations are mathematically similar to the helicoid, with removal of the helical boundary and jointing of the ends [3, 11]. This involves the use of hyperbolic sines and cosines in the equation, which can be expressed as: x = cos(α) sinh(v) sin(u) + sin(α) cosh(v) cos(u) y = cos(α) sinh(v) cos(u) + sin(α) cosh(v) sin(u) z = u cos(α) + v sin(α) (7) where u and v are the parameters specified by the user and α = π/2 is the value for a catenoid (α = 0 corresponds to a helicoid) [12, 13]. Finally, the Mobius strip was examined. This is one of the most commonly-seen and known minimal surface in mathematics and the basic surface related to many other mathematical structures, including knots and Klein bottles [14 16]. A common parametric expression for the Mobius strip is: ( x = cos(u) 1 + v ( )) u 2 cos 2 ( y = sin(u) 1 + v ( )) u 2 cos 2 ) (8) z = v 2 sin ( u 2 where u and v are the parameters specified by the user. Technical Report 2

3 In each of the following subsections, the basic MATLAB code and a visualization of the 3-D surface are given. The first five (gyroid, lidinoid, Schwarz P, Schwarz D, and Scherk) are defined above by 3-D implicit equations, requiring the use of the fimplicit3 tool to plot. The others (helicoid, catenoid, and Mobius strip) are defined according to parametric equations and are therefore easier to plot using a standard fsurf built-in function. In all cases, the basic equation definition is given, followed by limits on the values of the variables or parameters, the mesh density setting, and the face color for the plot. All of these settings are infinitely-adjustable within the domain of the defining functions. It is recommended that MATLAB version 2017a or higher be used. 2 Minimal Surface Cases 2.1 Gyroid Technical Report 3

4 2.2 Lidinoid Figure 1: Lidinoid Figure 2: Lidinoid Technical Report 4

5 2.3 Schwarz P Surface Figure 3: Schwarz P Figure 4: Schwarz P Technical Report 5

6 2.4 Schwarz D Surface Figure 5: Schwarz D Figure 6: Schwarz D Technical Report 6

7 2.5 Scherk Surface Figure 7: Scherk Figure 8: Scherk Technical Report 7

8 2.6 Helicoid Surface Figure 9: Helicoid Figure 10: Helicoid Technical Report 8

9 2.7 Catenoid Surface Figure 11: Catenoid Figure 12: Catenoid Technical Report 9

2.8 Mobius Strip Figure 13: Mobius Figure 14: Mobius 3 Summary and Conclusions Minimal surfaces are very important in mathematics, art, and design.

10 2.8 Mobius Strip Figure 13: Mobius Figure 14: Mobius 3 Summary and Conclusions Minimal surfaces are very important in mathematics, art, and design. Various forms of minimal surfaces can be identified; visualizing minimal surfaces is important because it allows users to view and handle physical objects representing various natural and mathematical phenomena. By studying minimal surfaces, structures that are otherwise difficult to represent mathematically can be graphically displayed and understood. References [1] A. Weyhaupt, Deformations of the gyroid and lidinoid minimal surfaces, Pacific Journal of Mathematics, vol. 235, no. 1, pp , [2] S. N. Krivoshapko and V. N. Ivanov, Minimal surfaces, in Encyclopedia of Analytical Surfaces, pp , Springer International Publishing, Technical Report 10

11 [3] H. Karcher and K. Polthier, Construction of triply periodic minimal surfaces, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, vol. 354, no. 1715, pp , [4] H. Karcher, The triply periodic minimal surfaces of alan schoen and their constant mean curvature companions, Manuscripta Mathematica, vol. 64, no. 3, pp , [5] D. M. Anderson, H. T. Davis, J. C. C. Nitsche, and L. E. Scriven, Periodic surfaces of prescribed mean curvature, in Springer Proceedings in Physics, pp , Springer Berlin Heidelberg, [6] A. J. Y. Hancco, G. A. Lobos, and V. R. Batista, Explicit minimal scherk saddle towers of arbitrary even genera in r3, Publicacions Matemàtiques, vol. 58, pp , [7] J. Perez and M. Traizet, The classification of singly periodic minimal surfaces with genus zero and scherk-type ends, Transactions of the American Mathematical Society, vol. 359, pp , [8] C. Mamaloukas, On determination of developable and minimal surfaces, International Journal of Pure and Applied Mathematics, vol. 39, pp , [9] W. Meeks and H. Rosenberg, The uniqueness of the helicoid, Annals of Mathematics, vol. 161, no. 2, pp , [10] R. Ghafouri and R. Bruinsma, Helicoid to spiral ribbon transition, Physical Review Letters, vol. 94, no. 13, [11] A. Boudaoud, P. Patrcio, and M. B. Amar, The helicoid versus the catenoid: Geometrically induced bifurcations, Physical Review Letters, vol. 83, no. 19, pp , [12] A. Ogawa, Helicatenoid, The Mathematica Journal, vol. 2, no. 2, [13] H. Rosenberg, Deformations of complete minimal surfaces, Transactions of the American Mathematical Society, vol. 295, no. 2, pp , [14] G. Pradisi and A. Sagnotti, Open string orbifolds, Physics Letters B, vol. 216, no. 1-2, pp , [15] E. L. Starostin and G. H. M. van der Heijden, The shape of a möbius strip, Nature Materials, vol. 6, no. 8, pp , [16] V. Petresin and L.-P. Robert, The double möbius strip studies, Nexus Network Journal, vol. 4, no. 2, pp , Technical Report 11

arxiv:math/ v1 [math.dg] 1 Jul 1993

arxiv:math/ v1 [math.dg] 1 Jul 1993 APPEARED IN BULLETIN OF THE AMERICAN MATHEMATICAL SOCIETY Volume 29, Number 1, July 1993, Pages 77-84 ADDING HANDLES TO THE HELICOID DAVID HOFFMAN, FUSHENG WEI, AND HERMANN KARCHER arxiv:math/9307226v1