- 追加された行はこの色です。
- 削除された行はこの色です。
#rst2hooktail_source
=============================================
流れ関数
=============================================
ベクトル場とそのポテンシャル関数に関する一つの応用として、流体力学における二次元平面流を考えます。このベクトル場 $\bm{V}$ は管状であるとします。（ $\nabla$ も二次元であることに注意して下さい。）
非圧縮流の連続の式
====================================================================
いま、流れ場 $\bm{V}$ の中に適当な流管を取ります。質量保存則より『流管のどの断面を取っても、単位時間に断面を過ぎる質量は等しい』と考えることが出来ます。これを式にすると次のようになります。流体の密度を $\rho$ とします。
<tex>
\nabla \cdot (\rho \bm{V})= 0 \tag{1}
</tex>
流管の考え方については ベクトル場の流線と流管_ を参照して下さい。式 $(1)$ は、流れ場が連続であれば必ず成り立つことが要請されるもので、 *連続の式* と呼ばれます。要するに質量保存則なわけですが、流体力学の大事な基礎方程式の一つです。特に、非圧縮流の場合、 $\rho =const.$ ですので、 $\rho$ を外に出して次のように簡単化できます。
<tex>
\nabla \cdot \bm{V} = 0 \tag{2}
</tex>
これは *非圧縮流の連続の式* ですが、同時に $\bm{V}$ が管状であることも示しています。質量保存の法則はどんな流体にも成り立ちますから、流体力学の立場では、管状ベクトル場とは非圧縮流の流れを意味すると考えても良さそうです。
二次元流と流れ関数
================================================================
前セクションの内容は三次元の流れでも成り立つものでしたが、ここでは特に二次元平面流を考えます。 $\nabla$ も二次元になることに注意して下さい。
<tex>
\nabla \cdot \bm{V} = \frac{\partial V_{x}}{\partial x} + \frac{\partial V_{y}}{\partial y} = 0 \tag{3}
</tex>
このとき、 $\bm{V}$ の成分を、スカラーポテンシャル $\Psi$ を用いて次のように表わすことができます。
<tex>
V_{x} = \frac{\partial \Psi}{\partial y}, \ \ \ \ V_{y} = - \frac{\partial \Psi}{\partial x} \tag{4}
</tex>
この $\Psi$ を *流れ関数* と呼びます。式 $(4)$ が自動的に式 $(3)$ を満たすことを確認してみて下さい。流れ関数は、次のようにベクトルポテンシャルの形で与えることもできます。この方が、すっきりと覚えやすいかも知れません。（管状ベクトル場にベクトルポテンシャルが存在することについては、 管状ベクトル場_ を参照して下さい。）
<tex>
\bm{V} = \nabla \times
\left( \begin{array}{c}
0 \\
0 \\
\Psi \\
\end{array}
\right) \tag{5}
</tex>
もしくは、次のように表現することも可能です。
<tex>
\bm{V} = \nabla \Psi \times
\left( \begin{array}{c}
0 \\
0 \\
1 \\
\end{array}
\right) \tag{6}
</tex>
流線と流れ関数の関係
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
ついでに、流れ関数の物理的意味も考えておきます。 $\bm{V}$ の中に置かれた曲線 $OP$ を単位時間当たりに通過する全流量 $Q$ を求めてみます。 $\bm{V}$ は二次元平面流なので、
<tex>
Q & = \int _{O}^{P} (\bm{V} \cdot \bm{n}) ds \\
&= \int _{O}^{P} (V_{x}n_{x}+V_{y}n_{y}) ds \\
&= \int _{O}^{P} \left( \frac{\partial \Psi}{\partial y}\frac{dy}{ds} +
\frac{\partial \Psi}{\partial x}\frac{dx}{ds} \right) ds \\
&= \int_{O}^{P}\frac{\partial \Psi}{\partial s}ds \\
& = \int_{O}^{P}d\Psi \\
&= \Psi (P) -\Psi (O) \tag{7}
</tex>
これより次のことが言えます。
.. admonition:: theorem
二次元非圧縮流の流れの中に二点 $O,P$ を取るとき、 $O,P$ を結ぶ任意の曲線を単位時間当たりに通過する流量は流れ関数の差 $\Psi (P) -\Psi (O)$ で表わされ、曲線の形にはよらない。
また、同一の流線上では流量が零なので、 $\Psi = const.$ が言えます。これも重要な性質です。流れ関数の物理的意味が、少し明らかになって来たでしょうか。
.. admonition:: theorem
同一の流線の上では、 $\Psi = const.$ が成り立ちます。逆に $\Psi = const.$ は、一本の流線を表わします。
ベクトル場をポテンシャル関数で表現する一例として、二次元非圧縮流の流れ関数を取り上げてみました。より詳しくは、きちんと流体力学の教科書を参考にして下さい。
ノイマン問題
======================================================================
流れ関数は、どのような境界条件を満たすのでしょうか？物理的状況を考えてみると、水の流れは、境界（つまり、水槽の壁や川底などの固体壁）を突き抜けて流れることはできません。境界では、必ず境界に沿った方向（接線方向）に流れ、境界の法線方向の速度成分は持たない、ということです。これを式で書くと次のように表わせます。
<tex>
\frac{\partial \psi}{\partial n}=0 \tag{8}
</tex>
ここで、もし流れ場が 層状ベクトル場_ だとすると、 $\nabla \times \bm{V}=\bm{0}$ が要請されます。流れ関数の表わす流れが二次元平面流であることに注意すると、この層状の条件は次のようになります。
<tex>
\nabla \times \bm{V} &=
\left( \begin{array}{c}
\frac{\partial}{\partial x} \\
\frac{\partial}{\partial y} \\
\frac{\partial}{\partial z} \\
\end{array}
\right)
\times
\left( \begin{array}{c}
\frac{\partial \psi}{\partial y} \\
-\frac{\partial \psi}{\partial x} \\
1 \\
\end{array}
\right) \\
&=
\left( \begin{array}{c}
0 \\
0 \\
-\frac{\partial^2 \psi}{\partial x^2}-\frac{\partial^2 \psi}{\partial y^2} \\
\end{array}
\right) \\
&= \bm{0} \tag{9}
</tex>
式 $(9)$ の $z$ 成分だけ考えて、次式に帰着します。ただし、 $\triangle_{2}$ は二次元ラプラシアンの意味だとします。
<tex>
\frac{\partial^2 \psi}{\partial x^2}+\frac{\partial^2 \psi}{\partial y^2}=\triangle_{2}\psi=0 \tag{10}
</tex>
式 $(8)(10)$ より、二次元層状ベクトル場の流れ関数は、ノイマン問題の解として与えられるということが分かりました。（ノイマン問題の定義は ラプラス場_ の最後に紹介しています。）このようにして、流れ場を求める問題が、微分方程式を境界条件の下に解くという問題に帰着してしまったわけです。流れの様子を生き生きと想像していた人にとって、もしかして微分方程式を解くという作業は、とても静的で無機質に思えるかも知れません。しかし、このようにして流れを微分方程式に帰着できるというのが、ポテンシャルを導入した旨味なのです。ポテンシャル論をもっと勉強していくと、ラプラス方程式に代表される様々な微分方程式を、様々な境界条件下で解くといった問題が主なトピックになってきます。
.. _ベクトル場の流線と流管: http://www12.plala.or.jp/ksp/vectoranalysis/VectorFieldFlux/
.. _管状ベクトル場: http://www12.plala.or.jp/ksp/vectoranalysis/Solenoidal/
.. _層状ベクトル場: http://www12.plala.or.jp/ksp/vectoranalysis/PotentialIrrotational/
.. _ラプラス場: http://www12.plala.or.jp/ksp/vectoranalysis/LaplacianField/
@@author:Joh@@
@@accept: @@
@@accept: 2006-10-11@@
@@category: ベクトル解析@@
@@id: StreamFunction@@
記事ソース/流れ関数 のバックアップの現在との差分(No.1)
