A boundary integral approach to the eigenvalue problem for the anisotropic bidomain operator with perfect contact conditions

Titel

Een randintegraalbenadering voor het eigenwaardeprobleem van de anisotrope bidomeinoperator met perfecte contactvoorwaarden.

1. Probleemstelling

Het artikel behandelt het eigenwaardeprobleem geassocieerd met de bidomeinoperator in een anisotroop, heterogeen domein dat de hartspier (myocardium) modelleert. Het domein is onderverdeeld in drie subregio's:

• $\Omega_I$: Het rechterventrikel.
• $\Omega_S$: Het septum (tussenschot).
• $\Omega_D$: Het linker ventrikel.

Kernuitdagingen:

• Anisotropie: De elektrische geleidbaarheid is richtingsafhankelijk (geleidend langs vezels vs. dwars op vezels). De vezelrichtingen veranderen per subdomein, wat leidt tot discontinuïteiten in de geleidbaarheidstensors op de interfaces.
• Transmissievoorwaarden: Er worden "perfecte contactvoorwaarden" aangenomen, wat betekent dat zowel het elektrische potentiaal als de stroomflux continu zijn over de interfaces, ondanks de discontinuïteit in de geleidbaarheid.
• Doel: Het berekenen van de eigenwaarden en eigenfuncties van dit systeem. Deze eigenfuncties zijn cruciaal voor theoretische analyses (zoals het bewijs van periodieke oplossingen) en numerieke benaderingen van de volledige bidomeinvergelijkingen.

2. Methodologie

De auteurs gebruiken potentiaaltheorie om het oorspronkelijke probleem van partiële differentiaalvergelijkingen (PDE's) te reduceren tot een systeem van randintegraalvergelijkingen (RIV's).

Stappen in de aanpak:

1. Formulering: Het probleem wordt gedefinieerd als een lineair systeem van elliptische vergelijkingen (Helmholtz-type) met anisotrope diffusie in elk subdomein, gekoppeld via transmissievoorwaarden op de interfaces ($\Gamma_{IS}$ en $\Gamma_{SD}$) en Neumann-voorwaarden (nul flux) op de externe rand.
2. Fundamentele Oplossing: Er wordt een expliciete uitdrukking afgeleid voor de fundamentele oplossing van de anisotrope Helmholtz-operator. Deze oplossing wordt uitgedrukt in termen van Besselfuncties van de tweede soort ($N_0$ en $N_1$), rekening houdend met de anisotrope schaling door de geleidbaarheidstensors.
3. Directe Formulering: De oplossing in elk domein wordt voorgesteld als een lineaire combinatie van enkelvoudige laag-potentialen (single-layer) en dubbelvoudige laag-potentialen (double-layer). De dichtheden van deze potentialen corresponderen met de onbekende waarden van het potentiaal en de stroomflux op de randen.
4. Grensvergelijkingen: Door de limiet te nemen wanneer een punt de rand nadert, worden de bekende sprongrelaties (jump relations) voor deze potentialen toegepast. Dit leidt tot een gekoppeld systeem van Fredholm-randintegraalvergelijkingen van de tweede soort.
5. Discretisatie: Het continue systeem wordt gediskretiseerd door de randen en interfaces in segmenten te verdelen. Speciale kwadratuurformules worden gebruikt om de zwakke singulariteiten (logaritmisch gedrag) in de kernen te behandelen.

3. Belangrijkste Bijdragen

• Dimensionaliteitsreductie: In plaats van het hele volume te diskretiseren (zoals bij de eindige-elementenmethode), reduceert deze methode het probleem tot de interfaces en de externe rand. Dit verlaagt de dimensionality en vermindert de rekenlast aanzienlijk.
• Expliciete Kernen: De auteurs leiden expliciete formules af voor de kernen van de integraalvergelijkingen, gebaseerd op Besselfuncties, specifiek voor de anisotrope bidomeincontext met verschillende vezeloriëntaties in de ventrikels en het septum.
• Nieuwe Hypothese over het Septum: Het artikel introduceert en onderzoekt het idee dat het septum fungeert als een functionele barrière voor golfvoortplanting, niet vanwege anatomische barrières, maar vanwege discontinuïteiten in de geleidbaarheidstensors op de interfaces, zelfs onder perfecte contactvoorwaarden.
• Numeriek Systeem: Er wordt een robuust numeriek schema ontwikkeld dat het eigenwaardeprobleem oplost door de rangverlies van het gediskretiseerde lineaire systeem te analyseren (via singuliere waarden of determinanten).

4. Resultaten

• Numerieke Validatie: De methode werd getest met specifieke parameters voor hartweefsel (geleidendheidswaarden voor intracellulair en extracellulair medium, vezeloriëntaties).
• Gedrag van Potentiaal: De numerieke resultaten tonen dat de elektrische potentiaal over de interfaces een sprong vertoont, ondanks de continuïteit van de stroomflux. Dit bevestigt dat de discontinuïteit in de geleidbaarheid de golfvoortplanting beïnvloedt.
• Schaal: De berekende potentiaalwaarden zijn zeer klein ($10^{-46}$ tot $10^{-47}$), wat consistent is met het passieve regime zonder externe stimulatie en de specifieke schaling van het geometrische model.
• Stabiliteit: Het systeem van integraalvergelijkingen (7 vergelijkingen voor 7 onbekenden) bleek goed oplosbaar en leverde gebonden, coherente oplossingen op die de verwachte kwalitatieve eigenschappen van het model weerspiegelen.

5. Betekenis en Conclusie

Dit werk biedt een efficiënt alternatief voor volumebaseerde methoden (zoals FEM) voor het analyseren van anisotrope hartproblemen met complexe interfacecondities.

• Efficiëntie: Door de dimensie te reduceren, wordt de complexiteit van het oplossen van grote lineaire systemen verminderd.
• Toepasbaarheid: De methode is bijzonder geschikt voor het bestuderen van de invloed van vezeloriëntatie en geleidbaarheidsdiscontinuïteiten op de elektrische activiteit van het hart.
• Toekomstperspectief: De auteurs wijzen erop dat deze benadering de basis legt voor het bestuderen van "imperfecte" transmissievoorwaarden (waarbij de stroom ook discontinu is), wat een volgende stap zou zijn in het modelleren van pathologische harttoestanden zoals bundeltakblokkades.

Samenvattend demonstreert dit artikel hoe potentiaaltheorie kan worden gebruikt om complexe anisotrope biofysische problemen in de cardiale elektrofyziologie te vereenvoudigen tot een numeriek hanteerbaar randintegraalprobleem.