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.