An Eulerian approach for simulations of wave propagation in multiphase, viscoelastic media is developed in the context of the Advection Upstream Splitting Method (AUSM). We extend the AUSM scheme to the five-equation model for simulations of interfaces between gases, liquids, and solids with constitutive relations appropriately transported. In this framework, the solid’s deformations are assumed to be infinitesimally small such that they can be modeled using linear viscoelastic models, e.g., generalized Zener. The Eulerian framework addresses the challenge of calculating strains, more naturally expressed in a Lagrangian framework, by using a hypoelastic model that takes an objective Lie derivative of the constitutive relation to transform strains into velocity gradients. Our approach introduces elastic stresses in the convective fluxes that are treated by generalizing AUSM flux-vector splitting (FVS) to account for the Cauchy stress tensor. We determine an appropriate discretization of non-conservative equations that appear in the five-equation multiphase model with AUSM schemes to prevent spurious oscillations at material interfaces. The framework’s spatial scheme is solution adaptive with a discontinuity sensor discriminating between smooth and discontinuous regions. The smooth regions are computed using explicit high-order central differences. At discontinuous regions (i.e., shocks, material interfaces, and contact surfaces), the convective fluxes are treated using a high-order Weighted Essentially Non-Oscillatory (WENO) scheme with $$\hbox {AUSM}^+$$ -up for upwinding. The framework is used to simulate one-dimensional (1D) and two-dimensional (2D) problems that demonstrate the ability to maintain equilibrium interfacial conditions and solve challenging multi-dimensional and multi-material problems.