  We present a solver of 3D two-fluid plasma model for the simulation of
short-pulse laser interactions with plasma. This solver resolves the equations
of the two-fluid plasma model with ideal gas closure. We also include the
Bhatnagar-Gross-Krook collision model. Our solver is based on PseudoSpectral
Time-Domain (PSTD) method to solve Maxwell's curl equations. We use a Strang
splitting to integrate Euler equations with source term: while Euler equations
are solved with a composite scheme mixing Lax-Friedrichs and Lax-Wendroff
schemes, the source term is integrated with a fourth-order Runge-Kutta scheme.
This two-fluid plasma model solver is simple to implement because it only
relies on finite difference schemes and Fast Fourier Transforms. It does not
require spatially staggered grids. The solver was tested against several
well-known problems of plasma physics. Numerical simulations gave results in
excellent agreement with analytical solutions or with previous results from the
literature.

Plasma Physics;Computational Physics