学際大規模情報基盤共同利用・共同研究拠点 平成 25 年度共同研究 最終報告書 2014 年 5 月

13-NA30

機械工学分野におけるシミュレーション科学の新展開
滝沢 寛之（東北大学大学院情報科学研究科）
概要 本研究の目的は、機械工学分野の先進的・萌芽的な数値シミュレーションを対象とし、スーパーコンピューターを有効利用した場合の性能向上とさらなる高度化・高性能化に必要な課題を明確にすることである。本研究を遂行することにより、機械工学分野におけるスーパーコンピューターの利用促進と、その結果として数値シミュレーションの高度化・高性能化を実現することを目指してきた。研究期間前半で検討した結果、すべての対象シミュレーションコードに関して、計算コストの高いルーチンの特定や、スーパーコンピューターへの移植に向けた課題の明確化を行うことができた。研究期間後半には、大幅な速度向上が期待できるシミュレーションコードを２つ選択し、それらの高速化支援を成功裏に行った。並列化を強く意識することなく書かれた萌芽的なシミュレーションコードであったが、高いノード性能を提供する SX-9 の特長を活かすことにより、ベクトル化とスレッド並列化によって大幅な速度向上を達成することができた。

発電機器の安全・安心設計に資する熱流動シ

1. 研究の目的と意義
機械工学分野では様々な数値シミュレーショ

ミュレーション

ン技術が研究されているにもかかわらず、そのソ安全性の高い発電機器を設計することは、電力フトウェア開発技術やスーパーコンピューター利の安定供給のために必要不可欠である。スーパー用技術が共有されているとは言い難い。また、スコンピューターを利用することにより安全でかつーパーコンピューターの利用技術は近年ますます高効率な発電機器が開発できる。

高度化しており、機械工学分野の研究者がスーパーコンピューターの演算能力を引き出すために要熱硬化性樹脂の機械特性評価のための分子シする労力は年々増大している。

ミュレーション

本研究の目的は、現在上記の理由からスーパー分子シミュレーションを用いて樹脂の機械特性コンピューターの本格的な利用に至っていない機を評価することは、樹脂の分子構造および損傷・械工学分野の先進的・萌芽的な数値シミュレーシ破壊メカニズムに関する知見を得ることが出来、ョンを対象とし、スーパーコンピューターを有効それらが機械的特性に及ぼす影響を予測すること利用した場合の性能向上とさらなる高度化・高性ができる。また、マルチスケールで巨視構造の設能化に必要な課題を明確にすることである。

計を行う場合には、分子スケールから議論を行う機械工学分野の数値シミュレーションの高速化ことで、材料設計へのフィードバックを行うこと支援のためには、実用的なシミュレーション規模も可能になる。

での性能評価が必要不可欠であり、そのための大規模計算資源が必要不可欠である。特に萌芽的研

バイパス遷移過程の解明のための数値流体シ

究の段階では、並列化を考慮せずにプログラム開ミュレーション発されてきたコードが多いため、SX-9 のように比境界層は層流状態から乱流状態に遷移するが、較的大きな計算ノードで実行できるシステムが有乱流状態になることで摩擦抵抗がひと桁増大する用であると期待される。

ため、その遷移過程を解明し遷移を遅らせる方法

本研究で対象とする４種類の数値シミュレーシを見つけることができればエネルギー資源の大幅ョンは、それぞれ以下の理由から学術的・社会的な節減となる。

意義の大きい研究開発に必要なものである。
1

学際大規模情報基盤共同利用・共同研究拠点 平成 25 年度共同研究 最終報告書 2014 年 5 月

プラズマによるナノ粒子創製プロセスシミュ

ため、実験で捉えることは難しく、計算機シミュレーションレーションが必須となる。

プラズマによるナノ粒子創製プロセス全体を包共同研究者の茂田は、「プラズマによるナノ粒子括する数値解析システムを構築できれば、ナノ粒創製プロセスシミュレーション」を担当する。プ子形成機構やプラズマ流の複雑熱流動場との干渉

ラズマによるナノ粒子創製プロセス全体を包括す現象を解明できるのみならず、コンピュータ上でる数値解析システムの構築は、時間・空間ともに仮想的な実験を低コストで繰り返すことができる

マルチスケールに渡る複雑な物理現象の再現の試ようになるため、工学的な意義は極めて大きい。

みに等しく、莫大な計算容量と計算時間を必要とする。

2. 当拠点公募型共同研究として実施した意義(1) 共同研究を実施した拠点名および役割分担(2) 共同研究分野研究代表者の滝沢と東北大学サイバーサイエン

超大規模数値計算系応用分野

スセンターは、スーパーコンピューター利用技術やアプリケーション高速化支援に取り組んでおり、 (3) 当公募型共同研究ならではという事項など
特にベクトル型スーパーコンピューター向けの高スーパーコンピューター利用技術の高度化によ速化支援ではこれまでに多くの実績を残している。 り、機械工学分野におけるシミュレーション科学
このため、本研究においても、対象とするシミュのさらなる発展には、構成拠点に所属する計算機レーションの高速化支援を主に担当する。

科学分野の研究グループ、研究者の協力が必要不

共同研究者の山本、佐野、古澤、三宅は、「発電可欠である。本研究では、先進的あるいは萌芽的機器の安全・安心設計に資する熱流動シミュレーな手法に基づく数値シミュレーションを対象としション」の研究開発を担当する。発電機器におけており、本研究を遂行することによってスーパーる破壊や損傷に至る熱流動を把握するためには複コンピューターの本格利用につながることを大い雑な非定常熱流動を数値計算で解明する以外に有に期待できる。また、研究代表者および共同研究

効な手段がなく、スーパーコンピューターの高い者が所属する東北大学工学部機械系と同大学サイ演算性能が必要である

バーサイエンスセンターとは物理的に近く、頻繁共同研究者の岡部と青柳は、「熱硬化性樹脂の機に議論することができることから、密な連携を期械特性評価のための分子シミュレーション」を担待できる。

当する。分子シミュレーションにより力学特性を本研究を遂行することにより、機械工学分野に求めるためには、ある程度大きい系でシミュレーおけるスーパーコンピューターの利用促進とそのションを行う必要がある。この場合、系の規模が結果として数値シミュレーションの高度化・高性大きくなるほど計算の高速化が強く望まれる。ゆ能化が期待でき、機械工学の研究者および計算機えに、大規模系の分子シミュレーションにおいて科学の研究者の双方にとって有意義な共同研究とスーパーコンピューターによる数値シミュレーシ考えている。東北大学サイバーサイエンスセンタ

ョンが必要となる。

ーでは、長期にわたって継続的に利用者のアプリ共同研究者の福西と伊澤は、「バイパス遷移過程ケーションの高速化支援を行ってきた実績がある。の解明のための数値流体シミュレーション」を担

本研究においてもこれまで蓄積してきたノウハウ当する。遷移過程の解明のためには、境界層の開を十分に活用して萌芽的研究を支援し、機械工学始点近傍において入り込んだ微小な変動が下流に分野での成果につながることが期待できる。

向かって流下する過程で境界層に与える影響について詳細に調べる必要がある。この変動は微小な2

学際大規模情報基盤共同利用・共同研究拠点 平成 25 年度共同研究 最終報告書 2014 年 5 月

3. 研究成果の詳細と当初計画の達成状況条件の場合について数値解析した。これは 25K の(1) 研究成果の詳細について過冷却状態に当たる条件であり、入口では乾燥し本研究で対象とする数値シミュレーションは、た蒸気が流入しているため、タービン三段流路のそれ自身が有意義な成果を期待できる重要な数値適当な位置で非平衡凝縮が起こる。図２(a)、２(b)シミュレーションである。本研究で行う高速化支にそれぞれ、計算により得られた全流路ミッドス

援と並行して、それぞれの分野でも成果を残してパンならびに三段動翼列流路断面における瞬間湿いる。以下、本研究で対象とする 4 つの研究分野

り度分布を示す。まず、図２(a)より一段目静翼列のそれぞれにおける研究成果を報告する。

下流域で凝縮が開始していることがわかる。凝縮により生成された液滴は下流に行くに従って成長発電機器の安全・安心設計に資する熱流動シし、三段目動翼下流域では最大 4.8％の湿り度が得ミュレーションられた。また、湿り度は各翼列の後流の影響を受まずは数値タービンの並列処理効率を検証するけながら増加していることも示されている。図２

ため、蒸気タービン最終三段静動翼列の非定常湿(b)には、三段目動翼断面においては湿り度が極端り蒸気流れについて全周を仮定して計算した場合に偏在化している様子が示されている。湿り度はの計算加速率を算出した。図１に計算に用いた計

断面中央域が全体的に高く、エンドウォールに行算格子を示す。翼列流路ブロック数は 551 ブロッ

くにしたがって減少しているが、エンドウォールクである。

付近ではかなり非定常性が強い分布になっている。
これは、前段に位置する、一段目静翼列、一段目動翼列、二段目静翼列、二段目動翼列、そして三段目静翼列から発生したコーナー渦や翼後流のすべてが干渉しながら三段目動翼列まで流れ込んでくるため、これらが湿り度の分布に影響した結果であると考えられる。実は、湿り度 4.8％は翼端近傍で得られており、湿り度がコーナー渦や後流によりかなり偏在化していることを示唆している結果である。蒸気の湿りは蒸気タービンの効率やタ

図 1 全周計算に用いた計算格子ービン翼の壊食に直接影響を与える因子である。

東北大学サイバーサイエンスセンターSX-9 上でその絶対値の評価ならびに最大値の位置を把握す

MPI に基づく並列計算を実施した。ただし、計算ることは、高効率化ならびに高信頼化を目指す次反復回数は 10000 回で最終解は求めていない。そ世代蒸気タービン設計に有益な知見を与えることの結果、１CPU の場合には、計算時間が 386429ができる。

秒であった者に対して、SX-9 の最大 CPU 数であ

今後はこれらについてさらなる改良を加えて、並る 64CPU の場合には、8627 秒となり、その比で列処理効率の向上を目指す。一方で、計算対象とある速度向上率は、44.8 倍となった。

して、より実際のタービンを模擬するため、蒸気次に、周方向の流動が周期的であると仮定して、タービン最終多段長翼列を通る非定常湿り蒸気流翼列流路 13 ブロックに簡略化した計算格子に対れの数値解析を実施することを計画している。

して非定常解が求まるまで計算した。具体的には、三菱重工が設計した蒸気タービン低圧三段の静動熱硬化性樹脂の機械特性評価のための分子シ

翼列を通る非定常三次元湿り蒸気流れに対して、ミュレーション

入口全圧 0.3296MPa、入口全温 384.3K の入口境界

炭素繊維強化プラスチック(CFRP)は軽量かつ高3

学際大規模情報基盤共同利用・共同研究拠点 平成 25 年度共同研究 最終報告書 2014 年 5 月

強度という優れた機械特性を有することから航空算時間が大幅に短縮され、大規模系の構造で短時宇宙機器に使用されている。特に航空宇宙分野に間に多くのサンプルを得ることが出来た。大規模適用される CFRP においてはマトリクス樹脂にエ系により詳細に分子構造を議論できるだけでなく、計算時間の短縮化により多くの計算結果（サンプル）を取り機械特性の傾向を比較・考察することが可能となった。図 3 および図 4 に引張シミュレーション結果を示す。このような大規模系および計算時間の短縮化は分子シミュレーションを用いたマトリクス樹脂の開発（機械特性の傾向調査など）には不可欠であり、今後の樹脂開発に果たす役割は大きいと考えられる。

(a) ミッドスパン(b) 三段目動翼流路断面
図 3 引張シミュレーション前

図 2 瞬間湿り度分布
ポキシ樹脂が用いられている。CFRP の機械特性にはマトリクス樹脂の機械特性が影響するため、最近では分子シミュレーションを用いたマトリクス樹脂の開発が注目されている。本研究では、分子シミュレーションを用いて硬化反応シミュレーションおよび引張シミュレーションを行うことにより、エポキシ樹脂の機械特性を調査した。分子シミュレーションは、大規模系（分子数が多い）であるほど高精度であるが、その反面一般的なパソコンでは計算時間がかかることが課題であった。
そこで、本研究ではスーパーコンピューターSX-9

図 4 引張シミュレーション後を用いて、大規模系で分子シミュレーションを行

オープンソースコードの LAMMPS に加えてオいエポキシ樹脂の機械特性を調査している。
リジナルのエポキシ樹脂分子シミュレーションプ
オープンソースコード LAMMPS を用い、スーログラムを開発中であり、主剤、硬化剤およびそパーコンピューターSX-9 を使用しエポキシ樹脂れらの配合比が硬化特性に与える影響を分子シミの引張シミュレーションを行った。その結果、計ュレーションにより調査し、高分子材料分野の権4

学際大規模情報基盤共同利用・共同研究拠点 平成 25 年度共同研究 最終報告書 2014 年 5 月

威ある学術雑誌に掲載された。また、その高性能は解析対象とする渦とその他のスケールの渦を抽化に向けた性能分析とサイバーサイエンスセンタ出するのに使用したバンドバスフィルタの中心波ーの環境への移植方法を検討した。

数の比であり、縦軸は渦の引き伸ばされやすさを表す伸長速度係数である。図 5 から、kr/kt = 0.5 のバイパス遷移過程の解明のための数値流体シとき係数の値が最大となり、kr/kt > 1 では負の値ミュレーションを取ることがわかる。このことは、慣性小領域に本研究の目的は、境界層遷移及び乱流素過程の動ある渦は、乱流レイノルズ数に関わらず、自身の力学明らかにすることにある。乱流が示す強い混2 倍程度大きいスケールの渦が作る速度場によっ合・拡散能力や高い熱伝達性能といった特徴的なて引き伸ばされる傾向があることに加えて、自身性質は、乱流中に存在する大小様々な渦の不規則より小さい渦は引き伸ばしには寄与しないことをで乱雑な運動によって、単位時間当たりの運動量意味している。

や熱の輸送量が極めて高い状態に保たれているお境界層遷移のシミュレーションでは，まず速度蔭である。すなわち、乱流遷移を含む乱流現象を変動の rms 値 に着目して噴流噴射地点下流に生解明するためには、渦の動力学に対する深い理解成された乱流領域を判別し，そのうえで乱流が発が求められる。

生しはじめる領域近傍に存在する渦構造を観察し本年度は以下の内容で研究を進めた。まず前者のた．その結果，流れ方向に並んだ符号の異なる縦境界層遷移については、バイバス遷移が起こる際渦の相互干渉が乱流の初生に関係することが強くに現れるストリーク構造とこれらが発達して形成示唆された．しかし，これらの結果は周期境界をされる局所的な乱流領域（乱流斑点）の関係に着課した比較的狭い計算領域のものであることに加目し、人為的に発生させた低速のストリークに対え，乱流の判別方法もまだ不十分であることから，して壁面から噴流を吹くことで流れ場がどのようさらに詳細な検討が必要である．今後は、スーパに変化するのかについて調べ、乱流斑点が生成しーコンピューターを活用して計算領域及び格子解やすい“外部刺激”の条件について探った。

像度を十分に確保することを予定しており、本研

また、後者の乱流の研究については、スーパーコ究の枠組みの中でその大規模化の支援を行っていンピューターを利用して大規模乱流の直接数値計く予定である。また、乱流シミュレーションに関算（DNS）データの解析を進めた。乱流の統計理しても、現在さらに詳細な検討を進めている段階

論は、低波数のエネルギー保有領域と高波数のエである。

ネルギー散逸領域の影響が無視できる慣性小領域についてのみ成立する。このような乱流の普遍的な統計量を議論するためには少なくとも 1,0243 個以上の格子解像度が必要となり、パソコンレベルではもはや扱いきれない。乱流計算とその解析にはスーパーコンピューターの利用が欠かせない状況にある。本研究では、渦度ベクトルに直交する面内にフーリエバンドパスフィルタを適用することで乱流場を渦運動のスケールに応じて複数の流れ場に分離する方法を新たに提案し、それぞれの流れ場の中から渦構造の抽出を行って、これらの渦構造の伸張と流れ場のスケールとの関係につい図 5 渦のスケールと伸長速度て評価した。その結果の一例を図 5 に示す。横軸
5

学際大規模情報基盤共同利用・共同研究拠点 平成 25 年度共同研究 最終報告書 2014 年 5 月

および計算コードを開発していく予定である。

(2) 当初計画の達成状況について
本共同研究では、機械工学分野の先進的、萌芽的な数値シミュレーションプログラムをスーパーコンピューターの環境に移植し、それらのシミュレーションを高性能化、大規模化することで、機械工学分野での新しい研究成果の創出を支援することを目指してきた。課題申請書には、当初計画が以下のように記述されている。
--- 抜粋ここから --本申請では、スーパーコンピューターの利用によってさらなる発展を期待できる機械工学分野の数値シミュレーションを選定し、どの程度1000

Temperature [K]11000の高度化・高性能化を期待できるのかを定量的

図 6 下流側から見た熱流動場に評価する。平成 25 年度前半には研究室で利用（中央断面における温度分布、9000 K の等温面、流線）されている以下の数値シミュレーションをスー
パーコンピューターで実行してその性能評価をプラズマによるナノ粒子創製プロセスシミュ行い、研究室のパソコンレベルの計算能力からレーションスーパーコンピューターレベルに移行すること

現在、プラズマトーチ内の３次元流動現象が明らで期待できる研究成果と、それに要する開発コかになりつつある。図 6 に本計算によって得られストの見積を行う。平成 25 年度後半には、そのた熱流動場の代表的な瞬間像を示す。温度分布や中でも特に高性能化が期待できるプログラムを流線から、プラズマトーチ内は３次元的で複雑な開発コストも考慮して選別し、東北大学サイバ熱流動場となっており、しばしば仮定される軸対ーサイエンスセンターの SX-9 システム向けに称２次元の場とは程遠いものとなっていた。高温並列化・最適化を行って機械工学分野での成果のプラズマ領域には、らせん状の流線が示す再循に貢献することを計画している。

環流が見られることが多い。これは導電率の高い--- ここまで ---プラズマ領域に Lorentz 力が発生し、流体運動を当初計画のとおり、平成 25 年度前半には対象シ駆動するためである。またその再循環流はトーチミュレーションコードの現状把握と性能解析を主上部から注入される低温のガス流にも影響を与え目的とし、共同研究者の希望するすべてのシミュるため、結果としてトーチ内の熱流動場に片寄りレーションコードをスーパーコンピューターで実が生じ、直進性の高いプラズマジェットも蛇行す行してコスト分析を行った。同年度後半には、そる。さらに高温領域と低温領域の界面が波打っての分析結果に基づいて、従来よりも大幅な高速化いる。これは異なる密度の流体が互いに異なる速を期待できるコードとして、共同研究者 伊澤のコ

度 で 運 動 し て い て 、 界 面 に お い てードおよび茂田のコードを選別し、その高速化支Kelvin-Helmholtz 不安定性が発現するためである。 援を行った。
今後は、非定常３次元的に流動するプラズマ場に茂田のシミュレーションコードで最も高コストおけるナノ粒子群の集団成長および輸送現象の反のサブルーチンでは、SOR 法による反復計算が行応非平衡性も組み込んだ、より高度な数理モデルわれている。SX-9 へ移植して高性能を達成するた6

学際大規模情報基盤共同利用・共同研究拠点 平成 25 年度共同研究 最終報告書 2014 年 5 月

めに、そのサブルーチンに対してベクトル化を促

ドの大規模な修正を伴う MPI 化までを行うことは進するコード修正を行った。カーネルループの改できなかったが、課題申請書に書かれた当初の期変により、スカラ変数の初期化、ループ中の分岐待どおり、比較的大きな計算ノードで実行できるの除去、作業配列導入、指示行挿入による自動ベ

SX-9 システムの活用により大幅な速度向上率をクトル化によりベクトル化を推進し、平均ベクト達成することができた。並列化を強く意識せずにル化率を大幅に押し上げることで 7.8 倍の速度向書かれたコードをベクトル化とスレッド並列化の上を実現した。これは、プラズマ近傍の渦の相互みで高速化したことを考えると、十分に高い速度作用やナノ粒子群の空間輸送といった現実的な物向上率であると考えている。今後は MPI 化によっ理を工学的に有意な時間スケールを以ってシミュて更なる高速化と大規模化の実現を目指す予定でレートするための大きな進展と言える。

ある。

一方、伊澤のシミュレーションコードの実行時

また、共同研究者山本らのコードや岡部らのコ間の大部分は高速フーリエ変換(FFT)に費やされードに関してはコードの改変作業を行うまでにはており、研究室のパソコンレベルでの計算では至らなかったものの、スーパーコンピューター利Intel Math Kernel Library(MKL)が使われている。

用のための技術支援などを行い、それぞれの研究HPC Challenge Benchmark の Global FFT において、分野において平成 25 年度の成果を出すための支地球シミュレータ 2 が好成績を残してきた例から援を行った。その支援を通じて、スーパーコンピもわかる通り、ベクトル型スーパーコンピュータ

ューター利用者と高速化支援者との間の連携の重ーには FFT の実行において高い実効性能を期待で要性を再確認するとともに、その価値観の違いにきる。このため、本共同研究では、伊澤のコード起因する連携の難しさも再認識できた。その結果の SX-9 向けチューニングを行った。FFT ライブラ

として得られた知見は、今後の計算機工学な研究

リを MKL から NEC 社製 ASL に変更するところからの方向性を検討する上で重要な示唆を与えるもの

移植作業を始め、ベクトル化やスレッド並列化をであると考えている。

進めることにより、従来の実行環境で達成していた性能の 76 倍の速度向上を達成することができ

5. 研究成果リストた。この飛躍的な速度向上の意義は、単なるプロ(1) 学術論文

グラム実行の高速化にとどまらず、乱流遷移を含

1.

む乱流現象の解明に関する研究の進捗を高速化す

Okabe, Ian A. Kinloch, Robert J. Young, "Investigationるという観点からも有意義であると考えている。

of the sp3 structure of carbon fibres using UV-Raman

本研究の場合、達成状況を定量的に評価すること

spectroscopy", TANSO, 259 (2013), 1-5.

は難しいが、対象アプリケーションをすべて解析

2. Tomonaga Okabe, Toshiki Sasayama, Jun Koyanagi,

し、選択された 2 つのシミュレーションについて"Micromechanical Simulation of Tensile Failure ofは大幅な性能向上を達成している。当初計画に記Discontinuous

載された予定通りに円滑に進捗し、速度向上率でComposites using Spring Element Model", Composites

は当初の期待以上の数値を達成できたことから、

Part A, 56 (2014), 64-71.

達成状況は良好であると自己評価している。

3. Jun Koyanagi, Yukihiro Sato, Tomonaga Okabe,Fumihiko Tanaka, Daniel Wolverson, Tomonaga

Fiber-Reinforced

PolymerMatrix

Satoru Yoneyama, "Numerical simulation of strain-rate
dependent transition of transverse tensile failure mode

4. 今後の展望

in fiber-reinforced composites", Composites Part A, 56

本共同研究では、当初計画の通り萌芽的シミュレーションコードの解析と高速化支援を成功裏(2014) 136–142.に行うことができた。この短い期間内では各コー

4. Fumihiko Tanaka, Tomonaga Okabe, Haruki Okuda
7

学際大規模情報基盤共同利用・共同研究拠点 平成 25 年度共同研究 最終報告書 2014 年 5 月

Ian A. Kinloch Robert J. Young, "Factors controlling

4. Xinrong Su, Satoru Yamamoto and Xin Yuan, On

the strength of carbon fibers in tension", Composites

the Accurate Prediction of Tip Vortex: Effect ofPart A, 57 (2014) 88–94.

Numerical Schemes, Proc. ASME Turbo Expo

5.

Y.Sato,T.Okabe,R.Higuchi,

K.Yoshioka,

2013, San Antonio, GT2013-94660, (2013),"Multiscale approach to predict crack initiation inunidirectional

off-axis

laminates",

CD-ROM(全文査読).

Advanced

5. Satoshi Miyake, Itsuro Koda, Satoru Yamamoto,

Composite Materials, accepted for the publication.

Yasuhiro Sasao, Kazuhiro Momma, Toshihiro

6. J. Watanabe, F. Tanaka, H. Okuda and T. Okabe,

Miyawaki and Hiroharu Ooyama, Unsteady Wake

"Tensile strength distribution of carbon fibers at shortand Vortex Interactions in 3-D Steam Turbinegauge

Low Pressure Final Three Stages, ASME Turbolength",

Advanced

Composite

Materials,

accepted for the publication.

Expo 2014, GT2014-25491, (2014), (掲載予定).

7. R. Higuchi, T. Okabe, Y. Ohtake, T. Honda, Y. Ueda

6. Masaya Shigeta, 3D fluid dynamic modelling ofand S. Ogihara, "Numerical simulation of damage

RF plasmas,（招待講演）, Proceedings of the 31stprogression and fracture in structures made of 3DInternational Conference on Phenomena in Ionized

woven

Gases, Granada, Spain, (July 14-19, 2013), USB

ceramic

matrixcomposites",

Advanced

Composite Materials, accepted for the publication.

Memory, 2 pages.
7. Masaya

Shigeta,

Numerical

Simulationof(2) 国際会議プロシーディングス

Multiple Thermal Plasma Plumes from a Jet Array,

1. Satoshi Miyake, Satoru Yamamoto, Yasuhiro

The 12th International Symposium on Fluid

Sasao, Kazuhiro Momma, Toshihiro Miyawaki

Control, Measurement and Visualization, Nara,and Hiroharu Ooyama, Unsteady Flow Effect onJapan, (November 18-23, 2013), USB Memory, 10

Nonequilibrium

Lowpages.

Pressure Steam Turbine Stages, Proc. ASME

8. Masaya

Condensation

in3-D

Shigeta,Modelsofnano-powder

Turbo Expo 2013, San Antonio, GT2013-94832,

formation in thermal plasma synthesis, Book of(2013), CD-ROM(全文査読).

Abstracts of 4th International Round Table on

2. Yasuhiro Sasao, Satoshi Miyake, Kenji Okazaki,
SatoruYamamotoandHiroharu

Thermal Plasmas for Industrial Applications:

Ooyama,Challenges

Eulerian-Lagragian Numerical Simulation of WetandOpportunities,

Marrakech,

Morocco, (March 3-7, 2014), p. 8.（招待講演）

Steam Flow through Multi-stage Steam Turbine,
Proc. ASME Turbo Expo 2013, San Antonio,(3) 国際会議発表

GT2013-95945, (2013), CD-ROM(全文査読).
3. Tadashi

Tanuma,Yasuhiro

Sasao,Satoru(4) 国内会議発表

Yamamoto, Yoshiaki Niizeki, Naoki Shibukawa,1.

山本悟、日本機械学会 2013 年度年次大会、ワ

Hiroshi Saeki, Aerodynamic Interaction Effects

ークショップ「CFD の産業活用における方向fromthe性」、マルチフィジックス CFD の産業応用、

Down-flow Type Exhaust Diffuser Performance in2013 年度年次大会 DVD 論文集、(2013-9)（招A Low Pressure Steam Turbine, Proc. ASME

待講演）.

UpstreamandDownstreamonTurbo Expo 2013, San Antonio, GT2013-95901,

2.

(2013), CD-ROM(全文査読).

笹尾泰洋、三宅哲、山本悟、蒸気タービン内部における液滴挙動のシミュレーション、日本機械学会 2013 年度年次大会 DVD 論文集、8

学際大規模情報基盤共同利用・共同研究拠点 平成 25 年度共同研究 最終報告書 2014 年 5 月(2013-9).
3.

三宅哲、笹尾泰洋、山本悟、蒸気タービン三段翼列を通る非平衡凝縮流れの三次元非定常解析、日本機械学会 2013 年度年次大会 DVD
論文集、(2013-9).

4.

岩崎俊樹、三宅哲、笹尾泰洋、山本悟、スタッキング翼からなるタービン多段静動翼列の非定常三次元流動解析、第９１期日本機械学会流体工学部門講演会 DVD 論文集、(2013-11).

5.

三宅哲、笹尾泰洋、山本悟、非平衡凝縮を伴う蒸気タービン三段静動翼列の三次元非定常流動解析,第２７回数値流体力学シンポジウ
ム, (2013-11), CD-ROM.

6.

吉川穣、茂田正哉、伊澤精一郎、福西祐, 「ストリーク構造と噴流の干渉による境界層の不安定化の促進」, 第 91 期流体工学部門講演会
講演論文集, (2013-11).

7.

廣田真人、茂田正哉、伊澤精一郎、福西祐, 「乱流の階層的渦構造の抽出とその運動に関する研究」, 第 27 回数値流体力学シンポジウム講演論文集, (2013-12).

8.

茂田正哉、
「熱プラズマを用いたナノ粒子群生成の諸現象に関する数理モデル」, 第 23 回日
本 MRS 年次大会(日本 MRS 学術シンポジウム)アブストラクト集, (2013-12), P-I10-008,CD-ROM.（招待講演）(5) その他（特許、プレス発表、著書等）
1.

山本悟, ガスタービン工学、日本ガスタービ
ン学会、(2013-8), 83-91(著書、分担執筆).

9