Two-derivative Runge–Kutta methods for PDEs using a novel discretization approach