Sugiyama, Daisuke*; Nakabayashi, Ryo*; Tanaka, Shingo*; Koma, Yoshikazu; Takahatake, Yoko
Journal of Nuclear Science and Technology, 58(4), p.493 - 506, 2021/04
Hashimoto, Shintaro; Sato, Tatsuhiko
EPJ Web of Conferences, 239, p.03015_1 - 03015_4, 2020/09
Particle transport simulation codes based on the Monte Carlo technique have been successfully applied to shielding calculations in accelerator facilities. Estimation of not only statistical uncertainties, which depend on the number of trials, but also systemic uncertainties, which are caused by uncertainty of total cross section models, is required to confirm the reliability of the simulation results. We evaluated unclear quantities of internal parameters included in the total cross section model by the KALMAN code, which is based on the least squares technique, comparing with experimental data of the total cross section. The uncertainties in the total cross sections obtained by the new model are comparable to the experimental errors. In the present study, the systematic uncertainty included in the simulation results can be estimated by performing the transport calculations with variation of the internal parameters within their unclear quantities.
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
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.
Sugiyama, Daisuke*; Nakabayashi, Ryo*; Koma, Yoshikazu; Takahatake, Yoko; Tsukamoto, Masaki*
Journal of Nuclear Science and Technology, 56(9-10), p.881 - 890, 2019/09
Riyana, E. S.; Okumura, Keisuke; Terashima, Kenichi
Proceedings of 27th International Conference on Nuclear Engineering (ICONE-27) (Internet), 4 Pages, 2019/05
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.
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.
Watanabe, Fumitaka; Okuno, Hiroshi
Proceedings of 18th International Symposium on the Packaging and Transport of Radioactive Materials (PATRAM 2016) (DVD-ROM), 9 Pages, 2016/09
This paper shows our calculations on the effects of a radiological release by assuming a severe accident in nuclear material transportation. Following recalculations of safe distance from the point of a severe accident during transportation of a transportation cask TN12 typically used in France containing spent nuclear fuel, and calculations to replicate the "Regulatory Guide: Emergency Preparedness for Nuclear Facilities", a similar calculation was made for a spent fuel transportation cask NFT-14P that was typically utilized in Japan instead of TN12. The safe distance was calculated to be about 30 m. The above calculations were made with the HotSpot codes which adopted the Gauss plume model and had been developed by the USA. Some additional calculations were made with EyesAct, which was developed and used in Japan, adopting also the Gauss plume model, to compare calculation results.
JAEA-Data/Code 2015-015, 162 Pages, 2015/10
MOSRA-SRAC is a lattice calculation module of the Modular code System for nuclear Reactor Analyses (MOSRA). This module performs the neutron transport calculation for various types of fuel elements including existing light water reactors, research reactors, etc. based on the collision probability method with a set of the 200-group cross-sections generated from the Japanese Evaluated Nuclear Data Library JENDL-4.0. It has also a function of the isotope generation and depletion calculation for up to 234 nuclides in each fuel material in the lattice. In these ways, MOSRA-SRAC prepares the burn-up dependent effective microscopic and macroscopic cross-section data to be used in core calculations.
Hayashi, Takao; Tobita, Kenji; Nishio, Satoshi; Ikeda, Kazuki*; Nakamori, Yuko*; Orimo, Shinichi*; DEMO Plant Design Team
Fusion Engineering and Design, 81(8-14), p.1285 - 1290, 2006/02
Neutron transport calculations were carried out to evaluate the capability of metal hydrides and borohydrides as an advanced shielding material. Some hydrides indicated considerably higher hydrogen content than polyethylene and solid hydrogen. The hydrogen-rich hydrides show superior neutron shielding capability to the conventional materials. From the temperature dependence of dissociation pressure, ZrH and TiH can be used without releasing hydrogen at the temperature of less than 640 C at 1 atm. ZrH and Mg(BH) can reduce the thickness of the shield by 30% and 20% compared to a combination of steel and water, respectively. Mixing some hydrides with F82H produces considerable effects in -ray shielding. The neutron and -ray shielding capabilities decrease in order of ZrH Mg(BH) and F82H TiH and F82H water and F82H.
Nagaya, Yasunobu; Okumura, Keisuke; Mori, Takamasa; Nakagawa, Masayuki
JAERI 1348, 388 Pages, 2005/06
To realize fast and accurate Monte Carlo simulation of neutron and photon transport problems, two vectorized Monte Carlo codes MVP and GMVP have been developed at JAERI. MVP is based on the continuous energy model and GMVP is on the multigroup model. Compared with conventional scalar codes, these codes achieve higher computation speed by a factor of 10 or more on vector supercomputers. Both codes have sufficient functions for production use by adopting accurate physics model, geometry description capability and variance reduction techniques. The first version of the codes was released in 1994. They have been extensively improved and new functions have been implemented. The major improvements and new functions are (1) capability to treat the scattering model expressed with File 6 of the ENDF-6 format, (2) time-dependent tallies, (3) reaction rate calculation with the pointwise response function, (4) flexible source specification, etc. This report describes the physical model, geometry description method used in the codes, new functions and how to use them.
Kishimoto, Katsumi; Arigane, Kenji*
JAERI-Tech 2005-016, 83 Pages, 2005/03
Revaluation to activation activity of reactor evaluated at the notification of dismantling submitted in 1997 was carried out in JRR-2 where decommissioning was advanced now. In the revaluation, estimation accuracy on neutron streaming at various horizontal experimental tubes was improved by applying 3 dimensional model to neutron transport calculation that had been carried out by 2 dimensional model, and calculating with TORT. As the result, excessive overestimations on horizontal experimental tubes and biological shield that had greatly contributed to total activation activity in evaluation at the notification of dismantling was revised, sum of their activation activities in the revaluation decreased to 1/18(case after 1 year from the permanent shutdown of reactor) of evaluation at the notification of dismantling, and the structural materials that had large activation activity were changed. By the above, it was shown that introducing 3 dimensional model was effective in evaluation on activation activity of the research reactor that had a lot of various experimental tubes.
Kosako, Kazuaki*; Yamano, Naoki*; Fukahori, Tokio; Shibata, Keiichi; Hasegawa, Akira
JAERI-Data/Code 2003-011, 38 Pages, 2003/07
The third revision of JENDL-3 (JENDL-3.3) was released in May 2002. The library is useful for many applications. For users' convenience, we have produced two JENDL-3.3 based libraries FSXLIB-J33 and MATXSLIB-J33 for transport calculation codes such as MCNP and ANISN. These two libraries are available on request.
JAERI-Research 2002-025, 34 Pages, 2002/11
no abstracts in English
Sato, Satoshi; Iida, Hiromasa; Nishitani, Takeo
Journal of Nuclear Science and Technology, 39(11), p.1237 - 1246, 2002/11
no abstracts in English
Journal of Nuclear Science and Technology, 39(3), p.256 - 263, 2002/03
Two vector computation algorithms; an odd-even sweep (OES) method and an independent sequential sweep (ISS) method, have been developed for the characteristics method to solve the neutron transport equation in a heterogeneous geometry. They realize long vector lengths without recursive operations for effective use of vector computers. Their efficiency has been investigated to a realistic fuel assembly calculation. For both methods, a vector computation is 15 times faster than a scalar computation. From a viewpoint of a comparison between the OES and ISS methods, the ISS method is superior to the OES method because the ISS method shows a faster convergence and saves a computer memory without reducing a computation speed.
Research Committee on Reactor Physics
JAERI-Review 2001-047, 180 Pages, 2002/02
Under the Research Committee on Reactor Physics, the Working Party on Reactor Physics of Accelerator-Driven System (ADS-WP) was set in July 1999 to review and investigate special subjects related to reactor physics research for the Accelerator-Driven Subcritical System (ADS).The ADS-WP, at the first meeting, discussed a task guideline of its activity for two years and decided to concentrate upon three subjects: (1) neutron transport calculations in high energy range, (2) static and kinetic (safety-related) characteristics of subcritical system, and (3) system design including ADS concepts and elemental technology developments required.The activity of ADS-WP continued from July 1999 to March 2001. In this duration, the members of ADS-WP met together four times and discussed the above subjects. In addition, the ADS-WP conducted a questionnaire on requests and proposals for the plan of Transmutation Physics Experimental Facility in the High-Intensity Proton Accelerator Project, which is a joint project between JAERI and KEK (High Energy Accelerator Research Organization).This report summarizes the results obtained by the above ADS-WP activity. The report will be useful to overview those results and moreover to set up a new guideline of future research activity in this field.
JAERI-Research 2001-051, 39 Pages, 2001/11
Fast computation of the characteristics method to solve the neutron transport equation in a heterogeneous geometry has been studied. Two vector computation algorithms; an odd-even sweep (OES) method and an independent sequential sweep (ISS) method have been developed and their efficiency to a typical fuel assembly calculation has been investigated. For both methods, a vector computation is 15 times faster than a scalar computation. From a viewpoint of comparison between the OES and ISS methods, the ISS method is superior to the OES method because the ISS method shows a faster convergence and saves a computer memory without reducing a computation speed. In the vector computation, a table-look-up method to reduce computation time of an exponential function saves only 20% of a whole computation time. Both the coarse mesh rebalance method and the Aitken acceleration method are effective as acceleration methods for the characteristics method, a combination of them saves 70-80% of outer iterations compared with a free iteration.
Sukegawa, Takenori; Hatakeyama, Mutsuo; Yanagihara, Satoshi
JAERI-Tech 2001-058, 81 Pages, 2001/09
In general, neutron transport and activation calculation codes are used for residual radioactive inventory estimation; however, it is essential to verify calculations by measurement results because of geometrical complexity of the reactor and so on. The comparison between measured and calculated radioactivity in the JPDR core components showed a relatively good agreement (factor of 2), and it was cleared that water content and weight ratio of steel bars to concrete materials significantly influenced the neutron flux distribution in the biological shield (factor of 2-10 error). The measured radioactivity inside of the reactor pressure vessel wall and at the inner part of the biological shield was compared in detail with the calculations to verify the methodology applied to calculations of radioisotope production. Then it was found that the radioactive inventory could be estimated accurately with combination of calculations and measurement of radioactivity in samples and dose rate distribution for planning of dismantling activities.