Kinetic derivation of thermal viscous models for nematic liquid crystal dynamics

Resumen Técnico: Derivación Cinética de Modelos Térmicos Viscosos para la Dinámica de Cristales Líquidos Nemáticos

Planteamiento del Problema
Trabajos recientes establecieron una variante compresible e invíscida de las ecuaciones de Leslie–Ericksen para cristales líquidos nemáticos derivada de la teoría cinética. Sin embargo, ese enfoque dependía de un cierre de momentos alrededor de un estado ya alineado, imponiendo efectivamente un comportamiento de orientación colectiva en lugar de permitir que este emergiera de la dinámica. Además, el modelo resultante era estrictamente invíscido, fallando en capturar los efectos disipativos, térmicos y de difusión de espín inherentes a los flujos reales de cristales líquidos. Este artículo aborda estas limitaciones mediante la derivación de una variante viscosa y compresible de las ecuaciones de Leslie–Ericksen que tiene en cuenta la emergencia de la alineación e incluye efectos térmicos y viscosos.

Metodología
Los autores desarrollan una teoría termodinámica macroscópica partiendo de una teoría cinética de fluidos ordenados. La metodología avanza a través de las siguientes etapas:

1. Marco Cinético: El punto de partida es la teoría cinética de fluidos ordenados propuesta en [7], la cual incorpora interacciones de alineación de campo medio mediante un potencial de Vlasov. El sistema asume una separación de escalas de tiempo donde la relajación orientacional ocurre más rápido que la relajación del momento traslacional. La función de distribución $f(x, p, \nu, \varsigma, t)$ depende de la posición, el momento lineal, la configuración microestructural (orientación $\nu$) y el momento conjugado $\varsigma$.
2. Operador de Colisión: Para hacer tratable el análisis, el operador de colisión de Boltzmann completo se sustituye por un único operador de relajación de tipo Bhatnagar–Gross–Krook (BGK). Este operador preserva las leyes de conservación y la estructura de producción de entropía de la teoría cinética completa.
3. Leyes de Balance: Los autores derivan las ecuaciones de balance macroscópicas para la masa, el momento lineal, el momento angular, la energía y la entropía mediante la toma de momentos de la ecuación cinética. Una característica clave es la identificación de un término de torque de campo medio en el balance del momento angular, el cual impulsa la emergencia de la alineación.
4. Estructura Termodinámica: La energía libre de Helmholtz se identifica como el potencial termodinámico. Los autores demuestran que la distribución de equilibrio (Maxwelliana) minimiza esta energía libre y derivan la distribución de Gibbs asociada para el perfil orientacional. Demuestran que, bajo potenciales de interacción específicos (tipos Kuramoto y Maier–Saupe), el sistema experimenta una transición de fase de un estado isotrópico a uno nemático a medida que la temperatura disminuye, con el torque de Vlasov desapareciendo en el equilibrio alineado.
5. Expansión de Chapman–Enskog: Se realiza una expansión de Chapman–Enskog de primer orden alrededor del equilibrio de Maxwell local. A diferencia de las derivaciones estándar, las derivadas temporales en la expansión se sustituyen utilizando las leyes de balance de orden cero (que incluyen el estrés de Ericksen y términos de esfuerzo de par), en lugar de las ecuaciones de Euler compresibles estándar. Esto asegura que la elasticidad orientacional se propague correctamente hacia las correcciones de primer orden.
6. Cierre Constitutivo: Los autores emplean el principio de maximización restringida de Rajagopal y Srinivasa. Combinan la tasa de producción de entropía derivada de la expansión cinética con la estructura termodinámica inferida de la energía libre de Helmholtz. Al maximizar la producción de entropía sujeta a las leyes de balance, derivan las relaciones constitutivas para el tensor de Cauchy, el tensor de esfuerzo de par y el flujo de calor.

Contribuciones Clave y Resultados

• Alineación Emergente: A diferencia de los modelos invíscidos previos que asumían la alineación, esta derivación permite que la transición isotrópico-nemática emerja naturalmente de la dinámica cinética a través del torque de Vlasov. El análisis confirma que el torque de campo medio desaparece en el equilibrio de Gibbs alineado, asegurando la consistencia con las leyes de balance macroscópicas.
• Efectos Viscosos y Térmicos: El artículo deriva un conjunto completo de ecuaciones constitutivas para un fluido nemático compresible, viscoso, conductor de calor y con difusión de espín. El tensor de estrés resultante incluye:
  • Una viscosidad de cizalladura Newtoniana.
  • Una viscosidad volumétrica isotrópica que surge del desequilibrio de energía traslacional-rotacional.
  • Un estrés viscoso anisotrópico acoplado al tensor de Ericksen y al gradiente de velocidad.
  • Un tensor de esfuerzo de par que incluye efectos de difusión de espín (un término micropolar de tipo Eringen) y acoplamiento al gradiente del director.
• Coeficientes de Transporte: Se derivan expresiones explícitas para los coeficientes de transporte (viscosidad de cizalladura, conductividad térmica, viscosidad volumétrica, viscosidad anisotrópica y viscosidad de espín) en términos del tiempo de relajación $\tau$ del BGK y el módulo de Frank $K$. La conductividad térmica satisface la relación de Eucken para gases diatómicos.
• Límite Incompresible: El marco produce naturalmente el contraparte incompresible mediante la imposición de una restricción de divergencia nula vía un multiplicador de Lagrange durante el procedimiento de maximización, sin requerir una derivación separada.
• Novedad Estructural: El tensor de estrés viscoso derivado está gobernado por un único parámetro de viscosidad (debido a la aproximación BGK isotrópica) y no recupera el tensor viscoso de Leslie con sus cinco coeficientes completos. Los autores señalan que recuperar el tensor de Leslie completo requeriría un operador BGK anisotrópico con múltiples escalas de tiempo de relajación, lo cual es objeto de trabajo continuo.

Significancia y Reivindicaciones
El artículo afirma generalizar la reciente teoría cinética invíscida de Farrell, Russo y Zerbinati [1] para dar cuenta de los efectos viscosos, térmicos y de difusión de espín. La principal significancia radica en proporcionar una derivación termodinámicamente consistente de estos efectos disipativos directamente desde la teoría cinética, cerrando la brecha entre las interacciones de partículas microscópicas y la mecánica de medios continuos macroscópica para fluidos ordenados.

Los autores enfatizan que su enfoque supera la limitación de los modelos previos al no imponer la alineación a priori; en su lugar, la alineación emerge de la minimización de la energía libre y del desvanecimiento del torque de Vlasov en el equilibrio. El sistema resultante (Ecuación 8.21) representa una variante compresible, viscosa, conductora de calor y con difusión de espín de las ecuaciones de Leslie–Ericksen. En el límite invíscido ($\tau \to 0$), el modelo se reduce al sistema invíscido derivado en [1].

El trabajo se presenta como un paso hacia una teoría cinética más completa de los cristales líquidos, con los autores señalando modestamente que la recuperación del tensor viscoso de Leslie completo con sus cinco viscosidades independientes sigue siendo un problema abierto que requiere un operador de colisión más complejo y anisotrópico.