一维热传导方程(第二类边界:诺伊曼条件)
热传导方程中的第二类边界条件(诺伊曼)是指边界热流已知的，它与第一类边界（固定温度）不同的是，第一类对于边界控制量是温度本身，第二类边界控制的是温度的变化率（梯度）,第二类边界条件直接控制了能量交换的速率，而不是温度本身。因此，它能更直接地描述“加热”或“冷却”的强度。
这是一个非常实际的问题。当边界条件给定为热流（第二类边界条件，诺伊曼条件）而非温度（第一类，狄利克雷条件）时，求解过程在数学上有所不同，但核心方法（分离变量法）依然有效。 下面我以一个一维有界区域 $0 < x < L$，左端绝热、右端有恒定热流输入为例，为你完整演示计算过程。
一、问题设定
考虑一根长度为$L$的均匀细杆，侧面绝热。其热传导方程为：
$$ u_t = \alpha u_{xx}， \quad 0 < x < L, t > 0 $$
边界条件（第二类）：左端绝热：$u_x(0,t)=0$（热流为零）；右端有恒定热流：$-ku_x(L,t)=q$（热流密度 $q$, 其中$q$是流入物体的热流密度（单位面积功率）。若 $q>0$,热量流入；若 $q < 0$, 热量流出。 $k$ 为导热系数）
初始条件：$$ u(x,0)=f(x)$$
二、计算过程
第一步：无量纲化:
为了简化，通常将热流边界条件转化为更标准的形式。令：
$$ u_x(L,t) = -\frac{q}{k} = \gamma$$
其中$\gamma$ 是常数（若 $q>0$，则 $\gamma < 0$，表示温度在右端下降），求解的问题为：
$$ \begin{cases} u_t = \alpha u_{xx}, & 0 < x < L, t >0 \\ u_x(0, t) = 0, & t \geq 0 \\ u_x(L,t) = \gamma, & t \geq 0 \\ u(x,0)=f(x), & 0 < x < L \end{cases}$$
第二步：边界条件齐次化
1. 特征值问题，针对齐次诺伊曼条件，即：对于齐次的第二类边界条件（即$u_x(0,t)=0,u_x(L,t)=0$）设$u(x,t)=X(x)T(t)$，代入方程得：
$$ X'' + \lambda X = 0, X'(0) = 0, X'(L) = 0$$
解特征值问题：
- 当$\lambda = 0$时,方程退化为：
- 当$\lambda > 0$时：$X_n(x)=\cos\left(\frac{n\pi x}{L}\right), \lambda_n = \left(\frac{\pi x}{L}\right)^2, n=1,2,3,\dots$
$$X'' = 0$$积分两次：
$$X'(x) = C_1$$
$$X(x)=C_1x + C_2$$
应用边界条件: $X'(0)=0 \implies C_1 = 0, X'(L)=0 \implies C_1=0$因此：
$$ X(x) = C_2 $$
其中$C_2$是任意非零常数。通常我们取最简单的形式，即$C_2=1$$，所以：
$$X_0(x)=1$$
2. 时间部分，对每个$\lambda_n$：
$$ T_n(t) = C_ne^{-\alpha \lambda_n t}$$
特别地，对于$\lambda_0=0$:
$$ T_0(t)=C_0$$
3.通解
$$u(x,t)=C_0 + \sum_{n=1}^{\infty}C_n \cos\left( \frac{n\pi x}{L}\right)e^{-\alpha\left(\frac{n\pi}{L}\right)^2t}$$
4. 利用初始条件,当$t=0$时
$$ f(x) = C_0 + \sum_{n=1}^{\infty}C_n \cos\left(\frac{n\pi x}{L} \right)$$
这是函数 $f(x)$ 在 $[0,L]$ 上的傅里叶余弦级数展开。系数为：
$$C_0 = \frac{1}{L}\int_{0}^L f(x)dx$$
$$C_n = \frac{2}{L}\int_{0}^L f(x)cos\left( \frac{n\pi x}{L} \right)dx, n=1,2,3,\dots$$
第三步：处理非齐次边界条件（恒定热流）
回到原问题$u_x(0,t)=0,u_x(L,t)=\gamma$,设：
$$u(x,t) = v(x,t) + \gamma x + \beta t$$
我们选择 $\beta$使得 $v$ 满足齐次边界条件。
$$u_x = v_x + \gamma$$
$$u_{xx} = v_{xx}$$
$$u_t = v_t + \beta $$
代入方程$u_t=\alpha u_{xx}$
$$v_t + \beta = \alpha v_{xx} \implies v_t = \alpha v_{xx} - \beta $$
三、完整结果（恒定热流输入）
对于问题：
$$ u_x(0,t) = 0, u_x(L,t) = \gamma, u(x,0)=f(x)$$
解为：
$$u(x,t) = \frac{1}{L} \int_{0}^L f(\xi)d\xi + \frac{\gamma}{L}t + \sum_{n=1}^{\infty} C_n \cos\left( \frac{n\pi x}{L}\right)e^{-\alpha (\frac{n\pi}{L})^2t}$$
其中：
$$C_n = \frac{2}{L}\int_{0}^L f(\xi) \cos\left( \frac{n\pi \xi}{L}\right)d\xi - \frac{2\gamma L}{n^2\pi^2}[(-1)^n-1]$$
四、两类边界的区别
| 第一类（固定温度） | 第二类（绝热） | |
|---|---|---|
| 边界行为 | 温度被强制固定 | 温度自由变化，但梯度为零 |
| 能量交换 | 热量可以自由进出 | 热量无法穿过边界 |
| 长期行为 | 温度 → 边界温度（通常非零） | 温度 → 初始平均值（常数） |
| 总热量 | 不守恒 | 守恒 |
五、最一般的问题表述
边界条件（第二类，诺伊曼型）：
$$u_x(c, t) = a, \quad u_x(d, t) = b$$
设原问题定义在 $x \in [c,d]$上
$$ u(x, t_0) = f(x) $$
$$ u_t = \alpha u_{xx},\quad c < x < d, \quad t > t_0$$
其中 $t_0$是任意给定的初始时刻（不一定为 0）。
统一求解公式
第 1 步：计算整体温度变化率
$$ \beta = \frac{b-a}{L}, \quad L = d-c$$
这个 $\beta$ 是平均温度随时间的变化率。
最终统一公式
$$ u(x,t) = \frac{\beta}{2}(x-c)^2 + a(x-c) + \beta(t-t_0) + C_0 + \sum_{n=1}^{\infty}C_n \cos\left( \frac{n\pi (x-c)}{L} \right) e^{-\alpha \left(\frac{n\pi}{L} \right)^2(t-t_0)}$$
其中
$$L = d - c, \quad \beta = \frac{b-a}{L}$$
$$ C_0 = \frac{1}{L} \int_{c}^{d} \left[ f(\xi) - \frac{\beta}{2}(\xi -c)^2 - a (\xi -c) \right] d\xi $$
$$ C_n = \frac{2}{L} \int_{c}^{d} \left[ f(\xi) - \frac{\beta}{2}(\xi -c)^2 - a (\xi -c) \right] \cos\left( \frac{n\pi(\xi -c)}{L}\right) d\xi $$
六、用傅里叶定律判断热流方向
对于一维问题，热流密度 $q$ 遵循傅里叶定律：
$$ q(x,t) = -k \frac{\partial u}{\partial x}(x,t)$$
其中$k>0$是导热系数。 $q>0$ 表示热量向$x$正方向（向右）流动。因此，要判断边界上的热流是流入还是流出杆，我们需要看边界上的温度梯度$u_x$。
- 左边界 ( $x=0$)：
- 如果$u_x(0,t)>0$：温度从左边界向右升高。那么$q(0,t)=-k \cdot u_x < 0$，即热流向左。由于左边界的外法向指向左，负的热流意味着热量从杆内流出到左端环境。
- 如果$u_x(0,t)<0$：温度从左边界向右降低。那么$q(0,t)=-k \cdot u_x > 0$，即热流向右。这表示热量从左端环境流入杆内。
- 右边界 ( $x=L$)：
- 如果$u_x(L,t)>0$：表示从 $x < L$ 到 $x=L$ 温度是升高的。那么$q(L,t) = -k \cdot u_x < 0$，即热流向左。由于右边界的外法向指向右，负的热流意味着热量从外界流入杆内。
- 如果$u_x(L,t)<0$：温度从右边界向左降低（即从内部到边界温度下降）。那么$q(L,t)=−k \cdot u_x>0$，即热流向右。这表示热量从杆内流出到右端环境。
例如：参数 $u_x(0,t)=20, u_x(1,t)=-30, u(x,0) = 10$描述的是：一个初始温度为 $10^{\circ}$的杆，其两个端面都被强制以恒定的速率向外散热（左端散热速率20，右端散热速率30）。因此，整个杆的温度会持续下降，且由于右端散热更快，杆的温度分布会呈现出右侧比左侧更冷的趋势。下图横坐标表示空间位置$x$，纵坐标表示时间$t$。可见随着时间的推进，杆的温度持续降低。
