セル型サブクリティカル・カセット設計

DT融合ハイブリッド炉における中性子束不均一性と材料劣化の同時緩和

要旨

私は、DT核融合ハイブリッド炉の実現性を左右する本質的課題は、単に「中性子を増やすこと」ではなく、不均一な14MeV中性子源を、空間的・時間的に制御された燃焼場と材料損傷場へ変換することにあると考える。トカマク型炉では、プラズマ形状、トロイダル幾何、インボード／アウトボード差、ダイバータ・ポート・遮蔽構造により、第一壁およびブランケットに入射する中性子束が不可避的に不均一となる。さらに融合‐核分裂ハイブリッドでは、DT中性子に加えて、サブクリティカル燃料カセット内で発生する核分裂中性子が重畳し、第一壁のdpa、He/H生成、ガンマ発熱、燃料燃焼分布を非線形に変化させる。既往研究でも、ハイブリッド炉の第一壁では純粋核融合炉よりもdpaやHe生成が増大し得ることが示されている。

本稿で私は、各ブランケットセルを「燃料集合体」ではなく、中性子スペクトルを局所的に整形する交換可能なサブクリティカル・カセットとして再定義する。そして、各セル前段に配置するBe/C/Gd中性子整形層を、入射中性子束、燃焼度、dpa、ガンマ発熱、局所未臨界度に応じて空間的に変化させる設計を提案する。結論を先取りすれば、燃焼ムラの平坦化はBe/C/Gd組成制御でかなり達成できるが、第一壁dpaの完全平坦化は整形層だけでは不可能である。dpaには、燃料側からの反射・二次中性子で制御できる成分と、プラズマから直接入射する14MeV中性子でほぼ制御不能な成分があるためである。ここを見誤ると、燃料燃焼は綺麗に平坦化できても、第一壁寿命は局所ピークで決まってしまう。

1. 序論：ハイブリッド炉の問題は「平均値」ではなく「局所値」である

融合‐核分裂ハイブリッド炉は、融合コアを外部中性子源とし、その周囲にサブクリティカル核分裂ブランケットを置く構想である。米国のハイブリッド炉検討報告でも、融合‐核分裂ハイブリッドは「中性子を発生する融合コア」と「それに駆動される核分裂ブランケット」から成るサブクリティカル系として整理されている。

しかし、私が問題にしたいのは、従来設計でよく用いられる「平均TBR」「平均エネルギー増倍率」「平均dpa」ではない。実機を壊すのは平均値ではなく、局所ピークである。高束セルでは燃料が先に燃え、局所発熱とFP蓄積が進み、Gdなどの吸収材も早く焼失する。一方、低束セルでは燃料が残り、燃焼度が遅れる。これにより、時間が経つほどパワー分布が再び歪む。さらに高束領域の第一壁は、局所dpa、He生成、熱応力、照射脆化が先行し、カセット全体の交換時期を支配する。

DT融合発電プラントでは、第一壁、燃料系、ブランケット、冷却、遮蔽、トリチウム管理が一体で設計される必要がある。IAEAは、融合発電プラントには真空容器、燃料サイクル、エネルギー抽出、冷却、遮蔽、トリチウム・放射化物の閉じ込めが必要であり、第一壁は高熱流束・中性子束に耐える必要があると整理している。

したがって、私の出発点は次の一点である。ハイブリッド炉の中性子設計は、平均増倍率の設計ではなく、セルごとの中性子インピーダンス整合問題である。

2. 現状の問題点

2.1 トカマクでは、第一壁中性子束は本質的に不均一である

トカマク炉の第一壁中性子束は、プラズマからの一次入射束だけでは決まらない。古典的だが今なお重要なJAERIのモンテカルロ計算は、トカマク第一壁の14MeV中性子束のポロイダル分布について、一次入射束と実際の束分布が大きく異なり得ること、また第一壁被覆材やブランケット反射を考慮する必要があることを示している。特に、実際の14MeV束は一次束に単純比例せず、非円形断面・複雑幾何ではモンテカルロ法が有効であると結論づけている。

この事実は、セル型カセット設計にとって決定的である。なぜなら、各セルの燃焼度やdpaを評価するとき、単純な幾何学的壁負荷だけでBe/C/Gd組成を決めると、実際の中性子応答を外すからである。各セルには、一次DT中性子、第一壁・被覆材で減衰した14MeV成分、ブランケット内で散乱・減速した成分、燃料から戻る核分裂中性子、捕獲ガンマが同時に入る。

2.2 ハイブリッド炉では、第一壁損傷が純核融合炉より複雑になる

融合‐核分裂ハイブリッド炉では、第一壁の前面からDT中性子が入り、背面側のサブクリティカル燃料から核分裂中性子とガンマ線が戻る。Khripunovの解析は、ハイブリッド第一壁材料では、同一DT中性子壁負荷でも純融合系よりdpaが約2倍、鋼材中のHe生成が約4〜5倍になり得ること、また核発熱が大きく増加し得ることを示している。

ここから私が導く結論は明確である。ハイブリッド炉では、第一壁は“核融合側の部品”ではなく、“融合中性子と核分裂中性子の両側から叩かれる境界部品”である。 したがって、燃料側に置く中性子整形層は、燃料燃焼だけでなく、第一壁へ戻る反射中性子・捕獲ガンマ・二次中性子を抑える役割を持たせなければならない。

2.3 Be/C/Gdは有効だが、役割を混同してはいけない

Beは14MeV領域で中性子増倍材として働き、低束セルで燃料反応率を引き上げる候補になる。2025年のCERN n_TOF提案でも、ブランケット中のLiおよびBe反応、特に 9Be(n,2n)^{9}\mathrm{Be}(n,2n)9Be(n,2n) 反応の断面積不確かさがTBR評価に大きく影響するため、高精度測定が重要とされている。

Cは、単なる構造材ではなく、スペクトルを鈍らせる「中性子粘性」として使うべきである。高速成分を減速し、Gdが効く熱・エピサーマル領域へ中性子を落とす。ただし、Cの使い方を誤ると、特定燃料の熱中性子核分裂を増やし、むしろ局所出力を上げる可能性がある。

Gdは極めて強い吸収材であり、特に 155Gd^{155}\mathrm{Gd}155Gd と 157Gd^{157}\mathrm{Gd}157Gd が重要である。OECD/NEAは、Gdが原子力産業で非常に大きな熱中性子吸収断面積を持つ吸収材として用いられ、燃焼性吸収材や停止毒として使われる一方、Gd断面積の不確かさが高燃焼度時の反応度評価に大きく効くと整理している。

したがって、私の設計原則は次の通りである。

低束セルにはBeを増やす。高束セルにはCを増やし、Cで減速された成分に対して少量Gdを効かせる。Gdを14MeV中性子の直撃位置に置いても効率は悪く、捕獲ガンマと局所発熱だけが増える。Gdは、Cまたは他の減速領域の後段に薄く分散配置すべきである。

3. 私が提案するセル型サブクリティカル・カセット

私は、ハイブリッドブランケットを連続した一体ブランケットではなく、ポロイダル方向・トロイダル方向に分割されたスペクトル制御セル群として設計する。各セルは、次のような機能層を持つ。

Plasma→First Wall / Armour→Be/C/Gd Shaping Layer→Subcritical Fuel→Tritium Breeder→Shield\text{Plasma} \rightarrow \text{First Wall / Armour} \rightarrow \text{Be/C/Gd Shaping Layer} \rightarrow \text{Subcritical Fuel} \rightarrow \text{Tritium Breeder} \rightarrow \text{Shield}Plasma→First Wall / Armour→Be/C/Gd Shaping Layer→Subcritical Fuel→Tritium Breeder→Shield

重要なのは、Be/C/Gd層を「固定材料」としないことである。各セル iii について、入射中性子束を sis_isi​、Be体積率を xBe,ix_{\mathrm{Be},i}xBe,i​、C体積率を xC,ix_{\mathrm{C},i}xC,i​、Gd等価体積率を xGd,ix_{\mathrm{Gd},i}xGd,i​ とし、

xBe,i+xC,i+xGd,i=1x_{\mathrm{Be},i}+x_{\mathrm{C},i}+x_{\mathrm{Gd},i}=1xBe,i​+xC,i​+xGd,i​=1

とする。ただし、Gdは実際にはGd2_22​O3_33​、Gd含有セラミック、または薄膜・分散相として数％以下に抑えるべきである。

私の基本制御則は、概念的には次のようになる。

xBe,i↓assi↑x_{\mathrm{Be},i} \downarrow \quad \text{as} \quad s_i \uparrowxBe,i​↓assi​↑xGd,i↑assi↑x_{\mathrm{Gd},i} \uparrow \quad \text{as} \quad s_i \uparrowxGd,i​↑assi​↑xC,i=1−xBe,i−xGd,ix_{\mathrm{C},i} = 1-x_{\mathrm{Be},i}-x_{\mathrm{Gd},i}xC,i​=1−xBe,i​−xGd,i​

すなわち、高束セルではBeを抜き、Cを厚くし、Gdを少量入れる。低束セルではBeを増やし、Gdを入れない。

これは単なるパワーフラットニングではない。私はこれを、中性子束の局所インピーダンス整合と呼ぶ。高束セルでは中性子を吸収・減速・後方散乱抑制し、低束セルでは中性子を増倍・有効利用する。これにより、サブクリティカル燃料の燃焼度分布と第一壁への戻り中性子成分を同時に制御する。

4. シミュレーションによる定量実証の設計

高忠実度実証には、3D CADベースの連続エネルギーモンテカルロ輸送計算が必要である。近年、OpenMC/DAGMCを用いたARC級炉の3D中性子輸送・放射化解析では、DPA評価には3Dモデルが必要であることが強調されている。

また、MCNP、Serpent、OpenMC、DAGMCなどの比較・ベンチマーク研究では、複雑な核融合幾何に対してCADベース輸送ワークフローが有効であり、今後の検証スイート統合が必要とされている。

私は、次の評価量をセルごとにタリーする。

燃料出力：

Pi(t)=∫Eϕi(E,t)∑mNm,i(t)σf,m(E)Ef,m dEP_i(t)=\int_E \phi_i(E,t)\sum_m N_{m,i}(t)\sigma_{f,m}(E)E_{f,m}\,dEPi​(t)=∫E​ϕi​(E,t)m∑​Nm,i​(t)σf,m​(E)Ef,m​dE

燃焼度：

Bi(T)=1MHM,i∫0TPi(t) dtB_i(T)=\frac{1}{M_{HM,i}}\int_0^T P_i(t)\,dtBi​(T)=MHM,i​1​∫0T​Pi​(t)dt

第一壁dpa：

Di(T)=∫0T∫EϕiFW(E,t)σdpa(E) dE dtD_i(T)=\int_0^T\int_E \phi^{FW}_i(E,t)\sigma_{\mathrm{dpa}}(E)\,dE\,dtDi​(T)=∫0T​∫E​ϕiFW​(E,t)σdpa​(E)dEdt

最適化関数：

J=wPCV(P)2+wBCV(B)2+wDCV(D)2+wQPF(q′′′)2+λkmax⁡(0,keff−klim)2+λTmax⁡(0,TBRreq−TBR)2J = w_P CV(P)^2 + w_B CV(B)^2 + w_D CV(D)^2 + w_Q PF(q''')^2 + \lambda_k \max(0,k_{\mathrm{eff}}-k_{\mathrm{lim}})^2 + \lambda_T \max(0,TBR_{\mathrm{req}}-TBR)^2J=wP​CV(P)2+wB​CV(B)2+wD​CV(D)2+wQ​PF(q′′′)2+λk​max(0,keff​−klim​)2+λT​max(0,TBRreq​−TBR)2

ここで、CVCVCV は変動係数、PFPFPF はピーキング係数である。私は、単純な「出力平坦化」ではなく、燃焼度、dpa、ガンマ発熱、TBR、未臨界余裕を同時に最適化する。

さらに、MOOSE/TMAP8のようなマルチフィジックス環境を接続する必要がある。2025年のMOOSEブランケットシミュレーション研究では、OpenMCに基づく中性子発熱、3D熱伝導、1D熱流動、3Dトリチウム輸送を連成し、NQA-1品質保証に沿う方向が示されている。

5. 概念実証シミュレーション：Be/C/Gd勾配設計の効果

以下の数値は、私が提案する設計思想を定量的に示すための簡易多群応答面モデルであり、最終設計値ではない。実機評価では、同じ最適化問題をOpenMC、MCNP、Serpent、FISPACT、ORIGEN、MOOSE等で置き換える必要がある。ただし、設計傾向を掴むには有効である。

前提として、24個のポロイダルセルを置き、入射DT中性子壁負荷のピーキング係数を2.09、変動係数を0.266とした。均一整形層をBe 30%、C 69.7%、Gd 0.3%相当とした場合、燃料出力分布もdpa分布も、ほぼ入射中性子束の不均一性をそのまま引き継ぐ。

表1　簡易応答面モデルによる最適化結果

ケース	燃料出力CV	燃料出力ピーキング	第一壁dpa CV	第一壁dpaピーキング	平均出力
均一Be/C/Gd層	0.266	2.09	0.266	2.09	1.000
燃料側整形層のみ最適化	0.082	1.28	0.231	1.93	0.993
第一壁近傍一体型整形を含む最適化	0.070	1.25	0.133	1.53	0.993

この結果から、私は三つの結論を得る。

第一に、燃料出力の平坦化はかなり可能である。 高束セルにCとGdを増やし、低束セルにBeを増やすだけで、燃料出力CVは0.266から0.07〜0.08程度まで低下する。

第二に、第一壁dpaの平坦化は難しい。 燃料側の整形層だけでは、dpa CVは0.266から0.231程度にしか下がらない。これは、第一壁dpaの大きな部分がプラズマから直接来る14MeV中性子で決まるためである。

第三に、第一壁近傍のC/Be構造を整形層に含めると、dpa CVは0.133程度まで下げられる。 ただし、これはプラズマ対向材、熱除去、トリチウム保持、スパッタリング、遠隔交換性と強く干渉するため、燃料側だけの問題として扱ってはいけない。

表2　最適化後の組成傾向

セル群	正規化入射束	Be体積率	C体積率	Gd等価体積率	設計意図
低束セル	0.70〜0.79	0.70〜0.75	0.25〜0.30	0	中性子増倍で燃焼を追いつかせる
中束セル	0.80〜1.05	0.25〜0.60	0.40〜0.75	0	出力とdpaのバランス調整
高束セル	1.18〜1.47	0.02〜0.05	0.93〜0.98	0.008〜0.022	減速・吸収・後方反射抑制

この組成分布は直感に合う。高束セルにBeを置くと、燃料出力と戻り中性子がさらに増える。したがって高束セルではBeを最小化し、Cでスペクトルを鈍らせ、Cで落ちた熱・エピサーマル成分をGdで刈り取る。低束セルではGdを入れず、Beを増やして燃焼を補償する。

6. 私の独自視点：dpa平坦化には「中性子整形層」だけでは足りない

ここが本稿の最も重要な主張である。

燃料燃焼ムラは、Be/C/Gdの局所組成最適化でかなり抑えられる。しかし、第一壁dpaの完全平坦化は、燃料側整形層だけでは原理的に限界がある。理由は、dpaが次の二成分で構成されるからである。

Di=Di,direct14MeV+Di,albedo/fissionD_i = D_{i,\mathrm{direct}}^{14MeV} + D_{i,\mathrm{albedo/fission}}Di​=Di,direct14MeV​+Di,albedo/fission​

Di,albedo/fissionD_{i,\mathrm{albedo/fission}}Di,albedo/fission​ は燃料側整形層である程度制御できる。しかし、Di,direct14MeVD_{i,\mathrm{direct}}^{14MeV}Di,direct14MeV​ はプラズマ形状、第一壁位置、遮蔽前の幾何でほぼ決まる。したがって、Be/C/Gd層だけでdpaを完全に平坦化しようとすると、過剰なGd、過剰なC、過剰な局所発熱、TBR低下、燃料利用低下を招く。

私の解決策は、二階建てである。

第一段階では、Be/C/Gd整形層で燃料燃焼と戻り中性子を平坦化する。

第二段階では、第一壁そのものについて、セル別の交換周期、犠牲層厚さ、RAFM/SiC/W/Be/Cの局所構造、冷却余裕、遠隔交換単位を変える。

つまり、燃焼平坦化は中性子工学で解き、dpa平坦化は中性子工学と保守工学の連成で解くべきである。ここを一つの最適化問題に押し込むことが、セル型サブクリティカル・カセット設計の本質である。

7. 実現化に向けた研究ロードマップ

私は、実証を五段階で進めるべきだと考える。

第一に、OpenMCまたはSerpentで単一セルのBe/C/Gd応答曲面を作る。ここでは、Be/C/Gd比率、層厚、燃料種類、冷却材、Gd配置位置を振り、燃料出力、dpa、ガンマ発熱、TBR、局所 keffk_{\mathrm{eff}}keff​ を評価する。

第二に、24〜96セル程度のトカマク3D部分モデルを作る。ここでは、現実のポート、ダイバータ、インボード遮蔽、アウトボード厚さ差を入れる。OpenMC/DAGMCのようなCADベース輸送を使うべきである。

第三に、燃焼計算を入れる。Gdは燃える。燃料も燃える。FPが貯まる。したがって、BOCだけで平坦でも意味がない。EOCまでの Pi(t)P_i(t)Pi​(t)、Bi(t)B_i(t)Bi​(t)、Di(t)D_i(t)Di​(t) を追う必要がある。

第四に、熱流動と構造を接続する。Gdで中性子を吸えば捕獲ガンマが出る。Cを増やせば熱伝導、照射寸法変化、トリチウム保持が問題になる。Beを増やせば増倍は有利だが、ヘリウム生成、膨張、毒性、製造・廃棄が重くなる。

第五に、カセット交換性を評価する。高束セルを低束セルと同じ寿命に合わせようとすると、組成が過激になる。私はむしろ、高束セルは短寿命交換、低束セルは長寿命使用と割り切る方が、全体最適になる可能性が高いと見る。

8. 結論

私は、DT融合ハイブリッド炉のブランケットを、均一な燃料領域ではなく、セルごとにスペクトルを設計するサブクリティカル・カセット集合体として扱うべきだと結論づける。

Be/C/Gd整形層の基本方針は明確である。

低束セルにはBeを増やす。高束セルにはCを増やし、Gdを少量入れる。GdはCによる減速後に効かせる。高束セルのBeは最小化する。燃料出力平坦化とdpa平坦化を同じ制御量で完全に達成しようとしない。

概念応答面シミュレーションでは、入射中性子束ピーキング2.09の条件に対し、Be/C/Gd最適化により燃料出力ピーキングを約1.25〜1.28まで低減できる可能性を示した。一方、第一壁dpaピーキングは燃料側整形層だけでは約1.93までしか下がらず、第一壁近傍一体型整形を含めても約1.53が残った。これは、dpa平坦化の限界が中性子整形層ではなく、プラズマ由来14MeV直接成分に支配されることを意味する。

したがって、私の最終提案は次である。

セル型サブクリティカル・カセットは、燃料燃焼を平坦化するためのBe/C/Gdスペクトル制御と、第一壁寿命を平坦化するための交換周期・犠牲層・局所構造設計を統合した「中性子‐材料‐保守連成システム」として設計すべきである。

この方向に進めば、DT融合ハイブリッド炉は、単なる中性子増倍率競争から、局所燃焼、局所損傷、局所保守を制御する実機工学へ進むことができる。私はここに、融合ハイブリッド炉実現化の最も現実的な研究テーマがあると考える。