热方程 · 量与起伏

留下的是平均温度

把两端保温,温度会慢慢变得均匀,却不一定降到零。改动平均温度 C,看看总量、能量和起伏中的哪一个跟着改变。

uₜ = κuₓₓ  E = ½∫u² dx  M = ∫u dx

温度剖面

实线:当前温差虚线:初值点线:当前平均值
总量 M
当前平均值 M/L
能量 E
起伏能量 ½∫(u−平均值)²
耗散 D = κ∫uₓ²
能量变化率 E′
总量变化率 M′
可衰减部分的上界

能量怎样下降

绝热时勾选,比较全能量与起伏;取消勾选,纵轴会放大到起伏的尺度。零端值时,可衰减部分就是全能量。

实线:可衰减部分虚线:第一模态速率给出的上界紫色:全能量 E

哪一模态先消失?

实心:当前模态能量轮廓:初始模态能量

每个非恒定模态的能量为 L aₙ²/4。横轴 n 是模态编号;常数背景另计入总能量,不画在这三个柱子里。

为什么绝热时要减去平均值?

绝热端没有热流穿过,M′=κ[uₓ(L)−uₓ(0)]=0,因此平均温度 C 保持不变。此时 E=LC²/2+起伏能量。指数估计控制的是后者;除非 C=0,否则总能量不会趋于零。

零端值时 E′=−D≤0,但 M′ 由两端斜率决定。有正有负的温差叠加时,不该把“能量下降”换成“带符号的总量单调下降”。零端值图中的起伏能量采用当前平均值计算,仅供比较;它不等于用于指数估计的 E。

kₙ=nπ/L,aₙ=Aₙe^(−κkₙ²t)
D=κL/2 Σ kₙ²aₙ²
零端值:E ≤ E(0)e^(−2κπ²t/L²)
绝热:起伏能量 ≤ 起伏初值 e^(−2κπ²t/L²)