A simulation framework for coherent X-ray imaging, based on scalar diffraction theory, is presented. It contains a core C++ library and an additional Python interface. A workflow is presented to include contributions of inelastic scattering obtained with Monte-Carlo methods. X-ray Talbot-Lau interferometry is the primary focus of the framework. Simulations are in agreement with measurements obtained with such an interferometer. Especially, the dark-field signal of densely packed PMMA microspheres is predicted. A realistic modeling of the microsphere distribution, which is necessary for correct results, is presented. The framework can be used for both setup design and optimization but also to test and improve reconstruction methods.