境界适合格子-名古屋大学.PPT

下　図 下　図 モンテカルロ衝突のアルゴリズム （「空衝突」導入） ? 全ての粒子についての計算は　 3. のみに ? 4., 5. は N×Pnull 個の粒子に 　 ついて計算 境界適合格子 ? 写像の概念を用いて、等間隔直交格子上（計算空間）で計算 を行う （詳細は後述） ? 格子にそこそこの柔軟性はあるが、特定の領域に格子点を 極端に集中させるのは困難な場合が多い マルチブロック、重なり格子 複雑な形状の場合には、複数の格子を組み合わせて 流れ場を覆う方法が一般的に用いられる マルチブロック法 重なり格子法 境界適合格子における変換式 (Appendix 1) という形の方程式は、計算空間 (x, h) では、 となる。 ここで 　　　　　　　　　　　　　　　である。 　 例えば、点 (i, j) での の差分式は などで数値的に評価しておく。 粒子コードに境界適合格子を適用するためには、 　 　　多数個の粒子量（ラグランジュ量）　と 　　非一様な格子点上の場の量（オイラー量） 　 を効率的に結びつける方法が必要である。 　 例えば、粒子量（ 　　　　）からソース項（　　　）を 格子点上でいったん求めれば、場の量に関する 方程式については Appendix 1 の変換式がその まま使用できる。 粒子コードへの適用 粒子コードへの適用（２） 境界適合格子　＋　粒子コード　の具体的な計算手順は、 　 １．粒子量（　　　 ）を格子点上に集めてソース項（　　 ）を計算する ２．格子点上で場の量（　　　 ）を進める　（計算空間、 Appendix 1） ３．格子点上の場の量（　　　 ）から、粒子に働く力　　　を計算する ４．粒子量（　　　 ）を進める　（物理空間） 下線部分の粒子量と格子点上の　　　　　　　　 「補間法」 場の量とのやりとりを行う方法　　　　　　　　　　「粒子局在化法」 通常の等間隔直交格子上での area-weighting method では、 粒子に働く場は格子点上の場 より、 で与えられる。 補間法（１）：　補間係数の求め方 補間法（２）：　補間係数の求め方 粒子局在化法 「粒子追跡アルゴリズム」 金星電磁圏シミュレーション 金星電磁圏のグローバルハイブリッド シミュレーション [Terada et al., 2002] 金星電離圏 * 非一様空間格子幅?時間ステップ幅 および衝突効果 寺田 直樹 名古屋大学太陽地球環境研究所 １．Introduction ?　しかし、粒子?ハイブリッドモデルでは 　　 直交格子を用いるのが一般的 非一様な格子分布の利用は、数値流体力学の分野では ２０年以上にわたり盛んに研究が行われている。 例２．衝撃波の反射解析 　 解像が必要な領域に計算格子を集中 　 複雑な形状の物体周りの場を解く （適切な境界条件） 例１．飛翔体周りの電磁環境の解析 講義の流れ 　? 非一様な格子分布　（テキスト p.233～） 　 を粒子?ハイブリッドモデルで利用する為の 手法を中心に説明。 　 　? 粒子数の動的変更　（テキスト p.237） 　? 非一様な Δｔ　（テキスト p.238～） 　? 衝突効果　（テキスト p.240～） 非一様な格子分布と併用されることがしばしばある ２．非一様な格子分布 （流体計算） ? 構造格子 (b) ? 非構造格子 (a) ? 構造格子と非構造格子の組み合わせ ? 解適合格子 格子の種類 非構造格子 ? 格子点の並びに規則性をもたない 　 ? 格子の自由度が高い 　 ? 有限要素法や有限体積法で使用される 　 ? 三角形（三次元では四面体）が最も多く用い られるが、四角形（三次元では六面体）や その他のセル形状やその混在もしばしば 用いられる 例．非構造格子を用いた デルタ翼周りの流れ場の計算 構造格子 ? 格子点の並びに規則性をもつ ? 差分法や有限体積法で使用される などがある 境界適合格子 マルチブロック 重なり格子 格子自由度の 低さを補う 例．境界適合格子を用いた 電極間プラズマの計算 非構造格子と構造格子の組み合わせ 解適合格子 流れの計算結果をもとに格子点を移動、 または格子点を追加?削除 三角形（三次元では四面体）セルを用いると 流速計算が多くなる ?　高レイノルズ数流れの計算では、高い 　　 空間分解能が必要となる物体近傍に 　　 のみ構造格子（またはプリズム形状の 　　 非構造格子）を用いる　→ 例 非構造格子 [Togashi et al., 2001] [Togashi et al