The mathematical optimisation of a batch cooling crystallization process is considered in this work. The objective is to minimize the standard deviation of the final crystal size distribution (CSD), which is an important feature in many industrial processes. The results with the problem written as a nonlinear programming and solved with the successive quadratic programming (SQP) coupled with the discretization of the control variable are compared with those obtained when SQP coupled with the parameterisation of the control variable is applied. Also it is proposed the implementation of the genetic algorithm (GA) coupled with parameterisation of the control variable. Extensive evaluations show that the SQP method is sensitive both to the parameterisation formulation and to the initial estimate. The solution with GA provided the control variable profile that leads to the minimum standard deviation of the final CSD. Nevertheless, it is a very time-consuming technique, which hampers its utilization in real time applications. However, its feature of global searching suggests its suitability in solving offline problems, in order to provide initial setup profiles. Bearing this in mind, it is proposed an algorithm which allows for the implementation of GA solution in a real time fashion, taking advantage of its robustness to find out the optimal solution. (c) 2005 Elsevier Ltd. All rights reserved.