Version 1.0 Revised on Nov. 19, 1992 IAEA Benchmark Problem Based on the Time-of-Flight Experiment on Lithium Oxide 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) Files: Data identification is made by file name, for example: LO20122.dat shows lithium-oxude, 20 cm-thick and 12.2 deg. WCT00.DAT ----- Source neutron spectrum LO05***.DAT ------ 50 mm-thick angular neutron flux LO20***.DAT ----- 200 mm-thick angular neutron flux LO40***.dat ----- 400 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 LO05***.DAT 4 0, 24.9, 42.8, 66.8 degrees LO20***.DAT 5 0, 12.2, 24.9, 42.8, 66.8 degrees LO40***.DAT 5 0, 12.2, 24.9, 41.8. 66.8 degrees *In the case of 50 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.8 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 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 ---------------------------------------------- 200 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 ---------------------------------------------- 400 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 Li2O atom density 10**22 atoms/cm**3 : Li-7= 5.338 Li-6= 0.4278 O = 2.883 Fe = 0.1557 Ni = 0.01837 Cr = 0.04223 Mn = 0.003309 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.085. 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 1 Calculational model for MCNP code FNS-TOF/31.4 CM(R)*40.0 CM(Z)-LI2O CYL./FIXED CONE JENDL-3 PT DETECTORS 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.85652-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 -40.0 2 PZ 0 3 CZ 31.4 4 PZ -70 5 CZ 4.9767 6 1 CZ 4.9850 7 2 CZ 5.0121 8 3 CZ 5.0744 9 4 CZ 5.2123 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 -60.00 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 4.600-05 6.959-05 8.399-05 1.727-04 1.554-04 1.932-04 2.512-04 3.390-04 4.520-04 5.189-04 6.487-04 7.450-04 8.417-04 9.862-04 1.215-03 1.421-03 1.591-03 1.810-03 1.933-03 2.203-03 2.353-03 2.397-03 2.468-03 2.501-03 1.255-03 1.238-03 1.243-03 1.162-03 1.164-03 1.193-03 1.185-03 1.448-03 1.718-03 2.186-03 1.653-03 1.040-03 9.251-04 8.761-04 8.385-04 7.505-04 7.832-04 7.368-04 6.683-04 6.248-04 5.473-04 5.120-04 4.928-04 4.124-04 3.956-04 4.339-04 4.881-04 6.196-04 1.714-04 2.017-04 2.109-04 2.134-04 2.369-04 2.547-04 2.576-04 2.685-04 2.849-04 2.881-04 3.445-04 5.177-04 5.260-04 5.417-04 7.421-04 7.535-04 7.649-04 9.229-04 9.597-04 9.782-04 5.512-03 8.224-03 8.337-03 6.554-02 1.400-01 1.424-01 1.440-01 1.456-01 1.480-01 1.454-01 2.054-02 2.095-02 C ********************************************* C * MATERIAL SPECIFICATION CARDS * C ********************************************* C ---- LI2O/SS COVER-------------------------- M1 3006.33C 4.30444-3 3007.33C 5.37069-2 8016.33C 2.90057-2 26000.33C 1.09763-3 28000.33C 1.29482-4 24000.33C 2.97692-4 25055.33C 2.33258-5 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 703.0000 1 148.7974 0 688.2151 1 298.0163 0 642.0212 1 477.4979 0 534.0523 1 676.2308 0 289.8329 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 200000 CTME 8 C *********************************************** C * PERIPHERAL CRADS * C *********************************************** PRDMP 50000 50000 LOST 10 10 PRINT