The growth rate and expansion history of the Universe can be measured from large galaxy redshift surveys using the Alcock-Paczynski effect. We validate the Redshift Space Distortion models used in the final analysis of the Sloan Digital Sky Survey (SDSS) extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 quasar clustering sample, in configuration and Fourier space, using a series of halo occupation distribution mock catalogues generated using the OuterRim N-body simulation. We test three models on a series of non-blind mocks, in the OuterRim cosmology, and blind mocks, which have been rescaled to new cosmologies, and investigate the effects of redshift smearing and catastrophic redshifts. We find that for the non-blind mocks, the models are able to recover fs8 to within 3 per cent and alpha(parallel to) and alpha(perpendicular to) to within 1 per cent. The scatter in the measurements is larger for the blind mocks, due to the assumption of an incorrect fiducial cosmology. From this mock challenge, we find that all three models perform well, with similar systematic errors on f sigma(8), alpha(parallel to), and alpha(perpendicular to) at the level of sigma f sigma(8) = 0.013, sigma(alpha parallel to) = 0.012, and sigma(alpha perpendicular to) = 0.008. The systematic error on the combined consensus is sigma f sigma(8) = 0.011, sigma(a) = 0.008, and sa. = 0.005, which is used in the final DR16 analysis. For baryon acoustic oscillation fits in configuration and Fourier space, we take conservative systematic errors of sigma(alpha parallel to) = 0.010 and sigma(alpha perpendicular to) = 0.007.