一、假设解的形式

热方程设$u(x,t)=X(x)T(t)$(空间函数 $X(x)$ 与时间函数 $T(t)$ 分离)

拉普拉斯方程设$\phi(x,y)=X(x)Y(y)$(x 方向与 y 方向函数分离)

1D Heat Equation

大部分Heat Equation是1D的,如杆的温度分布.

$$\frac{\partial u}{\partial t}=k\cdot \frac{\partial^2u}{\partial x^2}$$

$$\frac{\partial^2 u}{\partial x^2} = \frac{1}{k^2} \cdot \frac{\partial u}{\partial t}$$

2D Heat Equation

若温度在平面(如薄板)上变化,需考虑 x、y 两个空间方向,形式为:

$$\frac{\partial u}{\partial t} = k \cdot \left( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} \right)\$$

二、分解为ODEs

将假设解代入PDE,利用 “等式两边仅依赖不同变量” 推出两边等于同一常数 $\alpha$(称为 Eigenvalue)

热方程会分解为 $X^{”} – \alpha X = 0$ 和 $T’ – k^2\alpha T = 0$(k 为 Thermal Diffusivity)

三、Non-physical Solutions

分 $\alpha=0, \alpha>0, \alpha<0$ 三种情况求解 ODEs,结合物理意义筛选:

热方程中 $\alpha>0$ 会导致 $T(t)$ 指数增长(不符合温度衰减),排除;

拉普拉斯方程中半无限区域(如 $y \to \infty$)的 $\alpha>0$ 解无界(不符合稳态温度有限),排除

四、叠加形成级数解

收集 $\alpha<0$ 对应的物理解(称为 Eigenfunctions,如 $\sin(n\pi x/L), \cos(n\pi x/L)$,通过 Superposition Principle 组成 Series Solution


0 条评论

发表回复

Avatar placeholder

您的邮箱地址不会被公开。 必填项已用 * 标注