We demonstrate a fast, flexible, and accurate paraxial wave propagation model to serve as a forward model for propagation-based X-ray phase contrast imaging (XPCI) in parallel-beam or cone-beam geometry. This model incorporates geometric cone-beam effects into the multi-slice beam propagation method. It enables rapid prototyping and is well suited to serve as a forward model for propagation-based X-ray phase contrast tomographic reconstructions. Furthermore, it is capable of modeling arbitrary objects, including those that are strongly or multi-scattering. Simulation studies were conducted to compare our model to other forward models in the X-ray regime, such as the Mie and full-wave Rytov solutions.