Abstract: We develop a new high-order numerical method for continuum simulations of multimaterial phenomena in solids exhibiting elastic–plastic behavior using the diffuse interface numerical approximation. This numerical method ...