x=0 で 等温境界条件 u(0,t)=uをみたし、 x=1 で 断熱境界条件 ux(1,t)=0をみたす解 を求めるのに必要な道具
※数式がスマホで画面からはみ出る場合、横スクロールするかピンチインしてください。
目次
x=0で等温境界条件u(0,t)=uをみたし、x=1で断熱境界条件 ux(1,t)=0をみたす解を求めるのに必要な道具
熱伝導を表す偏微分方程式
u(x,t)=X(x)T(t)とおき\(u_t=\alpha^2u_{xx}\)にあてはめると
$$X(x)T'(t)=\alpha^2X”(x)T(t),\ \frac{T'(t)}{\alpha^2T(t)}=\frac{X”(x)}{X(x)}$$
物理的な制約(熱の出入りがない状況)からk<0と仮定する。
\(\frac{T'(t)}{\alpha^2T(t)}=\frac{X”(x)}{X(x)}=k\)とおくと\(T'(t)-k\alpha^2T(t)=0,\ X”(x)-kX(x)=0\)が得られる。
\(k=-\lambda^2\)とすれば\(T'(t)+\lambda^2\alpha^2T(t)=0,\ X”(x)+\lambda^2X(x)=0\)
一般解
$$T(t)=A_1e^{-\lambda^2\alpha^2t}$$
$$X(x)=B_1\ sin\lambda x+B_2\ cos\lambda x=0$$
\(u(x,t)=T(t)X(x)\)であるから、\(A=A_1B_1,\ B=A_1B_2\)と取り直せば
$$u(x,t)=A_1e^{-\lambda^2\alpha^2t}(B_1sin\lambda x+B_2cos\lambda x)=e^{-\lambda^2\alpha^2t}(Asin\lambda x+Bcos\lambda x)$$
重ね合わせの原理
線形同次微分方程式では\(x_1,x_2\)が解ならば、その1次結合\(c_1x_1+c_2x_2\)も解となる。
証明
Lはxに関数L[x]を対応させる写像 L[x]=0
\(x_1,\ x_2\)が同次微分方程式の解であるから\(L[x_1]=0,\ L[x_2]=0\) よって\(L[c_1x_1+c_2x_2]=c_1L[x_1]+c_2L[x_2]=0\)
x=0で等温境界条件u(0,t)=uをみたし、x=1で断熱境界条件 ux(1,t)=0をみたす解
u(0,t)=u_x(1,t)=0をみたす解
等温境界条件
$$u(0,t)=e^{-\lambda^2\alpha^2t}(Asin0+Bcos0)=e^{-\lambda^2\alpha^2t}B=0,\ B=0$$
断熱境界条件
\(u(x,t)=e^{-\lambda^2\alpha^2t}Asin\lambda x\)をxで偏微分すると\(u_x(x,t)=e^{-\lambda^2\alpha^2t}A\lambda\ cos\ \lambda x\)
$$u_x(1,t)=e^{-\lambda^2\alpha^2t}A\lambda\ cos\ \lambda=0,\ A=0の場合は自明解,\ cos\lambda=0,\ \lambda=\frac{\pi}{2}+n\pi=\frac{\pi+2n\pi}{2}=\frac{(2n+1)\pi}{2}$$
これらの解を重ね合わせて\(u(0,t)=u_x(1,t)=0\)をみたす解は
$$u(x,t)=\sum_{n=0}^\infty A_n \exp\left(-\left(\frac{(2n+1)\pi}{2}\right)^2\alpha^2t\right)sin\left(\frac{(2n+1)\pi x}{2}\right)$$
参考文献
検索しても計算過程が見つからない場合
検索しても計算過程が見つからない場合ココナラ を利用してみてはいかがでしょうか。
ココナラ 登録方法
会員登録しなくてもサービスの検索はできます。
サービスの購入・出品には会員登録が必要です。
スタディサプリ進路 社会人向け の 使い方
スタディサプリ進路 社会人向けで社会人が数学を学べる大学を検索してみます。
独学で大学数学の微分方程式を勉強しています!
- 1/sinxdxから1/tdtへの変形 トラクトリックス
- 1階非斉次線形微分方程式の一般解
- y’=(1+y)/sinxの解き方 変数分離形
- y’=(x^2-y^2)/2xyの解き方 同次形
- y’=2y/x-yの解き方 同次形
- dy/dx-(3/2)(y-a)^(1/3)=0 の一般解と、それらの解曲線の包絡線である特異解
- x=-e^2t∫3t^2e^(-2t)dt+te^2t∫3te^(-2t)dtの解き方 部分積分
- 特殊解を求めるのに 定数変化法 より クラメルの公式
- x”-x’=sint+2costの一般解(未定係数法)
- x”-2x’+5x=20cost, x(0)=x'(0)=0の解き方 初期値問題
- (t+2)x”-(2t+6)x’+(t+4)x=0 (x=e^t)の一般解 階数低下法
- (t^2+3t+4)x”+(t^2+t+1)x’-(2t+3)x=0 (x=e^(-t))の一般解 階数低下法
- 3つの関数の積の積分
- 4xy”+2y’+y=0の解き方(オイラーの微分方程式)
- 4xy”+2y’+y=0の解き方(フロベニウスの方法)
- x^2y”-2xy’+(x^2+2)y=0の解き方(フロベニウスの方法)
- 級数解法・フロベニウスの方法 使い分け
- x=0 で 等温境界条件 u(0,t)=uをみたし、 x=1 で 断熱境界条件 ux(1,t)=0をみたす解
- x=0 で 断熱境界条件 ux(0,t)=uをみたし、 x=1 で 等温境界条件 u(1,t)=0 をみたす解
- 4x^7y-28x^5y^3+28x^3y^5-4xy^7をr^nsinnθ,r^ncosnθの1次結合として表す。
- 4x^7y-28x^5y^3+28x^3y^5-4xy^7の最大値・最小値 ラプラシアン
- 独学で大学数学の積分因子を勉強しています!
- 完全微分形の一般解
- (x^2+3xy+2y^2)dy+(2x^2+3xy+y^2)dx=0の一般解 完全微分形