Instabilidade Jean Python 2D na simulação de galáxia espiral

Eu tenho ~ 1000 instantâneos de uma simulação de galáxia em espiral (das perturbações das nuvens moleculares t = 0 a t> 1Gyr. O que preciso fazer é determinar se qualquer posição em qualquer um dos momentos do instantâneo atende aos critérios de instabilidade do Jeans .

A simulação é 2D e o script de instantâneo pode ler a posição, velocidade e densidade da superfície ∑, mas não a temperatura. A simulação também não tem gás; em vez disso, a simulação foi modelada, dando a cada estrela individual massa extra que corresponde ao gás. Procurei on-line e vários sites afirmam a massa Jeans em 2d igual a c s 4 / G 2 * ∑. Mas não tenho idéia de como determinar a velocidade do som c s . Também não estou realmente interessado na massa de Jeans, mais ainda se alguma posição pode entrar em colapso para formar estrelas, ou seja, satisfazer a instabilidade de Jeans.$_{s}$$^4$$^2$$_{s}$

Primeiro um aviso, depois um hack:

1) Se você não executou esta simulação, se você não é um astrônomo teórico treinado de nível de graduação ou superior e se essa simulação não possui o tipo de dados, você precisa responder à pergunta que está fazendo , não deve confiar os resultados .

Demora anos de estudo para entender os efeitos numéricos estranhos (coisas que acontecem na simulação devido à maneira como ela é configurada e que nunca aconteceria na realidade), física de fundo e muitas suposições e aproximações que entram nesse tipo de coisa.

Costumamos dizer 'lixo dentro, lixo fora'. As respostas que você obtém disso provavelmente não serão muito esclarecedoras, porque dizem mais sobre a simulação do que o espaço real.

2) Dito isto, se você está apenas fazendo isso por diversão ou para aprender as habilidades, ÓTIMO! Continue jogando! E uma suposição muito básica que você pode fazer para obter alguns números (o que será bastante errado) é que todo gás no espaço é ~ 10K.