琴弦的振动是怎样的形式?声音在空气中的传播速度如何计算?太阳的寿命大约还剩 50 亿年?4 月 8 日,《张朝阳的物理课》线下第二课开讲,搜狐创始人、董事局主席兼 CEO 张朝阳与来自北京各高校的物理学子及物理学爱好者们一起探讨物理的奥秘。
课堂上,他陆续讲解了经典波动方程的推导与声速的计算,之后又与现场同学交流探讨了太阳内部的核反应、太阳寿命的估算,以及宇宙大爆炸理论的两个实验证据等物理问题。
推导经典波动方程,从琴弦振动开始
本节线下课,张朝阳先介绍了琴弦振动方程的推导。以琴弦初始位置建立 x 轴,弦的振动方向为 u 轴。采用平面运动模型并假设琴弦只能沿一个方向振动,这样就把三维问题简化成了二维问题。弦在某一时刻的位置用 u (x,t) 描述,它表示弦对平衡位置的偏离。
在课程中,假设弦的振动比较微小,u (x,t) 对 x 的偏导数远远小于 1。假设琴弦在初始位置的质量线密度为 μ,考虑在区间 (x,x+Δx) 上琴弦的微元的运动。忽略重力,那么这部分微元只受到弦自身张力的作用。
推导过程中,暂且假设弦的张力 T 处处相等。需要说明的是,因为弦只沿 u 方向振动,不沿 x 方向运动,所以微元受到的 x 方向的合力为零,于是张力在 x 轴上的分量 Tx 处处相等,而 u 方向的张力 Tu 不是处处相等的,从而 T 实际上不是处处相等的。不过,因为接下来的推导需要的是 Tu,当 u (x,t) 对 x 的偏导数远远小于 1 时,无论使用 T 还是 Tx 来估算 Tu,最低阶近似都是一样的,所以可以假设 T 为常数。而当 u (x,t) 对 x 的偏导数不是小量时,则须使用 Tx 为常数来推导,此时依然能得到严格的波动方程。
由于 u (x,t) 对 x 的偏导数远远小于 1,于是弦的切线倾角 θ 非常小,所以 Tu 约等于 T 乘以 θ,而 θ 可以用 u (x,t) 对 x 的偏导数来近似,于是微元受到的 u 方向的合力为:
利用牛顿第二定律,可得:
根据假设,T 为常数,不依赖于 x,对等式右边的偏导数做关于 x 的泰勒展开,并在消去一个 Δx 后让 Δx 趋向于 0,最后得到:
这就是琴弦的波动方程。由于时间关系,对于这个经典的微分方程,现场并未对其求解过程进行详细计算。张朝阳直接写出其通解的大致形式:
随后,他对这个通解作了进一步解释。以其中的 g (x-vt) 为例,假如函数 g (x) 的最大值点在 x0 处,那么在 t 时刻,g (x-vt) 的最大值点在 x-vt=x0 处,也就是 x=x0+vt 处,即最大值点在以速度 v 向右运动。换言之,g (x-vt) 表示向右以速度 v 传播的波。同理,f (x+vt) 表示以速度 v 向左传播的波。经过这一番解释,波动方程中 v 的物理含义就清晰了,它代表波的传播速度。
在这个例子的最后,他还简短地讨论了两端固定的琴弦的波动方程。它本质上是一个边值问题,需要利用傅里叶级数展开的方法来求解,展开的每一项都是满足边界条件的一个模式。
推导空气中声音的波动方程:纵波与密度波
在声波方程的推导中,为了突出基本原理而避免陷入太琐碎的数学细节,本次课程只考虑了一维的情形,并且作了和琴弦类似的假设:静态时处于 x 位置的空气质点在声音过来时的位置偏移 u (x,t) 比较小,并且对空间坐标 x 的偏导数远小于 1。
张朝阳介绍,声波是密度波,是一种纵波,质点的振动方向与传播方向平行,因此不像琴弦振动那样直观。不过,声波方程在数学推导上与琴弦的波动方程是比较类似的。考虑空气在 (x,x+Δx) 的“窄片”,“窄片”截面积是 A,厚度为 Δx。在水平方向上,“窄片”只受到气体的压力,根据牛顿第二运动定律可以得到:
将上式两端的 A 消掉,并且将 P 对 x 做泰勒展开,消掉 Δx 并忽略高阶项可以得到:
“是什么引起了压强的改变呢?答案是密度的改变。”因此,张朝阳转而考虑 P 对密度 ρ 的依赖关系。假设没有声波时空气压强是 P0,密度是 ρ0,那么
看到这里大家可能会疑问,这里 P 对密度的偏导数的准确定义是什么?根据理想气体物态方程,描述理想气体的状态需要两个物理量,而这两个量的选取具有任意性,因此必须明确这里 P 对 ρ 的偏导数是在保持哪个物理量不变的情况下的偏导数。对于这个问题,张朝阳后续又进行详细的解答。
接下来,开始推导 Δρ 和 u 的关系。考虑静态时处于 (x,x+Δx) 的空气“窄片”,横截面积为 A,声波过来时,“窄片”两个截面的偏移分别是 u (x,t) 和 u (x+Δx,t),偏移后的“窄片”质量不变,因此
截面积可以消掉。考虑 u (x+Δx,t) 的泰勒展开,有
代入前一个式子 (消掉 A),并展开可以得到
式中的删除线表示可以消掉或者忽略掉。左边的项和右边第一项互相抵消,由于 u 对 x 的偏导远小于 1,右边最后一项远小于其他项,因此可以忽略,于是:
根据前面推导的结果,有:
代入前面得到的关于 u 的方程,就得到了声音在空气中的波动方程:
声速计算要考虑绝热近似,牛顿曾经在这里栽过跟头
得到波动方程之后,张朝阳开始解释如何理解其中 P 对 ρ 的偏导数。首先,相对于声波导致的快速压缩和膨胀,空气的热传导过程是非常缓慢的。在空气的一个振动周期里,热的传导几乎可以忽略,因此可以把空气振动导致的压缩和膨胀看成是绝热过程,这就是绝热近似。
张朝阳介绍,牛顿曾经在这里栽过跟头。牛顿当时也得到了声波的波动方程,只不过因为对声音导致的空气压缩膨胀过程认识不足,错误地以为是等温过程。回到正题,在绝热近似下,波动方程中 P 对 ρ 的偏导数应该理解为绝热条件下的偏导数。
他进一步解释说,因为这里考虑的声音振幅很小,相对于空气的振动来说,空气微元内恢复平衡的速度非常快,因此可以假设在空间的每一处空气都满足理想气体状态方程。这个近似就是准静态近似。正因为这个近似,这里对声速的推导将不适用于特别稀薄的空气。
在准静态近似下,可以使用理想气体的平衡态性质进行推导。这样就有:
其中 U 是空气内能,因子 5/2 来自于能量均分定理。能量均分定理,每个自由度贡献 NkT / 2 的内能,空气中大部分是氮气和氧气,都是双原子分子,具有 3 个平动自由度和 2 个转动自由度。在室温下,振动自由度没有被激发。所以一共是 5 个自由度,故因子为 5/2。定义 γ 使得 γ-1=2/5。在这里,γ=8/5=1.4。于是
根据热力学第一定律
在绝热过程中有:
其中 C 为常数。由于密度和体积呈反比关系,所以 P 正比于 ρ^γ,设比例常数为 B,于是在绝热过程下
代入理想气体物态方程可得:
其中 m 上加一横表示空气的平均分子质量。所以,空气中声波的波动方程为:
在室温 20℃,即约 293 K,相应的声速为
其中 m_m 表示空气以 g 为单位的摩尔质量。这个数值结果与实验测得的声速非常接近。另外,还可以知道声速只与温度有关,而与声音的频率无关。
“这是一个令人震撼的结果。”张朝阳说,波动方程的出发点仅仅是牛顿第二运动定律,再考虑理想气体物态方程,就可以得到与实验值非常接近的计算值,不得不说计算的力量是强大的。
解释太阳能量源于核聚变,估算太阳寿命约为百亿年
在课程后半段的互动环节,有同学提问关于太阳寿命的问题,于是张朝阳对太阳的能量来源以及如何估算太阳寿命进行了解答。
他介绍,太阳内部有一个核心,太阳大部分的质量都集中在核心。核心内部温度非常高,约为 100 万开尔文。太阳的能量来源是核心处的核聚变,比如氢聚变成氦。由于氢核带正电,氢核之间存在库仑排斥力,因此存在势垒。但是强相互作用会使得氢核在靠得特别近时会释放能量形成新的原子核,也就是说在库仑势垒之内有一个深势阱。要想发生核聚变,氢核必须突破库仑势垒。静电势大小为
以原子核半径的量级 10^(-15) m 来估算库仑势垒高度,得到结果约为 10^(-13) 焦耳。再根据太阳的内核温度,利用 kT 可以估算内核粒子动能的量级为 10^(-16) 焦耳,这两个能量差了 3 个量级。
如果仅仅如此粗略地看,核聚变似乎不会发生。那为什么核聚变还是真实地发生了呢?可以从两个角度来理解:第一个角度是从经典的麦克斯韦速度分布出发,虽然内核粒子的平均动能突破不了库仑势垒,但是麦克斯韦速度分布允许更高动能粒子的存在,从而提供了穿越库仑势垒的可能;另一个角度是从量子隧穿出发,虽然根据经典力学,粒子无法穿过势垒,但是量子力学的隧穿效应保证了粒子实际上存在一定概率穿过库仑势垒。这就是对于为什么能发生核聚变的定性分析。
知道了太阳的能量来源,就可以估算太阳的寿命了。在太阳内部主要发生的是氢核聚变成氦核,只考虑最初和最末状态的话,就是四个氢核变成一个氦核,释放的能量为 ΔE=Δmc^2。于是,通过“质量亏损”可以计算得到:聚变产生一个氦核释放的能量为 0.4×10^(-11) 焦耳。
假设太阳最初全部由氢构成,那么可以根据太阳质量估算出有多少氢核。但是,不是所有氢都能聚变成氦,张朝阳在课上假设能变成氦的氢占比 10%,那么可以估算出太阳可以释放的能量为 0.12×10^(45) 焦耳。最后,再假设在太阳整个生命周期里表面温度近似不变,于是太阳将以恒定的温度以黑体辐射的形式向外辐射能量,这样就可以得到太阳寿命为
代入数值,并转化为以年为单位,可以得到太阳寿命约为 100 亿年。目前太阳已经燃烧了 50 亿年,剩余寿命大约还剩 50 亿年。在太阳的末期,太阳将会膨胀成为红巨星,会把周围的很多行星甚至可能包括地球都吞噬掉。
宇宙大爆炸理论的两个证据:哈勃红移和微波背景辐射
对于现场同学们感兴趣的宇宙大爆炸问题,张朝阳介绍了支持宇宙大爆炸理论的两个观测证据,分别是哈勃红移和微波背景辐射。
根据光的多普勒效应,以速度 u 远离地球而去的恒星,其光谱会发生偏移。(u<0 表示恒星向我们靠近),从这个恒星上发出频率为 ν0 的光,而被地球接收到的频率为
假如 u 远小于光速 c,那么进行泰勒展开并忽略高阶项,可以得到频率改变量 Δν=ν-ν0 满足
相应的波长变化量是
于是利用这个公式,根据恒星的光谱移动,就可以计算出恒星远离地球的速度。哈勃当年观测了很多恒星的频移,发现绝大部分都是红移,而且恒星远离地球的速度与恒星相对地球的距离成正比。尽管哈勃发现宇宙膨胀在先、大爆炸理论提出在后,然而从逻辑上来看,哈勃红移是大爆炸理论的一个结果,因此可以作为大爆炸理论的一个实验证据。
根据大爆炸的理论模型,宇宙从最初的极高温,经过膨胀而不断冷却,最后变成了现在这个样子。期间,物质的形成也是由宇宙温度降低所导致的。当宇宙温度降低时,夸克这些粒子会禁闭成重子。大部分重子是不稳定的,最后会变成质子和中子。中子寿命比较长,会衰变成质子和电子,外加一个反中微子。
不过这个过程是可逆的,最终会形成平衡,导致有一个不为零的中子丰度。当温度继续降低,直到光子能量不足以打破原子核时,这些中子、质子将会通过核聚变成为稳定的原子核,这个过程被称为“原初核合成”。但是因为温度依然很高,物质以离子形态存在,光子和这些带电粒子存在相互作用,从而导致宇宙是“不透明”的。
到了大约 38 万年的时候,温度已经进一步下降,电子与原子核开始结合成原子,并会释放出光子,这时候宇宙中的粒子大部分都是电中性的了,光子可以在宇宙中几乎畅通无阻。这个时期就是光退耦时期,光子不再与物质有大量的热交换,而主要以黑体辐射的形式存在,随着宇宙膨胀其温度不断下降。直到今天,这部分光子,也就是目前所说的微波背景辐射,对应的等效黑体温度大约为 2.7 开尔文。
介绍完黑体辐射的形成,张朝阳还根据黑体辐射的性质计算了微波背景辐射的峰值频率,约为 168 GHz。与哈勃红移一样,微波背景辐射的发现也比大爆炸的理论提出要早,不过同样可以将微波背景辐射看作是大爆炸理论的实验证据。