The free energy of a thermodynamic system is known to be a concave function of the temperature. This fact is used, together with some results due to Fisher, to deduce bounds on the internal energy of the two-dimensional Ising model, given reasonably accurate upper and lower bounds on the free energy. These free energy bounds are derived from renormalization group transformations. Unfortunately, the numerical accuracy of the bounds on the internal energy is poor. The reasons for the failure are discussed.