Enucleated incompressible red blood cells in shear flow: theoretical analysis of shape instabilities

Resumen Técnico: Glóbulos Rojos Desnucleados e Incompresibles en Flujo de Cizalla: Análisis Teórico de Inestabilidades de Forma

Planteamiento del Problema
Los glóbulos rojos (RBC, por sus siglas en inglés) son esenciales para el transporte de oxígeno, dependiendo de su capacidad para experimentar deformaciones significativas en el flujo. Mientras que los modelos de cuerpo rígido (por ejemplo, órbitas de Jeffery) describen con éxito la dinámica de los RBC en tasas de cizalla bajas (tumbling/giro, tank-treading/bandeo, swinging/oscilación), no logran explicar la emergencia de formas complejas de múltiples lóbulos y morfologías de ruptura de simetría (como los estomatocitos) observadas en tasas de cizalla intermedias ($100\text{--}1000 \text{ s}^{-1}$). Los mecanismos físicos que impulsan estas inestabilidades de forma —específicamente la interacción entre la elasticidad de la membrana, la resistencia a la flexión, la tensión de la membrana, el exceso de área y los cocientes de viscosidad— siguen siendo objeto de debate. Además, los marcos teóricos existentes suelen tratar al RBC como un elipsoide o cápsula estática, omitiendo los cambios geométricos específicos resultantes del proceso de maduración de la desnucleación, donde el núcleo de la célula es expulsado, alterando la relación entre el área de la membrana libre de tensiones y el volumen encerrado.

Metodología
Los autores emplean un marco teórico perturbativo que modela al RBC como una membrana cerrada, llena de fluido e incompresible, suspendida en un flujo de cizalla Newtoniano. El enfoque extiende el marco no lineal linealizado establecido previamente por Vlahovska et al. [26] al incorporar explícitamente dos áreas de exceso distintas:

1. $\Delta_0$ (Área de exceso libre de tensiones): Correspondiente a la forma elipsoidal del RBC previo a la desnucleación.
2. $\Delta$ (Área de exceso actual): Correspondiente a la forma de discoide tras la desnucleación, donde el volumen disminuye ($V < V_0$) mientras que el área superficial se conserva, lo que genera tensiones internas bloqueadas.

La expansión de la forma del RBC se realiza mediante armónicos esféricos ($Y_{l,m}$), separando la configuración actual $\rho$ y la configuración libre de tensiones $\rho_0$. Las ecuaciones de movimiento gobernantes se derivan del equilibrio entre los esfuerzos del fluido (flujo de Stokes) y las fuerzas de la membrana, las cuales se descomponen en contribuciones elásticas (modelo de Skalak), de flexión (energía de Helfrich) y de tensión de membrana. La tensión de la membrana ($\Sigma$) actúa como un multiplicador de Lagrange para imponer un área local constante y acopla los diferentes modos espaciales.

El análisis procede en dos etapas:

1. Sistema Reducido ($l=2$): El sistema se reduce primero a modos elipsoidales para revisar la dinámica de tumbling, tank-treading y swinging, contabilizando explícitamente el desajuste entre $\Delta$ y $\Delta_0$.
2. Sistema Perturbativo Completo ($l \geq 2$): Se introducen modos espaciales adicionales para derivar un criterio de inestabilidad lineal. La estabilidad del sistema se analiza examinando la tasa de crecimiento de los modos de orden superior ($l > 2$) como función de parámetros adimensionales: el número de flexión ($B$), el número capilar ($C$), el cociente de viscosidad ($\lambda$) y las áreas de exceso.

Contribuciones Clave y Resultados

• Extensión del Marco Teórico: El estudio generaliza el enfoque perturbativo existente para dar cuenta explícitamente de los cambios geométricos y de estado de tensión inducidos por la desnucleación. Al distinguir entre la forma libre de tensiones ($\Delta_0$) y la forma actual ($\Delta$), el modelo captura la "memoria" del proceso de maduración de la célula, una característica ausente en los modelos de vesículas fluidas.
• Papel del Área de Exceso en la Dinámica Reducida: En el subespacio $l=2$, los autores demuestran que el desajuste entre $\Delta$ y $\Delta_0$ renormaliza la constante elástica efectiva ($S$). Esto permite transiciones dinámicas (por ejemplo, de swinging a tumbling) impulsadas por cambios en el área de exceso de la forma de referencia, incluso a tasas de cizalla constantes. La orientación inicial (en el plano o fuera del plano) influye significativamente en las trayectorias resultantes.
• Criterio de Inestabilidad y Tensión de Membrana: La contribución principal es la derivación de un criterio de inestabilidad lineal para las transiciones de forma. El análisis revela que las inestabilidades de forma (la emergencia de formas de múltiples lóbulos) dependen críticamente de la tensión de la membrana $\sigma(t)$. Una inestabilidad ocurre cuando la tensión de la membrana cae por debajo de un umbral crítico ($\sigma < \sigma^*$), determinado por la relación entre la elasticidad de cizalla y la rigidez de flexión ($\alpha = B/C$).
  • La inestabilidad se inicia cuando la respuesta de flexión y elástica es insuficiente para contrarrestar el flujo externo y los esfuerzos internos.
  • El modo más inestable ($l^*$) y la tensión crítica dependen de $\alpha$, el cociente de viscosidad $\lambda$ y las áreas de exceso.
• Saturación y Redistribución: El estudio analiza el comportamiento a largo plazo tras una inestabilidad. El crecimiento de los modos inestables provoca un aumento en la tensión de la membrana, lo que eventualmente estabiliza el sistema. Este proceso conduce a una redistribución del área de exceso desde el modo dominante $l=2$ hacia modos de orden superior (por ejemplo, $l=3, 4$), resultando en una nueva configuración de forma estable y no elipsoidal.
• Diagramas de Fase: La integración numérica de las ecuaciones acopladas valida las predicciones analíticas. Los autores construyen diagramas de fase en el plano $(B, C)$, identificando regiones donde modos específicos ($l=3, 4, \dots$) se vuelven inestables. Estos diagramas muestran que el aumento del cociente de viscosidad $\lambda$ y el área de exceso $\Delta$ desestabiliza a los RBC en números capilares altos (elasticidad débil), promoviendo la emergencia de formas complejas.

Significancia
El artículo proporciona una explicación teórica para la emergencia de morfologías complejas de RBC, tales como estomatocitos y formas de trilobulados, bajo flujo de cizalla externo. Al vincular estas inestabilidades con las consecuencias mecánicas específicas de la desnucleación (el desajuste entre el área de exceso libre de tensiones y la actual), el trabajo ofrece una comprensión mecanicista de cómo los RBC transicionan de dinámicas elipsoidales simples a estados complejos de múltiples lóbulos. El marco identifica la tensión de la membrana como el motor clave que acopla los modos espaciales y destaca el papel crítico de la relación entre el módulo de flexión y el de cizalla, así como la viscosidad, en determinar el inicio de las inestabilidades de forma. Esta visión teórica cierra la brecha entre los modelos de cuerpo rígido y las simulaciones numéricas completas, ofreciendo una visión analítica de las condiciones bajo las cuales los RBC pierden su simetría elipsoidal.