Start with the Schwarzschild solution in the weak field approximation:
[tex](cd\tau)^2=(1-\frac{2\Phi}{c^2})(cdt)^2+(1-\frac{2\Phi}{c^2})^{-1}(dr)^2+...[/tex]
That should be:
[tex]{\color{red}(cd\tau)^2=(1+\frac{2\Phi}{c^2})(cdt)^2-(1+\frac{2\Phi}{c^2})^{-1}(dr)^2-...}[/tex]
See post 185 of this different thread
https://www.physicsforums.com/showthread.php?t=397403&page=12 where you gave the correct equation. This error in the signs propagates all the way through the rest of your calculations.
For the case [tex]dr=d\theta=d\phi=0[/tex] you get the well known relationship:
[tex]d\tau=\sqrt{1{\color{red}+}\frac{2\Phi}{c^2}}dt[/tex]
Writing the above for two different gravitational potentials [tex]\Phi_1[/tex] and [tex]\Phi_2[/tex] you obtain the well-known time dilation relationship:
[tex]\frac{d\tau_1}{d\tau_2}=\sqrt{\frac{1{\color{red}+}\frac{2\Phi_1}{c^2}}{1{\color{red}+}\frac{2\Phi_2}{c^2}}}[/tex]
At the Earth surface :
[tex]\Phi_1=-\frac{GM}{R}[/tex]
At the Earth center:
[tex]\Phi_2=-3/2\frac{GM}{R}[/tex]
Now, due to the fact that [tex]\frac{\Phi}{c^2}<<1[/tex] you can obtain the approximation:
[tex]\frac{d\tau_1}{d\tau_2}=1{\color{red}+}\frac{\Phi_1-\Phi_2}{c^2}=1{\color{red}+}\frac{GM}{2Rc^2}{\color{red}>}1[/tex]
So, [tex]f_1>f_2[/tex] where [tex]f_1[/tex] is the clock frequency on the Earth crust and [tex]f_2[/tex] is the frequency of the clock at the center of the Earth.
When corrected, your aproximation is in close agreement with the equations I gave and is on the right side of unity.
...