In this paper, we introduce a coarse-grained model of polymer crystallization using a multiphase-field approach. The model combines a multiphase-field method, Nakamura’s kinetic equation, and the equation of heat conduction for studying microstructural evolution of crystallization under isothermal and non-isothermal conditions. The multiphase-field method provides flexibility in adding any number of phases with different properties making the model effective in studying blends or composite materials. We apply our model to systems of neat PA6 and study the impact of initial distribution of crystalline grains and cooling rate on the morphology of the system. The relative crystallinity (conversion) curves show qualitative agreement with experimental data. We also investigate the impact of including carbon fibers on the crystallization and grain morphology. We observe a more homogeneous crystal morphology around fibers. This is associated with the higher initial volume fraction of crystal grains and higher heat conductivity of the fiber (compared to the polymer matrix). Additionally, we observe that the crystalline grains at the fiber surface grow perpendicular to the surface. This indicates that the vertical growth observed in experiments is merely due to geometrical constraints imposed by the fiber surface and neighbouring crystalline regions.