Particleincell methods have a long history in modelling of mantle convection, lithospheric deformation and crustal dynamics. They are primarily used to track material information, the strain a material has undergone, the pressuretemperature history of a certain material, or the amount of volatiles or partial melt present in a region. However, their efficient parallel implementation  in particular combined with adaptive meshes  is complicated due to the complex communication and frequent reassignment of particles to cells. Consequently, many scientific software packages accomplish this efficiency by designing particle methods for a single purpose, like the advection of scalar properties that do not evolve over time (e.g. chemical heterogeneities). Design choices for particle advection, data storage, and parallel communication are then optimized for this single purpose, making the code rigid to changing requirements. Here, I present the implementation of a flexible, scalable and efficient particleincell method for massively parallel finiteelement codes with adaptively changing meshes. A modular plugin structure allows for maximum flexibility of the generation of particles, the carried particle properties, the advection and output algorithms, and the projection of properties to the finiteelement mesh. I discuss the complexity of the implemented algorithms and present scaling tests ranging up to tens of thousands of cores and tens of billions of particles. Additionally, I discuss loadbalancing strategies like balanced repartitioning for particles in adaptive meshes with their strengths and weaknesses, and quantify sources of errors for the advection of particles. Finally, I show whole mantle convection models as application cases, compare our implementation to a modern advectionfield approach, and demonstrate under which conditions which method is more efficient. I implemented the presented methods in ASPECT (http://aspect.geodynamics.org), a freely available opensource community code for geodynamic simulations.
You are here
Implementing flexible and scalable particleincell methods for massively parallel numerical models
Research Group:
Rene Gassmoeller
Location:
A005
Schedule:
Thursday, September 14, 2017  11:00
Abstract:
Openings
 Public calls for academic personnel (Permanent positions)
 Professors (Temporary/Researchers/Visiting Professors)
 SISSA Mathematical Fellowships
 Post Doctoral Fellowships
 PhD Scholarships
 Call for Applications (PhD)
 Undergraduate Fellowships
 Postgraduate Fellowships
 Master of Science in Mathematics
 Marie SklodowskaCurie Grants
Upcoming events

Martijn Kool
Donaldson Thomas invariants on CalabiYau fourfolds
Thursday, February 22, 2018  14:30 to 15:30

Felice Iandoli
On the Cauchy problem for quasilinear Schrödingertype equations on the circle
Friday, February 23, 2018  14:00

Xiao Han
Principal fibrations over noncommutative spheres
Monday, February 26, 2018  16:00 to 17:00

Carlangelo Liverani
Mathematical Colloquium: Beyond averaging
Monday, March 12, 2018  16:00
Today's Lectures

Stefano Bianchini
09:00 to 11:00

Andrei Agrachev
11:00 to 13:00

Nicola Gigli
14:00 to 16:00
Recent publications

G. Cotti; B. Dubrovin; D. Guzzetti,Local moduli of semisimple Fro...

A. Michelangeli; A. Ottolini; R. Scandone,Fractional powers and singular...

R. Alicandro; G. Lazzaroni; M. Palombaro,Derivation of a rod theory fro...

G. Cotti; B. Dubrovin; D. Guzzetti,Isomonodromy Deformations at a...