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

Resumo Técnico: O Problema de Newton de Resistência Mínima com Densidade Não Constante

Enunciado do Problema
Este artigo investiga uma variante do problema clássico de Newton de resistência mínima, que busca a forma de um corpo sólido movendo-se a velocidade constante através de um fluido que oferece a menor resistência possível. Enquanto o modelo clássico assume um fluido homogêneo, este estudo aborda o cenário mais realista em que a densidade do fluido diminui exponencialmente com a altitude, uma condição característica da atmosfera terrestre. A densidade é modelada como $\rho(z) = \rho_0 e^{-cz}$, onde $z$ representa a altitude. O corpo move-se ao longo da coordenada $z$, e o funcional de resistência é derivado com base na lei de pressão do seno ao quadrado, assumindo ausência de atrito e que as partículas não colidem novamente com a superfície após o impacto.

O funcional de resistência resultante para uma superfície definida por um gráfico $u: \Omega \to \mathbb{R}$ é dado por:
$$E[u] = \int_{\Omega} \frac{e^{-u}}{1 + |Du|^2} \, dxdy$$
Ao contrário do funcional newtoniano clássico, a presença do peso exponencial $e^{-u}$ introduz uma mudança estrutural significativa na equação de Euler-Lagrange associada. O artigo afirma explicitamente que seu foco está na existência e regularidade de extremais clássicas $C^2$ (soluções radiais) e não na determinação de um minimizador global, observando que o problema de resistência mínima sem restrições no cenário clássico frequentemente carece de regularidade $C^2$.

Metodologia
A análise concentra-se em soluções radiais $u(r)$, onde a superfície é uma superfície de revolução em torno do eixo $z$. A equação de Euler-Lagrange para o funcional $E[u]$ reduz-se a uma equação diferencial ordinária (EDO) não linear de segunda ordem:
$$r(1 + 4u'^2 + 3u'^4) - 2u'(1 + u'^2) + 2ru''(3u'^2 - 1) = 0$$
O estudo aborda o problema de valor inicial com as condições $u(0) = 0$ e $u'(0) = 0$. Um desafio primário é a singularidade da equação em $r=0$, o que impede a aplicação direta dos teoremas clássicos de existência para EDOs.

Para superar isso, o autor emprega um argumento de ponto fixo baseado no Teorema do Ponto Fixo de Banach. A EDO é reformulada em uma equação integral envolvendo o inverso de uma função monótona específica $f(x) = x/(1+x^2)^2$. Ao definir um operador $T$ em uma bola fechada no espaço de Banach $C^1([0, R])$, o autor prova que $T$ é uma contração para um $R$ suficientemente pequeno, estabelecendo assim a existência local de uma solução $C^2$.

Para analisar o comportamento global e o domínio máximo da solução, o problema é transformado em uma análise de plano de fase. A curva é reparametrizada pelo comprimento de arco $s$, introduzindo um ângulo $\theta(s)$ tal que $u'(r) = \tan \theta$. Isso converte a EDO de segunda ordem em um sistema autônomo de equações de primeira ordem para as variáveis $x(s) = r$ e $\theta(s)$. O comportamento das soluções é estudado examinando as trajetórias no plano de fase $\Theta$, especificamente observando pontos de equilíbrio e o fluxo próximo a linhas críticas.

Contribuições e Resultados Principais

1. Existência e Regularidade Locais: O artigo prova (Teorema 2.1) que o problema de valor inicial com $u(0)=u'(0)=0$ admite uma solução $u \in C^2([0, R])$ para algum $R > 0$. Uma descoberta crucial é que a solução intercepta o eixo de rotação ortogonalmente ($u'(0)=0$). Isso contrasta agudamente com o modelo clássico de densidade constante, onde soluções radiais não podem atingir o eixo ortogonalmente a menos que a solução seja uma função constante trivial. A prova estabelece ainda que $u''(0) = 1/4$.

2. Monotonicidade e Convexidade: A análise demonstra que a inclinação $u'(r)$ é estritamente positiva e crescente no domínio máximo da solução. Consequentemente, a solução radial $u(r)$ é convexa.

3. Domínio Máximo Finito e Inclinação Crítica: O Teorema 2.2 estabelece que o domínio máximo da solução é um intervalo finito $[0, r_M)$. A solução termina em um ponto onde o coeficiente de $u''$ na EDO se anula, correspondendo a uma inclinação crítica:
   $$\lim_{r \to r_M} u'(r) = \frac{1}{\sqrt{3}}$$
   Cálculos numéricos indicam $r_M \approx 1.230$ e $u(r_M) \approx 0.228$. A análise do plano de fase revela que as trajetórias aproximam-se assintoticamente da linha $\theta = \pi/6$ (onde $\tan \theta = 1/\sqrt{3}$), e o domínio termina em um raio finito.

4. Estrutura do Plano de Fase: O estudo caracteriza o plano de fase do sistema autônomo associado, identificando pontos de equilíbrio e mostrando que o fluxo ao redor de pontos específicos consiste em círculos topológicos. Demonstra-se que as trajetórias que partem da origem terminam na linha horizontal $\theta = \pi/6$ com uma coordenada $x$ finita.

Significado e Afirmações
O artigo afirma que o decaimento exponencial da densidade do fluido altera fundamentalmente a geometria dos corpos de resistência mínima em comparação com o modelo homogêneo clássico. Especificamente, a existência de soluções radiais suaves ($C^2$) que interceptam o eixo de simetria ortogonalmente é um resultado novel. No modelo clássico, tal suavidade no eixo é impossível para soluções não triviais, que tipicamente requerem um disco de "frente plana" (uma descontinuidade na derivada) para minimizar a resistência.

O autor postula que essas soluções radiais suaves fornecem uma base teórica para formas ótimas em meios estratificados. A restrição geométrica de que a inclinação é limitada por $1/\sqrt{3}$ na borda do domínio é destacada como uma característica significativa para o projeto de corpos que atravessam a atmosfera. O artigo conclui que, embora essas extremais $C^2$ existam e sejam regulares, a questão de saber se elas representam minimizadores globais versus os minimizadores não suaves clássicos permanece um problema aberto para pesquisa futura. O trabalho é apresentado como uma investigação matemática rigorosa sobre a existência e propriedades dessas extremais específicas, sem propor novas aplicações experimentais ou reivindicar implementação de engenharia imediata.