2.15 Particles

ASPECT can, optionally, also deal with particles (sometimes called “tracers”). Particles can be thought of as point-like objects that are simply advected along with the flow. In other words, if $u\left(x,t\right)$ is the flow field that results from solving equations (1)–(2), then the $k$th particle’s position satisfies the equations

The initial positions of all particles also need to be given and are usually either chosen randomly, based on a fixed pattern, or are read from a file.

Particles are typically used to track visually where material that starts somewhere ends up after some time of a simulation. It can also be used to track the history of the volume of the fluid that surrounds a particle, for example by tracking how much strain has accumulated, or what the minimal or maximal temperature may have been in the medium along the trajectory of a particle. To this end, particles can carry properties. These are scalar- or vector-valued quantities that are attached to each particle, that are initialized at the beginning of a simulation, and that are then updated at each time step. In other words, if we denote by $p_{k,m}\left(t\right)$ the value of the $m$th property attached to the $k$th particle, then $p_{k,m}\left(t\right)$ will satisfy a differential equation of the form

The exact form of $g_m$ of course depends on what exactly a particular property represents. Like with compositional fields (see Section 2.7), it is possible to describe the right hand side $g_m$ in ways that also allows for impulse (delta) functions in time.

How particles are used in practice is probably best explained using examples. To this end, see in particular Section 5.2.5. All particle-related input parameters are listed in Section ??. The implementation of particles is discussed in great detail in [GHPB17].