主题
Search

热传导方程--圆盘


为了求解半径为 a=1 的二维圆盘上的热传导方程,尝试使用分离变量法:

 U(r,theta,t)=R(r)Theta(theta)T(t).
(1)

柱坐标系表示 thetar 项的拉普拉斯算子,得到:

 del ^2=(d^2R)/(dr^2)+1/r(dR)/(dr)+1/(r^2)(d^2Theta)/(dtheta^2),
(2)

因此,热传导方程变为:

 (RTheta)/kappa(dT)/(dt)=(d^2R)/(dr^2)ThetaT+1/r(dR)/(dr)ThetaT+1/(r^2)(d^2Theta)/(dtheta^2)RT.
(3)

等式两边同乘 r^2/RThetaT 得到:

 (r^2)/(kappaT)(dT)/(dt)=(r^2)/R(d^2R)/(dr^2)+r/R(dR)/(dr)+(d^2Theta)/(dtheta^2)1/Theta.
(4)

theta 项可以被分离。

 (d^2Theta)/(dtheta^2)1/Theta=-n(n+1),
(5)

其解为:

 Theta(theta)=Acos[sqrt(n(n+1))theta]+Bsin[sqrt(n(n+1))theta].
(6)

剩余部分变为:

 (r^2)/(kappaT)(dT)/(dt)=(r^2)/R(d^2R)/(dr^2)+r/R(dR)/(dr)-n(n+1).
(7)

等式两边除以 r^2 得到:

 1/(kappaT)(dT)/(dt)=1/R(d^2R)/(dr^2)+1/(rR)(dR)/(dr)-(n(n+1))/(r^2)=-1/(lambda^2),
(8)

这里选择了一个的分离常数,以便 t 部分保持有限:

 T(t)=Ce^(-kappat/lambda^2).
(9)

径向部分变为:

 1/R(d^2R)/(dr^2)+1/(rR)(dR)/(dr)-(n(n+1))/(r^2)+1/(lambda^2)=0
(10)
 r^2(d^2R)/(dr^2)+r(dR)/(dr)+[(r^2)/(lambda^2)-n(n+1)]R=0,
(11)

这就是球贝塞尔微分方程

考虑半径为 a 的圆盘,初始温度为 U(r,0)=0边界条件U(a,t)=1。那么解为:

 U(r,t)=1-2sum_(n=1)^infty(J_0((alpha_nr)/a))/(alpha_nJ_1(alpha_n))e^(-alpha_n^2kappat/a^2),
(12)

其中 alpha_n第一类贝塞尔函数 J_0(x) 的第 n零点 (Bowman 1958, pp. 37-39)。


另请参阅

热传导方程

使用 Wolfram|Alpha 探索

参考文献

Bowman, F. Introduction to Bessel Functions. New York: Dover, 1958.Carslaw, H. S. and Jaeger, J. C. "Some Two-Dimensional Problems in Conduction of Heat with Circular Symmetry." Proc. London Math. Soc. 46, 361-388, 1940.

请引用为

Weisstein, Eric W. "热传导方程--圆盘。" 来自 MathWorld--Wolfram Web 资源。 https://mathworld.net.cn/HeatConductionEquationDisk.html

主题分类