第章非圧縮性流れの解法C-YamamotoSanoLaboratory.pdf

MA C 第 章 非圧縮性流れの解法 型解法 e function metho d ３次元流れや乱流の計算には，流速と圧力の方程式を解く原始関数法 primitiv が適し vierStok esNS ている．この種の解法の基本は Na 方程式と連続方程式を解くものである．しかしながら， 実際にこれらの方程式を解こうとすると，NS 方程式から流速を求めることはできるが，連続方程式から残 りの変数，圧力を求めることができない．このことは，圧縮性流れの解法で連続方程式から密度，運動方程 式から流速，エネルギー方程式から例えば岐点内部エネルギーをそれぞれ求めることができるのとは対照 的である．このために非圧縮性流れの解法では何らかの戦略が必要になる．この章には最も一般的に用い C MA C NS られてきた MA 型の解法について述べる． 法では 方程式と，連続方程式に代わるものとして圧 力方程式が解かれる．他に 非圧縮性流れの有力な解法としては，連続方程式に擬似圧縮項を加えた方程式 から圧力を計算する擬似圧縮性法がある．ところで NS 方程式と圧力方程式を解く方法では，解が連続の条 C 件を満足しなくなる虞がある．MA 法ではこの難点を，いわゆる食違い格子を用いまた圧力の差分方程式 にある付加項を導入することによってみごとに克服している．なおこのアイディアは，曲線座標格子に拡 張することはできるが，任意形状要素を用いる有限要素法に拡張することは難しい．そのために有限要素 法では，NS 方程式と連続方程式をガレルキン法で同時に解く方法，ペナルティ関数法，ラグランジュ剰数 法，ソレノイダルの重み関数を用いる方法など各種の解法が考案されてきた． MA C 法 wW elc h MA Cmark er and cell Harlo によって提案された 法は ，はじめはラグランジュ座標系における マーカー粒子の動きを求めると同時にオイラー座標系における格子セルの方程式を解くものであった．そ の後マーカー粒子は自由表面にのみ置かれるようになり，通常は格子セルの方程式のみが解かれるように C なった．ここには格子セルの方程式の解き方について述べる．MA 法の基礎方程式は次の保存形で書かれ vierStok es た Na 方程式と圧力方程式である． u u d u u u u u u u g r u u u u r p r u u g g