Quadrature and Fast Summation for Boundary Integral Methods

QC 2026-09-16

Abstract

I denna avhandling utvecklar vi noggranna, stabila och effektiva algoritmer för två utmaningar inom randintegralmetoder: numerisk integration, eller kvadratur, av nästan singulära integraler och snabb summering av många parvisa interaktioner. Randintegralmetoder omformulerar linjära elliptiska partiella differentialekvationer till integralekvationer på randen, vilket minskar problemets dimension och möjliggör diskretiseringar av hög ordning för släta geometrier. Exempel på tillämpningar är elektrostatik, vågspridning vid en fix frekvens och rörelsen hos stela partiklar i viskösa vätskor.

De randintegraler som representerar lösningen, så kallade skiktpotentialer, är svåra att beräkna i punkter på eller nära randen eftersom deras integrander är singulära på randen och varierar snabbt i dess närhet. Standardkvadratur förlorar noggrannhet när evalueringspunkten, eller målpunkt, närmar sig randen, medan det blir kostsamt att bibehålla noggrannheten genom att öka antalet kvadraturpunkter. För målpunkter tillräckligt nära randen krävs därför specialkvadratur. Vi utvecklar en sådan metod för släta rotationssymmetriska ytor, med parametrar som väljs automatiskt för att uppfylla en given tolerans för kvadraturfelet. Metoden tillämpas på skiktpotentialer för de Laplace-, Helmholtz- och Stokesekvationer, vilka ligger till grund för tillämpningarna ovan. För vissa starkt singulära integraler kan stora termer nästan ta ut varandra vid beräkning mycket nära randen och orsaka en betydande precisionsförlust. Vi analyserar denna effekt och utvecklar en stabil omformulering som inte medför någon betydande extra beräkningskostnad.

Med utgångspunkt i dessa resultat konstruerar vi ett toleransstyrt kvadraturramverk som används i en randintegrallösare för Stokesflöde kring stela partiklar. För varje interaktion mellan en målpunkt och en partikel används förberäknade felindikatorer för att välja det minst kostsamma kvadraturalternativ som bedöms uppfylla toleransen. Den stabiliserade specialkvadraturen används därmed endast där den behövs.

Den andra utmaningen uppstår när standardkvadratur används för många målpunkter, eftersom bidraget från varje kvadraturpunkt på randen måste beräknas för varje målpunkt. Med N kvadraturpunkter på randen och ungefär lika många målpunkter kräver direkt summering O(N²) beräkningsarbete. Metoder för snabb summering minskar denna kostnad genom att undvika en explicit beräkning av samtliga parvisa interaktioner. Vi utvidgar en sådan adaptiv metod, kallad dual-space multilevel kernel splitting, så att den kan användas för interaktioner mellan punkter i rätblock med godtyckliga proportioner och periodiska randvillkor i en, två eller tre koordinatriktningar. Vid given noggrannhet minskar metoden beräkningskostnaden till O(N). Vi tillämpar metoden på interaktionssummor som uppstår för Laplace- och Stokesekvationerna.

Tillsammans bildar dessa metoder ett ramverk där interaktioner som kräver specialbehandling identifieras och beräknas noggrant och stabilt, medan övriga interaktioner beräknas effektivt med snabb summering. Metoderna utgör byggstenar för tillförlitliga och effektiva storskaliga simuleringar av suspensioner av stela partiklar, inklusive system med interaktioner mellan närliggande partiklar och periodiska randvillkor. Sådana simuleringar kan användas för att studera partiklars kollektiva rörelse och suspensioners reologi. Kvadraturmetoderna kan även användas för simuleringar av fibrer och filament inom biofysik på cellskala, medan den snabba summeringsmetoden är mer allmänt tillämpbar på problem med långräckviddiga interaktioner, däribland elektrostatiska beräkningar inom molekyldynamik.

Link to DiVA