Sphere–cylinder intersection

From HandWiki

In the theory of analytic geometry for real three-dimensional space, the curve formed from the intersection between a sphere and a cylinder can be a circle, a point, the empty set, or a special type of curve.

For the analysis of this situation, assume (without loss of generality) that the axis of the cylinder coincides with the z-axis; points on the cylinder (with radius [math]\displaystyle{ r }[/math]) satisfy

[math]\displaystyle{ x^2 + y^2 = r^2. }[/math]

We also assume that the sphere, with radius [math]\displaystyle{ R }[/math] is centered at a point on the positive x-axis, at point [math]\displaystyle{ (a, 0, 0) }[/math]. Its points satisfy

[math]\displaystyle{ (x-a)^2 + y^2 + z^2 = R^2. }[/math]

The intersection is the collection of points satisfying both equations.

Trivial cases

Sphere lies entirely within cylinder

If [math]\displaystyle{ a+R \lt r }[/math], the sphere lies entirely in the interior of the cylinder. The intersection is the empty set.

Sphere touches cylinder in one point

If the sphere is smaller than the cylinder ([math]\displaystyle{ R \lt r }[/math]) and [math]\displaystyle{ a+R = r }[/math], the sphere lies in the interior of the cylinder except for one point. The intersection is the single point [math]\displaystyle{ (r, 0, 0) }[/math].

Sphere centered on cylinder axis

If the center of the sphere lies on the axis of the cylinder, [math]\displaystyle{ a = 0 }[/math]. In that case, the intersection consists of two circles of radius [math]\displaystyle{ r }[/math]. These circles lie in the planes

[math]\displaystyle{ z = \pm\sqrt{R^2 - r^2}; }[/math]

If [math]\displaystyle{ r = R }[/math], the intersection is a single circle in the plane [math]\displaystyle{ z = 0 }[/math].

Non-trivial cases

Subtracting the two equations given above gives

[math]\displaystyle{ z^2 + (r^2 - R^2 + a^2) = 2ax. }[/math]

Since [math]\displaystyle{ x }[/math] is a quadratic function of [math]\displaystyle{ z }[/math], the projection of the intersection onto the xz-plane is the section of an orthogonal parabola; it is only a section due to the fact that [math]\displaystyle{ -r \lt x \lt r }[/math]. The vertex of the parabola lies at point [math]\displaystyle{ (-b, 0, 0) }[/math], where

[math]\displaystyle{ b = \frac{R^2 - r^2 - a^2}{2a}. }[/math]

Intersection consists of two closed curves

If [math]\displaystyle{ R \gt r + a }[/math], the condition [math]\displaystyle{ x \lt r }[/math] cuts the parabola into two segments. In this case, the intersection of sphere and cylinder consists of two closed curves, which are mirror images of each other. Their projection in the xy-plane are circles of radius [math]\displaystyle{ r }[/math].

Each part of the intersection can be parametrized by an angle [math]\displaystyle{ \phi }[/math]:

[math]\displaystyle{ (x,y,z) = \left(r\cos\phi,r\sin\phi,\pm\sqrt{2a(b + r\cos\phi)}\right). }[/math]

The curves contain the following extreme points:

[math]\displaystyle{ \left(-r, 0, \pm\sqrt{R^2 - (r+a)^2}\right);\quad \left(0, \pm r, \pm\sqrt{R^2 - (r-a)(r+a)}\right);\quad \left(+r, 0, \pm\sqrt{R^2 - (r-a)^2}\right). }[/math]

Intersection is a single closed curve

If [math]\displaystyle{ R \lt r + a }[/math], the intersection of sphere and cylinder consists of a single closed curve. It can be described by the same parameter equation as in the previous section, but the angle [math]\displaystyle{ \phi }[/math] must be restricted to [math]\displaystyle{ -\phi_0 \lt \phi \lt +\phi_0 }[/math], where [math]\displaystyle{ \cos\phi_0 = -b/r }[/math].

The curve contains the following extreme points:

[math]\displaystyle{ \left(-b, \pm\sqrt{r^2-b^2}, 0\right);\quad \left(0, \pm r, \pm\sqrt{R^2 - (r-a)(r+a)}\right);\quad \left(+r, 0, \pm\sqrt{R^2 - (r-a)^2}\right). }[/math]

Limiting case

Viviani's curve as intersection of a sphere and a cylinder

In the case [math]\displaystyle{ R = r + a }[/math], the cylinder and sphere are tangential to each other at point [math]\displaystyle{ (r, 0, 0) }[/math]. The intersection resembles a figure eight: it is a closed curve which intersects itself. The above parametrization becomes

[math]\displaystyle{ (x,y,z) = \left(r\cos\phi,r\sin\phi,2\sqrt{ar}\cos\frac{\phi}{2}\right), }[/math]

where [math]\displaystyle{ \phi }[/math] now goes through two full revolutions.

In the special case [math]\displaystyle{ a = r, R = 2r }[/math], the intersection is known as Viviani's curve. Its parameter representation is

[math]\displaystyle{ (x,y,z) = \left(r\cos\phi,r\sin\phi,R\cos\frac\phi2\right). }[/math]

The volume of the intersection of the two bodies, sometimes called Viviani's volume, is[1] [2] [3]

[math]\displaystyle{ V=2\int\int_{(x-a)^2+y^2\le r^2,x^2+y^2\le R^2} \sqrt{R^2-x^2-y^2}dxdy =\left(\frac{2\pi}{3}-\frac89\right)R^3. }[/math]

See also

References

  1. Lamarche, F.; Leroy, Claude (1990). "Evaluation of the volume of intersection of a sphere with a cylinder by elliptic integrals". Comput. Phys. Commun. 59 (2): 359–369. doi:10.1016/0010-4655(90)90184-3. Bibcode1990CoPhC..59..359L. 
  2. Viviani, V. (1692), """", Acta Eruditorum: 273–279 
  3. Woodhouse, Robert (1801). "VII. Demonstration of a theorem, by which such portions of the solidity of a sphere are assigned as admit an algebraic expression". Philos. Trans. R. Soc. Lond. 91: 153. doi:10.1098/rstl.1801.0009.