Analysis of the PHWR benchmark transient using a high-order spherical harmonics-nodal collocation method and treatment of rod cusping effect. (May 2021)