Collisions are an important dissipation mechanism in plasmas. When approximating collision operators numerically, it is important to preserve their mathematical structure in order to retain the laws of thermodynamics at the discrete level. This is particularly challenging when considering particle methods. A simple but commonly used collision operator is the Lenard–Bernstein operator, or its modified energy- and momentum-conserving counterpart. In this work, we present a macro-particle discretisation of this operator that is provably energy and momentum preserving.
In a recent publication (Toler et al., J. Plasma Phys., vol. 89, issue 2, 2023, p. 905890210), we demonstrated that for axisymmetric geometries, the Kapur–Rokhlin quadrature rule provided an efficient and high-order accurate method for computing the normal component, on the plasma surface, of the magnetic field due to the toroidal current flowing in the plasma, via the virtual-casing principle. The calculation was indirect, as it required the prior computation of the magnetic vector potential from the virtual-casing principle, followed by the computation of its tangential derivative by Fourier differentiation, to obtain the normal component of the magnetic field. Our approach did not provide the other components of the virtual-casing magnetic field. In this letter, we show that a more direct and more general approach is available for the computation of the virtual-casing magnetic field. The Kapur–Rokhlin quadrature rule accurately calculates the principal value integrals in the expression for all the components of the magnetic field on the plasma boundary, and the numerical error converges at a rate nearly as high as the indirect method we presented previously.
Linearized gravity around a rotating black hole or compact object introduces the concept of a gravitomagnetic field, which originates from the matter–current in the rotating object. Plasma in proximity to this object is subsequently subjected to motion guided by this gravitomagnetic term (where mass serves as the effective charge) in addition to the conventional magnetic field. Such an interplay of fields complicates the accessible plasma waves of the system and thus merits exploration to delineate this interplay, to identify any observable signatures of the gravitomagnetic field in weakly magnetized systems, and to motivate future numerical work in a fully relativistic setting where the effects may be stronger. In this work we analyse the dispersion of the upper and lower hybrid electrostatic waves in a plasma immersed in both magnetic and gravitomagnetic fields. In particular, we discuss the effective augmentation or cancellation of the two fields under the right conditions for the upper hybrid wave. In contrast, the lower hybrid wave experiences a frequency up-shift from the gravitomagnetic field regardless of whether it is parallel or antiparallel to the magnetic field for the studied field strengths.