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

Sintesi Tecnica: Il Problema di Newton della Resistenza Minima con Densità Non Costante

Enunciato del Problema
Questo articolo indaga una variante del classico problema di Newton della resistenza minima, che cerca la forma di un corpo solido in movimento a velocità costante attraverso un fluido che offre la minima resistenza. Mentre il modello classico assume un fluido omogeneo, questo studio affronta lo scenario più realistico in cui la densità del fluido diminuisce esponenzialmente con l'altitudine, una condizione caratteristica dell'atmosfera terrestre. La densità è modellata come $\rho(z) = \rho_0 e^{-cz}$, dove $z$ rappresenta l'altitudine. Il corpo si muove lungo la coordinata $z$, e il funzionale della resistenza è derivato sulla base della legge della pressione del seno al quadrato, assumendo assenza di attrito e che le particelle non rimbalzino sulla superficie dopo l'impatto.

Il funzionale della resistenza risultante per una superficie definita da un grafico $u: \Omega \to \mathbb{R}$ è dato da:
$$E[u] = \int_{\Omega} \frac{e^{-u}}{1 + |Du|^2} \, dxdy$$
A differenza del funzionale newtoniano classico, la presenza del peso esponenziale $e^{-u}$ introduce un cambiamento strutturale significativo nell'equazione di Eulero-Lagrange associata. L'articolo afferma esplicitamente che il suo focus è sull'esistenza e la regolarità degli estremali classici $C^2$ (soluzioni radiali) piuttosto che sulla determinazione di un minimizzatore globale, notando che il problema della resistenza minima non vincolato nell'impostazione classica spesso manca di regolarità $C^2$.

Metodologia
L'analisi si concentra sulle soluzioni radiali $u(r)$, dove la superficie è una superficie di rivoluzione attorno all'asse $z$. L'equazione di Eulero-Lagrange per il funzionale $E[u]$ si riduce a un'equazione differenziale ordinaria (ODE) non lineare del secondo ordine:
$$r(1 + 4u'^2 + 3u'^4) - 2u'(1 + u'^2) + 2ru''(3u'^2 - 1) = 0$$
Lo studio affronta il problema ai valori iniziali con le condizioni $u(0) = 0$ e $u'(0) = 0$. Una sfida primaria è la singolarità dell'equazione in $r=0$, che impedisce l'applicazione diretta dei teoremi classici di esistenza per le ODE.

Per superare questo ostacolo, l'autore impiega un argomento del punto fisso basato sul Teorema del Punto Fisso di Banach. L'ODE è riformulata in un'equazione integrale che coinvolge l'inversa di una specifica funzione monotona $f(x) = x/(1+x^2)^2$. Definendo un operatore $T$ su una palla chiusa nello spazio di Banach $C^1([0, R])$, l'autore dimostra che $T$ è una contrazione per un $R$ sufficientemente piccolo, stabilendo così l'esistenza locale di una soluzione $C^2$.

Per analizzare il comportamento globale e il dominio massimale della soluzione, il problema è trasformato in un'analisi del piano delle fasi. La curva è riparametrizzata mediante la lunghezza d'arco $s$, introducendo un angolo $\theta(s)$ tale che $u'(r) = \tan \theta$. Ciò converte l'ODE del secondo ordine in un sistema autonomo di equazioni del primo ordine per le variabili $x(s) = r$ e $\theta(s)$. Il comportamento delle soluzioni è studiato esaminando le traiettorie nel piano delle fasi $\Theta$, in particolare osservando i punti di equilibrio e il flusso vicino alle linee critiche.

Contributi e Risultati Chiave

1. Esistenza e Regolarità Locali: L'articolo dimostra (Teorema 2.1) che il problema ai valori iniziali con $u(0)=u'(0)=0$ ammette una soluzione $u \in C^2([0, R])$ per qualche $R > 0$. Una scoperta cruciale è che la soluzione interseca l'asse di rotazione ortogonalmente ($u'(0)=0$). Ciò contrasta nettamente con il modello classico a densità costante, dove le soluzioni radiali non possono raggiungere l'asse ortogonalmente a meno che la soluzione non sia una funzione costante banale. La prova stabilisce inoltre che $u''(0) = 1/4$.

2. Monotonia e Convessità: L'analisi dimostra che la pendenza $u'(r)$ è strettamente positiva e crescente nel dominio massimale della soluzione. Di conseguenza, la soluzione radiale $u(r)$ è convessa.

3. Dominio Massimale Finito e Pendenza Critica: Il Teorema 2.2 stabilisce che il dominio massimale della soluzione è un intervallo finito $[0, r_M)$. La soluzione termina in un punto in cui il coefficiente di $u''$ nell'ODE si annulla, corrispondente a una pendenza critica:
   $$\lim_{r \to r_M} u'(r) = \frac{1}{\sqrt{3}}$$
   Calcoli numerici indicano $r_M \approx 1.230$ e $u(r_M) \approx 0.228$. L'analisi del piano delle fasi rivela che le traiettorie si avvicinano asintoticamente alla retta $\theta = \pi/6$ (dove $\tan \theta = 1/\sqrt{3}$) e il dominio termina a un raggio finito.

4. Struttura del Piano delle Fasi: Lo studio caratterizza il piano delle fasi del sistema autonomo associato, identificando i punti di equilibrio e mostrando che il flusso attorno a punti specifici consiste in cerchi topologici. È dimostrato che le traiettorie che partono dall'origine terminano sulla retta orizzontale $\theta = \pi/6$ con una coordinata $x$ finita.

Significato e Affermazioni
L'articolo afferma che il decadimento esponenziale della densità del fluido altera fondamentalmente la geometria dei corpi a resistenza minima rispetto al modello omogeneo classico. In particolare, l'esistenza di soluzioni radiali lisce ($C^2$) che intersecano l'asse di simmetria ortogonalmente è un risultato nuovo. Nel modello classico, tale regolarità all'asse è impossibile per soluzioni non banali, che tipicamente richiedono un disco a "fronte piatto" (una discontinuità nella derivata) per minimizzare la resistenza.

L'autore ipotizza che queste soluzioni radiali lisce forniscano una base teorica per forme ottimali in mezzi stratificati. Il vincolo geometrico secondo cui la pendenza è limitata da $1/\sqrt{3}$ al bordo del dominio è evidenziato come una caratteristica significativa per la progettazione di corpi che attraversano l'atmosfera. L'articolo conclude che, sebbene questi estremali $C^2$ esistano e siano regolari, la questione se essi rappresentino minimizzatori globali rispetto ai minimizzatori non lisci classici rimane un problema aperto per la ricerca futura. Il lavoro è presentato come un'indagine matematica rigorosa sull'esistenza e le proprietà di questi specifici estremali, senza proporre nuove applicazioni sperimentali o rivendicare un'immediata implementazione ingegneristica.