On the development of OpenFOAM solvers for simulating MHD micropolar fluid flows with or without the effect of micromagnetorotation

Este trabajo presenta el desarrollo y validación de dos nuevos solvers de OpenFOAM, *epotMicropolarFoam* y *epotMMRFoam*, para simular flujos magnetohidrodinámicos de fluidos micropolares, demostrando que la inclusión del efecto de micromagnetorotación (MMR) es crucial para modelar con precisión la reducción de velocidad y la estabilización de vórtices en flujos sanguíneos bajo campos magnéticos, lo cual tiene importantes aplicaciones biomédicas.

Lo esencial

¡Claro que sí! Imagina que este artículo científico es como la historia de dos nuevos "superhéroes" del mundo de la simulación por computadora, creados para entender cómo se comporta la sangre cuando la sometemos a fuertes campos magnéticos.

Aquí tienes la explicación en español, sencilla y con analogías divertidas:

🩸 El Problema: La Sangre no es solo Agua

Imagina que la sangre es como una sopa de fideos. Los "fideos" son las células rojas (eritrocitos). En la física clásica, a veces tratamos la sangre como si fuera agua simple (un fluido "newtoniano"), donde todo fluye suavemente. Pero en realidad, la sangre es más compleja: es un fluido micropolar.

¿Qué significa esto? Imagina que cada célula roja no es solo un punto que se mueve, sino una pequeña hélice o un tornillo que, además de avanzar, puede girar sobre sí mismo. Cuando la sangre fluye, estas "hélices" giran y chocan entre sí.

🧲 El Nuevo Factor: La "Bailarina Magnética" (Micromagnetorrotación)

Aquí entra la magia del artículo. Cuando aplicas un imán fuerte a la sangre (como en una resonancia magnética o en terapias médicas), las células rojas, que tienen hierro, sienten el imán.

El artículo introduce un concepto llamado Micromagnetorrotación (MMR).

• La analogía: Imagina que tienes un grupo de bailarines (las células) girando en una pista. De repente, un director de orquesta (el campo magnético) les grita: "¡Dejen de girar y miren hacia mí!".
• El efecto: Las células intentan alinearse con el imán y dejan de girar libremente. Esto crea una especie de "fricción magnética" interna. El artículo descubre que si ignoramos este efecto, nuestros cálculos están muy equivocados. ¡Es como si el director de orquesta no existiera y los bailarines siguieran girando descontroladamente!

💻 Los Dos "Superhéroes" (Los Solvers)

Los autores crearon dos programas informáticos (llamados solvers) dentro de una herramienta gratuita llamada OpenFOAM (piensa en ella como un "Lego" gigante para construir simulaciones de fluidos):

1. epotMicropolarFoam: Es el primer héroe. Simula la sangre como un fluido con "hélices" que giran, bajo un campo magnético, pero ignora el efecto de que el imán las obligue a dejar de girar. Es como simular a los bailarines sin el director de orquesta.
2. epotMMRFoam: Es el héroe avanzado. Incluye todo lo anterior más la "Micromagnetorrotación". Aquí, el programa sabe que el imán está forzando a las células a alinearse y dejar de girar. Es el héroe que ve la realidad completa.

🔬 ¿Qué descubrieron? (Los Resultados)

Los científicos probaron estos programas en dos escenarios:

1. En una arteria recta (como un tubo de agua):

• Sin MMR: El imán apenas hacía nada. La sangre fluía casi igual que sin imán.
• Con MMR: ¡Cambio drástico! Cuando se activó el efecto de "alinear a las células", la velocidad de la sangre bajó hasta un 40% y el giro de las células casi desapareció (un 99.9%).
• La lección: Si no tienes en cuenta que el imán "congela" el giro de las células, subestimas mucho cómo el imán frena la sangre.

2. En un aneurisma (un globo en la arteria):
Imagina un aneurisma como un globo hinchado en medio de un tubo. Normalmente, el agua entra, se expande y crea remolinos o torbellinos dentro del globo (zonas de recirculación).

• Sin MMR: Los remolinos eran grandes y desordenados.
• Con MMR: ¡Los remolinos desaparecieron o se hicieron muy pequeños! El efecto de alinear las células con el imán actuó como un amortiguador o un "freno de estabilidad". La sangre se volvió más ordenada y estable dentro del globo.

🏥 ¿Por qué importa esto? (Aplicaciones Reales)

Esto no es solo teoría. Imagina que quieres enviar un medicamento a un tumor usando un imán para guiarlo (terapia dirigida).

• Si usas los programas viejos (sin MMR), pensarás que el imán no hace mucho.
• Si usas los nuevos programas (con MMR), verás que el imán puede frenar la sangre y controlar mejor dónde van las células o los medicamentos.

🚀 Conclusión Simple

Los autores crearon dos nuevas herramientas de computadora que nos dicen la verdad: cuando usas imanes con la sangre, no solo empujas las células, sino que también les quitas su capacidad de girar. Ignorar este "giro magnético" es como intentar predecir el clima sin tener en cuenta el viento.

Estos programas son un gran paso para mejorar tratamientos médicos, como la hipertermia magnética (calentar tumores con imanes) o la entrega precisa de fármacos, asegurando que los médicos entiendan realmente cómo se moverá la sangre bajo un imán.

Resumen técnico

Resumen Técnico: Desarrollo de Solvers OpenFOAM para Flujos MHD de Fluidos Micropolares con Efecto de Micromagnetorrotación

1. Planteamiento del Problema

Los fluidos que contienen partículas magnéticas, como la sangre (debido a la hemoglobina en los eritrocitos) y los ferrofluidos, experimentan un par magnético cuando se someten a un campo magnético externo. Este fenómeno, conocido como micromagnetorrotación (MMR), surge del desalineamiento entre la magnetización de las partículas y el campo magnético aplicado.

Aunque la teoría de fluidos micropolares (que considera la rotación de partículas individuales) es crucial para modelar la sangre, el efecto de la MMR a menudo se ignora en las simulaciones de flujo sanguíneo bajo campos magnéticos. La literatura existente carece de códigos numéricos de código abierto capaces de simular flujos magnetohidrodinámicos (MHD) de fluidos micropolares que incluyan tanto el acoplamiento de la fuerza de Lorentz como el término de par magnético (MMR) y la ecuación constitutiva de magnetización. Además, se ha observado experimentalmente que los campos magnéticos fuertes afectan la velocidad y la viscosidad de la sangre de una manera que la fuerza de Lorentz sola no puede explicar debido a la baja conductividad eléctrica de la sangre.

2. Metodología

Los autores desarrollaron dos nuevos solucionadores transitorios basados en la biblioteca de código abierto OpenFOAM (escrito en C++):

• epotMicropolarFoam: Diseñado para flujos MHD micropolares incompresibles y laminares sin el efecto de MMR.
• epotMMRFoam: Una extensión del anterior que incluye el efecto de MMR, resolviendo la ecuación de magnetización constitutiva y el par magnético en la ecuación de momento angular interno.

Aspectos Clave de la Implementación:

• Algoritmo: Ambos solucionadores utilizan el algoritmo PISO (Pressure-Implicit with Splitting of Operators) para el acoplamiento presión-velocidad.
• Aproximación MHD: Se adopta la aproximación de número de Reynolds magnético bajo ($Re_m \ll 1$), ignorando la ecuación de inducción magnética y utilizando una formulación de potencial eléctrico.
• Ecuaciones Gobernantes:
  • Se resuelve la ecuación de conservación de momento lineal, incluyendo el término de diferencia entre microrrotación y vorticidad ($2\mu_r \nabla \times (\omega - \Omega)$) y la fuerza de Lorentz.
  • Se resuelve la ecuación de cambio de momento angular interno. En epotMMRFoam, esta ecuación incluye el término de par magnético ($M \times H$) y se acopla con la ecuación de magnetización de Shizawa-Tanahashi.
• Validación: Los solucionadores se validaron contra soluciones analíticas para el flujo de Poiseuille MHD micropolar con diferentes valores de hematocrito (25% y 45%) e intensidades de campo magnético (1 T, 3 T y 8 T).
• Casos de Estudio: Se realizaron simulaciones en geometrías complejas:
  1. Flujo en una arteria 3D cilíndrica.
  2. Flujo a través de un aneurisma 2D fusiforme (con dilataciones del 166% y 200%).

3. Contribuciones Clave

• Desarrollo de Software: Creación de los primeros solucionadores de código abierto en OpenFOAM (epotMicropolarFoam y epotMMRFoam) capaces de manejar flujos MHD micropolares completos con el modelo de magnetización de Shizawa-Tanahashi.
• Inclusión de MMR: Implementación rigurosa del término de micromagnetorrotación, permitiendo estudiar cómo la alineación de las partículas magnéticas afecta la dinámica del flujo, algo que los modelos MHD clásicos no capturan.
• Validación Numérica: Demostración de una alta precisión (error < 2%) al comparar los resultados numéricos con soluciones analíticas para flujos de Poiseuille.
• Repositorio Público: Los códigos están disponibles públicamente bajo licencia GPLv3, facilitando su uso en investigación biomédica y de ingeniería.

4. Resultados Principales

• Precisión de Validación: Ambos solucionadores mostraron una excelente concordancia con las soluciones analíticas. Los errores máximos para la velocidad y la microrrotación no superaron el 2% en los casos de Poiseuille.
• Impacto de la MMR vs. Fuerza de Lorentz:
  • Sin MMR: El efecto del campo magnético en la sangre fue mínimo, incluso con campos fuertes (hasta 8 T), debido a la baja conductividad eléctrica de la sangre y el pequeño tamaño de los vasos. La fuerza de Lorentz por sí sola no alteró significativamente la velocidad o la microrrotación.
  • Con MMR: Se observaron reducciones drásticas en la velocidad (hasta un 40-42%) y en la microrrotación (hasta un 99.9%) bajo campos magnéticos intensos y altos valores de hematocrito.
• Mecanismo Físico: La inclusión de MMR alinea los eritrocitos con el campo magnético externo, "congelando" su rotación interna. Esto aumenta la resistencia efectiva al flujo y suprime la formación de núcleos de recirculación.
• Estudios de Aneurisma:
  • En aneurismas (166% y 200% de dilatación), la ausencia de MMR permitió la formación de grandes zonas de recirculación y vórtices intensos.
  • Con MMR, estas zonas de recirculación se suprimieron significativamente, los vórtices se volvieron más localizados y el flujo se estabilizó. Esto demuestra un efecto de amortiguamiento de cizalladura y estabilización del flujo.
• Dependencia del Hematocrito: Los efectos de la MMR fueron más pronunciados a mayor hematocrito (45%), ya que hay más partículas magnéticas (eritrocitos) presentes para interactuar con el campo.

5. Significado e Implicaciones

Este trabajo es fundamental para la ingeniería biomédica y la investigación de fluidos complejos:

• Aplicaciones Clínicas: Los solucionadores ofrecen herramientas precisas para estudiar aplicaciones como la hipertermia magnética (tratamiento de cáncer) y la administración dirigida de fármacos, donde el control preciso del flujo sanguíneo bajo campos magnéticos es crucial.
• Precisión en Modelado: Demuestra que ignorar la MMR en fluidos biológicos bajo campos magnéticos fuertes conduce a modelos inexactos que no capturan la física real de la interacción partícula-campo.
• Futuro: Aunque el estudio se limitó a flujos laminares y partículas rígidas, los solucionadores desarrollados sientan las bases para futuras investigaciones que incluyan flujos pulsátiles, viscosidad no newtoniana, deformabilidad de eritrocitos y geometrías específicas de pacientes.

En conclusión, el estudio confirma que el efecto de la micromagnetorrotación es dominante sobre la fuerza de Lorentz en el flujo sanguíneo bajo campos magnéticos, y proporciona la infraestructura computacional necesaria para explorar estas interacciones en escenarios biomédicos realistas.