Next Article in Journal
Robust-MBDL: A Robust Multi-Branch Deep-Learning-Based Model for Remaining Useful Life Prediction of Rotating Machines
Previous Article in Journal
Robustness of Real-World Networks after Weight Thresholding with Strong Link Removal
Previous Article in Special Issue
On the Structure of SO(3): Trace and Canonical Decompositions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameterizations of Delaunay Surfaces from Scratch

by
Clementina D. Mladenova
1 and
Ivaïlo M. Mladenov
2,*
1
Institute of Mechanics, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Bl. 4, 1113 Sofia, Bulgaria
2
Institute for Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, Tsarigradsko Chaussee 72, 1784 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(10), 1570; https://doi.org/10.3390/math12101570
Submission received: 30 April 2024 / Revised: 14 May 2024 / Accepted: 15 May 2024 / Published: 17 May 2024
(This article belongs to the Special Issue Differentiable Manifolds and Geometric Structures)

Abstract

:
Starting with the most fundamental differential-geometric principles we derive here new explicit parameterizations of the Delaunay surfaces of revolution which depend on two real parameters with fixed ranges. Besides, we have proved that these parameters have very clear geometrical meanings. The first one is responsible for the size of the surface under consideration and the second one specifies its shape. Depending on the concrete ranges of these parameters any of the Delaunay surfaces which is neither a cylinder nor the plane is classified unambiguously either as a first or a second kind Delaunay surface. According to this classification spheres are Delaunay surfaces of first kind while the catenoids are Delaunay surfaces of second kind. Geometry of both classes is established meaning that the coefficients of their fundamental forms are found in explicit form.

1. Introduction

By its very definition (cf. any book on differential geometry), the mean curvature H at some point of a given surface S in the three-dimensional Euclidean space is the arithmetic mean of the two principal curvatures (maximal and minimal ones) at that point.
In the most intuitive representation (the so-called Monge parameterization) of S , which relies on the Cartesian coordinates ( x , y , z ) in the space R 3 , the surface is given as a graph defined over some domain U of the plane, i.e.,
S = { ( x , y , z ) , z = z ( x , y ) , ( x , y ) U R 2 } .
Then, according to the formula available in the classical differential geometry (see [1,2,3,4]), the mean curvature is expressed in the form
H ( x , y ) = r ( 1 + q 2 ) 2 p q s + t ( 1 + p 2 ) ( 1 + p 2 + q 2 ) 3 / 2
in which
p = z x ( x , y ) , q = z y ( x , y ) r = 2 z 2 x ( x , y ) , s = 2 z x y ( x , y ) , t = 2 z 2 y ( x , y ) ·
Unfortunately, solving the equation
H ( x , y ) = c c o n s t
turns out to be a formidably complicated task!
In this situation the French astronomer and mathematician Charles-Eugène Delaunay immediately abandoned the general problem stated above and has restricted his considerations to axially-symmetric surfaces. In doing this he undertakes a genuinely geometric approach by showing that all profile curves (meridional sections) of the CMC surfaces of revolution (except for the plane) can be obtained as roulettes of the conic sections [5] (named after him Delaunay roulettes).
In fact, these Delaunay roulettes are nothing else, but traces of the foci of the non-degenerate conic sections when they roll (without slipping) along the axis of revolution (Figure 1).
In parallel, in a note to the same Delaunay’s paper, Sturm [6] had proposed another approach of variational nature in which these surfaces are characterized as surfaces of revolution having a minimal lateral surface area for a given (fixed in advance) volume.
More recent accounts on the subject can found in the paper by Eells [7] and in the handbook by Krivoshapko-Ivanov [8] (pp. 155–157).
Let us recall that in the case of a rotational CMC surface the principle curvatures are known as a parallel κ π and a meridional κ μ curvatures and therefore one should solve the equation
H : = κ π + κ μ 2 = c
where c R is an arbitrary constant.
This has been done already by Delaunay in 1841, who has found also the first analytical (integral) representation of the constant mean curvatures (CMC) surfaces of revolution [5]. The complete list of such surfaces consists of planes, cylinders, sequences of spheres, catenoids, unduloids, and nodoids (see Figure 2), known nowadays just as Delaunay surfaces.
In his paper Delaunay has introduced the constant mean curvature of the surfaces in question in the form
H = 1 2 1 R π + 1 R μ = 1 2 a
where by definition R π and R μ are the parallel and meridional curvature radii respectively and at the beginning the parameter a was just an arbitrary constant. Latter on this parameter has been recognized as the major semi-axis of the rolling conic.
In the variational approach respectively, a initially has played the role of the Lagrangian multiplier, but at the end it was also identified as a length of the semi-axis. In both considerations, during the processes of integrations another integration constant appears which was recognized readily as the length of the minor semi-axis of the rolling conic.
Viewed in this way, the situation is a little bit strange: the profile curve depends on two parameters, but the mean curvature of the generated rotational surface depends only on one of them!
In what follows the problem will be considered from a different viewpoint and the first step in this direction will be to rewrite the expression for the mean curvature H in another form.
This possibility is based on the obvious observation that any profile curve has a point which is situated at a maximal distance, say R from the symmetry axis. It is clear also that at such point the tangent to the curve is parallel to the symmetry axis of the surface and that its normal is parallel to the O X axis. From here, one can conclude that the radius of the parallel curvature R π is just R and that along the circle (which will be further referred as equator) with this radius the meridional curvature R μ is also constant (due to the basic assumption (4) about the surface). The ratio of these radii, R π ( R ) / R μ ( R ) will be denoted further on by ν , i.e.,
R π ( R ) R = ν R μ ( R ) .
So, we have in place of (5) the equation
H = 1 2 1 R π ( R ) + 1 R μ ( R ) = 1 2 1 R + ν R = ν + 1 2 R ·
The bottom line of all above is that the problem of describing the Delaunay surfaces with mean curvature specified by (4) transforms into that of finding the axially-symmetric surfaces which mean curvature is fixed by Formula (7).
Remark 1. 
The above reasoning is valid mutatis mutandis for the minimal distance r as well!
The actual organization of the paper is as follows. The main results are presented in the next Section 2 in which by relying on the first differential-geometric principles the fundamental equations behind the Delaunay surfaces are derived. Their solutions are formulated as two Propositions associated with the two different kinds of Delaunay surfaces and later on they are unified as Theorem 1 describing in explicit form all Delaunay surfaces. The above results are illustrated by various graphics which aims are to make them more clear.
As it is well known the intrinsic geometry of the surfaces in the Euclidean space is encoded into their fundamental forms which have been found explicitly in Section 3.
The Section 4 points out the importance of these surfaces for other branches of the sciences outside the mathematics. Finally, we have added an Appendix A in which we have collected all known (to the authors) explicit parameterizations of the Delaunay surfaces that can be useful for theoretical and applied purposes in the future.

2. The New Parameterizations

The starting point of our considerations will be some basic facts in the classical differential geometry which can be easily located in any of the above cited books [1,2,3,4].
Namely, let the axially-symmetric surface S R 3 is generated by revolving about the “vertical” axis O Z a curve  C lying in the X O Z -plane of the Cartesian coordinate system O X Y Z (Figure 3). If the generating curve  C is parameterized as
C = { ( x , z ) , z = z ( x ) , x 0 }
then the principal curvatures of S can be calculated by the well-known formulas
κ π ( x ) = z ˙ ( x ) x 1 + z ˙ 2 ( x ) 1 / 2 , κ μ ( x ) = z ¨ ( x ) 1 + z ˙ 2 ( x ) 3 / 2
where z ˙ ( x ) d z / d x , etc.
Besides, by making use of the explicit formulas in (9) it is straightforward to check directly that the above expressions satisfy the fundamental identity
κ μ ( x ) = d ( x κ π ( x ) ) d x ·
Combined with the defining Equation (4), the last equality produces (after integration) the solution
κ π ( x ) = c + m x 2
where m R is some integration constant. Inserting this back in (10) gives directly
κ μ ( x ) = c m x 2 ·
In this setting it remains only to fix the new integration constant m and for this purpose let us consider the last equality at x R . As it is well-known at such point of S we have (cf. (7))
c ν + 1 2 R and κ μ ( R ) = 1 R μ ( R ) = ν R π ( R ) = ν κ π ( R ) = ν R
and it is easy to conclude from (12) that
m 1 ν 2 R .
Having all this, we can write the final results for both curvatures as
κ π ( x ) = 1 + ν 2 R + ( 1 ν ) R 2 x 2 , κ μ ( x ) = 1 + ν 2 R ( 1 ν ) R 2 x 2 ·
The natural question that arises: for what these formulas are useful, can be answered immediately.
First, we can solve the first equation in (9) with respect to z ˙ ( x ) and by doing this we obtain the differential equation
d z ( x ) d x = ± x κ π 1 ( x κ π ) 2 ·
Second, by inserting there the expression for κ π ( x ) from (15) we end up with the explicit equation
d z ( x ) d x = ± ( ν + 1 ) x 2 ( ν 1 ) R 2 ( R 2 x 2 ) ( ν + 1 ) 2 x 2 ( ν 1 ) 2 R 2 ·
Let us replace x with another independent variable, say χ , via the formula
x = R χ .
Performing this change in (17) one can see that z is defined by the integral
z ( χ ) = ± R ( ν + 1 ) χ 2 ν + 1 d χ 1 χ 2 ( ν + 1 ) 2 χ 2 ( ν 1 ) 2 ·
One thing that this formula make clear is that the parameters R and ν which appear in the expression (7) of the mean curvature are separated multiplicatively and that the problem of finding the CMC surfaces reduces to the evaluation of the integral in (19).
Technically, this depends on the location of the roots of the quartic polynomial under the radical sign. Two of them are known, namely, these are χ 1 , 2 = ± 1 coming from the first factor in the quartic. It is easy to see that the other two roots can be easily found to be χ 3 , 4 = ± σ ( ν ) , where
σ ( ν ) = ν 1 ν + 1 2 .
Still on technical level one needs to know also if these roots are inside or outside of the interval [ 1 , + 1 ] on the real line. Some elementary analysis tells us that for ν ( , 0 ) , σ ( ν ) > 1 while for ν ( 0 , + ) , σ ( ν ) < 1 .
This situation is depicted by the graphics in Figure 4.
We will consider first the case when ν > 0 , which implies that the integral in (19) can be cast in the form
z ( χ ) = ± R χ 2 ν 1 ν + 1 d χ 1 χ 2 χ 2 σ ( ν ) ·
Taking into account that the uniformization of the integral
d χ 1 χ 2 χ 2 σ ( ν )
meaning its reduction to some trivial integrals, is ensured by the Jacobian elliptic function dn( u , k ), where u is the uniformization variable and k 2 = 1 σ ( ν ) is the so-called elliptical modulus, the result of the integration can be written in the form
z = ± R E a m u , k ( ν ) , k ( ν ) ν 1 ν + 1 u , k ( ν ) = 2 ν ν + 1 , ν R + .
Here a m ( u , k ) is the Jacobian amplitude function and R + denotes the set of all positive real numbers. More details about Jacobian elliptic function, elliptic integrals and related subjects can be found in any of the books [9,10] and the well-known Tables of elliptic integrals [11,12,13,14].
We should add also that all of the above mentioned differential geometry books [1,2,3,4] offer an implementation of the computer algebra systems for visualization of the analytical results and in this way provide a great help in increasing the reader’s understanding.
Figure 5 and Figure 6 below illustrate an unduloid and a nodoid generated via Formula (23).
Remark 2. 
Let us notice that in the special case when ν 1 the integral in (21) amounts to
z = ± R 1 χ 2 .
Taken together with (18) they represent a circle x 2 + z 2 = R 2 of radius R playing the role of the profile curve of the sphere with the same radius.
Our next task will be the case when ν < 0 , i.e., the parameter is assumed to be negative. Then, the relevant integral for z is of the form
z ( χ ) = ± r ( ν + 1 ) χ 2 ν + 1 d χ χ 2 1 ( ν 1 ) 2 ( ν + 1 ) 2 χ 2 , ν R .
It will be convenient to replace the negative number ν in the above integral with the positive number ν ˜ via the formula
ν = ν ˜ , ν ˜ R +
Performing the same change in (25) we end up with the problem of evaluation of the integral
z ( χ ) = ± r ν ˜ + 1 ( ν ˜ 1 ) χ 2 d χ χ 2 1 ( ν ˜ + 1 ) 2 ( ν ˜ 1 ) 2 χ 2 = ± r 1 ν ˜ 1 ν ˜ + 1 χ 2 d χ χ 2 1 1 ( ν ˜ 1 ν ˜ + 1 2 χ 2 ) , ν ˜ R + .
In this case the integral that has to be uniformized is
d χ χ 2 1 1 ( ν ˜ 1 ν ˜ + 1 2 χ 2 )
and this can be achieved introducing as uniformization variable u via the function
χ = 1 dn ( u , k ( ν ˜ ) ) , k ( ν ˜ ) = 1 ( ν ˜ 1 ν ˜ + 1 ) 2 = 2 ν ˜ ν ˜ + 1 ·
In these circumstances the evaluation of the integral in (27) produces at the very end the result
z ( u ) = ± r u ν ˜ + 1 ν ˜ 1 E a m ( u , k ( ν ˜ ) ) , k ( ν ˜ ) k 2 ( ν ˜ ) sn ( u , k ( ν ˜ ) ) cn ( u , k ( ν ˜ ) ) dn ( u , k ( ν ˜ ) ) .
An unduloid and a nodoid generated this time via the above formula are shown below in Figure 7 and Figure 8.
Remark 3. 
Let us notice that in the special case when ν ˜ 1 the integral in (28) amounts to
z = ± r d χ χ 2 1
and its uniformization is provided by χ = cosh ( u ) , so that we have also
z = r u u = z r ·
On the other side according to (18) we have as well
x = r χ = r cosh ( z r )
which is obviously the equation of the catenary with a focal distance r.
Actually, the results obtained above can be unified as the following propositions.
Proposition 1. 
With any point within the first quadrant in the ( R , ν ) plane is associated a unique Delaunay surface which profile curve in X O Z plane lies inside the strip ( 0 , R ˚ ) × R . The profile curves of the line R × { 1 } are semi-circles which radii coincide with the respective  R ˚ .
Their explicit parameterizations through the parameters in the specified range are given by the formulas
x = R dn u , k ( ν ) , z = ± R E a m u , k ( ν ) , k ( ν ) ν 1 ν + 1 u , k ( ν ) = 2 ν ν + 1
or
x = R χ , z = ± R 1 χ 2 , χ [ 0 , 1 ] , when ν 1 .
respectively.
Proposition 2. 
With any point within the first quadrant in the ( r , ν ˜ ) plane is associated a unique Delaunay surface which profile curve in X O Z plane lies outside the strip ( 0 , r ˚ ) × R . The profile curves of the line r × { 1 } are catenaries which focal distances coincide with the respective  r ˚ .
Parameterizations of these surfaces associated with the positive numbers ( r , ν ˜ ) are provided by the formulas
x = r dn ( u , k ( ν ˜ ) ) , z = ± r u ν ˜ + 1 ν ˜ 1 E a m ( u , k ( ν ˜ ) ) , k ( ν ˜ ) k 2 ( ν ˜ ) sn ( u , k ( ν ˜ ) ) cn ( u , k ( ν ˜ ) ) dn ( u , k ( ν ˜ ) )
or
x = r cosh z r , z [ , + ] , when ν ˜ 1 .
Definition 1. 
Surfaces described in Proposition  1 will be referred further on as Delaunay surfaces of first kind and those introduced in Proposition  2 as Delaunay surfaces of second kind.
Now we are ready to formulate our results as
Theorem 1. 
The analytical descriptions of all non-trivial Delaunay surfaces of the first, respectively the second kind, are given by Formulas (34) and (35) which are illustrated in Figure 9 and Figure 10.
All results presented up to now can be illustrated schematically as shown in Figure 11.

3. Geometry of Delaunay Surfaces

One of the most fundamental theorem in the classical differential geometry [1,2,3,4] states that any surface S
x ( u , v ) = ( x ( u , v ) , y ( u , v ) , z ( u , v ) ) , ( u , v ) D R 2
in the Euclidean space R 3 is defined up to a motion there (i.e., translation and rotation) by the coefficients of its first and second fundamental forms
I = E d u 2 + 2 F d u d v + G d v 2 , I I = d u 2 + 2 m d u d v + n d v 2 .
These coefficients are given by the formulas
E = E ( u , v ) = x u . x u , F = F ( u , v ) = x u . x v , G = G ( u , v ) = x v . x v = ( u , v ) = x u u . U , m = m ( u , v ) = x u v . U , n = n ( u , v ) = x v v . U
where U is the unit normal vector to S
U = U ( u , v ) = x u × x v | x u × x v | ·
In our settings and after some tedious computations we end up with the following results.
For the first kind Delaunay surfaces these coefficients turns out to be
E = 4 R 2 dn 2 ( u , k ( ν ) ) ( ν + 1 ) 2 , F 0 , G = R 2 dn 2 ( u , k ( ν ) ) = 2 R ν 1 ν + 1 + dn 2 ( u , k ( ν ) ) ν + 1 , m = 0 , n = R 1 2 ν ν + 1 sn 2 ( u , k ( ν ) ) .
Respectively, for the second kind Delaunay surfaces we have
E = 4 r 2 ( ν + 1 ) 2 dn 2 ( u , k ( ν ˜ ) ) , F 0 , G = r 2 dn 2 ( u , k ( ν ˜ ) ) = 2 r 1 ν ˜ ν ˜ + 1 dn 2 ( u , k ( ν ˜ ) ) ( ν ˜ + 1 ) dn 2 ( u , k ( ν ˜ ) ) , m = 0 , n = r ( 1 2 ν ˜ ν ˜ + 1 sn 2 u , k ( ν ˜ ) ) dn 2 ( u , k ( ν ˜ ) ) ·
In the above cited books on differential geometry [1,2,3,4] the reader will find also some classical formulas which describe the two principal types of curvatures at every point of the surface. These are the so-called Gaussian curvature K, and the mean (meaning average) curvature H, the latter being of immediate interest here, i.e.,
K = n m 2 E G F 2 , H = E n + G 2 F m 2 ( E G F 2 ) ·
It can be checked that the second formula above produces for the first kind of Delaunay surfaces the result
H = ν + 1 2 R
while for the second kind of Delaunay surfaces one obtains
H = 1 ν ˜ 2 r ·
However, despite that they seem different, one should keep in mind that according to its very definition (26), the parameter ν ˜ has just the opposite value of ν , so that the above formulas produce obviously absolutely identical results (as should be).
Remark 4. 
Having such complete and explicit description of the Delaunay surfaces opens the possibility to consider more subtle features of these surfaces like geodesics on them. Some preliminary results about geodesics on the unduloids can be found in [15].

4. Conclusions

Up to now we have been interested in purely mathematical questions like how to describe Delaunay surfaces starting with the first principles in differential geometry?
Actually, using the obtained results we have been able to characterize them globally by some intrinsic entities like the coefficients of their fundamental forms or (as we shall see in the Appendix A) via the intrinsic equation of their profile curves.
However, the historical development does not follow exactly the same route. The fundamental fact established by Sturm [6] about minimal property of the lateral area for a fixed volume was quite essential in the treatment of soap bubbles and liquid drops [16,17,18,19] as become evident by Plateau’s [20] experiments. In fact, Plateau is also the originator who coined the names Undulary and Nodary. More recently, these surfaces have been recognized to play the main role in the explanation of the myelin shapes [21], the “pearling” states [22], the capillarity models [23] and that of cells under compression in biology [24,25].
In view of their importance in science and technology we have collected the available parameterizations in the Appendix A below and we will be much obliged to the tempted readers for informing us about any other parameterization which they know but is not listed inside it.
Another our hope is that some of the talented readers of the Mathematics will succeed in finding the relationships between different parameters entering in these parameterizations.

Author Contributions

Conceptualization, C.D.M. and I.M.M.; methodology, C.D.M. and I.M.M.; formal analysis, C.D.M. and I.M.M.; investigation, C.D.M. and I.M.M.; writing—original draft preparation, C.D.M. and I.M.M.; writing—review and editing, C.D.M. and I.M.M.; visualization, C.D.M. and I.M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors would like to thank Vladimir I. Pulov (Varna Technical University) for the useful and constructive discussions on the subject, and to Venelin Chernogorov for his help with some of the figures presented in the paper. The authors are thankful also to all three reviewers for their valuable remarks and suggestions which definitely led to improvement of the presentation.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Here we collect the available parameterizations of the Delaunay surfaces starting with the seminal work of Delaunay on the subject. In all cases the original notation have been preserved.
  • 1841 Ch. Delaunay [5] (p. 313)
    x = a sin φ a tan φ 1 + α sin 2 φ + 0 φ a α d φ cos 2 φ 1 + α sin 2 φ + β z = a cos φ + a 1 + α sin 2 φ , φ R .
    Here a 0 , α > 1 and β —are three arbitrary real constants. The symmetry axis of the generated surfaces is O X and their mean curvatures are 1 2 a ·
  • 1886 A. Lane [26] (p. 133)
    Introducing
    f ( α ) = 0 α a 2 b 2 d α ( a 2 cos 2 α + b 2 sin 2 α ) 3 / 2
    in which a and b are semi-axes of the ellipse rolling on the line, he presents the respective roulette in the form
    x = f ( α ) + c 2 sin α cos α ( a 2 cos 2 α + b 2 sin 2 α ) 1 / 2 α sin α   y = ( a 2 cos 2 α + b 2 sin 2 α ) 1 / 2 α cos α , c 2 = a 2 b 2 .
    Delaunay or Sturm are not mentioned in the paper at all.
  • 1911 G. Loria [27] (p. 127)
    The symmetry axis of the surfaces is O X , a and e are real parameters and s is the natural parameter along the curve. Then
    y 2 = a 2 [ 1 2 e cos s a + e 2 ]
    is the equation of the profile curve in which e is the eccentricity of the rolling curve.
    Loria gives credit about this result to the relatively obscure paper published in 1868 by Spitzer [28], in which he does not mention either Delaunay or Sturm.
    However on the same page Loria presents also the natural equation of Delaunay curves in the form
    R ( s ) = a 1 2 e cos s a + e 2 e ( e cos s a ) ·
    The natural equation is listed here for later comparison.
  • 1969 M. Hady [29] (pp. 380–381)
    The profile curve with symmetry axis O Z is parameterized by the formulas
    x = u , z = ± ( H 0 u 2 + C ) d u u 2 ( H 0 u 2 + C ) 2
    in which H 0 and C are arbitrary real constants.
    Neither Delaunay, Sturm or some of the above mentioned people were cited in the paper.
  • 1980 K. Kenmotsu [30] (p. 150)
    The profile curve of CMC surfaces with symmetry axis O X is presented in the form
    X ( s ; H , B ) = 0 s ( 1 + B sin 2 H t ) ( 1 + B 2 + 2 B sin 2 H t ) 1 / 2 d t , 1 2 | H | 1 + B 2 + 2 B sin 2 H s 1 / 2
    in which s is the arc-length parameter, H denotes the fixed mean curvature and B is another real constant.
    An updated version of this work along nice pictures can be found in the book [31].
  • 2002 I. Mladenov [32] (pp. 22–24)
    The derived parameterization of the Delaunay surfaces is of the form
    x [ u , v ] = R E ( u , k ) ± r F ( u , k ) R k 2 sin u cos u Δ ( u , k ) u ( , ) y [ u , v ] = r Δ ( u , k ) cos v , z [ u , v ] = r Δ ( u , k ) sin v , v [ 0 , 2 π ] .
    in which R , r have the meaning used up to now, k 2 = R 2 r 2 R 2 , Δ ( u , k ) = 1 k 2 sin 2 u and ± signs refer to unduloids/nodoids.
  • 2002 I. Mladenov [33] (p. 329)
    The symmetry axis in this case is O X , K ( k ) and E ( k ) denote the complete elliptic integrals of the first, respectively second kind [11,12,13,14]. The parameterization is provided by the formulas
    x ( u ) = a ( 1 ε ) [ F ( 2 u π 4 , k ) + K ( k ) ] + a ( 1 + ε ) [ E ( 2 u π 4 , k ) + E ( k ) ] y ( u ) = a 1 + ε 2 + 2 ε sin u , k = 2 ε 1 + ε , ε R + .
    The coefficients of the first and the second fundamental forms for this parameterization are found to be
    E = a 2 , F = 0 , G = a 2 ( 1 + ε 2 + 2 ε sin u ) L = a ε ( ε + sin u ) 1 + ε 2 + 2 ε sin u , M = 0 , N = a ( 1 + ε sin u )
    and applying the formula from (42) one immediately obtains that
    H 1 2 a ·
  • 2003 I. Mladenov and J. Oprea [15] (p. 220)
    The parameterization derived in this paper has the form
    x ( u ) = a ( 1 ε ) [ F ( sn ( u , k ) , k ) + K ( k ) ] + a ( 1 + ε ) [ E ( sn ( u , k ) , k ) + E ( k ) ] y ( u ) = a ( 1 + ε ) dn ( u , k ) , ε = 1 b 2 a 2 , k = 2 ε 1 + ε
    while the coefficients of the first and the second fundamental forms were found to be
    E = 4 a 2 dn 2 ( u , k ) , F = 0 , G = a 2 ( 1 + ε ) 2 dn 2 ( u , k ) L = 4 a ε 1 + ε 1 2 1 + ε sn 2 ( u , k ) , M = 0 , N = a ( 1 + ε ) 1 2 ε 1 + ε sn 2 ( u , k ) .
  • 2004 I. Mladenov [34] (pp. 164–166)
    This paper presents the derivation of the conformal parameterizations for unduloids
    x ( u ) = dn ( u , k ) + k cn ( u , k ) , z ( u ) = E ( am ( u , k ) , k ) + k sn ( u , k ) k 2 = 1 m 2 , m 2 < 1
    and respectively, nodoids
    x ( u ) = dn ( u / k , k ) + k cn ( u / k , k ) , k 2 = 1 1 + m 2 , m 2 > 0 z ( u ) = E ( am ( u / k , k ) , k ) + k sn ( u / k , k ) 1 k 2 k u .
    An alternative approach for finding such parameterizations can be seen in [35].
  • 2008 I. Mladenov, M. Hadzhilazova, P. Djondjorov and V. Vassilev [36] (p. 278)
    The profile curve which symmetry axis is taken to be O Z , is parameterized via the natural parameter s in the form r ( s ) , z ( s ) . Relying on the intrinsic equation of the Delaunay curves
    κ ˙ = 2 ( λ κ ) λ κ μ ( 2 λ κ ) 2
    where κ ( s ) is the curvature of the profile curve which explicit form is found to be
    κ ( s ) = λ 1 4 λ μ + 1 4 λ μ sin ( 2 λ s ) 1 2 λ μ + 1 4 λ μ sin ( 2 λ s )
    and this leads to
    r ( s ) = 1 2 λ μ + 1 4 λ μ sin ( 2 λ s ) 2 λ z ( s ) = μ m ( λ , μ ) F ( λ s π 4 , k ) + m ( λ , μ ) λ E ( λ s π 4 , k ) .
    Here λ and μ are real parameters with specific ranges and the other ingredients are given by the expressions
    m ( λ , μ ) = 1 2 λ μ + 1 4 λ μ 2 , k = 2 1 4 λ μ 1 2 λ μ + 1 4 λ μ ·
  • 2014 E. Bendito, M. Bowick and A. Medina [37] (pp. 33–35)
    If the symmetry axis is O Z , roulettes of the ellipses were found to be described by the formulas
    x ( t ) = b a c cos ( t ) a 2 c 2 cos 2 ( t ) , c = a 2 b 2 , a > b z ( t ) = t 0 t a 2 c 2 cos 2 ( u ) d u c sin ( t ) a c cos ( t ) a 2 c 2 cos 2 ( t )
    while that ones of hyperbolas are given by
    x ( t ) = b c cosh ( t ) a c 2 cosh 2 ( t ) a 2 , c = a 2 + b 2 z ( t ) = t 0 t c 2 cosh 2 ( u ) a 2 d u c sinh ( t ) c cosh ( t ) a c 2 cosh 2 ( t ) a 2 ·
  • 2014 B. Athukorallage, T. Paragoda and M. Toda [38] (pp. 6–10)
    The parameterization of undulary with symmetry axis O X is cast into the form
    x ( t ) = b 2 a 0 t d u ( 1 + e cos u ) 1 e 2 cos 2 u y ( t ) = b 1 e cos u 1 + e cos u , e = a 2 b 2 a
    and that one corresponding to undulary is given by the formulas
    x ( t ) = b 2 a 0 t cos u d u ( e + cos u ) e 2 cos 2 u y ( t ) = b e cos u e + cos u ·
  • 2014 V. Pulov, M. Hadzhilazova and I. Mladenov [39] (p. 222)
    In this paper the reader will find the parameterization of the most interesting cases of Delaunay surfaces, namely—unduloids and nodoids, in terms of the complete set of the Weierstrassian elliptic functions, i.e., Weierstrass p-function ( u ) ( u ; g 2 , g 3 ) , σ ( u ) σ ( u ; g 2 , g 3 ) and ζ ( u ) ζ ( u ; g 2 , g 3 ) . Assuming that the O X -axis coincides with the axis of revolution and the generating curves of the surfaces lie in the X O Z -plane, the parameterization of the undularies (upper signs) and nodaries (lower signs) have the form
    x ± ( u ) = 2 c 2 ( 2 a 2 + b 2 ± c 2 ) J 1 ( u + C 1 ± ) c 2 ( 2 a 2 + b 2 ± c 2 ) 2 J 2 ( u + C 1 ± ) ( b 2 ± c 2 ) u + C 2 ± z ± ( u ) = c 6 ( u + C 1 ± ) + 10 a 2 5 b 2 3 c 2 6 ( u + C 1 ± ) 2 a 2 ± b 2 + 3 c 2
    where
    J 1 ( u ) = 1 ( u ) 2 ζ ( u ) u + ln σ ( u u ) σ ( u + u ) J 2 ( u ) = 1 2 ( u ) ( u ) J 1 ( u ) + 2 ( u ) u + ζ ( u u ) + ζ ( u + u ) .
    The above functions are built with the so-called invariants of the Weierstrassian functions given by the expressions
    g 2 g 2 ± = 4 3 ( a 4 a 2 b 2 + b 4 ) , g 3 g 3 ± = 4 27 ( 2 a 6 3 a 4 b 2 3 a 2 b 4 ± 2 b 6 )
    in which C 1 ± and C 2 ± are arbitrary real constants. Here u denotes the argument of ( · ) which produces the value 1 6 ( 2 a 2 b 2 3 c 2 ) for a > | b | , b R , c c ± = a + a 2 b 2 , and ( u ) d ( u ) / d u | u = u , etc. The mean curvature of the generated unduloids, respectively nodoids, is H = ± 1 2 a ·
  • 2017 M. Koiso and B. Palmer [40] (p. 14)
    Relying on variational arguments the authors have concluded that the surface of revolution with symmetry axis O Z will have CMC if its coordinates are given by the formulas
    r = cos θ ± cos 2 θ + H F 2 H , z = 1 2 H θ 0 θ 1 ± cos θ ^ cos 2 θ ^ + H F cos θ ^ d θ ^
    in which H denotes the mean curvature and F is the so-called flux parameter of the surface under consideration.

References

  1. Klimek, G.; Klimek, M. Discovering Curves and Surfaces with MAPLE; Springer: New York, NY, USA, 1997. [Google Scholar]
  2. Gray, A. Modern Differential Geometry of Curves and Surfaces with Mathematica®, 2nd ed.; CRC Press: Boca Raton, FL, USA, 1998. [Google Scholar]
  3. Rovenski, V. Geometry of Curves and Surfaces with Maple; Birkhäuser: New York, NY, USA, 2000. [Google Scholar]
  4. Oprea, J. Differential Geometry and Its Applications, 2nd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2003. [Google Scholar]
  5. Delaunay, C. Sur la surface de révolution dont la courbure moyenne est constante. J. Math. Pures Appl. 1841, 6, 309–314. [Google Scholar]
  6. Sturm, M. Note, à l’occasion de l’article précédent. J. Math. Pures Appl. 1841, 6, 315–320. [Google Scholar]
  7. Eells, J. The surfaces of Delaunay. Math. Intell. 1987, 9, 53–57. [Google Scholar] [CrossRef]
  8. Krivoshapko, S.; Ivanov, V. Encyclopedia of Analytical Surfaces; Springer: New York, NY, USA, 2015. [Google Scholar]
  9. Whittaker, E.; Watson, G. A Course of Modern Analysis, 4th ed.; Cambridge University Press: Cambridge, UK, 1927. [Google Scholar]
  10. Mladenov, I.; Hadzhilazova, M. The Many Faces of Elastica; Springer: Cham, Switzerland, 2017. [Google Scholar]
  11. Byrd, P.; Friedman, M. Handbook of Elliptic Integrals for Engineers and Scientists, 2nd ed.; Springer: New York, NY, USA, 1971. [Google Scholar]
  12. Janhke, E.; Emde, F.; Lösch, F. Tafeln Höherer Funktionen; Teubner Verlag: Stuttgart, Germany, 1960. [Google Scholar]
  13. Abramowitz, M.; Stegun, I.A. (Eds.) Handbook of Mathematical Functions; Dover: New York, NY, USA, 1972. [Google Scholar]
  14. Olver, F.; Lozier, D.; Boisvert, R.; Clark, C. The NIST Handbook of Mathematical Functions; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  15. Mladenov, I.; Oprea, J. Unduloids and their closed geodesics. Geom. Integr. Quant. 2003, 4, 206–234. [Google Scholar]
  16. Pulov, V.; Mladenov, I.M. Rotating liquid drops and Delaunay surfaces. J. Geom. Symmetry Phys. 2019, 54, 55–78. [Google Scholar] [CrossRef]
  17. Thompson, D.W. On Growth and Form; University Press: Cambridge, UK, 1917; Available online: https://www.gutenberg.org/files/55264/55264-h/55264-h.htm#fig63 (accessed on 5 March 2024).
  18. Isenberg, C. The Science of Soap Films and Soap Bubbles; Dover: New York, NY, USA, 1992. [Google Scholar]
  19. Oprea, J. The Mathematics of Soap Films: Explorations with Maple®; AMS: Providence, RI, USA, 2000. [Google Scholar]
  20. Plateau, J. Statique Expérimetal et Théorique des Liquides Soumis Aux Seules Forces Moléculares, 2 Volumes; Gauthier-Villars: Paris, France, 1873. [Google Scholar]
  21. Deuling, H.; Helfrich, W. A theoretical explanation for the myelin shapes of red blood cells. Blood Cells 1977, 73, 713–720. [Google Scholar]
  22. Bar-Ziv, R.; Moses, E. Instability and “pearling” states produced in tubular membranes by competition of curvature and tension. Phys. Rev. Lett. 1994, 73, 1392–1395. [Google Scholar] [CrossRef]
  23. Perera, B.; Paragoda, T.; Dharmasena, D. A survey of Delaunay surfaces with applications in capillary surfaces. J. Geom. Symmetry Phys. 2022, 64, 51–65. [Google Scholar] [CrossRef]
  24. Yoneda, M. Tension at the surface of sea-urchin egg: A critical examination of Cole’s Eexpariment. J. Exp. Biol. 1964, 41, 893–906. [Google Scholar] [CrossRef]
  25. Hadzhilazova, M.; Mladenov, I. Membrane approach to balloons and some related surfaces. Geom. Integr. Quant. 2006, 7, 176–186. [Google Scholar]
  26. Lane, A. Note on a Roulette. Am. J. Math. 1886, VIII, 132–137. [Google Scholar] [CrossRef]
  27. Loria, G. Spezielle Algebraische Transzendent Ebene Kurven: Theorie und Geschichte; Teubner Verlag: Leipzig, Germany, 1911. [Google Scholar]
  28. Spitzer, S. Merkwürdige Eigenschaft derjenigen Curve, welche von Brennpunkte einer Ellipse beschrieben wird, wenn diese auf einer Geraden rollt. Arch. Math. Phys. 1868, VIII, 235–238. [Google Scholar]
  29. Hady, M. Eine chrakterstische Eigenschaft der Sphäre as rotationsfläche. Math. Nachrichten 1969, 42, 379–385. [Google Scholar] [CrossRef]
  30. Kenmotsu, K. Surfaces of revolution with prescribed mean curvature. Tôhoku Math. J. 1980, 32, 147–153. [Google Scholar] [CrossRef]
  31. Kenmotsu, K. Surfaces with Constant Mean Curvature (Translated by Katsuhiro Molriya); AMS: Providence, RI, USA, 2003. [Google Scholar]
  32. Mladenov, I.M. Delaunay surfaces revisited. C. R. Acad. Sci. 2002, 55, 19–24. [Google Scholar]
  33. Mladenov, I.M. New solutions of the shape equation. Eur. Phys. J. B 2002, 29, 327–330. [Google Scholar] [CrossRef]
  34. Mladenov, I.M. Conformal immersions of Delaunay surfaces and their duals. Geom. Integr. Quant. 2004, 5, 158–168. [Google Scholar]
  35. Bracken, P. Delaunay surfaces expressed in terms of a Cartan moving frame. J. Appl. Anal. 2020, 26, 153–160. [Google Scholar] [CrossRef]
  36. Mladenov, I.; Hadzhilazova, M.; Djondjorov, P.; Vassilev, V. On the intrinsic equation behind the Delaunay surfaces. AIP Conf. Proc. 2008, 1079, 274–280. [Google Scholar]
  37. Bendito, E.; Bowick, M.; Medina, A. A natural parameterization of the roulettes of the conics generating the Delaunay surfaces. J. Geom. Symmetry Phys. 2014, 33, 27–45. [Google Scholar]
  38. Athukorallage, B.; Paragoda, T.; Toda, M. Roulettes of conics, Delaunay surfaces and applications. Surv. Math. Math. Sci. 2014, 4, 1–23. [Google Scholar]
  39. Pulov, V.; Hadzhilazova, M.; Mladenov, I.M. Delaunay surfaces in terms of Weierstrassian functions. In Proceedings of the International Conference on Integrability, Recursion Operators and Soliton Interactions, Sofia, Bulgaria, 29–31 August 2012; Aneva, B., Grahovski, G., Ivanov, R., Mladenov, D., Eds.; Avangard Prima: Sofia, Bulgaria, 2014; pp. 218–224. [Google Scholar]
  40. Koiso, M.; Palmer, B. Higher order variations of constant mean curvature surfaces. Calc. Var. 2017, 56, 159. [Google Scholar] [CrossRef]
Figure 1. Delaunay roulettes (the profile curves of the Delaunay surfaces, excluding the plane) are generated as traces of the foci of the non-degenerated conic sections by rolling them along the axis of revolution (the horizontal line).
Figure 1. Delaunay roulettes (the profile curves of the Delaunay surfaces, excluding the plane) are generated as traces of the foci of the non-degenerated conic sections by rolling them along the axis of revolution (the horizontal line).
Mathematics 12 01570 g001
Figure 2. The open parts of the cylinder, sphere, catenoid, unduloid and nodoid shown here are obtained by revolving the profile curves in Figure 1 about symmetry axis taken this time in vertical position.
Figure 2. The open parts of the cylinder, sphere, catenoid, unduloid and nodoid shown here are obtained by revolving the profile curves in Figure 1 about symmetry axis taken this time in vertical position.
Mathematics 12 01570 g002
Figure 3. Geometry of the generating profile curve C .
Figure 3. Geometry of the generating profile curve C .
Mathematics 12 01570 g003
Figure 4. Graphics of σ versus ν . The dashed lines indicate the asymptotic values of the function.
Figure 4. Graphics of σ versus ν . The dashed lines indicate the asymptotic values of the function.
Mathematics 12 01570 g004
Figure 5. Any point within the striped area where ν < 1 on the left produces (via Equation (23)) an unduloid shown on the right.
Figure 5. Any point within the striped area where ν < 1 on the left produces (via Equation (23)) an unduloid shown on the right.
Mathematics 12 01570 g005
Figure 6. With any point in the striped region ν > 1 in the left part is associated a nodoid (via Equation (23)) shown on the right side of the figure.
Figure 6. With any point in the striped region ν > 1 in the left part is associated a nodoid (via Equation (23)) shown on the right side of the figure.
Mathematics 12 01570 g006
Figure 7. With any point in the striped region ν < 1 in the left part is associated an unduloid (via Equation (30)) shown on the right hand side of the figure.
Figure 7. With any point in the striped region ν < 1 in the left part is associated an unduloid (via Equation (30)) shown on the right hand side of the figure.
Mathematics 12 01570 g007
Figure 8. Any point in the striped region ν > 1 on the left hand side generates a nodoid (via Equation (30)) shown on the right hand side of the figure.
Figure 8. Any point in the striped region ν > 1 on the left hand side generates a nodoid (via Equation (30)) shown on the right hand side of the figure.
Mathematics 12 01570 g008
Figure 9. An ensemble of undularies U ν in (left) and a circle C 1 along with an ensemble of nodaries N ν in and with a different scaling the same circle C 1 (in blue) on the right.
Figure 9. An ensemble of undularies U ν in (left) and a circle C 1 along with an ensemble of nodaries N ν in and with a different scaling the same circle C 1 (in blue) on the right.
Mathematics 12 01570 g009
Figure 10. Portions of a few undularies U ν out (left), along a nodary N ν out (right) in the presence of a catenary C ˜ 1 in both cases.
Figure 10. Portions of a few undularies U ν out (left), along a nodary N ν out (right) in the presence of a catenary C ˜ 1 in both cases.
Mathematics 12 01570 g010
Figure 11. Distribution of the Delaunay surfaces specified by the chosen range of ν R . The normal direction to the horizontal line is reserved for the first parameter R from the suggested classification of the Delaunay surfaces.
Figure 11. Distribution of the Delaunay surfaces specified by the chosen range of ν R . The normal direction to the horizontal line is reserved for the first parameter R from the suggested classification of the Delaunay surfaces.
Mathematics 12 01570 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mladenova, C.D.; Mladenov, I.M. Parameterizations of Delaunay Surfaces from Scratch. Mathematics 2024, 12, 1570. https://doi.org/10.3390/math12101570

AMA Style

Mladenova CD, Mladenov IM. Parameterizations of Delaunay Surfaces from Scratch. Mathematics. 2024; 12(10):1570. https://doi.org/10.3390/math12101570

Chicago/Turabian Style

Mladenova, Clementina D., and Ivaïlo M. Mladenov. 2024. "Parameterizations of Delaunay Surfaces from Scratch" Mathematics 12, no. 10: 1570. https://doi.org/10.3390/math12101570

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop