Improving the Accuracy and Efficiency of Hybrid Finite Element / Particle-In-Cell Methods for Modeling Geologic Processes

In many computational geodynamics models, properties such as density or viscosity are modeled with particles that carry the values of these properties and move with the underlying flow field. At each time step, the properties carried on the particles are interpolated onto the grid on which other quantities, such as the velocity, are computed. The most common interpolation method is to average the property values from the particles to the grid. More recently it has been shown that a Linear Least Squares (LLS) approximation to the particle values in each cell provides much greater accuracy. However, to achieve optimal accuracy, the LLS method requires the number of particles per cell (PPC) to continually increase as the grid size h goes to zero.

Here we demonstrate that a Quadratic Least Squares (QLS) interpolation requires only a small, fixed number of PPC, to achieve optimal accuracy as the cell size h goes to zero. In order to test the QLS method and compare it with the LLS method, we implemented both in the open source geodynamics code ASPECT and tested them on several 'classical' benchmarks, SolCx and SolKz, and on a new time-dependent benchmark in an annulus.