Coupled Continuous-Discontinuous Galerkin Finite Element Solver for Compound Flood Simulations

Este artigo apresenta um novo solver numérico que acopla métodos de Galerkin Contínuo e Descontínuo no modelo ADCIRC para simular com precisão e conservação os efeitos de inundações compostas resultantes da interação entre escoamento superficial e marés de tempestade, como observado no Furacão Harvey.

O essencial

Imagine que você está tentando prever como uma enchente vai se comportar quando duas forças da natureza se encontram: a água do mar subindo (maré de tempestade) e a chuva forte caindo do céu (escoamento superficial).

No passado, os cientistas tratavam esses dois eventos como se fossem vizinhos que nunca se falavam. Eles calculavam a maré, calculavam a chuva e depois apenas "somavam" os resultados, como se colocassem duas camadas de tinta uma sobre a outra. O problema é que, na vida real, a água da chuva e a água do mar se misturam, empurram e interagem de formas complexas. Somar os dois não dá o quadro completo.

Este artigo apresenta uma nova ferramenta computacional, uma espécie de "super simulador" chamado ADCIRC, que foi atualizado para entender essa interação complexa. Aqui está a explicação simplificada:

1. O Problema: A "Soma" Não Funciona

Pense em uma enchente composta como uma festa onde dois grupos de pessoas chegam ao mesmo tempo. Se você apenas contar quantas pessoas o Grupo A trouxe e quantas o Grupo B trouxe, você sabe o total, mas não sabe como eles vão se misturar na sala, quem vai empurrar quem, ou onde vai ficar mais apertado.

• O jeito antigo: Somar a maré + somar a chuva.
• O jeito novo: Simular como a chuva cai, corre pelas ruas, encontra o mar e como tudo isso se move junto em tempo real.

2. A Solução: Uma Equipe de Especialistas (DG e CG)

Para resolver isso, os autores criaram uma equipe de dois especialistas que trabalham juntos, cada um fazendo o que faz de melhor:

• O Especialista em "Contas" (Método DG - Galerkin Descontínuo):
  Imagine que este especialista é um contador rigoroso que não deixa nem uma gota de água escapar. Ele divide o mapa em pequenos quadrados (como um tabuleiro de xadrez) e garante que, em cada quadrado, a quantidade de água que entra é exatamente igual à que sai mais a que ficou. Ele é perfeito para lidar com a chuva e garantir que a conservação de massa (a água não desaparece magicamente) seja exata.

  • Analogia: É como ter um zelador em cada sala de um prédio que conta exatamente quantas pessoas entram e saem, garantindo que ninguém se perca.

• O Especialista em "Fluxo" (Método CG - Galerkin Contínuo):
  Este especialista é um engenheiro de tráfego que olha para o todo. Ele garante que a água se mova de forma suave e rápida entre os quadrados, sem criar "buracos" ou oscilações estranhas no mapa. Ele é muito eficiente para calcular a velocidade da água (momento).

  • Analogia: É como um maestro que garante que a orquestra toque em harmonia, sem que um instrumento desafine o outro.

A Magia: Eles acoplaram esses dois especialistas. O "Contador" cuida da água que entra (chuva) e garante que nada se perca. O "Engenheiro de Tráfego" cuida de como essa água se move. Juntos, eles simulam a enchente com muito mais precisão do que antes.

3. O Desafio da "Terra Seca" (Wetting and Drying)

Um dos maiores problemas em simular enchentes é lidar com o solo que está seco e que, de repente, fica molhado.

• O problema: Se a água da chuva começa a cair em um lugar seco, o computador precisa saber que aquele lugar agora é um rio. Se a água recuar, ele precisa saber que virou terra seca novamente.
• A solução: O novo sistema tem um "sensor de umidade" muito inteligente. Ele verifica a cada segundo se um pedaço de terra está molhado ou seco. Se a chuva cair em um lugar seco, o sistema "acorda" aquele pedaço de terra e começa a calcular o fluxo. Se a água secar, ele "desliga" o cálculo ali para economizar tempo e evitar erros.

4. A Prova de Fogo: Furacões Reais

Os autores testaram esse novo sistema com dois furacões famosos que atingiram os EUA:

• Furacão Harvey (2017): Este foi o caso perfeito para testar a "enchente composta". Harvey trouxe uma chuva histórica que inundou o interior, enquanto a maré subia no litoral. O novo modelo conseguiu prever com muito mais precisão como a chuva e a maré se juntaram para causar estragos, algo que os modelos antigos subestimavam.
• Furacão Ike (2008): Um caso de maré de tempestade muito forte. O novo modelo funcionou tão bem quanto os antigos para este caso, provando que a nova tecnologia não "quebrou" o que já funcionava.

5. Por que isso importa?

Antes, se você quisesse prever uma enchente causada por chuva e maré juntas, teria que fazer duas simulações separadas e tentar adivinhar como elas se somariam. Isso poderia levar a erros graves em planos de evacuação ou construção de diques.

Com este novo modelo:

1. Precisão: Eles veem a água real, não apenas uma soma teórica.
2. Segurança: Ajuda a prever melhor onde a água vai subir, salvando vidas e propriedades.
3. Eficiência: Eles conseguem fazer isso sem deixar o computador ficar lento demais, usando a força de supercomputadores de forma inteligente.

Em resumo: Os autores criaram um "casamento perfeito" entre dois métodos matemáticos. Um garante que a água não suma (contabilidade rigorosa), e o outro garante que a água flua suavemente (engenharia de tráfego). Juntos, eles nos dão uma visão muito mais clara e realista de como as enchentes compostas se formam, especialmente durante furacões.

Resumo técnico

Resumo Técnico: Solver Acoplado Contínuo-Descontínuo de Galerkin para Simulações de Inundações Compostas

1. Problema e Contexto

Ciclones tropicais recentes, como o Furacão Harvey (2017), demonstraram que a interação entre escoamento superficial (runoff) causado por chuvas intensas e a maré de tempestade (storm surge) pode amplificar drasticamente as inundações. Esses eventos de "inundação composta" não podem ser modelados corretamente pela simples superposição das fontes individuais.

O modelo numérico padrão para simulações de maré de tempestade, o ADCIRC (Advanced CIRCulation Model), utiliza o método de Galerkin Contínuo (CG) para discretização espacial. Embora eficiente, o método CG não garante conservação de massa local estrita. Além disso, a formulação do ADCIRC utiliza uma Equação de Continuidade de Onda Generalizada (GWCE) em vez da equação de continuidade primitiva, o que dificulta a incorporação direta e estável de fontes de massa variáveis no tempo e no espaço, como a precipitação.

2. Metodologia

Os autores desenvolveram um novo modelo numérico que acopla espaços de aproximação Contínuos (CG) e Descontínuos (DG) para resolver as Equações de Água Profunda (Shallow Water Equations - SWE) modificadas.

• Acoplamento Híbrido (DG-CG):
  • Equação de Continuidade: Discretizada utilizando o método Galerkin Descontínuo (DG). Isso garante a conservação local exata de massa, permitindo a inclusão direta de termos fonte (chuva) na equação sem induzir instabilidades numéricas significativas.
  • Equações de Momento: Discretizadas utilizando o método Galerkin Contínuo (CG) (herdado do ADCIRC). Isso preserva a eficiência computacional do modelo original para o cálculo de velocidades, mantendo o custo de solução baixo.
• Incorporação de Precipitação:
  • A equação de continuidade foi modificada para incluir um termo fonte $R$ (precipitação) que pode variar espacial e temporalmente.
  • O modelo suporta tanto modelos paramétricos de chuva (ex: R-CLIPER, baseado em parâmetros do ciclone) quanto dados observados/interpolados (formato GRIB2).
• Esquema de Molhamento e Secagem (Wetting and Drying):
  • Foi desenvolvido um esquema híbrido para lidar com regiões secas ($H=0$) e molhadas ($H>0$).
  • O estado de molhamento é gerenciado principalmente através das variáveis nodais do CG, com ajustes do DG para garantir a conservação de massa nos elementos. O algoritmo redistribui a massa e define velocidades zero em elementos parcialmente secos para manter a estabilidade.
• Implementação:
  • O solver foi integrado ao código ADCIRC.
  • Utiliza bases de polinômios lineares (Dubiner) para o DG e funções de base nodais para o CG.
  • O acoplamento entre os espaços descontínuo (elevação $\zeta$) e contínuo (velocidade $u$) é feito através de projeções e médias ponderadas por área para garantir a transferência de dados entre as malhas.

3. Contribuições Principais

1. Novo Esquema Numérico: Desenvolvimento de um solver SWE que combina a conservação local exata do DG (para elevação da superfície) com a eficiência do CG (para momento), otimizado para inundações compostas.
2. Modelagem de Chuva Direta: Capacidade de injetar chuva como um termo fonte na equação de continuidade primitiva, permitindo simulações realistas de inundações compostas em uma única malha computacional.
3. Integração no ADCIRC: Criação de uma estrutura modular para incorporar essa nova abordagem ao ecossistema existente do ADCIRC, minimizando o acoplamento complexo.
4. Validação Extensiva: Demonstração da robustez do método através de testes de convergência, casos de laboratório e simulações de grandes escalas de furacões reais.

4. Resultados e Validação

Os autores realizaram uma série de testes numéricos para validar o modelo:

• Testes de Convergência (Lynch e Gray): O esquema DG-CG demonstrou taxas de convergência ótimas (próximas a 2) para a elevação da superfície, comparáveis ao CG puro, apesar da natureza descontínua.
• Casos de Referência (Shinnecock Inlet, Quarter Annular): Os resultados mostraram diferenças mínimas em relação ao ADCIRC padrão em regiões estáveis, mas o DG-CG mostrou-se mais robusto em pontos de estrangulamento (choke points) onde o escoamento é altamente advectivo.
• Conservação de Massa e Chuva (Rain on a Hill): Em um teste de "lago em repouso" modificado com chuva, o modelo conservou a massa total com um erro relativo de apenas 0,07%, demonstrando a eficácia do termo fonte e do esquema de molhamento/secagem.
• Rio Neches (Hurricane Harvey): A comparação com dados observados e o modelo HEC-RAS mostrou que o solver DG-CG superou o ADCIRC padrão (CG) na previsão de picos de inundação, especialmente em locais onde o CG subestimava os níveis de água.
• Furacão Harvey (2017) - Escala Regional:
  • Simulações em grande escala (Golfo do México) mostraram que a inclusão da chuva no modelo DG-CG elevou os níveis de água após o pico da maré de tempestade, alinhando-se melhor com as marcas de água alta (High Water Marks) observadas.
  • O modelo capturou o efeito composto (maré + chuva) que o modelo CG puro, sem chuva, não conseguiu reproduzir.
• Desempenho Computacional:
  • Em execução serial, o solver DG-CG é cerca de 18% mais lento que o CG puro, devido ao custo adicional da integração de bordas e limitadores de inclinação no DG.
  • No entanto, o modelo escala bem em paralelo (escalabilidade forte), com a comunicação MPI sendo eficiente, tornando-o viável para simulações de grande escala em supercomputadores.

5. Significância e Conclusão

Este trabalho representa um avanço significativo na modelagem de inundações costeiras e compostas. Ao resolver a limitação de conservação de massa do método CG tradicional e permitir a integração direta de precipitação, o modelo DG-CG oferece uma ferramenta mais precisa para prever os impactos de ciclones tropicais que envolvem tanto marés de tempestade quanto inundações fluviais/pluviais.

A capacidade de usar modelos paramétricos de chuva (úteis para previsão) e dados observados (para hindcasting) torna o modelo versátil para aplicações de gestão de riscos e planejamento de resiliência costeira. Embora existam desafios remanescentes na estabilidade do esquema de molhamento/secagem em malhas extremamente refinadas e na inclusão de processos hidrológicos adicionais (infiltração, evapotranspiração), a abordagem proposta estabelece uma base sólida para a próxima geração de modelos de inundação acoplados.