Phase Diversity (PD) is a technique widely used to measure and compensate for Non-Common Path Aberrations (NCPA) in optical systems. Beyond its traditional applications, PD has found further uses in solar astrophysics, where it enhances image quality by processing multiple images to correct residual aberrations. In this work, we introduce a simulator designed to facilitate the development of PD algorithms, aimed at enabling researchers to produce higher-performance solutions of these algorithms.