An open-source computational framework for immersed fluid-structure interaction modeling using FEBio and MFEM

Este trabajo presenta un marco de código abierto para la simulación de interacción fluido-estructura inmersa que acopla las bibliotecas MFEM y FEBio para abordar desafíos computacionales en sistemas biológicos, como la dinámica de válvulas cardíacas, mediante un esquema monolítico implícito y métodos de estabilización avanzados.

Lo esencial

Imagina que el cuerpo humano es como una ciudad muy compleja donde circula un río de sangre. A veces, en ese río, hay "puertas" (las válvulas del corazón) que se abren y cierran para dejar pasar el agua. El problema es que estas puertas son muy finas, se doblan muchísimo y a veces chocan entre sí.

¿Cuál es el problema con los programas actuales?
Antes, los científicos usaban programas para simular esto como si fueran moldes de gelatina. Tenían que hacer que la "gelatina" (la malla de la computadora) se moviera exactamente igual que la puerta. Pero cuando la puerta se dobla mucho o se cierra, la gelatina se estira, se rompe y se hace un desastre. Para arreglarlo, tenían que volver a hacer el molde desde cero una y otra vez, lo cual es lento, costoso y a veces imposible de hacer.

¿Qué propone este nuevo trabajo?
Los autores (un equipo de científicos de hospitales y universidades) han creado un nuevo "sistema de simulación" de código abierto (gratis para todos) que funciona de una manera totalmente diferente. En lugar de usar moldes que se deforman, usan una técnica de "inmersión".

Aquí tienes la analogía principal:

La Analogía del "Tanque de Agua" vs. El "Molde"

Imagina que quieres estudiar cómo se mueve una hoja de papel en un río.

• El método antiguo (ALE): Era como intentar que el agua solo existiera dentro de un molde de plástico que cambiaba de forma cada vez que la hoja se movía. Si la hoja se doblaba mucho, el molde se rompía.
• El nuevo método (Inmersión): Es como llenar una piscina grande con agua (el fluido) y simplemente soltar la hoja dentro. No importa si la hoja se dobla, gira o toca el fondo; el agua ya está ahí, ocupando todo el espacio. La hoja "nada" dentro del agua sin necesidad de que el agua cambie de forma para acomodarla.

¿Cómo funciona este nuevo sistema? (La "Alianza de Gigantes")

Para que esto funcione rápido y con precisión, los autores unieron dos programas de computadora muy potentes, como si fueran dos superhéroes uniéndose para salvar el día:

1. MFEM (El Motor de Velocidad): Es un programa creado en un laboratorio nacional de EE. UU. Es como un coche de carreras de Fórmula 1. Es increíblemente rápido, puede usar múltiples procesadores a la vez (incluso tarjetas gráficas de videojuegos) y maneja el "agua" (el fluido) con una velocidad increíble.
2. FEBio (El Experto en Tejidos): Es un programa creado por universidades para estudiar biología. Es como un cirujano experto. Sabe exactamente cómo se comportan los tejidos blandos, cómo se estiran, cómo se doblan y cómo reaccionan cuando chocan entre sí.

La Magia:
Antes, tenías que elegir entre ser rápido (MFEM) o ser preciso con los tejidos (FEBio). Este nuevo trabajo crea un puente entre ambos.

• El "agua" se calcula con la velocidad del coche de carreras (MFEM).
• La "hoja" o la "válvula" se calcula con la precisión del cirujano (FEBio).
• Hablan entre ellos en tiempo real para ver cómo el agua empuja a la hoja y cómo la hoja frena el agua.

¿Por qué es importante esto?

1. Corazones de niños: Los niños tienen corazones pequeños y válvulas que crecen con ellos. Los médicos necesitan simular cómo se comportarán las válvulas a medida que el niño crece para poder repararlas mejor. Este sistema permite hacer esos cálculos complejos sin que la computadora se "rompa".
2. Válvulas artificiales: Ayuda a diseñar mejores válvulas artificiales que duren más y no se rompan.
3. Gratuito y Abierto: Al ser de "código abierto", cualquier científico en el mundo puede usarlo, mejorarlo y adaptarlo. No hay que pagar licencias costosas.

En resumen

Este paper presenta una nueva herramienta gratuita que permite a los científicos simular cómo se mueve la sangre a través de válvulas cardíacas que se doblan y chocan, sin tener que preocuparse por que la simulación se rompa.

Es como pasar de intentar dibujar un río en un papel que se arruga cada vez que mueves una piedra, a poner la piedra en una piscina real y observar cómo el agua fluye a su alrededor de forma natural, rápida y precisa. Esto ayudará a los médicos a entender mejor las enfermedades del corazón y a diseñar mejores tratamientos para pacientes, especialmente niños.

Resumen técnico

Aquí presento un resumen técnico detallado del artículo en español, estructurado según los puntos solicitados:

Resumen Técnico: Marco de Interacción Fluido-Estructura (IFE) Inmerso de Código Abierto

1. El Problema

La simulación de la interacción fluido-estructura (IFE) en sistemas biológicos, particularmente en la dinámica de válvulas cardíacas, presenta desafíos computacionales masivos. Estos incluyen:

• Grandes deformaciones y contacto: Las válvulas cardíacas experimentan deformaciones extremas y contacto (coaptación) en milisegundos, lo que provoca cambios topológicos en el dominio fluido.
• Limitaciones de los métodos tradicionales: Los métodos Arbitrary Lagrangian-Eulerian (ALE), que utilizan mallas acopladas a la interfaz, sufren de distorsión de malla severa y requieren remallado frecuente (remeshing) o algoritmos complejos de movimiento de malla, lo que los hace costosos y propensos a errores en problemas con grandes movimientos o contacto.
• Falta de software de código abierto robusto: Las herramientas existentes de código abierto a menudo carecen de modelos de mecánica de sólidos sofisticados (necesarios para tejidos biológicos) o no están optimizadas para arquitecturas de computación de alto rendimiento (HPC) y aceleradores de hardware (GPUs). Además, los modelos simplificados (rígidos o de orden reducido) no capturan la mecánica tisular necesaria para estudiar la progresión de enfermedades o la durabilidad de dispositivos.

2. Metodología

El trabajo presenta un nuevo marco de trabajo de código abierto que acopla estratégicamente dos bibliotecas de elementos finitos maduras:

• MFEM: Una biblioteca escalable, lista para GPU y de alto rendimiento desarrollada en el Laboratorio Nacional Lawrence Livermore, encargada de resolver la parte del fluido.
• FEBio: Un solucionador de elementos finitos no lineales especializado en biomecánica (desarrollado en la Universidad de Utah y Columbia), encargado de modelar el sólido inmerso con modelos constitutivos complejos.

Enfoque Numérico:

• Método de Dominio Ficticio / Multiplicador de Lagrange Distribuido (FD/DLM): El sólido se "sumerge" en un dominio fluido de fondo fijo. No se requiere que la malla del fluido se ajuste a la interfaz.
• Acoplamiento Monolítico Implícito: Se utiliza un esquema totalmente implícito para resolver el sistema de ecuaciones acopladas simultáneamente, lo que proporciona robustez para interacciones fuertemente acopladas típicas de aplicaciones cardiovasculares.
• Estabilización VMS (Multiescala Variacional): Se emplea el método VMS para mejorar la precisión y estabilidad en mallas subresueltas (comunes en métodos inmersos). Se introduce un factor de escala adicional ($s$) en los coeficientes de estabilización para manejar los grandes saltos de presión a través de estructuras delgadas inmersas (como las válvulas).
• Algoritmo de Solución: Se utiliza el método $\alpha$ generalizado de segundo orden para la integración temporal. El sistema no lineal se resuelve con el método de Newton. Una característica clave es la eliminación de los multiplicadores de Lagrange mediante una definición nodal, reduciendo el sistema lineal resultante al tamaño de un problema de fluido puro, disminuyendo así el costo computacional.

3. Contribuciones Clave

1. Integración Sinérgica: Creación de un plugin (MFEMiFSI) que combina la infraestructura de computación paralela y GPU de MFEM con las capacidades avanzadas de modelado de tejidos blandos (hiperelasticidad, viscoelasticidad, contacto) de FEBio.
2. Arquitectura Modular: El diseño permite la extensión fácil a formulaciones de sólidos adicionales y nuevas capacidades de modelado sin modificar el núcleo del marco, siguiendo el patrón de plugins de FEBio.
3. Paralelización y Escalabilidad: Implementación de paralelización de memoria distribuida (MPI) y una ruta clara hacia la aceleración por GPU, permitiendo simulaciones de alta fidelidad en supercomputadoras.
4. Solución a Problemas de Contacto: El marco evita los problemas de remallado de los métodos ALE, permitiendo simular ciclos cardíacos completos con contacto de válvulas de manera eficiente.

4. Resultados

El marco fue validado mediante una serie de problemas de prueba numérica:

• Sólido Anular Inmerso: Demostró tasas de convergencia óptimas (orden 2.0 para velocidad en norma $L_2$) comparado con soluciones analíticas.
• Válvula Cardíaca Idealizada (Abierta y Cerrada):
  • En el caso abierto, los resultados de desplazamiento de la punta de la hoja coincidieron con métodos de referencia ALE-FSI.
  • En el caso cerrado, se demostró la capacidad del método para capturar discontinuidades de presión y estados hidrostáticos, algo crítico para válvulas cerradas, gracias a la modificación de los coeficientes de estabilización.
• Hoja Flexible Oscilante: Simuló grandes deformaciones complejas impulsadas por fuerzas fluidas en un flujo transversal pulsátil, un escenario difícil para métodos ALE.
• Esfera en Caída Libre: Validó la formulación en 3D, mostrando una convergencia hacia la solución analítica de velocidad terminal con un error del 3.02% en la malla más fina.
• Válvula Semilunar 3D: Se realizó una simulación completa de una válvula cardíaca trivalva en un aorta idealizada. El modelo capturó con éxito la dinámica de vórtices complejos, la formación de chorros y las distribuciones de tensión principal en las hojas durante un ciclo cardíaco completo, demostrando la capacidad para manejar geometrías realistas y grandes deformaciones.

5. Significancia

Este trabajo llena un vacío crítico en la comunidad de biomecánica computacional al proporcionar una plataforma de código abierto que une la mecánica de fluidos de alto rendimiento con la mecánica de sólidos biológicos avanzados.

• Impacto Clínico: Permite estudios integrados que analizan simultáneamente el rendimiento hemodinámico y la respuesta mecánica a nivel de tejido, esencial para optimizar reparaciones de válvulas nativas, predecir la durabilidad de válvulas bioprotésicas y diseñar dispositivos de reemplazo transcatéter.
• Accesibilidad y Futuro: Al ser de código abierto y modular, fomenta la colaboración comunitaria. La infraestructura está preparada para aprovechar la computación exascale y la aceleración por GPU, lo que reducirá los tiempos de simulación y permitirá análisis más detallados y personalizados para pacientes pediátricos y adultos.
• Disponibilidad: El código del plugin y los problemas de ejemplo están disponibles públicamente en los repositorios de FEBio y GitHub, facilitando su adopción inmediata por investigadores y clínicos.