JustPIC.jl

Appearance

Velocity interpolation ​

Particle and marker advection needs velocities at particle positions, while geodynamic velocity fields are usually stored on staggered grids. JustPIC offers three interpolation paths for this step.

Linear ​

The default advection! path uses bilinear interpolation in 2D and trilinear interpolation in 3D. This is the cheapest and most general option, and it is the right default for most examples and tests.

LinP ​

advection_LinP! uses the linear-plus-pressure interpolation described by Pusok et al. 2017. The velocity at the m-th particle is given by

where is the bi- or trilinear interpolation from velocity nodes to the particle, is the interpolation from pressure nodes to the particle, and   is an empirical coefficient.

Modified Quadratic Spline ​

advection_MQS! uses the modified quadratic spline interpolation from Gerya et al. 2021. The scheme is designed so that velocity derivatives can be reconstructed from pressure-node locations where they are constrained by the continuity equation. Near boundaries, where the required stencil is unavailable, the implementation falls back to the linear interpolation path.

Example for the component in 2D:

Step 1: compute the normalized distances between the particle and the velocity node at the lower-left corner of the interpolation cell:

Step 2: lerp along on the bottom and top edges of the cell:

Step 3: add the quadratic correction, a second difference along over the three nodes closest to the particle. Which triplet that is depends on which half of the cell the particle sits in — the figure above shows the   case:

Step 4: lerp the corrected values along :

The component uses the same construction with the roles of the two directions exchanged.

Choosing a Scheme ​

  • Use advection! for the default linear interpolation.

  • Use advection_LinP! when matching the LinP reconstruction from the cited PIC literature.

  • Use advection_MQS! when the modified quadratic spline stencil is desired and the grid has enough interior support.

  • Use semilagrangian_advection!, semilagrangian_advection_LinP!, or semilagrangian_advection_MQS! for grid-field backtracking with the corresponding velocity reconstruction.