High Performance Finite Element Methods with Application to Simulation of Vertical Axis Wind Turbines and Diffusion MRI

• Finite element methods have been developed over decades, and together with the growth of computer power, they become more and more important in dealing with large-scale simulations in science and industry.The objective of this thesis is to develop high-performance finite element methods, with two concrete applications: computational fluid dynamics (CFD) with simulation of turbulent flow past a vertical axis wind turbine (VAWT), and computational diffusion magnetic resonance imaging (CDMRI). The thesis presents contributions in the form of both new numerical methods for high-performance computing frameworks and efficient, tested software, published open source as part of the FEniCS/FEniCS-HPC platform. More specifically, we have four main contributions through the thesis work.First, we develop a DFS-ALE method which combines the Direct finite element simulation method (DFS) with the Arbitrary Lagrangian-Eulerian method (ALE) to solve the Navier-Stokes equations for a rotating turbine. This method is enhanced with dual-based a posteriori error control and automated mesh adaptation. Turbulent boundary layers are modeled by a slip boundary condition to avoid a full resolution which is impossible even with the most powerful computers available today. The method is validated against experimental data with a good agreement.Second, we propose a partition of unity finite element method to tackle interface problems. In CFD, it allows for imposing slip velocity boundary conditions on conforming internal interfaces for a fluid-structure interaction model. In CDMRI, it helps to overcome the difficulties that the standard approaches have when imposing the microscopic heterogeneity of the biological tissues and allows for efficient solutions of the Bloch-Torrey equation in heterogeneous domains. The method facilitates a straightforward implementation on the FEniCS/ FEniCS-HPC platform. The method is validated against reference solutions, and the implementation shows a strong parallel scalability.Third, we propose a finite element discretization on manifolds in order to efficiently simulate the diffusion MRI signal in domains that have a thin layer or a thin tube geometrical structure. The method helps to significantly reduce the required simulation time, computer memory, and difficulties associated with mesh generation, while maintaining the accuracy. Thus, it opens the possibility to simulate complicated structures at a low cost, for a better understanding of diffusion MRI in the brain.Finally, we propose an efficient portable simulation framework that integrates recent advanced techniques in both mathematics and computer science to enable the users to perform simulations with the Cloud computing technology. The simulation framework consists of Python, IPython and C++ solvers working either on a web browser with Google Colaboratory notebooks or on the Google Cloud Platform with MPI parallelization.

• Finita elementmetoder har utvecklats under årtionden, och har, till- sammans med tillväxten i datorkraft, blivit allt viktigare för att utföra storskaliga simuleringar inom både akademin och industrin. Målet med denna avhandling är att utveckla finita elementmetoder med högprestanda, med särskilt fokus på två konkreta applikationer; beräknings- strömningsdynamik (eng. Computational Fluid Dynamics (CFD)) för simulering av turbulent flöde runt en vindturbin, och beräkningar inom diffusionsmagnetresonanstomografi (eng. Computational diffusion magnetic resonance imaging (CDMRI)). Denna avhandling innehåller bidrag till ovanstående områden i form av såväl nya numeriska metoder för högprestandaberäkningsramverk och testad effektiv programvara vilken publicerats som öppen källkod som del av plattformen FEniCS/FEniCS-HPC. Mer specifikt presenterar vi fyra huvudbidrag i detta avhandlingsarbete.Först utvecklar vi en DFS-ALE-metod som kombinerar Direkt Fini- ta Elementsimulering (DFS) med den Arbiträra Lagrange-Eulermetoden (ALE) för att lösa Navier-Stokes ekvationer för en roterande turbin. Vår metod är en förbättrad variant med dualbaserad a posteriori felkontroll och automatiserad adaptering av beräkningsnätet. Turbulenta gränsskikt modelleras med ett sliprandvillkor för att undvika full upplösning av problemet, vilket är omöjligt även med de mest kraftfulla datorer som finns att tillgå idag. Metoden valideras mot experimentell data, med god överensstämmelse.Därnäst föreslår vi en enhetspartitions finita element metod för att tackla interfaceproblem. Inom CFD möjliggör detta att påtvinga ett sliprandvillkor på konforma inre interface för en fluidstrukturinter-kationsmodell. Inom CDMRI bidrar det med att överkomma svårigheterna med att påtvinga mikroskopisk heterogenitet av den biologiska vävnaden, och möjliggör effektiv lösning av Bloch-Torrey ekvationen i heterogena domäner. Metoden gör det enklare att göra en rättfram implementering i FEniCS/FEniCS-HPC. Metoden valideras mot referenslösnignar, och implementationen visar på stark parallel skalning (eng. strong parallel scaling).Sedan föreslår vi en finita elementdisktretisering på mångfalder för att effektivt kunna simulera diffusions-MRI-signaler i områden med en tunn geometrisk struktur. Metoden bidrar med att signifikant reducera simuleringstiden, minnesåtgång och svårigheter associerade med genereringen av beräkningsnät, utan att påverka precisionen i beräkningarna. Detta öppnar för möjligheter att simulera komplicerade strukturer till låg kostnad, för att bättre förstå diffusionsmagnettomografi i hjärnan.Tilll sist föreslår vi ett effektivt portabelt simuleringsramverk som integrear nya avancerade tekniker inom både matematik och datave- tenskap för att möjliggöra för användaren att utföra simuleringar med datormolnberäkningsteknologin. Simuleringsramverket består av Python, IPython och C++-lösare som används tillsammans antingen i en webbläsare med Google Colaboration notebooks eller på Google Cloud-plattformen med MPI-parallellisering.

Nyckelord

High performance finite element method

computational diffusion MRI

turbulent flow

vertical axis wind turbine

Cloud computing.

högprestanda finita elementmetod

beräkningsdiffusionsmagnetresonanstomografi

turbulent flöde

vertikalaxlad vindturbin

datormolnberäkning

Tillämpad matematik och beräkningsmatematik

Applied and Computational Mathematics

Biological Physics

Biologisk fysik

Computer Science

Datalogi

Engineering Mechanics

Teknisk mekanik

Mathematics

Matematik

Biological and Biomedical Physics

Biologisk och biomedicinsk fysik

Publikations- och innehållstyp

vet (ämneskategori)

dok (ämneskategori)