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

Este artículo presenta un método de red de Boltzmann específico para el paciente que modela el flujo sanguíneo no newtoniano en aneurismas intracraneales embolizados con coils, tratando la masa de coils como un medio poroso inhomogéneo, validando así un flujo de trabajo para evaluar la hemodinámica post-tratamiento.

Lo esencial

La Gran Imagen: Reparar un Globo con Fugas

Imagina un punto débil y minúsculo en un globo de agua dentro de tu cabeza. En términos médicos, esto es un aneurisma intracraneal. Si estalla, provoca un tipo de accidente cerebrovascular peligroso.

Para repararlo, los médicos suelen utilizar una técnica llamada embolización con coils. Insertan un alambre diminuto en el globo y lo rellenan con una malla de espirales metálicas. Piensa en esto como meter una esponja en un cubo con fugas. El objetivo es frenar el agua que corre rápidamente (la sangre) dentro del globo para que deje de golpear las paredes débiles, dando al globo la oportunidad de sanar o cicatrizar.

El Problema: Adivinar el Resultado

La parte complicada es que el "globo con fugas" de cada paciente tiene una forma diferente, y cada médico rellena la "esponja" (los coils) de una manera ligeramente distinta. Antes de la cirugía, es muy difícil predecir exactamente cómo cambiará el flujo del agua una vez que los coils estén dentro. Actualmente, los médicos dependen en gran medida de su experiencia pasada para adivinar si el tratamiento funcionará.

La Solución: Un "Simulador de Vuelo" Digital

Los autores de este artículo crearon un programa informático que actúa como un simulador de vuelo para el flujo sanguíneo. En lugar de adivinar, pueden realizar una prueba virtual para ver exactamente cómo se comportará la sangre después de insertar los coils.

Así es como construyeron este simulador:

1. El Truco de la "Esponja Porosa"
Por lo general, para simular el agua fluyendo a través de un montón desordenado de alambres metálicos (los coils), una computadora tiene que dibujar cada alambre individual. Esto es como intentar contar cada grano de arena en un cubo; lleva una eternidad y requiere una supercomputadora.

Los autores encontraron una forma más inteligente. En lugar de dibujar cada alambre, trataron el montón de coils como una esponja.

• La Analogía: Imagina que estás vertiendo agua a través de un bosque denso. Podrías intentar calcular cómo fluye el agua alrededor de cada tronco de árbol individual (muy difícil). O, podrías simplemente decir: "Toda esta área es un pantano espeso y de movimiento lento" (mucho más fácil).
• La Ciencia: Utilizaron un modelo matemático llamado Navier-Stokes Promedio por Volumen (VANSE). Esto trata el área llena de coils como un "medio poroso" (una esponja) donde el agua se frena según lo apretada que esté la esponja.

2. El Factor de la "Sangre Inteligente"
La sangre no es como el agua; es "espesa" y cambia su forma de fluir dependiendo de la velocidad a la que se mueve (como el ketchup o la miel). El modelo informático tiene en cuenta este comportamiento "no newtoniano", asegurando que la simulación se sienta como sangre real, no solo como agua simple.

3. El "Motor Rápido" (Método de Boltzmann en Red)
Para que estos cálculos se ejecuten lo suficientemente rápido como para ser útiles, utilizaron un motor matemático específico llamado Método de Boltzmann en Red (LBM).

• La Analogía: Piensa en esto como un motor de videojuegos de alta velocidad. Mientras que otros métodos podrían intentar resolver la física de todo el océano de una vez, LBM divide el océano en pequeñas baldosas manejables y simula cómo las partículas rebotan entre ellas. Esto permite que la simulación se ejecute increíblemente rápido en las tarjetas gráficas modernas (GPU).

Lo Que Hicieron en el Estudio

El equipo tomó una tomografía computarizada real del aneurisma cerebral de un paciente y construyó un modelo 3D del mismo. Luego, ejecutaron dos tipos de simulaciones:

1. La Versión "Super Detallada": Modelaron cada alambre individual del coil (como contar cada árbol en el bosque).
2. La Versión "Esponja": Utilizaron su nuevo modelo de "medio poroso" (como tratar el bosque como un pantano).

Probaron esto con tres niveles diferentes de "atascamiento" (densidades de empaquetado del 15%, 20% y 25%), que son las cantidades que los médicos realmente utilizan en la cirugía.

Los Resultados: ¡La Esponja Funciona!

Los resultados fueron alentadores:

• Precisión: El modelo "Esponja" dio resultados casi idénticos al modelo "Super Detallado". El método "Esponja" fue mucho más rápido, pero no perdió la gran imagen.
• Reducción del Flujo: A medida que añadieron más coils (haciendo la esponja más densa), el flujo sanguíneo dentro del aneurisma disminuyó significativamente.
• Seguridad: El "esfuerzo cortante de la pared" (la fuerza de la sangre frotándose contra la pared débil) disminuyó drásticamente. En el aneurisma sin tratar, la pared estaba siendo golpeada con fuerza. Con los coils, la fuerza disminuyó aproximadamente un 40%, lo que sugiere un riesgo mucho menor de que el globo estalle.

La Conclusión

Este artículo presenta una nueva y más rápida forma de simular tratamientos de aneurismas cerebrales. Al tratar los coils metálicos como una "esponja" en lugar de contar cada alambre, los médicos podrían potencialmente ejecutar simulaciones específicas para cada paciente de manera rápida. Este flujo de trabajo permite una mejor evaluación de si un plan de tratamiento específico frenará con éxito el flujo sanguíneo y protegerá al paciente, alejándose de la pura adivinanza hacia decisiones basadas en datos.

Resumen técnico

Resumen Técnico: Un Método de Red de Boltzmann para el Flujo Sanguíneo No Newtoniano en Aneurismas Intracraneales Enrollados

Planteamiento del Problema
Los aneurismas intracraneales son una causa principal de accidente cerebrovascular hemorrágico, a menudo tratados mediante embolización con coils para inducir una reducción del flujo y una mitigación del esfuerzo cortante en la pared (WSS). Sin embargo, predecir los resultados hemodinámicos de tales tratamientos sigue siendo un desafío abierto. Una predicción precisa requiere geometría y parámetros específicos del paciente, particularmente la representación de la distribución inhomogénea de los coils dentro del aneurisma. Los enfoques tradicionales enfrentan una compensación: las simulaciones completamente resueltas de geometrías discretas de coils son computacionalmente costosas, mientras que los modelos simplificados a menudo fallan en capturar la heterogeneidad del empaquetamiento de los coils, lo cual es crucial para predicciones de flujo correctas. Además, la sangre es un fluido no newtoniano, lo que añade complejidad reológica a la simulación.

Metodología
Los autores proponen un flujo de trabajo que combina la geometría específica del paciente con un enfoque promediado en volumen para modelar el aneurisma enrollado como un medio poroso inhomogéneo. La metodología central involucra:

1. Ecuaciones Gobernantes: El modelo utiliza las Ecuaciones de Navier-Stokes Promediadas en Volumen (VANSE) para flujo no newtoniano. El dominio fluido se caracteriza por un campo de porosidad ($\phi$) que representa la fracción volumétrica del fluido. La ecuación de momento incluye un término de fuerza corporal ($f(u)$) que modela la resistencia de Darcy y Forchheimer para dar cuenta del medio poroso (los coils).
2. Reología: La viscosidad de la sangre se modela utilizando el modelo Carreau–Yasuda, que da cuenta del comportamiento de adelgazamiento por cizallamiento. La viscosidad dinámica depende de la tasa de deformación por cizallamiento.
3. Esquema Numérico: Las ecuaciones se resuelven utilizando un Método de Red de Boltzmann (LBM) adaptado al problema en una red D3Q27.
   • Manejo de Medios Porosos: El LBM incorpora porosidad espacialmente variable y una tasa de relajación variable para manejar la viscosidad no newtoniana y las fuerzas de arrastre poroso.
   • Implementación de Fuerzas: Se utiliza el esquema de forzamiento de Guo para incorporar las fuerzas de arrastre, que dependen cuadráticamente de la velocidad.
   • Solución de Velocidad: Para permeabilidad isotrópica, la velocidad se resuelve explícitamente. Para permeabilidad anisotrópica, se propone una iteración de punto fijo, aunque los autores notan que esto aumenta significativamente el tiempo computacional.
4. Geometría y Configuración de la Simulación:
   • Geometría: Se utilizó una geometría de aneurisma específica del paciente derivada de una tomografía computarizada (TC) anonimizada.
   • Representación de Coils: Los coils se simularon utilizando un modelo de varilla elástica discreta con densidades de empaquetamiento de 15%, 20% y 25%. Las estructuras discretas de coils resultantes se convirtieron en campos de porosidad continuos mediante convolución con ventana.
   • Condiciones de Contorno: Se aplicaron condiciones de no deslizamiento en las paredes (rebote a mitad de camino). La entrada fue impulsada por un perfil de velocidad resuelto en el tiempo derivado de un modelo de red arterial 1D-0D, mientras que la salida utilizó un método de extrapolación.
   • Implementación: Las simulaciones se implementaron utilizando el marco XLB, aprovechando la aceleración por GPU (Nvidia RTX A6000) para el cálculo paralelo.

Contribuciones Clave

• Enfoque de Modelado Híbrido: El artículo integra LBM para VANSE con reología no newtoniana para simular el flujo sanguíneo en aneurismas enrollados, tratando la masa de coils como un medio poroso en lugar de resolver cables individuales.
• Validación del Promedio: El estudio proporciona una comparación directa entre simulaciones de coils completamente resueltas y el enfoque de medio poroso promediado en volumen.
• Flujo de Trabajo Específico del Paciente: Los autores demuestran un pipeline completo desde la imagen médica (escáner TC) hasta la evaluación hemodinámica, incluyendo la generación de campos de porosidad a partir de despliegues de coils simulados.

Resultados

• Reducción del Flujo: Las simulaciones confirman que la inserción de coils reduce significativamente la velocidad del flujo dentro del saco del aneurisma mientras mantiene el flujo en el vaso padre. La reducción es más pronunciada con mayores densidades de empaquetamiento (hasta 25%).
• Validez del Modelo: Las simulaciones de medio poroso (promediadas en volumen) muestran un buen acuerdo con las simulaciones completamente resueltas. Si bien el enfoque de promedio pierde efectos locales a pequeña escala, el comportamiento general del flujo y las cantidades macroscópicas son consistentes.
• Esfuerzo Cortante en la Pared (WSS): El tratamiento reduce drásticamente la magnitud del WSS. En el aneurisma sin tratar, el WSS alcanzó aproximadamente 2 Pa, mientras que los aneurismas enrollados mostraron magnitudes por debajo de 1.5 Pa, lo que indica un riesgo reducido de ruptura. La diferencia en el WSS entre las tres densidades de empaquetamiento (15%, 20%, 25%) se encontró que era leve, sugiriendo que la distribución de la colocación de los coils puede ser más influyente que la densidad de empaquetamiento por sí sola.
• Métricas Cuantitativas: En toda la región del aneurisma, el modelo promediado en volumen predijo una reducción de la magnitud del flujo de aproximadamente el 50% y una reducción del WSS de aproximadamente el 40% en comparación con el caso sin tratar.

Significado y Afirmaciones
Los autores afirman que este flujo de trabajo permite una evaluación específica del paciente de los cambios en el flujo sanguíneo debido a posibles tratamientos con coils. El significado principal radica en la validez del enfoque de medio poroso promediado en volumen, que ofrece una alternativa computacionalmente eficiente a las simulaciones completamente resueltas sin sacrificar la precisión de las métricas de flujo globales.

El artículo concluye modestamente que, si bien tener en cuenta la permeabilidad anisotrópica podría teóricamente aumentar la precisión, la suposición isotrópica actual es probablemente suficiente para la práctica clínica, dadas las incertidumbres en la colocación exacta de los coils y la necesidad de una toma de decisiones rápida. Los autores identifican trabajos futuros como necesarios para la calibración detallada de parámetros, la cuantificación de incertidumbres y el análisis de extensiones de permeabilidad anisotrópica.