next up previous
Next: 3 現実的な恒星内部構造モデル Up: 8 常微分方程式の数値解法:境界値問題 Previous: 1 シューティング法

Subsections


2 境界値問題としてのLane-Emden 方程式

Lane-Emden方程式([*])を境界値問題として解いてみよう。 適当なフィッティングポイント$\xi_{\rm f}$を設け、 中心からフィッティングポイントまでの解は前節で解説した通りの方法で求め、 その解に、表面、即ち$\theta(\xi)=0$となる $\xi=\xi_{\rm s}$から 積分してきた解が、このフィッティングポイントで一致するようにすれば良い。

表面からの積分は$\xi_{\rm s}$の値を決めた上で行わなければならないのは言うまでもない が、微分方程式は2階だから、 $\xi=\xi_{\rm s}$での初期条件が2つ必要である。 その1つは $\theta(\xi)=0$だが、もう1つの$d\theta/d\xi$は最初はわからない。 そこで、$\xi_{\rm s}$とそこでの$d\theta/d\xi$の値をパラメーターとして与えてやれば、 外側の解をこの2つのパラメーターの関数として扱う事が出来る。 $y_1 = \theta$ $y_2 = \xi^2 d\theta/d\xi$とし、 $y_1^{{\rm (int)}}$ $y_2^{{\rm (int)}}$を中心からの積分で求めた $\xi=\xi_{\rm f}$での値、 $y_1^{{\rm (ext)}}$ $y_2^{{\rm (ext)}}$を表面からの積分で求めた $\xi=\xi_{\rm f}$での値としよう。内側の解はユニークに決まっている。 一方、外側の解は $x \equiv \xi_{\rm s}$と そこでの $z \equiv d\theta/d\xi$の値に依存している。 全体の解となるには、 $\xi=\xi_{\rm f}$

\begin{displaymath}
Y(x,z) \equiv
y_1^{{\rm (int)}}-y_1^{{\rm (ext)}}(x,z)
=0
\end{displaymath} (9)


\begin{displaymath}
Z(x,z) \equiv
y_2^{{\rm (int)}}-y_2^{{\rm (ext)}}(x,z)
=0
\end{displaymath} (10)

が満たされなければならない。

手順としては、 先ず、適当に選んだ$x_{\rm t}$$z_{\rm t}$に対する$Y$$Z$を 求める。 次に、$z$の値を変えずに、$x$$x_{\rm t}$から $x_{\rm t}+\delta x$に変えて、 $Y$及び$Z$の変化量を求める。 同様に、$x$の値を変えずに、$z$$z_{\rm t}$から $z_{\rm t}+\delta z$に変えて、 $Y$及び$Z$の変化量を求める。 これから

\begin{displaymath}
\partial Y/\partial x =
[Y(x_{\rm t}+\delta x,z_{\rm t})-Y(x_{\rm t},z_{\rm t})]/\delta x
\end{displaymath} (11)


\begin{displaymath}
\partial Y/\partial z =
[Y(x_{\rm t},z_{\rm t}+\delta z)-Y(x_{\rm t},z_{\rm t})]/\delta z
\end{displaymath} (12)

を数値的に求める。 今、$x_{\rm t}$$z_{\rm t}$の値が、真の解$x_0$$z_0$から $\Delta x$$\Delta z$ずれているとすると、1次の精度で、
\begin{displaymath}
Y(x_{\rm t},z_{\rm t}) +
(\partial Y/\partial x)\Delta x + (\partial Y/\partial z)\Delta z =0
\end{displaymath} (13)


\begin{displaymath}
Z(x_{\rm t},z_{\rm t}) +
(\partial Z/\partial x)\Delta x + (\partial Z/\partial z)\Delta z =0
\end{displaymath} (14)

となる。そこで、上の2式から$\Delta x$$\Delta z$を求め、 $x=x_{\rm t}+\Delta x$を改めて$x_{\rm t}$ $z=z_{\rm t}+\Delta z$を改めて$z_{\rm t}$として、 $\Delta x$$\Delta z$が充分小さくなるまで、逐次計算を繰り返せば良い。

1 問題

境界値問題として、Lane-Emden方程式を解け。



Jun Makino
平成15年6月5日