Multiharmonic algorithms for contrast-enhanced ultrasound

Resumo Técnico: Algoritmos Multiharmônicos para Ultrassom com Contraste

1. Problema e Contexto

O ultrassom com contraste (CEUS) utiliza microbolhas de gás para melhorar a resolução de imagens e terapias direcionadas (como entrega de fármacos). A interação entre as ondas acústicas e as microbolhas é altamente não linear, gerando harmônicos superiores que são cruciais para a imagem e terapia.

• Desafio Computacional: A modelagem tradicional no domínio do tempo (resolvendo equações diferenciais parciais acopladas a equações diferenciais ordinárias) exige passos de tempo extremamente pequenos devido às diferentes escalas de tempo da propagação da onda e da oscilação das microbolhas. Isso torna as simulações computacionalmente proibitivas para aplicações em tempo real.
• Objetivo: Desenvolver algoritmos eficientes baseados em uma discretização temporal que utilize expansões multiharmônicas (séries de Fourier truncadas) para capturar esses efeitos não lineares sem a necessidade de resolver o problema em cada passo de tempo fino.

2. Modelo Matemático

O sistema estudado acopla a Equação de Westervelt (para a pressão acústica $p$) com uma aproximação baseada em volume da Equação de Rayleigh-Plesset (para a variação de volume $v$ das microbolhas).

• Equação de Pressão (Westervelt): Inclui termos de amortecimento forte ($-b\Delta p_t$) e não linearidade quadrática ($\eta(p^2)_{tt}$).
• Equação de Volume (ODE): Descreve a oscilação das microbolhas acionada pela pressão, incluindo termos não lineares quadráticos e cúbicos.
• Condições de Contorno: Condições de absorção de tipo misto ($\beta p_t + \gamma p + \nabla p \cdot n = 0$) e soluções periódicas no tempo (regime de cavitação estável).

3. Metodologia

Os autores propõem uma abordagem baseada em Ansatz Multiharmônico, onde as soluções de pressão e volume são aproximadas por séries de Fourier truncadas com $N$ harmônicos:
$$p(x,t) \approx \sum_{m=0}^{N} \text{Re}\left( e^{i m \omega t} p_m(x) \right)$$
O trabalho explora duas formulações:

1. Campo Real: Mantém explicitamente a parte real da série, garantindo que a solução seja real. Isso gera um sistema acoplado complexo com termos de soma e diferença de frequências.
2. Campo Complexo: Abandona a restrição de realidade explícita (trabalhando apenas com o espectro de frequências positivas), resultando em um sistema triangular inferior que é computacionalmente mais eficiente, embora formalmente investigado.

Algoritmos Desenvolvidos:

• Algoritmo de Corte Multiharmônico: Derivação direta do sistema acoplado para os coeficientes de Fourier.
• Algoritmo Linearizado Multiharmônico: Uma abordagem iterativa onde os termos não lineares são tratados usando os harmônicos da iteração anterior. Isso transforma o problema não linear em uma sequência de problemas de Helmholtz lineares e equações algébricas, resolvíveis por substituição sucessiva.

4. Contribuições Principais

1. Existência de Soluções: Prova rigorosa da existência e unicidade de soluções periódicas no tempo para o sistema acoplado Westervelt-ODE, sob condições de pequenas perturbações (densidade de microbolhas e fonte acústica limitadas).
2. Derivação de Algoritmos: Formulação detalhada dos sistemas de equações para os coeficientes de Fourier ($p_m, v_m$) tanto para campos reais quanto complexos.
3. Análise de Erro e Convergência:
   • Estabelecimento de limites de erro para o algoritmo linearizado em termos do número de harmônicos ($N$) e do erro da iteração de ponto fixo.
   • Prova de que o esquema converge para a solução do sistema original quando $N \to \infty$.
4. Análise Numérica: Implementação via Método de Elementos Finitos (FEM) e comparação entre as formulações real e complexa, além de estudos sobre a influência do número de harmônicos e da presença de microbolhas.

5. Resultados Numéricos

• Eficiência Computacional: A formulação complexa (que resulta em um sistema triangular) mostrou-se numericamente equivalente à formulação real para os parâmetros testados, mas com menor custo computacional.
• Convergência Harmônica: Em cenários não lineares fortes, a inclusão de 3 a 5 harmônicos foi suficiente para capturar com precisão a dinâmica da onda, com erros residuais caindo abaixo da precisão da máquina para $N \ge 10$.
• Influência das Microbolhas:
  • As microbolhas introduzem um efeito de atenuação significativa na propagação da onda.
  • Elas alteram a velocidade de onda efetiva e aumentam o parâmetro de não linearidade efetivo.
  • A presença de microbolhas suprime rapidamente a amplitude dos harmônicos superiores em comparação com um meio sem bolhas, onde a decaimento é mais lento.
• Discretização Espacial: A análise mostrou que a discretização espacial (malha FEM) domina o erro em relação ao corte harmônico, e que o método é viável mesmo em laptops padrão sem GPUs dedicadas.

6. Significado e Impacto

Este trabalho preenche uma lacuna importante na modelagem de ultrassom com contraste ao fornecer:

• Fundamentação Teórica: A primeira prova rigorosa de bem-postura e convergência para algoritmos multiharmônicos aplicados a sistemas acoplados de ondas não lineares e dinâmica de microbolhas.
• Viabilidade Prática: Demonstra que algoritmos multiharmônicos são uma alternativa viável e eficiente aos métodos de domínio do tempo, permitindo simulações mais rápidas para aplicações clínicas em tempo real e planejamento terapêutico.
• Insights Físicos: Quantifica como a densidade de microbolhas e a amplitude da fonte afetam a geração de harmônicos e a atenuação, fornecendo diretrizes para a escolha do número de harmônicos necessários em diferentes cenários de imagem.

O artigo sugere futuras extensões para incluir sub-harmônicos (deformações não esféricas) e aplicações em problemas inversos (reconstrução de parâmetros a partir de sinais medidos).