This work investigates crystallization modeling by modifying an open-source computational fluid dynamics code OpenFOAM. The crystallization behavior of high-density polyethylene (HDPE) is implemented according to theoretical and experimental literature. A number of physical interdependencies are included. The cavity is modeled as deformable. The heat transfer coefficient in the thermal contact towards the mold depends on contact pressure. The thermal conductivity is pressure- and crystallinity-dependent. Specific heat depends on temperature and crystallinity. Latent heat is released according to the crystallization progress and temperature. Deviatoric elastic stress is evolved in the solidified material. The prediction of the cavity pressure evolution is used for the assessment of the solution quality because it is experimentally available and governs the residual stress development. Insight into the thermomechanical conditions is provided with through-thickness plots of pressure, temperature and cooling rate at different levels of crystallinity. The code and simulation setup are made openly available to further the research on the topic.