2.1.1 数学表述：矢量形式、势函数（Φ = −Gm/r）与泊松方程（∇²Φ = 4πGρ）

文档摘要

2.1.1 数学表述：矢量形式、势函数（Φ = −Gm/r）与泊松方程（∇²Φ = 4πGρ）在数值天体力学与计算宇宙学的工程实践中，牛顿万有引力理论远非教科书里几行优雅公式那般轻盈——它是一套必须被“编译”进内存、被离散化为矩阵、被迭代求解、被边界条件反复校验的可执行物理模型。当你在GPU集群上模拟一个包含十亿粒子的暗物质晕演化时，真正驱动每一帧时间步进的，并非开普勒的椭圆，而是$\nabla^2 \Phi = 4\pi G \rho$在三维网格上的有限差分离散；当你调试一个航天器轨道预报模块出现0.3 km/s的速度残差时，问题根源往往不是初始状态误差，而是势函数$\Phi(\mathbf{r}) = -G m / |\mathbf{r}|$在近场区未做正则化处理所引入的奇异性污染。