Barotropic-Baroclinic Splitting for Multilayer Shallow Water Models with Exchanges

Resumo Técnico: Divisão Barotrópica-Baroclínica para Modelos de Águas Rasas Multicamadas com Trocas

Enunciado do Problema
Os Modelos de Circulação Oceânica Global (OGCMs) enfrentam uma severa restrição computacional: o passo de tempo é limitado pela velocidade da onda mais rápida, especificamente as ondas de gravidade superficiais (aprox. 200 m/s), embora os principais fenômenos de interesse (estratificação e correntes) evoluam em escalas de tempo muito mais lentas (aprox. 3 m/s). Embora abordagens anteriores, como a aproximação de superfície rígida ou métodos de altura superficial implícita, tenham sido utilizadas, elas apresentam limitações. A decomposição barotrópica-baroclínica oferece uma solução ao separar as rápidas ondas de gravidade superficiais (barotrópicas) das ajustes baroclínicos mais lentos. No entanto, os métodos de divisão existentes frequentemente têm dificuldade em manter a consistência em estruturas multicamadas com trocas de massa, preservar propriedades de energia e garantir estabilidade sem perder a precisão, particularmente em coordenadas que seguem o terreno.

Metodologia
Os autores propõem uma análise numérica de uma divisão barotrópica-baroclínica dentro de um framework multicamadas não linear com trocas de massa, formulada como uma divisão de operador exata. O modelo baseia-se nas equações de Euler não viscosas com um traçador passivo em coordenadas de terreno ($\sigma$), assumindo densidade constante.

A metodologia central envolve a decomposição das equações de Euler médias por camada em dois operadores distintos:

1. Passo Barotrópico: Este passo lida com a evolução da altura da superfície livre ($h$) e da velocidade média da profundidade ($\bar{u}$). É formulado como um modelo de águas rasas de uma camada bem balanceado que inclui termos de pressão e termos de fonte de topografia, mas exclui o forçamento baroclínico. Crucialmente, este passo trata o transporte das velocidades individuais das camadas ($u_\alpha$) e dos traçadores ($T_\alpha$) na velocidade média $\bar{u}$.
2. Passo Baroclínico: Este passo gerencia as trocas verticais entre camadas e o ajuste das velocidades em relação à média. Ele resolve o desvio da velocidade média ($u_\alpha - \bar{u}$) e impõe a restrição de que as alturas das camadas permaneçam proporcionais à altura total ($h_\alpha = l_\alpha h$) via multiplicadores de Lagrange (termos de troca $G_{\alpha \pm 1/2}$). Este sistema é totalmente 3D, mas carece de termos de pressão e batimetria, permitindo um passo de tempo maior.

O esquema numérico emprega uma divisão de Lie com uma estratégia de subciclo:

• O passo baroclínico é resolvido em um grande passo de tempo $\Delta t$.
• O passo barotrópico é resolvido no mesmo intervalo usando $K$ pequenos subpassos $\delta t$ limitados pela velocidade da onda de gravidade.
• Um mecanismo de correção garante que a velocidade média prevista pelo passo barotrópico coincida com a média das velocidades atualizadas no passo baroclínico.

Principais Contribuições

• Formulação de Divisão de Operador Exata: O artigo reformula a divisão como uma soma de operador exata ($M_L = BT + BC$), esclarecendo a estrutura matemática do sistema multicamadas com trocas.
• Preservação de Propriedades Físicas: Os autores provam que a divisão preserva a conservação de energia total (ou dissipação) e satisfaz uma desigualdade de entropia discreta. Especificamente, o passo barotrópico satisfaz uma igualdade de energia, enquanto o passo baroclínico satisfaz uma desigualdade de energia devido à natureza dissipativa dos termos de troca.
• Princípios de Máximo Discreto: Demonstra-se que o esquema satisfaz princípios de máximo discreto para velocidade e concentrações de traçadores, garantindo que os valores permaneçam dentro de limites físicos.
• Bem-Balanceamento (Well-Balancing): Ao utilizar um solver bem-balanceado (reconstrução hidrostática) para o passo barotrópico, o esquema geral preserva o equilíbrio "lago em repouso" e, quando estendido com estratégias recentes de bem-balanceamento linear, mantém o equilíbrio geostrófico em malhas grosseiras.
• Análise de Hiperbolicidade: O artigo investiga a hiperbolicidade dos sistemas divididos, provando que o sistema barotrópico é hiperbólico e analisando os autovalores do sistema baroclínico, observando que ele permanece estritamente hiperbólico sob condições específicas em relação aos desvios de velocidade.

Resultados
Experimentos numéricos demonstram a eficácia do esquema proposto:

• Eficiência Computacional: Em simulações de baixo número de Froude (onde a razão entre a velocidade do fluxo e a velocidade da onda de gravidade é pequena), o método de divisão reduz drasticamente o custo computacional em comparação com esquemas não divididos. O ganho é proporcional ao número de Froude.
• Precisão: Apesar da divisão e do subciclo, o método mantém uma precisão comparável aos esquemas não divididos. Estudos de convergência contra soluções analíticas das equações de Euler hidrostáticas mostram convergência de primeira ordem.
• Características de Difusão: Em certos casos de teste (ex: cavidade impulsionada pelo vento), o método de divisão mostrou-se menos difusivo que o esquema não dividido, embora em outros (soluções estacionárias), tenha exibido uma difusão ligeiramente superior.
• Equilíbrio Geostrófico: A aplicação de estratégias de bem-balanceamento linear no passo barotrópico melhora significativamente a resolução de correntes de contorno oeste em integrações de longo prazo, capturando correntes mais estreitas e rápidas em comparação com esquemas padrão em malhas grosseiras.

Significância e Alegações
O artigo afirma que este trabalho fornece uma análise numérica rigorosa de uma divisão barotrópica-baroclínica que é simultaneamente eficiente computacionalmente e fisicamente consistente. A principal significância reside na capacidade de simular a circulação oceânica de longo prazo com grandes passos de tempo (limitados por ondas baroclínicas lentas em vez de rápidas ondas de gravidade) sem sacrificar a conservação de energia ou a preservação de equilíbrios críticos (lago em repouso e geostrófico). Os autores enfatizam que a estratégia de bem-balanceamento usada no passo barotrópico é herdada pelo esquema geral, permitindo simulações precisas de fluxos geostróficos em grades grosseiras. O método é apresentado como um framework robusto para aplicações tanto globais quanto regionais onde as ondas de superfície livre não são o foco principal.