A unified approach to optimal control of univariate and multivariate crystallization particulate processes with size-independent or/and size-dependent growth rate kinetics is developed by utilizing the minimum principle and the method of characteristics in conjunction with novel approximate integro-differential and ordinary differential equation computational schemes. The proposed theoretical approach leads to simple analytic solutions involving numerical boundary value problems with few unknown parameters. (C) 2012 Elsevier Ltd. All rights reserved.