Narukawa, Takafumi; Amaya, Masaki
Journal of Nuclear Science and Technology, 56(7), p.650 - 660, 2019/07
Asahi, Yuichi*; Grandgirard, V.*; Sarazin, Y.*; Donnel, P.*; Garbet, X.*; Idomura, Yasuhiro; Dif-Pradalier, G.*; Latu, G.*
Plasma Physics and Controlled Fusion, 61(6), p.065015_1 - 065015_15, 2019/05
The role of poloidal convective cells on transport processes is studied with the full-F gyrokinetic code GYSELA. For this purpose, we apply a numerical filter to convective cells and compare the simulation results with and without the filter. The energy flux driven by the magnetic drifts turns out to be reduced by a factor of about 2 once the numerical filter is applied. A careful analysis reveals that the frequency spectrum of the convective cells is well-correlated with that of the turbulent Reynolds stress tensor, giving credit to their turbulence-driven origin. The impact of convective cells can be interpreted as a synergy between turbulence and neoclassical dynamics.
Simonnet, M.; Barr, N.*; Drot, R.*; Le Naour, C.*; Sladkov, V.*; Delpech, S.*
Radiochimica Acta, 107(4), p.289 - 297, 2019/04
Negyesi, M.; Amaya, Masaki
Annals of Nuclear Energy, 114, p.52 - 65, 2018/04
Kitano, Akihiro; Nakajima, Ken*
Proceedings of 2018 International Congress on Advances in Nuclear Power Plants (ICAPP 2018) (CD-ROM), p.1205 - 1210, 2018/04
The feedback reactivity is taken into account in fast reactor core design especially in order to make the power coefficient negative, which is required to be confirmed in the operation. In the feedback reactivity experiment, the positive reactivity was inserted in the critical state at zero power, and the thermal data, such as reactor power and the R/V inlet temperature, was acquired until the power got stable by the feedback reactivity. In the conventional study, only two critical points in an experiment are available for evaluation of the feedback reactivity coefficients. This method needs three days for evaluation. The advanced method based on the inverse kinetics is newly applied in this work using the more extensive data. It is confirmed that this approach can evaluate the feedback reactivity coefficients in one experiment.
Negyesi, M.; Amaya, Masaki
Journal of Nuclear Science and Technology, 54(10), p.1143 - 1155, 2017/10
Miyahara, Naoya; Miwa, Shuhei; Nakajima, Kunihisa; Osaka, Masahiko
Proceedings of 2017 Water Reactor Fuel Performance Meeting (WRFPM 2017) (USB Flash Drive), 9 Pages, 2017/09
This paper presents the development of a reproductive experimental setup for FP release and transport and an analysis tool considering chemical reaction kinetics for the construction of the FP chemistry database. The performance test of the reproductive experimental setup TeRRa using CsI compounds show that TeRRa can reproduce well a FP chemistry-related behavior such as aerosol formation, growth and deposition behavior. An analytical tool has been developed based on the commercial ANSYS-FLUENT code. Some additional models was added to evaluate detailed FP chemistry during release and transport in this study. A test analysis simulating the CsI heating test in steam atmosphere was carried out to demonstrate the performance of the improved code. The result shows the appropriateness of the additional models.
Subekti, M.*; Kudo, Kazuhiko*; Nabeshima, Kunihiko; Takamatsu, Kuniyoshi
Atom Indonesia, 43(2), p.93 - 102, 2017/08
Reactor kinetics based on point kinetic model have been generally applied as the standard method for neutronics codes. As the central control rod (C-CR) withdrawal test has demonstrated in a prismatic core of HTTR, the transient calculation of kinetic parameter, such as reactivity and neutron fluxes, requires a new method to shorten calculation-process time. Development of neural network method was applied to point kinetic model as the necessity of real-time calculation that could work in parallel with the digital reactivity meter. The combination of TDNN and Jordan RNN, such as TD-Jordan RNN, was the result of the modeling approach. The application of TD-Jordan RNN with adequate learning, tested offline, determined results accurately even when signal inputs were noisy. Furthermore, the preprocessing for neural network input utilized noise reduction as one of the equations to transform two of twelve time-delayed inputs into power corrected inputs.
Asahi, Yuichi*; Latu, G.*; Ina, Takuya; Idomura, Yasuhiro; Grandgirard, V.*; Garbet, X.*
IEEE Transactions on Parallel and Distributed Systems, 28(7), p.1974 - 1988, 2017/07
High-dimensional stencil computation from fusion plasma turbulence codes involving complex memory access patterns, the indirect memory access in a Semi-Lagrangian scheme and the strided memory access in a Finite-Difference scheme, are optimized on accelerators such as GPGPUs and Xeon Phi coprocessors. On both devices, the Array of Structure of Array (AoSoA) data layout is preferable for contiguous memory accesses. It is shown that the effective local cache usage by improving spatial and temporal data locality is critical on Xeon Phi. On GPGPU, the texture memory usage improves the performance of the indirect memory accesses in the Semi-Lagrangian scheme. Thanks to these optimizations, the fusion kernels on accelerators become 1.4x - 8.1x faster than those on Sandy Bridge (CPU).
Kido, Kentaro; Hata, Kuniki; Maruyama, Yu; Nishiyama, Yutaka; Hoshi, Harutaka*
NEA/CSNI/R(2016)5 (Internet), p.204 - 212, 2016/05
Di Lemma, F. G.; Miwa, Shuhei; Osaka, Masahiko
JAEA-Review 2016-007, 27 Pages, 2016/03
During a nuclear power plant Severe Accident, complex boron melts can be formed, due to interaction of the control rods with the cladding materials. These can affect ultimately the source term assessment. This review will describe the results of previous studies on boron carbide/stainless steel/Zircaloy (BC/SS/Zry) melts, which will finally provide guidance for the needs of future experiments. This review showed that models for the behavior of complex BC/SS/Zry melts are limited, and unsuccessful in simulating core degradation, thus the improvement of the database for BC/SS/Zry melts is needed. Our experimental plan aims in providing thermodynamics and kinetics models for such melts, with the final aim of improving boron modelling in SA codes analysis and of understanding its effect on fission products behavior.
Kawaguchi, Munemichi; Doi, Daisuke; Seino, Hiroshi; Miyahara, Shinya
Proceedings of 23rd International Conference on Nuclear Engineering (ICONE-23) (DVD-ROM), 6 Pages, 2015/05
CONTAIN-LMR code is an integrated analysis tool to predict the consequence of severe accident in a liquid metal fast reactor. A sodium-concrete reaction is one of the most important phenomena, and Sodium-Limestone Concrete Ablation Model (SLAM) has been installed into the original CONTAIN code. The SLAM treats chemical reaction kinetics between the sodium and the concrete compositions mechanistically, the application is limited to the limestone concrete. In order to apply SLAM to the siliceous concrete which is an ordinary structural concrete in Japan, the chemical reaction kinetics model has been improved. The improved model was validated to analyze a series of sodium-concrete experiments which were conducted in Japan Atomic Energy Agency. It has been found that relatively good agreement between calculation and experimental results is obtained and the CONTAIN-LMR code has been validated with regard to the sodium-concrete reaction phenomena.
Takamatsu, Kuniyoshi; Tochio, Daisuke; Nakagawa, Shigeaki; Takada, Shoji; Yan, X.; Sawa, Kazuhiro; Sakaba, Nariaki; Kunitomi, Kazuhiko
Journal of Nuclear Science and Technology, 51(11-12), p.1427 - 1443, 2014/11
In a safety demonstration test involving a loss of both reactor reactivity control and core cooling, HTGRs such as the HTTR, which is the only HTGR in Japan, demonstrate that the reactor power would stabilize spontaneously. In the test at an initial power of 30%, when the insertion of all control rods was disabled and all gas circulators were tripped to reduce the coolant flow rate to zero, a reactor transient was initiated and examined. The results confirmed that the reactor power would decrease immediately and become effectively zero.
Idomura, Yasuhiro; Tokuda, Shinji; Kishimoto, Yasuaki
Nuclear Fusion, 45(12), p.1571 - 1581, 2005/12
Using a global gyrokinetic toroidal particle code, the toroidal electron temperature gradient driven (ETG) turbulence is studied in positive and reversed shear tokamaks. In the nolinear turbulent state, the ETG turbulence in the positive and reversed shear configurations show quite different structure formations. In the positive shear configuration, the ETG turbulence is dominated by streamers which have a ballooning type structure, and the electron temperature profile is quickly relaxed to the marginally stable state in a turbulent time scale. In the reversed shear configuration, quasi-steady zonal flows are produced in the regative shear region, while the positive shear region is characterized by streamers. Accordingly, the electron thermal diffusivity has a gap structure across the surface, and the gradinet is sustained above the marginal value for a long time in the quasi-steady phase. The results suggest a stiffness of the profile in positive shear tokamaks, and a possibility of the Te transport barrier in reversed shear tokamaks.
Purazuma, Kaku Yugo Gakkai-Shi, 81(8), p.581 - 592, 2005/08
A gyrokinetic particle simulation is a powerful tool in studying tokamak microturbulence. A method, which is a standard method in recent gyrokinetic particle simulations, dramatically improved an efficiency of a particle simulation by reducing a particle noise, and full torus turbulence simulations are enabled. In this paper, the method is reviewed, and issues in full torus gyrokinetic particle simulations are discussed.
Jolliet, S.*; Angelino, P.*; Bottino, A.*; Idomura, Yasuhiro; Villard, L.*
Theory of Fusion Plasmas, ISPP21, p.345 - 351, 2004/00
Global particle-in-cell (PIC) simulations are a very useful tool for studying the time evolution of turbulence induced by ion-temperature-gradient (ITG) instabilities. Unfortunately, the linear code LORB5 and its non-linear version ORB5 require high computational power. In order to study more sophisticated models, we need to optimize these codes. We will focus on LORB5, which uses a cylindrical grid (r,z) for solving the Vlasov equation and a (s,) grid for the Poisson equation. The approach presented in this work consists of implementing the gyrokinetic model using a single (s,) grid. Here is the straight-field-line poloidal coordinate. A method to avoid the singularity at the magnetic axis is presented, and a benchmark with the CYCLONE case is shown.
Okada, Michio*; Moritani, Kosuke; Goto, Seishiro*; Kasai, Toshio*; Yoshigoe, Akitaka; Teraoka, Yuden
Journal of Chemical Physics, 119(14), p.6994 - 6997, 2003/10
The oxidation of Cu(001) with hyperthermal O molecular beams was investigated using X-ray photoemission spectroscopy in conjunction with a synchrotron light source. The efficiency of oxidation is higher than that with ambient thermal O. Further oxidation under oxygen coverage larger than 0.5 ML occurs rather inefficiently even for the 2.3 eV beam irradiation. We found such slow oxidation of Cu corresponding to the initial stage of the CuO formation can be interpreted in terms of a collision-induced-adsorption mechanism. The kinetics of the dissociative adsorption is well described using the first order kinetics in a simple Langmuir-type adsorption model.
Norisue,Tomohisa*; Kida, Yusuke*; Masui, Naoki*; Tran-Cong-Miyata, Q.*; Maekawa, Yasunari; Yoshida, Masaru; Shibayama, Mitsuhiro*
Macromolecules, 36(16), p.6202 - 6212, 2003/08
The shrinking kinetics of poly(N-isopropylacrylamide) (PNIPA) gels has been studied for two types of PNIPA gels prepared by (i) copolymerization of constituent monomer and cross-linker and (ii) -ray irradiation in the PNIPA solutions in order to investigate the role of cross-linking on shrinking kinetics. The shrinking kinetics of the monomer cross-linked gels is quite similar to that of the polymer cross-linked gels. On the other hand, a significant difference was found when the microscopic structure and the dynamics were investigated by small-angle neutron scattering (SANS) and static/dynamic light scattering (SLS/DLS). The degree of built-in inhomogeneities and dynamic fluctuations were evaluated as a function of the cross-linking degree and the gel preparation temperature by intensity decomposition methods for both types of gels. It is concluded that the monomer cross-linked gels have extra built-in inhomogeneities due to the spatial distribution of crosslinks in addition to the frozen concentration fluctuations inherent in polymer gels.
Yoshigoe, Akitaka; Moritani, Kosuke; Teraoka, Yuden
Japanese Journal of Applied Physics, Part 1, 42(6B), p.3976 - 3982, 2003/06
It is well known that the thermal oxidation on Si(001) surface is an important reaction system to form of the gate-oxide films in MOSFET, since it is necessary to control the film thickness under a few nano-meter scale. Thus, we have studied the oxygen uptake and the Si oxidation states depending on the oxidation times by using the synchrotron radiation photoemission spectroscopy in 110Pa of O at the surface temperature from 870K to 1120K. We clarified the oxidation depending on the surface temperature was explained by the kinetics (Langumuir and auto-catalytic model). Using real time photoemission spectroscopy, we found that the Si sepcies was not formed at the initial oxidation stage.
Nagase, Fumihisa; Otomo, Takashi; Uetsuka, Hiroshi
Journal of Nuclear Science and Technology, 40(4), p.213 - 219, 2003/04
Isothermal oxidation tests in flowing steam were performed on low-Sn Zircaloy-4 cladding tubes over the wide temperature range from 773 to 1573 K in order to obtain oxidation kinetics applicable to various loss-of-coolant accident conditions of LWRs. The oxidation generally obeys a parabolic rate law for the examined time range up to 3600s at temperatures from 1273 K to 1573K, and for a limited time range up to 900s from 773 to 1253 K. A cubic rate law is preferable for evaluating the longer-term oxidation at 1253 K and below. The parabolic rate law constant and the cubic rate law constant for measured weight gain were evaluated at every examined temperature, and Arrhenius-type equations were determined in order to describe the temperature dependence of the rate constants. It was indicated that the change of the oxidation kinetics from the cubic to the parabolic rate and the discontinuities in the temperature dependence of the rate constants are caused by the monoclinic/tetragonal phase transformation of ZrO.