Formation of Implosion Singularities in 3D Compressible Navier-Stokes-Korteweg Equation

Resumen Técnico: Formación de Singularidades de Implosión en las Ecuaciones de Navier–Stokes–Korteweg de 3D Compresibles

Planteamiento del Problema
Este artículo investiga la formación de singularidades en tiempo finito para el sistema de Navier–Stokes–Korteweg (NSK) compresible tridimensional con coeficientes de viscosidad y capilaridad dependientes de la densidad. El sistema está definido por:
$$\begin{cases} \rho_t + \text{div}(\rho u) = 0, \\ (\rho u)_t + \text{div}(\rho u \otimes u) + \nabla P = \text{div}(2\mu(\rho)D(u)) + \nabla(\lambda(\rho)\text{div} u) + \text{div} K, \end{cases}$$
donde los coeficientes satisfacen las relaciones específicas de ley de potencia:
$$\mu(\rho) = \nu\rho^\alpha, \quad \lambda(\rho) = 2\nu(\alpha - 1)\rho^\alpha, \quad \kappa(\rho) = \varepsilon^2\alpha^2\rho^{2\alpha-3}.$$
Trabajos previos de Gu–Huang–Meng–Zhou [19] y Huang–Lei–Zhou [21] establecieron soluciones fuertes globales para datos iniciales arbitrariamente grandes (lejos del vacío) cuando el exponente $\alpha$ se encuentra en un rango adecuado específico (por ejemplo, $\alpha \in (\frac{9\sqrt{3}-4\sqrt{2}}{9\sqrt{3}-2\sqrt{2}}, 1)$ para 3D). La cuestión central abordada aquí es si las soluciones existen globalmente o si estallan en tiempo finito cuando $\alpha$ cae fuera de este rango, específicamente para exponentes positivos pequeños ($\alpha < 1/2$). En este régimen, la viscosidad volumétrica efectiva se vuelve negativa ($\lambda(\rho) + \frac{2}{3}\mu(\rho) < 0$), una característica asociada con medios acústicos fuera del equilibrio y materia activa, lo que desafía las teorías estándar de existencia global.

Metodología
Los autores emplean un enfoque perturbativo basado en la construcción de perfiles de implosión autosemejantes suaves derivados de las ecuaciones de Euler compresibles. La metodología procede a través de varias etapas clave:

1. Reformulación de la Velocidad Efectiva: Para manejar el tensor de tercer orden de Korteweg, el sistema se reformula utilizando una velocidad efectiva $v = u + d\rho^{\alpha-2}\nabla\rho$ (donde $d = \nu - \sqrt{\nu^2-\varepsilon^2}$). Esta transformación reestructura el sistema en una estructura parabólica-parabólica para $(\rho, v)$, revelando un término de difusión degenerado oculto $d\Delta(\rho^\alpha)$ en la ecuación de masa.
2. Escalamiento Autosemejante: El problema se transforma en coordenadas autosemejantes $(\tau, y)$ donde $\tau = -\log(T-t)/\Lambda$ y $y = x/(T-t)^{1/\Lambda}$. En estas coordenadas, los efectos viscosos y capilares aparecen como perturbaciones que decaen exponencialmente ($e^{-f_1\tau}$) siempre que el parámetro de escalamiento $\Lambda$ satisfaga una condición específica ($f_1 > 0$).
3. Análisis de Estabilidad de las Perturbaciones: El núcleo de la prueba consiste en analizar la estabilidad no lineal del sistema alrededor de un perfil de Euler autosemejante, suave y esféricamente simétrico $(\bar{S}, \bar{V})$. Las ecuaciones de perturbación se derivan y se construye un argumento de bootstrap para controlar la solución.
4. Estimaciones de Energía y Descomposición Espectral:
   • Estimaciones de Energía Ponderadas: Se establecen estimaciones de energía ponderadas de alto orden para controlar el decaimiento de las perturbaciones. Los autores utilizan una función de peso específica $\phi$ y una jerarquía de parámetros para gestionar la interacción entre transporte y disipación.
   • Argumento de Bootstrap: Se diseña un esquema de bootstrap riguroso para cerrar las estimaciones de las derivadas de orden inferior (específicamente el límite inferior de la densidad) y las de orden superior. Una innovación crítica es el uso de estructuras tanto de transporte como de disipación para establecer un límite inferior para la densidad, superando la dificultad de que la ecuación de la densidad posee una estructura de disipación además de la de transporte.
   • Control de Modos Inestables: El operador linealizado se descompone en subespacios estables e inestables. Los componentes inestables, que corresponden a un espacio de dimensión finita, se controlan mediante la selección de los coeficientes de los datos iniciales específicos. Esta selección asegura que la solución permanezca en una variedad de codimensión finita que conduce al estallido (blow-up).

Contribuciones Clave y Resultados
El artículo demuestra la existencia de datos iniciales suaves con densidad uniformemente separada del vacío que evolucionan hacia singularidades de implosión en tiempo finito para el sistema NSK 3D cuando $\alpha$ es suficientemente pequeño ($\alpha < 1/2$).

• Teorema Principal (Teorema 1.1): Para parámetros que satisfacen ya sea:
  • (P1) $1 < \gamma < 1 + 2/\sqrt{3}$, $0 < \alpha < 1/2$, y $f_1 > 0$; o
  • (P2) $\gamma \in (1, 1 + 2/\sqrt{3}) \setminus \mathcal{J}$ (donde $\mathcal{J}$ es un conjunto excepcional contable) y $0 < \alpha < \alpha^*(\gamma)$,
    existen datos iniciales $C^\infty$ $(\rho_0, u_0)$ tales que la solución correspondiente $(\rho, u)$ existe en $[0, T)$ y desarrolla una singularidad de implosión en $t=T$.
• Características de la Singularidad: En el tiempo de estallido $T$:
  • La densidad se vuelve infinita en el origen: $\lim_{t \to T^-} \rho(t, 0) = +\infty$.
  • La velocidad efectiva $v$ es ilimitada en cada vecindad del origen.
  • La solución converge, tras el reescalamiento, al perfil de implosión autosemejante prescrito.
• Corolario 1.1: El sistema original desarrolla una singularidad en tiempo finito, caracterizada específicamente por el estallido de la densidad ($\lim_{t \to T^-} \rho(t, 0) = +\infty$). Sin embargo, no se puede concluir directamente el estallido de la velocidad física $u$ a partir del estallido de la velocidad efectiva, ya que el término de gradiente $\nabla\rho$ también puede volverse singular cerca del origen, y no se puede descartar una posible cancelación en la relación $v = u + d\rho^{\alpha-2}\nabla\rho$.

Significado
El artículo afirma complementar la teoría de existencia global existente para el sistema NSK. Mientras que resultados previos demostraron la buena existencia global para un rango específico de $\alpha$, este trabajo identifica un mecanismo distinto de estallido en tiempo finito para el régimen de "$\alpha$ pequeño". Este régimen se caracteriza por una estructura de viscosidad volumétrica efectiva negativa, un fenómeno relevante para la física fuera del equilibrio. Los resultados demuestran que incluso con datos iniciales suaves lejos del vacío, la interacción entre la viscosidad dependiente de la densidad, la capilaridad y la viscosidad volumétrica negativa puede conducir a singularidades de implosión donde los campos de densidad y velocidad efectiva mismos estallan, en lugar de solo sus derivadas (formación de choques). La construcción se basa en la estabilidad de los perfiles de Euler autosemejantes, extendiendo el mecanismo de implosión conocido para las ecuaciones de Euler al entorno viscoso y capilar de NSK.