Envelope of Projectile Trajectories

Yes it is. To find the envelope, we just have to find the intersections between two trajectories associated to two slightly different angles. If we solve $$ x\tan\theta -\frac{gx^2}{2v^2}\sec^2(\theta) = x\tan(\theta+\varepsilon) -\frac{gx^2}{2v^2}\sec^2(\theta+\varepsilon) $$ we get $x=0$ or $$ x=\frac{2v^2}{g}\cdot\frac{\tan(\theta)-\tan(\theta+\varepsilon)}{\sec^2(\theta)-\sec^2(\theta+\varepsilon)}$$ and by letting $\varepsilon\to 0$ we get $x=\frac{v^2}{g}\cot(\theta)$, from which $y=\frac{v^2}{2g}\left(2-\frac{1}{\sin^2(\theta)}\right) $.

It follows that the equation of the envelope is given by: $$ y = \frac{v^2}{2g}\left(1-\left(\frac{gx}{v^2}\right)^2\right)=\frac{v^2}{2g}-\frac{g}{2 v^2}\,x^2$$ that clearly is a parabola with vertex in $\left(0,\frac{v^2}{2g}\right)$ through the points $\left(\pm\frac{v^2}{g} ,0\right)$.

Geometric sketch of the projectile motion

We may notice that the envelope and the trajectory with $\theta=\frac{\pi}{4}$ are homothetic, and the dilation ratio is just $2$. The vertices of the trajectories lie on an ellipse that is tangent to the envelope parabola, with centre at $\left(0,\frac{v^2}{4g}\right)$, a vertex in the origin and a vertex at $\left(\frac{v^2}{2g},\frac{v^2}{4g}\right)$.


@Jack provides a very nice and intuitive derivation of the envelope as the points of intersection of infinitely-close members of the curve family. The Wikipedia "Envelope" entry provides this less-illuminating abstraction:

The envelope of the family [of curves parameterized by $t$ is] the set of points for which $$F(t, x, y) = \frac{\partial F}{\partial t}(t,x,y) = 0 \tag{$\star$}$$ for some value of $t$ [...].

In $(\star)$, $F$ is the function that, when set equal to $0$, defines each curve in the family. For the question at hand, we have (with parameter $\theta$ instead of $t$) $$F(\theta,x,y) = -y + x\tan\theta -\frac{g x^2}{2v^2}\sec^2\theta \tag{1}$$ Therefore, differentiating with respect to $\theta$ gives $$\frac{\partial F}{\partial \theta}(\theta,x,y) =x\sec^2\theta -\frac{g x^2}{2v^2}\cdot 2 \sec^2\theta\tan\theta = \frac{x\sec^2\theta}{v^2} ( v^2 - g x \tan\theta) \tag{2}$$ Solving $\partial F/\partial \theta = 0$ for $\theta$ (noting that $\sec\theta$ never vanishes) gives $$\tan\theta = \frac{v^2}{g x} \qquad\text{so that}\qquad \sec^2\theta = 1 + \tan^2\theta = \frac{g^2 x^2 + v^4}{g^2 x^2}$$

Substituting into $(1)$, and setting $F=0$, we have

$$y = x\;\frac{v^2}{gx} - \frac{g x^2}{2 v^2}\;\frac{g^2 x^2 + v^4}{g^2 x^2} = \frac{v^2}{2g} - \frac{g x^2}{2v^2} = \frac{v^4-g^2x^2}{2gv^2}$$

which agrees with Jack's answer.