Non-trace full-F gyro-fluid interchange impurity advection

A full-F isothermal gyro-fluid model and code (which is based on the full distribution function F compared to only small fluctuations) is extended to handle self-consistent coupling of multiple quasi-neutral ion species via the polarisation equation in the long wavelength approximation. The numerical model is used to determine two-dimensional interchange driven 'blob' transport in a plasma with intrinsic impurity content for a range of impurity parameters. With the same model, the self-consistent advective interaction of a main plasma species blob with a non-trace impurity cloud is studied. For homogeneous impurity distributions an increased effective mass reduces blob transport, whereas it is found that localised impurity clouds can lead either to acceleration or slowing down of blob propagation depending on the alignment of the impurity density gradient during the acceleration phase of the main ion species blob.