核データ処理：
評価済み核データから断面積ライブラリへ
第45回炉物理夏期セミナー
2013年8月1日
琵琶湖 白浜荘

長家 康展
炉物理研究グループ
原子力基礎工学研究部門
日本原子力研究開発機構

1

本講義の狙い• 評価済み核データを理解すること– 評価済み核データはどのような形式で与えられるのか？

• 核データ処理を理解すること– 核計算コード用の断面積ライブラリはどのようにして作成されるのか？

2

アウトライン• 核データ処理について• 評価済み核データのフォーマット（ENDF-6）• ポイントワイズ断面積データの作成– 断面積の線形化– 共鳴断面積の再構成– ドップラー拡がり• 断面積ライブラリの作成– 連続エネルギーモンテカルロ用– 多群輸送計算コード用• 最近の核データ処理に関する話題– 一般化核データ（GND）– 核データ処理コードシステム– 国産核データ処理コードシステムの開発

3

核データ処理とは？
汎用目的の評価済み核データを、限定目的の核データセット（ライブラリ、データベース）へ変換する処理。
JEFF
ENDF/B
JENDL

核データ
処理

MVPライブラリ
MCNPライブラリ
SRACライブラリ
SLAROM-UＦライブラリ

評価済み核データ
ENDF

輸送計算コード用
多群断面積ライブラリ4

評価済み核データのフォーマット
• ENDF-6フォーマット
• 元々、米国の評価済み核データENDF/Bのために開発されたフォーマット
• ENDF/B-VI（バージョン）⇔ENDF-6 （フォーマット）• 米国断面積評価ワーキンググループ(CSEWG)によって管理• マニュアル ⇒ ENDF-102
 http://www.nndc.bnl.gov/endf
 ENDF-102, BNL-90365-2009 Rev.2
 最新版 ⇒ 2011 Dec

5

ENDF-6フォーマットのデータ構造• tape : 1つまたは複数のmaterial。計算機上の1つのファイル。
• material : 1つの核種もしくはcompound。
MAT番号で区別。
• file : 断面積、角度分布、エネルギー分布等の種類を表す。 MF番号で区別。
• section : 反応を表す。MT番号で区別。
• MAT番号のルール– 安定核種で最も軽い同位体をZ25。Z=原子番号。
– 質量数が1増えるごとに3プラス。
– 例：H-1 ⇒ 125, H-2 ⇒ 128, U-235 ⇒ 9228

6

代表的なMF番号• MF=1: コメント、ν値• MF=29: 共分散コメント• MF=2: 共鳴パラメータ • MF=30-36: 共分散データ• MF=3: 反応断面積• MF=4: 二次粒子の角度分布• MF=5: 二次粒子のエネルギー分布• MF=6: 二次粒子エネルギー角度分布• MF=7: 熱中性子散乱データ• MF=8: 放射性崩壊&FP収率データ• MF=9: 放射性崩壊の生成多重度• MF=10: 放射性崩壊の生成断面積• MF=11: 光子生成データコメント• MF=12-15: 光子生成データ

7

代表的なMT番号• MT=1: 全断面積
• MT=2: 弾性散乱
• MT=16: (n,2n)反応• MT=18: 核分裂• MT=51-90: 離散レベル非弾性散乱• MT=91: 連続レベル非弾性散乱• MT=102: 放射捕獲
• MT=151: 共鳴パラメータ• MT=452: 核分裂当りの平均全中性子数（ν値）8

ENDF-6フォーマットの基本データ構造
MAT MF MT line

ファイル1

ファイル2

ファイル3

tape id
Start of MF1, MT451 (コメント)開始...
SEND record
FEND record
Start of MF2, MT151 (共鳴パラメータ)...
SEND record
FEND record
Start of MF3, MT1 (全断面積)...
SEND record
Start of MF3, MT2 (弾性散乱断面積)...
FEND record
MEND record
TEND record
66カラム

0 0
1111 1451

11111 1 01111 1 099999
1111 2151
11111 21111 21111 300999991
11111 31111 30
21111 30 0-1 0099999
0
0

4 2 3

5カラム

9

ENDF-6フォーマット断面積の具体例
JENDL-4.0, Fe-56の(n,2n)断面積

MF
MAT MT

2.605600+4 5.545440+100002631 3 16
1
-1.120270+7-1.120270+7001112631 3 16
2
11
2
0
0
0
02631 3 16
31.140470+7 0.000000+0 1.170000+7 1.622410-2 1.200000+7 4.800450-22631 3 16
41.300000+7 2.138200-1 1.400000+7 3.891650-1 1.500000+7 5.134000-12631 3 1651.600000+7 5.817500-1 1.700000+7 6.107500-1 1.800000+7 6.118000-12631 3 1661.900000+7 5.977000-1 2.000000+7 5.759000-1
2631 3 16
72631 3 099999
11カラム
11カラム
11カラム
11カラム
11カラム
11カラム

11カラム実数表現•±1.234567±n•±1.23456±nn (nn≦38)•±1.23456789

単精度（32ビット精度）

10

ENDF-6フォーマットのデータ表現形式• テーブル形式– 離散化点を与える形式 ( E1 , σ 1 ), ( E2 , σ 2 ), 
– 汎用的、関数形に依存しない。

• 関数形式– 関数形を決めて、パラメータを与える形式– 例：ブライト・ウィーグナー一準位公式– パラメータが別のパラメータ（例えば、エネル
ギー）に依存する場合、そのパラメータはテーブル形式で与えられる。

11

テーブル形式のデータ表現
NBT(1)=3

NP = 10y(1)(ポイント数)
NR = 3(領域数)
ENDFに与えられている

NBT(2)=7

NBT(3)=10

領域1

領域2

領域3

INT(1)=5

INT(2)=2

INT(3)=2log-log内挿

線形内挿

線形内挿

y(6)

y(2)

y(7)
y(5)

y(9)

y(8)

y(4)

y(10)

y(3)
x(1) x(2)

x(3)
x(4)

x(5)
x(6)

x(7) x(8) x(9)x(10)◆点を与えて、線でつなぐ（内挿法を指定）

12

ENDF-6で用いられる内挿法• INT=1: ヒストグラム（一定）• INT=2: 線形内挿（linear-linear）• INT=3: y – ln(x)線形内挿（linear-log）• INT=4: ln(y) - x線形内挿（log-linear）• INT=5: ln(y) – ln(x)線形内挿（log-log）• INT=6: 荷電粒子断面積のための特別内挿法• INT=11-15: method of corresponding points（2次元関数内挿）• INT=21-25: ユニットベース内挿（2次元関数内挿）

13

ポイントワイズ断面積データの作成
ENDF
断面積の線形化（linearization）
共鳴断面積の再構成(reconstruction)
ドップラー拡がりの考慮(Doppler broadening)
PENDF• 連続エネルギー用MCライブラリ作成• 多群輸送コード用ライブラリ作成

共通の核データ処理

14

断面積データの線形化
JENDL-4.0, H-1, File 3
Total
ZA
AWR
1.001000+3 9.991673-10000 125 3 1
NR= 2
NP= 96 125 3 1
0.000000+0 0.000000+000NBT(1)= 30
INT(1)= 5 NBT(2)= 96
INT(2)= 200 125 3 1
1.000000-5 3.713628+1 2.000000-5 3.224498+1 5.000000-5 2.790478+1 125 3 1
E1
E2
E3
σ1
σ2
σ3

Elastic
ZA
AWR
1.001000+3 9.991673-10000 125 3
NR= 1
NP= 96 125 3
0.000000+0 0.000000+000NBT(1)= 96
INT(1)= 2 NBT(2)= 0
INT(2)= 000 125 31.000000-5 2.043634+1 2.000000-5 2.043634+1 5.000000-5 2.043634+1 125 3
E1
E2
E3
σ1
σ2
σ3
Capture

2
2
2
2

ZA
AWR
1.001000+3 9.991673-10000 125 3102
NR= 2
NP= 96 125 3102
2.224631+6 2.224631+600NBT(1)= 30
INT(1)= 5 NBT(2)= 96
INT(2)= 200 125 3102
1.000000-5 1.669994+1 2.000000-5 1.180864+1 5.000000-5 7.468441+0 125 3102
E1
E2
E3
σ1
σ2
σ3

15

線形化の必要性

16

区間二分法(interval halving)17

区間二分法(interval halving)18

区間二分法(interval halving)19

区間二分法(interval halving)

20

断面積データの線形化（利点）• ドップラー拡がり断面積をカーネル拡張法(kernel broadening)により精度よく計算できる。
• コード内での処理を線形内挿だけの処理に特化することができる。
• プロットが簡単。

21

共鳴断面積の再構成
軽核

重核

ファイル3(MF=3)で与えることが可能

J-4.0 U-235
MF=3

断面積は、共鳴パラメータ(MF=2, MT=151)で与えられる◆共鳴の再構成が必要

22

ENDF-6で用いられる共鳴公式• ブライト・ウィーグナー一準位公式– 表現が簡単、解析的に（ψχ法で）ドップラー拡がりを考慮可。負の断面積を与えることがある。

• ブライト・ウィーグナー多準位公式– 表現が複雑。負の断面積を与えない。多チャンネル
効果は無視。

• ライヒ・ムーア(Reich-Moore)– 2つの核分裂チャンネルをもつ多準位公式。

• 限定R行列(R-Matrix Limited)– 現在のところ、ENDF/B-VII.1のF-19でのみ使用され
ている。

23

ブライト・ウィーグナー一準位公式
ENDF-102における放射捕獲反応に対する定義式
NLS−1

σ n ,γ ( E ) = ∑ σ n,γ ( E ),
 =0

σ


n ,γ

Γnr ( E )Γγr
(E) = 2 ∑ g J ∑
,
2
2
k J
r =1 ( E − Er′ ) + Γr / 4

Γnr ( E ) =

π

NR J

P ( E )Γnr ( Er )Er′ = Er +P ( Er )

2mn A
, k=
E
 A +1S  ( Er ) − S  (E )2 P ( Er )2J +1
Γnr ( Er ), g J =2(2 I + 1)

24

部分波展開
入射中性子 ⇒ 平面波、部分波で展開∞ψ in = eikz = ∑ (2 + 1)i  j (kr ) Pl (cos θ ) =0

ラマーシュ, 「原子炉の初等理論」（上）, p.58, 第2－14図より引用

25

ブライト・ウィーグナー一準位公式
簡略化 ⇒ s波(  = 0 )のみ、孤立共鳴GGr

σ γr ( E ) = σ r
共鳴パラメータ

ERr

Γγr

Er

1
2
, xr = ( E − E r )2
E 1 + xr
Γr

Γr
GTr

GNr

AWR

AJ

2π  A + 1  Γnr 2 J + 1
σr =mn Er  A  Γr 2(2 I + 1)
2

2

SPI

26

ブライト・ウィーグナー一準位公式
ENDF/B-VII.0のPu-238

SLBW

分離共鳴

ZA
AWR
9.423800+4 2.360045+2
0
ZAI= 9.423800+4 1.000000+0 =ABN0
EL= 1.000000-5 2.000000+2 =EH LRU= 10
LFW= 1
LRF= 1

NIS= 1
NER= 2
NRO= 009434 2151
09434 2151
NAPS= 09434 2151

Subsection for 1st energy range

SPI= 0.000000+0 9.309000-1 =AP
NLS= 10009434 2151
AWRI= 2.360045+2 0.000000+0 =QX
LRX= 0 6*NRS= 96
NRS= 169434 2151
L= 0
-1.000000+1 5.000000-1 6.217000-2 1.581000-2 4.500000-2 1.360000-39434 2151
ER1
AJ1
GT1
GN1
GG1
GF1
-4.000000-1 5.000000-1 4.670000-2 3.400000-4 4.500000-2 1.360000-39434 2151
2.855000+0 5.000000-1 3.808600-2 7.470000-5 3.680000-2 1.211000-39434 2151
ER3
AJ3
GT3
GN3
GG3
GF3
J値

Γ

Γn

…

共鳴
エネルギー

Γγ

Γf

27

共鳴断面積の再構成
E3 = 2.855 eV
Γ = 3.8086e - 2
Γn = 7.47e - 5
Γγ = 3.68e - 2

28

共鳴断面積の再構成
E4 = 9.975 eV
Γ = 3.7213e - 2
Γn = 2.084e - 4
Γγ = 3.024e - 2

29

共鳴断面積の再構成

バックグラウンド断面積• 多準位効果を無視している影響• 考慮されていない共鳴の影響（ENDF/B-VII.0のPu-238では、バックグラウンド断面積は0）30

共鳴断面積の再構成

31

共鳴断面積の再構成

32

共鳴断面積の再構成

Γ3
2

Γ4
2

33

共鳴断面積の再構成

34

非分離共鳴データの処理• ENDF-6フォーマットの非分離共鳴領域⇒ 平均の共鳴幅、平均の共鳴間隔、
共鳴幅の分布（ある自由度のχ2分布）• 放射捕獲反応に対する断面積公式
NLS−1

σ n ,γ ( E ) = ∑ σ n,γ ( E ),
 =0

σ


n ,γ

2π 2 NJS g J Γn Γγ(E) = 2 ∑k J D , JΓ

無限希釈断面積
, J

2π 2 NJS g J Γ n, , J Γ γ , , J
Rγ , , J= 2 ∑k J D , JΓ , J

揺らぎ因子（ χ2分布の平均）• 自己遮へい因子、確率テーブルは別モジュールで計算

35

ドップラー拡がり
基礎式 ⇒ 反応率保存vσ (v, T ) = ∫ vrσ (vr ,0) P(V )dV

V

v
中性子

標的核

V
µ = cos θ

v r 相対速度v

マクスウェル・ガスモデル
β 
P(V )dV =   e − βV V 2 dφdµdV ,
π 
3/ 2

2

M
β=
2kT

vr = v + V − 2vVµ
2

2

36

SIGMA1 kernel broadening
ドップラー拡がり基礎式を変形[β ∞− β ( v −v )−β (v+v )
2
σ
dv
v
v
e
e
(
,
0
)
−
r r
r
∫
0
π

1
σ (v, T ) = 2
v

2

r

r

2

]

線形化

E − Ek
Ek +1 − E
1 2
σ ( E ,0) =σ k +1 +σ k , Er = mvr
Ek +1 − Ek
Ek +1 − Ek
2

F (a) =
n

2

π

a

∫ze
0

n −z2

dz

この積分値が計算できれば近似なく、σ (v, T ) が求まる。

37

ドップラー拡がりを考慮した断面積

JENDL-4.0
C-000 Elastic

38

ドップラー拡がりを考慮した断面積
JENDL-4.0
U-238 Capture

39

ドップラー拡がり（その他）• 実際のコード（SIGMA1, BROADR）では、断面積のドップラー拡がりに応じて、点を増やしたり、減らしたりして、内挿精度が悪くならないようにしている。
• 既にドップラー拡がりを考慮した温度（T1）から、更に別の温度（T2）へドップラー拡がりを考慮することができる。

40

連続エネルギーMC用断面積ライブラリ• 共通エネルギーグリッドの作成– 任意のエネルギー点における反応確率（σx/σt）が必要。
– エネルギー領域の検索が一度で済む。
– 断面積の線形化と共鳴の再構成と同時に行われる。

• NJOYコード– 室温における中性子エネルギー0.0253 eV– 1オーダー毎にその1倍、2倍、5倍のエネルギー点を追加する。

41

遷移確率データの作成
粒子輸送計算に必要なデータ• 断面積• 散乱後の粒子の飛行方向• 散乱後の粒子エネルギー

ENDFのファイル4, 5, 6をMC用のデータに変換
ENDF-6の微分断面積の形式

σ i ( µ , E → E ′) = σ ( E ) yi ( E ) f i ( µ , E → E ′) / 2π ,∫ dE ' ∫ dµ f (µ , E → E ′) = 1
ii はある生成物のインデックス

確率密度関数42

角度分布とエネルギー分布のタイプ• 角度分布のみ– ファイル4のみ f ( µ , E )– E’は運動学(kinematics)から計算– 弾性散乱(MT=2), 離散レベル非弾性散乱(MT=51-90)• 角度分布とエネルギー分布の積– ファイル4とファイル5 f ( µ , E ) × g ( E → E ′)– (n,2n)反応(MT=16), 核分裂反応(MT=18), 連続レベル
非弾性散乱(MT=91)• エネルギー角度分布– エネルギーと角度に相関– ファイル6 f ( µ , E → E ′)

43

ファイル4（角度分布）• ENDF-6フォーマットのフラグ(LTT番号)– LTT=0:等方– LTT=1:ルジャンドル展開形式– LTT=2:テーブル形式– LTT=3:ルジャンドルとテーブルの混合形式

44

ファイル5(エネルギー分布)• ENDF-6フォーマットのフラグ(LF番号)– LF=1:任意テーブル形式– LF=5:一般蒸発形式– LF=7:マクスウェル核分裂スペクトル– LF=9:蒸発スペクトル– LF=11:エネルギー依存ワットスペクトル– LF=12:エネルギー依存核分裂スペクトル(Madland and Nix fission spectrum)45

ファイル6(エネルギー角度分布)• ENDF-6フォーマットのフラグ(LAW番号)– LAW=1:未知の分布– LAW=1:連続エネルギー角度形式– LAW=2:離散二体散乱– LAW=3:等方離散放出– LAW=4:離散二体反応反跳分布– LAW=5:荷電粒子弾性散乱– LAW=6:N体位相空間分布– LAW=7:実験室系角度エネルギー則

46

モンテカルロ計算で必要なデータ• サンプリング ⇒ 乱数を使って、分布から具体的な値を決定する。
確率密度関数(PDF)
積算密度関数
(CDF)

p (x)

乱数 ξ

x

（サンプル値）x

C ( x) = ∫ p(t )dt , (a ≤ x < b)a

x = C −1 (ξ )• ENDF-6で用いられている分布形式– 関数形式 ⇒ 近似を用いることなく直接サンプリング– テーブル形式

47

関数形式からのサンプリング• マクスウェル分布(ファイル5, LF=5)
p ( E ′) =2

1

π T 3/ 2E ′e − E ′ / T−1  πE ′ = T − ln ξ 2 − (ln ξ 3 ) cos  ξ1 
 2 
• Kalbach系統式（ファイル6, LAW=1, LANG=2）af (µ ) =[cosh(aµ ) + r sinh( aµ )]2 sinh( a )近似を用いることなく、µのサンプリング可能

48

テーブル形式からのサンプリング
ヒストグラム
PDF

CDF

ξ

x

49

テーブル形式からのサンプリング
線形内挿
PDF

CDF

ξ

x

50

テーブル形式からのサンプリング
等確率ビン
PDF

等確率

CDF

ξ

等間隔

x

51

遷移確率データの作成（まとめ）• ENDF-6フォーマットにおける元データ– ファイル4(角度分布), ファイル5（エネルギー分布）、ファイル6 （エネルギー角度分布）– 分布 ⇒ 確率密度関数• MC計算用のデータに変換– 関数形式 ⇒ 直接サンプリング– テーブル形式 ⇒ グリッド点、確率密度関数、
積算密度関数のデータに変換

52

確率テーブルの作成• 非分離共鳴領域における自己遮へい効果• 確率テーブル法– 非均質効果、異核種間の共鳴干渉効果– MC法の計算過程において直接考慮• 確率テーブル法の仮定– 中性子は、どのエネルギーにも等確率入射– 共鳴は接近、多数の共鳴を飛び越して減速– 散乱後のエネルギーの断面積は統計的に取
扱可能

53

確率テーブルの作成
平均共鳴パラメータと分布則

I1 + I 3, P3∆E1 + ∆E3
∆E + ∆E4
σ 2 , P2 = 2∆EI∆E
σ 1 = 5 , P3 = 5
∆E5
∆E

σ3 =

ラダー

∆E1 ∆E2 ∆E3 ∆E4 ∆E5σ (barn)σj

σ (barn)

100
80

293K

60
40
20

Energy00

0.2 0.4 0.6 0.8 1
Cumulative probability

54

その他のデータ• 熱中性子散乱データ– ファイル7 ⇒ 熱中性子散乱則 S (α , β )– 遷移確率データ(PDF, CDF)へ変換• 反応率計算データ– 核発熱データ（Kerma), はじき出し断面積(DPA)• 光子生成データ– ファイル12(生成収率), ファイル13(光子断面積),14（光子角度分布）, 15（光子エネルギー分布）– ファイル3, 6でも表現

55

多群断面積ライブラリの作成
多群断面積

σ

i
, g

i
σ
∫ ( E )φ ( E )dE=g∫ φ ( E )dE, ∫ dE ≡ ∫g

E g −1Eg

dEg

σ ( E )φ ( E )dE ∫ dE ′∫ dµf ( E → E ′, µ ) P ( µ )∫=φ ( E )dE∫
feed関数i

σi
, g → g′



g

g

一般形

σ , gg′F( E )σ ( E )φ ( E )dE∫=∫ φ ( E )dE


g

g



ポイントワイズ断面積
を数値積分

56

多群断面積ライブラリの作成• 初期中性子束– 重み関数（スペクトル）をユーザ入力– マクスウェル ＋ 1/E + 核分裂スペクトル• 自己遮へい因子– Bondarenkoモデル ⇒⇒ φ ( E ) ∝W ( E )iσ⇒, g =
 +1[σ ( E ) + σ ]itσ i , g (σ 0 , T )f = i
σ  , g (σ 0 = ∞, T )i0W ( E )φ ( E ) ∝[Σ t ( E )]+1i
σ∫ ( E )φ ( E )dEg∫ φ ( E )dEg

バックグラウンド断面積と温度に対してテーブル化

57

最近の核データ処理に関する話題• 一般化核データ– GND = Generalized Nuclear Data

• 最近の核データ処理コード– NJOY2012– PREPRO2012– その他の核データ処理コード– 国産核データ処理コード

58

一般化核データ(GND)• ENDF-6フォーマット
– Fortran向き、テキストフォーマット
– 1960年代に設計、長い実績– フォーマットの制限（例、離散レベル非弾性散乱）、冗長なデータ、一貫性の欠如• GNDフォーマット
– XML言語– CSEWG, LLNLが中心に作成– OECD/NEAのWPEC SG38で検討中– 核データ管理パッケージfudge

59

GNDに関する論文

60

XML形式で表現された断面積例
JENDL-4.0, Fe-56の(n,2n)断面積(n,2n)反応
反応の定義 <reaction label="29" outputChannel="n[multiplicity:'2']+ Fe55 + gamma" date="1987-03-01" ENDF_MT="16">断面積の定義 <crossSection nativeData="linear"><linear xData="XYs" length="11" accuracy="0.001"><axes>
内挿法の定義<axis index="0" label="energy_in" unit="eV"
interpolation="linear,linear" frame="lab"/><axis index="1" label="crossSection" unit="b"frame="lab"/></axes><data> 11404700 0 1.17e7 0.0162241 1.2e7 0.0480045 1.3e7
断面積
数値
0.21382 1.4e7 0.389165 1.5e7 0.5134 1.6e7 0.58175 1.7e7 0.61075
データ 1.8e7 0.6118 1.9e7 0.5977 2e70.5759</data></linear></crossSection><outputChannel genre="NBody" Q="-11202700 eV">
二次中性子
<product name="n" label="n" multiplicity="2"
エネルギー
角度分布 ENDFconversionFlag="MF6">の定義<distributions nativeData="Legendre"><Legendre nativeData="LegendrePointwise">

61

NJOY
•LANLで開発•MCNP, PARTISN用の断面積ライブラリ作成
システム

ENDF
RECONR

線形化、再構成

BROADR

ドップラー拡がり

HEATR

核発熱定数

THERMR

熱中性子散乱多群定数作成GROUPR

PURR

確率的非分離処理

MATXS形式出力

MATXSR

ACER

ACE形式出力

多群断面積
ライブラリ

MCNP
ライブラリ

62

NJOY
• NJOY99(99.393)：2013年1月リリース– Fortran77– RSICC, NEA Data Bank, RISTから入手可• 最新版NJOY2012： 2013年2月リリース– Fortran90/95– R Matrix Limitedを処理可– LANLより直接入手, ソースも入手可
Issue Fee

Per Additional User (≥2)

Foreign Academic

$750

N/A

Foreign Government

$750

N/A

Foreign Commercial$2,000 (includes 1 user)

$500

63

PREPRO• D. E. Cullen氏により開発• IAEA核データ部で公開• 最新版PREPRO2012
ENDF/B-7.1が処理可
PREPROシステム

ユーザが作成熱中性子散乱THERM-J
確率テーブル

U3R-J

ENDF
LINEAR

線形化

RECENT

共鳴再構成

HEATR
(NJOY)MVPLIBMKSIGMA1
PENDF

ドップラー

MVP
ライブラリ

64

その他の核データ処理コード• AMPX – ORNL– SCALE用の核データ処理コードシステム• GALILEE – CEA– V1計画で自作の処理コードシステムを開発• CALENDF – Culham Center for Fusion
Energy, UK– J.-C. Sublet氏が開発• TNudy – 成均館大学校, South Korea– Rootシステム上の核データ処理システム

65

国産核データ処理コードシステムの開発• FRENDY計画
FRENDY＝FRom Evaluated Nuclear DatalibrarY (to any code/application)
NJOYの代替となる国産の核データ処理コード
システムの開発。
2013年よりJAEAにおいて開始。
オブジェクト指向、C++言語
今年度は、ACEフォーマットのライブラリを作成できることを目指す。（部分的に既存コードを再
利用）

66

FRENDY Endf6Parserクラスの利用例
JENDL-4.0

67

FRENDYによる共鳴の再構成
1.0E+4

Cross section (barn)

JENDL-4.0, U-235

Total

1.0E+3
1.0E+2
1.0E+1
1.0E+01.0E-2

1.0E+01.0E+2

1.0E+4

Energy (eV)

1.0E+6

68

まとめ• 核データ処理について解説– ENDF-6フォーマット、PENDFの作成、断面積ライ
ブラリ作成• 共分散データ処理（今回割愛）– 原子炉の設計精度評価に有用– ERRORJの開発、ORNLのTSUNAMIシステム• 核データ処理システムの現状– NJOY, PREPRO, etc.
– 国産核データ処理システムFRENDY69

テキスト正誤表• 2.2節、第1段落– (誤）「与えられる異なる」 → (正) 「与えられることになる」• 5節、(37)式– (誤） Cz → Ck （４ヶ所）• 6.3節、第２段落– (誤） σI → σ1 (下付き添え字を”I(アイ)”から”1(数字)”に変える。)• 7節、(66)式上– (誤） 中世子束 → 中性子束

70

関連図書
• ENDF-6 Formats Manual, CSEWG Document ENDF-102, BNL90365-2009 Rev.2 (2011).
• D. E. Cullen, “Nuclear Data Preparation,” Chap. 4, in
Handbook of Nuclear Engineering, Springer (2010).
• R. E. MacFarlane, “The NJOY Nuclear Data Processing
System, Version 2012,” LA-UR-12-27079 (2012).
• R. E. MacFarlane, A. C. Kahler, “Methods for Processing
ENDF/B-VII with NJOY,” Nucl. Data Sheets, 111, 2739-2890(2010).
• “Understanding NJOY,” LA-UR-00-1538 (2000).
http://t2.lanl.gov/nis/njoy/• “Introduction to ENDF format,” LA-UR-98-1779 (1998).
http://t2.lanl.gov/nis/endf/

71

終わり