Boundary integral methods for fast and accurate simulation of droplets in two-dimensional Stokes flow

Accurate simulation of viscous fluid flows with deforming droplets creates a number of challenges. This thesis identifies these principal challenges and develops a numerical methodology to overcome them. Two-dimensional viscosity-dominated fluid flows are exclusively considered in this work. Such flows find many applications, for example, within the large and growing field of microfluidics; accurate numerical simulation is of paramount importance for understanding and exploiting them. A boundary integral method is presented which enables the simulation of droplets and solids with a very high fidelity. The novelty of this method is in its ability to accurately handle close interactions of drops, and of drops and solid boundaries, including boundaries with sharp corners. The boundary integral method is coupled with a spectral method to solve a PDE for the time-dependent concentration of surfactants on each of the droplet interfaces. Surfactants are molecules that change the surface tension and are therefore highly influential in the types of flow problems which are considered herein. A method’s usefulness is not dictated by accuracy alone. It is also necessary that the proposed method is computationally efficient. To this end, the spectral Ewald method has been adapted and applied. This yields solutions with computational cost O(N log N ), instead of O(N^2), for N source and target points. Together, these innovations form a highly accurate, computationally efficient means of dealing with complex flow problems. A theoretical validation procedure has been developed which confirms the accuracy of the method.Att noggrant simulera viskösa flöden med deformerande droppar medför flera utmaningar. Denna avhandling identifierar de viktigaste utmaningarna och utvecklar numeriska metoder för att övervinna dem. Visköst dominerade tvådimensionella flöden studeras. Sådana flöden har många tillämpningar till exempel inom mikrofluidik och noggrann beräkning är av största vikt för att förstå och utnyttja dem. En randintegralsmetod som möjliggör simulering av droppar och fasta ränder med en mycket hög noggrannhet presenteras. Metoden särskiljer sig från andra genom dess förmåga att hantera nära samspel mellan droppar och förekomst av hörn på de fasta ränderna. Randin- tegralsmetoden är kopplad till en spektral metod som möjliggör inkluderandet av surfaktanter i flödesproblemet. Surfaktanter är molekyler som förändrar ytspänningen och de är därför betydelsefulla för de typer av flödesproblem som beaktas här. En metods användbarhet bestäms inte endast av dess noggrannhet. Det är också nödvändigt att den föreslagna metoden är effektiv. För detta ändamål har metoden spektral Ewald anpassats och tillämpats. Detta ger lösningar med beräkningskostnaden O(N log N ) istället för O(N^2), där N är antalet diskreta punkter i systemet. Tillsammans utgör dessa innovationer ett mycket noggrant, be- räkningseffektivt sätt att hantera komplexa flödesproblem. Ett teoretiskt valideringsförfarande har utvecklats som bekräftar metodens noggrannhet