Try without Evaluate:
TestParametricPlotC = ParametricPlot[{{{q[\[Phi]], q'[\[Phi]]} /. NumericalDiffequaWith1}, {{q[\[Phi]], q'[\[Phi]]} /. NumericalDiffequaWithout1}}, {\[Phi], 0, 50}, AxesLabel -> {q[\[Phi]], \[Phi]}, PlotRange -> All]