The Newton's problem assuming non-constant density of the fluid

Résumé technique : Le problème de Newton de la résistance minimale avec une densité non constante

Énoncé du problème
Cet article examine une variante du problème classique de Newton de la résistance minimale, qui cherche la forme d'un corps solide se déplaçant à vitesse constante dans un fluide offrant la résistance la plus faible. Alors que le modèle classique suppose un fluide homogène, cette étude aborde le scénario plus réaliste où la densité du fluide décroît exponentiellement avec l'altitude, une condition caractéristique de l'atmosphère terrestre. La densité est modélisée par $\rho(z) = \rho_0 e^{-cz}$, où $z$ représente l'altitude. Le corps se déplace le long de la coordonnée $z$, et le fonctionnel de résistance est dérivé sur la base de la loi de pression en sinus carré, en supposant l'absence de frottement et que les particules ne se heurtent pas à nouveau à la surface après l'impact.

Le fonctionnel de résistance résultant pour une surface définie par un graphe $u: \Omega \to \mathbb{R}$ est donné par :
$$E[u] = \int_{\Omega} \frac{e^{-u}}{1 + |Du|^2} \, dxdy$$
Contrairement au fonctionnel newtonien classique, la présence du poids exponentiel $e^{-u}$ introduit un changement structurel significatif dans l'équation d'Euler-Lagrange associée. L'article indique explicitement que son objectif porte sur l'existence et la régularité des extrémales classiques $C^2$ (solutions radiales) plutôt que sur la détermination d'un minimiseur global, notant que le problème de résistance minimale non contraint dans le cadre classique manque souvent de régularité $C^2$.

Méthodologie
L'analyse se concentre sur les solutions radiales $u(r)$, où la surface est une surface de révolution autour de l'axe $z$. L'équation d'Euler-Lagrange pour le fonctionnel $E[u]$ se réduit à une équation différentielle ordinaire (EDO) non linéaire du second ordre :
$$r(1 + 4u'^2 + 3u'^4) - 2u'(1 + u'^2) + 2ru''(3u'^2 - 1) = 0$$
L'étude traite du problème aux conditions initiales avec les conditions $u(0) = 0$ et $u'(0) = 0$. Un défi majeur réside dans la singularité de l'équation en $r=0$, qui empêche l'application directe des théorèmes classiques d'existence pour les EDO.

Pour surmonter cela, l'auteur emploie un argument de point fixe basé sur le théorème du point fixe de Banach. L'EDO est reformulée en une équation intégrale impliquant l'inverse d'une fonction monotone spécifique $f(x) = x/(1+x^2)^2$. En définissant un opérateur $T$ sur une boule fermée dans l'espace de Banach $C^1([0, R])$, l'auteur démontre que $T$ est une contraction pour un $R$ suffisamment petit, établissant ainsi l'existence locale d'une solution $C^2$.

Pour analyser le comportement global et le domaine maximal de la solution, le problème est transformé en une analyse dans le plan de phase. La courbe est reparamétrée par la longueur d'arc $s$, introduisant un angle $\theta(s)$ tel que $u'(r) = \tan \theta$. Cela convertit l'EDO du second ordre en un système autonome d'équations du premier ordre pour les variables $x(s) = r$ et $\theta(s)$. Le comportement des solutions est étudié en examinant les trajectoires dans le plan de phase $\Theta$, en particulier en observant les points d'équilibre et l'écoulement près des lignes critiques.

Contributions et résultats clés

1. Existence et régularité locales : L'article démontre (Théorème 2.1) que le problème aux conditions initiales avec $u(0)=u'(0)=0$ admet une solution $u \in C^2([0, R])$ pour un certain $R > 0$. Une découverte cruciale est que la solution intersecte l'axe de rotation orthogonalement ($u'(0)=0$). Cela contraste fortement avec le modèle classique à densité constante, où les solutions radiales ne peuvent atteindre l'axe orthogonalement à moins que la solution ne soit une fonction constante triviale. La preuve établit en outre que $u''(0) = 1/4$.

2. Monotonie et convexité : L'analyse démontre que la pente $u'(r)$ est strictement positive et croissante dans le domaine maximal de la solution. Par conséquent, la solution radiale $u(r)$ est convexe.

3. Domaine maximal fini et pente critique : Le Théorème 2.2 établit que le domaine maximal de la solution est un intervalle fini $[0, r_M)$. La solution se termine à un point où le coefficient de $u''$ dans l'EDO s'annule, correspondant à une pente critique :
   $$\lim_{r \to r_M} u'(r) = \frac{1}{\sqrt{3}}$$
   Des calculs numériques indiquent $r_M \approx 1,230$ et $u(r_M) \approx 0,228$. L'analyse du plan de phase révèle que les trajectoires s'approchent asymptotiquement de la droite $\theta = \pi/6$ (où $\tan \theta = 1/\sqrt{3}$), et que le domaine se termine à un rayon fini.

4. Structure du plan de phase : L'étude caractérise le plan de phase du système autonome associé, identifiant les points d'équilibre et montrant que l'écoulement autour de points spécifiques consiste en des cercles topologiques. Il est démontré que les trajectoires partant de l'origine se terminent sur la ligne horizontale $\theta = \pi/6$ avec une coordonnée $x$ finie.

Importance et affirmations
L'article affirme que la décroissance exponentielle de la densité du fluide modifie fondamentalement la géométrie des corps de résistance minimale par rapport au modèle homogène classique. Plus précisément, l'existence de solutions radiales lisses ($C^2$) qui intersectent l'axe de symétrie orthogonalement est un résultat nouveau. Dans le modèle classique, une telle régularité à l'axe est impossible pour des solutions non triviales, qui nécessitent généralement un disque à « face avant plate » (une discontinuité de la dérivée) pour minimiser la résistance.

L'auteur postule que ces solutions radiales lisses fournissent une base théorique pour les formes optimales dans des milieux stratifiés. La contrainte géométrique selon laquelle la pente est bornée par $1/\sqrt{3}$ au bord du domaine est soulignée comme une caractéristique significative pour la conception de corps traversant l'atmosphère. L'article conclut que, bien que ces extrémales $C^2$ existent et soient régulières, la question de savoir si elles représentent des minimiseurs globaux par rapport aux minimiseurs non lisses classiques demeure un problème ouvert pour des recherches futures. L'ouvrage est présenté comme une investigation mathématique rigoureuse sur l'existence et les propriétés de ces extrémales spécifiques, sans proposer de nouvelles applications expérimentales ni revendiquer une mise en œuvre immédiate en ingénierie.