Maeyama, Shinya*; Watanabe, Tomohiko*; Idomura, Yasuhiro; Nakata, Motoki*; Nunami, Masanori*
Computer Physics Communications, 235, p.9 - 15, 2019/02
We have implemented the Sugama collision operator in the gyrokinetic Vlasov simulation code, GKV, with an implicit time-integration scheme. The new method is versatile and independent of the details of the linearized collision operator, by means of an operator splitting, an implicit time integrator, and an iterative Krylov subspace solver. Numerical tests demonstrate stable computation over the time step size restricted by the collision term. An efficient implementation for parallel computation on distributed memory systems is realized by using the data transpose communication, which makes the iterative solver free from inter-node communications during iteration. Consequently, the present approach achieves enhancement of computational efficiency and reduction of computational time to solution simultaneously, and significantly accelerates the total performance of the application.
Kitamura, Yasunori*; Fukushima, Masahiro; Kitamura, Yasunori*
Annals of Nuclear Energy, 125, p.328 - 341, 2019/01
It has been taken for granted that the neutron correlation methods that employ two sets of neutron counting systems, e.g., the covariance-to-mean and the cross-correlation methods, are free from the count-loss effect for determination of the neutron decay constant. It was however found in the present study that these methods overestimate the neutron decay constant under high counting rate conditions. New formulae of these methods were hence obtained on the basis of a rigorous theoretical approach for treating the count-loss process. It is expected that the present formulae work better than conventional ones for determination of the neutron decay constant.
Yamaguchi, Masatake; Ebihara, Kenichi; Itakura, Mitsuhiro; Tsuru, Tomohito; Matsuda, Kenji*; Toda, Hiroyuki*
Computational Materials Science, 156, p.368 - 375, 2019/01
The segregation of multiple hydrogen atoms along aluminum (Al) grain boundaries (GBs) and fracture surfaces (FSs) was investigated through first-principles calculations considering the characteristics of GBs. The results indicate that hydrogen segregation is difficult along low-energy GBs. The segregation energy of multiple hydrogen atoms along GBs and FSs and the cohesive energy was obtained for three types of high-energy Al GBs. With increasing hydrogen segregation along the GBs, the cohesive energy of the GB decreases and approaches zero with no decrease in GB segregation energy. The GB cohesive energy decreases in parallel with the volume expansion of the region of low electron density along the GB.
Manabe, Kentaro; Matsumoto, Masaki*
Journal of Nuclear Science and Technology, 56(1), p.78 - 86, 2019/01
If an insoluble cesium-bearing particle is incorporated into the human body, the radioactivity will move as a single particle. In this case, it is impossible to estimate the number of disintegrations by considering the average behavior of countless nuclei. Then, a method was developed to simulate the behavior of the particle stochastically; and a biokinetic model was constructed to consider the characteristics of insoluble particles. Combination of the method and the model enables to estimate the number of disintegrations, and consequently the internal doses considering the stochastic behavior of the single cesium particle. We evaluated a probability density function of committed equivalent and effective doses and its 99th percentile value and arithmetic mean by repeating the above described procedure, and compared them to the reference values based on the existing models. As a result, the 99th percentile value of committed effective doses was 70 times the reference value when the number of incorporated particles was one, and consequently the dose level was quite low. When the exposure level is 1 mSv in committed effective dose, the uncertainty originating in the insolubility of cesium particles was negligibly small.
Kaneko, Masashi; Watanabe, Masayuki; Miyashita, Sunao*; Nakashima, Satoru*
Hyperfine Interactions, 239(1), p.20_1 - 20_10, 2018/12
We optimized a mixing ratio of exchange energy between pure DFT and exact Hartree-Fock using TPSS exchange-correlation functional to estimate the accurate coordination bonds in f-block complexes by numerically benchmarking with the experimental data of Mssbauer isomer shifts for Eu and Np nuclides. Second-order Douglas-Kroll-Hess Hamiltonian with segmented all-electron relativistically contracted basis set was employed to calculate the electron densities at Eu and Np nuclei, i.e. contact densities, for each five complexes for Eu(III) and Np(IV) systems. We compared the root mean square deviation values of their isomer shifts between experiment and calculation by changing the mixing ratio of Hartree-Fock exchange parameter from 0 to 100 % at intervals of 10 %. As the result, it was indicated that the mixing ratio of 30 and 60 % for Eu and Np benchmark systems, respectively, gives the smallest deviation values. Mulliken's spin population analysis indicated that the covalency in the metal-ligand bonds for both Eu and Np complexes decreases with increasing the Hartree-Fock exchange admixture.
Isotope News, (760), p.2 - 5, 2018/12
Recently, we proposed a new model for estimating the biological effectiveness for Boron Neutron Capture Therapy (BNCT) based on the absorbed dose distributions in cellular scale. The model quantitatively highlighted the indispensable need to consider the synergetic effect and the dose dependence of the biological effectiveness in the estimate of the therapeutic effect of BNCT. This paper reviews the basic features of the model.
Fukumitsu, Nobuyoshi*; Yamauchi, Yusuke*; Saptiama, I.*; Ariga, Katsuhiko*; Hatano, Kentaro*; Kumada, Hiroaki*; Fujita, Yoshitaka; Tsuchiya, Kunihiko
Isotope News, (760), p.15 - 18, 2018/12
Otosaka, Shigeyoshi; Sato, Yuhi*; Suzuki, Takashi; Kuwabara, Jun; Nakanishi, Takahiro
Journal of Environmental Radioactivity, 192, p.208 - 218, 2018/12
From August 2011 to October 2013, the concentration of iodine-129 (I) in the seabed sediment collected from 26 stations located within 160 km from the Fukushima Daiichi Nuclear Power Plant was estimated. The concentrations of I in seabed sediment off Fukushima in 2011 ranged between 0.02 and 0.45 mBq/kg. Although iodine is a biophilic element, the accident-derived radioiodine negligibly affects the benthic ecosystem. Until October 2013, a slightly increased activity of I in the surface sediment in the shelf-edge region (bottom depth: 200-400 m) was observed. The increase of the I concentrations in the shelf-edge sediments was affected by the (1) transport of I-bound particles from the land through rivers and (2) re-deposition of I desorbed from the contaminated coastal sediment to the shelf-edge sediments, which were considered to be dominant processes.
Nordlund, K.*; Zinkle, S. J.*; Sand, A. E.*; Granberg, F.*; Averback, R. S.*; Stoller, R. E.*; Suzudo, Tomoaki; Malerba, L.*; Banhart, F.*; Weber, W. J.*; et al.
Journal of Nuclear Materials, 512, p.450 - 479, 2018/12
Scientific understanding of any kind of radiation effects starts from the primary damage. We consider the extensive experimental and computer simulation studies that have been performed over the past several decades on what the nature of the primary damage is. We review both the production of crystallographic or topological defects in materials as well as radiation mixing, i.e. the process where atoms in perfect crystallographic positions exchange positions with other ones in non-defective positions. We also consider the recent effort to provide alternatives to the current international standard for quantifying this energetic particle damage, the Norgett-Robinson-Torrens displacements per atom (NRT-dpa) model for metals. We present in detail new complementary displacement production estimators ("athermal recombination corrected dpa": arc-dpa) and atomic mixing ("replacements per atom": rpa) functions that extend the NRT-dpa, and discuss their advantages and limitations.
Kusaka, Ryoji; Watanabe, Masayuki
Physical Chemistry Chemical Physics, 20(47), p.29588 - 29590, 2018/12
Mechanistic understanding of solvent extraction of uranyl ions (UO) by tributyl phosphate (TBP) will help improve the technology for the treatment and disposal of spent nuclear fuels. So far, it has been believed that uranyl ions in the aqueous phase are adsorbed to a TBP-enriched organic/aqueous interface, form complexes with TBP at the interface, and are extracted into the organic phase. Here we show that uranyl-TBP complex formation does not take place at the interface using vibrational sum frequency generation (VSFG) spectroscopy and propose an alternative extraction mechanism that uranyl nitrate, UO(NO), passes through the interface and forms the uranyl-TBP complex, UO(NO)(TBP), in the organic phase.
Kimura, Taiki*; Kaneko, Masashi; Watanabe, Masayuki; Miyashita, Sunao*; Nakashima, Satoru*
Dalton Transactions, 47(42), p.14924 - 14931, 2018/11
We demonstrated density functional calculations of Eu(III) and Am(III) complexes with pnictogen-donor (X) ligands, CH)X-CH-CH-X(CH) (X = N, P, As and Sb). We investigated the optimized structures of the cmoplexes and the Gibbs energy differences in the complex formation reactions. Those results indicated that the N- and P-donor ligands have Am(III) ion selectivity over Eu(III) ion, especially, the P-donor ligand showed the highest selectivity. The tendency of the Am(III)/Eu(III) selectivity by the pnictogen-dono ligands was found to be comparable to that of soft acid classification in hard and soft acids and bases rule. Mulliken's spin population analysis indicated that the bonding property between the metal ion and the pnictogen atoms correlated with the Am(III)/Eu(III) selectivity. In particular, the participation of f-orbital electrons of the metal ion in the covalency was indicated to have an important role for the selectivity.
Iwamoto, Hiroki; Stankovskiy, A.*; Fiorito, L.*; Van den Eynde, G.*
European Physical Journal; Nuclear Sciences & Technologies (Internet), 4, p.42_1 - 42_7, 2018/11
This paper presents a nuclear data sensitivity and uncertainty analysis of the effective delayed neutron fraction for critical and subcritical cores of the MYRRHA reactor using the continuous-energy Monte Carlo transport code MCNP. The sensitivities are calculated by the modified -ratio method proposed by Chiba. Comparing the sensitivities obtained with different scaling factors introduced by Chiba shows that a value of is the most suitable for the uncertainty quantification of . Using the calculated sensitivities and the JENDL-4.0u covariance data, the uncertainties for the critical and subcritical cores are determined to be 2.2 0.2% and 2.0 0.2%, respectively, which are dominated by delayed neutron yield of Pu and U.
Tsuda, Shuichi; Tanigaki, Minoru*; Yoshida, Tadayoshi; Saito, Kimiaki
Hoshasen, 44(3), p.109 - 118, 2018/11
JAEA has started to perform dose rate monitoring using a car-borne survey system KURAMA to rapidly produce the dose rate mappings of the deposited radionuclides in the environment after the Fukushima Dai-ichi Nuclear Power Plant accident. KURAMA is a car-borne survey system developed by Kyoto University to perform dose rate monitoring in a wide area in detail with rapidity. By improving KURAMA with continuous dose rate monitoring, the 2nd generation of KURAMA (KURAMA-II) succeeded in downsizing, durability and automated transmission of data so that enable detailed dose rate mapping in wide area in shorter period of time. This paper reports the radiation characteristics and the simulation analysis of KURAMA-II on the special issue of Hoshasen, the journal of Ionization Radiation Division in the Japan society of applied physics.
Kai, Takeshi; Yokoya, Akinari*; Fujii, Kentaro*; Watanabe, Ritsuko*
Hoshasen Kagaku (Internet), (106), p.21 - 29, 2018/11
It is thought to that the biological effects such as cell death or mutation are induced by complex DNA damage which are formed by several damage sites within a few nm. As the prediction of complex DNA damage at an electron track end, we report our outcomes. These results indicate that DNA damage sites comprising multiple nucleobase lesions with a single strand breaks can be formed by multiple collisions of the electrons within 1 nm. This multiple damage site cannot be processed by base excision repair enzymes. Pre-hydrated electrons can also be produced resulting in an additional base lesion over a few nm from the multi-damage site. This clustered damage site may be finally converted into a double strand break. These DSBs include another base lesion(s) at their termini that escape from the base excision process and which may result in biological effect. Our simulation is useful to reveal phenomena involved in radiation physico-chemistry as well as the DNA damage prediction.
Kobayashi, Hideki*; Suzuki, Rikie*; Yang, W.*; Ikawa, Hiroki*; Inoue, Tomoharu*; Nagano, Hirohiko; Kim, Y.*
Polar Data Journal (Internet), 2, p.14 - 29, 2018/11
The Arctic and boreal regions have been experiencing a rapid warming in the 21st century. It is important to understand the dynamics of boreal forest at the continental scale under the climate and environmental changes. While the role of understory vegetation in boreal forest ecosystems on carbon and nutrient cycling cannot be ignored, they are still one of least understood components in boreal ecosystems. Spectroscopic measurements of vegetation are useful to identify species and their biochemical characteristics. In this data paper, we present spectral reflectances of 44 typical understory formations and five 30-m long transects. The spectral reflectance covers the spectral region of visible, near infrared and shortwave infrared (350-2500 nm). For the transect measurements, we decided the length of transect at 30 m, similar to the scale of one pixel of a Landsat type satellite imagery. The photographs at all positions, where spectral reflectances were obtained, are included to understand the structure and status of each sample. The data set contains six dwarf shrubs (blueberry (), cowberry (), feltleaf willow (), young birch (), young aspen (), and young black spruce ()), two herbaceous (cottongrass () and marsh Labrador tea ()), three mosses (Sphagnum moss, splendid feather moss (), and polytrichum moss ()), and reindeer lichen (). Spectral reflectances from several non-vegetative such as snow, litter, and soil are also included. This spectral and photographic data set can be used for understanding the spectral characteristics of understory formations, designing newly planned spectral observations, and developing and validating the remote sensing methodology of large scale understory monitoring.
Uesawa, Shinichiro; Yamashita, Susumu; Shibata, Mitsuhiko; Yoshida, Hiroyuki
Proceedings of 11th Korea-Japan Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS-11) (Internet), 6 Pages, 2018/11
A dry method for fuel debris is proposed for decommissioning of TEPCO's Fukushima Daiichi NPS. We have been evaluating the air-cooling performance of the fuel debris in the dry method by using JUPITER. Because JUPITER can calculate relocation of the corium, it is expected to calculate thermal-hydraulic simulation of the air cooling of the fuel debris in the dry method based on the calculated debris position, shapes and composition with the relocation analysis. In this paper, the experiment of heat transfer and flow visualization of free convection adjacent to upward-facing horizontal heat transfer surface was performed to validate the calculation of the free convective heat transfer with JUPITER. In the experiment, the temperature distribution was measured with a thermocouple tree. In addition, the velocity distribution of free convection was visualized by a particle image velocimetry (PIV). In the comparison between the JUPITER and the experiment, the temperature distribution for the vertical direction in the quasi-steady state was fitted between the JUPITER and the experiment. The velocity distribution calculated with JUPITER was also in good agreement with the experimental result. Therefore, it is expected that JUPITER is a helpful numerical method to evaluate the air-cooling performance of the fuel debris in the dry method.
Yoshida, Hiroyuki; Uesawa, Shinichiro; Horiguchi, Naoki; Miyahara, Naoya; Ose, Yasuo*
Proceedings of 11th Korea-Japan Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS-11) (Internet), 6 Pages, 2018/11
The radioactive aerosol removal equipment is used as one of the safety systems of nuclear reactors. In this equipment, micro particles of aerosol remove through gas-liquid interfaces of two-phase flow. The mechanism related to the removal of micro particles through the gas-liquid interface is not clear, a numerical evaluation method of performance of aerosol removal equipment is not realized. Then, we have started to construct a numerical simulation method to simulate removal of micro particles through gas-liquid interfaces. In this simulation method, detailed two-phase flow simulation code TPFIT is used as the basis of this method. TPFIT adopts an advanced interface tracking method and can simulate interface movement and deformation directly. In addition, to simulate the movement of particles, the Lagrangian particle tracking method is incorporated. By combining the interface tracking method and the Lagrangian particle tracking method, the interaction between interfaces and micro particles can be simulated in detail. To solve the Lagrangian equations of particles, fluid properties and fluid velocity surrounding aerosol particles are evaluated by considering the relative position of particles and gas-liquid interface, to simulate particle movement near the interface. In this paper, outline and preliminary results of this simulation method are shown.
Yamashita, Susumu; Yoshida, Hiroyuki
Proceedings of 11th Korea-Japan Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS-11) (Internet), 5 Pages, 2018/11
To investigate detailed flow behaviors around the beam window of accelerator driven system (ADS), large scale simulation for unsteady thermal hydraulics around the beam window was performed using JUPITER. The input data, such as the beam window and nozzle, is designed as accurate as possible using the computer aided design software. As a result, the flow pattern around the beam window is quite different from previous results in which the steady flow is assumed. The flow pattern of the Lead-Bismuth Eutectic around the beam window and the exit of the nozzle were very complicated. According to complicated flow around those structures, the temperature distribution on the beam window is also complicated. Thus, it is confirmed that the complicated flow around structures will affect to the temperature distribution in the structures and the effect of flow field on the temperature must be evaluated.
Nagatake, Taku; Shibata, Mitsuhiko; Uesawa, Shinichiro; Ono, Ayako; Yoshida, Hiroyuki
Proceedings of 11th Korea-Japan Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS-11) (Internet), 3 Pages, 2018/11
In the Fukushima Daiichi Nuclear Power Plant accident, reactor cores were cooled by natural circulation due to pump trip. To investigate the accident progress of the Fukushima Daiichi Nuclear Power Plant, it is important to understand the thermal hydraulic behavior in reactor cores including fuel bundles. Flow rate inside cores was relatively low in the natural circulation conditions, then, thermal-hydraulic behavior in the fuel bundles was different from that in the normal operating conditions. To evaluate thermal hydraulic behavior under the accidental conditions, we are developing the numerical simulation codes named TPFIT and ACE3D. These codes are based on two-phase computational fluid dynamics and can simulate the two-phase flow inside fuel bundles including low flow rate condition. Before applying these codes to the thermal-hydraulic behavior, the applicability of these codes must be confirmed. Then, in this study, in order to obtain a validation data for TPFIT and ACE3D code, thermal hydraulic experiment was performed by using test section with a simulated fuel bundle with 44 unheated rods. In this simulated fuel bundle, there were wire mesh sensors, and void fraction distribution data inside the simulated fuel bundle under high pressure condition (max. 2.6 MPa) was obtained. The one of the advantage of wire mesh sensor is that a void fraction distribution of cross section at the same time can be measured. In this paper, void fraction distribution of two-phase flow in a simulated fuel bundle under high pressure condition are reported.
Ozu, Akira; Maeda, Makoto; Komeda, Masao; Toh, Yosuke; Koizumi, Mitsuo; Seya, Michio
Proceedings of 2017 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC 2017) (Internet), 4 Pages, 2018/11
A Differential Die-away Analysis (DDA) system using a compact pulsed neutron (14 MeV) generator has been newly developed for non-nuclear proliferation and nuclear security in the Japan Atomic Energy Agency (JAEA). The DDA system was designed to be able to detect a nuclear fissile material (Pu-239) of as low as 10 mg and to handle samples of a different volume: a vial bottle (20 mL), a pail container (20 L), through a Monte Carlo simulation. In the DDA system, the Fast Neutron Direct Interrogation (FNDI) technique, which utilizes fast neutrons for interrogation, was applied to measure the amount of fissile mass contained in the sample. The fundamental performance of the DDA system was investigated in the demonstration experiment. The simulation results show that the Pu-239 masses of less than 10 mg can be detected in the DDA system. The results of the experiment are discussed and compared with those of the simulation.