Unveiling the Role of Friction in Coarse-Grained Clay: A Hybrid Framework Integrating Long-Range Interactions and Granular Contact Mechanics

Resumo Técnico: Revelando o Papel do Atrito em Argila de Granulação Grossa

Definição do Problema
O comportamento da argila montmorilonita (MMT) é governado por uma complexa hierarquia de interações, que variam desde forças de hidratação em escala nanométrica até o travamento (jamming) e o intertravamento cinemático em escala mesoscópica. Embora a Dinâmica Molecular de Granulação Grossa (CGMD) tenha emergido como uma ferramenta vital para preencher a lacia entre simulações atomísticas e observações geotécnicas macroscópicas, os modelos tradicionais de CGMD sofrem de uma limitação crítica: a negligência histórica do atrito interparticular e do amortecimento viscoso localizado. Os modelos convencionais baseiam-se em potenciais efetivos derivados de estados de equilíbrio, que inerentemente suavizam o panorama de energia livre e omitem forças não conservativas. Essa omissão leva a comportamentos não físicos, como a incapacidade de capturar transições de jamming realistas e a falha em reproduzir a resistência ao cisalhamento e a integridade estrutural observadas em assembleias de argila macroscópicas. A extensão em que essa negligência influencia a precisão da simulação e o papel específico dos parâmetros dissipativos no controle da estrutura da argila permanecem questões não resolvidas.

Metodologia
Para abordar essas lacunas, os autores propõem um novo arcabouço híbrido de CGMD que acopla explicitamente interações moleculares de longo alcance com a mecânica de contato granular.

• Modelo de Interação: O arcabouço emprega uma abordagem de dois componentes. As interações de longo alcance (forças de van der Waals e repulsão de dupla camada eletrostática) são modeladas usando um potencial de Buckingham. As interações de curto alcance, dominadas por forças de hidratação e contato estérico, são modeladas usando um campo de força de contato Hertziano com amortecimento viscoelástico. Isso permite um controle sobre a sobreposição de partículas para mimetizar a compressibilidade da camada de hidratação, incorporando simultaneamente atrito tangencial e amortecimento viscoso.
• Construção do Sistema: Os plaquetas de MMT são modeladas como corpos rígidos hexagonais compostos por 271 esferas sobrepostas (diâmetro de 120 Å). O sistema consiste em 1000 dessas plaquetas dentro de uma caixa de simulação cúbica.
• Protocolo de Simulação: As simulações foram conduzidas no LAMMPS. O estudo utilizou um modelo de linha de base validado através de compressão isotrópica incremental (1 a 100 atm). Análises paramétricas foram realizadas para isolar os efeitos do amortecimento normal ($\eta_n$), da escala de amortecimento tangencial ($x_{\gamma,t}$) e dos coeficientes de atrito ($\mu$). Crucialmente, o estudo comparou simulações termostatizadas (NVT) contra protocolos não termostatizados (NVE) para distinguir entre o equilíbrio termodinâmico e o resfriamento (quenching) puramente mecânico. Testes de compressão uniaxial foram subsequentemente realizados para avaliar a resistência ao cisalhamento.

Resultos Principais

1. Validação do Modelo: O modelo híbrido de linha de base reproduziu com sucesso a compressibilidade macroscópica da argila. O índice de compressão simulado ($C_c$) alinhou-se com dados experimentais de bentonita e esmectita, mostrando uma resposta bimodal: uma alta compressibilidade ($C_c \approx 2,17$) em pressões baixas a médias e uma resposta rígida ($C_c \approx 0,39$) em altas pressões (>1000 kPa).
2. Papel do Amortecimento Viscoelástico: Estudos paramétricos revelaram que os coeficientes de amortecimento governam a evolução estrutural.
   • Amortecimento Normal: O aumento do amortecimento normal ($\eta_n$) suprimiu significativamente a densificação e o alinhamento. Aumentar $\eta_n$ de 0,2 para 2,0 quase dobrou a razão de vazios (de 1,03 para 1,97) e reduziu o parâmetro de ordem orientacional ($S$) de 0,86 para 0,04, efetivamente prendendo as plaquetas em configurações desorganizadas de alta razão de vazios.
   • Amortecimento Tangencial: Um alto amortecimento tangencial ($x_{\gamma,t} = 10$) inibiu a rotação e o deslizamento das partículas, impedindo o sistema de atingir o estado de equilíbrio denso e ordenado observado em níveis de amortecimento mais baixos ($x_{\gamma,t} = 1$ ou $5$).
3. Flutuações Térmicas vs. Aprisionamento Cinético: Comparações entre simulações termostatizadas e não termostatizadas demonstraram que depender apenas de forças de amortecimento sem um termostato atua como um processo de resfriamento (quenching) agressivo. Isso conduz o sistema a um estado de energia cinética próximo de zero, aprisionando as plaquetas em mínimos de energia locais e metaestáveis, impedindo a auto-organização em estados termodinamicamente estáveis e de energia livre minimizada.
4. Necessidade de Atrito: Testes de compressão uniaxial destacaram o papel crítico do atrito interparticular. Enquanto modelos com atrito ($\mu=0,1$) e sem atrito ($\mu=0$) exibiram módulos elásticos iniciais semelhantes, eles divergiram significativamente além de 10% de deformação. O modelo sem atrito exibiu um comportamento de escoamento fluido-like não físico e amolecimento por deformação (strain-softening) com uma baixa tensão de pico (~2 MPa), enquanto o modelo com atrito sustentou cargas significativamente maiores devido ao travamento das interfaces de deslizamento.

Significância e Alegações
O artigo alega apresentar um arcabouço híbrido pioneiro que unifica a mecânica de contato granular com potenciais moleculares bottom-up. Sua principal significância reside na resolução do problema de longo prazo do comportamento fluido-like não físico em simulações de argila, ao incorporar explicitamente o atrito interparticular e o amortecimento. O estudo estabelece que:

• A integridade mecânica macroscópica em modelos de argila de granulação grossa deriva diretamente da sinergia entre os panoramas de energia molecular e a dissipação de energia mesoscópica.
• O controle preciso dos parâmetros de temperatura e amortecimento é essencial para evitar o aprisionamento cinético artificial e garantir o realismo estrutural.
• O atrito explícito é indispensável para capturar a resistência ao cisalhamento e a estabilidade estrutural de assembleias de argila, oferecendo diretrizes essenciais para futuras simulações multiescala que visam unir descritores atomísticos à realidade complexa dos sistemas geotécnicos.