Ellipsoïde
En mathématiques, et plus précisément en géométrie euclidienne, un ellipsoïde est une surface du second degré de l'espace euclidien à trois dimensions. Il fait donc partie des quadriques, avec pour caractéristique principale de ne pas posséder de point à l'infini.
![](http://up.wiki.x.io/wikipedia/commons/thumb/b/b4/Ellipsoide.png/330px-Ellipsoide.png)
L'ellipsoïde admet un centre et au moins trois plans de symétrie. L'intersection d'un ellipsoïde avec un plan est une ellipse, un point ou l'ensemble vide.
L'équation d'un ellipsoïde centré à l'origine d'un repère orthonormé et aligné avec les axes du repère est de la forme
où a, b et c, appelés demi-axes de l'ellipsoïde, sont des paramètres strictement positifs.
Les ellipsoïdes dans l'histoire des sciences
modifierDepuis la fin du XVIIe siècle, les propriétés des ellipsoïdes ont fait l'objet d'intenses études par les mathématiciens et les physiciens en raison de leurs applications en physique céleste, en mécanique des fluides et plus récemment en physique nucléaire[1].
La thématique générale est l'étude de la forme d'équilibre des objets déformables en rotation. Selon les forces internes ou externes s'exerçant sur ces objets et leurs éventuels mouvements internes (écoulements, vortex), diverses formes d'équilibre et leur stabilité ont été étudiées par les plus grands mathématiciens. Ces formes à l'équilibre peuvent être des ellipsoïdes de révolution, mais leur stabilité nécessite la connaissance des propriétés des ellipsoïdes triaxiaux[2].
La forme de la Terre (et des planètes)
modifierLa recherche de la forme de la Terre, initiée par Newton est l'archétype de l'étude des corps déformables en rotation uniforme, dont la cohésion est assurée par les forces internes de gravitation, en l'absence de forces externes. Les résultats (proche de l'équilibre) de Newton furent développés par Maclaurin (1742) pour le calcul aux grandes vitesses de rotation. Jacobi (1834) a montré qu'au-delà d'un certain régime critique, les ellipsoïdes triaxiaux peuvent être des figures d'équilibre. Plus tard, dans les séquences de formes apparaissant à des vitesses de rotation de plus en plus élevées, des bifurcations ont été étudiées par Dirichlet, Dedekind. Riemann (1860) a généralisé les études au cas où la rotation n'est pas uniforme et où des mouvements internes (vortex) sont pris en compte. Poincaré et Cartan ont démontré qu'à partir de fréquences de rotation critiques, des bifurcations vers des formes non ellipsoïdales apparaissent (fission, formes asymétriques en poire).
Les gouttes d'eau (et autres objets macroscopiques)
modifierLes gouttes d'eau (ou tout autre liquide) sont des objets dont la cohésion est assurée par des forces de surface (tension superficielle). Mises en rotation, elles subissent des bifurcations entre différentes séquences de formes d'équilibre[3],[4], parmi lesquelles des ellipsoïdes de révolution ou des ellipsoïdes triaxiaux. Ces figures sont difficilement accessibles à l'expérience sur Terre, en raison de l'influence externe de la pesanteur. Il a fallu attendre les années 1990, avec des expériences en apesanteur[5] pour valider les calculs théoriques.
Les noyaux atomiques
modifierDans les années 1930, G.Gamow, puis N.Bohr et JA.Wheeler, pour modéliser le phénomène de fission nucléaire récemment découvert, ont développé le modèle de la goutte liquide du noyau atomique[6]. Dans ce modèle, les formes d'équilibre résultent de la compétition entre l'interaction nucléaire de courte portée et attractive (générant l'analogue d'une tension superficielle), l'interaction coulombienne à longue portée et répulsive et le cas échéant les forces centrifuges liées à la rotation. Plus tard, l'expérience et la théorie on montré que de nombreux noyaux ont des formes ellipsoïdales dans leur état fondamental[7].
Paramétrisations
modifierSelon le contexte et les propriétés que l'on veut étudier, différentes paramétrisations sont utilisées pour décrire un ellipsoïde.
Équation généralisée
modifierDans un repère cartésien en trois dimensions, l'équation d'une surface quadratique est où la matrice A est, par construction, une matrice symétrique réelle. D'après le théorème spectral, elle est diagonalisable et ses valeurs propres sont toutes réelles. Si ces trois valeurs propres sont strictement positives (ou strictement négatives), c'est-à-dire que A est de signature (3, 0) (ou (0, 3)), cette équation définit une quadratique type ellipsoïde. À condition éventuellement de changer tous les coefficients de l'équation par leur opposé, la matrice A est alors définie positive. Le déterminant de A n'étant pas nul, la quadratique possède un centre dont les coordonnées sont et son équation s'écrit sous la forme : avec
Si k est strictement positif, l'ellipsoïde (centré en v et arbitrairement orienté) est alors l'ensemble des points x vérifiant l'équation :
où A1 est réelle, définie positive.
De plus, les vecteurs propres de A1 définissent les axes de l'ellipsoïde et les valeurs propres de A1 sont égales à l'inverse du carré des demi-axes (c'est-à-dire 1/a2, 1/b2 et 1/c2)[8]. Les valeurs singulières de A1, étant égales aux valeurs propres, sont donc égales à l'inverse du carré des demi-axes.
Coordonnées sphériques
modifierUn ellipsoïde peut être paramétré de différentes manières. Une des possibilités, en choisissant l'axe z, est la suivante : où Les paramètres peuvent être vus comme des coordonnées sphériques. Pour un θ constant, nous obtenons une ellipse qui est l'intersection de l'ellipsoïde et d'un plan z = k. Le paramètre ϕ correspond alors à l'anomalie excentrique de cette ellipse. Seuls les ellipsoïdes de révolution possèdent une unique définition de la latitude réduite.
Coordonnées indicées
modifierLes propriétés d'un ellipsoïde sont invariantes par permutation des indices (voir infra), suivant Chandrasekhar, on peut définir un ellipsoïde :
, en ordonnant les demi-axes .
Coordonnées polaires réduites
modifierTout point d'une surface dans un espace à 3 dimensions est défini par 2 paramètres. La paramétrisation dite d'Hill-Wheeler[9] définit les trois demi-axes : où est lié à l'élongation de l'ellipsoïde (excentricité principale) et à son asymétrie (de prolate à oblate ). La symétrie par permutation des axes permet de limiter le domaine de variation de de 0 à 60°.
Espace projectif
modifierEn géométrie projective[réf. nécessaire], l'équation d'un ellipsoïde imaginaire est de la forme
L'équation genre ellipsoïde, cône imaginaire :
Cas particuliers
modifierEllipsoïde triaxial
modifierUn ellipsoïde est triaxial (ou scalène) si ses trois demi-axes sont différents.
Ellipsoïde de révolution
modifierDans le cas où seuls deux demi-axes sont égaux, l'ellipsoïde peut être engendré par la rotation d'une ellipse autour d'un de ses axes. Il s'agit d'un ellipsoïde de révolution, parfois appelé sphéroïde, permettant d'obtenir les miroirs elliptiques des projecteurs de cinéma et les ballons de rugby. On montre aussi que cette surface est optimale pour les dirigeables.
En prenant a = b, l'équation s'écrit :
On obtient un ellipsoïde de révolution d'axe Oz. En effet, les sections par les plans z = k sont des cercles d'axe Oz.
- Si , l’ellipsoïde est dit prolate (c'est-à-dire allongé).
- Si , l’ellipsoïde est une sphère.
- Si , l’ellipsoïde est dit oblate (c'est-à-dire aplati).
La méridienne dans le plan xOz que l'on obtient avec y = 0 est l'ellipse d'équation :
On remarque que l'on passe de l'équation de la méridienne à l'équation de la surface de révolution en remplaçant x2 par x2 + y2.
Ellipsoïde réciproque
modifierL'inversion géométrique d'un ellipsoïde de volume unité, définit l'ellipsoïde réciproque[10]
qui fait correspondre à un ellipsoïde prolate un ellipsoïde oblate (et vice versa) ; les propriétés des deux ellipsoïdes sont intrinsèquement liées (voir infra).
Propriétés de base
modifierVolume
modifierLe volume de l'espace délimité par un ellipsoïde est égal à :
Cette formule donne le volume d'une boule de rayon a dans le cas où les trois demi-axes sont de la même longueur.
Les volumes du plus grand parallélépipède rectangle inscrit et du plus petit parallélépipède rectangle circonscrit sont donnés par les formules suivantes :
Excentricité
modifierSi a ≥ b ≥ c (c'est-à-dire si a est la longueur du plus grand demi-axe et c est la longueur du plus petit demi-axe), on définit l'excentricité principale de l'ellipsoïde par la relation
Aire
modifierL'aire d'un ellipsoïde quelconque est donnée par la formule[11] (voir infra) :
où
et où F(ϕ , k) et E(ϕ , k) sont les intégrales elliptiques incomplètes de première et deuxième espèce respectivement (voir infra).
Propriétés développées
modifierCette section est basée sur le livre de référence de Chandrasekhar[2].
L'utilisation des ellipsoïdes dans la modélisation des systèmes physiques nécessite le calcul de grandeurs, notamment pour les comparer avec celles de la sphère de volume identique (que l'on supposera de volume unité pour simplifier les notations) ; ce sont notamment les moments de volume, les moments de surface et les moments angulaires :
où sont respectivement les éléments de volume, de surface et d'angle solide ; les moments sont normalisés à l'unité pour une sphère avec .
Notations
modifierPour unifier et simplifier les formules on utilisera les coordonnées indicées (voir supra) en ordonnant les demi-axes . Les trois excentricités de l'ellipsoïde sont
Elles ne sont pas indépendantes et seuls interviennent dans les résultats : leurs ratios ainsi qu' un angle , défini par
La plupart des propriétés des ellipsoïdes s'expriment en fonction des intégrales elliptiques de première et de seconde espèce :
Moments de volume
modifierLes moments de volumes sont les plus simples, par changement de variable on obtient : si n est impair, et si n est pair. Les moments d'inertie (n=2) sont particulièrement importants pour l'étude des ellipsoïdes en rotation ; le moment d'inertie d'un ellipsoïde de densité uniforme par rapport à un axe de ses axes de symétrie s'écrit donc :
.
Soit avec les notations non indicées, avec a sur l'axe Ox , avec b sur l'axe Oy et c sur l'axe Oz :
Les moments d'inertie non diagonaux sont tous nuls dans ce système d'axe :
Moments angulaires
modifierCes moments sont des propriétés statiques des ellipsoïdes ; ils diffèrent des moments angulaires (ou moments cinétiques) des solides en rotation. Les moments d'ordre 2 sont liés par l'équation de l'ellipsoïde : ; ils se déduisent de la relation : par le théorème d'Euler sur les fonctions homogènes, on en déduit :
Les moments angulaires d'un ellipsoïde sont fortement liés aux moments de surface de l'ellipsoïde réciproque[10] et interviennent dans le calcul de l'énergie potentielle (voir infra) des ellipsoïdes homogènes.
Moments de surface
modifierIls permettent d'analyser les déviations de la surface de l'ellipsoïde par rapport à celle de la sphère de même volume. Les plus utiles sont liés aux trois composantes de la normale à la surface : .
On en déduit l'aire de l'ellipsoïde (cf. supra les définitions de de k' ):
expression équivalente à celle ci-dessus en fonction des demi-axes a, b et c.
Énergie gravitationnelle et/ou coulombienne
modifierPour les ellipsoïdes dont la cohésion est assurée par la gravité ou pour les ellipsoïdes uniformément chargés, il est nécessaire d'évaluer le potentiel en un point par l'intégrale de volume :
.
Oliver Kellogg (en) a démontré que le potentiel peut s'exprimer en fonction des moments angulaires (cf. supra)[12] :
- .
Carlson a démontré plus précisément que l'énergie potentielle totale d'un ellipsoïde se déduit simplement de l'énergie potentielle d'une sphère équivalente[13] :
,
expression valable même pour les ellipsoïdes inhomogènes mais dont les équipotentielles sont des ellipsoïdes concentriques. La valeur moyenne de dans tout l'ellipsoïde est donc :
Applications et exemples
modifierLes propriétés des ellipsoïdes ont été étudiées depuis des siècles dans la recherche des formes que prennent les systèmes en rotation, Newton s'est intéressé en particulier à l'aplatissement des pôles (modèle ellipsoïdal de la Terre). Les mathématiciens les plus brillants ont contribué à l'étude des figures ellipsoïdales d'équilibre des systèmes en rotation[2] que l'on rencontre dans la nature.
Les corps célestes dont la cohésion est assurée par la gravitation ; en rotation autour d'un axe, les formes d'équilibre minimisent la somme de l'énergie potentielle (cf. ci- dessus) et de l'énergie cinétique de rotation (qui est fonction du moment d'inertie). Aux faibles vitesses de rotation, la figure d'équilibre est un ellipsoïde oblate ; à partir d'une vitesse de rotation critique[2], se produisent des bifurcations vers d'autres figures pour la plupart non ellipsoïdales (voire figure).
La figure ci-contre illustre les formes ellipsoïdales d'équilibre d'un corps céleste en rotation selon les modèles respectifs de Maclaurin et Jacobi : l'abscisse est proportionnelle à la vitesse de rotation angulaire, l'ordonnée donne les longueurs des demi-axes des ellipsoïdes ; l'axe de rotation est selon c.
Aux faibles vitesses de rotation (abscisse ≤ 0,3), les deux modèles prédisent un aplatissement progressif (forme oblate ). Pour les vitesses de rotation plus élevées (abscisse > 0,3), le modèle de Jacobi (lignes continues) diffère de celui Maclaurin (lignes pointillées) ; il prévoit une bifurcation vers une forme triaxiale , qui converge vers une forme prolate . Le modèle de Jacobi est plus réaliste car il prend en compte les mouvements internes (fluide visqueux).
Les deux modèles se limitent aux formes ellipsoïdales ; aux plus hautes vitesses de rotation des bifurcations apparaissent vers d'autres formes d'équilibre, comme par exemple des formes à deux lobes ou plus (formes de « poire » ou de « cacahuète ») ou des tores.
Ont été étudiés aussi les fluides dont la cohésion est assurée par la tension superficielle : par exemple les gouttes d'eau en impesanteur[4]. L'énergie potentielle à prendre en compte est la tension superficielle directement proportionnelle à l'aire de l'ellipsoïde , où T dépend de la nature du liquide (T = 8 × 10−3 J m−2 pour l'eau à température ordinaire). Comme pour les planètes, les gouttes d'eau ont des formes oblates aux faibles vitesses de rotation, qui bifurquent vers d'autres formes aux vitesses plus élevées.
Les noyaux atomiques peuvent être déformés[14] : soit par une brisure spontanée de symétrie dans leur état fondamental, soit lorsque leur énergie coulombienne (répulsive) atteint une valeur critique, soit enfin par rotation induite dans les réactions nucléaires. En première approximation, on peut considérer la matière nucléaire comme incompressible ; dans le modèle de la goutte liquide, l'énergie potentielle dépendante de la forme est donc dominée par la compétition entre une énergie de surface attractive et une énergie coulombienne répulsive. Dans leur état fondamental, la plupart des déformations des noyaux lourds sont de type ellipsoïdal ; pour les noyaux très lourds ou en rotation, ces formes deviennent instables ce qui conduit au phénomène de fission (induite ou spontanée).
Exemples
modifierCertains ellipsoïdes ont des propriétés spécifiques liées à leur domaine d'application :
- Ellipsoïde de Fresnel dans la propagation des ondes mécaniques ou électromagnétiques
- Ellipsoïde de Bessel en topographie
- Ellipsoïde diélectrique pour la propagation de la lumière dans les milieux non-isotropes
- Ellipsoïde de Poinsot en mécanique du solide
- Modèle ellipsoïdal de la Terre (ou ellipsoïde normale) dans l'étude de la forme de la Terre et des planètes ou étoiles.
Notes et références
modifier- Pour un rappel historique détaillé, voir (en) Subrahmanyan Chandrasekhar, « Ellipsoidal figures of equilibrium - Historical account », Comm. Pure Appl. Math., vol. XX, , p. 251-265 (DOI 10.1002/cpa.3160200203, lire en ligne).
- (en) Subrahmanyan Chandrasekhar, Ellipsoidal Figures of Equilibrium, New Haven (USA), Yale University Press, , 253 p. (ISBN 0-486-65258-0).
- (en) S. Chandrasekhar, « The stability of a rotating liquid drop », Proc. R. Soc. A, vol. 286, no 1404, , p. 1-26 (JSTOR 2415184).
- (en) R. A. Brown, L. E. Scriven et Michael James Lighthill, « The shape and stability of rotating liquid drops », Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences, vol. 371, no 1746, , p. 331-357 (DOI 10.1098/rspa.1980.0084).
- (en) T. G. Wang, A. V. Anilkumar, C. P. Lee et K. C. Lin, « Bifurcation of rotating liquid drops: results from USML-1 experiments in Space », Journal of Fluid Mechanics, vol. 276, , p. 389-403 (DOI 10.1017/S0022112094002612).
- Bernard Fernandez, De l'atome au noyau : une approche historique de la physique atomique et de la physique nucléaire, Ellipses, (ISBN 978-2-7298-2784-7, OCLC 69665126).
- « Une histoire de formes de noyaux vieille de 50 ans … enfin résolue ? », sur Ganil.spiral2.eu, (consulté le ).
- (en) Stephen P. Boyd (en), « Lecture 15. Symmetric matrices, quadratic forms, matrix norm, and SVD », sur Université Stanford, automne 2012-2013.
- (en) David Lawrence Hill et John Archibald Wheeler, « Nuclear constitution and the interpretation of fission phenomena », Physical Review, vol. 89, no 5, , p. 1102-1145 (DOI 10.1103/PhysRev.89.1102).
- (en) B. Remaud et G. Royer, « On the energy dependences of ellipsoidal leptodermous systems », Journal of Physics A (en): Mathematical and General, vol. 14, no 11, , p. 2897-2910 (DOI 10.1088/0305-4470/14/11/013).
- (en) « DLMF: § 19.33 », sur dlmf.nist.gov.
- O. D. Kellogg, Foundations of Potential theory, Springer, (ISBN 978-0-486-60144-1 et 978-1-4437-2153-0, OCLC 826373).
- (en) B. C. Carlson, « Ellipsoidal distributions of charge or mass », J. Math. Phys., vol. 2, no 3, , p. 441-450 (DOI 10.1063/1.1703729).
- (en) A. Bohr et B.R.Mottelson, Nuclear structure, vol. 2 - Nuclear deformations, World Scientific, (ISBN 978-981-02-3197-2, OCLC 38112813).