The Abraham general solvation model is used to predict the saturation solubility of crystalline non-electrolyte solutes in organic solvents. The derived equations take the form of log(CS/CW) = c + rR(2) + spi(2)(H) + a Sigma alpha(2)(H) + b Sigma beta(2)(H) + nuV(x) log(CS/CG) = c + rR(2) + spi(2)(H) + a Sigma alpha(2)(H) + b Sigma beta(2)(H) l log L-16 where C-S and C-W refer to the solute solubility in the organic solvent and water, respectively, C-G is a gas phase concentration, R-2 the solute's excess molar refraction, V-x the McGowan volume of the solute, Sigma alpha(2)(H) and Sigma beta(2)(H) are the measures of the solute's hydrogen-bond acidity and hydrogen-bond basicity, pi(2)(H) denotes the solute's dipolarity/polarizability descriptor, and L-16 is the solute's gas phase dimensionless Ostwald partition coefficient into hexadecane at 298 K. The remaining symbols in the above expressions are known coefficients, which have been determined previously for a large number of gas/solvent and water/solvent systems. Computations show that the Abraham general solvation model predicts the observed solubility behavior of pyrene, acenaphthene and fluoranthene to within an average absolute deviation of about +45%.