A route map for modeling pharmaceutical manufacturing processes utilizing morphological population balance (MPB) is presented in terms of understanding and controlling particle shape and size for optimizing the efficiency of both the manufacturing process and final properties of the formulated drugs. This is applied to batch cooling crystallization of the pharmaceutical compound ibuprofen from supersaturated ethanolic solutions in which the MPB is combined with the known crystal morphology and associated face-specific growth kinetics (Nguyen et al. CrystEngComm 2014, 16, 4568-4586) to predict the temporal evolution of the shape and size distributions of all crystals. The MPB simulations capture the temporal evolution of the size and shape of ibuprofen crystals and their distributions at each time instance during the crystallization processes. The volume equivalent spherical diameter and crystal size distribution converted from MPB simulation are validated against the experimental studies on the 1 L scale size (Rashid et al. Chem. Eng. Res. Des. 2012, 90, 158-161), confirming the promise of this approach as a powerful simulation, optimization, and control tool for the digital design of precision pharmaceutical processes and products with the desirable properties, with potential applications in crystallization personalized medicines.