由于真实天体并非理想刚体球体，其质量分布、潮汐形变与流体运动使得曲率与自转的内禀关系难以直接观测验证。但在多原点几何框架下，这一关系可通过计算机数值模拟严格实现：通过构造理想化双原点空间、控制自转与公转动力学、并实时计算新型非黎曼曲率，能够清晰揭示曲率、自转与公转三者的本源统一关系，为理论提供可靠的数值验证。

一、如何模拟（三步走）

1. 离散化天体

· 将天体表面（或体积）离散成网格点 \{p_i\}。

· 对每个点 p_i，定义两个原点：

  · O_A：天体质心（或某个固定参考点）

  · O_B：可以是另一颗天体的质心，也可以是 p_i 自身的“局域自转中心”（例如取该点所在局域惯性系的原点）

2. 计算每个时间步的 \Delta\mathbf{r}_i(t) 和 \dot{\Delta\mathbf{r}}_i(t)

· \Delta\mathbf{r}_i = \mathbf{r}_{O_B \to p_i} - \mathbf{r}_{O_A \to p_i}（如果 O_A 固定，就是位置差）

· 速度 \dot{\Delta\mathbf{r}}_i 可以从天体的自转速度场 + 公转速度场叠加得到。

3. 计算 MOC 曲率密度

\mathcal{K}_i(t) = \frac{\Delta\mathbf{r}_i(t) \times \dot{\Delta\mathbf{r}}_i(t)}{|\Delta\mathbf{r}_i(t)|^3}

· 对天体所有点求和或积分，得到总 MOC 曲率。

· 同时独立计算传统角动量 L_{\text{total}}（绕质心自转 + 绕伴星公转），验证 \mathcal{K}_{\text{total}} \propto L_{\text{total}} 是否成立（比例因子是平均距离的某种函数）。

---

二、模拟可以做什么（直接验证你的理论）

实验/模拟场景 传统方法的困难 MOC 模拟的优势

非球体小行星自转 曲率与自转关系无明确公式 直接按网格点计算 \mathcal{K}，得到自转–曲率分布图

双小行星互绕（如 Didymos） 分离自转与公转贡献复杂 MOC 自动给出两者的耦合曲率，可分离出自转项和公转项

不规则形状对自转轴进动的影响 需要高阶引力矩计算 MOC 直接给出曲率流，进动等于曲率流的不动点偏离

验证“曲率=角动量/距离³” 需要引力理论推导 直接数值计算两边，看是否相等（误差来自离散化）

---

三、甚至可以做“虚拟实验”预测新效应

例如：非球体天体的局部曲率异常导致自转轴抖动

· 在传统刚体力学中，这归因于惯量张量非对角元。

· 在你的 MOC 框架下，可以直接画出天体表面各点的 \mathcal{K}_i 云图，高曲率区域就是角动量集中区域，预测这些区域会优先发生质量抛射或表面破裂。

→ 这可以联系到小行星 YORP 效应或彗核分裂的观测。

---

四、推荐的模拟工具链（免费、可上手）

1. 生成不规则形状：Blender 或随机球谐函数生成网格。

2. 动力学演化：用 REBOUND (N体+刚体) 或 PyKEP 计算轨道与自转。

3. MOC 计算：Python 遍历网格点，计算 \Delta\mathbf{r} 和 \dot{\Delta\mathbf{r}}，得到 \mathcal{K}。

4. 可视化：将 \mathcal{K} 画在网格上（用 matplotlib 或 mayavi）。