Python 2D Jean Instability in Spiral Galaxy Simulation

Ho ~ 1000 istantanee di una simulazione di galassia a spirale (da perturbazioni della nuvola molecolare da t = 0 a t> 1 giro. Quello che devo fare è determinare se qualsiasi posizione in uno qualsiasi dei tempi dell'istantanea soddisfa i criteri di instabilità di Jeans .

La simulazione è 2D e lo script dell'istantanea può leggere posizione, velocità e densità della superficie ∑, ma non la temperatura. Anche la simulazione non ha gas, invece la simulazione è stata modellata dando ad ogni singola stella una massa extra che corrisponde al gas. Ho guardato online e diversi siti dichiarano la massa di Jeans in 2d pari a c / G * ∑. Ma non ho idea di come determinare la velocità del suono c . Inoltre, non mi interessa molto la massa dei Jeans, soprattutto se una posizione può collassare per creare stelle, ovvero soddisfare l'instabilità dei Jeans.$_{s}$$^4$$^2$$_{s}$

Prima un avvertimento, poi un trucco:

1) Se non hai eseguito questa simulazione, se non sei un astronomo teorico addestrato di livello universitario o superiore e se questa simulazione non ha il tipo di dati di cui hai bisogno per rispondere alla domanda che ti stai chiedendo non dovresti fidarti i risultati a tutti .

Ci vogliono anni di studio per comprendere gli strani effetti numerici (cose che accadono nella simulazione a causa del modo in cui è impostato che non accadrà mai nella realtà), fisica di fondo e un sacco di ipotesi e approssimazioni che vanno in questo tipo di cose.

Spesso diciamo "immondizia, immondizia". Le risposte che otterrai da questo probabilmente non saranno molto illuminanti perché dicono di più sulla simulazione che sullo spazio reale.

2) Detto questo, se lo stai facendo solo per divertimento o per imparare le abilità, GRANDE! Continuare a giocare! E un presupposto molto semplice che puoi fare per ottenere alcuni numeri (che sarà abbastanza sbagliato) è che tutto il gas nello spazio è ~ 10K.