Fluid-Solid Transition in Dense Granular Systems

Piaskowa gora nie jest fazą materii – dopóki jej nie uformujesz

Wyloju piasek z dłoni przypomina ciecz. Zatrzymując strumień i pozwalając mu się ussett, te same ziarnka teraz wytrzymują wagę jak ciało stałe. Nic chemicznego się nie zmieniło – nie utworzyły się żadnych wiązań między ziarnami, nie pojawił się żaden kryształowy wzór. Zmiana dotyczyła wyłącznie geometrii: ziarna skończyły się na przestawienie względem siebie. Jest to zamykanie, a jest to jeden z czystych przykładów systemu, który nabywa sztywności dzięki ograniczeniom wynikającym z pakowania, a nie ze względu na jakąkolwiek siłę przyciągającą.

Najprostszym laboratorium do jego badania nie jest piasek, ale 2D układ monodispersyjnych twardych dysków bezślizgowych, uwięzionych w pudełku. Powoli zmniejszaj pudełko (co odpowiada dodawaniu więcej dysków przy ustalonej wielkości pudełka) i śledź stopień wypełnienia φ – czyli procent powierzchni pudełka faktycznie pokrytej dyskami.

Krytyczna gęstość pakowania φ_J

Poniżej krytycznej gęstości, dyski mają wystarczająco dużo luzu, aby przesuwać się obok siebie pod wpływem jakiegokolwiek obciążenia – przepływy pakowania zachodzą, jednak bardzo powoli i nie generują żadnego stałego naprężenia odśrodkowego. Powyżej tej wartości, każdy dysk jest geometrycznie zablokowany przez sąsiednie dyski: nie ma ścieżki zmiany ułożenia, która nie wymagałaby nakładających się dysków, więc pakowanie opiera się na tarciu jak ciało stałe. Dla bezślizgłych, jednorodnych dysków w dwóch wymiarach ten próg znajduje się przy gęstości pakowania losowego:

φ_J ≈ 0,84 (2D, jednorodny, bezślizgały dyski – punkt J) φ < φ_J → pakowanie jest NIEZABLOKOWANE: przepływy, zerowy moduł sprężystości, zerowe ciśnienie φ > φ_J → pakowanie jest ZABLOKOWANE: skończony moduł sprężystości, skończone ciśnienie Nie jest to zaokrąglenie jakiejś gładkiej krzywej – dokładnie przy φ_J zarówno moduł sprężystości, jak i moduł ścisku przechodzą z dokładnie zerowego, a tuż nad nim rosną jako funkcja potęgowa nadmiarowej gęstości, φ − φ_J. Odpowiednia wartość dla 3D bezślizgłych sfer to φ_J ≈ 0,64, wartość, którą każdy student fizyki słyszał nazywaną „losowym pakowaniem w bliskie sąsiedztwo”. Obie wartości pochodzą od tej samej podstawowej idei: jammingu, a nie krystalizacji, uporządkowanej zgrupowania twardych cząstek.

φ_J ≈ 0.84   (2D, monodisperse, frictionless discs — point J)

φ < φ_J   →  packing is UNJAMMED: flows, zero shear modulus, zero pressure
φ > φ_J   →  packing is JAMMED: finite shear modulus, finite pressure

Izostatyczność: dlaczego przejście jest tak ostre

Ostrość punktu uwięzienia ma precyzyjne wyjaśnienie mechaniczne, oparte na argumentach Jamesa Clerk’a Maxwella dotyczącą sztywności: pakowanie N bezślizgłych dysków posiada 2N stopni swobody (w 2D) i potrzebuje co najmniej 2N sił kontaktowych, aby być mechanicznie sztywnym i statycznie określonym. Wprost na punkcie uwięzienia średnia liczba kontaktów na cząstkę, czyli liczba kontaktowa z, dokładnie mieści się w tym izostatycznym wartości:

z_iso = 2·d (d = wymiar przestrzenny: z_iso = 4 w 2D, z_iso = 6 w 3D) z < z_iso → niedostatecznie określone: istnieją "tryby plastyczne", struktura może się swobodnie deformować z ≥ z_iso → nadmiernie określone: brak możliwości swobodnej deformacji, pakowanie jest sztywne Wprost na φ_J pakowanie ma dokładnie z_iso kontaktów na cząstkę średnio, żadnego więcej, ani jednego mniej – ledwie wystarczająco, aby zamocować każdą cząsteczkę w miejscu. Dlatego uwięzienie wygląda jak prawdziwy punkt krytyczny, z nadmiarem liczby kontaktów Δz i nadmiarem modułu, które skalują się jako czyste prawa potęg φ - φ_J, oraz używa się do tego samego języka matematycznego, co w przypadku perkalizacji sieci sprężyn w artykule na tej stronie.

z_iso = 2·d          (d = spatial dimension: z_iso = 4 in 2D, z_iso = 6 in 3D)

z < z_iso  →  under-constrained: "floppy modes" exist, structure can deform for free
z ≥ z_iso  →  over-constrained: no free deformation, packing is rigid

Łączy oparte na siłowych łańcuchach: sztywne, ale nierównomierne

Po uwięzieniu, pakowanie nie rozkłada naskalanej obciążenia w taki sposób, jak robiłoby to ciągłe, sprężyste ciało stałe. Siły kontaktu skupiają się w cienkich, włóknistych łańcuchach siłowych, które wiją się przez pakowanie wzdłuż kontaktów, które zdarzają się być zorientowane w sposób korzystny, podczas gdy cząstki tuż obok tych łańcuchów przenoszą prawie żadnego obciążenia. Eksperymenty z dyskami fotoelastycznymi — gdzie naprężone dyski polaryzują światło i świecą pod skrzyżowanymi filtrami — sprawiają, że te łańcuchy są widoczne bezpośrednio, a to wyjaśnia, dlaczego materiały granulowane mogą wytrzymywać ogromne lokalnie obciążenia (jeden odcisk buta na piasku) podczas pozostawania globalnie miękkie i łatwo zaburzone w innych miejscach.

Faza utulzenia

Andrea Liu i Sidney Nagel zaproponowali w 1998 roku, że utulzenie powinno być postrzegane jako jeden narożnik znacznie szerszej, jednorodnej fazowej diagramy, z trzema osiami: gęstość (1/φ), temperatura i przyłożone naprężenie ślizgowe. Wzdłuż każdej osi, zwiększanie gęstości, obniżanie temperatury lub zmniejszanie naprężenia śliskiego prowadzi nieuporządkowany układ do tego samego sztywnego stanu utulzenia – dlatego też piany, emulsje, koloidalne szkła i opakowania granulowe, mimo że są to bardzo różne materiały fizyczne, dzielą ten sam krytyczny skalowanie w pobliżu ich odpowiednich punktów utulzenia. Granularna materia o temperaturze zerowej jest po prostu osią, na której fluktuacje termiczne są nieistotne, a układ geometryczny decyduje o wyniku.

Symulacja kompresji

Minimalna symulacja integruje siły nakładania się dysków miękkich (liniową lub Hertzowską odpychanie, które włącza się tylko wtedy, gdy dwa dyski się nakładają) z tłumioną dynamiką molekularną i powoli zmniejsza wymiar o partycji lub rośnie promień dysków, mierząc przy tym ciśnienie całkowite. Ciśnienie utrzymuje się na poziomie zerowym maszyny poniżej φ_J i gwałtownie wzrasta natychmiast, gdy dyski są zmuszone do nakładania się powyżej tego poziomu – to pojawienie się ciśnienia jest najczystszym sposobem na lokalizację φ_J numerycznie bez potrzeby liczenia kontaktów ręcznie.

// pairwise soft-disc repulsion, only active on overlap
for each pair (i, j):
  dr = distance(i, j);  overlap = (r_i + r_j) - dr;
  if (overlap > 0) {
    f = k * overlap;              // linear spring, k = stiffness
    applyForce(i, j, f, direction = (pos_i - pos_j) / dr);
  }
// pressure ≈ sum of overlap forces × contact distances / (2 × box area)

Często zadawane pytania

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Granular Jamming i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.