Locus problem on circle and parabola 
Circles are drawn through the vertex of a parabola $y^2 = 4ax$ to cut the parabola orthogonally at the other point. Find the locus of the centers of the circles.

What I did:
The tangent must surely be the diameter of the circle. Also, the line through the vertex perpendicular to the line joining the vertex to the point on parabola must intersect the circle at the other end of the diameter i.e. where the tangent meets the circle again. The center of the circle must be their midpoint

So, I take a parametric point $P(at^2, 2at)$ on the parabola. The equation of the line perpendicular to $VP$ going through $V$ is $2y = -tx$. It meets the tangent $ty = x + at^2$ at $G$ which is the other diametric end. The coordinates of $G$ are $\Large(\frac{-2at^2}{t^2+2}, \frac{at^3}{t^2+2})$.Now the midpoint of $GP$ is the center $L$ which comes out to be,
$$\large x = \frac{at^4}{2(t^2+2)}$$
$$\large y = \frac{3at^3+4at}{2(t^2+2)}$$ 
Now, I am not able to eliminate $t$ from these equations. How do I proceed?  
Note: This problem is from SL Loney Co-ordinate geometry which was first printed a 100 years ago. So, there surely must be a shorter solution. 
 A: We are given
\begin{gather}
\tag{1}\label{eq:1}
x = \frac{at^4}{2(t^2+2)}, \\
\tag{2}\label{eq:2}
y = \frac{3at^3+4at}{2(t^2+2)}.
\end{gather}
Combining \eqref{eq:1} and \eqref{eq:2},
$$
t\frac{y}{a} - \frac{x}{a} = \frac{2t^4 + 4t^2}{2(t^2+2)} = t^2,
$$
whence
\begin{equation}
\tag{3}\label{eq:3}
\left(t^2 + \frac{x}{a}\right)^2 = t^2\frac{y^2}{a^2}.
\end{equation}
Rewriting \eqref{eq:1},
$$
t^4 - 2\frac{x}{a}t^2 - 4\frac{x}{a} = 0,
$$
whence
\begin{equation}
\tag{4}\label{eq:4}
\left(t^2 - \frac{x}{a}\right)^2 = \frac{x^2}{a^2} + 4\frac{x}{a}.
\end{equation}
From \eqref{eq:3} and \eqref{eq:4},
$$
4t^2\frac{x}{a} = t^2\frac{y^2}{a^2} - \frac{x^2}{a^2} - 4\frac{x}{a},
$$
whence
\begin{equation}
\tag{5}\label{eq:5}
t^2\left(\frac{y^2}{a^2} - 4\frac{x}{a}\right) = \frac{x^2}{a^2} + 4\frac{x}{a}.
\end{equation}
With hindsight, this deduction could have been shortened, but I've preferred to keep it in a sequence that occurred fairly naturally. Anyway, we can now substitute \eqref{eq:5} in \eqref{eq:4}. How to do so, while making the least mess, is debatable. I decided it was about time to get rid of the denominators, thus
$$
(at^2 - x)^2 = x(x + 4a) = t^2(y^2 - 4ax),
$$
leading to this manageable mess,
\begin{align*}
4ax(y^2 - 4ax)^2 & = a^2[t^2(y^2 - 4ax)]^2 - 2ax(y^2 - 4ax)[t^2(y^2 - 4ax)] \\
& = a^2[x(x + 4a)]^2 - 2ax(y^2 - 4ax)[x(x + 4a)],
\end{align*}
whence
$$
\boxed{4(y^2 - 4ax)^2 = x(x + 4a)[a(x + 4a) - 2(y^2 - 4ax)]}
$$
Gathering all terms containing $y$ on the left hand side,
$$
4y^2(y^2 - 8ax) + 2x(x + 4a)y^2 = ax(x + 4a)[x + 4a + 8x] - 4(4ax)^2.
$$
Simplifying,
\begin{align*}
2y^2(2y^2 - 16ax + x^2 + 4ax) & = ax[(x + 4a)(9x + 4a) - 64ax] \\
& = ax(9x^2 - 24ax + 16a^2),
\end{align*}
and finally,
$$
\boxed{2y^2(2y^2 + x^2 - 12ax) = ax(3x - 4a)^2}
$$
A: Let us find the slope of the tangent:
$$y^2=4ax$$
$$2yy'=4a$$
$$y'=\frac{2a}y$$
Distances $OL$ and $LP$ are equal:
$$x_L^2+y_L^2=(x_P-x_L)^2+(y_P-y_L)^2\tag{1}$$
The point $P$ lies on the parabola:
$$y_P^2=4ax_P\tag{2}$$
Line $LP$ is tangent to parabola:
$$\frac{y_P-y_L}{x_P-x_L}=y'_P=\frac{2a}{y_P}\tag{3}$$
Introduce the following substitutions:
$$x_P-x_L=u\tag{4}$$
$$y_p-y_L=v\tag{5}$$
Replace (4) and (5) into (2) and (3):
$$(v+y_L)^2=4a(u+x_L)\tag{6}$$
$$\frac{v}{u}=\frac{2a}{v+y_L}\tag{7}$$
From (6) and (7) you get:
$$v^2=y_L^2-4ax_L\tag{8}$$
$$u=\frac{1}{2a}v(v+y_L)\tag{9}$$
Rearrange (1):
$$x_L^2+y_L^2=u^2+v^2=\frac{1}{4a^2}v^2(v+y_L)^2+v^2$$
$$4a^2\left(\frac{x_L^2+y_L^2}{v^2}-1\right)=(v+y_L)^2$$
$$\frac{4a^2}{v^2}(x_L^2+y_L^2-v^2)-v^2-y_L^2=2vy_L$$
$$\left(\frac{4a^2}{v^2}(x_L^2+y_L^2-v^2)-v^2-y_L^2\right)^2=4v^2y_L^2$$
Now replace (8) into (10) and the final result is the locus of points $L(x_L,y_L)$.
Hopefully I did not make any mistake along the way. Even if I did, I think that I have demonstrated the right way.
A: Simplify slightly by substituting $x = au/2$ and $y= av/2$. Then your parameterization yields
$$u(t^2+2)=t^4 \qquad v(t^2+2)=t(3t^2+4) \tag{1}$$
Squaring the second equation gives us a system in even powers of $t$, so define $s := t^2$, giving
$$u(s+2) = s^2 \qquad v^2(s+2)^2=s(3s+4)^2 \tag{2}$$
We can use the $u$ equation to whittle-down powers of $s$ in the $v$ equation. It'll be tedious, but completely mechanical. At least we can look forward to the fact that the end result will be linear in $s$.
$$\begin{align}
v^2(\;s^2+4s+4\;)&=s(\;9s^2+24s+16\;) \tag{3}\\
v^2(\;u(s+2)+4s+4\;)&=s(\;9\cdot u(s+2)+24s+16\;) \tag{4}\\
v^2(\;s(u+4)+2(u+2)\;)&=3s^2(3u+8)+2s(9u+8) \tag{5}\\
v^2(\;s(u+4)+2(u+2)\;)&=3(3u+8)\cdot u(s+2)+2s(9u+8) \tag{6}\\
s(\;v^2(u+4)-9u^2-42u-16\;) &= 6u(3u+8)-2v^2(u+2) \tag{7}
\end{align}$$
From here, we ---and by "we", I mean "the reader"--- can solve $(7)$ for $s$ and substitute back into the $u$ equation in $(2)$. Once the dust settles, "we" will have the $u$-$v$ equivalent of the target relation:
$$v^2 (\;2v^2+u^2-24 u\;) = u (\;3u-8\;)^2 \tag{$\star$}$$
