Решения геодезических уравнений в пространстве AdS3

Я борюсь с заданием от моего лектора, он попросил меня численно (математика) и аналитически вычислить решения для траектории точечной частицы в пространстве AdS3 с глобальными координатами, заданными метрикой

г с 2 "=" р 2 ( чушь 2 ( р ) г т 2 + г р 2 + грех 2 ( р ) г ф 2 )
Я рассчитал три геодезических уравнения движения

р ¨ + грех ( р ) чушь ( р ) ( т ˙ т ˙ ф ˙ ф ˙ ) "=" 0

ф ¨ + ткань ( р ) р ˙ ф ˙ "=" 0

т ¨ + танх ( р ) р ˙ т ˙ "=" 0

Я понимаю, что они связаны, но я, честно говоря, понятия не имею, как я найду р ( т ) , т ( т ) и ф ( т ) .

Во-первых, я думаю, что у вас могут отсутствовать некоторые множители 2 в ф и т уравнения.
Во-вторых, легко разрешимый случай, когда р фиксированный.
Привет Секавара, не могли бы вы уточнить?
О, я сделал несколько каракулей, так что могу ошибаться, но для ф и т я получил ф ¨ + 2 ткань ( р ) р ˙ ф ˙ "=" 0 и т ¨ + 2 танх ( р ) р ˙ т ˙ "=" 0 . И затем заметьте, что если мы возьмем р ˙ "=" 0 , то мы можем легко найти некоторые геодезические.
В таком случае, т "=" а + б т и ф "=" с ± б т .

Ответы (1)

Предполагая, что ваши уравнения верны (я не проверял), такие связанные ОДУ (и геодезические уравнения в целом) должны решаться численно. Хитрость заключается в использовании стандартных решателей Рунге-Кутты, наиболее удобно иметь дело с системами первого порядка. В общем, учитывая геодезические уравнения:

Икс ¨ а + Г б с а Икс ˙ б Икс ˙ с "=" 0 ,

мы можем записать их как систему первого порядка:

Икс ˙ а "=" в а , в ˙ а "=" Г б с а в б в с

Итак, для вашей системы выше определите Икс а "=" р ˙ а , у а "=" ф ˙ а , и г а "=" т ˙ а . Тогда в сочетании с оставшимися уравнениями у вас есть 6-мерная система первого порядка, и решения будут зависеть от задания 6 таких начальных условий.

Надеюсь это поможет.

Использование решателей RK вряд ли можно назвать «хитростью», так как они являются одними из самых популярных и часто являются «переходными», прежде чем прибегать к более сложным методам.
@JamalS Я думаю, ты слишком буквально воспринял мое предложение. Решатели RK изначально не работают с ОДУ 2-го порядка, «хитрость» заключается в том, чтобы написать ОДУ 2-го порядка как ОДУ первого порядка, а затем применить методы РК.
@Dr.IkjyotSinghKohli Я добился значительного прогресса в своей проблеме, однако у меня возникли проблемы с указанием начальных условий, которые приводят к стабильному решению в соответствии с используемым решателем. Спасибо
@Cyphox32 привет. Да. Это часто является проблемой в таких числовых приложениях. Возможно, вы захотите взглянуть на использование MATLAB/Octave/NumPY/SciPy в Python для таких вещей. В MATLAB есть очень мощные численные решатели: ODe23s, OdE45 и т. д., и с ними легче экспериментировать с начальными условиями.