Initialising ...

Initialising ...

Initialising ...

Initialising ...

Initialising ...

Initialising ...

Initialising ...

Kubo, Kotaro; Zheng, X.; Tanaka, Yoichi; Tamaki, Hitoshi; Sugiyama, Tomoyuki; Jang, S.*; Takata, Takashi*; Yamaguchi, Akira*

Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.308 - 315, 2020/10

Dynamic probabilistic risk assessment (PRA) is a method for improving the realism and completeness of conventional PRA. However, enormous calculation costs are incurred by these improvements. One solution is to select an appropriate sampling method. In this paper, we applied the Monte Carlo, Latin hypercube, grid-point, and quasi-Monte Carlo sampling methods to the dynamic PRA of a simplified accident sequence and compared the results for each method. Quasi-Monte Carlo sampling was found to be the most effective method in this case.

Hashimoto, Shintaro; Sato, Tatsuhiko

Journal of Nuclear Science and Technology, 56(4), p.345 - 354, 2019/04

Times Cited Count：1 Percentile：23.13(Nuclear Science & Technology)Particle transport simulations based on the Monte Carlo method have been applied to shielding calculations. Estimation of not only statistical uncertainty related to the number of trials but also systematic one induced by unclear physical quantities is required to confirm the reliability of calculated results. In this study, we applied a method based on analysis of variance to shielding calculations. We proposed random- and three-condition methods. The first one determines randomly the value of the unclear quantity, while the second one uses only three values: the default value, upper and lower limits. The systematic uncertainty can be estimated adequately by the random-condition method, though it needs the large computational cost. The three-condition method can provide almost the same estimate as the random-condition method when the effect of the variation is monotonic. We found criterion to confirm convergence of the systematic uncertainty as the number of trials increases.

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.

Maeda, Makoto; Komeda, Masao; Tobita, Hiroshi; Ozu, Akira; Kureta, Masatoshi; Bogucarska, T.*; Crochemore, J. M.*; Varasano, G.*; Pedersen, B.*

Dai-37-Kai Kaku Busshitsu Kanri Gakkai Nihon Shibu Nenji Taikai Rombunshu (CD-ROM), 7 Pages, 2017/02

JAEA and EC/JRC are carrying out collaborative research to develop NDA techniques that can be utilized for quantification of high radioactive special nuclear materials such as spent fuel and next generation minor actinide fuels. In the research, reliability of neutron transport codes is important because it is utilized for design and development of a demonstration system of next-generation Differential Die-away (DDA) technique in JAEA. In order to evaluate the reliability, actual neutron flux distribution in a sample cavity was examined in PUNITA device using JRC type DDA technique and JAWAS-T device using JAEA type DDA technique, and then the measurement results were compared with the simulation results obtained by the neutron transport codes. The neutron flux distribution in the target matrix was also examined in the PUNITA and compared with the simulation results. We report on the measurement and simulation results of the neutron flux distribution and evaluation results of the reliability of the neutron transport codes.

Maeda, Makoto; Komeda, Masao; Tobita, Hiroshi; Ozu, Akira; Kureta, Masatoshi; Bogucarska, T.*; Crochemore, J. M.*; Varasano, G.*; Pedersen, B.*

Proceedings of INMM 57th Annual Meeting (Internet), 9 Pages, 2016/07

The JAEA and EC/JRC have started collaborative research to develop a technique that can be utilized for quantification of high radioactive special nuclear materials such as next generation minor actinide fuels. In the study of a Differential Die-Away (DDA) technique, which is one of the techniques to be improved in the collaborative research, JRC type and JAEA type DDA techniques are compared. In the JRC type DDA technique, large amount of thermal neutron is generated using D-T neutron generator and graphite moderator to accomplish high detection sensitivity for small amount of fissile material. On the other hand, in JAEA type, relatively hard neutron spectrum and moderation of neutron in the target matrix are utilized to minimize position dependence of detection efficiency. Estimation of the neutron field is important to evaluate the performance of the system in DDA technique. The purpose of this study is to validate simulation results by experimental results and evaluate neutron flux distribution in the system by the simulation and the experiment. In this paper, we present the evaluation results of the neutron flux distributions in PUNITA which utilizes JRC type DDA technique and JAWAS-T which utilizes JAEA type DDA technique obtained by Monte Carlo simulation and activation method.

Fukushima, Masahiro; Kitamura, Yasunori; Kugo, Teruhiko; Okajima, Shigeaki

Journal of Nuclear Science and Technology, 53(3), p.406 - 424, 2016/03

Times Cited Count：8 Percentile：71.87(Nuclear Science & Technology)Tsutsumi, Masahiro; Tanimura, Yoshihiko

Nuclear Instruments and Methods in Physics Research A, 557(2), p.554 - 560, 2006/02

Times Cited Count：12 Percentile：65.95(Instruments & Instrumentation)A new cerium-doped LaCl(Ce) scintillator is evaluated with respect to the application in environmental -ray dosimetry and spectrometry. The scintillator is very attractive for -ray spectrometry in the case of high count rate, because it has excellent energy resolution and fast decay time. The performance characteristics of a scintillator with a 25.4 mm 25.4 mm LaCl(Ce) crystal are studied and compared to those of a NaI(Tl) scintillator with the same size crystal. Acquired pulse-height spectra are converted to dose rates by using the G(E) function method. Though the LaCl(Ce) crystal itself produces a rather high background in the crystal itself, the scintillator provides good energy information and dose-rate readings from low to high-level (several mGy/h) by subtracting the self-background. The properties of LaCl(Ce) scintillator suggest that the scintillator could be a promising candidate for monitoring at high-dose levels as in emergencies, as well as at ordinary levels of background radiation.

Asano, Yoshihiro; Sugita, Takeshi*; Suzaki,Takenori; Hirose, Hideyuki

Radiation Protection Dosimetry, 116(1-4), p.284 - 289, 2005/12

Times Cited Count：0 Percentile：0.01(Environmental Sciences)no abstracts in English

Kugo, Teruhiko; Mori, Takamasa

Proceedings of International Topical Meeting on Mathematics and Computation, Supercomputing, Reactor Physics and Nuclear and Biological Applications (M&C 2005) (CD-ROM), 10 Pages, 2005/09

A new deterministic transport code based on the method of characteristics (MOC) has been developed for heterogeneous transport calculations in core design of innovative reactors which have complex structures. We have investigated the capability of the MOC code for general geometry with an unstructured geometry PWR problem. The comparison of the results with accurate Monte Carlo calculation results by GMVP has confirmed that the MOC code produces satisfactory results and has a capability to treat unstructured geometry.

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.

Nagaya, Yasunobu; Mori, Takamasa

Journal of Nuclear Science and Technology, 42(5), p.428 - 441, 2005/05

Times Cited Count：50 Percentile：95.28(Nuclear Science & Technology)A new method to estimate a change in the effective multiplication factor due to the perturbed fission source distribution has been proposed for Monte Carlo perturbation calculations with the correlated sampling and differential operator sampling techniques. The method has been implemented into the MVP code for verification. Simple benchmark problems have been set up for fast and thermal systems and the applicability of the method has been verified with the problems. In consequence, it has been confirmed that the method is very effective to estimate the change. It has been also shown that there are some cases where the perturbed source effect is significant and the change in reactivity cannot be estimated accurately without taking the effect into account. Even in such cases, the new method can estimate the perturbed source effect and the estimation of the change in reactivity has been remarkably improved.

Matsubayashi, Masahito; Hibiki, Takashi*; Mishima, Kaichiro*; Yoshii, Koji*; Okamoto, Koji*

Nuclear Instruments and Methods in Physics Research A, 533(3), p.481 - 490, 2004/11

Times Cited Count：3 Percentile：25.71(Instruments & Instrumentation)The validity of a fast neutron radiography quantification method, the -scaling method, which was originally proposed for thermal neutron radiography was examined with Monte Carlo calculations and experiments conducted at the YAYOI fast neutron source reactor. Water and copper were selected as comparative samples for a thermal neutron radiography case and a dense object, respectively. Although different characteristics on effective macroscopic cross-sections were implied by the simulation, the -scaled experimental results with the fission neutron spectrum cross-sections were well fitted to the measurements for both the water and copper samples. This indicates that the -scaling method could be successfully adopted for quantitative measurements in fast neutron radiography.

Mori, Takamasa; Nagaya, Yasunobu; Okumura, Keisuke; Kaneko, Kunio*

JAERI-Data/Code 2004-011, 119 Pages, 2004/07

The 2nd version of code system, LICEM-2, has been developed to produce neutron cross section libraries for the MVP continuous energy Monte Carlo code from an evaluated nuclear data library in the ENDF format. The code system can process nuclear data in the latest ENDF-6 format and produce cross section libraries for MVP's capability of transport calculation at arbitrary temperature. By using the present system, MVP neutron cross section libraries have been prepared from the latest evaluations of JENDL, ENDF/B and JEFF data bases. This report describes the specification of MVP neutron cross section library, the details of each code in the code system, how to use them and MVP neutron cross section libraries produced with the code system.

Sato, Satoshi; Kawasaki, Nobuo*; Kume, Etsuo; Nishitani, Takeo

Journal of Nuclear Science and Technology, 41(Suppl.4), p.62 - 65, 2004/03

no abstracts in English

Oguri, Tomomi*; Nakashima, Hiroshi; Tanaka, Susumu; JAERI-University Collaboration Group for Accelerator Shielding Study

Journal of Nuclear Science and Technology, 41(Suppl.4), p.54 - 57, 2004/03

In shielding design of proton accelerator facilities, it applies simplified formulae and Monte Carlo calculation methods to calculation radiation streaming of access way and duct. A radiation streaming experiment was carried out in order to validate accuracy of intermediate energy region in labyrinth of access way connected to the 1st heavy ion room in TIARA at JAERI. In this paper, it is described comparing of measured data and calculations by MCNP-4B Monte Carlo code and simplified formulae. In Monte Carlo calculations, HILO86 constant set and LA150 was used for cross section data. The results of this experimental analysis show the accuracy of 50% errors for Monte Carlo calculations overall and a factor of 3 for simplified formulae. So applicability for these calculation methods in region of intermediate energy is validated.

Takamatsu, Kuniyoshi; Shimakawa, Satoshi; Nojiri, Naoki; Fujimoto, Nozomu

JAERI-Tech 2003-081, 49 Pages, 2003/10

In the case of evaluations for the highest temperature of the fuels in the HTTR, it is very important to expect the power density distributions accurately; therefore, it is necessary to improve the analytical model with the neutron diffusion and the burn-up theory. The power density distributions are analyzed in terms of two models, the one mixing the fuels and the burnable poisons homogeneously and the other modeling them heterogeneously. Moreover these analytical power density distributions are compared wtih the ones derived from the gross -ray measurements and the Monte Carlo calculational code with continuous energy. As a result the homogeneous mixed model isn't enough to expect the power density distributions of the core in the axial direction; on the other hand, the heterogeneous model improves the accuracy.

Sato, Satoshi

JAERI-Research 2003-014, 223 Pages, 2003/09

no abstracts in English

Aoki, Kosuke*; Izumi, Yoshinobu*; Nishijima, Shigehiro*; Okuno, Kiyoshi; Koizumi, Norikiyo

IEEE Transactions on Applied Superconductivity, 13(2), p.1744 - 1747, 2003/06

Times Cited Count：1 Percentile：12.74(Engineering, Electrical & Electronic)The strand position in CIC conductor consisting of 1,152(34446) strands was analytically evaluated using Monte Carlo method. During the conductor fabrication, the conduit was compressed with the cable to fix the size and shape from one end. This makes the strands to be stretched along the axis. In the calculation, such manufacturing process was simulated. In addition, the contact energy between strands and the strain energy in the strands are considered. The calculation results show that the strands are moved by compressing the conduit. By this calculation, not only strand positions but also the distribution of contact stress between strands could be evaluated.

Kuroishi, Takeshi; Nomura, Yasushi

Journal of Nuclear Science and Technology, 40(6), p.433 - 440, 2003/06

Times Cited Count：1 Percentile：11.39(Nuclear Science & Technology)Effective source acceleration method is studied in criticality safety analysis for realistic spent fuel transport cask. Various axial burnup profiles based on in-core flux measurements are proposed in the OECD/NEA/BUC benchmark Phase II-C. In some cases, calculations by ordinary Monte Carlo method show very slow convergence of fission source distribution, and unacceptably large skipped cycles are needed. The matrix eigenvector calculation that has been developed and incorporated in the ordinary Monte Carlo calculation to improve the slow convergence is applied to the benchmark. The efficiency of this method depends on the precision of matrix elements. In a certain stage of insufficient convergence of fission source distribution, especially for this benchmark of very slow convergence, more acceleration procedure causes anomalous results because of large statistical fluctuations of matrix elements corresponding to low source levels. Therefore, we propose effective source acceleration method with less calculation time than increasing histories for the estimation of matrix elements.

Masaki, Kei; Sugiyama, Kazuyoshi*; Tanabe, Tetsuo*; Goto, Yoshitaka*; Tobita, Kenji; Miyo, Yasuhiko; Kaminaga, Atsushi; Kodama, Kozo; Arai, Takashi; Miya, Naoyuki

Nihon Genshiryoku Gakkai Wabun Rombunshi, 2(2), p.130 - 139, 2003/06

no abstracts in English