# Spacecraft path with polar coordinates

There is a circular gate rotating at a constant angular speed of ω. The circular gate has a tunnel across its diameter. The mission is to pass through the gate. (That is, come in one side of the gate, travel the whole diameter, and exit at the other side.)

Also, craft is neutrally buoyant, the craft has four thrusters (in the i, -i, j, and -j directions), and sudden movements cause "slosh" in the fuel tank that destroys the craft.

My thought was to describe the position of the craft in the complex plane $$z=r \cdot e^{i\theta}$$ where r and θ depend on time t. Let's say at time 0 the craft is in front of the tunnel opening.

Differentiating once, I get $$v = (\dot r +i r \dot \theta)e^{i\theta}$$
which describes a velocity of $$\dot r$$ in the radial direction and a velocity of $$r \dot \theta$$ in the perpendicular direction.

Differentiating again, I get $$a = (\ddot r - r {\dot \theta}^2 )+i(2\dot r \dot \theta + r \ddot \theta)] e^{i\theta}$$ which has the familiar expressions for acceleration in the radial and perpendicular directions.

Now, $$\dot \theta = \omega$$ which is constant. We also want to move smoothly, so we want the radial speed to be constant as well. Thus, our acceleration is $$a = [(-r \omega^2) + i(2\dot r \dot \theta + r \ddot \theta)]e^{i\theta}.$$

Now, all that matters is that the craft gets through the tunnel, so I set $$\dot r \equiv -1$$ until the craft gets to the center and then $$\dot r \equiv 1$$ after we get to the center.

Thus, the acceleration before getting to the center is $$a = [(-r {\dot \theta}^2)+i(-2\omega)]e^{i\theta}.$$

Since the thrusters are in the i and j directions, I multiply the expression out, substitute cosθ = x/r, and sinθ = y/r, and r = √x^2+y^2 to obtain $$\left(-\omega^2 \cdot x+2\omega \cdot \frac{y}{\sqrt{x^2+y^2}}\right)+i \left(-\omega^2 \cdot y - 2\omega \frac {x}{\sqrt{x^2+y^2}} \right)$$.

Of course, the real part represents the horizontal acceleration and the imaginary part gives the vertical.

The thrust would then be given by the mass times the acceleration.

This is apparently incorrect, since I have run it on a simulation. Does anyone see where I went wrong?

Last edited: