Plasma turbulence within cometary plasma environments

We present a numerical work in which the interaction between a comet and the solar wind is studied in 2D in the plane perpendicular to the solar wind mean field direction. Our simulations are conducted with the hybrid Particle-in-Cell (PIC) code Menura that allows for the injection of a turbulent solar wind [1]. First, we consider the case of laminar solar wind and we present a study on the equivalent Mach number of the two-ion-species (cometary and solar wind) plasma surrounding the comet. We develop an expression for the Mach number having suitable limits in the two asymptotic cases of infinite cometary and solar wind ion density; our expression is derived by extending previous studies on bi-ion plasma models [2]. Through numerical simulations in which the cometary activity is varied, we show how our Mach number is able to unambiguously describe the existence and location of the cometary shock. Second, we compare two runs, one with a laminar and one with a turbulent solar wind in the case of moderate cometary activity. We divide the simulation domain into the regions upstream and downstream the cometary shock. We analyze how plasma turbulence properties are affected by the passage through the shock in the case of a turbulent solar wind. Then, we divide the downstream region into three different regions identified by different solar wind-to-cometary ion density ratios. We study the downstream turbulence properties in the case of laminar and turbulent impinging solar wind and how they vary in those regions.