Update-Magnitude State Redistribution (UM-SRD): A Shut-off Extension of Weighted SRD for Cut-Cell Methods

Resumo Técnico: Redistribuição de Estado por Magnitude de Atualização (UM-SRD)

Enunciado do Problema
Os métodos de células cortadas permitem que esquemas de volumes finitos lidem com geometrias complexas sem malhas adaptadas ao corpo, mas introduzem desafios severos de estabilidade devido a células cortadas arbitrariamente pequenas. Os métodos de Redistribuição de Estado (SRD), especificamente o framework de SRD ponderado desenvolvido por Berger & Giuliani (2024), abordam esses problemas de estabilidade fundindo células pequenas em vizinhanças, evoluindo médias na malha fundida e redistribuindo os resultados de volta para a malha original.

No entanto, uma limitação crítica da formulação de SRD ponderado existente é a falta de um mecanismo de "desligamento". Diferentemente dos métodos clássicos de redistribuição de fluxo, onde o pós-processamento naturalmente desaparece quando os desequilíbrios de fluxo são zero, o SRD ponderado aplica um operador de redistribuição não trivial em cada passo de tempo, independentemente da magnitude da atualização de volume finito. Essa ativação contínua leva a duas deficiências práticas:

1. Falha em preservar estados estacionários: Mesmo quando o equilíbrio de fluxo é exatamente zero, o SRD substitui valores pontuais por médias ponderadas por volume, impedindo a convergência para precisão de máquina e introduzindo erros de $O(1)$ em regiões estáticas.
2. Dissipação desnecessária: Em regiões suaves com atualizações pequenas, mas não nulas, a aplicação contínua da redistribuição introduz difusão numérica proporcional ao gradiente da vizinhança, independentemente da dinâmica física real.

Metodologia
Este artigo introduz a Redistribuição de Estado por Magnitude de Atualização (UM-SRD), uma extensão de desligamento do SRD ponderado projetada para fundir o operador SRD com o operador identidade com base na magnitude local da atualização de volume finito.

O método opera dentro do framework padrão de três estágios de passo de tempo:

1. Atualização de Volume Finito: Calcular o estado pré-redistribuição $U^*$ usando o esquema base.
2. Agregação de Vizinhança: Calcular médias ponderadas de vizinhança $\hat{Q}^*_j$ em vizinhanças fundidas $M_j$.
3. Redistribuição Fundida: Em vez de aplicar o operador SRD $S_j$ diretamente, o UM-SRD aplica uma combinação convexa:
   $$R_j(s_j) = (1 - s_j)I + s_j S_j$$
   onde $I$ é o operador identidade e $s_j \in [0, 1]$ é um parâmetro de mistura.

O parâmetro de mistura $s_j$ é determinado por uma função suave de desligamento, definida localmente $h(\eta_j)$:
$$s_j = \frac{\eta_j^p}{\eta_j^p + \tau^p}$$
Aqui, $\eta_j$ é um indicador normalizado da magnitude máxima de atualização absoluta dentro da vizinhança $M_j$ (especificamente $\eta_j = \frac{\max \|\Delta U_i\|}{\epsilon + \max \|\Delta U_k\|}$).

• Regime de Atualização Zero: Se todas as atualizações na vizinhança estiverem na precisão de máquina ($\Delta U_i \approx 0$), então $\eta_j = 0$, $s_j = 0$, e o método reduz-se exatamente à identidade ($U^{n+1} = U^*$).
• Regime Ativo: Se ocorrerem atualizações significativas (por exemplo, perto de descontinuidades ou quando a rigidez de células pequenas está ativa), $\eta_j \approx 1$, $s_j \approx 1$, e o método comporta-se como SRD ponderado completo.

Contribuições Analíticas Principais
O artigo fornece uma análise rigorosa do UM-SRD para um modelo de advecção linear unidimensional com uma única célula cortada pequena:

• Conservação: O operador fundido é provado ser conservativo por construção para qualquer $s_j \in [0, 1]$.
• Diminuição da Variação Total (TVD): Sob a mesma condição CFL do esquema de upwind de primeira ordem base, o UM-SRD é provado ser TVD na malha de células cortadas original. Isso depende da combinação convexa da identidade e do operador TVD SRD.
• Preservação de Estado Estacionário: É provado que o UM-SRD preserva exatamente os estados estacionários do esquema base. Se a atualização de volume finito for zero, $s_j=0$, e nenhuma redistribuição ocorre, evitando o desvio espúrio inerente ao SRD padrão.
• Precisão: Em regiões suaves, o UM-SRD é mostrado como uma perturbação do esquema base. Se o parâmetro de mistura escala apropriadamente com o tamanho da malha (ou através do indicador normalizado usado na prática), a contribuição adicional ao erro de truncamento local é de ordem superior ao esquema base de primeira ordem. Consequentemente, o UM-SRD retém a taxa de convergência de primeira ordem do método de upwind subjacente.

Resultados Numéricos
Experimentos numéricos em problemas modelo unidimensionais e bidimensionais confirmam as previsões teóricas:

• Convergência: O UM-SRD exibe convergência de primeira ordem nas normas $L_1$ e $L_\infty$ para advecção suave, com erros indistinguíveis do esquema base e do SRD ponderado padrão (diferenças $< 0,01\%$).
• Comportamento de Desligamento: Em testes com pulsos compactos, o parâmetro de mistura $s_j$ cai exatamente para zero quando o pulso está ausente, recuperando o esquema base, e sobe para $\approx 0,99$ quando o pulso interage com a célula pequena.
• Estabilidade: O UM-SRD estabiliza o esquema base sob passos de tempo de malha completa (onde o CFL local na célula pequena é 10 vezes o limite de estabilidade), prevenindo a divergência exponencial observada no esquema de upwind base.
• Propriedades Não-Oscilatórias: Testes com dados descontínuos mostram que o UM-SRD permanece não-oscilatório, correspondendo ao comportamento do SRD padrão perto da célula pequena.
• Preservação de Estado Estacionário: Experimentos confirmam que, enquanto o SRD padrão introduz erros de $O(1)$ em perfis estáticos, o UM-SRD mantém o perfil inicial com precisão de máquina.
• Extensão 2D: O algoritmo foi aplicado com sucesso a uma malha cartesiana 2D com uma fronteira embutida, demonstrando convergência de primeira ordem sem modificação na lógica de mistura.

Significado
O artigo afirma que o UM-SRD oferece uma modificação prática e mínima à infraestrutura existente de SRD ponderado. Ao reutilizar as estruturas de vizinhança e pesos estabelecidas e adicionar apenas um mecanismo de mistura local, ele dota o SRD com uma propriedade de atualização que desaparece. Isso combina os benefícios de estabilidade da redistribuição de estado para células cortadas arbitrariamente pequenas com as propriedades de preservação de precisão dos métodos de redistribuição de fluxo, permitindo especificamente a preservação exata de estados estacionários e reduzindo a dissipação desnecessária em regiões suaves. Os autores observam que, embora a análise atual seja restrita à advecção linear escalar 1D, a estrutura algorítmica é independente da dimensão e pronta para integração em códigos de células cortadas de produção.