A numerical study on the coefficient of restitution of wet collisions

Dit onderzoek gebruikt SPH-simulaties om te bepalen dat de terugkaatsingscoëfficiënt bij natte botsingen afhangt van het Stokes-getal en een dimensieloze filmdikte, waarbij twee verschillende regimes met specifieke machtwet-exponenten worden geïdentificeerd.

De kern

De Kunst van het Droge en Natte Springen: Een Simpele Uitleg

Stel je voor dat je een balletje laat vallen. Als je het op een droge vloer laat vallen, springt het een beetje terug. Dat is makkelijk te begrijpen. Maar wat gebeurt er als je datzelfde balletje laat vallen op een dun laagje water of olie? Dan is het verhaal veel ingewikkelder.

Deze wetenschappelijke studie, gedaan door onderzoekers uit India en Engeland, probeert precies dat uit te leggen: hoeveel energie een balletje verliest als het in een nat laagje landt.

Hier is de uitleg, vertaald naar alledaags taalgebruik:

1. Het Probleem: Waarom is nat vallen lastig?

Wanneer een balletje op een droge ondergrond landt, is het verlies aan energie vrijwel altijd hetzelfde. Maar als er een laagje vloeistof tussen zit, verandert alles.

• De Smeertheorie (Lubricatie): Bij lage snelheden werkt het als een smerende laag. Het balletje moet de vloeistof "wegduwen" voordat het de grond raakt. Dat kost energie, net als wanneer je met je hand door honing probeert te duwen.
• De Inertie (De "Spat"): Bij hogere snelheden is het niet meer alleen maar smeren. De vloeistof moet heel snel opzij. Hierdoor ontstaat er een drukgolf onder het balletje, alsof je met een snelheid van 100 km/u op een wateroppervlak landt.

De onderzoekers wilden weten: Is er één simpele regel die voorspelt hoe hoog het balletje terugsprikt, of hangt het van meer dingen af?

2. De Methode: Een Digitale Zandbak

In plaats van duizenden balletjes fysiek te laten vallen (wat veel tijd en geld kost), gebruikten de onderzoekers een computerprogramma (genaamd SPH).

• De Analogie: Denk aan dit programma als een gigantische, digitale zandbak. In plaats van zandkorrels, gebruiken ze duizenden kleine "deeltjes" om water of olie te simuleren.
• De Slimme Truc: Om de computer niet te laten crashen, hebben ze het balletje niet volledig in detail nagebootst. Ze simuleerden alleen het onderste deel van het balletje dat met het water in contact komt, terwijl de rest van het balletje "onzichtbaar" maar wel zwaar bleef. Dit bespaart enorm veel rekenkracht, net als het bouwen van een maquette van een auto in plaats van de hele fabriek.

3. De Ontdekking: Twee Werelden

De onderzoekers lieten balletjes van verschillende maten vallen in vloeistof van verschillende diktes. Ze ontdekten dat er twee verschillende werelden zijn, afhankelijk van hoe dik het laagje vloeistof is ten opzichte van het balletje.

Stel je twee situaties voor:

• Wereld 1: Het "Dikke" Laagje (Kleine balletjes of dik water)
  Hier is het laagje vloeistof relatief dik. Het balletje moet veel water verplaatsen.

  • Wat gebeurt er? Het balletje maakt grote, chaotische bewegingen in het water (zoals kleine draaikolken of wervelingen). Het is alsof je probeert te springen in een badkuip die halfvol is; je maakt veel waterbewegingen die energie kosten.
  • De Regel: De snelheid waarmee het balletje valt (de "Stokes-getal") maakt hier weinig uit. Het belangrijkste is hoe dik het laagje is. De energie gaat vooral verloren door die wervelingen.

• Wereld 2: Het "Dunne" Laagje (Grote balletjes of heel dun water)
  Hier is het laagje vloeistof heel dun, bijna als een film.

  • Wat gebeurt er? Het balletje duwt het water heel snel weg, maar er is geen ruimte voor grote wervelingen. Het is meer als het wegduwen van een dunne laagje olie op een tafel.
  • De Regel: Hier telt de snelheid van het balletje wel degelijk mee. Hoe harder het valt, hoe meer energie er verloren gaat door wrijving.

4. De Grote Regel (De "Wiskundige Formule")

Voorheen dachten wetenschappers dat er maar één regel was die alles beschreef. Deze studie bewijst dat dat niet klopt.
Ze hebben een nieuwe formule bedacht die twee verschillende "modi" heeft:

1. Als het laagje dik is ten opzichte van het balletje, geldt de ene formule.
2. Als het laagje dun is, geldt de andere formule.

Het is alsof je twee verschillende regels hebt voor het rijden:

• Op een sneeuwweg (dik laagje) telt vooral hoe zwaar je auto is en hoe je remt (wervelingen).
• Op een natte asfaltweg (dun laagje) telt vooral hoe hard je rijdt (wrijving).

5. Waarom is dit belangrijk?

Deze kennis is niet alleen leuk voor de wetenschap, maar ook heel nuttig voor de praktijk:

• Industrie: Bedrijven die poeders verwerken (zoals in de farmacie of bouw) moeten weten hoe korrels met elkaar botsen als ze nat zijn.
• Natuur: Het helpt bij het begrijpen van modderstromen of hoe zand wordt getransporteerd in rivieren.
• Computersimulaties: Als je een computerprogramma maakt dat duizenden balletjes simuleert (bijvoorbeeld voor een video-game of een fabrieksontwerp), kun je nu veel nauwkeuriger voorspellen hoe ze zich gedragen zonder dat de computer uren moet rekenen.

Kortom:
De onderzoekers hebben ontdekt dat "nat vallen" niet één ding is. Het hangt af van de verhouding tussen de dikte van het water en de grootte van het balletje. Soms is het chaos (wervelingen), en soms is het gewoon wrijving. Met hun nieuwe regels kunnen we dit nu veel beter voorspellen.

Technische samenvatting

Titel: Een numeriek onderzoek naar de terugkaatsingscoëfficiënt van natte botsingen

Auteurs: Abhishek Kumar Singha, Christopher Robert Kit Windows-Yule, Prapanch Nair
Affiliaties: IIT Delhi (India) en Universiteit van Birmingham (VK)

---

1. Probleemstelling

Het modelleren van granulaire stromingen (zoals puinstromen, coatingprocessen en sedimenttransport) vereist een nauwkeurige schatting van energie-dissipatie tijdens botsingen. Voor droge botsingen wordt de terugkaatsingscoëfficiënt (COR, Coefficient of Restitution) vaak als een materiaalconstante behandeld die slechts zwak afhankelijk is van de impactsnelheid.

Bij natte botsingen (waarbij een interstitiële vloeistoffilm aanwezig is) verandert dit beeld fundamenteel. De aanwezigheid van de vloeistof introduceert extra dissipatiemechanismen zoals viskeuze stroming, inertiale drukopbouw en vervorming van het vrij oppervlak.

• Huidige beperking: De meeste macroscopische modellen gaan er ten onrechte van uit dat de natte COR uitsluitend wordt bepaald door het Stokes-getal (St).
• Het gat in de kennis: Experimentele studies (zoals Gollwitzer et al.) hebben aangetoond dat de COR niet alleen van St afhangt, maar ook van de verhouding tussen de vloeistoffilm-dikte en de deeltjesdiameter ($\delta/D$). Eerdere studies hebben de onafhankelijke invloed van deze geometrische verhouding niet volledig onderzocht, vooral niet in het regime van matige tot hoge Weber-getallen (waar traagheidseffecten domineren boven oppervlaktespanning).

2. Methodologie

De auteurs gebruiken Incompressible Smoothed Particle Hydrodynamics (ISPH) om de botsing van een stijf bolletje met een dunne vloeistoffilm te simuleren.

• Simulatiekader: Een mesh-vrij, updated-Lagrangiaans raamwerk dat de drukgedreven dissipatie en vloeistof-vaste stof interacties oplost.
• Vloeistofmodel: Gebruik van de Cummins-Rudman projectie-methode voor incompressibiliteit en het Cleary-Monaghan viskeuze model. Oppervlaktespanning wordt verwaarloosd omdat het onderzoek zich richt op het regime van matige tot hoge Weber-getallen ($O(10) - O(10^2)$).
• Stijf lichaam dynamiek: Het vaste deeltje wordt gemodelleerd als een cluster van Lagrange-deeltjes, maar om rekentijd te besparen, wordt alleen een afgeknotte bolvormige schil (dikte 4 deeltjesbreedtes) gebruikt in het gebied waar interactie met de vloeistof optreedt. De totale massa, zwaartepunt en traagheidsmoment worden echter berekend voor de volledige bol.
• Botsingsmodel:
  • De botsing met de bodem wordt gedetecteerd en de snelheid wordt aangepast op basis van een droge COR ($e_r = 0.98$).
  • De tijdstap wordt adaptief verkleind wanneer het deeltje dicht bij de wand komt om numerieke stabiliteit te garanderen.
• Validatie: Het model is gevalideerd tegen experimentele data van Gollwitzer et al. (botsing van glazen bollen op siliconenolie). De simulaties tonen een uitstekende overeenkomst (RMSE = 0.00524) bij een resolutie van ongeveer 178.000 vloeistofdeeltjes.

3. Belangrijkste Bijdragen

1. Ontwikkeling van een schaalwet: De auteurs identificeren dat de natte COR niet alleen een functie is van het Stokes-getal, maar ook van een dimensieloze filmdikte $\gamma = \delta/D$.
2. Identificatie van twee regimes: Door systematisch variaties in deeltjesdiameter, impactsnelheid, viscositeit en filmdikte, worden twee distincte schaalregimes ontdekt met verschillende machts-wet exponenten.
3. Efficiënt algoritme: Een nieuwe methode voor het modelleren van stijve lichamen in SPH die de rekentijd aanzienlijk verlaagt zonder nauwkeurigheid te verliezen.
4. Fysische interpretatie: Een kwalitatieve verklaring voor de twee regimes gebaseerd op de evolutie van kinetische energie en secundaire stromingsverschijnselen (zoals wervels).

4. Resultaten

De studie toont aan dat de COR ($e$) een machtswet volgt van de vorm $e \sim St^a \gamma^b$. De data vertoont twee duidelijk gescheiden regimes:

• Regime R1 (Lager $\gamma$, grotere deeltjes):

  • De COR hangt sterk af van zowel het Stokes-getal als de dimensieloze filmdikte.
  • Fitted vergelijking: $e \approx 0.17 \cdot St^{0.17} \cdot \gamma^{-0.16}$.
  • Fysiek: De dissipatie wordt voornamelijk bepaald door viskeuze weerstand en traagheidseffecten in de primaire stroming. De kinetische energie in de vloeistof dissipeert monotoon na contact met de bodem.

• Regime R2 (Hoger $\gamma$, kleinere deeltjes):

  • De afhankelijkheid van het Stokes-getal is verwaarloosbaar zwak; de COR wordt gedomineerd door de geometrische verhouding $\gamma$.
  • Fitted vergelijking: $e \approx 0.4 \cdot St^{0.01} \cdot \gamma^{-0.22}$.
  • Fysiek: Secundaire stromingsverschijnselen (zoals wervelvorming en golfbewegingen) spelen een grotere rol in de dissipatie. Dit blijkt uit een tweede piek in de kinetische energie-evolutie van de vloeistof, wat wijst op complexere stromingspatronen.

De overgang tussen deze twee regimes wordt gedefinieerd door een lijn in de log-log ruimte van $St$ en $\gamma$.

5. Betekenis en Conclusie

Dit onderzoek weerlegt de aanname dat het Stokes-getal voldoende is om natte botsingen te karakteriseren. De introductie van de dimensieloze filmdikte $\gamma$ is essentieel voor een nauwkeurige beschrijving, vooral wanneer de verhouding tussen filmdikte en deeltjesgrootte varieert.

• Toepassing: De ontwikkelde schaalwetten kunnen direct worden toegepast in Discrete Element Method (DEM) simulaties en andere statistische benaderingen van granulaire stromingen om energie-verliezen bij natte botsingen realistischer te modelleren.
• Fysisch inzicht: Het onderzoek benadrukt dat de overgang van viskeus-gedreven naar traagheid-gedreven dissipatie (en de rol van secundaire stromingen) cruciaal is voor het begrijpen van het terugkaatsingsgedrag van natte deeltjes.

Samenvattend biedt deze studie een robuust numeriek raamwerk en een nieuwe, tweedelige schaalwet die de complexiteit van natte botsingen in het matige tot hoge Weber-getal regime beter beschrijft dan eerdere modellen.