An Integrated Open Source Software System for the Generation and Analysis of Subject-Specific Blood Flow Simulation Ensembles

Resumen Técnico: Un Sistema Integrado de Software de Código Abierto para la Generación y Análisis de Conjuntos de Simulaciones de Flujo Sanguíneo Específicos del Sujeto

Planteamiento del Problema
El análisis hemodinámico es crítico para el diagnóstico de enfermedades cardiovasculares (por ejemplo, aneurismas, estenosis) y la investigación de parámetros como la turbulencia y el esfuerzo cortante en la pared. Si bien el análisis específico del sujeto requiere combinar imágenes médicas no invasivas (MRI) con Dinámica de Fluidos Computacional (CFD) para resolver las ecuaciones de Navier-Stokes, los flujos de trabajo actuales están fragmentados. Los investigadores a menudo dependen de múltiples herramientas separadas para la segmentación de imágenes, el modelado de superficies 3D, la configuración de simulaciones y el postprocesamiento. Las soluciones existentes, ya sean comerciales (por ejemplo, Ansys, SimScale) o de código abierto (por ejemplo, SimVascular, CRIMSON, OpenFOAM), generalmente carecen de soporte integrado para la generación y el análisis de conjuntos de simulaciones —conjuntos de simulaciones con parámetros variables—. En consecuencia, los investigadores recurren frecuentemente a scripts personalizados, lo que limita la reproducibilidad, la generalización y la accesibilidad para quienes no son expertos en informática o análisis numérico.

Metodología
Los autores abordan estas brechas integrando el marco de simulación Lattice Boltzmann de alto rendimiento OpenLB en el marco de análisis visual interactivo Voreen. Esta integración crea un sistema unificado y de código abierto donde todos los pasos del flujo de trabajo —desde el preprocesamiento de imágenes médicas hasta la ejecución de simulaciones y el análisis de conjuntos— ocurren dentro de un único entorno de software.

• Arquitectura Central: El sistema utiliza OpenLB como un módulo backend flexible para simulaciones de multifísica que emplean el Método Lattice Boltzmann (LBM) para aproximar las ecuaciones de Navier-Stokes incompresibles. OpenLB maneja la resolución numérica en arquitecturas de CPU/GPU. Voreen actúa como frontend, proporcionando una interfaz gráfica de usuario (GUI) personalizable para la gestión del flujo de trabajo, la renderización volumétrica y el análisis interactivo.
• Integración del Flujo de Trabajo: El sistema lee los parámetros de simulación desde archivos XML de parámetros de Voreen. Soporta la importación de archivos geométricos (por ejemplo, .stl) y datos médicos volumétricos (por ejemplo, 4D-PC-MRI).
• Generación de Conjuntos: La herramienta facilita la creación de conjuntos de simulaciones permitiendo a los usuarios configurar parámetros variables (por ejemplo, condiciones de contorno, velocidades de entrada) y ejecutar múltiples miembros de simulación en paralelo.
• Técnicas de Análisis: Para el análisis de conjuntos, el sistema emplea incrustaciones 2D del espacio de similitud. Utiliza una medida de similitud de campo vectorial (propuesta por Leistikow et al.) para comparar los miembros de la simulación entre sí o contra datos medidos. Estas comparaciones se proyectan utilizando Escalamiento Multidimensional (MDS) para visualizar la disimilitud a lo largo del tiempo.

Contribuciones Clave
El artículo presenta tres contribuciones principales:

1. Un Sistema Integrado: Un sistema novedoso, extensible y basado en una GUI amigable que combina OpenLB y Voreen para configurar, realizar y analizar simulaciones 3D de Lattice Boltzmann sin requerir software externo ni transferencias manuales de datos.
2. Casos de Uso Demostrados: La efectividad del sistema se valida mediante tres flujos de trabajo específicos orientados a la investigación desarrollados en colaboración con expertos en el dominio:
   • Caso de Uso A: Configuración de simulaciones basadas en archivos de geometría preexistentes (STL).
   • Caso de Uso B: Generación de simulaciones específicas del sujeto a partir de datos de imágenes médicas (MRI), incluyendo segmentación, creación de máscaras y configuración de simulaciones utilizando la velocidad medida como condiciones de contorno de entrada.
   • Caso de Uso C: Análisis visual de los conjuntos de simulación resultantes para derivar conclusiones, específicamente comparando miembros de simulación contra mediciones utilizando gráficos de similitud y visualización de diferencias voxel a voxel.
3. Evaluación por Expertos: El software fue evaluado mediante sesiones de taller con investigadores de los dominios de CFD e imágenes cardiovasculares. La retroalimentación se utilizó para identificar fortalezas de usabilidad (por ejemplo, configuración rápida de condiciones de contorno) y áreas de mejora (por ejemplo, necesidad de una colocación de límites más flexible y flujos de trabajo especializados de imágenes médicas).

Resultados
Los autores aplicaron la herramienta a tres escenarios del mundo real para demostrar sus capacidades:

• Simulación Basada en Geometría: Configuración y visualización exitosa de resultados de simulación in-situ para una geometría STL dada, permitiendo actualizaciones en tiempo real de la magnitud de la velocidad y los campos de presión.
• Simulación Específica del Sujeto: Procesamiento de datos 4D-PC-MRI de la arteria pulmonar de una rata. El sistema realizó la segmentación de imágenes (utilizando un algoritmo de caminante aleatorio), generó un modelo 3D y ejecutó un conjunto de simulación con tres miembros inicializados por datos de velocidad medidos.
• Análisis de Conjuntos: Visualización del conjunto utilizando gráficos de similitud 2D y comparaciones directas voxel a voxel. En un caso específico, el análisis reveló que una simulación destinada a replicar el flujo de la arteria pulmonar de una rata no logró capturar una estructura de vórtice crítica observada en las mediciones de MRI, clasificando correctamente la simulación como inadecuada.

Significado y Afirmaciones
El artículo afirma que, al integrar capacidades de CFD y MRI en una única herramienta de código abierto, el sistema proporciona una comprensión más completa de los parámetros hemodinámicos. Aborda específicamente la brecha en la investigación de simulaciones médicas respecto a la naturaleza engorrosa de generar y analizar conjuntos de simulaciones. Los autores afirman que su sistema facilita un análisis más preciso de biomarcadores hemodinámicos al permitir que investigadores con experiencia mínima previa en programación configuren, ajusten y visualicen datos de flujo complejos de manera interactiva. El trabajo se posiciona como un paso hacia la mejora de la reproducibilidad y la generalización en estudios hemodinámicos al reemplazar scripts personalizados con una aplicación unificada y documentada.