Brinkman's term for the steady-state Navier-Stokes equations with prescribed flux rate or pressure drop in perforated pipes

Resumen Técnico: Homogeneización de las ecuaciones de Navier-Stokes en estado estacionario en tuberías perforadas con condiciones de contorno mixtas

Planteamiento del problema
Este trabajo investiga el comportamiento asintótico del movimiento estacionario de un fluido viscoso e incompresible dentro de una tubería distorsionada que contiene un gran número de pequeñas partículas sólidas. El dominio $\Omega_\varepsilon$ es un subconjunto perforado de un dominio acotado y admisible $\Omega \subset \mathbb{R}^3$, obtenido mediante la eliminación de $N(\varepsilon)$ partículas de diámetro $\varepsilon^3$ separadas por una distancia mutua $\varepsilon$. Esta configuración corresponde al régimen crítico ($\alpha = 3$) en la teoría de homogeneización, donde el tamaño de la partícula escala con el cubo de la distancia entre partículas.

El flujo del fluido está gobernado por las ecuaciones de Navier-Stokes en estado estacionario con condiciones de contorno mixtas. Específicamente:

1. Se imponen condiciones de no deslizamiento en las paredes físicas de la tubería ($\Gamma_W$) y en las superficies de las partículas ($\partial K_n^\varepsilon$), donde las partículas pueden moverse con velocidades constantes prescritas $\mu_n^\varepsilon$.
2. Se aplican condiciones de presión de Bernoulli en la entrada ($\Gamma_I$) y la salida ($\Gamma_O$). Estas involucran la derivada normal de la velocidad y la presión de Bernoulli ($\Phi = p + \frac{1}{2}|u|^2$).
3. Se analizan dos escenarios distintos:
   • Flujo neto prescrito: El flujo total $F$ a través de la salida es fijo, mientras que la caída de presión es una constante desconocida determinada por la solución.
   • Caída de presión prescrita: La caída de presión entre la entrada y la salida es fija, mientras que el flujo resultante es una cantidad desconocida.

El principal desafío matemático radica en derivar las ecuaciones efectivas (homogeneizadas) para el límite cuando $\varepsilon \to 0$ sin imponer supuestos de pequeñez sobre la fuerza externa, la tasa de flujo o la caída de presión. Esto se complica debido a que los estándares estimaciones de energía para las ecuaciones de Navier-Stokes a menudo no logran proporcionar cotas uniformes en estos problemas de valores en la frontera mixtos sin restricciones adicionales.

Metodología
Los autores emplean el método de la energía dentro del marco de la teoría de homogeneización, siguiendo la estrategia general de Tartar. La demostración procede en tres etapas principales:

1. Cotas uniformes (estimaciones independientes de $\varepsilon$):
   Establecer cotas uniformes para las soluciones $(u_\varepsilon, \Phi_\varepsilon)$ es el paso más crítico y novedoso. A diferencia de los casos estándar donde las simples estimaciones de energía bastan, los autores no pueden confiar únicamente en probar la ecuación de cantidad de movimiento con la propia solución debido a las condiciones de contorno mixtas y la caída de presión desconocida.

   • En su lugar, utilizan un argumento de contradicción basado en la ley de Bernoulli para soluciones de las ecuaciones de Euler estacionarias.
   • Suponen que las normas de las soluciones crecen indefinidamente cuando $\varepsilon \to 0$, normalizan la secuencia y muestran que el límite satisface las ecuaciones de Euler.
   • Al analizar el comportamiento de la frontera de la presión de Bernoulli límite en las paredes de la tubería (utilizando resultados de regularidad local y la geometría específica de la entrada/salida), demuestran que la presión límite debe ser constante en las paredes. Esto conduce a una contradicción con las condiciones de flujo o de caída de presión, probando así la acotación uniforme de la secuencia.
   • Para el problema de flujo prescrito, se realiza un supuesto geomético adicional de que las partículas no se acumulan cerca de los cilindros de entrada y salida para facilitar el análisis de la frontera de la presión límite.

2. Construcción de correctores y operadores de restricción:
   Para manejar el dominio perforado, los autores construyen correctores de Stokes $(\omega_{n,1}^\varepsilon, \omega_{n,2}^\varepsilon)$ en el exterior de las partículas y en las regiones anulares que las rodean. Estos correctores se utilizan para definir un operador de restricción $R_\varepsilon$.

   • A diferencia de trabajos previos en el régimen subcrítico donde bastaban los correctores capacitarios escalares, el régimen crítico requiere estos correctores de Stokes de valor vectorial para dar cuenta de la capacidad relativa no nula del enjambre de partículas.
   • El operador de restricción $R_\varepsilon(\phi)$ mapea las funciones de prueba del dominio homogeneizado $\Omega$ al dominio perforado $\Omega_\varepsilon$ modificándolas cerca de las partículas para satisfacer las condiciones de contorno cero, preservando al mismo tiempo la propiedad de divergencia nula.

3. Paso al límite:
   Utilizando las cotas uniformes y el operador de restricción, los autores pasan al límite en la formulación débil de las ecuaciones de Navier-Stokes.

   • El término convectivo no lineal converge mediante convergencia fuerte en $L^4$ (garantizada por las cotas uniformes y las inclusiones compactas).
   • La interacción entre el fluido y las partículas genera un término adicional en la ecuación límite. Este término surge del comportamiento límite de la suma de las fuerzas de arrastre ejercidas por las partículas, caracterizado por las matrices de resistencia de Stokes $G_n$.

Contribuciones clave y resultados
El artículo establece dos teoremas principales (Teoremas 2.4 y 2.5) que describen el límite homogeneizado del sistema:

• Emergencia del término de Brinkman: Las ecuaciones homogeneizadas en el dominio límite $\Omega$ son las ecuaciones de Navier-Stokes aumentadas por un término de Brinkman. La ecuación efectiva toma la forma:
  $$-\Delta u + (\nabla \times u) \times u + \nabla \Phi + Gu - J = f$$
  donde:

  • $Gu$ representa un término de fricción proporcional a la velocidad del fluido $u$, que surge del arrastre colectivo de las partículas (siendo $G$ el límite de la suma escalada de las matrices de resistencia de Stokes).
  • $J$ es un término fuente que surge del flujo límite de las propias partículas en movimiento.
  • Las condiciones de contorno en la entrada y la salida siguen siendo de tipo mixto involucrando la presión de Bernoulli.

• Extensión al régimen crítico con datos no nulos: Los resultados extienden los hallazgos de homogeneización previos (específicamente en [59] para el régimen subcrítico) al régimen crítico ($\alpha=3$). Crucialmente, estos resultados se mantienen sin ningún supuesto de pequeñez sobre la fuerza externa, la tasa de flujo o la caída de presión. Este es un avance significativo, ya que muchos resultados de homogeneización para Navier-Stokes dependen de supuestos de datos pequeños para asegurar la unicidad o la existencia de cotas.

• Manejo de partículas en movimiento: El marco permite acomodar partículas con velocidades constantes no nulas $\mu_n^\varepsilon$, siempre que su energía cinética esté uniformemente acotada (Supuesto 2.11). Esto permite un modelado más realista de las interacciones fluido-partícula en comparación con los modelos de partículas fijas.

Significancia
El artículo afirma que su contribución más significativa es la derivación rigurosa de la ley de Brinkman para las ecuaciones de Navier-Stokes en estado estacionario en el régimen crítico bajo datos no restringidos. Los autores enfatizan que la derivación de las cotas uniformes mediante el argumento de contradicción que involucra la ley de Bernoulli es un asunto delicado, particularmente porque la estructura compleja de las perforaciones debe tenerse en cuenta para asegurar que el argumento se mantenga al nivel de $\varepsilon$.

Al construir exitosamente un verdadero operador de restricción que no modifica las funciones de prueba que se anulan en las partículas y al manejar la convergencia del término de presión en un entorno polidisperso con velocidades de partículas no nulas, este trabajo proporciona una base matemática robusta para modelar el flujo de fluidos en tuberías perforadas donde la densidad de partículas es lo suficientemente alta como para inducir un arrastre de tipo Brinkman, pero no tanto como para resultar en la ley de Darcy. Los resultados son directamente aplicables a campos como la hemodinámica, la ingeniería de petróleo y la aerodinámica, donde el flujo de fluido a través de arreglos de obstáculos sólidos es un fenómeno fundamental.