Barotropic-Baroclinic Splitting for Multilayer Shallow Water Models with Exchanges

Resumen Técnico: División Barotrópica-Baroclínica para Modelos de Aguas Someras Multicapa con Intercambios

Planteamiento del Problema
Los Modelos de Circulación Oceánica Global (OGCM, por sus siglas en inglés) enfrentan una severa restricción computacional: el paso de tiempo está limitado por la velocidad de la onda más rápida, específicamente las ondas de gravedad superficiales (aprox. 200 m/s), a pesar de que los fenómenos primarios de interés (estratificación y corrientes) evolucionan en escalas de tiempo mucho más lentas (aprox. 3 m/s). Aunque se han utilizado enfoques previos como la aproximación de tapa rígida o métodos de altura superficial implícitos, estos presentan limitaciones. La descomposición barotrópica-baroclínica ofrece una solución al separar las rápidas ondas de gravedad superficial (barotrópicas) de los ajustes verticales más lentos (baroclínicos). Sin embargo, los métodos de división existentes suelen tener dificultades para mantener la consistencia en marcos multicapa con intercambios de masa, preservar las propiedades de energía y asegurar la estabilidad sin perder precisión, particularmente en coordenadas que siguen el terreno.

Metodología
Los autores proponen un análisis numérico de una división barotrópica-baroclínica dentro de un marco no lineal multicapa con intercambios de masa, formulado como una división de operadores exacta. El modelo se basa en las ecuaciones de Euler sin viscosidad con un trazador pasivo en coordenadas de seguimiento del terreno ($\sigma$), asumiendo densidad constante.

La metodología central consiste en descomponer las ecuaciones de Euler promediadas por capa en dos operadores distintos:

1. Paso Barotrópico: Este paso gestiona la evolución de la altura de la superficie libre ($h$) y la velocidad promediada por la profundidad ($\bar{u}$). Se formula como un modelo de aguas someras de una sola capa bien equilibrado que incluye términos de presión y términos de fuente de topografía, pero excluye el forzamiento baroclínico. Crucialmente, este paso trata el transporte de las velocidades individuales de las capas ($u_\alpha$) y los trazadores ($T_\alpha$) a la velocidad media ($\bar{u}$).
2. Paso Baroclínico: Este paso gestiona los intercambios verticales entre capas y el ajuste de las velocidades respecto a la media. Resuelve la desviación de la velocidad media ($u_\alpha - \bar{u}$) y garantiza la restricción de que las alturas de las capas permanezcan proporcionales a la altura total ($h_\alpha = l_\alpha h$) mediante multiplicadores de Lagrange (términos de intercambio $G_{\alpha \pm 1/2}$). Este sistema es totalmente 3D pero carece de términos de presión y batimetría, lo que permite un paso de tiempo mayor.

El esquema numérico emplea una división de Lie con una estrategia de subciclos:

• El paso baroclínico se resuelve en un paso de tiempo grande $\Delta t$.
• El paso barotrópico se resuelve en el mismo intervalo utilizando $K$ subpasos pequeños $\delta t$ restringidos por la velocidad de la onda de gravedad.
• Un mecanismo de corrección asegura que la velocidad media predicha por el paso barotrópico coincida con la media de las velocidades individuales actualizadas en el paso baroclínico.

Contribuciones Clave

• Formulación de División de Operadores Exacta: El artículo reformula la división como una suma de operadores exactos ($M_L = BT + BC$), aclarando la estructura matemática del sistema multicapa con intercambios.
• Preservación de Propiedades Físicas: Los autores demuestran que la división preserva la conservación de la energía total (o disipación) y satisface una desigualdad de entropía discreta. Específicamente, el paso barotrópico satisface una igualdad de energía, mientras que el paso baroclínico satisface una desigualdad de energía debido a la naturaleza disipativa de los términos de intercambio.
• Principios de Máximo Discreto: Se demuestra que el esquema satisface principios de máximo discreto para la velocidad y las concentraciones de trazadores, asegurando que los valores permanezcan dentro de límites físicos.
• Equilibrio Bien Balanceado: Al utilizar un resolvedor bien equilibrado (reconstrucción hidrostática) para el paso barotrópico, el esquema general preserva el equilibrio de "lago en reposo" y, al extenderse con estrategias recientes de equilibrio bien balanceado lineal, mantiene el equilibrio geostrófico en mallas gruesas.
• Análisis de Hiperbolicidad: El artículo investiga la hiperbolicidad de los sistemas divididos, demostrando que el sistema barotrópico es hiperbólico y analizando los autovalores del sistema baroclínico, observando que este permanece estrictamente hiperbólico bajo condiciones específicas respecto a las desviaciones de velocidad.

Resultados
Los experimentos numéricos demuestran la eficacia del esquema propuesto:

• Eficiencia Computacional: En simulaciones de bajo número de Froude (donde la relación entre la velocidad del flujo y la velocidad de la onda de gravedad es pequeña), el método de división reduce drásticamente el costo computacional en comparación con los esquemas no divididos. La ganancia es proporcional al número de Froude.
• Precisión: A pesar de la división y los subciclos, el método mantiene una precisión comparable a la de los esquemas no divididos. Los estudios de convergencia frente a soluciones analíticas de las ecuaciones de Euler hidrostáticas muestran una convergencia de primer orden.
• Características de Difusión: En ciertos casos de prueba (por ejemplo, cavidad impulsada por el viento), se encontró que el método de división es menos difusivo que el esquema no dividido, aunque en otros (soluciones estacionarias), exhibió una difusión ligeramente mayor.
• Equilibrio Geostrófico: La aplicación de estrategias de equilibrio bien balanceado lineal en el paso barotrópico mejora significamente la resolución de las corrientes de frontera oeste en integraciones de largo plazo, capturando corrientes más estrechas y rápidas en comparación con los esquemas estándar en mallas gruesas.

Significado y Reivindicaciones
El artículo sostiene que este trabajo proporciona un análisis numérico riguroso de una división barotrópica-baroclínica que es tanto computacionalmente eficiente como físicamente consistente. El significado principal radica en la capacidad de simular la circulación oceánica a largo plazo con pasos de tiempo grandes (restringidos por ondas baroclínicas lentas en lugar de ondas de gravedad rápidas) sin sacrificar la conservación de la energía o la preservación de equilibrios críticos (lago en reposo y geostrófico). Los autores enfatizan que la estrategia de equilibrio bien balanceado utilizada en el paso barotrópico es heredada por el esquema general, permitiendo simulaciones precisas de flujos geostróficos en mallas gruesas. El método se presenta como un marco robusto tanto para aplicaciones globales como regionales donde las ondas de la superficie libre no son el foco principal.