October 30th, 2018

Гравитационный потенциал тонкого диска

Предыдущий пост: https://dims12.livejournal.com/365433.html

Задачи: https://nplus1.ru/news/2018/10/24/physics

Исходное обсуждение: https://polenadisto.livejournal.com/768774.html

Программа для расчёта находится здесь: https://github.com/dims12/DiskPotential

Я вычислил гравитационный потенциал тонкого диска и он получился таким



Видно, что я был неправ, предполагая, что ближайшая к диску эквипотенциальная поверхность будет прямоугольной. В то же время, видно, что эквипотенциальные поверхности наклонены к диску под очень небольшим углом и наклоняются под 45 градусов только на 70-80 единицах от центра диска. Иными словами, на значительной площади диска поле будет однородным и мало отличимым от земного.

Программа написана на Питоне в двух версиях: для Numpy и для Tensorflow. Обе версии эквивалентны друг другу, но вторая версия при наличии GPU работает в десятки раз быстрее.