This function computes the specific internal energy of the fluid, but neglects the (small) influence of the pressure term p/d.