A Volume of Fluid Immersed Boundary Method for Industrial Polymer Mixing

Este artículo presenta un nuevo solucionador acoplado por bloques de Volumen de Fluido con Frontera Inmersa (BC-VOF-IB) implementado en OpenFOAM que supera las inestabilidades numéricas causadas por altos contrastes de viscosidad para simular con precisión la mezcla de polímeros con superficie libre en extrusores industriales parcialmente llenos.

Lo esencial

Imagina que estás intentando mezclar un tanque gigante de miel espesa y pegajosa con aire dentro de una máquina giratoria. Esto es esencialmente lo que sucede en la mezcla industrial de polímeros, donde empresas como Pirelli necesitan combinar plástico fundido con aditivos para fabricar neumáticos, dispositivos médicos o piezas de automóviles. El objetivo es lograr que todo se mezcle perfectamente para que el producto final sea resistente y uniforme.

Sin embargo, simular este proceso en una computadora es una pesadilla para matemáticos e ingenieros. Aquí está el porqué, y cómo este artículo lo resuelve, utilizando analogías simples:

El Problema: La Lucha entre "Miel Espesa vs. Aire Delgado"

En estas máquinas, tienes dos fluidos muy diferentes:

1. Fusión de Polímero: Extremadamente espesa, pegajosa y de movimiento lento (como miel fría).
2. Aire: Muy delgado y de movimiento rápido.

Cuando intentas simular cómo interactúan estos dos dentro de una máquina con tornillos giratorios, los programas informáticos estándar se confunden. Es como intentar calcular el movimiento de un caracol y un coche de carreras en la misma pista usando el mismo conjunto de reglas. La computadora intenta dar pasos diminutos, diminutos, para evitar que el "caracol" (el plástico espeso) se mueva demasiado rápido, lo que hace que la simulación sea increíblemente lenta; a veces tarda días en completar unos pocos segundos de mezcla en tiempo real.

Además, las máquinas tienen partes complejas y giratorias (tornillos) que se mueven dentro de un contenedor fijo. Tradicionalmente, para simular esto, tienes que construir una malla digital (una cuadrícula de cajas diminutas) que envuelva perfectamente los tornillos giratorios. A medida que los tornillos giran, esta cuadrícula debe remodelarse constantemente, lo cual es como intentar tejer un suéter mientras la persona que lo lleva está corriendo un maratón. Es desordenado, difícil y propenso a errores.

La Solución: Una Nueva "Cuadrícula Inteligente" y un "Enfoque de Equipo"

Los autores de este artículo desarrollaron una nueva forma de ejecutar estas simulaciones utilizando un software llamado OpenFOAM. Combinaron dos técnicas poderosas:

1. El Método de Frontera Inmersa (El Truco de la "Pared Fantasma")
En lugar de remodelar la cuadrícula para que se ajuste a los tornillos giratorios, mantuvieron la cuadrícula fija y rígida (como un bloque sólido de hielo). Luego le dijeron a la computadora: "Oye, hay un tornillo giratorio dentro de este bloque de hielo".

• La Analogía: Imagina una piscina con una cuadrícula fija de baldosas en el fondo. En lugar de mover las baldosas para ajustarlas a un nadador, simplemente le dices al agua: "No pases a través del nadador". La computadora utiliza matemáticas para crear una "pared fantasma" alrededor del tornillo, forzando al fluido a fluir alrededor de él sin necesidad de reconstruir nunca la cuadrícula. Esto hace que manejar formas complejas y en movimiento sea mucho más fácil.

2. El Método de Volumen de Fluido (VOF) (El Truco de la "Pintura de Rastreo")
Para ver dónde termina el plástico espeso y dónde comienza el aire, utilizan una "pintura" que llena las celdas.

• La Analogía: Imagina que la cuadrícula de la computadora es un tablero de ajedrez 3D. Algunos cuadros son 100% plástico, otros son 100% aire, y algunos son una mezcla. La computadora rastrea cuánta "pintura de plástico" hay en cada cuadro para ver la superficie del líquido.

3. El Esquema Acoplado por Bloques (El "Huddle de Equipo")
Este es el avance más importante. En las simulaciones estándar, la computadora resuelve la velocidad del fluido en las direcciones X, Y y Z una por una, como tres personas hablando por turnos. Cuando el fluido es súper espeso (como el polímero), este enfoque de "tomar turnos" hace que la simulación se bloquee o se ralentice hasta un arrastre porque el fluido espeso acopla todas las direcciones juntas de manera estrecha.

Los autores cambiaron esto a un enfoque Acoplado por Bloques.

• La Analogía: En lugar de que tres personas tomen turnos, todas se reúnen en un círculo y resuelven el problema juntas exactamente al mismo tiempo. Al tratar el movimiento en todas las direcciones como un solo equipo gigante e interconectado, la computadora puede manejar la enorme diferencia entre el plástico espeso y el aire delgado sin quedarse atascada.

Los Resultados: De Horas a Minutos

El equipo probó su nuevo método en dos escenarios:

1. Un Canal con Forma de Hueso de Perro: Un caso de prueba donde el plástico espeso se inyecta en un canal estrecho y sinuoso.

   • Método Antiguo: El programa informático estándar se bloqueó o tardó 7 horas en simular unos pocos segundos porque se vio obligado a dar pasos diminutos.
   • Nuevo Método: Su nuevo método de "Huddle de Equipo" completó el mismo trabajo en solo 16 minutos y no se bloqueó, incluso cuando el plástico se volvió extremadamente espeso.

2. Máquinas Industriales Reales: Simularon extrusoras de tornillo único y tornillo doble del mundo real (las máquinas utilizadas para fabricar gránulos de plástico).

   • Mostraron con éxito cómo el plástico llena la máquina, cómo se acumula la presión y cómo el aire es expulsado.
   • Demostraron que su método de "Pared Fantasma" funciona tan bien como el antiguo y difícil método de remodelar la cuadrícula, pero mucho más rápido y más fácil de configurar.

¿Qué Sigue?

El artículo concluye que este es un gran paso adelante para la industria. Cierra la brecha entre las matemáticas académicas y las necesidades reales de la fábrica. Sin embargo, los autores señalan que su modelo actual asume que la temperatura se mantiene constante (isotérmico). En la realidad, mezclar plástico genera calor, lo cual cambia qué tan espeso es el plástico. Agregar efectos de temperatura y comportamientos de plástico más complejos y "elásticos" son los siguientes pasos para la investigación futura.

En resumen: Construyeron una forma más rápida y estable de observar cómo el plástico espeso se mezcla con el aire en máquinas giratorias dentro de una computadora, convirtiendo un proceso que antes tardaba horas en uno que toma minutos, sin necesidad de reconstruir el mundo digital cada vez que un tornillo gira.

Resumen técnico

Resumen Técnico: Un Método de Frontera Inmersa con Volumen de Fluido para la Mezcla Industrial de Polímeros

Enunciado del Problema
La simulación numérica de procesos de mezcla de polímeros presenta desafíos significativos debido a la compleja interacción entre la reología no newtoniana, la dinámica de superficies libres multifásicas y geometrías complejas en movimiento. Las mezcladoras industriales, como los extrusores de husillo simple (SSE) y de husillo doble (TSE), a menudo operan bajo condiciones de llenado parcial donde coexisten fundidos de polímero altamente viscosos con aire. Simular estos escenarios requiere capturar la interfaz entre fases mientras se manejan fronteras sólidas en rotación.

Los enfoques estándar de Dinámica de Fluidos Computacional (CFD) enfrentan dos obstáculos principales en este contexto:

1. Complejidad Geométrica: La generación de mallas adaptadas al cuerpo (conformes) para husillos en rotación con geometrías complejas es computacionalmente costosa y a menudo requiere técnicas de malla dinámica que pueden conducir a distorsión de la malla.
2. Inestabilidad Numérica: El contraste extremo de viscosidad entre los fundidos de polímero y el aire (que a menudo supera de $10^5$ a $10^8$) causa inestabilidades numéricas severas en los solucionadores estándar de Volumen de Fluido (VOF) segregados. Específicamente, el tratamiento explícito del gradiente de velocidad transpuesto en el término de difusión viscosa impone restricciones severas al paso de tiempo ($\Delta t \propto h^2$), haciendo que las simulaciones de flujos de alta viscosidad sean computacionalmente prohibitivas o inestables.

Metodología
Este trabajo propone un marco numérico novedoso que integra un enfoque de captura de interfaz de Volumen de Fluido (VOF) con un método de Frontera Inmersa (IB) no conforme, implementado dentro de la biblioteca de Volúmenes Finitos de OpenFOAM.

• Enfoque de Frontera Inmersa (IB): El método utiliza una cuadrícula de fondo fija donde las fronteras sólidas (por ejemplo, husillos en rotación) se representan mediante superficies STL trianguladas. Se emplea un operador de interpolación de Mínimos Cuadrados Ponderados (WLS) para imponer condiciones de frontera (Dirichlet para la velocidad, Neumann para la fracción volumétrica) en las celdas intersectadas por la frontera inmersa. Esto elimina la necesidad de generación de mallas conformes y permite el manejo de líneas de contacto móviles a lo largo de geometrías complejas y en rotación.
• Formulación VOF: La interfaz se rastrea resolviendo una ecuación de transporte para la fracción volumétrica de fase ($\alpha$). El método emplea el algoritmo MULES (Limitador Universal Multidimensional con Solución Explícita) para mantener la nitidez y acotación de la interfaz, incluyendo un término de velocidad compresiva para contrarrestar la difusión numérica.
• Esquema Acoplado por Bloques (BC): Para abordar los problemas de estabilidad que surgen de los altos contrastes de viscosidad, los autores introducen un esquema acoplado por bloques para la ecuación de momento. A diferencia de los solucionadores segregados estándar que tratan las componentes de velocidad secuencialmente, el enfoque BC trata el término completo de difusión viscosa $\nabla \cdot (\mu(\nabla \mathbf{u} + \nabla^\top \mathbf{u}))$ de forma implícita. Esto implica acoplar las tres componentes de velocidad dentro de un único sistema lineal. La discretización del término de gradiente tangencial utiliza un método de Área Finita (FA) combinado con interpolación WLS en plantillas extendidas (que involucran vecinos puntuales en lugar de solo vecinos de cara), adaptado de aplicaciones de elasticidad [11].
• Algoritmo de Solución: El solucionador emplea un algoritmo PIMPLE modificado (combinando SIMPLE y PISO) para manejar el acoplamiento presión-velocidad. Las restricciones IB se integran en el sistema algebraico, modificando los coeficientes de la matriz y los términos fuente para satisfacer las condiciones de frontera en las superficies inmersas mientras se mantiene la conservación de masa.

Contribuciones Clave

1. Desarrollo de Solucionadores BC-VOF-IB: El artículo introduce el solucionador UCInterIbFoam, que combina el método IB no conforme con el enfoque VOF acoplado por bloques. Esta es la primera implementación de tal método acoplado dentro del marco OpenFOAM-10 para flujos bifásicos.
2. Mejora de la Estabilidad mediante Acoplamiento por Bloques: El trabajo demuestra que el tratamiento totalmente implícito de la difusión viscosa relata significativamente las restricciones de estabilidad del paso de tiempo. Esto permite la simulación de flujos de polímeros de alta viscosidad con relaciones de viscosidad de hasta $10^8$ sin que la simulación falle debido a violaciones del número de Courant.
3. Manejo de Líneas de Contacto Móviles: El marco aborda con éxito las dificultades geométricas y algorítmicas de las líneas de contacto móviles donde la interfaz bifásica interseca fronteras inmersas, un requisito crítico para dispositivos de mezcla llenados parcialmente.

Resultados Numéricos
Los métodos propuestos fueron validados y probados en dos categorías de problemas:

• Caso de Referencia (Moldeo por Inyección de Hueso de Perro):

  • Se realizó una comparación entre solucionadores segregados (interFoam), conformes acoplados por bloques (UCInterFoam) y no conformes acoplados por bloques (UCInterIbFoam) en un caso de inyección de fluido de ley de potencia de alta viscosidad.
  • El solucionador segregado estándar no logró producir resultados estables dentro de plazos razonables debido a restricciones severas del paso de tiempo ($\Delta t \approx 10^{-6}$ s).
  • Los solucionadores acoplados por bloques lograron soluciones estables con pasos de tiempo de hasta $10^{-3}$ s, reduciendo el tiempo computacional de 7 horas a 16 minutos (conforme) y 58 minutos (no conforme).
  • La solución IB no conforme mostró acuerdo cualitativo con la solución conforme, con discrepancias menores en los picos de viscosidad cerca de restricciones geométricas, atribuidas a la naturaleza no conforme de la cuadrícula.

• Aplicaciones Industriales (SSE y TSE):

  • Se realizaron simulaciones en Extrusores de Husillo Simple (SSE) simplificados y Extrusores de Husillo Doble (TSE) realistas bajo condiciones tanto de llenado completo como de llenado parcial (alimentación por escasez).
  • Los solucionadores capturaron con éxito la evolución de la superficie libre, las relaciones de llenado y los campos de presión.
  • Los resultados se alinearon con las expectativas físicas: los escenarios de llenado completo ($Q_{in} > Q_0$) exhibieron gradientes de presión negativos desde la entrada hasta la salida, mientras que los escenarios de llenado parcial ($Q_{in} < Q_0$) mostraron aumentos de presión cerca de la salida.
  • El método IB demostró ser efectivo para manejar la compleja geometría de engranaje interpenetrado del TSE sin requerir regeneración dinámica de la malla.

Significado y Afirmaciones
El artículo afirma que el marco BC-VOF-IB propuesto representa un paso significativo hacia la superación de la brecha entre la investigación académica en CFD y las demandas del procesamiento industrial de polímeros. Al superar las inestabilidades numéricas asociadas con los altos contrastes de viscosidad y eliminar la necesidad de una compleja generación de mallas adaptadas al cuerpo para geometrías en rotación, el método permite predicciones físicamente consistentes de los campos de velocidad y presión en dispositivos de mezcla industrialmente relevantes y llenados parcialmente.

Los autores reconocen que el trabajo actual se limita a procesos isotérmicos y reología newtoniana generalizada. Declaran que los desarrollos futuros deben incluir efectos térmicos (ecuación de energía) y modelos constitutivos viscoelásticos para capturar completamente la física de la mezcla industrial de polímeros. Sin embargo, el marco actual establece una base rigurosa para estas extensiones y proporciona una herramienta robusta para simular flujos de superficie libre en geometrías complejas y en movimiento.