环形区域拉普拉斯方程傅里叶级数解
题目
问题 2. 在环形区域 { (r,θ):1<r≤2, −π≤θ<π} \{(r, \theta) : 1 < r \leq 2, \, -\pi \leq \theta < \pi\} { (r,θ):1<r≤2,−π≤θ<π} 中求解
Δu=0,\Delta u = 0,Δu=0,
u∣r=1=0,u|_{r=1} = 0,u∣r=1=0,
u∣r=2=θ2.u|_{r=2} = \theta^2.u∣r=2=θ2.
解应以适当的傅里叶级数形式表示。
解答
该问题为拉普拉斯方程在环形区域上的 Dirichlet 边值问题。区域为内半径 r=1 r=1 r=1、外半径 r=2 r=2 r=2 的圆环,边界条件为:
- 在内边界 r=1 r=1 r=1 上,u=0 u = 0 u=0;
- 在外边界 r=2 r=2 r=2 上,u=θ2 u = \theta^2 u=θ2.
拉普拉斯方程在极坐标下为:
Δu=∂2u∂r2+1r∂u∂r+1r2∂2u∂θ2=0. \Delta u = \frac{\partial^2 u}{\partial r^2} + \frac{1}{r} \frac{\partial u}{\partial r} + \frac{1}{r^2} \frac{\partial^2 u}{\partial \theta^2} = 0. Δu=∂r2∂2u+r1∂r∂u+r21∂θ2∂2u=0.
由于边界条件关于 θ \theta θ 为偶函数(θ2 \theta^2 θ2 是偶函数),解中仅含余弦项。使用分离变量法,假设解的形式为 u(r,θ)=R(r)Θ(θ) u(r, \theta) = R(r) \Theta(\theta) u(r,θ)=R(r)Θ(θ)。分离变量后,得到角度部分的解为 Θn(θ)=cos(nθ) \Theta_n(\theta) = \cos(n\theta) Θn(θ)=cos(nθ),径向部分的解为:
- 当 n=0 n = 0 n=0 时,R0(r)=A0+B0lnr R_0(r) = A_0 + B_0 \ln r R0(r)=A0+B0lnr;
- 当 n≥1 n \geq 1 n≥1 时,Rn(r)=Anrn+Bnr−n R_n(r) = A_n r^n + B_n r^{-n} Rn(r)=Anrn+Bnr−n.
因此,通解为:
u(r,θ)=A0+B0lnr+∑n=1∞(Anrn+Bnr−n)cos(nθ). u(r, \theta) = A_0 + B_0 \ln r + \sum_{n=1}^{\infty} (A_n r^n + B_n r^{-n}) \cos(n\theta). u(r,θ)=A0+B0lnr+n=1∑∞(Anrn+Bnr−n)cos(nθ).
应用边界条件
-
内边界条件 u∣r=1=0 u|_{r=1} = 0 u∣r=1