The main objective of this research is to control the evolution of the Crystal Size Distribution (CSD) in a batch crystallization process. We are led to study a numerical resolution of the population balance and propose an algorithm for fast and accurate simulation. This method was extended to the case of the two-dimensional crystal by using a moving mesh. We studied the problem of controllability of the system from its discretized model and then from the continuous model. To compensate the uncertainty of the model parameters, we added the second state feedback control to ensure the tracking of the desired CSD in presence of parameter uncertainty. We constructed an observer who provides us with on-line estimation of state variables. These state variables estimated are used in the control law synthesis