Я бы сделал примерно такую систему уравнений.
r = sqrt(t)
phi = a·r
t — параметр, условное «время»; phi — полярный угол, r — длина радиус-вектора.
Ну и, соотвественно, x = r cos phi, y = r sin phi.
В общем, радиус (ну или угол) должен увеличиваться со скоростью квадратного корня.
В этом деле есть физический смысл — это решение дифура r′(t)=1/r. Только двоечку и константу интегрирования упустил, ибо они нам как бы не нужны. Метод не точный, но если посмотреть на длину дуги спирали, там самый большой член квадратичный.
Если нужен СОВСЕМ стабильный зазор (например, расположить по спирали какие-то кружочки), у меня есть рекуррентный алгоритм.
Как написать алгоритм спирали?