返回首页
最新
我正在开发一个名为 Astroworld 的模块化 N 体模拟器,使用 Python 编写。最初它是一个太阳系可视化工具,但我最近将其重构为一个通用引擎,能够将物理法则与行星数据解耦。技术亮点包括:
- **辛积分法**:使用速度 Verlet 积分器来保持长期的能量守恒(在稳定系统中,$\Delta E/E \approx 10^{-8}$)。
- **无关架构**:可以通过轨道元素(开普勒元素)或状态矢量输入任何系统。我已经用它验证了超紧凑系统(如 TRAPPIST-1)和长期扰动(如行星九假说)的稳定性。
- **验证**:包括90多个物理测试,其中包括使用施瓦兹希尔德度规修正的水星相对论进动。
**行星九实验**:我进行了一个1万年的模拟,以跟踪像 Sedna 这样的远日天体(TNO)在近日点参数($\omega$)中的差异信号。结果(约 $0.002^{\circ}$)为引擎的精度提供了很好的合理性检查,因为这一效应是世俗的,需要数百万年才能完全显现。
**技术栈**:使用 NumPy 进行向量化,Matplotlib 进行二维分析,以及 Plotly 进行交互式三维轨迹展示。我目前正在开发一个实时三维渲染层。希望能收到关于高偏心率轨道的积分器稳定性反馈,或者关于实现更复杂引力势的建议。
独立开发者,白天是IT经理。创建这个是为了自己解决问题——非常希望能得到反馈。