Trkitlts Vletlen pakolsok Balla Pter Bevezet Alapfogalmak Szimulcik
Térkitöltés Véletlen pakolások Balla Péter
• Bevezető • Alapfogalmak • Szimulációk Súrlódó golyók pakolása, <z> viselkedése • Kísérleti eredmények, a Delaunay-dekompozíció ▪ Kísérleti technika ▪ A rendeződés geometriája ▪ A modell ▪ Eredmények, összehasonlítás a modellel • Összefoglalás
Bevezetés • • • Szemcsés anyagok (golyók), T=0 „Jamming” Nincs tökéletes térkitöltés Erős függés a minta preparációjától A mechanikai stabilitás feltétele, kontaktusszám (2 d, d+1) Izostatikusság, hiperstatikusság (JÖVŐ HÉT!)
Alapfogalmak • • • Kérdés: dobozba véletlenszerűen ejtett gömbök a térfogat hány százalékát foglalják el? Pakolási hányados: A zöldséges dilemmája: mi a lehetséges maximális térkitöltés? (Kepler-sejtés 1611; Hilbert 18. problémája 1900; Hales 1998)
Maximális térkitöltés
Szimulációk Súrlódó és súrlódásmentes pakolások • • MD szimuláció, N=20000 golyó d átmérő, m tömeg, periodikus határfeltétel, 20× 50 golyó 0, 02<φi<0, 3, innen a gravitáció hatása alatt esnek (forgás modellezve) Kezdeti kinetikus energia: 20 -100 mgd Összenyomhatóság keménység (mg/d egységekben): 2 E 5<kn<2 E 9 (kemény gömb) Súrlódás: 0<μ<10 Sebességek aránya egy centrális ütközés alatt: ε=0, 26; 0, 50; 0, 88
<z> és a pakolási hányados
P(z) és a keménygömb-átmenet
Függés a preparációtól
Kísérleti eredmények Röntgen-CT
Kísérleti eredmények • • Azonos méretű műanyaggolyók (2%-on belül) hengeralakú dobozban Pozíciómérés hibája: 0, 1%
A modell • • • Az elemi építőkövek feltérképezése A térfogatfluktuációk feltérképezése, térfogateloszlás meghatározása A kísérleti eredmények és a modell összevetése
A rendeződés geometriája • • • Cél: a rendeződés jellemzése egyszerű geometriai megfontolások alapján Kötött párok (küszöbtávolság: r=1. 0 d… 1. 11 d) Közös szomszéd (számuk n) Közös szomszédok eloszlása (rögzített r-re, kül. n-ek részvételi aránya) ▪ ha r<1, 118 d, akkor n max. 5 Lokális jellemzés ▪ ρ erősen korellál a fenti eloszlással Diéderszög
A rendeződés geometriája
A rendeződés geometriája: tetraéderes szerkezet
Következmények • • A rendeződés tetraéderek formájában történik Az eredmény lényegében független a küszöbtávolságtól (r=1. 0 d… 1. 11 d) A probléma kezelése: Delaunay-dekompozíció, tetraéderekre Minimális, véletlenszerűen elhelyezett tetraéderek, amik nem lógnak bele egymás körülírt gömbjébe (nincs szükség küszöbtávra!) Egyértelműen azonosít egy N-részecskés rendszert egy térkitöltő tetraéderrendszerrel (<f>: egy golyóra eső átlagos tetraéderszám, T tetraéderek száma): <f>=14… 15, 53 (szoros pakolás…„granuláris gáz”) mechanikailag stabil, egyforma gömbökből álló rendszerben, a gravitáció hatása alatt <f>~14, 5 (az egyensúlyi kontaktusszámmal ekvivalens)
A tetraéderek elrendeződése • Lokálisan lehetséges legsűrűbb elrendeződés, tökéletes d-élhosszú tetraéderek, ez a dekmpozícióban szereplő minimálisan elérhető elemi térfogat (Roger, geometriai alsó határ): • Térkitöltés (az elemi tetraéderek kitöltik a térfogatot): • A fenti feltételekkel, és az elemi térfogatok között térfogatcserét megengedve, a termodinamikai határesetben:
A tetraéderek elrendeződése • A mechanikai stabilitás: • A geometriai kényszerek: • A térkitöltés: • Nincs szabad paraméter!
Referenciák • Szimulációk Geometry of frictionless and frictional sphere packings Leonardo E. Silbert et al. (2006) • Kísérlet és Delaunay-dekompozíciós model Volume fluctuations and geometrical constraints in granular packs Tomaso Aste (2006)
Összefoglalás • • Szimulációk (a kontaktusszám, izostatikusság, súrlódás): izostatikusság: csak súrlódásmentes esetben, súrlódással mindig hiperstatikus szerkezet (bonyolultabb modell szükségeltetik) igaz a keménygömb-határesetben is <z> folytonosan változik a súrlódás bekapcsolásával a részletek erős előélet-függést mutatnak A kísérlet és a modell (a térkitöltés geometriája): a rendeződés lokálisan tetraéderek formájában történik a kísérletekkel jó egyezésben leírható a Delaunay-dekompozícióval, amely eleget tesz az alábbi (természetes) feltételeknek: mechanikai stabilitás geometriai kényszerek térkitöltés
Köszönöm a figyelmet!
- Slides: 23