lo que vemos en tu simulación son los primeros momentos del inicio del planeta, y además vemos lo que pasa en el exterior, pero no lo que pasa en el interior.
Pero DCM, se supone que la simulación está recreando el interior de una esfera de radio \( R \) con partículas dentro! Mira lo que pasa con una esfera hueca de radio \( 200 \) unidades:

video de la simulación.

Como podéis ver en el interior la materia de las capas más interiores (color verde) se desplazan hacia la superficie y la materia de las capas superiores (color azul) se desplazan hacia el centro. Llegando a crearse una esfera hueca de un radio menor al inicial. La materia está siendo comprimida.
Por curiosidad, ¿cómo has conseguido que las partículas formen una esfera hueca? ¿Has metido las fórmulas físicas (¿cuáles y dónde?) y ha ocurrido eso de forma espontánea, es decir, sin ninguna otra acción por tu parte, o cómo?
Hola DaniM, por favor revisa el mensaje siguiente, para entender lo que pasa en la esfera maciza de densidad de material homogénea.
https://foro.rinconmatematico.com/index.php?topic=123119.msg513712#msg513712 Pongo sólo el esquema

En resumen el hecho que las esferas huecas de espesor finito tenga un campo gravitatorio no nulo , en su interior provoca que la esfera maciza sea inestable y comience el proceso de generación de la esfera hueca final que te describo.
Quiero compartir que esta simulación no la he realizado yo, sino que me ha ayudado un joven físico (Álvaro), al que estoy muy agradecido.
La simulación que he compartido es muy simple pues sólo se ha tenido en cuenta la dirección y sentido de la gravedad que está actuando sobre cada partícula que formaría la esfera maciza de densidad homogénea, pero nada sobre su magnitud y cuando las partículas llegan a un radio determinado, donde se produciría el choque entre los 2 frentes de partículas, hacemos desaparecer las partículas.
Esto está realizado para mostrar a nivel cualitativo y no cuantitativo, lo que ocurría en esta esfera maciza y cómo se transformaría en la esfera hueca.
Tienes que tener claro una cosa; queremos saber qué ocurre dentro de la esfera maciza con una densidad homogénea de materia, pero ésta no está compuesta de un número de partículas finitas, en mi simulación se ha sustituido esta densidad de materia por un número finito de partículas , formando capas, para ver el movimiento de las capas. De forma que viendo cómo se comportan las partículas testigo de cada capa vemos que ocurre con el material de esa capa.
Las partículas son como testigos para ver qué ocurre con cada capa de material.
DaniM, tu simulación es buenísima, con los conocimientos que hasta hoy conocemos, y por ello lo que se obtiene es que después de un tiempo en una región del espacio con masas finitas todas se forman en una única de un diámetro determinado y una masa suma de todas las partículas. Pero una vez que se ha formado una esfera maciza de una radio R, y una masa M, con una densidad \( \rho \) homogénea, no se ha tenido en cuenta la gravedad que existe dentro de esta esfera maciza de densidad homogénea, actuando sobre las capas de material de ella.
Para hacer una simulación real de esto se debe utilizar las ecuaciones de la gravedad en cada punto de la esfera maciza de densidad homogénea.
Las ecuaciones de la gravedad dentro de una esfera maciza de densidad homogénea que he obtenido son las siguientes:
Ecuaciones

La variable r es la distancia al punto de cálculo de la gravedad, desde la superficie.
M es la masa total de la esfera maciza.
Indicar que la gravedad positiva va en sentido hacia el centro de la esfera y negativa hacia el exterior de la esfera.
Importante: Como no es posible crear en el ordenador una esfera maciza de densidad homogénea, sería conveniente simular una región del espacio esférica de radio R (sería el radio de la esfera maciza homogénea), en la que la gravedad que exista sea la de las ecuaciones anteriores y en ésta ponemos partículas de masas todas iguales testigos, por capas, poniendo de 5 a 10 capas, distribuidas a lo largo del radio R, siendo la masa de la esfera la suma de las masas de todas las partículas. La distribución de las partículas es importante que se haga de la siguiente forma:
esquema de distribución de piezas.

En el caso real, nos imaginamos que unas capas de material chocarán con otras sin problemas, pero aquí como tenemos un nº finito de partículas, si no las alineamos en la misma dirección, los choques entre capas no se podrán simular.
Importante: Las ecuaciones de arriba son para una esfera maciza de densidad homogénea, pero al actuar la gravedad las capas inferiores se irán hacia la superficie y entonces esta esfera maciza se transformará en una esfera hueca de un gran espesor (casi el radio R) al principio y este espesor se irá reduciendo con el tiempo.
En definitiva las ecuaciones en este caso ya no serán las de las esferas macizas sino las de las esferas huecas de espesor que exista en cada instante.
Esto se complica enormemente pero podemos suponer que las ecuaciones de la esfera maciza son vigentes durante el desplazamiento de las capas hacia la superficie, excepto cuando ya todas las capas hayan chocado unas con otras (sería la situación en la que la partícula más interior choca con la partícula más exterior). Esta situación se producirá a una distancia, R1, del centro de la esfera.
En esta circunstancias tendremos una esfera hueca de un espesor muy pequeño respecto a R1.
En este momento las ecuaciones se transforman a las ecuaciones de una esfera hueca de muy pequeño espesor en comparación con el Radio y podemos simplificar :
r' es la distancia de las partículas al centro de la esfera:
para \( r'>R1 \) \( g= \frac{GM}{2R1^2} \) (positiva hacia el centro de la esfera)
para \( r'< R1 \) \( g= - \frac{GM}{2R1^2} \) (negativa hacia la superficie de la esfera)
DaniM muchas gracias por todo, por tu enorme aportación.
Si deseas seguir viendo qué pasa dentro de las esferas macizas tú tienes la decisión de ir adelante o parar. Lo que hagas está bien hecho.
Cualquier aclaración o duda que tengáis ya sabéis donde estoy.
DaniM muchas gracias por todo.
Un saludo.