Simulaciones del centro galáctico

Uso de simulaciones hidrodinámicas para comprender la entrada de gas en los ~200 pc internos de la Vía Láctea provocada por la fusión con una galaxia satélite.

CÓDIGO: órbitas periódicas en potenciales no axisimétricos

Código C++ usado en Gallego & Cuadra 2017 para estimar órbitas periódicas a 200 pc del centro galáctico, donde domina el potencial de la barra. Estas órbitas se usaron como condiciones iniciales para una simulación hidrodinámica con Gadget-2. El método para identificar órbitas periódicas se llama “técnica de disparo” (shooting technique). Parte de una suposición para la velocidad inicial en el espacio de fases x−y, en el que la componente x de la velocidad es cero sobre el semieje mayor de la órbita. Luego, las ecuaciones de movimiento se integran con un método de Runge–Kutta de cuarto orden en el sistema de referencia de la barra en rotación. La suposición de la velocidad inicial se optimiza hasta obtener una órbita periódica con la precisión relativa deseada. Descarga el código aquí: find_periodic_orbits.cpp