Mixing of miscible liquids: Dimensionless scaling for intermediate-to-large density differences in a stirred tank

Este estudo utiliza simulações numéricas de um tanque agitado com líquidos miscíveis para demonstrar que, embora o tempo de mistura corra positivamente com o número de Richardson, uma lei de escala exponencial derivada baseada nos números de Potência, Froude e Richardson colapsa com sucesso todos os dados em uma única curva mestra para diferenças de densidade intermediárias a grandes.

O essencial

Imagine que você está tentando misturar um tanque gigante de dois líquidos diferentes: um xarope pesado e grosso no fundo e um suco mais leve e fino no topo. Você solta uma pá giratória gigante (um rotor) para misturá-los.

No mundo real, esta é uma tarefa comum em fábricas que produzem tudo, desde medicamentos até tratamento de águas residuais. Mas aqui está o problema: como os líquidos têm pesos diferentes (densidades), o pesado quer ficar no fundo e o leve quer flutuar no topo. Isso cria uma "luta" entre a pá giratória tentando misturá-los e a gravidade tentando mantê-los separados.

Este artigo é como uma história de detetive onde cientistas usaram simulações computacionais poderosas para descobrir exatamente quanto tempo leva para misturar perfeitamente esses dois líquidos e como prever esse tempo sem precisar construir um tanque físico e realizar testes caros a cada vez.

O Cenário: Uma Cozinha de Teste Digital

Os pesquisadores construíram uma versão virtual de um tanque de mistura industrial padrão.

• O Tanque: É um grande cilindro com paredes e quatro aletas verticais (defletores) para impedir que o líquido apenas gire em círculo como um rio preguiçoso.
• A Pá: Uma lâmina giratória no meio.
• Os Líquidos: Eles simularam uma mistura 50/50 de um líquido pesado e um líquido mais leve. Eles não usaram produtos químicos reais; apenas trataram-nos como fluidos "pesados" e "leves" com a mesma espessura (viscosidade).
• O Método: Em vez de usar equações matemáticas padrão, eles usaram um truque inteligente chamado Método de Boltzmann de Rede. Pense nisso como simular o líquido não como uma massa contínua, mas como bilhões de pequenas bolas de bilhar invisíveis quicando e colidindo. Isso permitiu que eles vissem exatamente como a turbulência (o redemoinho caótico) se comportava.

A Grande Pergunta: Quão Rápido Podemos Misturar?

O objetivo principal era encontrar uma "fórmula mágica" para prever o tempo de mistura.

• As Variáveis: Eles alteraram duas coisas principais:
  1. A velocidade de rotação da pá (número de Reynolds): Girar mais rápido geralmente significa mais turbulência e mistura mais rápida.
  2. Quão diferentes são os pesos (número de Richardson): Se os líquidos têm quase o mesmo peso, misturam-se facilmente. Se um é muito mais pesado, a gravidade luta contra a mistura, criando camadas difíceis de romper.

A Descoberta: A Batalha "Gravidade vs. Rotação"

Os pesquisadores encontraram alguns padrões interessantes:

1. Quando a Gravidade Não Importa (Mesmo Peso):
   Se os dois líquidos têm exatamente o mesmo peso, o tempo de mistura é surpreendentemente consistente. Não importa quão rápido você gire a pá (dentro de uma certa faixa), o "tempo de mistura adimensional" (uma maneira chique de dizer "quantas voltas da pá são necessárias") permanece constante em cerca de 20 voltas. É como uma regra da natureza: uma vez que a água está agitada o suficiente, girá-la mais rápido não faz com que ela misture mais rápido em termos de voltas da pá.

2. Quando a Gravidade Luta de Volta (Pesos Diferentes):
   Assim que os líquidos têm pesos diferentes, o pesado quer ficar no fundo. Quanto maior a diferença de peso, mais difícil é misturar.

   • A Tendência: Quanto mais diferentes os pesos, mais tempo leva para misturar.
   • O Giro Surpreendente: Se você mantiver a "diferença de peso" constante e apenas girar a pá mais rápido, o tempo de mistura nem sempre diminui. Às vezes, girar mais rápido na verdade faz com que demore mais para atingir um ponto específico de mistura.
   • Por quê? Imagine que o líquido pesado é como um cobertor grosso. Se você girar a pá muito rápido, cria muita energia, mas o líquido pesado forma uma "tampa" estável que o líquido mais leve não consegue penetrar. A energia é desperdiçada ao redemoinhar a camada superior enquanto a camada inferior permanece selada. É como tentar mexer uma panela de sopa onde os vegetais pesados se assentaram em um bloco sólido no fundo; girar a colher mais rápido apenas salpica o caldo no topo sem quebrar o bloco de vegetais.

A Solução: Uma Nova "Curva Mestra"

A maior conquista da equipe foi criar uma única fórmula simples que combina todos esses fatores. Eles descobriram que, se você olhar para o tempo de mistura através da lente de três números específicos (Potência, Froude e Richardson), todos os seus dados desordenados colapsam em uma única curva suave e exponencial.

Pense assim: Antes, os engenheiros tinham que adivinhar ou realizar centenas de testes para ver como um novo líquido se misturaria. Agora, eles têm uma "receita". Se você disser a eles a diferença de peso e a velocidade planejada de rotação, esta fórmula prevê o tempo de mistura com alta precisão.

A Conclusão

O artigo conclui que, para esses tanques industriais específicos:

• A turbulência é fundamental: Uma vez que o líquido está totalmente agitado, o comportamento da mistura é previsível.
• A gravidade é o chefe: Se os líquidos têm densidades diferentes, a gravidade cria uma "estratificação" (camadas) que resiste à mistura.
• Mais rápido nem sempre é melhor: Em sistemas com grandes diferenças de densidade, simplesmente aumentar a velocidade do motor não garante uma mistura mais rápida; às vezes, isso apenas cria uma separação mais estável.

Os autores fornecem esta nova fórmula para ajudar os engenheiros a projetar melhores processos de mistura sem precisar construir protótipos caros primeiro. Eles planejam testar esta fórmula em diferentes formas de tanque e tipos de pá no futuro, mas, por enquanto, funciona perfeitamente para os tanques padrão que simularam.

Resumo técnico

Resumo Técnico: Escalonamento Adimensional para Diferenças de Densidade Intermediárias a Grandes em um Tanque Agitado

Declaração do Problema
A mistura de líquidos miscíveis é uma operação unitária crítica em processos industriais, variando da produção farmacêutica ao tratamento de águas residuais. Embora a mistura turbulenta em sistemas monofásicos e homogêneos seja bem caracterizada, a presença de diferenças de densidade entre os fluidos introduz estratificação impulsionada pela flutuabilidade. Este fenômeno cria regiões estratificadas verticalmente onde as forças gravitacionais contrabalançam o transporte convectivo, potencialmente suprimindo a mistura vertical e aumentando o tempo necessário para alcançar a homogeneidade. Embora existam correlações empíricas para mistura de traçadores monofásicos e casos estratificados específicos, um framework geralmente válido que considere consistentemente os efeitos gravitacionais, inerciais, viscosos e de flutuabilidade através de diferentes regimes de fluxo e geometrias permanece elusivo. Este estudo aborda o desafio de prever tempos de mistura em tanques agitados com bafles contendo dois líquidos miscíveis com diferenças de densidade intermediárias a grandes.

Metodologia
Os autores empregaram Dinâmica dos Fluidos Computacional (CFD) usando o Método de Boltzmann em Rede (LBM) para simular a mistura turbulenta em um tanque agitado com bafles. O sistema consistiu em uma razão volumétrica de 50/50 de dois líquidos miscíveis com viscosidades dinâmicas idênticas, mas densidades variáveis ($\rho_l$ e $\rho_h$), onde o fluido mais leve foi inicialmente estratificado acima do fluido mais pesado.

• Abordagem Numérica: As simulações resolveram as equações de Navier-Stokes incompressíveis usando um esquema de rede D3Q19. Para lidar com altos números de Reynolds ($Re > 10^5$), um modelo de Simulação de Grandes Vórtices (S-L LES) Smagorinsky-Lilly foi implementado para modelar a turbulência sub-grade.
• Transporte Escalar: Um modelo de transporte escalar baseado em fluxo de massa foi adotado para rastrear a concentração do líquido mais leve. Este modelo acoplou o campo escalar à equação de momento através de uma força de corpo gravitacional dependente da concentração, permitindo a simulação de efeitos de flutuabilidade decorrentes de variações locais de densidade.
• Espaço de Simulação: Um estudo sistemático de parâmetros foi conduzido variando a velocidade de rotação do agitador ($N$) e a densidade da fase mais leve ($\rho_l$). Isso resultou em uma ampla gama de números adimensionais: números de Reynolds ($Re$) de 50.000 a 125.000 e números de Richardson ($Ri$) de 0 a 14.
• Validação: A precisão numérica foi verificada através de uma análise do Índice de Convergência de Malha (GCI), que estimou a incerteza numérica em aproximadamente 2,39% na malha de produção. Adicionalmente, os resultados da simulação foram validados contra dados experimentais de traçadores da literatura, mostrando boa concordância na resposta da sonda e nas definições de tempo de mistura.

Principais Resultados
O estudo quantificou o tempo de mistura ($t_m$), definido como o tempo necessário para alcançar 90% de homogeneidade, e expressou-o em forma adimensional ($t^*_m = N \cdot t_m$).

1. Base Homogênea ($Ri = 0$): Na ausência de diferenças de densidade, o tempo de mistura adimensional colapsou para um valor constante de $t^*_m \approx 20$ em toda a faixa de números de Reynolds investigada. Isso confirma que o fluxo é totalmente turbulento e que a mistura inercial domina, consistente com a literatura estabelecida para tanques com bafles.
2. Efeito da Estratificação ($Ri > 0$): À medida que o número de Richardson aumentou, indicando maiores diferenças de densidade, o tempo de mistura adimensional aumentou sistematicamente. Os dados seguiram uma tendência exponencial clara em todos os números de Reynolds.
3. Lei de Escalonamento Proposta: Os autores derivaram uma curva mestra que colapsa todos os dados de simulação em uma única relação de escalonamento exponencial:
   $$t^*_m = 20 \cdot e^{0.58 \cdot N_p \cdot Fr \cdot Ri}$$
   Onde $N_p$ é o número de potência, $Fr$ é o número de Froude e $Ri$ é o número de Richardson. O pré-fator de 20 representa a base inercial, enquanto o termo exponencial captura o efeito estabilizante da flutuabilidade.
4. Comportamento Não Monotônico: Ao analisar os tempos de mistura dimensionais ($t_m$) em $Ri$ constante, o estudo observou uma dependência não monotônica em relação a $Re$ para números de Richardson mais altos ($Ri \geq 8$). Embora o aumento de $Re$ (e, portanto, a velocidade do agitador) inicialmente reduza o tempo de mistura ao aumentar a turbulência, aumentos adicionais levam a uma estratificação de densidade mais forte (devido ao acoplamento de parâmetros na configuração da simulação), o que suprime o transporte vertical. Esta interação resulta em um tempo de mistura mínimo em números de Reynolds intermediários.
5. Estrutura do Fluxo: A análise do campo de velocidade revelou que, em altos $Ri$, a fase mais leve pode ficar "selada" da entrada de potência do agitador, com estruturas verticais aparecendo principalmente no líquido mais denso. A análise de Densidade Espectral de Potência (PSD) indicou que, embora a flutuabilidade afete a organização espacial do fluxo, a dinâmica impulsionada pelo agitador ainda domina a cascata de energia.

Significado e Alegações
O artigo alega fornecer uma correlação preditiva compacta para tempos de mistura turbulenta em sistemas inicialmente estratificados por densidade dentro do regime totalmente turbulento. A lei de escalonamento derivada, baseada na influência combinada dos números de Potência, Froude e Richardson, oferece uma expressão analítica simples para engenheiros e especialistas em aplicações.

Os autores observam que, até onde lhes é conhecido, nenhuma correlação comparável para sistemas estratificados no regime totalmente turbulento foi previamente relatada. O significado deste trabalho reside em sua capacidade de otimizar o projeto de processos de mistura, reduzindo a dependência de testes reais caros. O modelo proposto permite a estimativa de tempos de mistura em sistemas onde os efeitos de flutuabilidade são significativos, capturando a competição complexa entre turbulência aprimorada e estratificação fortalecida. Os autores sugerem modestamente que trabalhos futuros envolverão testar esta correlação em diferentes geometrias de tanque e tipos de agitadores, incluindo agitadores angulados, para avaliar a robustez do escalonamento proposto.