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

Technische Samenvatting: Het Newton-probleem van minimale weerstand met niet-constante dichtheid

Probleemstelling
Dit artikel onderzoekt een variant van Newton's klassieke probleem van minimale weerstand, dat de vorm zoekt van een vast lichaam dat met constante snelheid door een vloeistof beweegt en de minste weerstand ondervindt. Waar het klassieke model een homogene vloeistof veronderstelt, behandelt deze studie het realistischere scenario waarbij de vloeistofdichtheid exponentieel afneemt met de hoogte, een voorwaarde die kenmerkend is voor de aardatmosfeer. De dichtheid wordt gemodelleerd als $\rho(z) = \rho_0 e^{-cz}$, waarbij $z$ de hoogte voorstelt. Het lichaam beweegt langs de $z$-coördinaat en het weerstandsfunctionaal wordt afgeleid op basis van de wet voor druk naar het kwadraat van de sinus, onder de aanname van geen wrijving en dat deeltjes na impact niet opnieuw met het oppervlak botsen.

Het resulterende weerstandsfunctionaal voor een oppervlak gedefinieerd door een grafiek $u: \Omega \to \mathbb{R}$ wordt gegeven door:
$$E[u] = \int_{\Omega} \frac{e^{-u}}{1 + |Du|^2} \, dxdy$$
In tegenstelling tot het klassieke Newton-functionaal introduceert de aanwezigheid van de exponentiële weging $e^{-u}$ een significante structurele verandering in de bijbehorende Euler-Lagrange-vergelijking. Het artikel stelt expliciet dat de focus ligt op het bestaan en de regulariteit van klassieke $C^2$-extrema (radiale oplossingen) en niet op de bepaling van een globaal minimizer, met de opmerking dat het onbeperkte probleem van minimale weerstand in de klassieke setting vaak ontbreekt aan $C^2$-regulariteit.

Methodologie
De analyse richt zich op radiale oplossingen $u(r)$, waarbij het oppervlak een omwenteling is rond de $z$-as. De Euler-Lagrange-vergelijking voor het functionaal $E[u]$ reduceert tot een niet-lineaire gewone differentiaalvergelijking (ODE) van de tweede orde:
$$r(1 + 4u'^2 + 3u'^4) - 2u'(1 + u'^2) + 2ru''(3u'^2 - 1) = 0$$
De studie behandelt het beginwaardeprobleem met voorwaarden $u(0) = 0$ en $u'(0) = 0$. Een primaire uitdaging is de singulariteit van de vergelijking bij $r=0$, wat de directe toepassing van klassieke ODE-bestaansstellingen verhindert.

Om dit te overwinnen, past de auteur een vastpuntargument toe gebaseerd op de Banach-vastpuntstelling. De ODE wordt herschreven tot een integraalvergelijking die de inverse bevat van een specifieke monotoone functie $f(x) = x/(1+x^2)^2$. Door een operator $T$ te definiëren op een gesloten bal in de Banachruimte $C^1([0, R])$, bewijst de auteur dat $T$ een contractie is voor een voldoende kleine $R$, waardoor het lokale bestaan van een $C^2$-oplossing wordt vastgesteld.

Om het globale gedrag en het maximale domein van de oplossing te analyseren, wordt het probleem getransformeerd tot een fasevlakanalyse. De kromme wordt gereparametriseerd door booglengte $s$, waarbij een hoek $\theta(s)$ wordt geïntroduceerd zodat $u'(r) = \tan \theta$. Dit zet de ODE van de tweede orde om in een autonoom stelsel van eerste-orde vergelijkingen voor de variabelen $x(s) = r$ en $\theta(s)$. Het gedrag van de oplossingen wordt bestudeerd door de trajecten in het fasevlak $\Theta$ te onderzoeken, met name met het oog op evenwichtspunten en de stroming in de buurt van kritieke lijnen.

Belangrijkste Bijdragen en Resultaten

1. Lokaal Bestaan en Regulariteit: Het artikel bewijst (Stelling 2.1) dat het beginwaardeprobleem met $u(0)=u'(0)=0$ een oplossing $u \in C^2([0, R])$ toestaat voor een zekere $R > 0$. Een cruciale bevinding is dat de oplossing het rotatieas orthogonaal snijdt ($u'(0)=0$). Dit staat in schril contrast met het klassieke model met constante dichtheid, waarbij radiale oplossingen de as niet orthogonaal kunnen bereiken tenzij de oplossing een triviale constante functie is. Het bewijs stelt verder vast dat $u''(0) = 1/4$.

2. Monotonie en Convexiteit: De analyse toont aan dat de helling $u'(r)$ strikt positief en stijgend is in het maximale domein van de oplossing. Bijgevolg is de radiale oplossing $u(r)$ convex.

3. Eindig Maximaal Domein en Kritieke Helling: Stelling 2.2 stelt vast dat het maximale domein van de oplossing een eindig interval $[0, r_M)$ is. De oplossing eindigt op een punt waar de coëfficiënt van $u''$ in de ODE verdwijnt, wat overeenkomt met een kritieke helling:
   $$\lim_{r \to r_M} u'(r) = \frac{1}{\sqrt{3}}$$
   Numerieke berekeningen geven aan dat $r_M \approx 1.230$ en $u(r_M) \approx 0.228$. De fasevlakanalyse onthult dat de trajecten asymptotisch naderen tot de lijn $\theta = \pi/6$ (waar $\tan \theta = 1/\sqrt{3}$), en dat het domein eindigt bij een eindige straal.

4. Fasevlakstructuur: De studie karakteriseert het fasevlak van het bijbehorende autonome stelsel, identificeert evenwichtspunten en toont aan dat de stroming rond specifieke punten uit topologische cirkels bestaat. Er wordt aangetoond dat trajecten die vanuit de oorsprong vertrekken, eindigen op de horizontale lijn $\theta = \pi/6$ met een eindige $x$-coördinaat.

Beteekenis en Beweringen
Het artikel beweert dat de exponentiële afname van de vloeistofdichtheid de geometrie van lichamen met minimale weerstand fundamenteel verandert in vergelijking met het klassieke homogene model. Specifiek is het bestaan van gladde ($C^2$) radiale oplossingen die de symmetrieas orthogonaal snijden een nieuw resultaat. In het klassieke model is dergelijke gladheid aan de as onmogelijk voor niet-triviale oplossingen, die doorgaans een "vlak-front" schijf vereisen (een discontinuïteit in de afgeleide) om de weerstand te minimaliseren.

De auteur stelt dat deze gladde radiale oplossingen een theoretische basis vormen voor optimale vormen in gestratificeerde media. De geometrische beperking dat de helling begrensd is door $1/\sqrt{3}$ aan de rand van het domein, wordt benadrukt als een significant kenmerk voor het ontwerp van lichamen die door de atmosfeer bewegen. Het artikel concludeert dat hoewel deze $C^2$-extrema bestaan en regulier zijn, de vraag of ze globale minimizers vertegenwoordigen versus de klassieke niet-gladde minimizers een open probleem blijft voor toekomstig onderzoek. Het werk wordt gepresenteerd als een rigoureuze wiskundige verkenning van het bestaan en de eigenschappen van deze specifieke extrema, zonder nieuwe experimentele toepassingen voor te stellen of directe technische implementatie te claimen.