Hydrodynamic Behavior of Non-spherical Particles in Confined Vertical Flows: A Resolved CFD-DEM Study

Este estudio emplea simulaciones CFD-DEM resueltas para demostrar que los nódulos polimetálicos no esféricos experimentan una resistencia significativamente mayor y velocidades terminales reducidas en comparación con esferas de volumen equivalente debido a la asimetría de la estela inducida por la forma, al tiempo que revela cómo el tamaño de la partícula y el confinamiento gobiernan comportamientos distintos de la varianza de la resistencia durante el transporte hidráulico vertical.

Lo esencial

Imagina que estás intentando mover un montón de rocas extrañas y deformes (llamadas nódulos polimetálicos) desde el fondo del océano hasta un barco utilizando un gigantesco tubo vertical. Así es como funciona la minería en aguas profundas. La gran pregunta para los ingenieros es: ¿Cómo se comportan estas rocas extrañas cuando el agua corre hacia arriba por el tubo para transportarlas?

La mayoría de los modelos informáticos utilizados para diseñar estos sistemas tratan cada roca como si fuera una canica perfecta y lisa. Pero en la realidad, estas rocas son rugosas, irregulares y no se parecen en nada a las canicas. Este artículo pregunta: ¿Funciona realmente el tratar una roca deforme como si fuera una canica lisa, o nos está dando la respuesta equivocada?

Para averiguarlo, los investigadores construyeron una simulación informática súper detallada (como un motor de física de un videojuego de alta tecnología) que no utiliza atajos. En lugar de adivinar cómo el agua empuja la roca, calcularon el empuje del agua sobre cada bulto y grieta de la roca.

Aquí está lo que descubrieron, explicado de forma sencilla:

1. La "Roca Deforme" frente a la "Canica Lisa"

Cuando los investigadores dejaron caer estas rocas deformes en agua quieta para ver qué tan rápido se hundían, las rocas deformes cayeron aproximadamente un 28% más lento que una canica lisa del mismo peso y tamaño exactos.

• La Analogía: Imagina que intentas nadar a través del agua. Si eres un delfín suave y aerodinámico, te deslizas fácilmente. Si eres un trozo de madera rugoso y dentado, atrapas más agua en tu camino hacia abajo.
• Por qué sucede: Las rocas deformes tienen un "área frontal" mayor (atrapan más agua) y crean una estela asimétrica y desordenada detrás de ellas (como un rastro caótico de burbujas). Este arrastre adicional las ralentiza significamente.
• El Detalle: Aunque caen más lento, la fuerza total que las empuja hacia arriba (flotabilidad) es la misma que la de la canica. Simplemente necesitan moverse más lento para equilibrar esa fuerza.

2. El "Atasco" en la Tubería

A continuación, simularon el agua corriendo hacia arriba por el tubo para transportar estas rocas. Observaron dos tamaños: rocas "Pequeñas" y rocas "Grandes".

• Las Canicas Lisas: Cuando la velocidad del agua aumentaba, las canicas lisas se comportaban de manera predecible. A velocidades bajas, se tambaleaban y se asentaban. A velocidades altas, subían a toda velocidad en línea recta, como coches incorporándose a una autopista.
• Las Rocas Deformes: Estas eran mucho más caóticas.
  • A velocidades bajas: ¡Las pequeñas rocas deformes ni siquiera lograban subir por el tubo! Se quedaban flotando cerca del fondo, girando y tambaleándose en su lugar, incapaces de vencer la gravedad. Las canicas lisas, sin embargo, lograban subir.
  • A velocidades altas: Incluso cuando el agua era lo suficientemente rápida como para transportarlas, las rocas deformes tardaban más tiempo en llegar a la parte superior y se movían siguiendo una trayectoria mucho más errática y giratoria. Eran como un grupo de personas intentando subir una escalera mecánica mientras giran en círculos, mientras que las canicas lisas simplemente subían corriendo en línea recta.

3. El Efecto "Peonza"

La mayor diferencia fue cómo rotaban las rocas.

• Canicas Lisas: La mayoría de las veces solo subían. No giraban mucho.
• Rocas Deformes: Debido a que son rugosas, el agua al golpearlas las hacía girar salvajemente. Este giro (rotación) estaba estrechamente vinculado a su movimiento hacia arriba y hacia abajo.
• La Analogía: Piensa en una canica lisa como una bala disparada por un arma: va recta. Piensa en la roca deforme como un bumerán o un trompo lanzado a un túnel de viento. Gira, cambia de dirección y se retuerce constantemente debido a su forma. Este giro crea una "fricción" adicional con el agua, dificultando su transporte.

4. Las "Fluctuaciones de Fuerza" (El Viaje Accidentado)

Los investigadores midieron el "empuje" (fuerza de arrastre) que el agua ejercía sobre las rocas.

• Rocas Pequeñas: Ya fueran lisas o deformes, el empuje era relativamente constante.
• Rocas Grandes: Aquí es donde la cosa se volvió salvaje.
  • Canicas Lisas Grandes: El empuje variaba un poco a medida que el agua pasaba rápidamente a su alrededor, creando un patrón predecible de "baches" en la fuerza.
  • Rocas Deformes Grandes: El empuje era extremadamente impredecible. Debido a que las rocas giraban y cambiaban de orientación respecto al agua, la fuerza sufría picos repentinos. Era como conducir un coche en una carretera lisa (canicas lisas) frente a conducir un coche en una carretera donde los baches cambian cada segundo dependiendo de la inclinación del coche (rocas deformes).

La Conclusión

El estudio concluye que, si bien se puede usar un modelo de canica lisa para tener una idea general de cómo se comportarán estas rocas, este omite los detalles.

• Si utilizas un modelo de canica lisa, podrías pensar que las rocas subirán por el tubo más rápido y con mayor facilidad de lo que realmente lo harán.
• Las rocas deformes necesitan más velocidad de agua para ponerse en movimiento y, una vez que se mueven, son menos estables y más difíciles de controlar porque giran y se tambalean.

En resumen: La naturaleza es desordenada. No puedes simplemente pretender que una roca dentada es una esfera perfecta si quieres diseñar una máquina que realmente funcione. La "rugosidad" añade mucho arrastre y caos adicionales que los modelos simples ignoran.

Resumen técnico

Resumen Técnico: Comportamiento Hidrodinámico de Partículas No Esféricas en Flujos Verticales Confinados

Planteamiento del Problema
La recuperación de nódulos polimetálicos (PMN) de las llanuras abisales de aguas profundas depende de sistemas hidráulicos que transportan partículas gruesas, irregulares y no esféricas a través de elevadores verticales. La eficiencia de estos sistemas depende de la predicción de la captación de partículas, la estabilidad de la suspensión y el comportamiento de transporte bajo una fuerte fuerza gravitatoria y confinamiento geométrico. Las prácticas de diseño actuales suelen depender de modelos de Dinámica de Fluidos Computacional y Método de Elementos Discretos (CFD-DEM) no resueltos que utilizan correlaciones de arrastre empíricas calibradas para partículas esféricas o parametrizadas mediante descriptores escalares como la esfericidad. Estos enfoques fallan en flujos confinados (donde las relaciones de diámetro partícula-tubo son $d/D \ge 0.2$) porque no pueden capturar el arrastre dependiente de la orientación, las estelas asimétricas desviadas por la pared o el acoplamiento rotacional-traslacional inherente a las partículas no esféricas. Además, sigue sin estar claro si la sustitución de las morfologías realistas de los PMN por esferas de volumen equivalente preserva las métricas de transporte a nivel de conjunto o sesga las predicciones de los umbrales de captación y las caídas de presión.

Metodología
Este estudio emplea un marco CFD-DEM totalmente resuelto para investigar la hidrodinámica de los PMN no esféricos en tuberías verticales representativas de los elevadores de minería en aguas profundas.

• Marco Numérico: El enfoque acopla las ecuaciones de Navier-Stokes incompresibles (resueltas mediante OpenFOAM) con las ecuaciones de Newton para el movimiento de partículas discretas (resueltas mediante LIGGGHTS).
• Acoplamiento Fluido-Sólido: Se utiliza un método de Frontera Inmersa (IB) para resolver las interacciones fluido-sólido directamente en la interfaz sin leyes de arrastre empíricas. La condición de no deslizamiento se impone mediante un enfoque de fuerza directa y las fuerzas hidrodinámicas se calculan integrando el tensor de esfuerzo del fluido sobre la superficie de la partícula.
• Representación de Partículas: Las geometrías complejas de los PMN, reconstruidas a partir de tomografías computarizadas de rayos X (CT), se modelan mediante una aproximación de multiesferas (ensamblajes rígidos de subpartículas esféricas superpuestas). Esto permite una representación de alta fidelidad de las superficies irregulares manteniendo la viabilidad computacional.
• Verificación: El marco fue verificado mediante datos experimentales para la sedimentación de una sola esfera (Ten Cate et al., 2002) y el benchmark de "drafting-kissing-tumbling" (DKT) para múltiples esferas (Glowinski et al., 2001; Sharma y Patankar, 2005). Además, una aproximación de multiesferas de una esfera fue validada contra una esfera lisa para confirmar la precisión de la representación geométrica.
• Configuración de la Simulación: Las simulaciones se realizaron para PMN y sus esferas de volumen equivalente con diámetros efectivos de 20 mm y 54 mm (diámetros de volumen equivalente de 16.4 mm y 44 mm) en una tubería de diámetro $D=200$ mm. El estudio cubrió números de Reynolds de partícula ($Re_p$) de 98 a 2904 y relaciones de confinamiento ($d/D$) de 0.082 y 0.22. El comportamiento de transporte se analizó a través de velocidades de fluido que variaron desde $1.0u_t$ hasta $3.0u_t$ (donde $u_t$ es la velocidad de sedimentación terminal).

Resultos Clave

1. Sedimentación y Mejora del Arrastre:

   • Los PMN no esféricos sedimentan entre un 27 y un 29% más lento que las esferas de volumen equivalente, a pesar de tener la misma masa y flotabilidad.
   • Esta reducción en la velocidad terminal corresponde a una mejora del coeficiente de arrastre de 1.8 a 2.0 veces la de las esferas.
   • El aumento del arrastre surge de dos factores: (1) un área frontal instantánea mayor (hasta 1.5 veces el área de volumen equivalente debido a la relación de la esfera circunscrita) y (2) efectos inducidos por la morfología, incluyendo estructuras de estela asimétricas y el acoplamiento rotacional-traslacional.
   • Crucialmente, la fuerza de arrastre terminal permanece igual al peso sumergido para ambas formas; la menor velocidad de los PMN es el resultado de generar la resistencia hidrodinámica requerida a velocidades de deslizamiento menores.

2. Tiempo de Residencia y Regímenes de Transporte:

   • El comportamiento de transporte exhibe una transición impulsada por la velocidad, desde un movimiento intermitente dominado por la sedimentación hacia una captación estable dominada por la convección.
   • A velocidades bajas ($u_f = 1.0u_t$), los PMN pequeños permanecen en un estado de suspensión marginal cerca de la entrada, exhibiendo un movimiento oscilatorio sin transporte ascendente neto, mientras que las esferas muestran dispersión vertical.
   • A medida que la velocidad aumenta, las distribuciones de tiempo de residencia se estrechan para todas las partículas. Sin embargo, los PMN exhiben consistentemente tiempos de residencia medios más altos y mayores varianzas en comparación con las esferas, lo que se atribuye al mayor arrastre y al fuerte acoplamiento rotacional-traslacional.
   • A altas velocidades ($u_f = 3.0u_t$), todas las partículas alcanzan distribuciones espacialmente homogéneas, aunque los PMN mantienen un movimiento rotacional mejorado (velocidades angulares de hasta 80 rad/s) y trayectorias helicoidales.

3. Estadísticas y Fluctuaciones de la Fuerza de Arrastre:

   • La fuerza de arrastre normalizada media ($\hat{f}$) permanece cerca de la unidad en todas las condiciones, confirmando un equilibrio de estado estacionario entre el arrastre y el peso sumergido.
   • Partículas Pequeñas ($d/D = 0.082$): Tanto las esferas como los PMN exhiben distribuciones de fuerza de arrastre estrechas. Para los PMN pequeños, la reorientación rápida promedia eficazmente las fluctuaciones del área proyectada dependientes de la orientación, haciendo que su comportamiento sea cuantitativamente similar al de las esferas.
   • Partículas Grandes ($d/D = 0.22$):
     • Las esferas grandes muestran un ensanchamiento monotónico de la distribución de arrastre con el aumento de la velocidad, impulsado por la inestabilidad de la estela y los efectos de confinamiento.
     • Los PMN grandes exhiben distribuciones significativamente más amplias y sesgadas con colas de alto arrastre extendidas. Esto es impulsado por la superposición de los efectos de la historia de la estela y la fuerza dependiente de la orientación. A diferencia de las partículas pequeñas, las escalas de tiempo rotacionales de los PMN grandes son comparables a las escalas de tiempo de evolución de la estela, lo que impide el promedio de las variaciones del área proyectada.

Significancia y Conclusiones
El estudio elucida los mecanismos por los cuales la forma de la partícula influye en el arrastre, la dinámica de la estela y el acoplamiento rotacional-traslacional en flujos verticales confinados. La contribución principal es el establecimiento de los límites cuantitativos para la aplicación de modelos esféricos de volumen equivalente al transporte de PMN no esféricos.

• Similitud Cualitativa: La física del transporte subyacente (progresión desde regímenes de sedimentación hacia regímenes dominados por la convección) es cualitativamente similar tanto para los PMN como para las esferas. Esto sugiere que los comportamientos de transporte de primer orden pueden capturarse utilizando modelos esféricos con leyes de arrastre calibradas.
• Divergencia Cuantitativa: Existen diferencias cuantitativas significativas. Los PMN exhiben tiempos de residencia entre un 40 y un 90% más largos y coeficientes de arrastre de 1.8 a 2.0 veces mayores que las esferas. Estas diferencias son críticas para predecir las velocidades mínimas de suspensión y las caídas de presión en los elevadores de minería en aguas profundas.
• Implicaciones de Modelado: Aunque las aproximaciones esféricas pueden ser suficientes para el modelado de orden reducido del comportamiento del conjunto, fallan en capturar la intermitencia específica y las fluctuaciones de fuerza asociadas con partículas grandes y no esféricas en relaciones de alto confinamiento. El enfoque CFD-DEM resuelto presentado aquí proporciona una base necesaria para el desarrollo de correcciones específicas de forma para tales modelos.

Los autores señalan que estos hallazgos se aplican a conjuntos diluidos y sirven como entradas de calibración; la extrapolación a lodos densos que involucren sedimentación impedida o flujo de tapón permanece fuera del alcance de este trabajo.