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

Este estudo emprega simulações CFD-DEM resolvidas para demonstrar que nódulos polimetálicos não esféricos experimentam um arrasto significativamente aumentado e velocidades terminais reduzidas em comparação com esferas de volume equivalente devido à assimetria da esteira induzida pela forma, ao mesmo tempo em que revela como o tamanho da partícula e o confinamento governam comportamentos distintos de variância de arrasto durante o transporte hidráulico vertical.

O essencial

Imagine que você está tentando mover uma pilha de rochas estranhas e irregulares (chamadas nódulos polimetálicos) do fundo do oceano para um navio usando um canudo gigante e vertical. É assim que funciona a mineração em águas profundas. A grande questão para os engenheiros é: Como essas rochas de formatos estranhos se comportam quando a água corre para cima pelo canudo para carregá-las?

A maioria dos modelos computacionais usados para projetar esses sistemas trata cada rocha como se fosse uma bola de gude perfeita e lisa. Mas, na realidade, essas rochas são irregulares, cheias de calos e não se parecem nada com bolas de gude. Este artigo pergunta: Tratar uma rocha irregular como uma bola de gude lisa realmente funciona ou nos dá a resposta errada?

Para descobrir, os pesquisadores construíram uma simulação computacional super detalhada (como um motor de física de um videogame de alta tecnologia) que não utiliza atalhos. Em vez de apenas adivinhar como a água empurra a rocha, eles calcularam o empuxo da água em cada saliência e fenda da rocha.

Aqui está o que eles descobriram, explicado de forma simples:

1. A "Rocha Irregular" vs. A "Bola de Gude Lisa"

Quando os pesquisadores jogaram essas rochas irregulares em água parada para ver o quão rápido elas afundavam, as rochas irregulares caíram cerca de 28% mais devagar do que uma bola de gude lisa do exato mesmo peso e tamanho.

• A Analogia: Imagine tentar nadar através da água. Se você for um golfinho liso e aerodinâmico, você desliza facilmente. Se você for um pedaço de madeira irregular e serrilhado, você "prende" mais água no seu caminho para baixo.
• Por que isso acontece: As rochas irregulares têm uma "área frontal" maior (elas capturam mais água) e criam um rastro de esteira assimétrico e caótico atrás delas (como um rastro caótico de bolhas). Esse arrasto extra as retarda significamente.
• O Detalhe: Embora caiam mais devagar, a força total que as empurra para cima (empuxo) é a mesma da bola de gude. Elas apenas precisam se mover mais devagar para equilibrar essa força.

2. O "Engarrafamento" no Cano

Em seguida, eles simularam a água correndo para cima pelo cano para carregar essas rochas. Eles observaram dois tamanhos: rochas "Pequenas" e rochas "Grandes".

• As Bolas de Gude Lisas: Quando a velocidade da água aumentava, as bolas de gude lisas se comportavam de forma previsível. Em velocidades baixas, elas oscilavam e se acomodavam. Em velocidades altas, elas subiam rapidamente em linha reta, como carros entrando em uma rodovia.
• As Rochas Irregulares: Estas eram muito mais caóticas.
  • Em velocidades baixas: As pequenas rochas irregulares nem sequer conseguiam subir o cano! Elas pairavam perto do fundo, girando e oscilando no lugar, incapazes de vencer a gravidade. As bolas de gude lisas, no entanto, conseguiam subir.
  • Em velocidades altas: Mesmo quando a água era rápida o suficiente para carregá-las, as rochas irregulares levavam mais tempo para chegar ao topo e moviam-se em um caminho muito mais errático e giratório. Eram como um grupo de pessoas tentando subir uma escada rolante enquanto giram em círculos, enquanto as bolas de gude lisas estavam apenas correndo direto para cima.

3. O Efeito "Pião"

A maior diferença foi como as rochas rotacionavam.

• Bolas de Gude Lisas: Elas apenas subiam. Não giravam muito.
• Rochas Irregulares: Por serem irregulares, a água batendo nelas as fazia girar descontroladamente. Esse giro (rotação) estava intimamente ligado ao seu movimento para cima e para baixo.
• A Analogia: Pense em uma bola de gude lisa como uma bala disparada por uma arma — ela vai direto. Pense na rocha irregular como um bumerangue ou um pião lançado em um túnel de vento. Ela torce, vira e muda de direção constantemente devido à sua forma. Esse giro cria um "atrito" extra com a água, tornando o transporte mais difícil.

4. As "Flutuações de Força" (A Viagem Turbulenta)

Os pesquisadores mediram o "empurrão" (força de arrasto) que a água dava às rochas.

• Rochas Pequenas: Independentemente de serem lisas ou irregulares, o empurrão era relativamente constante.
• Rochas Grandes: Aqui é onde a coisa ficou selvagem.
  • Grandes Bolas de Gude Lisas: O empurrão variava um pouco conforme a água passava por elas, criando um padrão previsível de "solavancos" na força.
  • Grandes Rochas Irregulares: O empurrão era extremamente imprevisível. Como as rochas giravam e mudavam de forma em relação à água, a força dava picos repentinos. Era como dirigir um carro em uma estrada lisa (bolas de gude lisas) versus dirigir um carro em uma estrada onde os solavancos mudam a cada segundo, dependendo de como o carro está inclinado (rochas irregulares).

A Conclusão

O estudo conclui que, embora você possa usar um modelo de bola de gude lisa para ter uma ideia aproximada de como essas rochas se comportarão, ele perde os detalhes.

• Se você usar um modelo de bola de gude lisa, pode pensar que as rochas subirão o cano mais rápido e com mais facilidade do que realmente subirão.
• As rochas irregulares precisam de mais velocidade de água para começar a se mover e, uma vez em movimento, são menos estáveis e mais difíceis de controlar porque giram e oscilam.

Em resumo: A natureza é bagunçada. Você não pode simplesmente fingir que uma rocha irregular é uma esfera perfeita se quiser projetar uma máquina que realmente funcione. A "irregularidade" adiciona muito arrasto extra e caos que os modelos simples ignoram.

Resumo técnico

Resumo Técnico: Comportamento Hidrodinâmico de Partículas Não Esféricas em Fluxos Verticais Confinados

Definição do Problema
A recuperação de nódulos polimetálicos (PMNs) de planícies abissais profundas depende de sistemas hidráulicos que transportam partículas grosseiras, irregulares e não esféricas através de colunas de elevação (risers) verticais. A eficiência desses sistemas depende da previsão de arraste de partículas, estabilidade de suspensão e comportamento de transporte sob forte força gravitacional e confinamento geométrico. As práticas de projeto atuais dependem frequentemente de modelos de Dinâmica de Fluidos Computacional - Método de Elementos Discretos (CFD-DEM) não resolvidos, que utilizam correlações de arrasto empíricas calibradas para partículas esféricas ou parametrizadas por descritores escalares como a esfericidade. Essas abordagens falham em fluxos confinados (onde a razão diâmetro-partícula/diâmetro-tubo $d/D \ge 0,2$) porque não conseguem capturar o arrasto dependente da orientação, os rastros (wakes) assimétricos defletidos pela parede ou o movimento acoplado rotacional-translacional inerente às partículas não esféricas. Além disso, permanece incerto se a substituição das morfologias realistas dos PMNs por esferas de volume equivalente preserva as métricas de transporte em nível de conjunto ou enviesa as previsões de limiares de arraste e quedas de pressão.

Metodologia
Este estudo emprega uma estrutura CFD-DEM totalmente resolvida para investigar a hidrodinâmica de PMNs não esféricos em tubos verticais representativos de colunas de elevação de mineração submarina.

• Estrutura Numérica: A abordagem acopla as equações de Navier-Stokes incompressíveis (resolvidas via OpenFOAM) com as equações de Newton para o movimento de partículas discretas (resolvidas via LIGGGHTS).
• Acoplamento Fluido-Sólido: Um método de Fronteira Imersa (IB) é utilizado para resolver as interações fluido-sólido diretamente na interface sem leis de arrasto empíricas. A condição de não-deslizamento é imposta via uma abordagem de força direta, e as forças hidrodinâmicas são computadas integrando o tensor de tensão do fluido sobre a superfície da partícula.
• Representação de Partículas: Geometrias complexas de PMNs, reconstruídas a partir de tomografias computadorizadas (CT) de raios-X, são modeladas usando uma aproximação de multiesferas (conjuntos rígidos de subpartículas esféricas sobrepostas). Isso permite uma representação de alta fidelidade de superfícies irregulares mantendo a tratabilidade computacional.
• Verificação: A estrutura foi verificada contra dados experimentais para sedimentação de esfera única (Ten Cate et al., 2002) e o benchmark de drafting-kissing-tumbling (DKT) para múltiplas esferas (Glowinski et al., 2001; Sharma e Patankar, 2005). Adicionalmente, uma aproximação de multiesferas de uma esfera foi validada contra uma esfera lisa para confirmar a precisão da representação geométrica.
• Configuração da Simulação: As simulações foram conduzidas para PMNs e suas esferas de volume equivalente com diâmetros efetivos de 20 mm e 54 mm (diâmetros de volume equivalente de 16,4 mm e 44 mm) em um tubo de diâmetro $D=200$ mm. O estudo cobriu números de Reynolds de partícula ($Re_p$) de 98 a 2904 e razões de confinamento ($d/D$) de 0,082 e 0,22. O comportamento de transporte foi analisado através de velocidades de fluido variando de $1,0u_t$ a $3,0u_t$ (onde $u_t$ é a velocidade de sedimentação terminal).

Principais Resultados

1. Sedimentação e Aumento de Arrasto:

   • PMNs não esféricos sedimentam 27–29% mais devagar do que esferas de volume equivalente, apesar de possuírem massa e flutuabilidade idênticas.
   • Essa redução na velocidade terminal corresponde a um aumento do coeficiente de arrasto de 1,8 a 2,0 vezes em relação às esferas.
   • O aumento do arrasto decorre de dois fatores: (1) uma área frontal instantânea maior (até 1,5 vezes a área de volume equivalente devido à razão da esfera circunscrita) e (2) efeitos induzidos pela morfologia, incluindo estruturas de rastro assimétricas e acoplamento rotacional-translacional.
   • Crucialmente, a força de arrasto terminal permanece igual ao peso submerso para ambas as formas; a menor velocidade dos PMNs é resultado da geração da resistência hidrodinâmica necessária em velocidades de deslizamento (slip velocities) menores.

2. Tempo de Residência e Regimes de Transporte:

   • O comportamento de transporte exibe uma transição impulsionada pela velocidade de um movimento intermitente, dominado pela sedimentação, para um arraste estável, dominado pela convecção.
   • Em baixas velocidades ($u_f = 1,0u_t$), pequenos PMNs permanecem em um estado de suspensão marginal próxima à entrada, exibindo movimento oscilatório sem transporte ascendente líquido, enquanto as esferas mostram dispersão vertical.
   • À medida que a velocidade aumenta, as distribuições de tempo de residência estreitam-se para todas as partículas. No entanto, os PMNs exibem consistentemente tempos de residência médios e variâncias mais elevados em comparação com as esferas, atribuídos ao maior arrasto e ao forte acoplamento rotacional-translacional.
   • Em altas velocidades ($u_f = 3,0u_t$), todas as partículas alcançam distribuições espacialmente homogêneas, embora os PMNs mantenham movimento rotacional aumentado (velocidades angulares de até 80 rad/s) e trajetórias helicoidais.

3. Estatísticas de Força de Arrasto e Flutuações:

   • A força de arrasto normalizada média ($\hat{f}$) permanece próxima da unidade em todas as condições, confirmando um equilíbrio de estado estacionário entre o arrasto e o peso submerso.
   • Partículas Pequenas ($d/D = 0,082$): Tanto as esferas quanto os PMNs exibem distribuições de força de arrasto estreitas. Para pequenos PMNs, a reorientação rápida efetivamente neutraliza as flutuações da área projetada dependentes da orientação, tornando seu comportamento quantitativamente semelhante ao das esferas.
   • Partículas Grandes ($d/D = 0,22$):
     • Grandes esferas mostram um alargamento monotônico da distribuição de arrasto com o aumento da velocidade, impulsionado pela instabilidade do rastro e efeitos de confinamento.
     • Grandes PMNs exibem distribuições significativamente mais amplas e assimétricas com caudas de alto arrasto estendidas. Isso é impulsionado pela superposição de efeitos de histórico de rastro e força dependente da orientação. Diferente das partículas pequenas, as escalas de tempo rotacionais dos grandes PMNs são comparáveis às escalas de tempo de evolução do rastro, impedindo a neutralização das variações da área projetada.

Significância e Conclusões
O estudo elucida os mecanismos pelos quais a forma da partícula influencia o arrasto, a dinâmica do rastro e o acoplamento rotacional-translacional em fluxos verticais confinados. A principal contribuição é estabelecer os limites quantitativos para a aplicação de modelos esféricos de volume equivalente ao transporte de PMNs não esféricos.

• Similaridade Qualitativa: A física de transporte subjacente (progressão de regimes dominados por sedimentação para regimes dominados por convecção) é qualitativamente semelhante tanto para PMNs quanto para esferas. Isso sugere que comportamentos de transporte de primeira ordem podem ser capturados usando modelos esféricos com leis de arrasto calibradas.
• Divergência Quantitativa: Existem diferenças quantitativas significativas. Os PMNs exibem tempos de residência 40–90% maiores e coeficientes de arrasto 1,8–2,0 vezes superiores aos das esferas. Essas diferenças são críticas para prever as velocidades mínimas de suspensão e as quedas de pressão em colunas de elevação de mineração submarina.
• Implicações de Modelagem: Embora aproximações esféricas possam ser suficientes para a modelagem de ordem reduzida do comportamento de conjunto, elas falham em capturar a intermitência específica e as flutuações de força associadas a partículas não esféricas grandes em razões de confinamento elevadas. A abordagem CFD-DEM resolvida apresentada aqui fornece uma base necessária para o desenvolvimento de correções específicas de forma para tais modelos.

Os autores observam que estas descobertas se aplicam a conjuntos diluídos e servem como entradas de calibração; a extrapolação para polpas densas envolvendo sedimentação dificultada ou fluxo de plugue permanece fora do escopo deste trabalho.