境界适合格子.ppt

下　図 下　図 モンテカルロ衝突のアルゴリズム （「空衝突」導入） ? 全ての粒子についての計算は　 3. のみに ? 4., 5. は N×Pnull 個の粒子に 　 ついて計算 境界適合格子 ? 写像の概念を用いて、等間隔直交格子上（計算空間）で計算 を行う （詳細は後述） ? 格子にそこそこの柔軟性はあるが、特定の領域に格子点を 極端に集中させるのは困難な場合が多い マルチブロック、重なり格子 複雑な形状の場合には、複数の格子を組み合わせて 流れ場を覆う方法が一般的に用いられる マルチブロック法 重なり格子法 境界適合格子における変換式 (Appendix 1) という形の方程式は、計算空間 (x, h) では、 となる。 ここで 　　　　　　　　　　　　　　　である。 　 例えば、点 (i, j) での の差分式は などで数値的に評価しておく。 粒子コードに境界適合格子を適用するためには、 　 　　多数個の粒子量（ラグランジュ量）　と 　　非一様な格子点上の場の量（オイラー量） 　 を効率的に結びつける方法が必要である。 　 例えば、粒子量（ 　　　　）からソース項（　　　）を 格子点上でいったん求めれば、場の量に関する 方程式については Appendix 1 の変換式がその まま使用できる。 粒子コードへの適用 粒子コ