Solid phase solubility is a fundamental parameter in the design of crystallization processes. The development and optimization of crystallization processes requires screening of numerous solvent systems for which the solubility of the compound of interest has to be measured as a function of temperature and solvent composition. Tools that quickly estimate the solubility in different solvents can be very useful in the initial phases of the solvent system selection process. In this paper, we report our experience applying two thermodynamic models in the solubility estimation of pharmaceutical compounds: the NRTL-SAC method (Chen and Song, 2004, Ind Eng Chem Res 43: 8354) which provides a correlative and predictive model from limited solubility measurements, and the COSMO-SAC (Lin and Sandler, 2002, Ind Eng Chem Res 41: 899) method which predicts solubility from ab initio calculations. These theoretical methods, coupled with rapid experimental measurement for verification, provide a powerful solubility screening protocol for the development of crystallization processes.