A Generalised Curvilinear Coordinate system-based Patch Dynamics Scheme in Equation-free Multiscale Modelling

Resumo Técnico: Um Esquema de Dinâmica de Patches Baseado em um Sistema de Coordenadas Curvilíneas Generalizadas em Modelagem Multiescala Sem Equação

Enunciado do Problema
O artigo aborda os desafios computacionais inerentes à modelagem multiescala de sistemas complexos, especificamente aqueles regidos por equações de convecção-difusão-reação (CDR) lineares, heterogêneas e não estacionárias. Esquemas numéricos tradicionais frequentemente exigem simulações em níveis finos sobre grandes domínios espaciais e intervalos de tempo longos, levando a custos computacionais e requisitos de memória proibitivos. Embora métodos de multiescala "sem equação" (equation-free), como o esquema de dinâmica de patches, tenham sido desenvolvidos para unir as escalas microscópica e macroscópica simulando apenas frações do domínio espaço-temporal, as implementações existentes são amplamente restritas a domínios retangulares com grades uniformes e patches retangulares. Essa limitação dificulta sua aplicação em problemas da vida real caracterizados por geometrias irregulares, requisitos de grade não uniformes (por exemplo, para capturar camadas limite) e domínios físicos complexos.

Metodologia
Os autores propõem uma nova extensão do esquema de dinâmica de patches que integra coordenadas curvilíneas ortogonais generalizadas. A metodologia central envolve:

1. Transformação de Coordenadas: O domínio físico ($\Omega_p$), que pode ser não retangular ou exigir espaçamento de grade não uniforme, é mapeado para um domínio computacional retangular ($\Omega_c$) com uma malha uniforme via uma transformação suave e um-para-um. Essa transformação converte a equação CDR heterogênea original em uma forma onde o termo de derivada mista desaparece (condição ortogonal), facilitando a computação em uma grade estruturada.
2. Estrutura de Dinâmica de Patches: O esquema opera no domínio computacional usando uma combinação do esquema gap-tooth e integração projetiva grosseira.
   • Microssimulação: O problema microscópico é resolvido apenas dentro de pequenos e esparsos patches centrados em pontos de grade macroscópicos.
   • Elevação (Lifting) e Restrição: Um operador de elevação reconstrói as condições iniciais microscópicas dentro dos patches a partir do estado macroscópico usando interpolação polinomial bi-quadrática. Inversamente, um operador de restrição faz a média da solução microscópica sobre cada patch para atualizar o estado macroscópico.
   • Condições de Borda: Condições de contorno de Neumann são impostas nas bordas dos patches usando interpolação polinomial de Lagrange baseada em pontos de grade macroscópicos vizinhos (um stencil de 9 pontos).
   • Integração Projetiva: A derivada temporal macroscópica é estimada a partir de curtas microssimulações, e a solução é extrapolada para frente sobre um grande passo de tempo macroscópico usando o método de Euler progressivo (forward Euler).
3. Tratamento de Heterogeneidade: O esquema acomoda tensores de difusão dependentes do espaço, velocidades convectivas e termos de fonte dependentes de espaço-tempo. Ele é projetado para problemas com coeficientes suaves e não oscilatórios.

Principais Contribuições

• Estrutura Generalizada: Esta é a primeira aplicação de um sistema de coordenadas curvilíneas ortogonais generalizadas dentro do esquema de dinâmica de patches, permitindo o tratamento de tanto grades não uniformes em domínios retangulares quanto domínios físicos não retangulares.
• Patches Não Uniformes e Não Retangulares: O artigo introduz o conceito de configurações de patches não uniformes e não retangulares no domínio físico, que são mapeadas para patches retangulares uniformes no domínio computacional.
• Formulação Explícita: Uma representação explícita do esquema é fornecida para problemas CDR lineares, heterogêneos e não estacionários em duas dimensões, incluindo a derivação dos coeficientes transformados para as equações governantes.
• Implementação Algorítmica: Os autores detalham o algoritmo completo, incluindo o tratamento específico das condições de borda, elevação e extrapolação projetiva no contexto de coordenadas curvilíneas.

Resultos
O esquema proposto foi validado através de dois casos de teste numéricos distintos:

1. CDR Heterogêneo Dominado por Convecção em uma Grade Esticada:

   • Um problema com um domínio retangular apresentando altos gradientes (camadas limite) perto das fronteiras superior e direita foi resolvido.
   • Grades Esticadas: Usando um parâmetro de estiramento $\lambda = 0,1$, o esquema alcançou uma precisão significativamente maior comparado a grades uniformes ($\lambda = 0$), particularmente em regiões de alto gradiente.
   • Desempenho: O esquema de dinâmica de patches demonstrou precisão superior (erros até 13 vezes menores em casos específicos) em comparação com uma simulação de domínio total (FDS) usando o método das linhas. Também exigiu menos tempo computacional (redução de aproximadamente 40%) e memória (redução de aproximadamente 10%) utilizando menos pontos de rede.
   • Independência de Grade: As soluções mostraram-se independentes da grade em uma resolução de $20 \times 20$ grades macro com 8.500 passos de tempo.

2. Difusão Não Assimétrica em um Anel:

   • Um problema de difusão em um domínio anular (geometria não retangular) foi resolvido usando um sistema de coordenadas curvilíneas ajustado ao corpo (body-fitted).
   • Precisão: Os resultados numéricos mostraram excelente concordância com a solução analítica.
   • Convergência: O esquema demonstrou convergência conforme os níveis de grade e de tempo macroscópicos eram refinados.
   • Condições de Contorno: O método lidou com sucesso tanto com condições de contorno periódicas (na direção angular) quanto com condições de contorno de Dirichlet (na direção radial).

Significância e Alegações
O artigo afirma que o esquema proposto oferece uma alternativa robusta e eficiente para problemas multiescala complexos. Sua significância reside na:

• Flexibilidade Geométrica: Supera a limitação de domínios retangulares, permitindo a simulação precisa de sistemas com geometrias complexas e regiões localizadas de alto gradiente sem a necessidade de uma resolução de grade fina global.
• Eficiência Computacional: Ao restringir as microssimulações a pequenos patches e utilizar a integração projetiva, o esquema reduz significativamente o tempo computacional e o uso de memória, mantendo alta precisão.
• Versatilidade: O método é aplicável tanto a problemas "Tipo-A" (onde detalhes de microescala são necessários apenas em regiões específicas, como camadas limite) quanto a problemas "Tipo-B" (onde modelos de microescala são necessários em todo o domínio devido à ausência de uma descrição macroscópica).

Os autores observam que a implementação atual é um desenvolvimento inicial focado em equações CDR lineares, não estacionárias e heterogêneas com coeficientes suaves. Eles reconhecem que trabalhos futuros são necessários para estender o framework para problemas com coeficientes altamente oscilatórios, onde erros de ressonância provenientes de condições de borda artificiais podem exigir regiões de amortecimento (buffer regions) ou núcleos de média.