Théorie des corps élancés

La théorie de Munk (Munk theory) ou théorie des corps élancés (slender bodies theory) est un modèle visant à expliquer le comportement aérodynamique d’un corps élancé en incidence dans des fluides non visqueux. Plutôt que de se concentrer sur la trajectoire des particules du fluide, cette théorie modélise les échanges de quantité de mouvement entre le corps élancé et le fluide traversé.

Si ce bandeau n'est plus pertinent, retirez-le. Cliquez ici pour en savoir plus.

La mise en forme de cet article est à améliorer (février 2021).

La mise en forme du texte ne suit pas les recommandations de Wikipédia : il faut le « wikifier ».

Si ce bandeau n'est plus pertinent, retirez-le. Cliquez ici pour en savoir plus.

Cet article n’est pas rédigé dans un style encyclopédique (février 2021).

Vous pouvez améliorer sa rédaction !

Émise en 1921 par Max M. Munk dans un rapport historique^[1] du National Advisory Committee for Aeronautics^[2], la théorie de Munk conduira à d'importants progrès, notamment dans le calcul de la stabilité aérodynamique des fusées et des missiles en déplacement dans l'atmosphère.

Pour démontrer sa théorie, Munk se base sur l'observation que les dirigeables, alors largement utilisés, doivent déplacer une masse d’air considérable pour se frayer un passage à travers l’atmosphère^[3]. En s'appuyant sur la deuxième loi de Newton, il estime les forces impliquées dans ce déplacement massif d'air, permettant ainsi de modéliser les efforts aérodynamiques en jeu.

Principes et hypothèses

Munk, en s’en tenant à l’observation des échanges de quantité de mouvement entre les dirigeables et l’air environnant, a dégagé une méthode de calcul, simple, des efforts de l’air déplacé en réaction à ce déplacement, appliquée aux dirigeables.

Notion de corps élancé

La première hypothèse de Munk fut de considérer un dirigeable comme un corps élancé, à savoir un corps dont la pente locale est suffisamment faible, c'est-à-dire un corps présentant de faibles variations de diamètre par rapport à sa longueur, (sa surface extérieure se présentant toujours avec un angle assez faible par rapport à la vitesse du fluide loin du corps)^[4]. Une autre définition du corps élancé est qu'un corps est dit élancé si l’évolution de sa section frontale est très progressive le long de son grand axe, de sorte qu’il ne provoque aucune sous ou survitesse de l’écoulement à sa surface^[5]. Puis Munk a divisé l’air environnant ce corps élancé en tranches d’air normales à son axe, cet air étant par hypothèse considéré comme non visqueux.

Ceci posé, il s’est intéressé à ce qui se passe dans l’air présent dans l’une de ces tranches élémentaires lorsque le corps s’avance d’une longueur élémentaire $d x$ correspondant à son épaisseur. Munk considère que les mouvements d’air dus à l’avancée du corps se font uniquement à l’intérieur de la tranche d’air élémentaire considérée, c'est-à-dire sous forme d’un écoulement 2D dans cette tranche (ceci suppose implicitement que le flux d'air généré par chaque tranche n'affecte pas les autres tranches, ceci est obtenu si le corps est suffisamment élancé et que ses sections ne varient pas abruptement^[6])^,^[7]^,^[8].

À l'issue de ces réflexions, Munk a calculé que l’effort normal élémentaire suscité par l’incidence sur la tranche du corps dessinée ci-dessus est^[9] :

dF_{N}=q\,\mathrm {d} S\,\sin(2\alpha )

libellé où $q$ est la pression dynamique de l’écoulement, $d S$ est la variation de section du corps élancé sur l’épaisseur de la tranche, et $α$ l’incidence.

Pour les petites incidences (inférieures à 10°), le sinus peut évidemment être pris comme égal à l'angle en radians, ce qui donne :

dF_{N}=2\;q\;\mathrm {d} S\;\alpha .

Exploitations pratiques de ce résultat

On peut remarquer que lorsque la section du corps reste constante (comme la partie cylindrique médiane d’un dirigeable ou le fuselage de fusée), la force normale de toutes les tranches concernées est nulle. Cette nullité de la portance des parties cylindriques d’un corps est paradoxale pour le néophyte ; par exemple sur l'image ci-contre, il n’est pas naturel de penser que le fuselage cylindrique, bien que placé en incidence par rapport au vent, ne développe pas de portance.

Pourtant ce constat est assez bien vérifié par les mesures en soufflerie, même si les parties cylindriques des fuselages produisent quand-même, aux petits angles d’incidence, une faible portance normale linéaire^[10]^,^[11] (voir plus bas à ce sujet la section déroulante "La vraie valeur du $C_{n\alpha }$ des parties cylindriques" ainsi que les articles Stabilité aérodynamique de la fusée et Rapport des Barrowman).

Une deuxième remarque importante est que lorsque l’accroissement $dS$ de la section élémentaire du corps élancé est positif, la force normale élémentaire est aussi positive (ce qui correspond à une force vers le haut sur les images ci-contre où les corps ont une incidence vers le haut). Dans le cas inverse (accroissement $d S$ négatif) la force normale élémentaire est négative (ce qui est très contre-intuitif).

Il résulte de tout cela que l’intégration des efforts normaux élémentaires au long d’une ogive de fusée, par exemple, résulte en une force normale positive (image ci-dessus). Ce qui implique que lorsqu'une fusée prend de l’incidence, la force normale qui naît sur l’ogive crée un moment d’instabilité par rapport au centre des masses de la fusée (si celui-ci est situé, comme généralement, très en arrière de l'ogive).

De même, la partie avant d'un dirigeable (son ogive) développe pour une incidence positive (vers le haut) une force normale positive (vers le haut).

Par contre, la partie arrière d'un dirigeable, du fait que ses sections normales décroissent, développe une portance négative (une déportance). Cette déportance de la partie arrière du dirigeable, en combinaison avec la portance de la partie avant, va créer un moment déstabilisateur (parfois nommé moment de Munk), moment qui va tendre à augmenter l'incidence du dirigeable (c.-à-d. que si une perturbation atmosphérique place le dirigeable en légère incidence, cette incidence ne va faire que croître jusqu'à placer le dirigeable en travers du vent de sa vitesse). Sur l'image ci-contre, la portance des tranches normales d'un corps à double pointes est représentée par la courbe bleue d'après la théorie des corps élancés. En rouge à marques carrées est la courbe des portances des tranches normales réellement mesurées. Les deux vecteurs bleus représentent la portance de l'avant-corps et la déportance de l'arrière-corps tirées de ces mesures (par intégration). On note que cette dernière déportance mesurée est nettement plus faible que la déportance prédite par la théorie des corps élancés. Appliquée à un corps élancé se terminant par une pointe arrière, la théorie des corps élancés promet en effet une déportance de l'arrière-corps de même module que la portance d'avant-corps. Ces deux forces se composent alors pour créer :

d'une part une résultante nulle (le corps ne développe donc théoriquement aucune portance lorsqu'il est placé en incidence) ;
d'autre part le fameux moment de déstabilisation déjà évoqué (ou moment de Munk). Le centre de portance aérodynamique^[12] d'un tel corps élancé à pointe arrière se trouve projeté à l'infini avant. Nous verrons plus loin que, d'après les mesures en soufflerie, ce centre de portance, ne se trouve pas à l'infini avant mais quand-même très en avant des corps (dans l'image ci-contre, portance de l'avant-corps et déportance de l'arrière-corps sont composés graphiquement pour dessiner le Centre de portance aérodynamique mesuré du corps). La même composition de ces forces résulte également en une portance (ce qui est vérifié par la portance que développent les dirigeables lorsqu'ils prennent de l'incidence).

Les résultats de la théorie des corps élancés furent très tôt utilisés pour des corps autrement véloces que les dirigeables : les fusées (voir l’article Rapport des Barrowman). Dans l’un de ses textes, H. Julian Allen (qui allait étendre la théorie des corps élancés de Munk aux grandes incidences et aux fluides visqueux) précise d’ailleurs à ce sujet, citant Tsien : « les antiques théories potentielles donnant la répartition de la portance sur les dirigeables sont encore applicables à de tels corps même dans un flux supersonique, pourvu que ces corps soient élancés et à faible incidence^[13]. »

Le comportement très contre intuitif du cône de rétreint.

Ainsi qu'il a été dit plus haut, lorsqu'un élément d'un corps montre une évolution négative de ses sections (accroissement $dS$ négatif lorsque l'on va vers l'arrière), comme par exemple un cône de rétreint (c.-à-d. un cône dont la pointe est en arrière), cet élément placé en incidence va développer un $C_{n\alpha }$ négatif. Ce comportement du cône de rétreint est très contre intuitif dans la mesure où l'on perçoit ce cône comme se mettant en travers du vent et donc comme un élément stabilisateur (comme un genre d'empennage). Or (schéma ci-contre) le cône de rétreint placé en incidence génère bien, non pas une portance, mais une déportance (cette déportance se conjuguant, par exemple, avec la portance de l'ogive pour créer un fort moment déstabilisateur, appelé moment de Munk).

L'apport de l'expérience

Imaginons un cône de rétreint articulé autour d'un axe passant par sa pointe. Plaçons-le en incidence comme sur le schéma ci-contre. Cette situation de départ posée, la question est : comment va se comporter le cône ? Il y a deux façons contraires de raisonner :

Raisonnement intuitif :

Le cône s'est mis en travers du vent. En conséquence, il va être entraîné par le vent et augmenter son incidence.

Raisonnement éclairé par la théorie des corps élancés :

Nous savons qu'en application de la théorie des corps élancés, ce cône va développer une déportance (c'est la force bleue symbolisée sur le schéma). Cette force tend à ramener le cône à sa position neutre (sans incidence). Ce raisonnement est très contre intuitif.

Conditions de l'expérience : l'image ci-contre à gauche montre le dispositif expérimental installé dans le flux d'air d'une soufflerie artisanale. Il faut noter que la grande section du cône est recouverte par une calotte sphérique dont le centre est situé sur l'axe de rotation du cône : ainsi l'effet du jeu des pressions sur la calotte sphérique ne participe pas à la rotation.

Expérience du cône de rétreint articulé autour de sa pointe.

Verdict de l'expérience : au vu de la vidéo ci-contre, on remarque bien que chaque mise en incidence du cône est suivie d'un retour au neutre. Le phénomène auquel on est confronté, cependant, c'est que tous ces retours au neutre sont plus forts qu'il serait nécessaire ; ce qui fait que le cône, après un certain nombre d'oscillations, fini par attendre les limites de l'expérience (il finit par percuter de part et d'autre les butées qui ont été préalablement installées, à 25° de chaque côté de l'axe).
On peut donc dire qu'il existe une stabilité statique du cône de rétreint (une tendance instantanée à revenir au neutre), mais que la trop grande force des retours au neutre le met en situation d'instabilité dynamique.

Le cône de rétreint en soufflerie... avec amortissement.

Par chance, il est facile d’annihiler cette trop grande force des retours au neutre : on peut amortir visqueusement la rotation du cône en faisant tourner l'axe de celui-ci dans de la graisse de vaseline : L'énergie trop forte des retours au neutre est alors diminuée : le mouvement est amorti. On observe alors que le cône de rétreint se maintient face au vent de la soufflerie, malgré les turbulences de cette soufflerie (vidéo ci-contre).

Causes de l'instabilité dynamique

Il est facile de comprendre que sur chaque point du cône, la combinaison des vitesses (vitesse de l'écoulement d'air et vitesse géométrique propre du point du cône lors de son retour au neutre) rend les efforts aérodynamiques des retours au neutre plus forts qu'il ne faudrait et en tout cas plus fort que les efforts aérodynamiques qui naîtraient autour du cône immobilisé en rotation (autrement dit, le retour au neutre du cône, qui se fait, rappelons-le contre le vent, augmente la force du vent, donc les efforts de retour au neutre)^[14]^,^[15].

Inversement, les efforts aérodynamiques de retour au neutre d'une ogive, du fait de la même combinaison des vitesses (vitesse de l'écoulement et vitesse propre des points de l'ogive) sont toujours amortis. Il en est de même pour les retours au neutre suscités par les ailerons d'empennage, qui sont eux aussi amortis par le même phénomène de combinaison des vitesses. Cela s'observe à l'amortissement des oscillations d'une fusée lors de son vol, ou même à l'amortissement des oscillations d'une simple girouette.

Confrontation de la théorie avec les mesures expérimentales

Bien que le raisonnement de Munk fasse appel à beaucoup d’hypothèses simplificatrices, force est de constater que les mesures en soufflerie donnent des résultats tout à fait en accord avec la théorie des corps élancés (surtout pour les avant-corps, nous le verrons à l'instant).

Ainsi, l'image ci-contre compare la distribution des forces normales mesurée sur un modèle de dirigeable à 10° d'incidence (en vert clair) avec la distribution déterminée par la Théorie des corps élancés (en rouge). Pour la partie avant du dirigeable, l'accord est assez bon (sauf au point d'arrêt).

De sorte que les calculs de stabilité aérodynamique des fusées sont toujours basés sur ladite théorie (et même, pour les fusées d'amateurs, dans sa version la plus simple). Deux exemples parlants, à ce sujet, sont la force normale développée par une ogive de fusée et le point d’application de cette force normale : les fuséistes amateurs^[16] considèrent la force normale développée par l’ogive comme valant toujours $2\;q\;\alpha \;S_{og}$ (si $q$ est la pression dynamique de l'écoulement, $\alpha$ l'incidence de la fusée en radians et $S_{og}$ la section maximale de l'ogive. Ce qui signifie que le coefficient de force normale $C_{n}$ de l'ogive est $2$ (en référence à la section maximale de l'ogive et par radian), ceci pour des ogives de toutes les formes, y compris l'ogive plate (ou l'absence d'ogive qui donne une tête plate à la fusée) (voir plus bas).

De même, le point d'application de la force normale développée par une ogive est remarquablement bien déterminé par la théorie des corps élancés (la position de cette force normale est donnée dans l'image ci-contre).

D'une façon générale, la position relative de la force normale développée par l'ogive (position nommée XCPA) est donnée par l'équation :

{\frac {XCPA_{og}}{L_{og}}}=1-{\frac {Vol_{og}}{L_{og}S_{og}}}

Vol_{og}

,

S_{og}

et

L_{og}

étant le volume, la section maximale et la longueur de l'ogive (cette équation n'est pas donnée dans le rapport des Barrowman mais dans un rapport suivant de James S. Barrowman seul^[17]^,^[18]).

En ce qui concerne les arrière-corps, les prédictions de la théorie des corps élancés se montrent moins précises (on l'observe sur la première image de cette section où la déportance arrière est, en ordre de grandeur, la moitié de la portance avant alors que la Théorie des corps élancés la prédit égale à la portance avant). La raison en est que ladite Théorie ne prend pas en compte l’épaisseur de la couche Limite sur la partie arrière des corps profilés (couche limite qui devient très importante, en comparaison avec son épaisseur sur la partie avant)^[19]. De sorte que le moment de déstabilisation d'un corps profilé est plus faible que prédit par la théorie.

Une façon d'estimer la déstabilisation d'un corps profilé en incidence est de s'intéresser à son Centre de portance aérodynamique (nommé CPA sur l'image ci-contre à droite et Centre de portance sur l'image ci-contre à gauche)^[20] ^,^[21] : on observe alors que ce CPA est toujours très en avant des corps profilés à faible incidence, parfois même à plus d'une longueur de corps en avant de leur pointe avant (voir à ce sujet l'article Corps de moindre traînée).

L'image de gauche, d'après Hoerner, donne un panorama des Centres de portance (ou CPA) de divers types de corps profilés courants.

Cas particuliers de corps non élancés

L'ogive plate (l'absence d'ogive) ou l'ogive hémisphérique ne sont en rien des corps élancés^[22]. Si ces deux ogives très particulières peuvent être considérées comme développant, en ordre de grandeur, un coefficient de force normale d'environ 2 (par radian et en référence à leur section maximale, soit la section de la partie cylindrique qui les suit) (image ci-contre à gauche)^[23], le point d'application de cette force normale ne peut être donné par la théorie de Munk.
Pour l'ogive hémisphérique, la logique, comme les données expérimentales recueillies en soufflerie (marques bleu clair sur l'image ci-contre à droite), indiquent que le point d'application de la portance de l'ogive hémisphérique se place en son centre géométrique (donc au centre de sa section maximale)^[24]^,^[23].

Quant à l'ogive en tête plate, les mesures en soufflerie placent le point d'application de sa portance notablement en arrière d'elle-même. En tant que disque présenté à l’écoulement, en effet, la tête plate a ceci de particulier qu’elle ne présente aucune surface orientée de façon à transformer une sur- ou sous-pression en force normale. En elle-même, elle ne peut donc développer aucune portance normale. Par contre, toujours pour les faibles incidences, la partie cylindrique du fuselage qui la suit est l’objet d’une forte portance normale^[23]^,^[25] (ce qui donne un $C nα$ de l’ordre de 2,5 par radian, en référence à la section du fuselage cylindrique)^[23]^,^[26]^,^[27]. Plus curieusement encore, le point d’application de cette portance normale peut être en arrière de la tête plate de 0,5 à 2,5 diamètres de fuselage.

La vraie valeur du

C nα

des parties cylindriques

Comme il a été précisé plus haut, selon les enseignements de la Théorie des corps élancés, les parties cylindriques d'un corps placé à faible incidence ne développent aucune portance^[28]. Les fuséistes amateurs se satisfont de cette nullité théorique en s'appuyant sur les hypothèses suivantes :

Cette portance est faible (par rapport à celle de l'ogive, par exemple).
S'il y a une portance du fuselage cylindrique d'une fusée, cette portance est répartie (en général) de part et d'autre du Centre des Masses de la fusée, ce qui fait qu'elle n'a pas d'action sur la stabilité de l'engin^[29].

Cependant, les mesures en soufflerie font état d'une certaine portance de ces parties cylindriques ; c'est ce que l'on voit sur le graphe ci-contre où la pente générale des courbes ne peut être attribuée qu'à l'augmentation de l'élancement L/D des parties cylindriques de corps ogivo-cylindriques.

S. F. Hoerner écrit dans son ouvrage Lift : « Dans une combinaison ogive plus cylindre, l'arrière-corps [c.-à-d. le cylindre] ne devrait théoriquement pas produire de portance, en particulier à de petits angles d'attaque. Il existe cependant des preuves cohérentes selon lesquelles ces corps ont des gradients de portance linéaire plus grands que ceux indiqués par [la théorie de Munk]. Une explication peut être trouvée dans la croissance de la couche limite ( [son épaisseur de] déplacement) […]. »
Et Hoerner ajoute, en renvoyant, entre autres, vers sa fig. 15 qui montre une partie cylindrique longue de 5 diamètres^[30] : « Cet incrément de portance est de l’ordre de 10% de la portance de l’ogive. »

Le même Hoerner, en s’appuyant sur un nombre assez limité de marques expérimentales, donne alors la courbe du gradient de portance linéaire de corps ogivo-cylindriques selon leur élancement total ^[31]. En linéarisant cette courbe entre les élancements totaux 3 et 14 on obtient un $C nα$ (ou gradient de portance linéaire de la partie cylindrique du fuselage) de 0,038 par diamètre de longueur (en référence à la section du fuselage et par radian).

Sur le graphe de Darling ci-contre, la droite rouge a une pente de 0,047 par diamètre : cette pente^[32] est bien dans le même ordre de grandeur que la linéarisation (0,038) de la courbe d'Hoerner.

Évolution de la théorie des corps élancés

S'appuyant sur les travaux d'autres chercheurs, H. Julian Allen (bientôt accompagné de Edward W. Perkins)^[33] a mis en forme une extension aux fluides visqueux (ou réels) de la Théorie des corps élancés de Munk. Dans cette extension, on considère que chaque tranche normale du corps développe, en plus de la force normale dégagée par la Théorie de Munk, une force normale tourbillonnaire naissant du $C_{x}$ traversier appliqué à cette tranche (cette tranche étant donc assimilée à une tranche de cylindre de même élancement que le corps, au premier régime dans les cas les plus simples (image de droite), la vitesse de cette même tranche étant $V\sin(\alpha )$ si $V$ est la vitesse de l'écoulement et $\alpha$ l'incidence du corps. De fait, les observations d'écoulements sur les corps élancés (image de gauche) dévoilent l'existence de deux nappes tourbillonnaires, ces tourbillons étant relâchés périodiquement dans le sillage desdits corps^[34].

Dans le cas simple d'une fusée ogivo-cylindrique en subsonique, par exemple, le $C_{n}$ total peut alors être approché par la formule :

C_{n}=\sin(2\alpha )+1,2\;\eta \;{\frac {S_{p}}{S_{og}}}\sin ^{2}(\alpha )

Dans ce libellé, le premier terme après le signe égal est bien sûr celui déterminé par la théorie des corps élancés ;

\alpha

est l'angle d'incidence (en radians);

\eta

le coefficient donné par la courbe de Goldstein ou celle du NACA RM A50L07, l'élancement du cylindre étant pris comme celui de la fusée^[35]^,^[36] ;

S_{p}

est la surface projetée du corps sur un plan parallèle à son grand axe ;

et

S_{og}

est la surface maximale de l'ogive qui sert ici de référence au

C_{n}

.

Le coefficient

1,2

apparaissant dans ce libellé est le

C_{x}

traversier du cylindre infini en régime sous-critique (ou Premier régime). Il peut être adopté tant que le Nombre de Reynolds diamétral et le Mach transverse, tous deux basés sur la même vitesse transverse

V\sin(\alpha )

, sont inférieurs à respectivement 180 000 et Mach 0,3 ou 0,4 (selon la précision requise)^[37]. C'est cette plage de

C_{x}=1,2

(en fonction du Reynolds et du Mach) que l'on voit sur le graphe 1.

Graphe 1

L'extension d'Allen et Perkins permet de prédire, avec prudence, le comportement, pour toutes les incidences possibles, de corps de révolution ou même de sections non circulaires^[38].

Aux vitesses supersoniques et pour les élancements supérieurs à 6, l'effet de l'élancement du corps peut être négligé, c.-à-d. qu'on peut prendre $\eta =1$ (mais le $C_{x}$ traversier du cylindre infini est très variable avec le Mach et bien supérieur à $1,2$ ^[39].)

Notes et références

↑ Les contributions de Max M. Munk ont été jugées si remarquables par le National Advisory Committee for Aeronautics (NACA, le prédécesseur de la NASA) qu'en 1925 cet organisme a publié un article de Joseph Sweetman Ames intitulé A Résumé of the Advances in Theoretical Aerodynamics Made by Max M. Munk. [1]
↑ (en) Max M. Munk, « THE AERODYNAMIC FORCES ON AIRSHIP HULLS », NACA Report, n^o 184,‎ 1921 (lire en ligne)
↑ Le Zeppelin Hindenburg déplaçait jusqu’à 248 tonnes d’air à la vitesse maximale de 135 km/h.
↑ Munk a d’ailleurs proposé des corrections à ses résultats pour les corps qui seraient moins élancés, corrections que la postérité n’a pas retenues.
↑ « Si un dirigeable très allongé est en mouvement horizontal de translation dans l'air avec une légère incidence, la composante du mouvement de l'air dans la direction de l'axe du navire peut être négligée. L'air livre passage au dirigeable en se déplaçant perpendiculairement à son axe, mais pas en se déplaçant le long du dirigeable. L'air situé dans un plan vertical perpendiculaire au mouvement reste donc dans ce plan, de sorte que le mouvement dans chaque plan perpendiculaire au mouvement peut être considéré comme bidimensionnel. », Munk, THE AERODYNAMIC FORCES ON AIRSHIP HULLS.
↑ Extension of Lighthill’s slender-body theory to moderate aspect ratios, Zhanle Yu, Christophe Eloy, HAL, [2]
↑ On peut admettre cette simplification plus facilement si l’on considère que la pente très faible du corps élancé fait que celui-ci se glisse facilement sous les particules d’air en les écartant (diamétralement) mais sans les entraîner en avant.
↑ (en) H. Julian Allen et Edward W. Perkins, « A study of effects of viscosity on flow over slender inclined bodies of revolution », NACA Report 1048,‎ 1951 (lire en ligne [PDF])
↑ Dans ce résultat n'est pas tenu compte du coefficient multiplicateur $(k 2 - k 1)$ ; $k 1$ est le coefficient de masse additive pour les déplacements longitudinaux du corps et $k 2$ le coefficient de masse additive pour les déplacements normaux du corps dans le plan de l'incidence. De fait, en 1977, Jorgensen écrira : « Pour les élancements L/D suffisamment grands, il est d'usage de négliger l'effet des masses additives [à savoir le coefficient multiplicateur $(k 2 - k 1)$ utilisé par Munk dans sa démonstration]. » (en) Howard Jorgensen, « Prediction of static aerodynamic characteristics for slender bodies alone and with lifting surfaces to very high angles of attack », NASA TR R-474,‎ 1977 (lire en ligne)
↑ Le qualificatif linéaire signifie que cette portance est proportionnelle à l'angle d'incidence.
↑ Cette faible portance normale linéaire du fuselage cylindrique est généralement négligée par les amateurs fuséistes.
↑ Le centre de portance aérodynamique ou (CPA) est le point d'où il est possible de tenir un corps en incidence sans faire usage d'un moment (en usant uniquement d'une force).
↑ ESTIMATION OF THE FORCE AND MOMENTS ACTING ON INCLINED BODIES OF REVOLUTION OF HIGH FINENESS RATIO, H. Julian Allen, RM A 9126, [3]
↑ Cet effet est évoqué p. 32 du texte de Planète-Sciences "Le vol de la fusée, Stabilité et Trajectographie", Juillet 2008 [4], sauf que dans le cas du cône de rétreint, le $C_{n\alpha }$ est négatif, c.-à-d. qu'il va susciter un coefficient d'amortissement négatif, donc un anti-amortissement.
↑ James Barrowman, qui a présenté ces calculs d'amortissement pour la première fois dans The practical calculation of the aerodynamic characteristics of slender finned vehicles, [5] écrit : "Quand l'incidence varie, la combinaison de la vitesse locale [géométriquement due à cette variation d'incidence] [...] et de la vitesse de l'écoulement produit un angle d'attaque local [différent de l'angle d'attaque pour un corps fixe] [...]".
↑ En France, ces fuséistes amateurs n'ont le droit de lancer leurs fusées à feu que sous le contrôle de l'association Planète Sciences.
↑ Barrowman 1967, p. 29
↑ On la trouve aussi dans le Bonnet et Luneau : A. BONNET et J. LUNEAU, AÉRODYNAMIQUE, THÉORIES DE LA DYNAMIQUE DES FLUIDES, CÉPADUÈS-éditions
↑ Hoerner écrit d’ailleurs p 19-15 de son ouvrage Lift : « Il apparaît que [l'épaisseur] de la couche limite le long de l’arrière-corps croît avec l’élancement L/D du corps, réduisant par là même les forces normales négatives qui se développent le long de cet arrière-corps […]. En fait, nous pourrions conjecturer que le moment [déstabilisant] pourrait être seulement la moitié du moment théorique pour un rapport D/L ~ zéro [soit un élancement infini]. »
↑ Ce CPA est appelé par Hoerner et les anglo-saxons, Neutral point, soit Point Neutre : c’est le point d’où l’on peut tenir le corps, à une incidence donnée, sans déployer aucun moment ; c’est donc aussi le point de l’axe où s’applique la résultante des forces aérodynamiques sur le corps à cette incidence.
↑ Si l'on tient le corps (par une articulation sphérique) en un point de son axe situé en aval de son centre de portance aérodynamique, il aura tendance à prendre de l'incidence ; par contre, si on le tient en amont de ce centre de portance aérodynamique, il fera girouette et manifestera une tendance à se replacer à l'incidence nulle.
↑ La pente de leur surface est infinie pour l'ogive plate et trop forte, même à l'écart du point d'arrêt pour l'ogive hémisphérique.
↑ ^{a b c et d} HANDBOOK OF BLUNT-BODY AERODYNAMICS VOL. I, STATIC STABILITY, by J. A. Darling, AD-776-586, Naval Ordinance Laboratory, Silver Spring, MD, DEC 73 [6]
↑ Le point d'application des forces de pressions d'un hémisphère ne peut être qu'en son centre géométrique.
↑ Darling note : "La théorie [des corps élancés] ne s'applique pas très bien aux corps d'élancement faible (L/D < 4) parce que l'abruptesse de leur avant-corps exacerbe leur effet aérodynamique sur la partie [cylindrique] qui les suit."
↑ On peut estimer qu’il y a report du $C nα$ de la tête plate sur la partie cylindrique qui la suit.
↑ Comme on le voit sur cette image il se forme en aval de la tête plate (sans incidence) une zone d’eau morte (zone décollée tourbillonnaire) de ~1,5 diamètres.
↑ Dans la présente réflexion, nous n'évoquerons que les cylindres à base circulaire mais les cylindres à base quelconque sont également justiciables de la théorie des corps élancés...
↑ LE VOL DE LA FUSÉE 2008
↑ Hoerner 1985, p. 19-7
↑ Hoerner 1985, p. 19-9 fig.16
↑ C'est le $C nα$ par diamètre de longueur des parties cylindriques, par radian et en référence section du fuselage.
↑ CHARACTERISTICS OF FLOW OVER INCLINED BODIES OF REVOLUTION, by H. Julian Allen and Edward W. Perkins, NACA RM A50L07 [7]
↑ Ce lâché de tourbillons se fait symétriquement pour les faibles incidences $\alpha$ puis dissymétriquement pour les incidences plus fortes, cette dissymétrie s'approchant de celle existant à l'aval du cylindre dans un écoulement perpendiculaire à son axe, avec création d'un Nombre de Strouhal.
↑ Jorgensen 1977, p. 9
↑ Pour la compréhension de cette courbe de Goldstein, voir l'article Ventilation de l'aval du cylindre
↑ Ces deux conditions sont évidemment respectées plus facilement pour les incidences $\alpha$ modérées.
↑ Jorgensen 1977, p. 11
↑ Drag of circular cylinders for a wide range of reynolds numbers and ach numbers, F. E. Gowen & E. W. Perkins, NACA TN 2960, 1953 https://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19930084018.pdf

Bibliographie

S. F. Hoerner, Résistance à l'avancement dans les fluides, Paris, Gauthier-Villars, 1965 (OCLC 727875556, ASIN B07B4HR4HP).
(en) S. F. Hoerner, Fluid-dynamic drag : theoretical, experimental and statistical information, 1992 (OCLC 228216619).
(en) S. F. Hoerner and H. V. Borst, Fluid-Dynamic Lift, practical information on Aerodynamic and hydrodynamic lift, Liselotte A. Hoerner ed., 1985 ([8]).
(en) James S. Barrowman, The theoretical prediction of the center of pressure, Apogeerockets, 1966 (lire en ligne).
Planète Sciences, Le vol de la fusée : Stabilité et Trajectographie, Version 2.0 - juillet 2008, Planète-Sciences, CNES, 2008 (lire en ligne)
(en) Howard Jorgensen, Prediction of static aerodynamic characteristics for slender bodies alone and with lifting surfaces to very high angles of attack, NASA, 1977 (lire en ligne)

Voir aussi

Articles connexes

[1] Les contributions de Max M. Munk ont été jugées si remarquables par le National Advisory Committee for Aeronautics (NACA, le prédécesseur de la NASA) qu'en 1925 cet organisme a publié un article de Joseph Sweetman Ames intitulé A Résumé of the Advances in Theoretical Aerodynamics Made by Max M. Munk. [1]

[2] (en) Max M. Munk, « THE AERODYNAMIC FORCES ON AIRSHIP HULLS », NACA Report, n^o 184,‎ 1921 (lire en ligne)

[3] Le Zeppelin Hindenburg déplaçait jusqu’à 248 tonnes d’air à la vitesse maximale de 135 km/h.

[4] Munk a d’ailleurs proposé des corrections à ses résultats pour les corps qui seraient moins élancés, corrections que la postérité n’a pas retenues.

[5] « Si un dirigeable très allongé est en mouvement horizontal de translation dans l'air avec une légère incidence, la composante du mouvement de l'air dans la direction de l'axe du navire peut être négligée. L'air livre passage au dirigeable en se déplaçant perpendiculairement à son axe, mais pas en se déplaçant le long du dirigeable. L'air situé dans un plan vertical perpendiculaire au mouvement reste donc dans ce plan, de sorte que le mouvement dans chaque plan perpendiculaire au mouvement peut être considéré comme bidimensionnel. », Munk, THE AERODYNAMIC FORCES ON AIRSHIP HULLS.

[6] Extension of Lighthill’s slender-body theory to moderate aspect ratios, Zhanle Yu, Christophe Eloy, HAL, [2]

[7] On peut admettre cette simplification plus facilement si l’on considère que la pente très faible du corps élancé fait que celui-ci se glisse facilement sous les particules d’air en les écartant (diamétralement) mais sans les entraîner en avant.

[AllenPerkins-8] (en) H. Julian Allen et Edward W. Perkins, « A study of effects of viscosity on flow over slender inclined bodies of revolution », NACA Report 1048,‎ 1951 (lire en ligne [PDF])

[9] Dans ce résultat n'est pas tenu compte du coefficient multiplicateur $(k 2 - k 1)$ ; $k 1$ est le coefficient de masse additive pour les déplacements longitudinaux du corps et $k 2$ le coefficient de masse additive pour les déplacements normaux du corps dans le plan de l'incidence. De fait, en 1977, Jorgensen écrira : « Pour les élancements L/D suffisamment grands, il est d'usage de négliger l'effet des masses additives [à savoir le coefficient multiplicateur $(k 2 - k 1)$ utilisé par Munk dans sa démonstration]. » (en) Howard Jorgensen, « Prediction of static aerodynamic characteristics for slender bodies alone and with lifting surfaces to very high angles of attack », NASA TR R-474,‎ 1977 (lire en ligne)

[10] Le qualificatif linéaire signifie que cette portance est proportionnelle à l'angle d'incidence.

[11] Cette faible portance normale linéaire du fuselage cylindrique est généralement négligée par les amateurs fuséistes.

[12] Le centre de portance aérodynamique ou (CPA) est le point d'où il est possible de tenir un corps en incidence sans faire usage d'un moment (en usant uniquement d'une force).

[EstimAllen-13] ESTIMATION OF THE FORCE AND MOMENTS ACTING ON INCLINED BODIES OF REVOLUTION OF HIGH FINENESS RATIO, H. Julian Allen, RM A 9126, [3]

[14] Cet effet est évoqué p. 32 du texte de Planète-Sciences "Le vol de la fusée, Stabilité et Trajectographie", Juillet 2008 [4], sauf que dans le cas du cône de rétreint, le $C_{n\alpha }$ est négatif, c.-à-d. qu'il va susciter un coefficient d'amortissement négatif, donc un anti-amortissement.

[15] James Barrowman, qui a présenté ces calculs d'amortissement pour la première fois dans The practical calculation of the aerodynamic characteristics of slender finned vehicles, [5] écrit : "Quand l'incidence varie, la combinaison de la vitesse locale [géométriquement due à cette variation d'incidence] [...] et de la vitesse de l'écoulement produit un angle d'attaque local [différent de l'angle d'attaque pour un corps fixe] [...]".

[16] En France, ces fuséistes amateurs n'ont le droit de lancer leurs fusées à feu que sous le contrôle de l'association Planète Sciences.

[p29-17] Barrowman 1967, p. 29

[18] On la trouve aussi dans le Bonnet et Luneau : A. BONNET et J. LUNEAU, AÉRODYNAMIQUE, THÉORIES DE LA DYNAMIQUE DES FLUIDES, CÉPADUÈS-éditions

[19] Hoerner écrit d’ailleurs p 19-15 de son ouvrage Lift : « Il apparaît que [l'épaisseur] de la couche limite le long de l’arrière-corps croît avec l’élancement L/D du corps, réduisant par là même les forces normales négatives qui se développent le long de cet arrière-corps […]. En fait, nous pourrions conjecturer que le moment [déstabilisant] pourrait être seulement la moitié du moment théorique pour un rapport D/L ~ zéro [soit un élancement infini]. »

[20] Ce CPA est appelé par Hoerner et les anglo-saxons, Neutral point, soit Point Neutre : c’est le point d’où l’on peut tenir le corps, à une incidence donnée, sans déployer aucun moment ; c’est donc aussi le point de l’axe où s’applique la résultante des forces aérodynamiques sur le corps à cette incidence.

[21] Si l'on tient le corps (par une articulation sphérique) en un point de son axe situé en aval de son centre de portance aérodynamique, il aura tendance à prendre de l'incidence ; par contre, si on le tient en amont de ce centre de portance aérodynamique, il fera girouette et manifestera une tendance à se replacer à l'incidence nulle.

[22] La pente de leur surface est infinie pour l'ogive plate et trop forte, même à l'écart du point d'arrêt pour l'ogive hémisphérique.

[DarlingHémisph-plate-23] {a b c et d} HANDBOOK OF BLUNT-BODY AERODYNAMICS VOL. I, STATIC STABILITY, by J. A. Darling, AD-776-586, Naval Ordinance Laboratory, Silver Spring, MD, DEC 73 [6]

[24] Le point d'application des forces de pressions d'un hémisphère ne peut être qu'en son centre géométrique.

[25] Darling note : "La théorie [des corps élancés] ne s'applique pas très bien aux corps d'élancement faible (L/D < 4) parce que l'abruptesse de leur avant-corps exacerbe leur effet aérodynamique sur la partie [cylindrique] qui les suit."

[26] On peut estimer qu’il y a report du $C nα$ de la tête plate sur la partie cylindrique qui la suit.

[27] Comme on le voit sur cette image il se forme en aval de la tête plate (sans incidence) une zone d’eau morte (zone décollée tourbillonnaire) de ~1,5 diamètres.

[28] Dans la présente réflexion, nous n'évoquerons que les cylindres à base circulaire mais les cylindres à base quelconque sont également justiciables de la théorie des corps élancés...

[29] LE VOL DE LA FUSÉE 2008

[p19-79-30] Hoerner 1985, p. 19-7

[p19-9fig16-31] Hoerner 1985, p. 19-9 fig.16

[32] C'est le $C nα$ par diamètre de longueur des parties cylindriques, par radian et en référence section du fuselage.

[33] CHARACTERISTICS OF FLOW OVER INCLINED BODIES OF REVOLUTION, by H. Julian Allen and Edward W. Perkins, NACA RM A50L07 [7]

[34] Ce lâché de tourbillons se fait symétriquement pour les faibles incidences $\alpha$ puis dissymétriquement pour les incidences plus fortes, cette dissymétrie s'approchant de celle existant à l'aval du cylindre dans un écoulement perpendiculaire à son axe, avec création d'un Nombre de Strouhal.

[p9Jorgensten-35] Jorgensen 1977, p. 9

[36] Pour la compréhension de cette courbe de Goldstein, voir l'article Ventilation de l'aval du cylindre

[37] Ces deux conditions sont évidemment respectées plus facilement pour les incidences $\alpha$ modérées.

[p11Jorgensten-38] Jorgensen 1977, p. 11

[39] Drag of circular cylinders for a wide range of reynolds numbers and ach numbers, F. E. Gowen & E. W. Perkins, NACA TN 2960, 1953 https://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19930084018.pdf

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]

[16]

[17]

[18]

[19]

[20]

[21]

[22]

[23]

[24]

[25]

[26]

[27]

[28]

[29]

[30]

[31]

[32]

[33]

[34]

[35]

[36]

[37]

[38]

[39]