射影法

射影法（しゃえいほう、）とは、非圧縮性流れの問題を時間発展的に解くときの手法のひとつである。投影法（とうえいほう）、プロジェクションなどとも言う。
ここで”射影”とは、回転成分および発散成分の双方を含むベクトル場から発散を消し回転成分を取り出す操作のこと。
1967年にによって、非圧縮性ナビエ・ストークス方程式を解くために導入された
。この方法は、圧力勾配項とそれ以外の項を分離して計算するため、分離型解法とか演算子分割とか多段法などと呼ばれる。
１時刻の計算は２段階に分けられる。第一段階では、圧力以外の効果による速度増分を計算し、仮の速度場（速度予測子などと呼ばれる）を得る。つづく第二段階で、速度の発散成分を消す。すなわちdivergence-free の空間へ射影する。実際の射影プロセスは反復解法によりポアソン式を解くことであり、同時に圧力場が得られる。
ヘルムホルツ・ホッジ分解.
任意のベクトル場は回転成分（）と流入出成分（に分解できる。（ヘルムホルツ・ホッジ分解と呼ばれる。）
つまり、
と分解できる。ここで formula_1 は回転成分であり formula_2 である。また formula_3 はスカラー場である（任意の流入出成分は対応するスカラー場の勾配で表される）。
よって両辺の発散をとると
という、formula_3に関するポアソン方程式の問題に帰着する。
任意ベクトル場 formula_5 が与えられたとき、上記方程式の formula_3 は一意に定められる。
すなわち、発散ゼロ場が
として得られる。
このようにポアソン式を解くことで圧力場と発散ゼロの速度場を得る。
ベクトルの直交成分の一つを消す操作であることから射影法と呼ばれる。
コリンの射影法.
非圧縮性すなわち密度一定のナビエストークス方程式（流体の運動量保存則）は次の通り。
コリンの原論文では、この運動量式に仮の流速 formula_7 を導入し、圧力勾配項と他の項を分ける。
 {\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} \\
つまり
となるような仮の流速 formula_7 を考える。ここで上付きの formula_9 は、時刻を表す。
すると、まず第一段階として、圧力勾配項を無視した(1)式により
から仮の速度場 formula_7 が得られる。
この仮速度と本来の速度 formula_11 との関係が(2)式である。
この第一段階の粘性と第二段階の圧力項を順次適用するやり方は、などと呼ばれる。
この第二段階で時刻 n+1 での流速を計算するためには、右辺第2項で時刻n+1での圧力を必要とする。ここで、この式の両辺の発散を取れば、左辺はformula_12となり、よって
というポアソン方程式を反復計算でとくことで、時刻 n+1 の圧力 formula_13 を得る。
フラクショナルステップ.
射影法の計算を実現するには速度と圧力を分ける必要がある。また、移流項、粘性項、外力項もそれぞれに分けて多段処理とした方がポアソン式の反復計算の収束はより確実であり、計算量も大きくは増えない。よって実際の計算は例えば以下のように多段分割される。
このような多段分割は、フラクショナルステップ、 Fractional step projection などとも呼ばれる。
日本国内ではMAC法およびSMAC法との対比から、N番目ステップの速度から圧力効果無視の速度予測子を得てN+1番ステップの圧力で補正しN+1ステップの非圧縮速度場を得る、すなわち圧力の時間差分のみ陰的のスキームを指して特にフラクショナルステップ法と呼ぶ。