User:Peter Mercator/Draft for ellipsoid

From Wikipedia, the free encyclopedia
Jump to: navigation, search


This page is about the geometric three dimensional shape. For other uses, see Peter Mercator/Draft for ellipsoid (disambiguation).
Tri-axial ellipsoid with distinct semi-axes a,b and c
Ellipsoids of revolution (spheroid) with a pair of equal semi-axes (a) and a distinct third semi-axis (c) which is an axis of symmetry. The ellipsoid is oblate or prolate as c is less than or greater thana.

An ellipsoid is a closed quadric surface that is a three dimensional analogue of an ellipse. The standard equation of an ellipsoid centered at the origin of a Cartesian coordinate system is

{x^2 \over a^2}+{y^2 \over b^2}+{z^2 \over c^2}=1,

The points (a,0,0), (0,b,0) and (0,0,c) lie on the surface and the line segments from the origin to these points are called the semi-principal axes of length a, b, c. They correspond to the semi-major axis and semi-minor axis of the appropriate ellipses.

There are four distinct cases of which one is degenerate:

  • a>b>c\quad: \,\!   Tri-axial or (rarely) scalene ellipsoid;
  • a=b>c\quad: \,\!   oblate ellipsoid of revolution ( oblate spheroid );
  • a=b<c\quad:\,\!   prolate ellipsoid of revolution ( prolate spheroid );
  • a=b=c\quad: \,\!   the degenerate case of a Sphere;

Mathematical literature often uses 'ellipsoid' in place of 'tri-axial ellipsoid'. Scientific literature (particularly geodesy) often uses 'ellipsoid' in place of 'ellipsoid of revolution' and only applies the adjective 'tri-axial' when treating the general case. Older literature uses 'spheroid' in place of 'ellipsoid of revolution'.

Any planar cross section passing through the center of an ellipsoid forms an ellipse on its surface: this degenerates to a circle for sections normal to the symmetry axis of an ellipsoid of revolution (or all sections when the ellipsoid degenerates to a sphere.)

Generalised equations[edit]

More generally, an arbitrarily oriented ellipsoid, centered at v, is defined by the equation

(\mathbf{x-v})^\mathrm{T} A^{-1} (\mathbf{x-v}) = 1,

where A is a positive definite matrix and x, v are vectors.

The eigenvectors of A define the principal directions of the ellipsoid and theeigenvalues of A are the inverse squares of the semi-axes: {1/a^2}, {1/b^2}  and {1/c^2}. An invertible linear transformation applied to a sphere produces an ellipsoid, which can be brought into the above standard form by a suitable rotation, a consequence of thepolar decomposition (also, see spectral theorem). If the linear transformation is represented by asymmetric 3-by-3 matrix, then the eigenvectors of the matrix are orthogonal (due to the spectral theorem) and represent the directions of the axes of the ellipsoid: the lengths of the semiaxes are given by the eigenvalues. The singular value decomposition and polar decomposition are matrix decompositions closely related to these geometric observations.

Parameterization[edit]

The surface of the ellipsoid may be parameterized in several ways. One possible choice which singles out the 'z'-axis is:

\begin{align}
x&=a\,\cos u\cos v,\\
y&=b\,\cos u\sin v,\\
z&=c\,\sin u;\end{align}\,\!
where

-{\pi}/{2}\leq u\leq+{\pi}/{2},
\qquad
-\pi\leq v\leq+\pi.\!\,\!

There is no simple interpretation of this parameterization. For constant u, that is on the ellipse which is the intercept with a constant z plane, v then plays the role of theeccentric anomaly for that ellipse. For constant v on a plane through the Oz axis the parameter u plays the same role for the ellipse of intersection. Two other similar parameterizations are possible, each with their own interpretations. Only on an ellipse of revolution can a unique definition of reduced latitude be made.

Volume and surface area[edit]

Volume[edit]

The volume of an ellipsoid is given by the formula

 V= \frac{4}{3}\pi abc =   \frac{4}{3}\pi\sqrt{\det(A)} .\,\!

Note that this equation reduces to that of the volume of a sphere when all three elliptic radii are equal, and to that of an oblate or prolate spheroidwhen two of them are equal. The volumes of the maximum inscribed and minimum circumscribed boxes are respectively:

V_{\rm max}=\frac{8}{3\sqrt 3} abc, \qquad V_{\rm min}= 8abc.

The volume of an ellipse of dimension higher than 3 can be calculated using the product of square roots of the eignvalues of A^{-1} and the dimensional constant given for thevolume of a hypersphere. One can also define ellipsoids in higher dimensions, as the images of spheres under invertible linear transformations. The spectral theorem can again be used to obtain a standard equation akin to the one given above.

Surface area[edit]

The surface area of a general (tri-axial) ellipsoid is[1][2]


S=2\pi c^2+ \frac{2\pi ab}{\sin\phi}
\left(E(\phi,k)\, \sin^2\phi  + F(\phi,k)\, \cos^2\phi \right),
where
 
\cos\phi = \frac{c}{a}, \qquad
k^2 =\frac{a^2(b^2-c^2)}{b^2(a^2-c^2)}, \qquad
a\ge b \ge c,

and F(φ,k), E(φ,k) are incomplete elliptic integrals of the first and second kind respectively.[1]

The surface area of an ellipsoid of revolution (or spheroid) may be expressed in terms ofelementary functions:

S_{\rm oblate} =  2\pi a^2\left(1+\frac{1-e^2}{e}\tanh^{-1}e\right)
\quad\mbox{where}\quad e^2=1-\frac{c^2}{a^2}\quad(c<a).
S_{\rm prolate} =  2\pi a^2\left(1+\frac{c}{ae}\sin^{-1}e\right)
\quad\qquad\mbox{where}\;\quad e^2=1-\frac{a^2}{c^2}\quad(c>a).

In both cases e may again be identified as the eccentricity of the ellipse formed by the cross section through the symmetry axis. (See ellipse). Derivations of these results may be found at Mathworld

Approximate formula[edit]

S\approx 4\pi\!\left(\frac{ a^p b^p+a^p c^p+b^p c^p }{3}\right)^{1/p}.\,\!

Here p ≈ 1.6075 yields a relative error of at most 1.061% [3]; a value of p = 8/5 = 1.6 is optimal for nearly spherical ellipsoids, with a relative error of at most 1.178%.

In the "flat" limit of c much smaller than a, b, the area is approximately 2πab.

Dynamical properties[edit]

The mass of an ellipsoid of uniform density ρ is:

m = \rho V = \rho \frac{4}{3} \pi abc\,\!

The moments of inertia of an ellipsoid of uniform density are:

I_{\mathrm{xx}} =  \frac{1}{5} m( b^2+c^2),\qquad
I_{\mathrm{yy}} =  \frac{1}{5} m(c^2+a^2),\qquad
I_{\mathrm{zz}} =  \frac{1}{5} m(a^2+b^2),
 I_{\mathrm{xy}}=  I_{\mathrm{yz}} = I_{\mathrm{zx}} =0.\,\!

For a=b=c these moments of inertia reduce to those for a sphere of uniform-density.

Artist's conception of Haumea, a tri-axialdwarf planet, with its two moons

Tri-axial ellipsoids and cuboids rotate stably along their major or minor axes, but not along their median axis. This can be seen experimentally by throwing an eraser with some spin. In addition, moment of inertia considerations mean that rotation along the major axis is more easily perturbed than rotation along the minor axis.[4]

One practical effect of this is that triaxial astronomical bodies such as Haumea generally rotate along their minor axes (as does the Earth, which is merely oblate); in addition, because oftidal locking, tri-axial moons in synchronous orbit such asMimas orbit with their major axis aligned radially to their planet.

A relaxed ellipsoid, that is, one in hydrostatic equilibrium, has an oblatenessa−c directly proportional to its mean density and mean radius. Ellipsoids with a differentiated interior—that is, a denser core than mantle—have a lower oblateness than a homogeneous body. Over all, the ratio (b–c)/(a−c) is approximately 0.25, though this drops for rapidly rotating bodies.[5]

Fluid properties[edit]

The ellipsoid is the most general shape for which it has been possible to calculate the creeping flow of fluid around the solid shape. The calculations include the force required to translate through a fluid and to rotate within it. Applications include determining the size and shape of large molecules, the sinking rate of small particles, and the swimming abilities ofmicroorganisms.[6]


See also[edit]

References[edit]

  1. ^ F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, editors, 2010, NIST Handbook of Mathematical Functions (Cambridge University Press), available on line at http://dlmf.nist.gov/19.33 (see next reference).
  2. ^ NIST (National Institute of Standards and Technology) athttp://www.nist.gov
  3. ^ [http://www.numericana.com/answer/ellipsoid.htm#thomsen Final answers] by Gerard P. Michon (2004-05-13). See Thomsen's formulas and Cantrell's comments. CAUTION: Note that the labelling of the semi-axes differs to the notation used in this article.
  4. ^ Goldstein, H G (1980). Classical Mechanics, (2nd edition) Chapter 5.
  5. ^ "Shapes of the Saturnian Icy Satellites". 
  6. ^ Dusenbery, David B. (2009).Living at Micro Scale, Harvard University Press, Cambridge, Mass. ISBN 978-0-674-03116-6.

External links[edit]