Background: Crystallization is a widely used unit operation in the pharmaceutical industry where high aspect ratio crystals are often produced. To avoid undesired shape i.e. high aspect ratio crystals according to downstream operations, appropriate level of understanding process is necessary. Breakage often occurs during the process and gives rise to change the Crystal shape. A well-adjusted, detailed model can be usable tool for investigation and optimization of control. Method: A two-dimensional population balance model of continuous cooling crystallization, involving nucleation, growth of two characteristic crystal facets and random binary breakage of high aspect ratio crystals is developed. The randomness of breakage is described by beta distribution of broken fractions. The population balance model is reduced into a closed moment equation model for the joint moments of the two size variables of crystals by means of which the influence of parameters of breakage on the dynamic and steady state behavior of crystallizer is analyzed according to understand the impact of complex kinetic on crystal shape. Results: Analysis of the process provides information of breakage influence on aspect ratio of crystal. It is shown that the mean value of crystal width increases as the intensity of breakage along the crystal length resulted in increasing number of crystals. (C) 2013 The Institution of Chemical Engineers. Published by Elsevier B.V. All rights reserved.