miniaturadeimagen| Una simulación de n cuerpos de la formación cosmológica de un cúmulo de galaxias en un universo en expansión. En física y astronomía, una simulación de n cuerpos es una simulación de un sistema dinámico de partículas, generalmente bajo la influencia de fuerzas físicas, como la gravedad (ver problema de los n cuerpos para otras aplicaciones). Las simulaciones de n cuerpos son herramientas ampliamente utilizadas en astrofísica, desde la investigación de la dinámica de sistemas de pocos cuerpos como el sistema Tierra-Luna-Sol hasta la comprensión de la evolución de la estructura a gran escala del universo. En cosmología física, las simulaciones de n cuerpos se utilizan para estudiar procesos de formación de estructuras no lineales, como filamentos y halos de galaxias, debido a la influencia de la materia oscura. También se utilizan simulaciones directas de n cuerpos para estudiar la evolución dinámica de los cúmulos estelares.
Naturaleza de las partículas Las 'partículas' tratadas por la simulación pueden corresponder o no a objetos físicos de naturaleza definida por partículas. Por ejemplo, una simulación de n cuerpos de un cúmulo de estrellas podría estar compuesto de n partículas, por lo que cada partícula representa una estrella y tiene un significado físico. Por otro lado, en el caso de una simulación de una nube de gas, resulta imposible representar cada átomo o molécula de gas por una partícula, ya que esto requeriría del orden de 1023 partículas por cada mol de material (ver constante de Avogadro). Debido a limitaciones de la computación actual, una sola 'partícula' debe representar una cantidad mucho mayor de gas (a menudo implementado utilizando Hidrodinámica de Partículas Suavizadas). Esta cantidad no necesita tener ningún significado físico, sino que debe elegirse como un compromiso entre precisión y requisitos informáticos manejables.
Simulación de Materia Oscura La materia oscura toma un papel importante en la formación de galaxias. La ecuación de Boltzmann sin colisiones describe la evolución de la densidad f de las partículas de materia oscura respecto al tiempo (en el espacio de fases):
En esta ecuación, es la velocidad y Φ es el potencial gravitacional obtenido por la ecuación de Poisson. Estas dos ecuaciones juntas se resuelven en un universo que se rige por las ecuaciones de Friedmann tras determinar las condiciones iniciales de las partículas de materia oscura. El método convencional para inicializar posiciones y velocidades implica mover las partículas dentro de una red cartesiana uniforme o una configuración de partículas similar al vidrio. Esto se puede conseguir mediante una aproximación de teoría lineal o una teoría de perturbaciones de bajo orden.
Simulaciones de n cuerpos gravitacionales directos miniaturadeimagen| Simulación de n cuerpos de 400 objetos con parámetros similares a los de los planetas del sistema solar . En las simulaciones gravitacionales directas de n cuerpos, estos se encuentran bajo la mutua influencia de sus fuerzas gravitacionales. Las respectivas ecuaciones se integran numéricamente sin aproximaciones ni simplificaciones. Estos cálculos se utilizan en casos donde las interacciones entre objetos individuales, como estrellas o planetas, son importantes para la evolución del sistema.
Las primeras simulaciones gravitacionales directas de n cuerpos fueron llevadas a cabo por Erik Holmberg en el Observatorio de Lund en 1941, determinando las fuerzas entre estrellas en galaxias influidas mutuamente por su gravitación mediante la equivalencia matemática entre la propagación de la luz y la interacción gravitacional. Colocando bombillas en las posiciones de las estrellas y midiendo los flujos de luz direccionales con una célula fotoeléctrica, las ecuaciones de movimiento se pueden integrar con una complejidad de O(N). Las primeras simulaciones puramente por cálculos fueron realizadas por Sebastian von Hoerner en el Astronomisches Rechen-Institut de Heidelberg, Alemania. Sverre Aarseth de la Universidad de Cambridge (Reino Unido) ha dedicado toda su vida científica al desarrollo de una serie de códigos de n cuerpos altamente eficientes para aplicaciones astrofísicas que utilizan pasos de tiempo adaptativos (jerárquicos), un esquema vecino de Ahmad-Cohen y la regularización de encuentros cercanos. La regularización es un truco matemático para eliminar la singularidad en la ley de gravitación newtoniana de dos partículas que se acercan arbitrariamente entre sí. Los códigos de Sverre Aarseth se utilizan para estudiar la dinámica de cúmulos estelares, sistemas planetarios y núcleos galácticos.
Simulaciones de relatividad general Muchas simulaciones son lo suficientemente grandes como para que los efectos de la relatividad general sean significativos al establecer una cosmología de Friedmann-Lemaitre-Robertson-Walker. Esto se incorpora en la simulación como una medida evolutiva de la distancia (o factor de escala) en un sistema de coordenadas en movimiento, lo que hace que las partículas se desaceleren en las coordenadas en movimiento (así como debido al desplazamiento al rojo de su energía física). Sin embargo, las contribuciones de la relatividad general y la velocidad finita de la gravedad pueden ignorarse, ya que las escalas de tiempo dinámicas típicas son largas en comparación con el tiempo de cruce de la luz para la simulación, y la curvatura espacio-temporal inducida por las partículas y las velocidades de las partículas son pequeñas. Las condiciones de contorno de estas simulaciones cosmológicas suelen ser periódicas (o toroidales), de modo que un borde del volumen de simulación coincide con el borde opuesto.
Optimizaciones de cálculo Las simulaciones de N cuerpos son simples en principio, porque implican simplemente integrar las 6 N ecuaciones diferenciales ordinarias que definen los movimientos de las partículas en la gravedad newtoniana . En la práctica, el número N de partículas involucradas suele ser muy grande (las simulaciones típicas incluyen muchos millones, la simulación Millennium incluyó diez mil millones) y el número de interacciones entre partículas que necesitan calcularse aumenta del orden de N2, y por lo tanto, La integración de las ecuaciones diferenciales puede resultar costosa o incluso imposible desde el punto de vista computacional. Por lo tanto, comúnmente se utilizan una serie de mejoras.
La integración numérica generalmente se realiza en pequeños intervalos de tiempo utilizando un método como la integración por salto . Sin embargo, toda integración numérica conduce a errores. Los pasos más pequeños dan menos errores pero se ejecutan más lentamente. La integración de Leapfrog es aproximadamente de segundo orden en el paso de tiempo, otros integradores, como los métodos de Runge-Kutta, pueden tener una precisión de cuarto orden o mucho mayor.
Uno de los refinamientos más simples es que cada partícula lleva consigo su propia variable de paso de tiempo, de modo que las partículas con tiempos dinámicos muy diferentes no tienen que evolucionar hacia adelante al ritmo del tiempo más corto.
Existen dos esquemas de aproximación básicos para disminuir el tiempo de cálculo de tales simulaciones. Estos pueden reducir la complejidad computacional a O (N log N) o mejor, con pérdida de precisión.
Métodos de árbol En los métodos de árbol, como la simulación de Barnes-Hut, generalmente se usa un octree para dividir el volumen en celdas cúbicas y solo las interacciones entre partículas de celdas cercanas deben tratarse individualmente; Las partículas en células distantes pueden tratarse colectivamente como una única partícula grande centrada en el centro de masa de la célula distante (o como una expansión multipolar de bajo orden). Esto puede reducir drásticamente el número de interacciones de pares de partículas que deben calcularse. Para evitar que la simulación se vea saturada por el cálculo de las interacciones entre partículas, las celdas deben refinarse a celdas más pequeñas en partes más densas de la simulación que contienen muchas partículas por celda. Para simulaciones donde las partículas no están distribuidas uniformemente, los métodos de descomposición de pares bien separados de Callahan y Kosaraju producen O óptimo ( n registro n ) tiempo por iteración con dimensión fija.
Método de malla de partículas Otra posibilidad es el método de malla de partículas en el que el espacio se discretiza en una malla y, a los efectos de calcular el potencial gravitacional, se supone que las partículas están divididas entre los vértices circundantes de 2x2 de la malla. Encontrar la energía potencial Φ es fácil, porque la ecuación de Poissondonde G es la constante de Newton y es la densidad (número de partículas en los puntos de la malla), es trivial de resolver usando la transformada rápida de Fourier para ir al dominio de la frecuencia donde la ecuación de Poisson tiene la forma simple
dónde es el número de onda comóvil y los sombreros denotan transformadas de Fourier. Desde , el campo gravitacional ahora se puede encontrar multiplicando por y calcular la transformada inversa de Fourier (o calcular la transformada inversa y luego utilizar algún otro método). Dado que este método está limitado por el tamaño de la malla, en la práctica se utiliza una malla más pequeña o alguna otra técnica (como la combinación con un árbol o un algoritmo simple partícula-partícula) para calcular las fuerzas a pequeña escala. A veces se utiliza una malla adaptativa, en la que las celdas de la malla son mucho más pequeñas en las regiones más densas de la simulación.
Optimizaciones de casos especiales Se utilizan varios algoritmos de perturbación gravitacional diferentes para obtener estimaciones bastante precisas de la trayectoria de los objetos en el sistema solar .
La gente suele decidir poner un satélite en una órbita congelada . La trayectoria de un satélite que orbita cerca de la Tierra se puede modelar con precisión a partir de la órbita elíptica de dos cuerpos alrededor del centro de la Tierra y agregando pequeñas correcciones debido al achatamiento de la Tierra, la atracción gravitacional del Sol y la Luna y la resistencia atmosférica., etc. Es posible encontrar una órbita congelada sin calcular la trayectoria real del satélite.
La trayectoria de un pequeño planeta, cometa o nave espacial de largo alcance a menudo se puede modelar con precisión partiendo de la órbita elíptica de dos cuerpos alrededor del Sol y agregando pequeñas correcciones de la atracción gravitacional de los planetas más grandes en sus órbitas conocidas.
Algunas características de las trayectorias a largo plazo de un sistema de partículas se pueden calcular directamente. No es necesario calcular la trayectoria real de ninguna partícula en particular como paso intermedio. Tales características incluyen la estabilidad de Lyapunov, el tiempo de Lyapunov, varias medidas de la teoría ergódica, etc.
Sistemas de dos partículas Aunque hay millones o miles de millones de partículas en las simulaciones típicas, normalmente corresponden a una partícula real con una masa muy grande, normalmente 109 masas solares. Esto puede introducir problemas con interacciones de corto alcance entre las partículas, como la formación de sistemas binarios de dos partículas. Como se supone que las partículas representan una gran cantidad de partículas de materia oscura o grupos de estrellas, estos binarios no son físicos. Para evitar esto, se utiliza una ley de fuerza newtoniana suavizada, que no diverge como el radio del cuadrado inverso en distancias cortas. La mayoría de las simulaciones implementan esto de forma bastante natural ejecutando las simulaciones en celdas de tamaño finito. Es importante implementar el procedimiento de discretización de tal manera que las partículas siempre ejerzan una fuerza evanescente sobre sí mismas.
Suavizamiento