Cómo modelar la formación de planetas en nuestro sistema solar

Modelar la formación de planetas es un problema extremadamente difícil y se realiza por etapas. Para las primeras etapas necesitas hidrocódigos. Para las últimas etapas, necesitas códigos de N cuerpos. Las cosas entre las personas todavía están tratando de averiguar qué hacer. Estoy más familiarizado con los códigos N-body utilizados para la acumulación de etapa tardía. Los principales actores son Swift y Symba de Hal Levison y Mercury de John Chambers. Esos están escritos en Fortran. La mayoría de los códigos hyrdo que conozco también están escritos en Fortran o C. Python es un lenguaje horrible para el trabajo de simulación numérica, especialmente para este tipo de problemas que realmente requieren muchos cálculos rápidos de números.

En cuanto a la gravedad dependiente del tiempo, puede modificar cualquiera de los anteriores para incluirlo siempre que se sienta cómodo editando el código Fortran.

No creo que la haya. Demasiadas partículas Hace años escribí una simulación de pozo de gravedad y cuando se aplicó la gravedad a todos los objetos del sistema, comenzó a arrastrarse a unos 35 objetos más o menos. No lo he probado en hardware moderno, pero es difícil creer que una computadora personal pueda manejar la carga necesaria para que realmente “funcione” como imaginas. Hay cosas como el algoritmo de bandada de pájaros que te permitirán manejar más objetos, pero en realidad ya no estás hablando de la gravedad. Y, después de todo esto, solo estaba tratando de hacer que los “planetas” giraran alrededor del “sol”, en lugar de modelar la formación de un sistema solar. ¡Buena suerte sin embargo!