A Variational Formulation for Deformable Particle Simulations and its Level Set Discrete Element Method Implementation

Dit artikel presenteert een variatieformulering voor de Discrete Element Method (DEM) die, via een energiegebaseerde benadering en Level Set-technieken, de robuustheid en schaalbaarheid van klassieke DEM combineert met een efficiënte modellering van deformatie in deeltjessystemen.

De kern

Stel je voor dat je een bak vol met knikkers hebt. Als je die bak schudt, bewegen de knikkers als harde, onbuigzame balletjes. Ze botsen tegen elkaar, rollen en stuiteren, maar ze vervormen nooit. Ze blijven perfect rond. Dit is hoe computers vaak simuleren hoe zand, grind of poeders zich gedragen: als een hoopje harde balletjes.

Maar in het echte leven zijn deeltjes niet altijd perfect hard. Een korrel zand kan een beetje platdrukken, een stukje schuim kan vervormen, of een biologisch weefsel kan buigen. De oude methoden om dit te simuleren waren ofwel te simpel (alsof je deeltjes als onbreekbare stenen ziet) ofwel veel te zwaar voor de computer (alsof je elke korrel zand uit miljoenen atomen opbouwt, wat duizenden jaren rekentijd kost).

De oplossing uit dit papier: De "Slimme Vouwknikkers"

De auteurs van dit artikel, Thomas Henzel en Konstantinos Karapiperis, hebben een nieuwe manier bedacht om deze deeltjes te simuleren. Ze noemen het een Deformable DEM (Discrete Element Method).

Hier is hoe het werkt, vertaald naar alledaagse taal:

1. De "Vouwbare" Deeltjes

In plaats van een deeltje te zien als een stijf balletje, zien ze het als een slimme, vouwbare knikker.

• Oude manier: Een deeltje heeft alleen een positie (waar het is) en een draaiing (hoe het ligt).
• Nieuwe manier: Het deeltje heeft diezelfde posities, maar krijgt ook een paar "magische knoppen" (wiskundige variabelen) waarmee je het kunt laten buigen, rekken of plompen.

Stel je voor dat je een deeltje hebt dat je kunt veranderen van een perfect bolletje naar een platte koekje. In de computercode gebeurt dit niet door het deeltje in duizenden stukjes te hakken, maar door een paar vooraf ingestelde "buigpatronen" te gebruiken. Het is alsof je een origami-vogel hebt die je niet opnieuw vouwt, maar die je met één handbeweging kunt laten uitrekken of samendrukken volgens een vast patroon.

2. De "Level Set" (Het Onzichtbare Net)

Hoe weet de computer nu hoe het deeltje eruitziet als het vervormt?
Ze gebruiken een techniek die ze Level Set noemen.

• De analogie: Stel je voor dat je een deeltje in een bak met water legt. Het wateroppervlak is je "Level Set". Als het deeltje vervormt, verandert de vorm van het wateroppervlak.
• In de computer wordt dit gedaan met een digitaal raster (een soort 3D-scherm). Het deeltje is het gebied waar de "waarde" negatief is, en de lucht is waar de waarde positief is. Als het deeltje vervormt, schuift dit onzichtbare net mee.
• Het slimme aan hun methode is dat ze dit net niet elke seconde opnieuw hoeven te berekenen (wat heel langzaam is). Ze gebruiken een slimme truc: ze "trekken" het net terug naar de oorspronkelijke vorm, kijken hoe het eruitzag, en passen dat toe op de nieuwe vorm. Dit is als het updaten van een Google Maps-kaart: je hoeft niet de hele wereld opnieuw te tekenen, je past alleen de straten aan die je hebt veranderd.

3. Waarom is dit zo cool?

De auteurs hebben bewezen dat hun methode twee grote problemen oplost:

1. Snelheid: Het is net zo snel als het simuleren van harde balletjes. Je kunt dus miljoenen deeltjes simuleren zonder dat je computer ontploft.
2. Realisme: Het deeltjes kunnen nu echt vervormen. Als je een hoopje zand samendrukt, worden de korrels plat. In de oude methoden bleven ze rond en gleedden ze over elkaar. In hun nieuwe methode plompen ze, wat zorgt voor een veel realistischer beeld van hoe materialen zich gedragen.

De "Wiskundige Zwaartekracht" (Variational Formulation)

Het papier zit vol met zware wiskunde over "variaties" en "energie". In het kort:
Ze hebben een nieuwe wet bedacht die zegt: "De natuur zoekt altijd de weg van de minste weerstand."
Ze hebben deze wet zo aangepast dat hij niet alleen geldt voor het bewegen van deeltjes, maar ook voor het vervormen ervan. Hierdoor kunnen ze precies berekenen hoeveel energie er in een vervormde korrel zit en hoe die korrel terugveert. Het is alsof je een elastiekje hebt: je weet precies hoe hard je moet trekken om het te rekken, en hoe hard het terugveert.

Samenvattend

Dit artikel introduceert een nieuwe manier om te simuleren hoe zand, poeder of biologische weefsels zich gedragen.

• Vroeger: Deeltjes waren als harde stenen (te simpel) of als duizenden atomen (te traag).
• Nu: Deeltjes zijn als slimme, vouwbare knikkers die snel en realistisch vervormen.

Het is een beetje alsof je van een poppenkast met stijve poppen overstapt naar een poppenkast met poppen van zachte, rekbaar rubber. Je kunt ze nog steeds gooien en laten stuiteren, maar nu kun je ze ook een beetje platdrukken, waardoor ze veel natuurlijker reageren op de wereld om hen heen. Dit helpt ingenieurs en wetenschappers om betere medicijnen, sterkere bouwmaterialen en betere voorspellingen voor aardverschuivingen te maken.

Technische samenvatting

Titel: Een Variationale Formulering voor Simulaties van Vervormbare Deeltjes en de Implementatie via de Level Set Discrete Element Method (LS-DEM)

Auteurs: Thomas Henzel en Konstantinos Karapiperis (EPFL, Zwitserland)

---

1. Het Probleem

Discrete deeltjessystemen (zoals granulaire materialen, bodems, biologische weefsels en farmaceutische poeders) vertonen mechanisch gedrag dat wordt bepaald door interacties tussen deeltjes én de vervorming van de deeltjes zelf.

• Beperkingen van bestaande methoden:
  • Classieke DEM (Discrete Element Method): Gaat uit van starre deeltjes. Hoewel contactwetten lokale vervorming (zoals Hertziaanse contacten) benaderen, kunnen ze geen bulk-vervorming (globale vormverandering van het deeltje) modelleren.
  • FEM (Finite Element Method): Kan vervorming nauwkeurig oplossen, maar is computatief te duur voor grote systemen met duizenden deeltjes vanwege de noodzaak van mesh-hergeneratie en complexe contactproblemen.
  • Bestaande vervormbare DEM-aanpakken:
    • FEM-DEM koppeling: Te duur voor schaalbaarheid.
    • Aggregaten van kleine starre deeltjes: Vereisen veel sub-deeltjes voor realisme, wat de rekentijd explodeert.
    • Vereenvoudigde modalen: Vaak beperkt tot specifieke geometrieën, gebrek aan fysische consistentie of thermodynamische grondslagen.

Er is een behoefte aan een methode die bulk-vervorming mogelijk maakt met een rekentijd vergelijkbaar met starre DEM, terwijl de geometrische flexibiliteit en schaalbaarheid behouden blijven.

---

2. Methodologie

De auteurs presenteren een unificerend raamwerk dat de klassieke DEM uitbreidt met een reduced-order beschrijving van elastische deeltjesvervorming, gebaseerd op een variational formulering en geïmplementeerd binnen de Level Set DEM (LS-DEM).

A. Variationale Formulering (Theorie)

Het model is gebaseerd op het Lagrange-d'Alembert-principe, dat de dynamica van het systeem beschrijft door de variatie van de actie te balanceren met het virtuele werk van niet-conservatieve krachten (zoals wrijving).

• Toestandsvariabelen: Naast de gebruikelijke translatie ($X$) en rotatie ($\Theta$) worden modale vervormingscoördinaten ($\varepsilon_\alpha$) geïntroduceerd.
• Kinematica: De vervorming van een deeltje wordt uitgedrukt als een lineaire combinatie van vooraf gedefinieerde vormfuncties (modes) $\Phi_\alpha$:
  $$u_\varepsilon(x, t) = \varepsilon_\alpha(t) \Phi_\alpha(x)$$
  Hierbij zijn de vormfuncties vast in de referentieconfiguratie, maar evolueren de amplitudes $\varepsilon_\alpha$ in de tijd.
• Energie: Het systeem omvat kinetische energie (translatie, rotatie, vervorming) en elastische potentiële energie. De kinetische energie van vervorming leidt tot een generaliseerde massamatrix ($M_{\alpha\beta}$), en de potentiële energie tot een generaliseerde stijfheid ($K_{\alpha\beta}$) of een niet-lineair potentiaal $U_\varepsilon$.
• Besturingsvergelijkingen: Door de variatie te nemen, worden de bewegingsvergelijkingen afgeleid. Deze vormen een gekoppeld systeem van gewone differentiaalvergelijkingen (ODE's) voor translatie, rotatie en de modale coördinaten:
  $$M_\varepsilon \ddot{\varepsilon} + K_\varepsilon \varepsilon = F_\varepsilon$$
  Waarbij $F_\varepsilon$ de generaliseerde krachten zijn die voortvloeien uit contactinteracties.

B. Numerieke Implementatie (LS-DEM)

De methode is geïmplementeerd in het LS-DEM-raamwerk, waar deeltjesgeometrie wordt gedefinieerd door een level set-functie op een rooster en oppervlaktenodes.

• Level Set Evolutie: De grootste uitdaging is het bijwerken van de level set-functie bij vervorming zonder de nauwkeurigheid te verliezen. De auteurs vergelijken verschillende strategieën:
  1. Herberekening: Accuraat maar duur.
  2. Advectie-oplossing: Vereist het oplossen van een PDE, wat complex is.
  3. Semi-Lagrangiaans (Gekozen methode): Gebruikt de bekende verplaatsingsvelden om punten terug te projecteren naar de onvervormde configuratie ($x_{ref} = x - u$) en interpolatie toe te passen. Dit is zeer efficiënt en vereist geen globale PDE-oplossing.
  4. Interpolatie van voorgecomputeerde levels: Zeer snel, maar vereist veel geheugen.
• Krachtenberekening: Contactkrachten worden berekend op de nodes van het ene deeltje tegen de level set van het andere. Deze krachten worden vervolgens geprojecteerd op de modale coördinaten om de generaliseerde krachten $F_\varepsilon$ te bepalen.

C. Modale Extractie

Voor complexe deeltjes zijn analytische vormenfuncties vaak niet beschikbaar. De auteurs stellen een strategie voor om numerieke modes te extraheren uit high-fidelity FEM-simulaties (bijv. via Proper Orthogonal Decomposition/SVD) of experimenten. Deze modes kunnen zowel lineair als niet-lineair zijn (bijv. door veranderende contactcondities).

---

3. Belangrijkste Resultaten

De methode werd gevalideerd aan de hand van verschillende voorbeelden:

1. Analytische Bending Mode: Een slanke balk onder driedpunt-bending. De simulatie toonde uitstekende overeenkomst met de analytische oplossing voor de kracht-verplaatsingsrespons en de vervormingsvorm.
2. Numerieke Compaction Mode: Een bolvormig deeltje onder hydrostatische compressie (geïnduceerd door niet-lineaire contactcondities).
   • Een vervormingsmode werd geëxtraheerd uit een FEM-simulatie.
   • De vervormbare LS-DEM kon de contactvlak-vervlakking en volumetrische verdichting nauwkeurig nabootsen.
   • De kracht-verplaatsingsrespons stemde perfect overeen met de FEM-resultaten.
3. Systeem-schaal Simulatie: Een pakket van 27 bollen onder isotrope compressie.
   • De simulatie resulteerde in een kracht-verplaatsingsrespons die perfect overeenkwam met een semi-analytische afleiding gebaseerd op de energiebalans van een enkel deeltje.
   • Dit bewijst dat de methode schaalbaar is en correcte macroscopisch gedrag voorspelt.

---

4. Bijdragen en Significantie

De belangrijkste bijdragen van dit werk zijn:

• Unificerend Variational Raamwerk: Voor het eerst wordt de vervormbaarheid van deeltjes in DEM systematisch afgeleid uit het Lagrange-d'Alembert-principe, wat een fysisch consistente basis biedt voor energie- en impulsbehoud.
• Efficiëntie en Schaalbaarheid: De methode behoudt de rekentijd van de orde van grootte van starre DEM, terwijl het bulk-vervorming mogelijk maakt. Dit maakt simulaties van grote systemen (duizenden deeltjes) met vervorming haalbaar.
• Geometrische Flexibiliteit: Door gebruik te maken van Level Sets en modale decompositie, is de methode niet beperkt tot specifieke vormen (zoals bollen) en kan het complexe topologieën en willekeurige vervormingsmodi hanteren.
• Niet-lineariteit: Het raamwerk kan niet-lineaire effecten (zowel materiaal- als contact-gedreven) accommoderen door het gebruik van niet-quadratische energiefunctionaliteiten.
• Praktische Toepasbaarheid: De methode biedt een brug tussen de nauwkeurigheid van FEM en de snelheid van DEM, wat essentieel is voor het modelleren van materialen zoals nacre, biomassa, en geavanceerde granulaire metamaterialen.

Conclusie:
De voorgestelde Deformable LS-DEM is een krachtige, uitbreidbare en computatie-efficiënte tool voor het modelleren van de mechanica van deeltjessystemen waarbij de vervorming van individuele korrels een cruciale rol speelt in het macroscopische gedrag. Het overbrugt de kloof tussen starre contactmechanica en continuümmechanica.