On the dynamics of two connected objects located at different circular orbits

Question

I'm just thinking about the qualitative and quantitative implications of a hypothetical orbital problem, but what I ended up does not make sense as it should.

In particular, assume that two objects are placed at two close circular orbits, and they are also connected to each other by a rigid tether-like stuff, whose length is $l$, as the following sketch depicts.

I assume that the distance of the upper (resp. lower) object from the center of the Earth is $r'$ (resp. $r$). What I am interested is the study of the overall stability of this system in view of the overtake rate of the object on the lower orbit (we know that the lower-orbit object is faster, thereby overtaking its peer over the passage of time).

Here is my analysis:

The overtake can be written as

$\Delta x = (v_{d}-v_{u})t.$

Given the gravitational parameter $\mu$, we have

$v_d = \sqrt{\dfrac{\mu}{r}},$

$v_u = \sqrt{\dfrac{\mu}{r + l \sin \theta}},$

which end up with

$\Delta x = \sqrt{\mu} (\dfrac{1}{\sqrt{r}} - \dfrac{1}{\sqrt{r + l \sin \theta}})t$.

I'm particularly interested in the $\dot{\Delta x} = 0$. Thus given $\theta = \theta(t)$, the total differentiation role reads

$\dot{\Delta x} = \dfrac{\partial \Delta x}{\partial t}\dfrac{\partial t}{\partial t}+ \dfrac{\partial \Delta x}{\partial \theta}\dfrac{\partial \theta}{\partial t}.$

In particular, one obtains

$\dfrac{\partial \Delta x}{\partial t} = \sqrt{\mu} (\dfrac{1}{\sqrt{r}} - \dfrac{1}{\sqrt{r + l \sin \theta}}),$

and

$\dfrac{\partial \Delta x}{\partial \theta} = \sqrt{\mu}(\dfrac{1}{2}l\dot{\theta}\cos \theta (r + l \sin \theta)^{-\dfrac{3}{2}}) t$.

Thus, applying $\dot{\Delta x} = 0$ condition yields the following dynamics corresponding to $\theta$.

$\dot{\theta} = \dfrac{2(r + l \sin \theta)^{\dfrac{3}{2}})[\sqrt{\dfrac{1}{r + l \sin \theta}}-\sqrt{\dfrac{1}{r}}]}{lt\cos \theta}$

Now, time for numerically integrating the equation above. In particular, the following Python script

import numpy as np
import math
from scipy.integrate import odeint
import matplotlib.pyplot as plt

def f(s,t):
    l = 20
    h = s[0]
    r = 1000
    dhdt = 2*((l*math.sin(h)+r)**(1.5))*((1/math.sqrt(l*math.sin(h)+r))-(1/r))/(t*l*math.cos(h))
    return dhdt

t = np.linspace(0.1,2000)
s0=[20]

s = odeint(f,s0,t)

plt.plot(t,s[:,0],'r--', linewidth=2.0)

plt.show()

generates

The dynamics look reasonable since one expects that $\theta$ eventually converges to the equilibrium $\theta = 0$. However, that $t$, jumping up and down in the denominator of $\dot{\theta}$, looks pretty weird because it literally implies that the rate of $\theta$ is infinitely large when the system is initiated! Did I mess during my computations?!

Any comments to clarify the situation are highly appreciated.

Answer 1

I'm afraid that your whole approach is wrong due to starting with contradictory assumptions. If the two bodies are connected by a rigid tether, then in general, they will not move in circular orbits, since in addition to Earth's gravity, forces from the tether will act on them. (As uhoh pointed out, there are some configuration where the bodies do move in circles.)

$\Delta x = (v_{d}-v_{u})t$ would be correct only if $v_d - v_u$ stayed constant, but even in your post, this difference changes with time, since it depends on $\theta$, which changes with time. $v_d - v_u$ is the momentary rate of change of $\Delta x$, so what you should have instead is $\dot{\Delta x} = v_d − v_u$.

The formulas for $v_d$ and $v_u$ are not correct because the bodies would not in general move in circular orbits.

And if $\dot{\Delta x} = 0$ over an interval of time, then $\dot{\theta}=0$ as well, since $\Delta x = l\cos\theta$.