5.3 燃料账本：静能、动能与总能量

5.3 燃料账本：静能、动能与总能量

本站摘要：把能量账分成三栏——静能 mc 平方、动能 (γ−1)mc 平方、总能量γmc 平方。本节给一吨重的掠光号逐速记账，算出"从 0.99c 提到 0.999c 比从零加速到 0.99c 还贵两倍"的账目，并验证低速极限退回经典动能、恒力下加速度跌到牛顿值的三百五十分之一。

5.2 开出质能账单，本站逐笔登记。账本同时回答第二章的伏笔——为什么有质量的物体永远到不了光速。

燃料账本

给一吨（1000 千克）的列车记账，能量单位用焦耳，另附参照物让数字有体感。三栏关系一目了然：总能量恒等于静能加动能，静能是"存在本身"的存款，动能才是燃料烧出来的部分。

# 演算 1：一吨列车的能量三栏账
import math
c = 299_792_458
m = 1000.0                      # 1 吨
E0 = m * c**2                   # 静能
print(f"静能 E0 = {E0:.4e} 焦耳")     # 输出: 8.9876e+19 焦耳

def gamma(beta):
    return 1 / math.sqrt(1 - beta**2)

GLOBAL_ELEC = 1.0e20            # 全球年发电量量级（约 28000 太瓦时）
for beta in [0.6, 0.8, 0.9, 0.99, 0.999]:
    g = gamma(beta)
    Ek = (g - 1) * E0
    print(f"v={beta}c: γ={g:8.3f}  动能={Ek:.3e} 焦耳  相当于全球年发电 {Ek/GLOBAL_ELEC:.2f} 倍")
# 输出:
# v=0.6c:  γ=   1.250  动能=2.247e+19 焦耳  相当于全球年发电 0.22 倍
# v=0.8c:  γ=   1.667  动能=5.992e+19 焦耳  相当于全球年发电 0.60 倍
# v=0.9c:  γ=   2.294  动能=1.163e+20 焦耳  相当于全球年发电 1.16 倍
# v=0.99c: γ=   7.089  动能=5.472e+20 焦耳  相当于全球年发电 5.47 倍
# v=0.999c:γ=  22.366  动能=1.920e+21 焦耳  相当于全球年发电 19.20 倍
​

一吨物质提到 0.9c 就要烧掉相当于全球一年的发电量——能量账单第一次让人对"光速很近但够不着"有了肌肉记忆。

提速的边际成本

账本最刺眼的一栏是边际成本：速度每前进一小步，账单如何暴涨。

# 演算 2：每一步提速的燃料开销（单位 mc²）
import math
def gamma(beta):
    return 1 / math.sqrt(1 - beta**2)

steps = [(0.0, 0.6), (0.6, 0.8), (0.8, 0.9), (0.9, 0.99), (0.99, 0.999)]
cum = 0.0
for a, b in steps:
    cost = gamma(b) - gamma(a)           # 动能增量（以 mc² 计）
    cum += cost
    print(f"{a}c → {b}c: 新增 {cost:7.3f} mc², 累计 {cum:7.3f} mc²")
# 输出:
# 0c → 0.6c: 新增   0.250 mc², 累计   0.250 mc²
# 0.6c → 0.8c: 新增   0.417 mc², 累计   0.667 mc²
# 0.8c → 0.9c: 新增   0.628 mc², 累计   1.294 mc²
# 0.9c → 0.99c: 新增   4.795 mc², 累计   6.089 mc²
# 0.99c → 0.999c: 新增  15.277 mc², 累计  21.366 mc²
​

读表：最后一步（0.99c 到 0.999c）的开销，是此前全部旅程累计的三倍半还多。而且趋势没有尽头——γ在 c 处发散，把有质量的物体推到光速需要无穷能量。第二章埋的伏笔至此结案：光速上限不是交警开出的罚单，是燃料账本的数学必然。

低速对账与力的旧公式

低速区新账本必须与牛顿账本吻合。以巡航客机（每秒 250 米）为例核对；再看恒力下场（力等于动量对时间的变化率，动量换成γ乘 m 乘 v 之后力与加速度的关系被γ三次方拉开）：

# 演算 3：低速极限与恒力加速
import math
c = 299_792_458
m, v = 1000.0, 250.0            # 1 吨，250 m/s 巡航

# ① 经典动能 vs 相对论动能
Ek_newton = 0.5 * m * v**2
beta = v / c
g = 1 / math.sqrt(1 - beta**2)
# 数值稳定写法：γ-1 = β²γ²/(γ+1)，避免两个近似 1 的数相减
Ek_rel = m * c**2 * beta**2 * g**2 / (g + 1)
print(f"经典动能 {Ek_newton:.6e} J, 相对论动能 {Ek_rel:.6e} J")
print(f"相对差 {abs(Ek_rel-Ek_newton)/Ek_newton:.2e}")
# 输出: 经典动能 3.125000e+07 J, 相对论动能 3.125000e+07 J
# 相对差 5.22e-13 —— 万亿分之一量级，日常世界两种账本无法区分

# ② 恒力下的加速度衰减：a = F/(γ³m)（纵向）
def gamma(b):
    return 1 / math.sqrt(1 - b*b)
for beta in [0.1, 0.6, 0.9, 0.99]:
    print(f"v={beta}c: 加速度只剩牛顿值的 1/{gamma(beta)**3:.0f}")
# 输出:
# v=0.1c: 加速度只剩牛顿值的 1/1
# v=0.6c: 加速度只剩牛顿值的 1/2
# v=0.9c: 加速度只剩牛顿值的 1/12
# v=0.99c: 加速度只剩牛顿值的 1/356
​

演算 3 的①里藏着一个数值计算的真坑：若直接写"γ减 1"，两个极其接近 1 的数相减会把 16 位有效数字毁掉一大半，算出的相对差会假报成万分之几——这是浮点灾难性消元的典型现场。改用恒等式γ减 1 等于β平方γ平方除以 (γ 加 1) 后，真实差距显现为 5 乘 10 的负 13 次方量级：低速世界里，牛顿账本精确得令人尊敬。

②则解释了加速器工程师的日常：同样的推力，0.99c 时提速效果只有牛顿预期的三百五十六分之一——不是机器坏了，是账本的γ三次方在起作用（动量对速度的斜率正是 mγ³，5.2 推导时已见过它）。

⚠️ 记账误区：把"静能"当成可以随意支取的油箱。静能确实巨大，但普通过程只允许动用质量差额——化学反应改变质量小数点后十位，核反应到第三四位；要全额支取，只有正反物质湮灭（迄今无法大规模储存反物质）。账本上的存款与柜台的流通性是两回事。

💡 换个单位看世界：粒子物理里能量直接用电子伏特、质量也用电子伏特每 c 平方——"质子质量 938 兆电子伏"说的是静能。单位合并到这种程度，正是质能同一性的行业日常。6.1 节的对撞机账目将全程用这套单位。

本节要点回顾

• 三栏账：静能 mc²、动能 (γ−1)mc²、总能量γmc²，总恒等于静能加动能
• 体感数字：一吨列车到 0.9c 需约 1.16 乘 10 的 20 次方焦耳，约全球一年发电量
• 边际暴涨：0.99c 到 0.999c 一步花费 15.3 份静能，超过此前全部累计
• 上限结案：γ在 c 处发散，加速到光速需无穷能量——燃料账本层面的证明
• 低速对账：客机速度下新旧动能相差约 5 乘 10 的负 13 次方，牛顿账本是极好近似
• 恒力衰减：纵向加速度按γ三次方缩水，0.99c 时只剩三百五十六分之一

三栏账齐了。最后一笔：能量与动量合成直角三角形——货运部的镇店之宝。