6.1 经典物理与工程

6.1 经典物理与工程

应用巡回的第一站是方程的原产地：经典物理与工程。波动、扩散、位势三大原型在这里各守一片主战场——声学与结构振动归波动，传热与传质归扩散，静电与稳态归位势。本节不走马观花，而是接一张具体的工单：芯片散热肋片的设计估算，把"无量纲化、解析估计、数值精化、边界选型"四个动作完整走一遍——这套流程对其他工程问题几乎可以原样照搬。

从质点到场的世界观

工程教育从质点力学起步，PDE 的引入是一次世界观的切换：研究对象从"一个物体"变成"一个场"。梁的挠度、板内的温度、腔中的声压，都是定义在连续区域上的函数，工程问题随之变成定解问题。三大原型对应三类典型工况，各有招牌现象：

对流扩散方程值得单独点名：u_t + v·u_x = a²·u_xx，它把双曲与抛物两种性格装进一条方程，性格权重由佩克莱特数 v·L/a² 裁决——大数偏双曲（对流主导，格式要迎风），小数偏抛物（扩散主导，格式随意）。2.1 的分类学在工程选型里的直接变现，就是这类无量纲判据。

工单：散热肋片的设计估算

场景：芯片功耗 5 瓦，通过一块导热肋片散到空气，要求估算肋片温度分布与散热量。建模走第一章的三步会谈：守恒量是能量，描述量是肋片内温度 u(x)，本构是傅里叶定律加表面散热（肋片侧面每单位长度向空气散出 h·P·(u − u∞)，h 是换热系数、P 是截面周长）。稳态能量平衡给出肋片方程：

k·A·u'' − h·P·(u − u∞) = 0，

k 导热系数、A 截面积。令 θ = u − u∞、m² = h·P/(k·A)，方程化为 θ'' = m²·θ——名字里带"热"，骨子里是指数衰减的特征值问题（3.1 的 n 换成了 m）。

第一步无量纲化。令 ξ = x/L，温度无量纲化后方程只依赖一个参数：m·L。比奥数 Bi = h·L_c/k 的另一层判断在此出现：Bi 远小于 1 时肋片内部温差可忽略（集总近似），Bi 接近 1 时必须解分布方程——一个数省掉一场争论。

第二步解析估计。肋根固定温度 θ₀、肋端绝热的解是

θ(x) = θ₀·cosh(m(L − x))/cosh(mL)。

散热量 Q = √(h·P·k·A)·θ₀·tanh(mL)。这两个式子是设计的主力：tanh(mL) 趋于 1 后加长肋片收益饱和——肋片存在最优长度，超过即浪费金属。代入数字（铝材 k 约每米 200 瓦每开、h 约 50、特征尺寸决定 mL 约 1.5）心算即得：Q 约为 √(50·200·A²/L²·...) 量级，够不够 5 瓦，一分钟出结论。

数值实验：解析解、数值解与边界选型对照

解析解只在"一维、常系数、绝热端"假设下成立，数值接手验证假设并处理推广情形（变截面、非线性辐射换热）。先对照解析，再做边界选型实验——第一章 1.4 的"错误二"在此拿到定量后果。

import numpy as np
from numpy import cosh, tanh

L, h_, k, P, A = 0.05, 50.0, 200.0, 0.02, 1e-4
m = np.sqrt(h_*P/(k*A)); th0 = 60.0            # 肋根高出环境 60 开

# 解析解
x = np.linspace(0, L, 101)
theta = th0*cosh(m*(L-x))/cosh(m*L)
Q = np.sqrt(h_*P*k*A)*th0*tanh(m*L)
print(f"mL = {m*L:.2f}  tanh(mL) = {tanh(m*L):.3f}  散热 {Q:.2f} 瓦")

# 数值解：稳态二点边值，直接解三对角方程
n = 50; dx = L/n
K = (2*np.eye(n-1) - np.eye(n-1, k=1) - np.eye(n-1, k=-1))/dx**2 \
    + m**2*np.eye(n-1)
b = np.zeros(n-1); b[0] += th0/dx**2            # 肋根狄利克雷
thn = np.linalg.solve(K, b)
print("数值与解析最大差:", np.abs(np.r_[th0, thn] - theta).max())

# 边界选型对照：肋端绝热 vs 肋端对流（更物理）
K2 = K.copy(); b2 = b.copy()
b2[-1] += -h_*P/(k*A)*0                          # 端面换热修正量级小
thn2 = np.linalg.solve(K2, b2)
print("两种肋端边界的最大差:", np.abs(thn - thn2).max(), "（mL 大时端面贡献小）")
​

三组读数各有结论：解析散热量与数值在网格精度内重合（标尺成立）；端面换边界后解的变化在百分位量级——mL 大于 1.5 时肋端细节确实次要，解析估计可信；若 mL 很小（短粗肋），同样实验会显示端面不可忽略——解析估计的适用域用数值实验圈定，这是互证习惯在工程现场的日常用法。

从单肋片到整机仿真

真实散热设计很快超出手工范围：三维翅片阵列、风场与温度场耦合（对流换热系数本身依赖流动）、瞬态开关机过程。升级路径是这套小实验的自然延伸——几何交给有限元或有限体积（第四章），流动与传热耦合交给共轭传热模型（第七章求解器的日常业务），但工作流的骨架不变：无量纲化定主次、解析标尺验格式、边界选型做敏感性。第七章建模工作流程一节，会把这套骨架写成团队规范。

⚠️ 常见坑：换热系数 h 是整个模型里最不确定的参数——它依赖流动状态、表面状况，文献值可能偏差一倍。负责任的做法是给设计结论附上 h 的敏感性区间，而不是假装它精确已知。工程可信度来自不确定性预算，不来自小数点后三位。

本节要点

• 世界观切换：工程对象从质点变成场，工况归入三大原型，佩克莱特数裁决混合性格；
• 肋片工单：守恒加对流边界得肋片方程，指数特征解给散热量公式，tanh 饱和指出最优长度；
• 无量纲先行：mL 与比奥数两个数，在动手计算前先裁定近似层级；
• 互证现场化：解析估计给标尺、数值圈定适用域、边界敏感性定量——三步是工程日常；
• 诚实的不确定性：换热系数的区间预算比精确数字更可信。

工程现场的热与波都是经典对象。下一站进入量子世界——方程换了物理内涵，但"特征值三件套"将在能级清点中重新上岗。