Version 1.0 Revised on Nov. 19, 1992 IAEA Benchmark Problem Based on the Time-of-Flight Experiment on Graphite Slabs at FNS/JAERI Yukio OYAMA and Hiroshi MAEKAWA Fusion Reactor Physics Laboratory Department of Reactor Engineering Japan Atomic Energy Research Institute Tokai-mura, Naka-gun, Ibaraki-ken, 319-11 Japan (Phone) 81-292-82-6075 (for Y.O.) or -6015 (for H.M.) (Telefax) 81-292-82-5709 (e-mail) oyama@fnshp.tokai.jaeri.go.jp (for Y.O.) maekawa@fnshp.tokai.jaeri.go.jp (for H.M.) 1. Numerical Data and Their Format a) Datafiles: Data identification is made by file name, for example: C20122.dat shows carbon, 20 cm-thick and 12.2 deg. ACT00.DAT ----- Source neutron spectrum C05***.DAT ------ 50.6mm-thick angular neutron flux C20***.DAT ----- 202.4 mm-thick angular neutron flux C40***.DAT ----- 404.8 mm-thick angular neutron flux ENERGY.DAT ---------- Boundary energy b) Units of the data Source spectrum [n/sr/lethargy/source] Angular flux [n/sr/cm2/lethargy/source] c) Data format for experimental data Comment 1 line 20A4 Energy (mid-point) [MeV] I = 1, 150 6E12.4 Angular flux [See above] I = 1, 150 6E12.4 Error (fraction) [%] I = 1, 150 6E12.4 d) Number of data set wct00.DAT 1 0 degree C05***.DAT 4 0, 24.9, 42.8, 66.8 degrees C20***.DAT 5 0, 12.2, 24.9, 42.8, 66.8 degrees C40***.DAT 5 0, 12.2, 24.9, 42.8, 66.8 degrees *In the case of 50.6 mm-thick and 12.2 degree, it was revealed that the measured flux was contaminated by some of direct neutrons from the target and too difficult to analysis. Therefore the data of 50.6 mm-thick and 12.2 degree is not included in the file. e) Data format for boundary energy (ENERGY.DAT) Comment 1 line 20A4 Boundary energy [MeV] I = 1, 151 8F9.5 2. Flight Path and Effective Measured Area The flight path and effective measured area are summarized as follows. The meaning of them is described in the references. These data are useful for Monte Carlo calculations. Table Flight path and measured area 50.6 mm-thick Assembly ---------------------------------------------- Angle Flight Path [cm] Measured Area [cm**2] ---------------------------------------------- 0.0 738 85.88 24.9 740 86.34 41.8 744 87.26 66.8 753 89.33 ---------------------------------------------- 202.4 mm-thick Assembly ---------------------------------------------- Angle Flight Path [cm] Measured Area [cm**2] ---------------------------------------------- 0.0 723 82.42 12.2 724 82.65 24.9 726 83.11 41.8 732 84.49 66.8 746 87.72 ---------------------------------------------- 404.8 mm-thick Assembly ---------------------------------------------- Angle Flight Path [cm] Measured Area [cm**2] ---------------------------------------------- 0.0 703 77.81 12.2 704 78.04 24.9 708 78.96 41.8 716 80.81 66.8 736 85.41 ---------------------------------------------- 3. Calculational Model The calculational model is of course depend on the code used. Basic data are as follows: Radius of assembly : 31.5 cm Thickness of assembly : 5, 20 and 40 cm Distance between target and assembly : 20 cm C atom density 10**22 atoms/cm**3 : 8.280 Area of detector : (pai)x2.54**2 [cm**2] 4. How to Calculate the Angular Flux (1) Neutron source You can use the data of "WCT00.DAT" as the source. Before you start the calculation you should interpolate the spectrum to adjust it to the group structure used. It is notable that the integrated source spectrum multiplied by 4x(pai) is not unity (1.0) but about 1.055. The source neutrons are generated isotropically at the target position. (2) Two-demensional discrete ordinate code such as DOT3.5 You can easily understand from the reference 4) how to do. It is important to average over the measured area when the calculated angular flux is compared with the measured one. (3) Monte Carlo calculation You can see a sample of method in the reference 1). 5. Comparison Two types of comparison will be done for this benchmark problem. (1) To compare the measured and calculated angular fluxes directly in graph. (2) To compare the integrated flux over energy regions in C/E values (Ratio of calculated to experimental values); 6. References When you publish the results using this numerical data, you should refer at least following references 1) and 2). 1) Oyama Y., Maekawa H.: "Measurement and Analysis of an Angular Neutron Flux on a Beryllium Slab Irradiated with Deuteron-Tritium Neutrons," Nucl. Sci. Eng., 97, 220-234 (1987). 2) Oyama Y., Yamaguchi S., Maekawa H.: "Experimental Results of Angular Neutron Spectra Leaking from Slabs of Fusion Reactor Candidate Materials (I)," JAERI-M 90-092 (1990). 3) Oyama Y., Maekawa H.: "Measurements of Angle-Dependent Neutron Spectra from Lithium-Oxide Slab Assemblies by Time-of-Flight Method," JAERI-M 83-195 (Nov. 1983). 4) Oyama Y., Yamaguchi S., Maekawa H.: "Analysis of Time-of-Flight Experiment on Lithium-Oxide Assemblies by a Two-Dimensional Transport Code DOT3.5," JAERI-M 85-031 (March 1985). 5) Oyama Y., Maekawa H.: "Spectral Measurement of Angular Neutron Flux on the Restricted Surface of Slab Assemblies by the Time-of-Flight Method," Nucl. Instr. Methods, A245 173-181 (1986). 6) Oyama Y., Yamaguchi S., Maekawa H.: "Measurements and Analyses of Angular Neutron Flux Spectra on Graphite and Lithium-Oxide Slabs Irradiated with 14.8 MeV Neutrons," J. Nucl. Sci. Technol., 25, 419-428 (1988). 7) Oyama Y.: "Experimental Study of Angular Neutron Flux Spectra on a Slab Surface to Assess Nuclear Data and Calculational Methods for a Fusion Reactor Design," JAERI-M 88-101 (June 1988). Attachment MCNP INPUT sample FNS-TOF/31.4 CM(R)*20.24CM(Z)-C CYL./FIXED CONE BIAS/JENDL-3 C ************************************** C * CELL CARAD * C ************************************** C ******* EXTERNAL VOID **************** 1 0 -4 : +4 -2 +3 : +10 : +2 -10 +5 +6 +7 +8 +9 C ******* SOURCE VACUUM REGION ********* 2 0 -3 +4 -1 C ******* MATERIAL REGION ************** 3 1 8.694-2 -3 +1 -2 C ******* DETECTOR VACUUM REGION ******* 4 0 +2 -10 -5 : +2 -10 -6 : +2 -10 -7 : +2 -10 -8 : +2 -10 -9 C ------- THE FOLLOWING IS A BLANK DELIMETER C ************************************** C * SURFACE CARD * C ************************************** 1 PZ -20.24 2 PZ 0 3 CZ 31.4 4 PZ -60 5 CZ 5.122 6 1 CZ 5.129 7 2 CZ 5.143 8 3 CZ 5.186 9 4 CZ 5.284 10 SO 1000 C ------- THE FOLLOWING IS A BLANK DELIMETER C ************************************** C * MODE CARD - NEUTRON ONLY * C ************************************** MODE 0 C ************************************** C * TRANSFORMATION CARDS * C * ROTATION ABOUT THE Y AXIS BY THETA* C ************************************** *TR1 0 0 0 12.2 90 102.2 90 0 90 77.8 90 12.2 +1 *TR2 0 0 0 24.9 90 114.9 90 0 90 65.1 90 24.9 +1 *TR3 0 0 0 41.8 90 131.8 90 0 90 48.2 90 41.8 +1 *TR4 0 0 0 66.8 90 156.8 90 0 90 23.2 90 66.8 +1 C ************************************** C * CELL PARAMETER CARDS * C ************************************** IN 0 1 1 1 C **************************************************** C * SOURCE SPECIFICATION CARDS * C * SRC1=POINT ISOTROPIC OPTION * C * SDIR DIRC. BIASING - HEIGHT REDUCTION CONSIDERED* C * SI(ENG.) AND SP(PROB.) TAKEN FROM BETOF SOURCE * C * EXPT. DATA * C **************************************************** SRC1 0 0 -40.24 2 1.0 SDIR 0 0 1 1 0.5372231 SI 4.6308-02 5.2474-02 5.9461-02 6.7378-02 7.6349-02 8.6515-02 9.8035-02 1.1109-01 1.2588-01 1.4264-01 1.6163-01 1.8315-01 2.0754-01 2.3517-01 2.6649-01 3.0197-01 3.4217-01 3.8774-01 4.3936-01 4.9786-01 5.6415-01 6.3927-01 7.2438-01 8.2084-01 9.3013-01 1.0540+00 1.1943+00 1.3533+00 1.5335+00 1.7377+00 1.8498+00 1.9691+00 2.0961+00 2.2313+00 2.3752+00 2.5284+00 2.6914+00 2.8650+00 3.0498+00 3.2465+00 3.4559+00 3.6787+00 3.9160+00 4.1686+00 4.4374+00 4.7236+00 5.0282+00 5.3525+00 5.6978+00 6.0652+00 6.4564+00 6.8728+00 7.3161+00 7.7879+00 8.2902+00 8.8249+00 9.3940+00 9.9999+00 1.0157+01 1.0317+01 1.0480+01 1.0645+01 1.0812+01 1.0983+01 1.1156+01 1.1331+01 1.1510+01 1.1691+01 1.1875+01 1.2062+01 1.2252+01 1.2445+01 1.2641+01 1.2840+01 1.3042+01 1.3248+01 1.3456+01 1.3668+01 1.3883+01 1.4102+01 1.4324+01 1.4550+01 1.4779+01 1.5012+01 1.5248+01 1.5488+01 1.5732+01 1.5980+01 1.6231+01 1.6487+01 SP 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2.1438-04 9.1563-04 1.5388-03 2.1627-03 2.5871-03 2.9698-03 3.3277-03 3.7304-03 3.9914-03 4.1831-03 4.3067-03 2.2786-03 2.1176-03 2.1039-03 2.0513-03 2.1117-03 2.0411-03 2.0591-03 2.1096-03 2.1281-03 2.4053-03 2.0935-03 1.8313-03 1.8008-03 1.7868-03 1.6668-03 1.6064-03 1.5991-03 1.6805-03 1.6701-03 1.5432-03 1.2904-03 1.2595-03 1.0310-03 9.2870-04 8.9403-04 8.8635-04 8.7646-04 9.5094-04 2.2846-04 2.2733-04 2.9803-04 3.5343-04 4.2411-04 4.5049-04 4.8501-04 5.7358-04 6.6128-04 7.2272-04 1.0292-03 1.1708-03 1.3518-03 1.5306-03 1.6062-03 1.6985-03 1.7860-03 1.8887-03 2.6840-03 4.0301-03 5.7679-03 1.2112-02 2.2014-02 5.0048-02 1.4342-01 1.9945-01 2.1655-01 1.6617-01 9.8846-02 2.5913-02 9.2423-03 1.2696-03 C ********************************************* C * MATERIAL SPECIFICATION CARDS * C ********************************************* C ---- CARBON ------------------------------- M1 6012.33C 1.0 C DRXS C ********************************************** C * TALLY SPECIFICATION CARDS * C ********************************************** FC5 --- FLUXES AT 5 PT DTS(TH=0.0, 12.2, 24.9, 41.8, 66.8 DEG) F5 0.0 0 723.0000 1 152.999 0 707.649 1 305.672 0 658.514 1 487.901 0 545.688 1 685.675 0 293.881 1 DD 0.5 100 E0 4.6308-02 5.2474-02 5.9461-02 6.7378-02 7.6349-02 8.6515-02 9.8035-02 1.1109-01 1.2588-01 1.4264-01 1.6163-01 1.8315-01 2.0754-01 2.3517-01 2.6649-01 3.0197-01 3.4217-01 3.8774-01 4.3936-01 4.9786-01 5.6415-01 6.3927-01 7.2438-01 8.2084-01 9.3013-01 1.0540+00 1.1943+00 1.3533+00 1.5335+00 1.7377+00 1.8498+00 1.9691+00 2.0961+00 2.2313+00 2.3752+00 2.5284+00 2.6914+00 2.8650+00 3.0498+00 3.2465+00 3.4559+00 3.6787+00 3.9160+00 4.1686+00 4.4374+00 4.7236+00 5.0282+00 5.3525+00 5.6978+00 6.0652+00 6.4564+00 6.8728+00 7.3161+00 7.7879+00 8.2902+00 8.8249+00 9.3940+00 9.9999+00 1.0157+01 1.0317+01 1.0480+01 1.0645+01 1.0812+01 1.0983+01 1.1156+01 1.1331+01 1.1510+01 1.1691+01 1.1875+01 1.2062+01 1.2252+01 1.2445+01 1.2641+01 1.2840+01 1.3042+01 1.3248+01 1.3456+01 1.3668+01 1.3883+01 1.4102+01 1.4324+01 1.4550+01 1.4779+01 1.5012+01 1.5248+01 1.5488+01 1.5732+01 1.5980+01 1.6231+01 1.6487+01 FQ0 E F C ********************************************** C * ENERGY AND THERMAL CARDS * C ********************************************** ERGN 0 16.5 0 C *********************************************** C * PROBLEM CUTOFF CARDS * C *********************************************** CUTN 0 4.6308-02 -10 -0.01 NPS 100000 CTME 10 C *********************************************** C * PERIPHERAL CRADS * C *********************************************** PRDMP 50000 50000 LOST 10 10 PRINT