If your convex set $C$ is a linear subspace, then you can take some orthonormal basis of it, say $e_i, i=1,..,k\leq n$ and your projection must satisfy $x-y\perp e_i, i=1,2,..,k$. Therefore, if $y$ is represented in the mentioned basis as $y=\sum\limits_{i=1}^{k}{a_ie_i}$, then we must have $(x,e_i)-a_i=0$, i.e $a_i=(x,e_i)\Rightarrow y=\sum\limits_{i=1}^{k}{(x,e_i)e_i}$.
If $C$ is just a convex set, then I guess an algorithmic approach would be to draw "circles" with center $x$ and a radius equal to $\|x-y1\|$ for some $y_1\in C$. Then start decreasing the radius until you do not have an intersection with $C$. So this is a mediocre method to locate the projection $y$.