A carrier-wave factored one-way Navier--Stokes method for boundary-layer instability modelling

Resumen Técnico: M-OWNS para el Modelado de la Inestabilidad de la Capa Límite

Planteamiento del Problema
Predecir la transición laminar-turbulenta requiere rastrear la evolución de las perturbaciones mediante la amplificación en flujos base en desarrollo. Mientras que el análisis de autovalores locales (LST) es computacionalmente económico pero limitado por las suposiciones de flujo paralelo, y los resolvedores de Navier–Stokes armónicos linealizados (LHNS) globales son generales pero prohibitivamente costosos debido a las restricciones de memoria e inversión de matrices, los métodos de marcha espacial ofrecen un punto medio. Sin embargo, la marcha espacial estándar enfrenta una dificultad fundamental: el operador de Navier–Stokes linealizado soporta modos acústicos que se propagan hacia arriba (upstream), lo que puede desestabilizar una marcha descendente (downstream) ingenua.

Las soluciones existentes presentan compensaciones (trade-offs):

• Ecuaciones de Estabilidad Parabolizadas (PSE): Utilizan el factor de onda portadora para desplazar el espectro descendente hacia el origen, permitiendo una resolución más gruesa. Sin embargo, dependen de estrategias de regularización (por ejemplo, omitir gradientes de presión en la dirección de la corriente o amortiguamiento implícito) para suprimir los modos ascendentes, lo que introduce errores característicos. Además, la PSE sufre de una restricción de tamaño de paso mínimo; refinar el paso para capturar perturbaciones multimodales o no modales a menudo vuelve a capturar los modos ascendentes inestables, desestabilizando la marcha.
• Navier–Stokes de Una Sola Vía (OWNS): Elimina por completo los modos ascendentes mediante una proyección recursiva (OWNS-R), reteniendo el gradiente de presión en la dirección de la corriente y evitando los errores de regularización de la PSE. Sin embargo, la OWNS no emplea el factor de onda portadora. En consecuencia, debe resolver los autovalores en sus números de onda completos desde el origen, lo que requiere una resolución fina en la dirección de la corriente. Esta resolución fina amplía la ventana espectral del integrador, demandando una proyección más compleja para suprimir un rango más amplio de modos ascendentes, lo cual está limitado por la acumulación de error numérico en los parámetros de la recursión.

Metodología: M-OWNS
El artículo presenta M-OWNS, un método de marcha espacial que integra el factor de onda portadora de la PSE en el marco de la recursión de Navier–Stokes de una sola vía (OWNS-R).

1. Factor de Onda Portadora: El estado de la perturbación $\phi$ se descompone en una envolvente de variación lenta $\tilde{\phi}$ y una onda portadora de variación rápida $\exp(i\Theta_0)$, donde $d\Theta_0/dx = \alpha_0$. Esto desplaza el espectro del operador espacial por un escalar $\alpha_0$.
2. Proyección de Una Sola Vía: El operador de proyección $P$ de la OWNS-R se aplica al sistema factorizado. Los autores demuestran que el mismo operador de proyección que paraboliza el sistema no factorizado también paraboliza el sistema de onda portadora factorizado, eliminando los modos que se propagan hacia arriba sin requerir la modificación de los parámetros de proyección.
3. Selección de la Portadora ($\alpha_0$):
   • Variante Iterativa: $\alpha_0$ se actualiza en cada estación utilizando el cociente de Rayleigh, convergiendo hacia el autovalor dominante (similar a la PSE).
   • Variante de Portadora Fija: $\alpha_0$ se establece como un valor representativo derivado del flujo base (por ejemplo, el extremo de la rama de vorticidad/entropía $\alpha_1 = (\omega - \beta W_\infty)/U_\infty$). Esta variante no requiere iteraciones y mantiene el mismo costo por paso que la OWNS no factorizada.

Contribuciones Clave

• Ventaja de Resolución Espectral: Un análisis de primer orden demuestra que, para cualquier autofunción excitada cuyo autovalor se encuentre más cerca del número de onda de la portadora $\alpha_0$ que del origen, el sistema factorizado resuelve el modo con un tamaño de paso en la dirección de la corriente más grueso que el sistema no factorizado. El factor desplaza efectivamente el "disco de resolución" del integrador desde el origen hacia el grupo de los modos excitados.
• Desacoplamiento de Resolución y Regularización: Al combinar el factor con la proyección de una sola vía, M-OWNS logra los beneficios de resolución de la PSE (pasos más gruesos) sin las restricciones de tamaño de paso o los errores de regularización (por ejemplo, gradientes de presión omitidos) que afectan a la PSE.
• Robustez en Regímenes Multimodales: Se muestra que el método es efectivo tanto para perturbaciones de un solo modo (donde la iteración es óptima) como para perturbaciones multimodales o de banda ancha (donde se prefiere una portadora fija para evitar problemas de convergencia del cociente de Rayleigh asociados con autofunciones no ortogonales).

Resultados y Validación
El método fue validado en capas límite de placa plana incomprimibles y subsónicas, perturbaciones de flujo cruzado tridimensional sobre un cilindro con barrido y una capa límite hipersónica de Mach 4.5. Se probaron cuatro configuraciones de excitación: excitación de entrada de autofunción, succión/inyección en la pared, excitación de flujo libre multimodal y excitación de entrada aleatorizada.

• Eficiencia de Resolución: En todos los casos de prueba, M-OWNS capturó amplitudes de perturbación, campos de radiación acústica y secuencias de sincronización modal con resoluciones en la dirección de la corriente más gruesas que la OWNS no factorizada.
  • Para el caso del cilindro con barrido, la OWNS no factorizada requirió $N_x \approx 8000$ pasos para la convergencia, mientras que la M-OWNS de portadora fija convergió con $N_x \approx 1000$ (una reducción de 8 veces), y la M-OWNS iterativa convergió con $N_x \approx 150$ (una reducción de 47 veces).
  • Para el caso de la entrada aleatorizada hipersónica (la prueba multimodal más estricta), la OWHS no factorizada requirió $N_x \approx 6000$, mientras que la M-OWNS de portadora fija convergió con $N_x \approx 3000$ (una reducción de 2 veces).
• Exactitud: Los resultados de M-OWNS fueron validados contra PSE, OWNS no factorizada y resolvedores LHNS totalmente elípticos. En el caso de succión/inyección en la pared, los resultados de M-OWNS coincidieron con la referencia LHNS, mientras que la OWNS no factorizada a la misma resolución subestimó las amplitudes.
• Costo: El costo computacional total (total de resoluciones lineales) se redujo por factores de dos a ocho en los casos de prueba. La variante de portadora fija logró estas reducciones con costos por paso idénticos a la OWNS no factorizada, ya que elimina la sobrecarga iterativa.

Significancia
El artículo afirma que M-OWNS ofrece una ventaja distintiva de resolución numérica y eficiencia en el modelado de la evolución del estado de las perturbaciones e inestabilidades. Al desplazar la región espectral resuelta hacia el grupo de autovalores físicamente relevantes, el método permite pasos en la dirección de la corriente significativamente más gruesos sin sacrificar la exactitud ni introducir los errores de regularización inherentes a la PSE. La variante de portadora fija se destaca como particularmente significativa, ya que logra estas ganancias sin aumentar el costo computacional por paso ni requerir convergencia iterativa, lo que la convierte en una alternativa robusta tanto para problemas de inestabilidad de un solo modo como para problemas de inestabilidad complejos y multimodales. Los autores señalan que el método es aplicable a flujos generales donde el espectro excitado se agrupa cerca de un número de onda específico, proporcionando una solución práctica a las restricciones de resolución de la OWNS estándar.