Python 2D Jean không ổn định trong mô phỏng thiên hà xoắn ốc

Tôi có ~ 1000 ảnh chụp mô phỏng thiên hà xoắn ốc (từ nhiễu loạn phân tử đám mây t = 0 đến t> 1Gyr. Điều tôi cần làm là xác định xem bất kỳ vị trí nào tại bất kỳ thời điểm chụp nhanh nào đều đáp ứng tiêu chí cho tính không ổn định của Jeans .

Mô phỏng là 2D và tập lệnh chụp nhanh có thể đọc vị trí, vận tốc và mật độ bề mặt, nhưng không phải là nhiệt độ. Mô phỏng cũng không có khí, thay vào đó mô phỏng được mô hình hóa bằng cách cho mỗi khối sao riêng lẻ tương ứng với khí. Tôi đã xem trực tuyến và một số trang web cho biết khối lượng Jeans trong 2ngày bằng c / G *. Nhưng tôi không biết làm thế nào để xác định tốc độ âm thanh c . Ngoài ra tôi không thực sự quan tâm đến khối lượng Jeans, hơn nữa liệu bất kỳ vị trí nào cũng có thể sụp đổ để tạo ra các ngôi sao, tức là thỏa mãn sự bất ổn của Jeans.$_{s}$$^4$$^2$$_{s}$

Đầu tiên là cảnh báo, sau đó là hack:

1) Nếu bạn không chạy mô phỏng này, nếu bạn không phải là nhà thiên văn học lý thuyết được đào tạo ở trình độ đại học trở lên và nếu mô phỏng này không có loại dữ liệu bạn cần trả lời câu hỏi bạn không nên tin tưởng kết quả ở tất cả .

Phải mất nhiều năm nghiên cứu để hiểu được các hiệu ứng số kỳ lạ (những điều xảy ra trong mô phỏng vì cách nó được thiết lập sẽ không bao giờ xảy ra trong thực tế), vật lý nền và vô số giả định và xấp xỉ đi vào loại công cụ này.

Chúng ta thường nói 'rác vào, rác ra'. Các câu trả lời bạn nhận được từ điều này có lẽ sẽ không được chiếu sáng lắm vì chúng nói nhiều về mô phỏng hơn là không gian thực.

2) Điều đó nói rằng, nếu bạn chỉ làm điều này cho vui hoặc để học các kỹ năng, TUYỆT VỜI! Tiếp tục chơi! Và một giả định rất cơ bản mà bạn có thể đưa ra để rút ra một số con số (điều đó sẽ khá sai) là tất cả khí trong không gian là ~ 10K.