A Lattice Boltzmann Method for Non-Newtonian Blood Flow in Coiled Intracranial Aneurysms

Este artigo apresenta um método de Lattice Boltzmann específico para o paciente que modela o fluxo sanguíneo não newtoniano em aneurismas intracranianos coilados, tratando a massa de coils como um meio poroso inhomogêneo, validando assim um fluxo de trabalho para avaliar a hemodinâmica pós-tratamento.

O essencial

A Visão Geral: Consertando um Balão Vazando

Imagine um ponto fraco minúsculo em um balão de água dentro da sua cabeça. Em termos médicos, isso é um aneurisma intracraniano. Se ele estourar, causa um tipo perigoso de derrame.

Para consertar isso, os médicos frequentemente usam uma técnica chamada embolização com coils. Eles inserem um fio minúsculo no balão e o preenchem com uma malha de espirais metálicas. Pense nisso como enfiar uma esponja em um balde vazando. O objetivo é desacelerar a água correndo (o sangue) dentro do balão, para que ela pare de bater nas paredes fracas, dando ao balão a chance de cicatrizar ou se recuperar.

O Problema: Adivinhar o Resultado

A parte complicada é que o "balão vazando" de cada paciente tem um formato diferente, e cada médico enche a "esponja" (os coils) de uma maneira ligeiramente diferente. Antes da cirurgia, é muito difícil prever exatamente como o fluxo de água mudará assim que os coils estiverem dentro. Atualmente, os médicos dependem fortemente de sua experiência passada para adivinhar se o tratamento funcionará.

A Solução: Um "Simulador de Voo" Digital

Os autores deste artigo criaram um programa de computador que atua como um simulador de voo para o fluxo sanguíneo. Em vez de adivinhar, eles podem executar um teste virtual para ver exatamente como o sangue se comportará após a inserção dos coils.

Veja como eles construíram esse simulador:

1. O Truque da "Esponja Porosa"
Geralmente, para simular a água fluindo através de uma pilha bagunçada de fios metálicos (os coils), um computador precisa desenhar cada fio individualmente. Isso é como tentar contar cada grão de areia em um balde; leva uma eternidade e requer um supercomputador.

Os autores encontraram uma maneira mais inteligente. Em vez de desenhar cada fio, eles trataram a pilha de coils como uma esponja.

• A Analogia: Imagine que você está derramando água através de uma floresta densa. Você poderia tentar calcular como a água flui ao redor de cada tronco de árvore individual (muito difícil). Ou, você poderia apenas dizer: "Toda esta área é um pântano espesso e de movimento lento" (muito mais fácil).
• A Ciência: Eles usaram um modelo matemático chamado Navier-Stokes Média de Volume (VANSE). Isso trata a área preenchida por coils como um "meio poroso" (uma esponja), onde a água desacelera com base no quão apertada está a esponja.

2. O Fator "Sangue Inteligente"
O sangue não é como água; é "espesso" e muda como flui dependendo da velocidade (como ketchup ou mel). O modelo de computador leva em conta esse comportamento "não newtoniano", garantindo que a simulação pareça sangue real, e não apenas água comum.

3. O "Motor Rápido" (Método de Lattice Boltzmann)
Para fazer esses cálculos rodarem rápido o suficiente para serem úteis, eles usaram um motor matemático específico chamado Método de Lattice Boltzmann (LBM).

• A Analogia: Pense nisso como um motor de videogame de alta velocidade. Enquanto outros métodos podem tentar resolver a física de todo o oceano de uma vez, o LBM divide o oceano em pequenos blocos gerenciáveis e simula como as partículas saltam entre eles. Isso permite que a simulação rode incrivelmente rápido em placas gráficas modernas (GPUs).

O Que Eles Fizeram no Estudo

A equipe pegou uma tomografia computadorizada real do aneurisma cerebral de um paciente e construiu um modelo 3D dele. Em seguida, executaram dois tipos de simulações:

1. Versão "Super Detalhada": Eles modelaram cada fio individual do coil (como contar cada árvore na floresta).
2. Versão "Esponja": Eles usaram seu novo modelo de "meio poroso" (como tratar a floresta como um pântano).

Eles testaram isso com três níveis diferentes de "apertamento" (densidades de preenchimento de 15%, 20% e 25%), que são as quantidades que os médicos realmente usam na cirurgia.

Os Resultados: A Esponja Funciona!

Os resultados foram encorajadores:

• Precisão: O modelo "Esponja" deu resultados quase idênticos ao modelo "Super Detalhado". O método "Esponja" foi muito mais rápido, mas não perdeu a visão geral.
• Redução do Fluxo: À medida que adicionavam mais coils (tornando a esponja mais densa), o fluxo sanguíneo dentro do aneurisma desacelerou significativamente.
• Segurança: A "tensão de cisalhamento na parede" (a força do sangue esfregando contra a parede fraca) caiu dramaticamente. No aneurisma não tratado, a parede estava sendo batida com força. Com os coils, a força caiu cerca de 40%, sugerindo um risco muito menor de o balão estourar.

A Conclusão

Este artigo apresenta uma maneira nova e mais rápida de simular tratamentos de aneurisma cerebral. Ao tratar os coils metálicos como uma "esponja" em vez de contar cada fio, os médicos poderiam potencialmente executar simulações específicas para cada paciente rapidamente. Esse fluxo de trabalho permite uma melhor avaliação de se um plano de tratamento específico desacelerará com sucesso o fluxo sanguíneo e protegerá o paciente, afastando-se do puro palpite em direção a decisões baseadas em dados.

Resumo técnico

Resumo Técnico: Um Método de Boltzmann em Rede para Fluxo Sanguíneo Não Newtoniano em Aneurismas Intracranianos Enrolados

Declaração do Problema
Os aneurismas intracranianos são uma causa principal de acidente vascular cerebral hemorrágico, frequentemente tratados por meio de embolização com coils para induzir redução do fluxo e mitigação da tensão de cisalhamento na parede (WSS). No entanto, prever os resultados hemodinâmicos de tais tratamentos permanece um desafio em aberto. Uma previsão precisa requer geometria e parâmetros específicos do paciente, particularmente a representação da distribuição heterogênea dos coils dentro do aneurisma. Abordagens tradicionais enfrentam um compromisso: simulações totalmente resolvidas de geometrias discretas de coils são computacionalmente caras, enquanto modelos simplificados frequentemente falham em capturar a heterogeneidade do empacotamento dos coils, o que é crucial para previsões corretas do fluxo. Além disso, o sangue é um fluido não newtoniano, adicionando complexidade reológica à simulação.

Metodologia
Os autores propõem um fluxo de trabalho que combina geometria específica do paciente com uma abordagem de média volumétrica para modelar o aneurisma enrolado como um meio poroso inhomogêneo. A metodologia central envolve:

1. Equações Governantes: O modelo utiliza as Equações de Navier-Stokes de Média Volumétrica (VANSE) para fluxo não newtoniano. O domínio do fluido é caracterizado por um campo de porosidade ($\phi$) representando a fração volumétrica do fluido. A equação de momento inclui um termo de força de corpo ($f(u)$) modelando o arrasto de Darcy e Forchheimer para contabilizar o meio poroso (os coils).
2. Reologia: A viscosidade do sangue é modelada usando o modelo Carreau–Yasuda, que accounts para o comportamento de afinamento por cisalhamento. A viscosidade dinâmica depende da taxa de deformação por cisalhamento.
3. Esquema Numérico: As equações são resolvidas usando um Método de Boltzmann em Rede (LBM) adaptado ao problema em uma rede D3Q27.
   • Tratamento de Meios Porosos: O LBM incorpora porosidade espacialmente variável e uma taxa de relaxamento variável para lidar com a viscosidade não newtoniana e as forças de arrasto poroso.
   • Implementação de Força: O esquema de força de Guo é utilizado para incorporar as forças de arrasto, que dependem quadraticamente da velocidade.
   • Solução de Velocidade: Para permeabilidade isotrópica, a velocidade é resolvida explicitamente. Para permeabilidade anisotrópica, uma iteração de ponto fixo é proposta, embora os autores observem que isso aumenta significativamente o tempo computacional.
4. Geometria e Configuração da Simulação:
   • Geometria: Uma geometria de aneurisma específica do paciente, derivada de uma tomografia computadorizada (TC) anonimizada, foi utilizada.
   • Representação dos Coils: Os coils foram simulados usando um modelo de haste elástica discreta com densidades de empacotamento de 15%, 20% e 25%. As estruturas discretas de coils resultantes foram convertidas em campos de porosidade contínuos por meio de convolução com janela.
   • Condições de Contorno: Condições de não-deslizamento foram aplicadas nas paredes (reflexão de meio caminho). A entrada foi impulsionada por um perfil de velocidade resolvido no tempo, derivado de um modelo de rede arterial 1D-0D, enquanto a saída utilizou um método de extrapolação.
   • Implementação: As simulações foram implementadas usando o framework XLB, aproveitando a aceleração por GPU (Nvidia RTX A6000) para computação paralela.

Principais Contribuições

• Abordagem de Modelagem Híbrida: O artigo integra LBM para VANSE com reologia não newtoniana para simular o fluxo sanguíneo em aneurismas enrolados, tratando a massa de coils como um meio poroso em vez de resolver fios individuais.
• Validação da Média: O estudo fornece uma comparação direta entre simulações de coils totalmente resolvidas e a abordagem de meio poroso de média volumétrica.
• Fluxo de Trabalho Específico do Paciente: Os autores demonstram um pipeline completo, desde imagens médicas (TC) até avaliação hemodinâmica, incluindo a geração de campos de porosidade a partir de implantações de coils simuladas.

Resultados

• Redução do Fluxo: As simulações confirmam que a inserção de coils reduz significativamente a velocidade do fluxo dentro do saco do aneurisma, mantendo o fluxo no vaso pai. A redução é mais pronunciada com maiores densidades de empacotamento (até 25%).
• Validade do Modelo: As simulações de meio poroso (média volumétrica) mostram boa concordância com simulações totalmente resolvidas. Embora a abordagem de média perca efeitos locais de pequena escala, o comportamento geral do fluxo e as quantidades macroscópicas são consistentes.
• Tensão de Cisalhamento na Parede (WSS): O tratamento reduz drasticamente a magnitude da WSS. No aneurisma não tratado, a WSS atingiu aproximadamente 2 Pa, enquanto aneurismas enrolados mostraram magnitudes abaixo de 1,5 Pa, indicando um risco reduzido de ruptura. A diferença na WSS entre as três densidades de empacotamento (15%, 20%, 25%) foi encontrada como pequena, sugerindo que a distribuição da colocação dos coils pode ser mais influente do que a densidade de empacotamento isoladamente.
• Métricas Quantitativas: Em toda a região do aneurisma, o modelo de média volumétrica previu uma redução na magnitude do fluxo de aproximadamente 50% e uma redução na WSS de aproximadamente 40% em comparação com o caso não tratado.

Significado e Alegações
Os autores afirmam que este fluxo de trabalho permite uma avaliação específica do paciente das mudanças no fluxo sanguíneo devido a tratamentos potenciais com coils. O significado principal reside na validade da abordagem de meio poroso de média volumétrica, que oferece uma alternativa computacionalmente eficiente às simulações totalmente resolvidas, sem sacrificar a precisão das métricas globais de fluxo.

O artigo conclui modestamente que, embora levar em conta a permeabilidade anisotrópica pudesse teoricamente aumentar a precisão, a suposição isotrópica atual é provavelmente suficiente para a prática clínica, dadas as incertezas na colocação exata dos coils e a necessidade de tomada de decisão rápida. Os autores identificam trabalhos futuros como necessários para calibração detalhada de parâmetros, quantificação de incertezas e análise de extensões de permeabilidade anisotrópica.