Counter differential rigid-rotation equilibrium of electrically non-neutral two-fluid plasma with finite pressure

We derive the two-dimensional counter-differential rotation equilibria of two-component plasmas, composed of both ion and electron ( $e^-$ ) clouds with finite temperatures, for the first time. In the equilibrium found in this study, as the density of the $e^{-}$ cloud is always larger than that of the ion cloud, the entire system is a type of non-neutral plasma. Consequently, a bell-shaped negative potential well is formed in the two-component plasma. The self-electric field is also non-uniform along the $r$ -axis. Moreover, the radii of the ion and $e^{-}$ plasmas are different. Nonetheless, the pure ion as well as $e^{-}$ plasmas exhibit corresponding rigid rotations around the plasma axis with different fluid velocities, as in a two-fluid plasma. Furthermore, the $e^{-}$ plasma rotates in the same direction as that of $\boldsymbol {E \times B}$ , whereas the ion plasma counter-rotates overall. This counter-rotation is attributed to the contribution of the diamagnetic drift of the ion plasma because of its finite pressure.

A solver based on pseudo-spectral analytical time-domain method for the two-fluid plasma model

A number of physical processes in laser-plasma interaction can be described with the two-fluid plasma model. We report on a solver for the three-dimensional two-fluid plasma model equations. This solver is particularly suited for simulating the interaction between short laser pulses with plasmas. The fluid solver relies on two-step Lax–Wendroff split with a fourth-order Runge–Kutta scheme, and we use the Pseudo-Spectral Analytical Time-Domain (PSATD) method to solve Maxwell’s curl equations. Overall, this method is only based on finite difference schemes and fast Fourier transforms and does not require any grid staggering. The Pseudo-Spectral Analytical Time-Domain method removes the numerical dispersion for transverse electromagnetic wave propagation in the absence of current that is conventionally observed for other Maxwell solvers. The full algorithm is validated by conservation of energy and momentum when an electromagnetic pulse is launched onto a plasma ramp and by quantitative agreement with wave conversion of p-polarized electromagnetic wave onto a plasma ramp.