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

Resumen Técnico: El Problema de Newton de Resistencia Mínima con Densidad No Constante

Enunciado del Problema
Este artículo investiga una variante del problema clásico de resistencia mínima de Newton, que busca la forma de un cuerpo sólido que se mueve a velocidad constante a través de un fluido ofreciendo la menor resistencia posible. Mientras que el modelo clásico asume un fluido homogéneo, este estudio aborda el escenario más realista donde la densidad del fluido disminuye exponencialmente con la altitud, una condición característica de la atmósfera terrestre. La densidad se modela como $\rho(z) = \rho_0 e^{-cz}$, donde $z$ representa la altitud. El cuerpo se mueve a lo largo de la coordenada $z$, y el funcional de resistencia se deriva basándose en la ley de presión del seno al cuadrado, asumiendo ausencia de fricción y que las partículas no vuelven a colisionar con la superficie tras el impacto.

El funcional de resistencia resultante para una superficie definida por un gráfico $u: \Omega \to \mathbb{R}$ viene dado por:
$$E[u] = \int_{\Omega} \frac{e^{-u}}{1 + |Du|^2} \, dxdy$$
A diferencia del funcional clásico de Newton, la presencia del peso exponencial $e^{-u}$ introduce un cambio estructural significativo en la ecuación de Euler-Lagrange asociada. El artículo establece explícitamente que su enfoque reside en la existencia y regularidad de extremales clásicas $C^2$ (soluciones radiales) en lugar de la determinación de un minimizador global, señalando que el problema de resistencia mínima sin restricciones en el contexto clásico a menudo carece de regularidad $C^2$.

Metodología
El análisis se centra en soluciones radiales $u(r)$, donde la superficie es una revolución alrededor del eje $z$. La ecuación de Euler-Lagrange para el funcional $E[u]$ se reduce a una ecuación diferencial ordinaria (EDO) no lineal de segundo orden:
$$r(1 + 4u'^2 + 3u'^4) - 2u'(1 + u'^2) + 2ru''(3u'^2 - 1) = 0$$
El estudio aborda el problema de valor inicial con las condiciones $u(0) = 0$ y $u'(0) = 0$. Un desafío principal es la singularidad de la ecuación en $r=0$, lo que impide la aplicación directa de los teoremas clásicos de existencia para EDO.

Para superar esto, el autor emplea un argumento de punto fijo basado en el Teorema del Punto Fijo de Banach. La EDO se reformula en una ecuación integral que involucra la inversa de una función monótona específica $f(x) = x/(1+x^2)^2$. Al definir un operador $T$ sobre una bola cerrada en el espacio de Banach $C^1([0, R])$, el autor demuestra que $T$ es una contracción para un $R$ suficientemente pequeño, estableciendo así la existencia local de una solución $C^2$.

Para analizar el comportamiento global y el dominio máximo de la solución, el problema se transforma en un análisis del plano de fases. La curva se reparametriza mediante la longitud de arco $s$, introduciendo un ángulo $\theta(s)$ tal que $u'(r) = \tan \theta$. Esto convierte la EDO de segundo orden en un sistema autónomo de ecuaciones de primer orden para las variables $x(s) = r$ y $\theta(s)$. El comportamiento de las soluciones se estudia examinando las trayectorias en el plano de fases $\Theta$, específicamente observando los puntos de equilibrio y el flujo cerca de las líneas críticas.

Contribuciones y Resultados Clave

1. Existencia Local y Regularidad: El artículo demuestra (Teorema 2.1) que el problema de valor inicial con $u(0)=u'(0)=0$ admite una solución $u \in C^2([0, R])$ para algún $R > 0$. Un hallazgo crucial es que la solución intersecta el eje de rotación ortogonalmente ($u'(0)=0$). Esto contrasta marcadamente con el modelo clásico de densidad constante, donde las soluciones radiales no pueden alcanzar el eje ortogonalmente a menos que la solución sea una función constante trivial. La prueba establece además que $u''(0) = 1/4$.

2. Monotonía y Convexidad: El análisis demuestra que la pendiente $u'(r)$ es estrictamente positiva y creciente en el dominio máximo de la solución. En consecuencia, la solución radial $u(r)$ es convexa.

3. Dominio Máximo Finito y Pendiente Crítica: El Teorema 2.2 establece que el dominio máximo de la solución es un intervalo finito $[0, r_M)$. La solución termina en un punto donde el coeficiente de $u''$ en la EDO se anula, correspondiente a una pendiente crítica:
   $$\lim_{r \to r_M} u'(r) = \frac{1}{\sqrt{3}}$$
   Cálculos numéricos indican $r_M \approx 1.230$ y $u(r_M) \approx 0.228$. El análisis del plano de fases revela que las trayectorias se aproximan asintóticamente a la línea $\theta = \pi/6$ (donde $\tan \theta = 1/\sqrt{3}$), y el dominio termina en un radio finito.

4. Estructura del Plano de Fases: El estudio caracteriza el plano de fases del sistema autónomo asociado, identificando puntos de equilibrio y mostrando que el flujo alrededor de puntos específicos consiste en círculos topológicos. Se demuestra que las trayectorias que parten del origen terminan en la línea horizontal $\theta = \pi/6$ con una coordenada $x$ finita.

Significado y Afirmaciones
El artículo afirma que el decaimiento exponencial de la densidad del fluido altera fundamentalmente la geometría de los cuerpos de resistencia mínima en comparación con el modelo homogéneo clásico. Específicamente, la existencia de soluciones radiales suaves ($C^2$) que intersectan el eje de simetría ortogonalmente es un resultado novedoso. En el modelo clásico, tal suavidad en el eje es imposible para soluciones no triviales, las cuales típicamente requieren un disco de "frente plano" (una discontinuidad en la derivada) para minimizar la resistencia.

El autor postula que estas soluciones radiales suaves proporcionan una base teórica para formas óptimas en medios estratificados. La restricción geométrica de que la pendiente está acotada por $1/\sqrt{3}$ en el borde del dominio se destaca como una característica significativa para el diseño de cuerpos que atraviesan la atmósfera. El artículo concluye que, aunque estas extremales $C^2$ existen y son regulares, la cuestión de si representan minimizadores globales frente a los minimizadores no suaves clásicos permanece como un problema abierto para futuras investigaciones. El trabajo se presenta como una investigación matemática rigurosa sobre la existencia y propiedades de estas extremales específicas, sin proponer nuevas aplicaciones experimentales ni afirmar una implementación ingenieril inmediata.