  The r-process nuclei are robustly synthesized in the material ejected during
a neutron star binary merger (NSBM), as tidal torques transport angular
momentum and energy through the outer Lagrange point in the form of a vast
tidal tail. If NSBM are indeed solely responsible for the solar system r-
process abundances, a galaxy like our own would require to host a few NSBM per
million years, with each event ejecting, on average, about 5x10^{-2} M_sun of
r-process material. Because the ejecta velocities in the tidal tail are
significantly larger than in ordinary supernovae, NSBM deposit a comparable
amount of energy into the interstellar medium (ISM). In contrast to extensive
efforts studying spherical models for supernova remnant evolution, calculations
quantifying the impact of NSBM ejecta in the ISM have been lacking. To better
understand their evolution in a cosmological context, we perform a suite of
three-dimensional hydrodynamic simulations with optically-thin radiative
cooling of isolated NSBM ejecta expanding in environments with conditions
adopted from Milky Way-like galaxy simulations. Although the remnant morphology
is highly complex at early times, the subsequent radiative evolution that
results from thermal instability in atomic gas is remarkably similar to that of
a standard supernova blast wave. This implies that sub-resolution supernova
feedback models can be used in galaxy-scale simulations that are unable to
resolve the key evolutionary phases of NSBM blast waves. Among other
quantities, we examine the radius, time, mass and kinetic energy content of the
NSBM remnant at shell formation as well as the momentum injected to the ISM. We
find that the shell formation epoch is attained when the swept-up mass is about
10^3 M_sun, at this point the mass fraction of r-process material is
drastically enhanced up to two orders of magnitude in relation to a solar
metallicity ISM.
