跳到主要内容

東京工業大学 環境・社会理工学院 土木・環境工学系 2021年8月実施 午前 基礎科目 問題2

Author​

思齐塾, 祭音Myyura

Description​

  1. 0≤x≤π,t≥00 \le x \le \pi, t \ge 0 で定義された次の偏微分方程式 (1) を,与えられた境界条件・初期条件 (2) のもと,変数分離法により解け.
∂u(x,t)∂t=4∂2u(x,t)∂x2(1)\frac{\partial u(x, t)}{\partial t} = 4 \frac{\partial^2 u(x, t)}{\partial x^2} \tag{1}
u(0,t)=0,u(π,t)=0,u(x,0)=2sin⁡3x−4sin⁡5x(2)u(0, t) = 0, \quad u(\pi, t) = 0, \quad u(x, 0) = 2 \sin 3x - 4 \sin 5x \tag{2}

题目描述​

在区域

0≤x≤π,t≥00\leq x\leq\pi,\qquad t\geq0

上,使用分离变量法求解热方程

∂u(x,t)∂t=4∂2u(x,t)∂x2,(1)\frac{\partial u(x,t)}{\partial t} =4\frac{\partial^2u(x,t)}{\partial x^2}, \tag{1}

并满足齐次固定端边界条件与初始条件

u(0,t)=0,u(π,t)=0,u(x,0)=2sin⁡3x−4sin⁡5x.(2)u(0,t)=0,\qquad u(\pi,t)=0,\qquad u(x,0)=2\sin3x-4\sin5x. \tag{2}

Kai​

2.​

境界条件は u(0,t)=u(π,t)=0u(0,t)=u(\pi,t)=0 と読む。 u(x,t)=X(x)T(t)u(x,t)=X(x)T(t) とおくと、

XT′=4X′′TXT'=4X''T

より

T′4T=X′′X=−λ\frac{T'}{4T}=\frac{X''}{X}=-\lambda

である。したがって

X′′+λX=0,T′+4λT=0.X''+\lambda X=0,\qquad T'+4\lambda T=0.

境界条件 X(0)=X(π)=0X(0)=X(\pi)=0 を満たす非自明解は

λ=n2,Xn(x)=sin⁡nx(n=1,2,…)\lambda=n^2,\qquad X_n(x)=\sin nx\qquad(n=1,2,\ldots)

であり、対応する時間因子は Tn(t)=e−4n2tT_n(t)=e^{-4n^2t} である。よって一般に

u(x,t)=∑n=1∞cne−4n2tsin⁡nx.u(x,t)=\sum_{n=1}^{\infty}c_ne^{-4n^2t}\sin nx.

初期条件

u(x,0)=2sin⁡3x−4sin⁡5xu(x,0)=2\sin3x-4\sin5x

から、直交性により c3=2, c5=−4c_3=2,\ c_5=-4 、その他の cnc_n は 00 である。したがって

u(x,t)=2e−36tsin⁡3x−4e−100tsin⁡5x\boxed{u(x,t)=2e^{-36t}\sin3x-4e^{-100t}\sin5x}

を得る。各項は ut=4uxxu_t=4u_{xx} を満たし、 x=0,πx=0,\pi で消え、 t=0t=0 で所定の初期関数になるので、すべての条件を満たしている。