TOPICS
Search

Minimal Surface of Revolution


Calculus of variations can be used to find the curve from a point (x_1,y_1) to a point (x_2,y_2) which, when revolved around the x-axis, yields a surface of smallest surface area A (i.e., the minimal surface). This is equivalent to finding the minimal surface passing through two circular wire frames. The area element is

 dA=2piyds=2piysqrt(1+y^('2))dx,
(1)

so the surface area is

 A=2piintysqrt(1+y^('2))dx,
(2)

and the quantity we are minimizing is

 f=ysqrt(1+y^'^2).
(3)

This equation has f_x=0, so we can use the Beltrami identity

 f-y_x(partialf)/(partialy_x)=a
(4)

to obtain

 ysqrt(1+y^('2))-y^'(yy^')/(sqrt(1+y^('2)))=a
(5)
 y(1+y^('2))-yy^('2)=asqrt(1+y^('2))
(6)
 y=asqrt(1+y^('2))
(7)
 y/(sqrt(1+y^('2)))=a
(8)
 (y^2)/(a^2)-1=y^('2)
(9)
 (dx)/(dy)=1/(y^')=a/(sqrt(y^2-a^2))
(10)
 x=aint(dy)/(sqrt(y^2-a^2))=acosh^(-1)(y/a)+b
(11)
 y=acosh((x-b)/a),
(12)

which is called a catenary, and the surface generated by rotating it is called a catenoid. The two constants a and b are determined from the two implicit equations

y_1=acosh((x_1-b)/a)
(13)
y_2=acosh((x_2-b)/a),
(14)

which cannot be solved analytically.

CatenoidXAxis

The general case is somewhat more complicated than this solution suggests. To see this, consider the minimal surface between two rings of equal radius y_0. Without loss of generality, take the origin at the midpoint of the two rings. Then the two endpoints are located at (-x_0,y_0) and (x_0,y_0), and

 y_0=acosh((-x_0-b)/a)=acosh((x_0-b)/a).
(15)

But cosh(-x)=cosh(x), so

 cosh((-x_0-b)/a)=cosh((-x_0+b)/a).
(16)

Inverting each side

 -x_0-b=-x_0+b,
(17)

so b=0 (as it must by symmetry, since we have chosen the origin between the two rings), and the equation of the minimal surface reduces to

 y=acosh(x/a).
(18)

At the endpoints

 y_0=acosh((x_0)/a),
(19)

but for certain values of x_0 and y_0, this equation has no solutions. The physical interpretation of this fact is that the surface breaks and forms circular disks in each ring to minimize area. Calculus of variations cannot be used to find such discontinuous solutions (known in this case as Goldschmidt solutions). The minimal surfaces for several choices of endpoints are shown above. The first two cases are catenoids, while the third case is a Goldschmidt solution.

To find the maximum value of x_0/y_0 at which catenary solutions can be obtained, let p=1/a. Then (17) gives

 y_0p=cosh(px_0).
(20)

Now, denote the maximum value of x_0 as x_0^*. Then it will be true that dx_0/dp=0. Take d/dp of (20),

 y_0=sinh(px_0)(x_0+p(dx_0)/(dp)).
(21)

Now set dx_0/dp=0

 y_0=x_0sinh(px_0^*).
(22)

From (20),

 py_0^*=cosh(px_0^*).
(23)

Take (23) ÷ (22),

 px_0^*=coth(px_0^*).
(24)

Defining u=px_0^*,

 u=cothu.
(25)

This has solution u=1.1996789403.... From (22), y_0p=coshu. Divide this by (25) to obtain y_0/x_0=sinhu, so the maximum possible value of x_0/y_0 is

 (x_0)/(y_0)=cschu=0.6627434193....
(26)

Therefore, only Goldschmidt ring solutions exist for x_0/y_0>0.6627....

The surface area of the minimal catenoid surface is given by

 A=2(2pi)int_0^(x_0)ysqrt(1+y^('2))dx,
(27)

but since

y=sqrt(1+y^('2))a
(28)
=acosh(x/a),
(29)
A=(4pi)/aint_0^(x_0)y^2dx
(30)
=4piaint_0^(x_0)cosh^2(x/a)dx
(31)
=4piaint_0^(x_0)1/2[cosh((2x)/a)+1]dx
(32)
=2pia[int_0^(x_0)cosh((2x)/a)dx+int_0^(x_0)dx]
(33)
=2pia[a/2sinh((2x)/a)+x]_0^(x_0)
(34)
=pia^2[sinh((2x)/a)+(2x)/a]_0^(x_0)
(35)
=pia^2[sinh((2x_0)/a)+(2x_0)/a].
(36)
CatenoidSolution

Some caution is needed in solving (◇) for a. If we take x_0=1/2 and y_0=1 then (◇) becomes

 1=acosh(1/(2a)),
(37)

which has two solutions: a_1=0.2350... ("deep"), and a_2=0.8483... ("flat"). However, upon plugging these into (◇) with x_0=1/2, we find A_1=6.8456... and A_2=5.9917.... So A_1 is not, in fact, a local minimum, and A_2 is the only true minimal solution.

The surface area of the catenoid solution equals that of the Goldschmidt solution when (◇) equals the area of two disks,

 pia^2[sinh((2x_0)/a)+(2x_0)/a]=2piy_0^2
(38)
 a^2[2sinh((x_0)/a)cosh((x_0)/a)+(2x_0)/a]-2y_0^2=0
(39)
 a^2[cosh((x_0)/a)sqrt(cosh^2((x_0)/a)-1)+(x_0)/a]-y_0^2=0.
(40)

Plugging in

 (y_0)/a=cosh((x_0)/a)
(41)
 (y_0)/asqrt(((y_0)/a)^2-1)+cosh^(-1)((y_0)/a)-((y_0)/a)^2=0.
(42)

Defining

 u=(y_0)/a
(43)

gives

 usqrt(u^2-1)+cosh^(-1)u-u^2=0.
(44)

This has a solution u=1.2113614259. The value of x_0/y_0 for which

 A_(catenary)=A_(2 disks)
(45)

is therefore

 (x_0)/(y_0)=((x_0)/a)/((y_0)/a)=(cosh^(-1)((y_0)/a))/((y_0)/a)=(cosh^(-1)u)/u=0.5276973967.
(46)

For x_0/y_0 in (0.52770,0.6627), the catenary solution has larger area than the two disks, so it exists only as a local minimum.

There also exist solutions with a disk (of radius r) between the rings supported by two catenoids of revolution. The area is larger than that for a simple catenoid, but it is a local minimum. The equation of the positive half of this curve is

 y=c_1cosh(x/(c_1)+c_3).
(47)

At (0,r),

 r=c_1cosh(c_3).
(48)

At (x_0,y_0),

 y_0=c_1cosh((x_0)/(c_1)+c_3).
(49)

The area of the two catenoids is

A_(catenoids)=2(2pi)int_0^(x_0)ysqrt(1+y^('2))dx
(50)
=(4pi)/(c_1)int_0^(x_0)y^2dx
(51)
=4pic_1int_0^(x_0)cosh^2(x/(c_1)+c_3)dx.
(52)

Now let u=x/c_1+c_3, so du=dx/c_1

A=4pic_1^2int_(c_3)^(x_0/x_1+c_3)cosh^2udu
(53)
=4pic_1^21/2int_(c_3)^(x_0/x_1+c_3)[cosh(2u)+1]du
(54)
=2pic_1^2[1/2sinh(2u)+u]_(c_3)^(x_0/x_1+c_3)
(55)
=2pic_1^2{1/2sinh[2((x_0)/(c_1)+c_3)]-1/2sinh(2c_3)+(x_0)/(c_1)}
(56)
=pic_1^2{sinh[2((x_0)/(c_1)+c_3)]-sinh(2c_3)+(2x_0)/(c_1)}.
(57)

The area of the central disk is

 A_(disk)=pir^2=pic_1^2cosh^2c_3,
(58)

so the total area is

 A=pic_1^2{sinh[2((x_0)/(c_1)+c_3)]+[cosh^2c_3-sinh(2c_3)]+(2x_0)/(c_1)}.
(59)

By Plateau's laws, the catenoids meet at an angle of 120 degrees, so

tan30 degrees=[(dy)/(dx)]_(x=0)
(60)
=[sinh(x/(c_1)+c_3)]_(x=0)
(61)
=sinhc_3=1/(sqrt(3))
(62)

and

 c_3=sinh^(-1)(1/(sqrt(3))).
(63)

This means that

cosh^2c_3-sinh(2c_3)=[1+sinh^2c_3]-2sinhc_3sqrt(1+sinh^2c_3)
(64)
=(1+1/3)-2(1/(sqrt(3)))sqrt(1+1/3)
(65)
=4/3-2/(sqrt(3))2/(sqrt(3))=0,
(66)

so

 A=pic_1^2{sinh[2((x_0)/(c_1)+c_3)]+(2x_0)/(c_1)}.
(67)

Now examine x_0/y_0,

 (x_0)/(y_0)=((x_0)/(c_1))/((y_0)/(c_1))=((x_0)/(c_1))/(cosh((x_0)/(c_1)+c_3))=usech(u+c_3),
(68)

where u=x_0/c_1. Finding the maximum ratio of x_0/y_0 gives

 d/(du)((x_0)/(y_0))=sech(u+c_3)-utanh(u+c_3)sech(u+c_3)=0
(69)
 utanh(u+c_3)=1,
(70)

with c_3=sinh^(-1)(1/sqrt(3)) as given above. The solution is u=1.0799632187, so the maximum value of x_0/y_0 for two catenoids with a central disk is y_0=0.4078241702.

If we are interested instead in finding the curve from a point (x_1,y_1) to a point (x_2,y_2) which, when revolved around the y-axis (as opposed to the x-axis), yields a surface of smallest surface area A, we proceed as above. Note that the solution is physically equivalent to that for rotation about the x-axis, but takes on a different mathematical form. The area element is

 dA=2pixds=2pixsqrt(1+y^('2))dx
(71)
 A=2piintxsqrt(1+y^('2))dx,
(72)

and the quantity we are minimizing is

 f=xsqrt(1+y^'^2).
(73)

Taking the derivatives gives

(partialf)/(partialy)=0
(74)
d/(dx)(partialf)/(partialy^')=d/(dx)((xy^')/(sqrt(1+y^('2)))),
(75)

so the Euler-Lagrange differential equation becomes

 (partialf)/(partialy)-d/(dx)(partialf)/(partialy^')=d/(dx)((xy^')/(sqrt(1+y^('2))))=0.
(76)
 (xy^')/(sqrt(1+y^('2)))=a
(77)
 x^2y^('2)=a^2(1+y^('2))
(78)
 y^('2)(x^2-a^2)=a^2
(79)
 (dy)/(dx)=a/(sqrt(x^2-a^2))
(80)
 y=aint(dx)/(sqrt(x^2-a^2))+b=acosh^(-1)(x/a)+b.
(81)

Solving for x then gives

 x=acosh((y-b)/a),
(82)

which is the equation for a catenary. The surface area of the catenoid product by rotation is

A=2piintxsqrt(1+y^('2))dx=2piintxsqrt(1+(a^2)/(x^2-a^2))dx
(83)
=2piintx/(sqrt(x^2-a^2))sqrt((x^2-a^2)+a^2)dx
(84)
=2piint(x^2dx)/(sqrt(x^2-a^2))
(85)
=[x/2sqrt(x^2-a^2)+(a^2)/2ln(x+sqrt(x^2-a^2))]_(x_1)^(x_2)
(86)
=1/2[x_2sqrt(x_2^2-a^2)-x_1sqrt(x_1^2-a^2)+a^2ln((x_2+sqrt(x_2^2-a^2))/(x_1+sqrt(x_1^2-a^2)))].
(87)

Isenberg (1992, p. 80) discusses finding the minimal surface passing through two rings with axes offset from each other.


See also

Minimal Surface, Surface of Revolution

Explore with Wolfram|Alpha

References

Arfken, G. Mathematical Methods for Physicists, 3rd ed. Orlando, FL: Academic Press, pp. 931-937, 1985.Goldstein, H. Classical Mechanics, 2nd ed. Reading, MA: Addison-Wesley, p. 42, 1980.Isenberg, C. The Science of Soap Films and Soap Bubbles. New York: Dover, pp. 79-80 and Appendix III, 1992.

Referenced on Wolfram|Alpha

Minimal Surface of Revolution

Cite this as:

Weisstein, Eric W. "Minimal Surface of Revolution." From MathWorld--A Wolfram Web Resource. https://mathworld.wolfram.com/MinimalSurfaceofRevolution.html

Subject classifications