Drastic reduction of computing time in dynamic simulation for hydrogen isotope distillation columns

木下 正弘

Fusion Technology, 9, p.492 - 498, 1986/00

水素精留塔の動的シミュレーションにおいては、基本となる常微分方程式が非常に大きな次元及び「堅さ」を有するため、長大な計算時間が大きな問題となることが多い。著者は、現有の積分アルゴリズムをレビューし、そのうちのいくつかをテストすることにより、安定性及び簡便さの点からBallard-Brosilowのアルゴリズムを選定した。このアルゴリズムを用いると、各ステップにおいて、線形の三項方程式を解くこととスカラー沸点計算が必要であるに過ぎない。いくつかの数値実験をもとに、時間きざみの大きさをシミュレーションの過程で修正してゆく手順を提案した。1つの定常状態から他の定常状態に移る過程を対象とする典型的な例においては、例えば改良オイラー法に比べると約百分の一の計算時間において十分な精度で数値解が得られた。

改訂版ERATO-JによるMHD安定性解析

徳田 伸二; 常松 俊秀; 安積 正史; 滝塚 知典; 奈良岡 賢逸*; 竹田 辰興

JAERI-M 9899, 56 Pages, 1982/01

JAERI-M-9899.pdf:1.34MB

線型理想MHD安定性解析コード、ERATO-J、の4点の改良を行なった。(i)高精度写像モジュールを開発。この写像法により固有値の収束性が改良された。(ii)中間および高トロイダルモード安定性解析に不可欠な準モード表示法の追加。(iii)真空磁場計算へのベクトル・ポテンシャル法の採用。(iv)大次元疎行列の圧縮。以上の改良により精度及び計算効率が飛躍的に向上した。この改良版ERATO-Jを用いて、n=1内部キンク・モードおよびバルーニング・モードの安定性を行った。N=1内部キンク・モードに対する第2安定領域が見つかった。バルーニング・モードの成長率へのトロイダルモード数およびシアの影響が明らかになった。