🔬 condensed matter

Criticality of the viscous to inertial transition near jamming in non-Brownian suspensions

Utilizando simulações de Método de Elementos Discretos, este estudo demonstra que a transição de espessamento por cisalhamento viscoso-inercial em suspensões não-Brownianas densas e sem atrito próximas ao empacotamento é governada por uma escala de comprimento microestrutural divergente que faz com que a taxa de cisalhamento característica da transição desapareça à medida que a fração volumétrica se aproxima do limite de empacotamento.

Autores originais: Nishanth Murugan, Donald Koch, Sarah Hormozi

Publicado 2026-07-16

📖 5 min de leitura🧠 Leitura aprofundada

CC BY 4.0

Autores originais: Nishanth Murugan, Donald Koch, Sarah Hormozi

Artigo original sob licença CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta é uma explicação gerada por IA do artigo abaixo. Não foi escrita nem endossada pelos autores. Para precisão técnica, consulte o artigo original. Ler aviso legal completo

Imagine um mundo onde as regras de como as coisas fluem mudam dependendo da velocidade com que você as empurra. Este é o fascinante reino da dinâmica de fluidos, focando especificamente em "suspensões" — misturas onde minúsculas partículas sólidas estão flutuando em um líquido, como areia na água ou gotas de chocolate no leite. Normalmente, pensamos nos fluidos como sendo ou escorregadios (como a água) ou espessos (como o mel), mas essas misturas repletas de partículas podem ser complicadas. Às vezes, se você mexer lentamente, elas agem como um líquido normal. Mas se você mexer rápido, elas podem subitamente tornar-se rígidas, quase como um sólido. Esse fenômeno é chamado de "engrossamento por cisalhamento" (shear thickening). Os cientistas estão obcecados em entender exatamente quando e por que essa mudança acontece, porque isso afeta tudo, desde o processamento de cimento e cerâmica até a previsão de eventos naturais perigosos, como deslizamentos de terra e avalanches. O grande mistério tem sido: qual é o exato "limite de velocidade" onde o fluido decide deixar de ser escorregadio e começar a agir como um sólido? É um número fixo ou muda dependendo de quão crowded (congestionadas) as partículas estão?

Neste estudo, os pesquisadores Nishanth Murugan, Donald Koch e Sarah Hormozi, da Universidade de Cornell, usaram simulações computacionais poderosas para resolver uma peça deste quebra-cabeça. Eles focaram em um tipo específico de mistura: suspensões densas de partículas que não se grudam (atrito zero/frictionless) e são grandes demais para serem agitadas pelo calor (não-Brownianas). Pense nelas como uma pista de dança lotada onde os dançarinos estão usando roupas escorregadias e não podem bater uns nos outros. A equipe queria encontrar o momento preciso em que a pista de dança muda de um deslize suave para uma multidão caótica e agitada.

Suas simulações revelaram que não existe um único limite de velocidade universal para essa transição. Em vez disso, o "limite de velocidade" depende inteiramente de quão compacta a pista de dança está. À medida que as partículas se aproximam do ponto onde estão completamente travadas e não conseguem se mover (um estado chamado de "jamming" ou entupimento), a velocidade necessária para desencadear o comportamento de endurecimento cai para quase zero. Em outras palavras, quanto mais congestionado o sistema fica, mais fácil é transformar o fluido em um sólido.

Os autores descobriram que isso acontece devido a uma estrutura oculta e crescente dentro da mistura. À medida que as partículas se aglomeram, elas começam a formar grandes clusters temporários que agem como uma única unidade massiva. O tamanho desses clusters cresce cada vez mais à medida que o sistema se aproxima do jamming. Os pesquisadores descobriram que a transição de "escorregadio" para "rígido" é, na verdade, desencadeada pela inércia (o momento) desses gigantescos clusters, e não apenas pelas partículas individuais. Eles provaram isso ao mostrar que o tamanho desses clusters e a velocidade com que o fluido endurece seguem o mesmo padrão matemático.

Curiosamente, o artigo argumenta contra a ideia de que essa transição é controlada pelo tamanho de uma única partícula ou por uma força fixa. Em sistemas onde as partículas possuem atrito (como areia áspera), a transição ocorre em uma velocidade diferente e constante. Mas para essas partículas suaves e escorregadias, a transição é "crítica", o que significa que está profundamente ligada ao congestionamento do sistema. O estudo sugere que, conforme você se aproxima do ponto de jamming, uma "escala de comprimento" (uma medida de quão grandes esses clusters ficam) cresce infinitamente, tornando o sistema incrivelmente sensível até ao menor empurrão.

Os pesquisadores utilizaram um método chamado Método de Elementos Discretos (DEM), que é como um videogame super preciso onde eles calculam as forças entre milhões de partículas. Eles simularam versões 2D e 3D desses sistemas lotados. Seus resultados mostraram que o "número de Stokes" (uma forma sofisticada de medir o equilíbrio entre a viscosidade do fluido e o momento das partículas) cai para zero conforme a fração volumétrica se aproxima do limite de jamming. Eles descobriram que o tamanho dos clusters cresce de acordo com uma lei de potência específica, com um expoente de cerca de 0,79 em 2D e 0,84 em 3D.

Ao testar como os resultados mudavam quando tornavam a caixa de simulação menor ou maior, eles confirmaram que esses gigantescos clusters realmente existem e são a chave para a transição. Se a caixa de simulação fosse menor que o cluster, os resultados mudavam, provando que o cluster estava tentando "atravessar" as paredes. Isso confirmou que a transição é um evento não-local, o que significa que o que acontece em uma parte do fluido depende do que está acontecendo longe dali, conectado por essas estruturas massivas e invisíveis.

Em resumo, o artigo sugere que, para partículas suaves e congestionadas, a mudança de líquido para sólido não é sobre os dançarinos individuais; é sobre a formação de uma multidão massiva e invisível que cresce até tomar conta de toda a sala. Essa percepção ajuda os cientistas a entender que, próximo ao ponto de jamming, as regras de fluxo tornam-se críticas e dependem dessas estruturas crescentes, abrindo caminho para melhores teorias sobre como esses materiais complexos se comportam.

Resumo Técnico: Criticidade da Transição Viscosa para Inercial próximo ao Entupimento em Suspensões Não-Brownianas

Definição do Problema
A reologia de suspensões densas não-brownianas é crítica para compreender fenômenos naturais (ex: deslizamentos de terra, lamaçais) e processos industriais (ex: lama de cimento, cerâmicas). Uma característica central desses sistemas é a transição de um regime viscoso, caracterizado por uma reologia Newtoniana ( $\Sigma \propto \dot{\gamma}$ ), para um regime inercial, caracterizado por uma reologia Bagnoldiana ( $\Sigma \propto \dot{\gamma}^2$ ). Esta transição é tradicionalmente atribuída ao aumento do domínio da inércia das partículas sobre as forças hidrodinâmicas.

No entanto, o mecanismo físico que desencadeia essa transição e a taxa de cisalhamento característica na qual ela ocorre permanecem em debate. A literatura experimental e numérica apresenta resultados conflitantes em relação ao número de Stokes crítico ( $St_{v \to i}$ ), uma medida adimensional da razão entre as tensões inerciais e viscosas. Em suspensões sem atrito (ex: esferas de poliestireno com repulsão eletrostática), $St_{v \to i}$ foi observado variando com a fração volumétrica ( $\phi$ ) e aproximando-se de zero conforme o sistema se aproxima da fração volumétrica de entupimento (jamming) ( $\phi_m$ ). Por outro lado, em suspensões com atrito (ex: partículas de PMMA rugosas), $St_{v \to i}$ parece ser independente de $\phi$ e permanece finito. A falta de consenso sobre o mecanismo físico, particularmente o papel da microestrutura em sistemas sem atrito, motiva este estudo.

Metodologia
Os autores utilizam simulações de Método de Elementos Discretos (DEM) para investigar suspensões densas, não-brownianas e sem atrito através de uma ampla gama de taxas de cisalhamento e frações volumétricas.

Modelo de Força: A dinâmica é governada pelas equações de movimento de Newton, contabilizando duas forças primárias: forças de lubrificação hidrodinâmica ( $F^L_{ij}$ ) decorrentes do movimento relativo no intervalo interparticular, e forças de repulsão de contato ( $F^R_{ij}$ ).
Aproximação Sem Atrito: Para modelar a barreira de repulsão eletrostática encontrada em experimentos (que impede o contato direto), a força de repulsão é modelada como uma mola de Hooke ( $F^R_{ij} = K\delta n$ ) que é ativada quando a separação superficial $h$ cai abaixo de um limiar mínimo $h_{min}$ .
Configuração da Simulação: As simulações são conduzidas em domínios 2D e 3D usando condições de contorno de Lees-Edwards para simular o fluxo de cisalhamento. O sistema utiliza uma distribuição de partículas polidispersa (16%) para evitar a cristalização. As forças de lubrificação são calculadas assumindo uma separação mínima de $h_{min}$ quando as partículas tentam se aproximar mais, efetivamente restringindo o movimento normal.
Parâmetros: O estudo explora o espaço de parâmetros definido pela fração volumétrica ( $\phi$ ), o número de Stokes ( $St = 4\rho a^2 \dot{\gamma} / \eta_f$ ) e a razão de rigidez adimensional ( $\Gamma$ ). As simulações são realizadas em um regime onde $\Gamma \gg 1$ , garantindo que a barreira de repulsão domine sobre as forças de lubrificação em pequenas separações.

Principais Contribuições e Resultados

Criticidade da Transição: As simulações confirmam que, para suspensões sem atrito, o número de Stokes característico para a transição para o espessamento por cisalhamento inercial ( $St_{v \to i}$ ) não é constante. Em vez disso, ele diminui conforme a fração volumétrica aumenta, aproximando-se de zero conforme $\phi \to \phi_m$ . A relação segue uma lei de potência: $St_{v \to i} \sim (\phi_m - \phi)^\beta$ .
Estrutura de Escalonamento (Scaling Framework): Para explicar esse comportamento, os autores propõem uma estrutura de escalonamento baseada na hipótese de que uma escala de comprimento microestrutural diverge conforme o sistema se aproxima do entupimento. Eles introduzem um número de Stokes modificado, $St_{Cluster} = B St (\phi_m - \phi)^{-\beta}$ , que escala os efeitos inerciais em relação a um tamanho de cluster crescente $L_c \sim (\phi_m - \phi)^{-\beta/2}$ .
Colapso de Dados: Utilizando esta estrutura de escalonamento, os dados reológicos (viscosidade de cisalhamento relativa $\eta_r$ $η_{r}$ vs. $St$) para várias frações volumétricas colapsam em uma única curva mestra em simulações 2D e 3D.
- Em 2D, o expoente é $\beta \approx 0.79$ .
- Em 3D, o expoente é $\beta \approx 0.84$ .
- Estes valores alinham-se com observações experimentais que sugerem $St_{v \to i} \sim (\phi_m - \phi)^{-1}$ .
Dependência do Tamanho do Sistema: Para validar a existência da escala de comprimento divergente proposta, os autores realizaram testes de dependência do tamanho do sistema. Observaram que, quando o tamanho da caixa de simulação é menor que o tamanho característico do cluster, a viscosidade aumenta devido ao cluster interagir consigo mesmo através de fronteiras periódicas. Ao ajustar os dados de viscosidade a uma lei de potência dependente do tamanho do sistema, eles extraíram com sucesso uma escala de comprimento de cluster $L_c$ que diverge conforme $\phi \to \phi_m$ com um expoente consistente com o escalonamento reológico ( $L_c \sim (\phi_m - \phi)^{-\beta/2}$ ).

Significância e Alegações
O artigo afirma que a transição do espessamento por cisalhamento viscoso para o inercial em suspensões sem atrito é um fenômeno crítico impulsionado por uma escala de comprimento microestrutural divergente (clusters) conforme o sistema se aproxima do entupimento.

Mecanismo: A criticidade da transição é atribuída ao crescimento destes clusters, que altera a escala na qual os efeitos inerciais se tornam dominantes. À medida que o tamanho do cluster diverge, a taxa de cisalhamento característica necessária para desencadear efeitos inerciais desaparece.
Reologia Não-Local: A existência de uma escala de comprimento divergente implica que a reologia destas suspensões é não-local; a transmissão de tensão depende da formação e destruição destes clusters de grande escala.
Contraste com Sistemas com Atrito: Os autores observam que este comportamento crítico é específico para sistemas sem atrito. Em suspensões com atrito, onde os contatos entre partículas são governados pela rugosidade superficial em vez de repulsão de longo alcance, a transição parece ser independente da fração volumétrica, sugerindo que o mecanismo microestrutural difere fundamentalmente.
Direção Futura: O trabalho motiva o desenvolvimento de teorias não-locais que deem conta da transmissão de tensão através destes clusters crescentes. Também destaca a necessidade de novos estudos sobre as leis de escala para tamanhos de cluster em suspensões com atrito para estabelecer uma teoria unificada através de diferentes regimes de fluxo.

Mais como este