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

Technische Samenvatting: De rol van wrijving in grofkorrelige klei ontrafeld

Probleemstelling
Het gedrag van montmorilloniet (MMT) klei wordt beheerst door een complexe hiërarchie van interacties, variërend van nanoschaal hydratatiekrachten tot mesoschaal jamming en kinematische vergrendeling. Hoewel Coarse-Grained Molecular Dynamics (CGMD) een vitale methode is geworden om de kloof tussen atomistische simulaties en macroscopische geotechnische observaties te overbruggen, lijden traditionele CGMD-frameworks aan een kritieke beperking: de historische verwaarlozing van interpartikel wrijving en gelokaliseerde viskeuze demping. Conventionele modellen vertrouwen op effectieve potentialen afgeleid van evenwichtstoestanden, die inherent het vrije-energieprofiel afvlakken en niet-conserverende krachten weglaten. Deze omissie leidt tot onfysisch gedrag, zoals het onvermogen om realistische jamming-transities te vangen en het falen om de afschuifsterkte en structurele integriteit te reproduceren die worden waargenomen in macroscopische kleiassemblages. De mate waarin dit gebrek aan inzicht de simulatienauwkeurigheid beïnvloedt en de specifieke rol van dissipatieve parameters bij het bepalen van de kleifabriek blijven onopgeloste vragen.

Methodologie
Om deze hiaten aan te pakken, stellen de auteurs een nieuw hybride CGMD-framework voor dat langafstands moleculaire interacties expliciet koppelt aan granulaire contactmechanica.

• Interactiemodel: Het framework maakt gebruik van een duaal componentenmodel. Langafstandsinteracties (van der Waals en elektrostatische dubbellaag-repulsie) worden gemodelleerd met een Buckingham-potentiaal. Kortafstandsinteracties, gedomineerd door hydratatiekrachten en sterische contacten, worden gemodelleerd met een Hertziaans contactkrachtveld met visco-elastische demping. Dit maakt gecontroleerde partikeloverlap toe om de compressibiliteit van de hydratatieschil na te bootsen, terwijl het tegelijkertijd tangentiële wrijving en viskeuze demping incorporeert.
• Systeemconstructie: De MMT-plaatjes worden gemodelleerd als hexagonale rigide lichamen bestaande uit 271 overlappende sferen (diameter 120 Å). Het systeem bestaat uit 1000 van dergelijke plaatjes binnen een kubische simulatiebox.
• Simulatieprotocol: Simulaties werden uitgevoerd met LAMMPS. De studie maakte gebruik van een basismodel dat gevalideerd is via incrementele isotrope compressie (1 tot 100 atm). Parametrische analyses werden uitgevoerd om de effecten van normale demping ($\eta_n$), tangentiële dempingsschaling ($x_{\gamma,t}$) en wrijvingscoëfficiënten ($\mu$) te isoleren. Cruciaal was de vergelijking tussen thermostatted (NVT) simulaties en unthermostatted (NVE) protocollen om onderscheid te maken tussen thermodynamisch evenwicht en puur mechanische quenching. Vervolgens werden uniaxiale compressietests uitgevoerd om de afschuifsterkte te evalueren.

Belangrijkste Resultaten

1. Modelvalidatie: Het hybride basismodel reproduceerde succesvol de macroscopische compressibiliteit van klei. De gesimuleerde compressie-index ($C_c$) stemde overeen met experimentele gegevens voor bentoniet en smectiet, waarbij een bimodal respons werd getoond: een hoge compressibiliteit ($C_c \approx 2,17$) bij lage- tot gemiddelde druk en een stijvere respons ($C_c \approx 0,39$) bij hoge druk (>1000 kPa).
2. Rol van Visco-elastische Demping: Parametrische studies toonden aan dat dempingscoëfficiënten de structurele evolutie beheersen.
   • Normale Demping: Verhoogde normale demping ($\eta_n$) onderdrukte densificatie en uitlijning aanzienlijk. Het verhogen van $\eta_n$ van 0,2 naar 2,0 verdubbelde bijna de porositeit (van 1,03 naar 1,97) en verminderde de oriëntatie-orde-parameter ($S$) van 0,86 naar 0,04, waardoor de plaatjes effectief gevangen werden in gedesorganiseerde configuraties met een hoog porievolume.
   • Tangentiële Demping: Hoge tangentiële demping ($x_{\gamma,t} = 10$) remde de rotatie en het glijden van de deeltjes, waardoor het systeem niet de dichte, geordende evenwichtstoestand kon bereiken die werd waargenomen bij lagere dempingsniveaus ($x_{\gamma,t} = 1$ of $5$).
3. Thermische Fluctuaties vs. Kinetische Trapping: Vergelijkingen tussen thermostatted en unthermostatted simulaties toonden aan dat het uitsluitend vertrouwen op dempingskrachten zonder een thermostaat fungeert als een agressief quenching-proces. Dit drijft het systeem naar een staat met bijna nul kinetische energie, waardoor de plaatjes gevangen raken in lokale, metastabiele energiedips en de zelfassemblage naar thermodynamisch stabiele, vrije-energie-geminimaliseerde toestanden wordt voorkomen.
4. Noodzaak van Wrijving: Uniaxiale compressietests onderstreepten de cruciale rol van interpartikel wrijving. Hoewel wrijvende ($\mu=0,1$) en wrijvingsloze ($\mu=0$) modellen vergelijkbare initiële elasticiteitsmoduli vertoonden, divergeerden ze significant voorbij 10% rek. Het wrijvingsloze model vertoonde onfysisch fluïde-achtig vloeien en rek-verweking met een lage piekspanning (~2 MPa), terwijl het wrijvende model aanzienlijk hogere belastingen kon dragen dank aan de vergrendeling van glijdende interfaces.

Betekenis en Claims
Het artikel claimt een baanbrekend hybride framework te presenteren dat granulaire contactmechanica verenigt met bottom-up moleculaire potentialen. De primaire betekenis ligt in het oplossen van het langdurige probleem van onfysisch fluïde-achtig gedrag in kleisimulaties door expliciet interpartikel wrijving en demping op te nemen. De studie stelt vast dat:

• Macroscopische mechanische integriteit in CG-kleimodellen direct voortkomt uit de synergie tussen moleculaire energielandschappen en mesoscopische energie-dissipatie.
• Precieze controle over temperatuur- en dempingsparameters essentieel is om kunstmatige kinetische trapping te voorkomen en structurele realisme te waarborgen.
• Expliciete wrijving onmisbaar is voor het vastleggen van de afschuifsterkte en structurele stabiliteit van kleiassemblages, wat essentiële richtlijnen biedt voor toekomstige multi-schaal simulaties die streven naar het overbruggen van de kloof tussen atomistische beschrijvers en de complexe realiteit van geotechnische systemen.