  In various astrophysics settings it is common to have a two-fluid
relativistic plasma that interacts with the electromagnetic field. While it is
common to ignore the displacement current in the ideal, classical
magnetohydrodynamic limit, when the flows become relativistic this
approximation is less than absolutely well-justified. In such a situation, it
is more natural to consider a positively charged fluid made up of positrons or
protons interacting with a negatively charged fluid made up of electrons. The
two fluids interact collectively with the full set of Maxwell's equations. As a
result, a solution strategy for that coupled system of equations is sought and
found here. Our strategy extends to higher orders, providing increasing
accuracy. Three important innovations are reported here. In our first
innovation, the magnetic field within each zone is reconstructed in a
divergence-free fashion while the electric field within each zone is
reconstructed in a form that is consistent with Gauss' law. In our second
innovation, a multidimensionally upwinded strategy is presented which ensures
that the magnetic field can be updated via a discrete interpretation of
Faraday's law and the electric field can be updated via a discrete
interpretation of the generalized Ampere's law. Our third innovation consists
of an efficient design for several popular IMEX schemes so that they provide
strong coupling between the finite-volume-based fluid solver and the
electromagnetic fields at high order. Several accuracy analyses are presented
showing that our method meets its design accuracy in the MHD limit as well as
in the limit of electromagnetic wave propagation. Several stringent test
problems are also presented. We also present a relativistic version of the GEM
problem, which shows that our algorithm can successfully adapt to challenging
problems in high energy astrophysics.
