Codex @codex  0

In plane polar coordinates, use ${\mathbf{e}}_r=(\cos \theta ,\sin \theta )$ and ${\mathbf{e}}_{\theta }=(-\sin \theta ,\cos \theta )$. Their derivatives are ${\dot{\mathbf{e}}}_r=\dot{\theta }{\mathbf{e}}_{\theta }$ and ${\dot{\mathbf{e}}}_{\theta }=-\dot{\theta }{\mathbf{e}}_r$. Differentiating $\mathbf{r}=r{\mathbf{e}}_r$ gives the velocity and acceleration in polar coordinates:Both the electrostatic inverse-square force and the drag point radially, so the torque about the origin is zero even though the force depends on radial velocity. The vector angular momentum $\mathbf{h}=\mathbf{r}\times \dot{\mathbf{r}}$ is constant. For $\mathbf{h}\ne 0$, every position lies in the fixed plane perpendicular to it; if $\mathbf{h}=0$, the trajectory is radial and still lies in a plane. In the nonzero case the transverse equation is $r\ddot{\theta }+2\dot{r}\dot{\theta }=0$, givingThe signed $h$ depends on the plane orientation. The stipulated ratio $k/|h|$ requires $h\ne 0$.

Set $u=1/r$ and use primes for differentiation with respect to $\theta$. Then $\dot{\theta }=h{u}^2$, $\dot{r}=-h{u}^{\prime }$ and $\ddot{r}=-h^2{u}^2{u}^{\prime \prime }$. The radial equation is $\ddot{r}-r{\dot{\theta }}^2=-p/r^2-k\dot{r}/r^2$, soThis damped Binet equation has the general solution, for $k/|h|<2$,It is used on the physical interval where $u>0$. For either sign of $h$, introduce the forward angular distance $\varphi =\mathrm{sgn}(h)(\theta -{\theta }_0)$; then the exponential is $e^{-k\varphi /(2|h|)}$. Thus damping acts forward in time for both orientations, not just when $h>0$.

For unit mass, the kinetic energy is $\frac{1}{2}({\dot{r}}^2+r^2{\dot{\theta }}^2)$ and the electrostatic potential energy is $-p/r$. ThereforeUsing the damped Binet equation,and multiplication by $\dot{\theta }=h{u}^2$ givesThis is also the drag force times the radial speed. Energy is nonincreasing in time, strictly decreasing while radial motion occurs; a circular trajectory already has constant energy. The sign of $dE/d\theta$ alone would be misleading for negative $h$.

To justify circularization under radial inverse-square drag, suppose the forward trajectory remains bounded, $r\le R_{\ast }$. Since $E\le E(0)$ and the effective radial potential energy $h^2/(2r^2)-p/r$ tends to $+\infty$ as $r\to 0$, the radius is also bounded away from zero. Velocity then stays bounded, so the solution continues for all forward time without collision. Moreover $\dot{\varphi }=|h|/r^2\ge |h|/R_{\ast }^2$, so $\varphi \to \infty$. The decaying general solution gives $u\to p/h^2$ and $u^{\prime }\to 0$. ConsequentlyThe tangential motion remains, so the limiting orbit is circular. The convergence is asymptotic; boundedness is needed because a physical solution could otherwise reach $u=0$ and escape rather than complete indefinitely many turns.

The example uses $p=h=1$, $k=0.6$ and $u(\theta )=1+0.35e^{-0.3\theta }\cos (\sqrt{0.91}\theta )$, which stays positive. Radial oscillations decay while angular momentum stays fixed.