Spectral interpolation in semi-implicit semi-Lagrangian methods for shallow water equations on the sphere

Resumen Técnico: Interpolación Espectral en Métodos Semi-Implicitos Semi-Lagrangianos para las Ecuaciones de Agua Somera en la Esfera

Enunciado del Problema
Los métodos semi-implícitos semi-lagrangianos (SISL) son el estándar para los modelos operativos globales de clima y tiempo porque permiten pasos de tiempo mayores que los permitidos por la condición de Courant-Friedrichs-Lewy (CFL). Sin embargo, existe una inconsistencia numérica significativa en los esquemas operativos actuales. Mientras que la discretización espacial de las derivadas de las variables prognósticas emplea típicamente métodos de precisión espectral (como las Transformadas de Armónicos Esféricos [SHT] o el método de la Esfera de Doble Fourier [DFS]), el paso de interpolación requerido para evaluar variables en puntos de salida fuera de la red depende de técnicas de bajo orden, como polinomios de Lagrange de producto tensorial cúbicos o quínticos.

Esta discrepancia genera una diferencia en la precisión numérica. Si bien la interpolación de bajo orden a menudo se justifica por preocupaciones sobre el costo computacional (con una complejidad de $O(n^2)$ frente a $O(n^4)$ para la interpolación espectral) y límites de error teóricos para la advección simple, existe evidencia numérica limitada en la literatura que justifique esta compensación para sistemas complejos como las Ecuaciones de Agua Somera (SWE). Además, se sabe que la interpolación de bajo orden introduce difusión numérica excesiva y errores de fase, degradando potencialmente la precisión y las propiedades de conservación a lo largo de largos tiempos de integración.

Metodología
Este trabajo presenta la primera investigación numérica de la interpolación de precisión espectral dentro de esquemas SISL para las SWE en una esfera. Los autores se basan en el modelo SISL basado en DFS desarrollado por Yoshimura [56], que utiliza una base de Fourier estable y un esquema SETTLS de dos niveles temporales.

La innovación metodológica central es la sustitución de la interpolación de Lagrange estándar por un esquema de interpolación espectral basado en la misma base de Fourier DFS utilizada para las derivadas espaciales. Para abordar la barrera del costo computacional, los autores emplean la Transformada Rápida de Fourier No Uniforme (NUFFT).

• Mecanismo de Interpolación: Las variables prognósticas (específicamente los términos $R_v$ y $R_h$ en el sistema semi-implícito) se expanden utilizando la base DFS de regularidad parcial. Estas expansiones se evalúan luego en los puntos de salida fuera de la red.
• Reducción de Complejidad: La evaluación directa de series de Fourier bivariadas en $O(n^2)$ puntos de salida requeriría $O(n^4)$ operaciones. Al reformular las series en forma de exponencial compleja y utilizar la NUFFT (específicamente la biblioteca FINUFFT), la complejidad computacional se reduce a $O(n^2 \log^2(1/\epsilon) + n^2 \log n)$, donde $\epsilon$ es la tolerancia de precisión. Esto hace que el paso de interpolación sea computacionalmente consistente con el resto del esquema espectral SISL.
• Implementación: Los autores implementaron esta interpolación espectral en dos contextos:
  1. Un modelo SISL basado en DFS (utilizando la base de Fourier estable de Yoshimura).
  2. Un modelo SISL basado en SHT (utilizando redes de Gauss-Legendre), donde el método DFS se adapta para generar interpolantes de Fourier bivariados a partir de los puntos de latitud no equidistantes.

Contribuciones y Resultados Clave
Los autores evaluaron los nuevos esquemas de interpolación espectral frente a la interpolación de Lagrange de bajo orden estándar (LAG) utilizando una serie de casos de prueba SWE estándar (Williamson TC1, TC2, TC5, TC6 y el caso de prueba Galewsky). Los resultados demuestran varias ventajas distintas:

1. Precisión y Convergencia:

   • TC1 (Advección de una campana coseno): La interpolación DFS espectral mantuvo la misma tasa de convergencia espacial que LAG, pero logró errores aproximadamente dos órdenes de magnitud menores. A lo largo de largos tiempos de integración (10 revoluciones), el método DFS redujo los errores en tres órdenes de magnitud en comparación con LAG.
   • TC2 (Flujo Geostrofico): Mientras que la interpolación LAG no logró mostrar la convergencia temporal de segundo orden esperada (probablemente debido a errores espaciales dominantes), los esquemas de interpolación DFS alcanzaron la precisión teórica de segundo orden.
   • TC5 y TC6 (Flujo sobre Montañas y Ondas de Rossby-Haurwitz): El método DFS produjo consistentemente errores menores en los campos de altura en comparación con LAG, particularmente en resoluciones más finas y pasos de tiempo más pequeños. En TC6, el método DFS preservó mejor la forma de los valles de las ondas, coincidiendo más estrechamente con las soluciones de referencia eulerianas de alto orden.

2. Propiedades de Conservación:

   • La interpolación espectral mejoró significativamente la conservación de masa y energía. En TC1, el error absoluto de masa después de 120 días fue de $1.21 \times 10^{-11}$ para DFS frente a $4.34 \times 10^{-4}$ para LAG.
   • En TC2 y TC5, el método DFS mantuvo errores estables de conservación de masa y energía a lo largo del tiempo, mientras que los errores de LAG tendieron a aumentar u oscilar más significativamente.

3. Difusión Numérica y Dinámica de Pequeña Escala:

   • Caso de Prueba Galewsky: Esta prueba involucra la evolución de un chorro no lineal hacia turbulencia. El análisis espectral del espectro de energía cinética reveló que la interpolación LAG causó una disipación rápida de energía de alto número de onda. Por el contrario, los esquemas de interpolación DFS preservaron la cascada de energía, permitiendo que el espectro de alto número de onda se acercara a la escala teórica $k^{-5/3}$ observada en flujos rotantes turbulentos. La interpolación LAG no logró capturar este comportamiento, indicando una amortiguación numérica excesiva.

4. Eficiencia Computacional:

   • Las comparaciones de tiempo real mostraron que para requisitos de baja precisión, la interpolación LAG es más rápida. Sin embargo, a medida que aumenta la precisión requerida (tolerancias de error más estrictas), el método de interpolación DFS se vuelve más eficiente, alcanzando la precisión objetivo en menos tiempo que el enfoque LAG. Esto confirma que la interpolación espectral acelerada por NUFFT es viable para aplicaciones de alta precisión.

Significado
El artículo afirma que incorporar interpolación de precisión espectral en esquemas SISL resuelve la discrepancia numérica entre la discretización espacial y la interpolación de advección. El estudio demuestra que este enfoque no es meramente una mejora teórica, sino que ofrece beneficios tangibles:

• Elimina la "brecha de precisión" entre los cálculos de derivadas espectrales y los pasos de interpolación.
• Reduce la difusión numérica, preservando así estructuras de pequeña escala y características ondulatorias a lo largo de largos períodos de integración.
• Mejora la conservación de masa y energía, lo cual es crítico para simulaciones climáticas a largo plazo.
• Logra estas mejoras sin incurrir en costos computacionales prohibitivos, gracias a la aceleración NUFFT.

Los autores concluyen que, si bien la interpolación de bajo orden ha sido un compromiso práctico, el nuevo enfoque espectral ofrece una alternativa superior para simulaciones de alta precisión y a largo plazo de la dinámica atmosférica, motivando su extensión a núcleos dinámicos tridimensionales más complejos.