The implementation, optimization, and performance of DFT-D, including the effects of solvation, has been tested on applications of polar processes in solution, where dispersion and hydrogen bonding is known to be involved. Solvent effects are included using our ab initio continuum solvation strategy, COSab, a conductor-like continuum solvation model, modified for ab initio in the quantum chemistry program GAMESS. Structure and properties are investigated across various functionals to evaluate their ability to properly model dispersion and solvation effects. The commonly used S22 set with accurate interaction energies of organic complexes has been used for parametrization studies of dispersion parameters and relevant solvation parameters. Dunning’s correlation consistent basis sets, cc-pVnZ (n = D, T), are used in the optimization, together with the Grimme B97-D exchange-correlation functional. Both water (ε = 78.4) and ether (ε = 4.33) environments are considered. Optimized semiempirical dispersion correction parameters and solvent extent radii are proposed for several functionals. We find that special parametrization of the semiempirical dispersion correction when used together in the DFT-D/COSab approach is not necessary. The global performance is quite acceptable in terms of chemical accuracy and suggests that this approach is a reliable as well as economical method for evaluation of solvent effects in systems with dispersive interactions. The resulting theory is applied to a group of push−pull pyrrole systems to illustrate the effects of donor/acceptor and solvation on their conformational and energetic properties.