Oscillatory behaviour of the RBF-FD approximation accuracy under increasing stencil size

Titel: Oscillerend gedrag van de RBF-FD benaderingsnauwkeurigheid bij toenemende stencil-grootte

Auteurs: Andrej Kolar-Požun, Mitja Jančič, Miha Rot en Gregor Kosec (Jožef Stefan Instituut & Universiteit van Ljubljana)

1. Probleemstelling

Bij het oplossen van partiële differentiaalvergelijkingen (PDV's) op verspreide knooppunten met de Radial Basis Function generated Finite Difference (RBF-FD) methode, is de keuze van de stencil-grootte ($n$, het aantal dichtstbijzijnde buren) een kritieke parameter.

• Context: De auteurs focussen op Polyharmonic Splines (PHS) met monomiale augmentatie. Een groot voordeel van PHS is dat deze geen vormparameter (shape parameter) vereist, wat de optimalisatie vereenvoudigt.
• De Uitdaging: Hoewel de nauwkeurigheidsorde in PHS RBF-FD wordt bepaald door de graad van de augmentatiepolynomen, heeft de keuze van de stencil-grootte een aanzienlijke invloed op de daadwerkelijke nauwkeurigheid.
• Observatie: De auteurs stellen vast dat de fout in de oplossing niet monotoon afneemt bij het vergroten van de stencil, maar een oscillerend gedrag vertoont. Dit betekent dat er specifieke stencil-groottes zijn waarbij de fout minimaal is, gevolgd door perioden van hogere fouten. Het doel is om deze minima te begrijpen en te voorspellen.

2. Methodologie

De studie is uitgevoerd met de volgende opzet:

• Testprobleem: De Poisson-vergelijking ($\nabla^2 u = f$) op een cirkelvormig domein.
• Analytische oplossing: $u(x, y) = \sin(\pi x)\sin(\pi y)$, zodat de exacte fout kan worden berekend.
• Discretisatie: Het domein is gediscretiseerd met een afstand $h = 0.01$. De Laplace-operator is benaderd met RBF-FD, gebruikmakend van radiale kubieken ($\phi(r) = r^3$) als PHS, geaugmenteerd met monomen tot graad $m=3$.
• Variabele: De stencil-grootte $n$ (aantal buren) werd systematisch gevarieerd.
• Foutmeting: Er werd gekeken naar zowel de maximale absolute fout ($e_{max}$) als de gemiddelde absolute fout ($e_{avg}$) voor zowel de oplossing als de Laplace-benadering.
• Controle-experimenten: Om uit te sluiten dat de oscillaties veroorzaakt werden door numerieke artefacten, werd getest met verschillende lineaire oplosmethoden (SparseLU, BiCGSTAB) en met verschillende discretisatiedichtheden. Ook werd onderzocht of rand-effecten de oorzaak waren door de stencil-grootte bij de randen vast te houden.

3. Belangrijkste Resultaten

• Oscillaties: De fout ($e_{max}$ en $e_{avg}$) toont duidelijke lokale minima en maxima bij specifieke stencil-groottes (bijv. minima bij $n=28$ en $n=46$). Het verschil in fout tussen een lokaal maximum en minimum kan meer dan een orde van grootte bedragen.
• Oorzaakanalyse:
  • De oscillaties zijn niet veroorzaakt door de keuze van de lineaire oplosmethode.
  • Ze zijn niet het gevolg van een te grove discretisatie (de vorm van de oscillatie blijft behouden bij verfijning, alleen de schaal verandert).
  • Ze zijn niet primair veroorzaakt door de problematische rand-stencils.
• Ruimtelijke Foutverdeling: De sleutel tot het inzicht ligt in de tekens van de puntsgewijze fout ($e^{\pm}$).
  • Bij fout-maxima heeft de fout over het hele domein hetzelfde teken (bijv. overal positief of overal negatief).
  • Bij fout-minima wisselt het teken van de fout binnen het domein; er zijn gebieden met positieve en negatieve fouten die elkaar gedeeltelijk opheffen.
• Nieuwe Kwantiteit ($\delta N^{\pm}$): Gebaseerd op deze observatie introduceerden de auteurs een nieuwe maatstaf:
  $$\delta N^{\pm}_{poiss} = \frac{1}{N_{int}} \left( |\{x_i : e > 0\}| - |\{x_i : e < 0\}| \right)$$
  Deze grootheid representeert het gemiddelde teken van de fout.
  • Bij fout-maxima is $|\delta N^{\pm}| \approx 1$ (dominant één teken).
  • Bij fout-minima nadert $|\delta N^{\pm}|$ de waarde 0 (balans tussen positieve en negatieve fouten).
  • De grootheid $\delta N^{\pm}_{lap}$ (gebaseerd op de Laplace-fout) gedraagt zich als een "gegladde" versie en heeft een nulpunt nabij de optimale stencil-groottes.

4. Bijdragen en Significatie

• Inzicht in Mechanisme: Het artikel legt een directe link tussen de oscillatie in de globale nauwkeurigheid en de ruimtelijke verdeling van de fouttekens. Het toont aan dat optimale nauwkeurigheid optreedt wanneer de fouten zich intern in het domein opheffen (wisselende tekens) in plaats van zich op te bouwen (zelfde teken).
• Praktische Indicator: De geïntroduceerde grootheid $\delta N^{\pm}$ fungeert als een indicator voor lokale optimaliteit. Als deze waarde dicht bij nul ligt, is de stencil-grootte waarschijnlijk optimaal.
• Beperking en Toekomst: De huidige indicator vereist nog de analytische oplossing om te worden berekend, wat in praktische toepassingen vaak niet beschikbaar is. De auteurs zien dit als een eerste stap naar een indicator die zonder de exacte oplossing kan werken.
• Impact: Het begrijpen en voorspellen van deze optimale stencil-groottes kan de nauwkeurigheid van RBF-FD methoden aanzienlijk verbeteren zonder extra rekenkosten (zoals het verhogen van de discretisatie-dichtheid of de orde van augmentatie).

Conclusie

De auteurs tonen aan dat de nauwkeurigheid van PHS RBF-FD methoden niet lineair verbetert met een grotere stencil, maar oscilleert. Dit gedrag is inherent aan de methode en niet aan het specifieke probleem. Door de correlatie tussen de oscillatie en de verdeling van fouttekens te identificeren, bieden ze een nieuw perspectief voor het selecteren van de meest efficiënte stencil-grootte, wat essentieel is voor het optimaliseren van meshvrije methoden.