A system, methods, and user-friendly program product to calculate global energy utility targets and define optimal driving force distribution for a process or cluster of processes under all possible process changes and streams specific minimum temperature approach values, simultaneously, and without enumeration, are provided. The program product can utilize stream-specific minimum temperature approach values ΔTmini, where the superscript i represents the specific hot stream, as the optimization parameters instead of the single global ΔTmin currently used, in addition to identifying the optimal operating conditions. The program product can define optimal process conditions and an optimal driving force distribution in heat recovery systems, and can produce an optimal Pareto-curve that shows the rigorous trade off between energy cost and capital cost for any energy recovery system.