Mixing induced by microswimmers as probed by mutual information

Technische Samenvatting: Door microswimmers geïnduceerde menging, onderzocht middels wederzijdse informatie

Probleemstelling
Vloeistofmenging op microschaal is een cruciaal proces in zowel biologische systemen (bijv. nutriëntenopname, biogene menging in oceanen) als kunstmatige toepassingen (bijv. microreactoren, lab-on-a-chip). In tegenstelling tot macroscopische menging worden stromingen op microschaal beheerst door lage Reynoldsgetallen, wat ze tijdsomkeerbaar maakt; roeren alleen is daarom onvoldoende, en menging is afhankelijk van het samenspel tussen advectie en moleculaire diffusie. Terwijl conventionele metrieken voor mengingsefficiëntie bestaan (bijv. $L_2$-normen, Lyapunov-exponenten), zijn deze vaak afhankelijk van specifieke initiële patronen of lokale rek (stretching). De auteurs beogen de menging die wordt geïnduceerd door microswimmers te onderzoeken met behulp van wederzijdse informatie (mutual information), een informatietheoretische maatstaf die eerder alleen gevalideerd is in vereenvoudigde modelsystemen (zoals Couette-stroming), om de toepasbaarheid ervan te bepalen voor complexe, actieve vloeistofstromingen gegenereerd door swimmers.

Methodologie
De studie maakt gebruik van een twee-dimensionaal (2D) squirmer-model om niet-deformeerbare, krachtvrije circulaire deeltjes te beschrijven die zwemmen in een Stokes-stroming.

• Stromingsmodel: Het stroomveld wordt gegenereerd door een tangentiële slipsnelheid op het oppervlak van de swimmer, gekarakteriseerd door modi $B_n$. De auteurs behouden de eerste twee modi: $B_1$ (bron-dipool, bepaalt de translationale snelheid) en $B_2$ (krachtdipool/stresslet, bepaalt het type swimmer: pusher, puller of stresslet). De squirmer-parameter wordt gedefinieerd als $\beta = B_2/B_1$.
• Opsluiting (Confinement): Simulaties worden uitgevoerd in een 2D vierkant domein met perfect-slip (stress-vrije) randvoorwaarden. De geconfineerde Green-functie wordt geconstrueerd met de methode van beelden en berekend in het Fourier-domein om numerieke stabiliteit te garanderen, inclusief spectrale regularisatie om eindige grootte van de swimmers te representeren.
• Kwantificering van menging: Mengingsefficiëntie wordt gekwantificeerd door de wederzijdse informatie $I[x_0; x_t]$ tussen de initiële positie van een deeltje $x_0$ en de uiteindelijke positie $x_t$ op tijdstip $t$. Dit wordt berekend als het verschil tussen de Shannon-entropie van de uiteindelijke positieverdeling en de conditionele entropie gegeven de initiële positie.
• Transportdynamica: De evolutie van de waarschijnlijkheidsdichtheidsfunctie (PDF) van tracer-deeltjes wordt numeriek opgelost met behulp van de advectie-diffusie (Fokker-Planck) vergelijking. De conditionele PDF wordt verkregen via een pseudo-spectrale methode.
• Parameters: De studie vergelijkt verschillende swimmer-configuraties (geordende vs. ongeordende posities/oriëntaties), verschillende swimmer-typen (stresslets, pushers, pullers, bron-dipolen) onder een vaste totale energie-dissipatie beperking, en systemen met en zonder hydrodynamische interacties (HIs).

Belangrijkste Resultaten

1. Impact van Ruimtelijke Organisatie:

   • Wanorde bevordert menging: Zeer symmetrische configuraties (bijv. geordende posities met identieke oriëntaties) onderdrukken grootschalig transport door stroomdomeinen te creëren die gescheiden worden door stagnatielijnen. Het introduceren van wanorde in de oriëntaties en posities van de swimmers doorbreekt deze symmetrieën, wat resulteert in meer onregelmatige snelheidsvelden die de stretching en dispersie van tracers versterken.
   • Aggregatie onderdrukt menging: Sterke aggregatie van swimmers in een gelokaliseerd gebied concentreert de geïnduceerde stroom, waardoor de buitenste regio's grotendeels onveranderd blijven. Naarmate de ruimtelijke verdeling van de swimmers breder wordt, neemt de mengingsefficiëntie toe.

2. Optimale Swimmer-type en Dissipatie:

   • Niet-monotone afhankelijkheid: Bij een vaste energie-dissipatiesnelheid hangt de mengingsefficiëntie niet-monotoon af van de squirmer-parameter $|\beta|$.
   • Trade-off: Er is een afweging tussen het genereren van sterke, lang reikende shear-stromingen (gedomineerd door de stresslet-component, $|\beta| \to \infty$) en het bereiken van ruimtelijke dispersie via de translationale beweging van de swimmer (gedomineerd door de bron-dipool component, $\beta=0$).
   • Optimaal eindige waarde: De optimale mengingsefficiëntie vindt plaats bij een eindige waarde van $|\beta|$, waarbij een balans wordt gevonden tussen de lang reikende stroom van stresslets en de translationale beweging van de swimmer. Pure stresslets ($\beta \to \infty$) zijn niet-motiel in afwezigheid van HIs en zijn daarom minder effectief voor globale menging dan swimmers met een eindige $|\beta|$.

3. Hydrodynamische Interacties (HIs):

   • Wanneer HIs worden meegerekend, presteren pushers ($\beta < 0$) beter dan pullers ($\beta > 0$). Dit wordt toegeschreven aan de neiging van pullers om aggregaten te vormen of een polaire uitlijning te vertonen, wat de mengingsefficiëntie vermindert.
   • De inclusie van HIs vlakt de afhankelijkheid van de efficiëntie van $|\beta|$ af en verschuift het optimale swimmer-type naar grotere waarden van $|\beta|$.

4. Mengregimes en Systeemgrootte:

   • Het verval van de wederzijdse informatie vindt plaats in drie duidelijke stadia:
     1. Diffusie-gedomineerd: Een initiële fase waarin menging pure diffusie-schaling volgt.
     2. Advectie-versterkt: Een intermediaire fase waarin advectie de menging versnelt; de vervalsnelheid is grotendeels onafhankelijk van de systeemgrootte voor een vaste swimmer-dichtheid.
     3. Relaxatiefase: De laatste fase die wordt gecontroleerd door de systeemgrootte, waarbij tracers interageren met de grenzen.
   • Grotere systemen behouden de wederzijdse informatie langer, en de relaxatietijd schaalt met de systeemgrootte, hoewel de afhankelijkheid zwakker is dan de $L^2$-schaling die verwacht wordt bij pure effectieve diffusie, vanwege de uitgebreide reikwijdte van de advectieve stromen.

Betekenis en Claims
Dit artikel toont aan dat wederzijdse informatie een robuuste, universele en aannamevrije maatstaf is voor het kwantificeren van mengingsefficiëntie in complexe actieve stromingen, en breidt de geldigheid ervan uit voorbij vereenvoudigde shear-stromingen. De studie onthult dat de door microswimmers geïnduceerde menging wordt beheerst door een fundamentele trade-off: het maximaliseren van menging vereist een balans tussen het genereren van sterke shear-stromingen (via dipolaire componenten) en de ruimtelijke dispersie die wordt geboden door de translationale beweging van de swimmer. Verder benadrukken de resultaten dat wanorde in de organisatie van swimmers en collectieve hydrodynamische interacties (die pushers bevoordelen boven pullers) cruciale factoren zijn bij het bepalen van de efficiëntie van biogene menging. Het kader biedt een kwantitatieve basis voor het begrijpen van transportfenomenen in actieve vloeistoffen, met potentiële relevantie voor ciliaire stromingen, microbiële suspensies en intracellulair transport.