Yokoyama, Kenji; Lahaye, S.*
Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.109 - 116, 2020/10
CEA/DEN/DM2S/SERMA and JAEA/NSEC are working on benchmarks for burnup, isotopic concentrations and decay heat calculations in the collaboration framework between both organisms. Both actors of this benchmark are independently developing their own simulation code systems for computing quantities of interest in nuclear fuel cycle domain: MENDEL in CEA and MARBLE in JAEA. The purpose of the benchmark is to verify each system by comparing both calculation results on specific applications. MENDEL uses a several solvers for the resolution of Bateman equation. Runge-Kutta method or Chebyshev Rational Approximation method (CRAM) are used for irradiation computations. An analytical solver can also be used for decay calculations. MARBLE can use Krylov subspace method or CRAM method. As the first phase of the benchmark, we compared the calculated results of decay heat and isotropic concentrations following by a Pu-239 fast fission pulse. We applied nuclear data from three libraries: (1) JEFF-3.1.1, (2) JENDL/DDF-2015 + JENDL/FPY-2011, and (3) ENDF/B-VII.1. Nuclear data and burnup chain were generated from these libraries independently on each system. We confirmed that the results for both systems were in very good agreement with each other. Numerical results were also compared to experimental data. As the second phase of the benchmark, we are proceeding with a burnup calculation benchmark of MENDEL and MARBLE using the nuclear data and burnup chain provided by ORLIBJ33, which is a set of cross-section data based on JENDL-3.3 for ORIGEN-2 code system. We will also compare with calculation results by the ORIGEN-2 code with ORLIBJ33. Since the series of ORLIB, that is, ORLIBJ32, ORLIBJ33, and ORLIBJ40, have been widely used especially in Japan for many years, the comparison with ORLIB is effective for confirming the performance of MENDEL and MARBLE.
Ishitsuka, Etsuo; Nakashima, Koki*; Nakagawa, Naoki*; Ho, H. Q.; Ishii, Toshiaki; Hamamoto, Shimpei; Takamatsu, Kuniyoshi; Kenzhina, I.*; Chikhray, Y.*; Matsuura, Hideaki*; et al.
JAEA-Technology 2020-008, 16 Pages, 2020/08
As a summer holiday practical training 2019, the feasibility study for nuclear design of a nuclear battery using HTTR core was carried out, and the U enrichment and burnable poison of the fuel, which enables continuous operation for 30 years with thermal power of 5 MW, were studied by the MVP-BURN. As a result, it is clear that a fuel with U enrichment of 12%, radius of burnable poison and natural boron concentration of 1.5 cm and 2wt% are required. As a next step, the downsizing of core will be studied.
Harada, Hideo; Takayama, Naoki; Komeda, Masao
Journal of Physics Communications (Internet), 4(8), p.085004_1 - 085004_17, 2020/08
A new convention of epithermal neutron spectrum is formulated for improving accuracy of resonance integrals. The new type function is proposed as an approximating function of epithermal neutron spectrum based on calculations by the state-of-art Monte Carlo code MVP-3. Bias effects on determination of resonance integrals due to utilizing approximating functions of the traditional types and the new type are compared. The other bias effect is also investigated, which is caused by neglecting position dependence of a neutron spectrum inside an irradiation capsule. For demonstrating the bias effects due to these assumptions on neutron spectrum quantitatively in a practical case, the thermal neutron-capture cross section and resonance integral of Cs measured at a research reactor JRR-3 are re-evaluated. A superior property of the proposed new convention is discussed. The experimental method is proposed to determine the new shape factor introduced in the convention by a combinational use of triple flux monitors (Au, Co and Zr), and its analytical methodology is formulated.
Ikeuchi, Hirotomo; Yano, Kimihiko; Washiya, Tadahiro
Journal of Nuclear Science and Technology, 57(6), p.704 - 718, 2020/06
To suggest efficient process of the fuel debris treatment after the retrieval from the Fukushima Daiichi Nuclear Power Plant (1F), thorough investigation is indispensable on potential source of U in the fuel debris. Estimation on the fuel debris accumulated in the reactor pressure vessel is specifically important due to its limited accessibility. The present study aims to estimate the chemical forms of U in the in-vessel fuel debris, especially in the minor phases such as metallic phases, by performing the thermodynamic calculation considering the material relocation and changing environment during the accident progression in the 1F Unit 2. Input conditions for the thermodynamic calculation such as composition, temperature, and oxygen amount were assumed mainly based on the results of severe accident analysis. The chemical form of U varied depending on the local amount of Fe and O. In regions of low steel content, the U-containing metallic phase was dominated by -(Zr,U)(O), while regions of high steel content were dominated by Fe(Zr,U) (Laves phase). A few percent of U was transferred to the metallic phases under reducing conditions, raising challenging issues on the chemical removal of nuclear material from fuel debris.
Riyana, E. S.; Okumura, Keisuke; Terashima, Kenichi; Matsumura, Taichi; Sakamoto, Masahiro
Mechanical Engineering Journal (Internet), 7(3), p.19-00543_1 - 19-00543_8, 2020/06
Ono, Masato; Fujiwara, Yusuke; Matsumoto, Tetsuro*; Iigaki, Kazuhiko
Nippon Genshiryoku Gakkai Wabun Rombunshi, 19(2), p.110 - 120, 2020/06
Integrity confirmation for buildings against collisions of projectiles has been conducted to evaluate collisions between a projectile with simple shape and a wall using empirical formulas. It is a matter of fact, there is a possibility that structures with complex shape such as stack may collide with a reactor building. However, there were not so many studies of collisions between structures with complex shape and buildings in the literature. Impact evaluation was carried out using reactor building and stack with real shape and adequate physical property. It was found that ceiling of reactor building was not damaged by the collision, confirming that there was no effect inside of reactor building.
Nagasumi, Satoru; Matsunaka, Kazuaki*; Fujimoto, Nozomu*; Ishii, Toshiaki; Ishitsuka, Etsuo
JAEA-Technology 2020-003, 13 Pages, 2020/05
The influence of the control rod model on the nuclear characteristics of the HTTR has been evaluated, by creating detailed control rod model, in which geometric shape was close to that of the actual control rod structure, in MVP code. According to refinement of the control rod model, the critical control rod position was 11 mm lower than that of the conventional model, and this was close to the measured value of 1775 mm. The reactivity absorbed by the shock absorber located at the tip of the control rod was 0.2%k/k, and this was 14 mm difference at the critical control rod position. Considering the effect of refinement of the control rod and the effect of the shock absorber, the correction amount for the analysis value in SRAC code due to the shape effect of the control rod, is -0.05%k/k in reactivity, and -3 mm in the critical control rod position at low temperature criticality.
Sugino, Kazuteru; Takino, Kazuo
JAEA-Data/Code 2019-011, 110 Pages, 2020/01
A deterministic discrete ordinates method (SN method) transport calculation code for three-dimensional hexagonal geometry has been developed as the MINISTRI code (Ver. 7.0). MINISTRI is based on the triangle-mesh finite difference method, which can perform neutron transport calculations with high accuracy for cores of fast power reactors and assemblies of the Russian BFS critical facility. The present study has derived a proper scheme for remarkably improving the convergence of MINISTRI by investigating the issue of previous MINISTRI (Ver. 1.1), which sometimes plays a poor convergence performance in calculations for large-scale power reactor cores. The verification test of improved MINISTRI has been carried out for various cores by setting the reference result as the multi-group Monte-Carlo calculation with the same cross-sections as used in MINISTRI. As a result, it is found that the agreements are within 0.1% for eigenvalues and within 0.7% for power distributions. Thus, the satisfying accuracy of MINISTRI has been confirmed. In order to reduce the calculation time, the initial diffusion calculation scheme and the parallel processing have been implemented. As a result, the calculation time is reduced to the approximately one tenth compared with previous MINISTRI. Furthermore, adoption of the treatment of the anisotropic cell streaming effect, preparation of the perturbation calculation tool, implementation of the function for specification of the triangle-mesh-wise material and merging of the hexagonal-mesh calculation code MINIHEX have been carried out. Thus, the versatility of MINISTRI has been enhanced.
Riyana, E. S.; Okumura, Keisuke; Terashima, Kenichi
Proceedings of 27th International Conference on Nuclear Engineering (ICONE-27) (Internet), 4 Pages, 2019/05
JAEA-Data/Code 2018-018, 103 Pages, 2019/03
The latest available thermodynamic data were critically reviewed and the selected values were included into the JAEA-TDB for performance assessment of geological disposal of high-level radioactive and TRU wastes. This critical review specifically addressed thermodynamic data for (1) a zirconium-hydroxide system through comparison of thermodynamic data selected by the Nuclear Energy Agency within the Organisation for Economic Co-operation and Development (OECD/NEA), (2) complexation of metal ions with isosaccharinic acid based on the latest review papers. Furthermore, the author performed (3) tentative selection of thermodynamic data on ternary complexes among alkaline-earth metal, uranyl and carbonate ions, and (4) integration with the latest version of JAEA's thermodynamic database for geochemical calculations. The internal consistency of the selected data was checked by the author. Text files of the updated and integrated thermodynamic database have been prepared for geochemical calculation programs of PHREEQC and Geochemist's Workbench.
Miyahara, Naoya; Miwa, Shuhei; Horiguchi, Naoki; Sato, Isamu*; Osaka, Masahiko
Journal of Nuclear Science and Technology, 56(2), p.228 - 240, 2019/02
In order to improve LWR source term under severe accident conditions, the first version of a fission product (FP) chemistry database named "ECUME" was developed. The ECUME is intended to include major chemical reactions and their effective kinetic constants for representative SA sequences. It is expected that the ECUME can serve as a fundamental basis from which FP chemical models in the SA analysis codes can be elaborated. The implemented chemical reactions in the first version were those for representative gas species in Cs-I-B-Mo-O-H system. The chemical reaction kinetic constants were evaluated from either literature data or calculated values using ab-initio calculations. The sample chemical reaction calculation using the presently constructed dataset showed meaningful kinetics effects at 1000 K. Comparison of the chemical equilibrium compositions by using the dataset with those by chemical equilibrium calculations has shown rather good consistency for the representative Cs-I-B-Mo-O-H species. From these results, it was concluded that the present dataset should be useful to evaluate FP chemistry in Cs-I-B-Mo-O-H system under LWA SA conditions.
Ho, H. Q.; Honda, Yuki; Goto, Minoru; Takada, Shoji
Annals of Nuclear Energy, 112, p.42 - 47, 2018/02
Sugawara, Takanori; Katano, Ryota; Tsujimoto, Kazufumi
Annals of Nuclear Energy, 111, p.449 - 459, 2018/01
This study aims to review the ADS design based on the outcome for the last dozen years and to investigate the impact of impurities in the transmutation cycle on the ADS neutronics design. The impact of impurities in the transmutation cycle is investigated for the reviewed reference design. For the uranium from the partitioning, the accompaniment of 20 wt.% U against the Pu weight is acceptable although the MA transmutation amount will be decreased slightly. For the rare earth (RE) from the partitioning, the accompaniment of 5 wt.% RE against the MA weight is allowable. In the reprocessing, the decontamination factor, DF=10 for RE is enough from the viewpoint of the neutronics design. The impact of the fuel composition accuracy is also investigated. The uncertainty of the ZrN ratio against the MA fuel should be less than 0.2% to minimize a surplus proton beam current due to the uncertainty.
Di Lemma, F. G.; Yamashita, Shinichiro; Miwa, Shuhei; Nakajima, Kunihisa; Osaka, Masahiko
Energy Procedia, 127, p.29 - 34, 2017/09
Chemical effects of molybdenum (Mo) and boron (B), which were considered to form compounds with Cs, on the Cs chemisorption were predicted using a chemical equilibrium calculation. It is seen that CsMoO were formed in the chemisorbed compounds. On the other hand, little effects were observed for B. The results suggest that the effects of Mo should be considered for further experimental investigation.
Ikeuchi, Hirotomo; Piluso, P.*; Fouquart, P.*; Excoffier, E.*; David, C.*; Brackx, E.*
Proceedings of 8th European Review Meeting on Severe Accident Research (ERMSAR 2017) (Internet), 12 Pages, 2017/05
no abstracts in English
Ohgama, Kazuya; Ikeda, Kazumi*; Ishikawa, Makoto; Kan, Taro*; Maruyama, Shuhei; Yokoyama, Kenji; Sugino, Kazuteru; Nagaya, Yasunobu; Oki, Shigeo
Proceedings of 2017 International Congress on Advances in Nuclear Power Plants (ICAPP 2017) (CD-ROM), 10 Pages, 2017/04
Nagaya, Yasunobu; Okumura, Keisuke; Sakurai, Takeshi; Mori, Takamasa
JAEA-Data/Code 2016-018, 421 Pages, 2017/03
In order to realize fast and accurate Monte Carlo simulation of neutron and photon transport problems, two Monte Carlo codes MVP (continuous-energy method) and GMVP (multigroup method) have been developed at Japan Atomic Energy Agency. The codes have adopted a vectorized algorithm and have been developed for vector-type supercomputers. They also support parallel processing with a standard parallelization library MPI and thus a speed-up of Monte Carlo calculations can be achieved on general computing platforms. The first and second versions of the codes were released in 1994 and 2005, respectively. They have been extensively improved and new capabilities have been implemented. The major improvements and new capabilities are as follows: (1) perturbation calculation for effective multiplication factor, (2) exact resonant elastic scattering model, (3) calculation of reactor kinetics parameters, (4) photo-nuclear model, (5) simulation of delayed neutrons, (6) generation of group constants, etc. This report describes the physical model, geometry description method used in the codes, new capabilities and input instructions.
Inaba, Yoshitomo; Nishihara, Tetsuo
Annals of Nuclear Energy, 101, p.383 - 389, 2017/03
In order to ensure the thermal integrity of fuel in High Temperature Gas-cooled Reactors (HTGRs), it is necessary that the maximum fuel temperature in normal operation is to be lower than a thermal design target. In the core thermal-hydraulic design of block-type HTGRs, the maximum fuel temperature should be evaluated considering data such as thermal power, core geometry, power density and neutron fluence distributions, and core coolant flow distribution. The fuel temperature calculation code used in the design stage of the High Temperature engineering Test Reactor (HTTR) presupposes to run on UNIX systems, and its operation and execution procedure are complicated and are not user-friendly. Therefore, a new fuel temperature calculation code named FTCC which has a user-friendly system such as a simple and easy operation and execution procedure, was developed. This paper describes calculation objects and models, basic equations, improvement points from the HTTR design code in FTCC, and the result of a validation calculation with FTCC. The calculation result obtained by FTCC provides good agreement with that of the HTTR design code, and then FTCC will be used as one of the design codes for HTGRs. In addition, the effect of cooling forms on the maximum fuel temperature is investigated by using FTCC. As a result, it was found that the effect of center hole cooling for hollow fuel compacts and gapless cooling with monolithic type fuel rods on reducing the temperature is very high.
Shinohara, Masanori; Ishitsuka, Etsuo; Shimazaki, Yosuke; Sawahata, Hiroaki
JAEA-Technology 2016-033, 65 Pages, 2017/01
To reduce the neutron exposure dose for workers during the replacement works of the startup neutron sources of the High Temperature Engineering Test Reactor, calculations of the exposure dose in case of temporary neutron shielding at the bottom of fuels handling machine were carried out by the PHITS code. As a result, it is clear that the dose equivalent rate due to neutron radiation can be reduced to about an order of magnitude by setting a temporary neutron shielding at the bottom of shielding cask for the fuel handling machine. In the actual replacement works, by setting temporary neutron shielding, it was achieved that the cumulative equivalent dose of the workers was reduced to 0.3 man mSv which is less than half of cumulative equivalent dose for the previous replacement works; 0.7 man mSv.