  A self-consistent model for the simulation of Ge-rich Ge$_2$Sb$_2$Te$_5$
phase change memories is presented. Combining the multi-phase field model and a
phase-aware electro-thermal solver, it reproduces the multi-physics behavior of
the material. Simulations of memory operations are performed to demonstrate its
ability to reproduce experimental observations.
