JAXA スーパーコンピュータシステム利用成果報告(2022 年 2 月～2023 年 1 月)

将来輸送技術の研究(エンジン流路形状の研究)

報告書番号：R22JG3104
利用分野：研究開発
URL：https://www.jss.jaxa.jp/ar/j2022/20698/
責任者
沖田耕一, 研究開発部門第四研究ユニット
問い合せ先
宇宙航空研究開発機構研究開発部門第四研究ユニット 佐藤 茂
sato.shigeru@jaxa.jp(sato.shigeru@jaxa.jp)
メンバ
佐藤 茂, 福井 正明, 井上 拓, 宗像 利彦, 高橋 正晴

事業概要
再使用型宇宙推進機関であるスクラムジェットエンジンに就いて, その内部形状がエンジン性能に及ぼす影響を CFD 援用にて空気力学的に解明し, 同エンジンの設計に資すること.
角田宇宙センターに蓄積されている豊富なエンジン試験結果からエンジン性能向上に繋がる因子を抽出するためエンジン試験結果と対比する CFD を行い, またエンジン性能向上に繋がるものとして考案されたエンジン形態を CFD で試すこと.
参考 URL:https://jaxa.repo.nii.ac.jp/?action=pages_view_main&active_action=repository_view_main_item_det
ail&item_id=9276&item_no=1&page_id=13&block_id=21

JAXA スーパーコンピュータを使用する理由と利点
角田宇宙センターではスクラムジェットエンジンの研究を進めており, ラムジェットエンジン試験設備(RJTF)を用いてエンジン性能試験を多数重ねて来ている. その過程では, 特に飛行マッハ 6 条件
下ではエンジンの内部形状の一部の違いがエンジン性能に大きな影響を与えることが分かっている.
エンジンの主要素であるインレット, 分離部, ストラット等の内部形状の違いがエンジン性能にどの様に影響するのか, 当角田センターに蓄積されている豊富な実験データに基づき各形態のシミュレーションを CFD にて遂行すると共に, また試験未実施の形態についてもシミュレーションを行っている.
エンジン内部形状の空気力学的効果を CFD 援用により追求し, 実験データとの対比を重ね, 空力的効果を体系化し, スクラムジェットエンジンの設計の判断材料を構築する.

- 319 -JAXA スーパーコンピュータシステム利用成果報告(2022 年 2 月～2023 年 1 月)

今年度の成果角田宇宙センターに設置されているラムジェット試験設備(RJTF)にてエンジン試験実施済み形態と対比すべく, 改良型のエンジン形態の仮想実験を進めている. この形態は, スクラムエンジン流路の中央に設置される支柱(ストラット)の後縁を絞り込み且つ切り落としたボートテイル形のもので, エンジン試験実施済みのエンジンと基本寸法は同一ながらエンジン推力性能の改善を図ったものである.
図 1 にエンジンの概要を, 図 2 には各形態を示す. この両者に就いてエンジン内部諸量の違いを比較した. 比較は 3 次元燃焼条件にて行い, エンジン内の衝撃波等の影響を調べることで行った.
計算には汎用熱流体解析ソルバ Fluent を用いた. 最少格子幅はストラット前縁周辺等の 0.1mm である. 計算ではエンジン左右中央断面を対称面とする右舷半裁を対象とし, 格子数は 503 万, 対流項評価の数値流束は AUSM+, 制限関数は二次精度, 時間積分は陽的解法, 乱流モデルは k-ω である. 壁面は断熱壁と仮定した. 気流条件は, 流入マッハ数 5.3 等 RJTF エンジン試験条件に揃えている. 機体下面境界層に相当する設備境界層も込みである. 燃焼計算に関しては, 反応モデルは Petersen andHanson の論文(1999 年)を元に水素-酸素の反応式を Fluent に読み込ませたものを用いている. 燃焼モデルはスーパーコンピュータ"JSS3"を用いた. 角田からは遠隔利用である.
1)次に図 3 を示し, 燃焼下の OH 分布に就いて述べる. この図 3 は燃焼条件下のストラット後流域における OH 分布を示すものである. 図 3a が 5/5 高さストラット形態, 図 3b がボートテイル形ストラット形態である.
両形態とも後流域の天板側に OH の濃度の高い領域が現れている. ここは天板と側板の隅に出来ている剥離領域である. また, カウルからの斜め衝撃波がストラット後縁を通過した直ぐの下流でもOH の発生が見られる. ボートテイル形ストラット形態の方が顕著である. 5/5 高さストラット形態(M6S43)に比べ, 赤色の領域が広い. OH の発生が多く, 燃焼が相対的に活発である様子が現れている.
実際の所, 燃焼条件の計算値を用いて求めた積分推力では, 燃料流量 48g/s で 5/5 高さストラット形態が 715N, ボートテイル形ストラット形態が 825N となっており, ボートテイル形ストラット形態の方が高い値を示している. この OH 分布の違いが性能の裏付けとなっている.
何れにしても, カウル衝撃波が反応を促進している様子が分かる. ボートテイル形ストラット形態
では 5/5 高さストラット形態に比し, カウル衝撃波が余り攪乱されておらず, 影響がより下流に及んでいることにより反応がより促進されるのではないかと考えている.
このボートテイル形ストラット形態では, 前年度の報告でも記した通り, 燃料当量比分布が改善され同時に性能が向上しているが, そのことは OH 分布の違いからも裏付けられる. こちらの形態に優位性がある.
2)格子依存性を見る為に当計算に用いた格子に対して細密化した格子 A(40x106)と細密化した格子B(10x106)を用いて気流条件の計算を行った. 各々は当計算に用いた格子 C(5x106)に対して 8 倍, 2 倍である. 各格子にて算出した圧力分布を比較した. その結果, 若干の違いが見出された.
違いが顕著な箇所は, 後ろ向き段差の周辺であり, 側板とストラットとに挟まれた狭い通路の区間である. この区間では下流に向かって, 衝撃波及び膨張波の入射及び反射が繰り返され, それらが折り重なっており, 加えて後ろ向き段差部からの剥離と再循環域を含み, それ故に計算による再現が困難な区間である. この区間の流れ場をより精度高く再現するには, 更なる研究が必要である.
以上から, 細密化格子よりも当計算で種々の結果を導いて来た格子 C が適切であることが分かる.
気流計算でも燃焼計算でも格子 C は実験結果と良好な一致を得ている. それらの結果は過去に発表済- 320 -JAXA スーパーコンピュータシステム利用成果報告(2022 年 2 月～2023 年 1 月)

みである. 当計算で用いた格子は精度としても, 計算時間や費用の点でも妥当なものである.
エンジン内の空気力学や燃焼学等の諸現象の中から推力発生の源泉を見出し, これを敷衍することを究極の目標と考える. 源泉を抑えることで, 設計に生かすことが可能となる.

図 1: 供試体エンジンの概要 試験台には天地逆に設置, 5/5 高さストラット装備した姿

図 2: CFD 比較向けストラット二形態 a)5/5 高さストラット形態は試験済み形態, b)ボートテイ
ル形ストラット形態は改良型仮想形態(単位:mm)- 321 -JAXA スーパーコンピュータシステム利用成果報告(2022 年 2 月～2023 年 1 月)

図 3: 燃焼条件下におけるストラット後流域の OH 分布

成果の公表-査読なし論文
佐藤 茂 福井正明 宗像利彦 渡邉孝宏 髙橋正晴 井上 拓
スクラムジェットエンジン性能向上に関する試み―ストラット形状による燃料当量比分布改善案
第 54 回流体力学講演会/第 40 回航空宇宙数値シミュレーション技術シンポジウム(令和 4 年 6 月7 月盛岡)論文集, 令和 5 年 2 月, 宇宙航空研究開発機構特別資料
宇宙航空研究開発機構 JAXA-SP-22-007-口頭発表1)佐藤 茂 福井正明 宗像利彦 渡邉孝宏 髙橋正晴, 井上 拓
スクラムジェットエンジン性能向上に関する試み―ストラット形状による燃料当量比分布改善案
第 54 回流体力学講演会/第 40 回航空宇宙数値シミュレーション技術シンポジウム 2022, 令和 4 年6 月-7 月, 盛岡

JSS 利用状況⚫ 計算情報プロセス並列手法

熱流体解析汎用ソフト FLUENT 内部の並列処理に依拠

スレッド並列手法

熱流体解析汎用ソフト FLUENT 内部の並列処理に依拠

プロセス並列数

4 - 36

１ケースあたりの経過時間

1680 時間- 322 -JAXA スーパーコンピュータシステム利用成果報告(2022 年 2 月～2023 年 1 月)

⚫ JSS3 利用量
総資源に占める利用割合※1 (%)： 0.08

内訳
計算資源
計算システム名CPU 利用量(コア・時)資源の利用割合※2 (%)

TOKI-SORA

0.00

0.00

TOKI-ST

208,947.47

0.21

TOKI-GP

0.00

0.00

TOKI-XM

0.00

0.00

TOKI-LM

648.650.04

TOKI-TST

0.24

0.00

TOKI-TGP

0.00

0.00

TOKI-TLM

0.00

0.00

ファイルシステム資源
ファイルシステム名

ストレージ割当量(GiB)

資源の利用割合※2 (%)

/home

53.55

0.05

/data 及び/data2

446.240.00/ssd

100.950.01

アーカイバ資源
アーカイバシステム名

利用量(TiB)

資源の利用割合※2 (%)

J-SPACE

0.270.00※1 総資源に占める利用割合：3 つの資源(計算,ファイルシステム,アーカイバ)の利用割合の加重平均※2 資源の利用割合：対象資源一年間の総利用量に対する利用割合⚫ ISV 利用量
ISV ソフトウェア資源

ISV ソフトウェア(合計)

利用量(時)

資源の利用割合※2 (%)

5,611.30

3.90※2 資源の利用割合：対象資源一年間の総利用量に対する利用割合- 323 -