We explore the rich and complex dynamics that a phenomenological equation of state (EoS) with non-convex regions in the pressure-density plane may develop as a result of genuinely relativistic effects. We study the parameter space of the EoS to ensure its causality and thermodynamical consistency. We approximate the non-conventional dynamics developed in the evolution of relativistic blast waves by means of a high order shock capturing scheme.