returns internal energy derivative (d2u/dd2)@T=const
function d2ud2T input EoS.HelmholtzDerivs f; output Der2EnergyByDensity2 d2ud2T; end d2ud2T;