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.