Selecting suitable solvents for the crystallization of pharmaceuticals can be a challenging task, given the vast number of solvents that we can choose from.To simplify this problem, we can use principles from solid-liquid equilibrium (SLE) alongside established thermodn. models to identify promising candidates prior to conducting experimentsThis study addresses the batch cooling crystallization of flurbiprofen - a non-steroidal anti-inflammatory drug (NSAID) used to treat arthritis - using a simple model framework implemented within MATLAB.The Apelblat equation is employed to describe the thermophys. behavior of flurbiprofen over a wide temperature range (283.15-323.15 K) in twelve (12) solvents: three alkanes (n-hexane, n-heptane, n-octane); two (iso-Pr, methyl-tert-butyl) ethers; five alcs. (n-propanol, isopropanol, n-butanol, isobutanol, isopentanol); an ester (iso-Pr acetate); and a nitrile (acetonitrile).Moreover, we have used green metrics (E-factor, Scope 1 and 2 carbon emissions) with established process economics models to determine the most promising solvent (n-propanol) for an environmentally friendly and economical manufacturing process.