  In this work we develop a high-resolution mapped-grid finite volume method
code to model wave propagation in two dimensions in systems of multiple
orthotropic poroelastic media and/or fluids, with curved interfaces between
different media. We use a unified formulation to simplify modeling of the
various interface conditions -- open pores, imperfect hydraulic contact, or
sealed pores -- that may exist between such media. Our numerical code is based
on the Clawpack framework, but in order to obtain correct results at a material
interface we use a modified transverse Riemann solution scheme, and at such
interfaces are forced to drop the second-order correction term typical of
high-resolution finite volume methods. We verify our code against analytical
solutions for reflection and transmission of waves at a material interface, and
for scattering of an acoustic wave train around an isotropic poroelastic
cylinder. For reflection and transmission at a flat interface, we achieve
second-order convergence in the 1-norm, and first-order in the max-norm; for
the cylindrical scatterer, the highly distorted grid mapping degrades
performance but we still achieve convergence at a reduced rate. We also
simulate an acoustic pulse striking a simplified model of a human femur bone,
as an example of the capabilities of the code. To aid in reproducibility, at
the web site http://dx.doi.org/10.6084/m9.figshare.701483 we provide all of the
code used to generate the results here.
