热传导方程的第二边值问题

第二边值问题,也称为 Newmann 边值问题。它是函数在边界上没有变化量的条件,在热传导方程的实际问题里,就是在边界上没有热量的交换。在数学上的表示,是函数在边界上的导数为 \(0\)。

第二边值问题的分离变量法,与第一边值的分离变量法,思想是一样的。就是将假设未知函数可以写成变量分离的形式,然后代入到方程里去,从而求出解的方法。

1,考虑初边值问题:

\begin{cases}u_t=c^2u_{xx},\quad & 0<x<L, t>0\\ u(x,0)=\phi(x),& 0<x<L\\ \frac{\partial u}{\partial x}(0,t)=0,\frac{\partial u}{\partial x}(L,t)=0,&t>0\end{cases}

这样的边界条件称为第二边界条件,也称为 Newman 边界条件。

2,分离变量法:我们设 \(u(x,t)=T(t)X(x)\),代入方程, 我们得到

\[T'(t)X(x)=c^2T(x)X^{\prime\prime}(x)\]

两边除以 \(c^2T(t)X(x)\),得到 \[\frac{T'(t)}{c^2T(t)}=\frac{X^{\prime\prime}(x)}{X(x)}\]

因为左边只跟 \(t\) 有关,右边只跟 \(x\) 有关,要使得它们相等,只能是两边都是常数,我们设这个常数为 \(-\lambda\),也就是

\[\frac{T'(t)}{c^2T(t)}=\frac{X^{\prime\prime}(x)}{X(x)}=-\lambda\]

现在我们得到了两个方程

\[\frac{T'(t)}{c^2T(t)}=-\lambda\]

\[\frac{X^{\prime\prime}(x)}{X(x)}=-\lambda\]

3,解 \(X(x)\)。由原方程的边界条件,我们得到

\begin{cases}X^{\prime\prime}(x)+\lambda X(x)=0\\ X'(0)=0, X'(L)=0\end{cases}

现在我们考虑三种情况:

(1)\(\lambda<0\),由二阶常微分方程理论,我们得到

\[X(x)=C_1e^{\sqrt{-\lambda}x}+C_2e^{-\sqrt{-\lambda}x}\]

代入边界条件,我们得到 \(C_1=C_2=0\),所以 \(\lambda <0\) 不是特征值。

(2)\(\lambda=0\),我们得到 \[X^{\prime\prime}(x)=0\]

积分两次,得到 \[X(x)=Ax+B\]

代入边界条件,得到 \(B=0\),所以 \(X=A\) 是特征函数。

(3)\(\lambda>0\),由二阶常微分方程理论,我们得到

\[X(x)=C_1\cos\sqrt{\lambda}x+C_2\sin\sqrt{\lambda}x\]

代入边界条件,我们得到

\[c_2=0, \sqrt{\lambda}=\frac{n\pi}{l}, \lambda=\left(\frac{n\pi}{l}\right)^2, n=1,2,\cdots\]

我们记 \[X_0=A, X_n=\cos\frac{n\pi x}{l}\]

它们是特征函数。

4,解 \(T(t)\)。\(T(t)\) 满足的方程是

\[\frac{T'(t)}{c^2T(t)}=-\lambda\]

当 \(\lambda=0\) 时,\[\frac{T_0′(t)}{T_0(t)}=0\]

两边积分,得到 \[\ln |T_0(t)|=C\] 再取对数,得到 \(T_0(t)=C\)

当 \(\lambda=\left(\frac{n\pi}{l}\right)^2\) 时,

\[\frac{T_n'(t)}{c^2T_n(t)}=-\left(\frac{cn\pi}{l}\right)^2\]

由常微分方程的分离变量法,它的解为

\[T_n(t)=A_ne^{-\left(\frac{n\pi}{l}\right)^2t}\]

5,原方程的解:由齐次线性微分方程的叠加原理,\(u(x,t)\) 的表达式为

\begin{align*}u(x,t)&=A_0T_0(t)X_0(x)+\sum_{n=1}^{\infty}A_nT_n(t)X_n(x)\\ &=A_0+\sum_{n=1}^{\infty}A_ne^{-\left(\frac{cn\pi}{l}\right)^2t}\cos\frac{n\pi x}{l}\end{align*}

代入初值条件,

\[u(x,0)=\phi(x)=A_0+\sum_{n=1}^{\infty}A_n\cos\frac{n\pi x}{l}\]

等式右边是傅里叶余弦级数,所以它的系数是

\[A_0=\frac{1}{L}\int_0^L\phi(x)dx,\]

\[A_n=\frac{2}{L}\int_0^L\phi(x)\cos\frac{n\pi x}{l}dx\]

所以原方程的解为

\[u(x,t)=\frac{1}{L}\int_0^L\phi(x)dx+\sum_{n=1}^{\infty}A_ne^{-\left(\frac{cn\pi}{l}\right)^2t}\cos\frac{n\pi x}{l}\]

其中

\[A_n=\frac{2}{L}\int_0^L\phi(x)\cos\frac{n\pi x}{l}dx\]