Skip to main content
Physics LibreTexts

9.6: Position in a Parabolic Orbit

  • Page ID
    8183
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    When a “long-period” comet comes in from the Oort belt, it typically comes in on a highly eccentric orbit, of which we can observe only a very short arc. Consequently, it is often impossible to determine the period or semi major axis with any degree of reliability or to distinguish the orbit from a parabola. There is therefore frequent occasion to have to understand the dynamics of a parabolic orbit.

    We have no mean or eccentric anomalies. We must try to get \(v\) directly as a function of \(t\) without going through these intermediaries.

    The angular momentum per unit mass is given by Equation 9.5.28a:

    \[h = r^2 \dot{v} = \sqrt{2G \textbf{M} q}, \label{9.7.1} \tag{9.7.1}\]

    where \(v\) is the true anomaly and \(q\) is the perihelion distance.

    But the Equation to the parabola (see Equation 2.4.16) is

    \[r = \frac{2q}{1 + \cos v}, \label{9.7.2} \tag{9.7.2}\]

    or (see section 3.8 of Chapter 3), by making use of the identity

    \[\cos v = \frac{1-u^2}{1+u^2} , \quad \text{where} \quad u = \tan \frac{1}{2} v , \label{9.7.3a,b} \tag{9.7.3a,b}\]

    the Equation to the parabola can be written

    \[ r = q \sec^2 \frac{1}{2} v . \label{9.7.4} \tag{9.7.4}\]

    Thus, by substitution of Equation \ref{9.7.4} into \ref{9.7.1} and integrating, we obtain

    \[q^2 \int_0^v \sec^4 (\frac{1}{2} v) dv = \sqrt{2G\textbf{M} q} \int_T^t dt. \label{9.7.5} \tag{9.7.5}\]

    Upon integration (drop me an email if you get stuck!) this becomes

    \[u + \frac{1}{3} u^3 = \frac{\sqrt{\frac{1}{2}G \textbf{M}}}{q^{3/2}} (t - T) . \label{9.7.6} \tag{9.7.6}\]

    This Equation, when solved for \(u\) (which, remember, is \(\tan \frac{1}{2} v\)), gives us \(v\) as a function of \(t\). As explained at the end of section 9.5, if \(q\) is in astronomical units and \(t − T\) is in sidereal years, and if the mass of the comet is negligible compared with the mass of the Sun, this becomes

    \[u + \frac{1}{3} u^3 = \frac{π \sqrt{2}(t - T)}{q^{3/2}} \label{9.7.7} \tag{9.7.7}\]

    or \[3u + u^3 - C = 0, \quad \text{where} \quad C = \frac{π \sqrt{18} (t-T)}{q^{3/2}}. \label{9.7.8a,b} \tag{9.7.8a,b}\]

    There is a choice of methods available for solving Equation \ref{9.7.8a,b}, so it might be that the only difficulty is to decide which of the several methods you want to use! The constant \(\frac{1}{3}C\) is sometimes called the “parabolic mean anomaly”.

    Method 1: Just solve it by Newton-Raphson iteration. Thus \(f = 3u + u^3 − C = 0\) and \(f^\prime = 3(1 + u^2 )\), so that the Newton-Raphson \(u = u − f / f^\prime\) becomes

    \[u = \frac{2u^3 + C}{3(1 + u^2)}, \label{9.7.9} \tag{9.7.9}\]

    which should converge quickly. For economy, calculate \(u^2\) only once per iteration.

    Method 2:

    Let \[u = x-1/x \quad \text{and} \quad C = c-1/c. \label{9.7.10a,b} \tag{9.7.10a,b}\]
    Then Equation 9.7.8a becomes

    \[x = c^{1/3} . \label{9.7.11} \tag{9.7.11}\]

    Thus, as soon as \(c\) is found, \(x\), \(u\) and \(v\) can be calculated from Equations 9.7.11, 10a, and 3a or b, and the problem is finished – as soon as \(c\) is found!

    So, how do we find c? We have to solve Equation 9.7.10b.

    Method 2a:

    Equation 9.7.10b can be written as a quadratic Equation:

    \[c^2 - Cc - 1 = 0. \label{9.7.12} \tag{9.7.12}\]

    Just be careful that you choose the correct root; you should end with \(v\) having the same sign as \(t − T\).

    Method 2b:

    Let \[ C = 2 \cot 2 \phi \label{9.7.13} \tag{9.7.13}\]

    and calculate \(\phi\). But by a trigonometric identity,

    \[2 \cot 2 \phi = \cot \phi - 1/\cot \phi \label{9.7.14} \tag{9.7.14}\]

    so that, by comparison with Equation 9.7.10b, we see that

    \[c = \cot \phi . \label{9.7.15} \tag{9.7.15}\]

    Again, just make sure that you choose the right quadrant in calculating \(\phi\) from Equation \ref{9.7.13}, so as to be sure that you end with \(v\) having the same sign as \(t − T\).

    Method 3.

    I am told that Equation 9.7.8 has the exact analytic solution

    \[u = \frac{1}{2} w^{\frac{1}{3}} - 2w^{-\frac{1}{3}}, \label{9.7.16} \tag{9.7.16}\]

    where \[w = 4C + \sqrt{64 + 16C^2}. \label{9.7.17} \tag{9.7.17}\]

    I haven’t verified this for myself, so you might like to have a go.

    Example: Solve the Equation \(3u + u^3 = 1.6\) by all four methods. (Methods 1, 2a, 2b and 3.)

    Example: A comet is moving in a parabolic orbit with perihelion distance \(0.9 \ \text{AU}\). Calculate the true anomaly and heliocentric distance 20 days after perihelion passage. (A sidereal year is 365.25636 days.)

    Exercise: Write a computer program that will return the true anomaly as a function of time, given the perihelion distance of a parabolic orbit. Test it with your answer for the previous example.


    This page titled 9.6: Position in a Parabolic Orbit is shared under a CC BY-NC 4.0 license and was authored, remixed, and/or curated by Jeremy Tatum via source content that was edited to the style and standards of the LibreTexts platform.