1998 Fiscal Year Annual Research Report
Project/Area Number |
09680474
|
Research Institution | TOKYO INSTITUTE OF TECHNOLOGY |
Principal Investigator |
矢部 孝 東京工業大学, 大学院・総合理工学研究科, 教授 (60016665)
|
Co-Investigator(Kenkyū-buntansha) |
青木 尊之 東京工業大学, 原子炉工学研究所, 助教授 (00184036)
|
Keywords | メソスケール / ボルツマン / ブラソク / 数値解法 / CIP / Euler法 / 粒子法 / スプライン法 |
Research Abstract |
ミクロな現象を見る計算手法としては分子動力学法、マクロの現象を見るためには流体計算手法が種々提案されているが、本研究ではこの中間領域のメソスケールの現象を追跡する計算法を追求した。この目的のために、計算グリッドが少なくても高精度の計算が可能なCIP法を用いたここで用いたボルツマン方程式は空間及び速度空間の6次元の方程式であり、時間発展はこの6次元位相空間での移流を解くことになる。このため、移流方程式に対して非常に有力なCIP法を採用した。ここでは、CIP法の1次元手法を6回繰り返して用いる簡便な手法を採用して超多次元空間の計算法を開発した。この手法をチェックするために、まず、位相空間2次元と4次元のブラソフ方程式を解き、線形ランダウ減衰や二流体不安定性に適用した。 従来用いられている粒子法およびスプライン法とも比較をしたが、粒子法ではCIP法の10倍の粒子数を用いなければ同等の結果を得られないことがわかった。また、たとえ10倍の粒子数を用いたとしても、メソスケールに適用できるほど十分な分布関数の形状が得られないことが分かった。またスプライン法では、分布が負になったりスパイク状の不安定性が現れやすい。この比較により、CIP法の有用性を示すことができた。 速度空間グリッドを変化させ、どのくらい少ないグリッドまで可能かをチェックした結果、速度空間グリッドが10以下であっても充分な精度を得ることができることがわかった。最終的に6次元のコードを作成し、実計算を行ったが、Alpha-CPUの1台のパソコンでも計算可能であり、16×16×16×8×8×8の6次元グリッドでのランダウ減衰のシミュレーションが約7時間で終了した。
|
Research Products
(8 results)
-
[Publications] 田島雅和,渡辺和義,矢部孝: "CIP法による小型船体のスラミング現象のシミュレーション" 日本機械学会論文集(B編). 64巻618号. 63-70 (1998)
-
[Publications] 矢部 孝: "固体・液体・気体の汎用シミュレーシュン手法(解説)" シミュレーション. 17巻・3号. 7-13 (1998)
-
[Publications] 矢部 孝: "パルスレーザー加工・溶接現象のシミュレーション(解説)" レーザー研究. 26巻11号. 788-792 (1998)
-
[Publications] Y.Ogata,T.Yabe: "Shock Capturing with Improved Numerical Viscosity In Primitive Euler Representation" Comput.Phys.Commun.(印刷中).
-
[Publications] S.Yoon,T.Yabe: "The Unified Simulation for Incompressible and Compressible Flow by Predicter-Corrector Scheme based on CIP Method" Conput.Phys.Commun.(印刷中).
-
[Publications] T.Yabe,Y.Zhang,F.Xiao: "A Numerical Procedure-CIP-to Solve All Phases of Matter Together" Proc.16th Int.Conf.Num.Meth.Fluid Dynamics. 439-457 (1998)
-
[Publications] T.Yabe: "Mathematical Modeling of Weld Phenomena (分担)" The University Press,Cambridge, 683 (1998)
-
[Publications] T.Yabe: "Computational Fluid Dynamics Review 1998 (分担)" World Scientific, 1130 (1998)