A stabilized finite element method for a flow problem arising from 4D flow magnetic resonance imaging

Resumo Técnico: Um Método de Elementos Finitos Estabilizado para um Problema de Fluxo Decorrente de 4D Flow MRI

Enunciado do Problema
O artigo aborda o desafio de avaliar a qualidade dos dados de Ressonância Magnética (RM) de fluxo 4D (4D Flow MRI) e reconstruir biomarcadores hemodinâmicos, especificamente a pressão, a partir de medições de velocidade ruidosas. Embora o 4D Flow MRI seja o padrão-ouro não invasivo para quantificar o fluxo sanguíneo cardiovascular, técnicas agressivas de subamostragem utilizadas para reduzir os tempos de aquisição introduzem artefatos e ruído. Além disso, algoritmos de reconstrução padrão frequentemente ignoram as leis físicas que regem a dinâmica de fluidos (especificamente a conservação de massa e as equações de Navier-Stokes), levando a erros sistemáticos em quantidades derivadas, como a pressão.

Os autores modelam o campo de velocidade físico real $\mathbf{u}$ como a soma da velocidade reconstruída por RM $\mathbf{u}_m$ (tratada como um dado) e um componente de ruído desconhecido $\mathbf{w}$. Ao substituir essa decomposição nas equações de Navier-Stokes incompressíveis, o problema é reformulado para resolver a velocidade do ruído $\mathbf{w}$ e a pressão do fluido $p$. Esta abordagem transforma o problema de reconstrução de pressão em um problema de valor de contorno onde os termos de fonte dependem dos dados ruidosos de MRI $\mathbf{u}_m$.

Metodologia
Os autores propõem um método de elementos finitos (FEM) estabilizado, projetado para lidar com as características específicas deste problema reformulado:

1. Derivação do Modelo: Partindo das equações de Navier-Stokes, os autores derivam um problema do tipo Oseen modificado (versão linearizada) para a velocidade do ruído $\mathbf{w}$ e a pressão $p$. A formulação leva em conta o divergente dos dados de MRI $\mathbf{u}_m$, que é geralmente não nulo devido ao ruído.
2. Interpolação de Ordem Igual: Para acomodar a necessidade prática de eficiência e a natureza do processamento de dados de MRI, o método emprega espaços de elementos finitos de ordem igual tanto para velocidade quanto para pressão (por exemplo, $P_k$ para ambos). Esta escolha viola a condição padrão de estabilidade inf-sup (Ladyzhenskaya-Babuška-Brezzi).
3. Estratégia de Estabilização: Para restaurar a estabilidade e permitir o uso de pares de ordem igual, os autores introduzem uma formulação estabilizada inspirada em "Métodos de Elementos Finitos Estabilizados Incomuns" (referenciando especificamente o trabalho em [5]). A estabilização inclui:
   • Um termo de mínimos quadrados sobre o resíduo da equação de momento.
   • Um termo que lida com a derivada convectiva.
   • Um termo de estabilização de pressão.
   • Um parâmetro $\tau_T$ dependente do tamanho do elemento $h_T$, viscosidade $\mu$ e o parâmetro do passo de tempo $\sigma$.
4. Linearização: A análise teórica é conduzida em uma versão linearizada do problema (equação de Oseen modificada). O caso não linear (Navier-Stokes perturbado) é abordado numericamente usando iterações de Picard.

Principais Contribuições e Resultados Teóricos

• Bem-posedness (Bem-posto): Os autores provam que o problema de Oseen modificado linearizado admite uma solução única sob condições específicas, notadamente uma restrição no parâmetro $\sigma$ em relação ao gradiente dos dados de MRI ($\sigma > 4\rho \|\nabla \mathbf{u}_m\|_\infty$).
• Estabilidade e Estimativas de Erro: Para o esquema de elementos finitos estabilizado proposto, os autores estabelecem:
  • Coercividade: A forma bilinear satisfaz uma estimativa de estabilidade em uma norma de energia específica.
  • Convergência Ótima: Eles derivam estimativas de erro de ordem ótima tanto para a velocidade quanto para a pressão. Os limites de erro dependem do tamanho da malha $h$, da regularidade da solução e da qualidade de aproximação do campo de velocidade convectivo.
  • Erro de Consistência: O artigo reconhece que o esquema não é estritamente consistente devido à aproximação do termo convectivo no contexto não linear, mas demonstra que esse erro de consistência pode ser limitado de forma ótima.
• Reconstrução de Pressão: O método fornece um procedimento alternativo e não invasivo para a reconstrução de pressão ao resolver o problema inverso de estimar simultaneamente o ruído e a pressão.

Resultados Numéricos
Os autores validam o método através de três experimentos numéricos utilizando a biblioteca FEniCS:

1. Fluxo Kovasznay (Teste Acadêmico): Usando uma solução analítica, o método demonstra taxas de convergência ótimas para elementos lineares ($k=1$), quadráticos ($k=2$) e cúbicos ($k=3$) através de vários valores de viscosidade. Os resultados sustentam-se tanto para o modelo linearizado quanto para o problema de Navier-Stokes perturbado não linear.
2. Velocidade Aleatória em um Domínio Curvo: Simulando uma configuração de 4D Flow MRI com dados de velocidade constantes por partes em uma malha quadrilátera, o método reconstrói com sucesso o campo de pressão. Embora algum desfoque seja observado (atribuído a erros de interpolação entre os dados constantes por partes e os elementos finitos contínuos), a recuperação da pressão permanece robusta.
3. Velocidade de Navier-Stokes Interpolada: Em um cenário realista, os dados de velocidade são gerados a partir de uma simulação de Navier-Stokes de baixa resolução, interpolados e, então, utilizados para reconstruir a pressão. A pressão reconstruída e os campos de velocidade corrigidos mostram alta concordância com uma solução de referência de alta resolução, com erros máximos em torno de 1,6% em seções transversais específicas.

Significância e Alegações
O artigo afirma que o método de elementos finitos estabilizado proposto oferece um arcabouço matemático robusto para:

• Avaliação de Qualidade: Estimar a qualidade dos campos de velocidade medidos ao quantificar o componente de "ruído" $\mathbf{w}$.
• Reconstrução de Pressão: Fornecer uma alternativa baseada na física para métodos de reconstrução de pressão puramente fenomenológicos, evitando assim procedimentos invasivos.
• Eficiência: Permitizar o uso de pares de elementos finitos de ordem igual computacionalmente eficientes sem sacrificar a estabilidade, o que é crucial para problemas de fluxo 3D complexos.

Os autores observam que, embora o caso linearizado seja totalmente analisado, a extensão para problemas totalmente dependentes do tempo e a busca por discretizações livres de divergência ou robustas à pressão permanecem como questões abertas para pesquisas futuras. O trabalho posiciona-se como um passo em direção ao cálculo mais confiável de biomarcadores hemodinâmicos a partir de dados clínicos de 4D Flow MRI.