Body-fitted tracking of 2d open curves with a level set based mesh evolution method

Resumo Técnico: Rastreamento de Curvas Abertas 2D Ajustadas ao Corpo com um Método de Evolução de Malha Baseado em Level Set

Definição do Problema
O artigo aborda o desafio numérico de rastrear o movimento de curvas abertas (ou coleções de curvas abertas disjuntas) em um espaço bidimensional. Embora exista uma literatura significativa para o rastreamento de domínios ou fronteiras fechadas, o tratamento de objetos abertos — curvas com extremidades estritamente contidas dentro do domínio computacional — permanece difícil. Métodos puramente Lagrangianos (baseados em partículas) frequentemente lutam com robustez sob grandes deformações e mudanças topológicas, exigdo reamostragem frequente. Por outro lado, métodos puramente Eulerianos (malha fixa, por exemplo, o Level Set padrão) carecem de uma discretização explícita ajustada ao corpo, o que torna difíceis as computações geométricas e mecânicas precisas (como a avaliação de curvatura ou a resolução de problemas de valor de contorno na curva) sem aproximações. O objetivo é conciliar esses requisitos: manter uma representação malhada precisa para cálculos, permitindo simultaneamente uma evolução robusta e arbitrária da forma da curva.

Metodologia
A estratégia proposta emprega uma abordagem híbrida que acopla duas representações complementares da curva $\Gamma$ em cada iteração do processo de evolução:

1. Representação de Malha Explícita: A curva é discretizada como uma coleção de arestas pertencentes a uma malha triangular conforme do domínio computacional total. Isso permite computações geométricas e mecânicas de alta precisão (por exemplo, Método de Elementos Finitos) diretamente na curva.
2. Representação de Level Set Implícita: A curva é capturada simultaneamente usando duas funções escalares auxiliares, $\phi$ e $\psi$, definidas em todo o domínio.
   • $\phi$ define uma curva fechada $\tilde{\Gamma}$ (o zero level set) que estende a curva aberta $\Gamma$.
   • $\psi$ delimita o segmento aberto $\Gamma$ dentro de $\tilde{\Gamma}$ (especificamente, $\Gamma = \{x \mid \phi(x)=0, \psi(x)<0\}$).
   • As extremidades $\Sigma$ são definidas onde $\phi(x)=0$ e $\psi(x)=0$.
   • Uma condição de ortogonalidade ($\nabla\phi \cdot \nabla\psi = 0$) é imposta próximo às extremidades para evitar colisões indesejadas dos level sets.

Fluxo de Trabalho Algorítmico
A evolução procede via um loop iterativo (Algoritmo 1) que alterna entre essas representações:

1. Malha para Level Set: A partir da curva malhada atual, dois level sets são gerados. Isso envolve o cálculo de distâncias não assinadas e distâncias assinadas usando o Método de Marcha Rápida (Fast Marching Method) para construir $\phi$ (estendendo a curva para um loop fechado) e $\psi$ (definindo o segmento aberto).
2. Cálculo de Velocidade: O campo de velocidade $V$ que impulsiona o movimento é calculado. Esta etapa utiliza a malha explícita para quantidades geométricas precisas (ex: curvatura) ou para resolver problemas de valor de contorno (ex: dinâmica de fluidos ou elasticidade) via Método de Elementos Finitos.
3. Evolução do Level Set: As equações de evolução para $\phi$ e $\psi$ (equações do tipo advecção) são resolvidas sobre um passo de tempo $\Delta t$ usando o Método das Características (integração Runge-Kutta). Isso atualiza a curva implicitamente, lidando com grandes deformações e mudanças topológicas de forma robusta.
4. Level Set para Malha: As novas funções de level set são convertidas de volta em uma malha triangular de alta qualidade, ajustada ao corpo. Isso envolve uma etapa de extração de isosuperfície (marching triangles) para criar uma malha inicial de baixa qualidade, seguida por operadores de remalhagem local (divisão de aresta, colapso, troca, realocação de vértice) para restaurar a qualidade dos elementos enquanto preserva a discretização explícita da curva.

Principais Contribuições

• Acoplamento Inovador: O artigo introduz um framework algorítmico específico que alterna perfeitamente entre uma malha ajustada ao corpo (para precisão) e uma representação de dois level sets (para robustez), especificamente para curvas abertas.
• Geração de Dois Level Sets: Detalha um método para gerar o par de funções de level set ($\phi, \psi$) a partir de uma malha de linha, garantindo a condição de ortogonalidade nas extremidades e a extensão correta da curva aberta em um loop fechado.
• Estratégia de Remalhagem: Descreve um procedimento de remalhagem orientado à qualidade que extrai a curva dos level sets e regenera uma malha de alta qualidade, um passo crítico para manter a estabilidade numérica em simulações físicas.
• Framework de Otimização de Forma: O método é integrado em um contexto de otimização de forma usando o método de variação de contorno de Hadamard, derivando derivadas de forma para funcionais envolvendo tanto o volume da curva quanto suas extremidades.

Resultados e Aplicações
A metodologia é validada através de diversas simulações físicas e problemas de otimização:

• Advecção Analítica: Um teste de benchmark envolvendo uma curva sofrendo grande deformação e retornando ao seu estado inicial demonstra alta precisão (erro de distância de Hausdorff $\approx 10^{-3}$) apesar do estiramento extremo.
• Roll-up de Folha de Vórtice: O método simula a dinâmica de uma folha de vórtice 2D, onde a velocidade é determinada por uma integral fracamente singular sobre a curva. Os resultados mostram boa concordância com benchmarks Lagrangianos existentes.
• Otimização de Perímetro Anisotrópico: O framework otimiza curvas abertas para minimizar o comprimento anisotrópico (métrica de Riemann) na presença de campos de vento e restrições topográficas, encontrando com sucesso caminhos ótimos que desviam contra o vento.
• Manufatura Aditiva: O método otimiza a trajetória do laser para Fusão por Feixe de Elétrons. O objetivo é fundir uma região alvo específica em um leito de pó enquanto evita fundir fora dela, sujeito a equações de difusão de calor. A otimização molda com sucesso a trajetória do laser para corresponder à geometria alvo.
• Reconstrução de Fratura: Um problema inverso é resolvido para reconstruir uma linha de falha subterrânea a partir de medições de deslocamento de superfície. O método minimiza a diferença de mínimos quadrados entre os deslocamentos observados e previstos, recuperando com sucesso a inclinação de um padrão de falha conhecido.

Significância e Alegações
O artigo afirma oferecer uma metodologia robusta e ajustada ao corpo que supera as limitações de abordagens puramente Lagrangianas ou Eulerianas para curvas abertas. Ao combinar a precisão geométrica de representações de malha com a flexibilidade topológica dos métodos de Level Set, o trabalho permite computações mecânicas e geométricas precisas mesmo durante evoluções de grande escala.

Os autores são modestos quanto ao escopo, observando que:

• A implementação atual é 2D, embora a metodologia seja projetada com extensões para 3D em mente (que é um trabalho em andamento).
• Mudanças topológicas (ramificações, junções) são tratadas naturalmente pelo aspecto do Level Set, mas não são o foco principal dos exemplos apresentados.
• A implementação 2D pode ser menos eficiente do que solvers planares especializados, mas serve como uma prova de conceito para uma tecnologia destinada a aplicações 3D (ex: otimização de cascas, telas eletromagnéticas).

O trabalho é apresentado como um passo preliminar para um framework geral de superfícies 3D abertas, com os resultados 2D atuais servindo para validar os componentes algorítmicos centrais.