射影法のソースを表示

{{otheruses|数値流体力学の解法|その他|投影法}}
'''射影法'''（しゃえいほう、{{lang-en|projection method}}）とは、[[非圧縮性流れ]]の問題を時間発展的に解くときの手法のひとつである。'''投影法'''（とうえいほう）、プロジェクションなどとも言う。

ここで”射影”とは、回転成分および発散成分の双方を含むベクトル場から発散を消し回転成分を取り出す操作のこと。

[[1967年]]に{{日本語版にない記事リンク|アレクサンダー・コリン|en|Alexandre Chorin}}によって、非圧縮性[[ナビエ・ストークス方程式]]を解くために導入された
<ref>
{{Citation 
 | surname1 = Chorin
 | given1 = A. J.
 | title = The numerical solution of the Navier-Stokes equations for an incompressible fluid
 | journal = Bull. Am. Math. Soc. 
 | volume = 73
 | year = 1967
 | pages = 928&ndash;931
 | url =http://math.berkeley.edu/~chorin/chorin67.pdf
}}</ref> 
<ref>
{{Citation 
 | surname1 = Chorin
 | given1 = A. J.
 | title = Numerical Solution of the Navier-Stokes Equations
 | journal = Math. Comp.
 | volume = 22
 | year = 1968
 | pages = 745&ndash;762
 | url =
 | doi=10.1090/s0025-5718-1968-0242392-2
}}</ref>。この方法は、圧力勾配項とそれ以外の項を分離して計算するため、分離型解法とか演算子分割とか多段法などと呼ばれる。

１時刻の計算は２段階に分けられる。第一段階では、圧力以外の効果による速度増分を計算し、仮の速度場（速度予測子などと呼ばれる）を得る。つづく第二段階で、速度の発散成分を消す。すなわちdivergence-free の空間へ射影する。実際の射影プロセスは反復解法によりポアソン式を解くことであり、同時に圧力場が得られる。

== ヘルムホルツ・ホッジ分解 ==
{{Main|ヘルムホルツ・ホッジ分解}}

任意のベクトル場は回転成分（{{日本語版にない記事リンク|管状回転ベクトル場|en|Solenoidal vector field}}）と流入出成分（{{日本語版にない記事リンク|渦なし場|en|irrotational field}}に分解できる。（[[ヘルムホルツ・ホッジ分解]]と呼ばれる<ref>{{cite book | title = A Mathematical Introduction to Fluid Mechanics | author1 = Chorin, A. J. | author2 = J. E. Marsden | edition = 3rd | publisher = [[Springer Science+Business Media|Springer-Verlag]] | year = 1993 | isbn = 0-387-97918-2}}</ref>。）

つまり、
{{indent|<math>
  \mathbf{u} = \mathbf{u_s} + \mathbf{u_r} = \mathbf{u_s} + \nabla \phi
</math>}}
と分解できる。ここで <math>\mathbf{u_s}</math> は回転成分であり <math>\nabla \cdot \mathbf{u_s} = 0</math> である。また <math>\phi</math> はスカラー場である（任意の流入出成分は対応するスカラー場の勾配で表される）。

よって両辺の発散をとると
{{indent|<math>
  \nabla \cdot \mathbf{u} = \nabla \cdot \mathbf{u_s} + \nabla \cdot \nabla \phi = \nabla^2 \phi
</math>}}
という、<math>\phi</math>に関する[[ポアソン方程式]]の問題に帰着する。


任意ベクトル場 <math>\mathbf{u}</math> が与えられたとき、上記方程式の <math>\phi</math> は一意に定められる。

すなわち、発散ゼロ場が
{{indent|<math>
  \mathbf{u_s} = \mathbf{u} - \nabla \phi
</math>}}
として得られる。


このようにポアソン式を解くことで圧力場と発散ゼロの速度場を得る。

ベクトルの直交成分の一つを消す操作であることから射影法と呼ばれる。

== コリンの射影法 ==
非圧縮性すなわち密度一定のナビエストークス方程式（流体の運動量保存則）は次の通り。
{{indent|<math>
  \frac {\partial \mathbf{u}} {\partial t} + (\mathbf{u}\cdot\nabla)\mathbf{u} = - \frac {1}{\rho} \nabla p + \nu \nabla^2 \mathbf{u}
</math>}}

コリンの原論文では、この運動量式に仮の流速 <math>\mathbf{u}^*</math> を導入し、圧力勾配項と他の項を分ける。
{{indent|<math>
\begin{align}
  \frac {\partial \mathbf{u}} {\partial t} &= \frac {\mathbf{u}^{n+1} - \mathbf{u}^n} {\Delta t} = \frac {\mathbf{u}^{n+1} - \mathbf{u}^*} {\Delta t} + \frac {\mathbf{u}^* - \mathbf{u}^n} {\Delta t} \\
  &= -(\mathbf{u}\cdot\nabla)\mathbf{u} - \frac {1}{\rho} \nabla p + \nu \nabla^2 \mathbf{u}
\end{align}
</math>}}
つまり
{{indent|
<math>  
\begin{align}
  (1) \qquad& \frac {\mathbf{u}^* - \mathbf{u}^n} {\Delta t} = -(\mathbf{u}\cdot\nabla)\mathbf{u} + \nu \nabla^2 \mathbf{u} \\
  (2) \qquad& \frac{\mathbf{u}^{n+1} - \mathbf{u}^*} {\Delta t} = - \frac {1}{\rho} \nabla p
\end{align}
</math>}}
となるような仮の流速 <math>\mathbf{u}^*</math> を考える。ここで上付きの <math>n</math> は、時刻を表す。

すると、まず第一段階として、圧力勾配項を無視した(1)式により
{{indent|<math>
  \mathbf{u}^* = \mathbf{u}^n + \Delta t \left\{ -( \mathbf{u} \cdot \nabla ) \mathbf{u} + \nu \nabla^2 \mathbf{u} \right\} 
</math>}}
から仮の速度場 <math>\mathbf{u}^*</math> が得られる。

この仮速度と本来の速度 <math>\mathbf{u}^{n+1}</math> との関係が(2)式である。
{{indent|<math>
  \mathbf{u}^{n+1} = \mathbf{u}^* - \frac {\Delta t}{\rho} \nabla p ^{n+1}
</math>}}

この第一段階の粘性と第二段階の圧力項を順次適用するやり方は、{{日本語版にない記事リンク|演算子分割法|en|operator splitting}}などと呼ばれる。

この第二段階で時刻 n+1 での流速を計算するためには、右辺第2項で時刻n+1での圧力を必要とする。ここで、この式の両辺の発散を取れば、左辺は<math>\nabla\cdot \mathbf{u}^{n+1} = 0</math>となり、よって
{{indent|<math>
  \nabla ^2 p^{n+1} = \frac {\rho} {\Delta t} \nabla\cdot \mathbf{u}^*
</math>}}
というポアソン方程式を反復計算でとくことで、時刻 n+1 の圧力 <math>\nabla ^2 p^{n+1}</math> を得る。

== フラクショナルステップ ==

射影法の計算を実現するには速度と圧力を分ける必要がある。また、移流項、粘性項、外力項もそれぞれに分けて多段処理とした方がポアソン式の反復計算の収束はより確実であり、計算量も大きくは増えない。よって実際の計算は例えば以下のように多段分割される。

# まず、任意の時刻の即度場は、射影後の非圧縮場かそれに準じた分布とする。
# 移流項を計算する。多段法の場合、移流項は射影直後の非圧縮に極力近いベクトル場に適用するのが望ましい。
# この段階を'''予測子'''（よそくし）段階と呼ぶ。
# スキームによってはこの中間段階の流速に対しても射影操作する。
# 粘性項等や外力項を計算する。
# この段階はときに'''修正子'''（しゅうせいし）段階と呼ばれる。（とくに陰的な圧力項反映の後）
# 最終の射影操作を実施。これをもって次の時刻の速度場とし、1ステップの計算が完了する。
このような多段分割は、フラクショナルステップ、 Fractional step projection などとも呼ばれる。

日本国内では[[MAC法]]およびSMAC法との対比から、N番目ステップの速度から圧力効果無視の速度予測子を得てN+1番ステップの圧力で補正しN+1ステップの非圧縮速度場を得る、すなわち圧力の時間差分のみ陰的のスキームを指して特にフラクショナルステップ法と呼ぶ。

== 参考文献 ==
{{reflist}}

{{デフォルトソート:しやえいほう}}
[[Category:数値流体力学]]
[[Category:数値解析]]