I am trying to find the equation of a parabola given the following conditions.
Given:
- circle with radius $R_1$ centered at $(0, R_1 + R_2)$
- $R_2$ is given
- $\frac{R_1}{R_2} = 0.382$
- $R_3 = R_2 * \sqrt{16}$
- $L$ is a given constant which is the length of the conic section from the vertex to points $D$ and $E$
- parabola does not have a vertex at the origin
- parabola is tangent to circle at some point to the right of the x-axis, '$A$'
- parabola has points: $B(x_A, y_A)$, $C(x_A, -y_A)$, $D(x_D, R_3)$, $E(x_, -R_3)$
- points B and C are symmetric wrt. $x$ axis
- points D and E are symmetric wrt. $x$ axis
My approach:
Use the general equation of a horizontal conic section: $x = ay^2 + by + c$
By substituting in points B and C or D and E, I find that: $b = 0$. Therefore, $$ x = a y^2 + c $$
I want $\frac{dy}{dx}$, so through implicit differentiation I get: $$dx = 2a*ydy$$ $$\frac{dy}{dx} = \frac{1}{2ay}$$
Next, I attempted to express $x_n$ and $y_n$ in terms of a common unknown. $$x_n = R_1 * \sin{\theta_N}$$ $$y_n = R_2 + (R_1 - R_1 * \cos{\theta_N})$$
Substituting $y_n$ into the equation for $\frac{dy}{dx}$, $$\frac{dy}{dx} = \tan{\theta_N} = \frac{1}{2a*(R_2 + (R_1 - R_1\cos{\theta_N}))}$$
Substituting $y_n$ and $x_n$ into the equation of a parabola yields, $$R_1 * \sin{\theta_N} = a*(R_2 + (R_1 - R_1 * \cos{\theta_N}))^2 + c$$
Substituting $y_n$ into the equation for $\frac{dy}{dx}$, $$\frac{dy}{dx} = \tan{\theta_N} = \frac{1}{2a*(R_2 + (R_1 - R_1\cos{\theta_N}))}$$
Since $\frac{R_2}{R_1} = 0.382$, $R_2 = 2.6178R_1$ $$\frac{dy}{dx} = \tan{\theta_N} = \frac{1}{2a*(3.6178R_1 - R_1\cos{\theta_N})}$$ $$R_1 * \sin{\theta_N} = a*(3.6178R_1 - R_1 * \cos{\theta_N})^2 + c$$
At this point, I am not sure how to proceed.
Context:
This problem is used in finding the diverging section of a bell nozzle for rocket engines. I am attempting to find the equations of the plots on page 80 of Rocket Propulsion Elements 7th edition by Sutton and Biblarz.
A more general form of this problem would have left the $16$ in the equation $R_3 = R_2 * \sqrt{16}$ in variable form. This would have allowed the plots in RPE to be generated. $16$ in this context represents the area ratio between the exit plane and the throat.