Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Idomura, Yasuhiro; Dif-Pradalier, G.*; Garbet, X.*; Sarazin, Y.*; Tore Supra Teams*
Physics of Plasmas, 30(4), p.042508_1 - 042508_18, 2023/04
Times Cited Count:2 Percentile:29.50(Physics, Fluids & Plasmas)Two time slices in linear and saturated ohmic confinement (LOC and SOC) regimes in a Tore Supra ohmic L-mode discharge are analyzed using nonlinear global full- gyrokinetic simulations, and qualitative features of the LOC-SOC transition are reproduced. The exhaust of carbon impurity ions is caused by fast ion mixing, which is driven by the toroidal field stress. The intrinsic rotation develops in the opposite direction between the LOC and SOC phases, which is characterized by the different features of the mode asymmetry between trapped electron modes in the LOC phase and ion temperature gradient driven modes in the SOC phase, leading to the change of the profile shear stress. The energy fluxes of electrons and deuterium ions are dominant in the LOC and SOC phases, respectively, and the ratio of the energy confinement time between two phases agree with the experimental value.
Asahi, Yuichi; Fujii, Keisuke*; Heim, D. M.*; Maeyama, Shinya*; Garbet, X.*; Grandgirard, V.*; Sarazin, Y.*; Dif-Pradalier, G.*; Idomura, Yasuhiro; Yagi, Masatoshi*
Physics of Plasmas, 28(1), p.012304_1 - 012304_21, 2021/01
Times Cited Count:6 Percentile:36.28(Physics, Fluids & Plasmas)This article demonstrates a data compression technique for the time series of five dimensional distribution function data based on Principal Component Analysis (PCA). Phase space bases and corresponding coefficients are constructed by PCA in order to reduce the data size and the dimensionality. It is shown that about 83% of the variance of the original five dimensional distribution can be expressed with 64 components. This leads to the compression of the degrees of freedom from to
. One of the important findings - resulting from the detailed analysis of the contribution of each principal component to the energy flux - deals with avalanche events, which are found to be mostly driven by coherent structures in the phase space, indicating the key role of resonant particles.
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
Times Cited Count:5 Percentile:27.82(Physics, Fluids & Plasmas)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.
Donnel, P.*; Garbet, X.*; Sarazin, Y.*; Asahi, Yuichi; Wilczynski, F.*; Caschera, E.*; Dif-Pradalier, G.*; Ghendrih, P.*; Gillot, C.*
Plasma Physics and Controlled Fusion, 61(1), p.014003_1 - 014003_11, 2019/01
Times Cited Count:13 Percentile:61.32(Physics, Fluids & Plasmas)Poloidal asymmetries of the plasma flow are known to play a role in neoclassical transport. According to conventional neoclassical theory, the level of poloidal asymmetry of the electric potential is expected to be very small. In the present work, a general framework for the generation of axisymmetric structures of potential by turbulence is presented. Zonal flows, geodesic acoustic modes and convective cells are described by a single model. This is done by solving the gyrokinetic equation coupled to the quasi-neutrality equation. This calculation provides a predictive calculation of the frequency spectrum of flows given a specified forcing due to turbulence. It also shows that the dominant mechanism comes from zonal flow compression at intermediate frequencies, while ballooning of the turbulence Reynolds stress appears to be the main drive at low frequency.
Asahi, Yuichi*; Grandgirard, V.*; Idomura, Yasuhiro; Garbet, X.*; Latu, G.*; Sarazin, Y.*; Dif-Pradalier, G.*; Donnel, P.*; Ehrlacher, C.*
Physics of Plasmas, 24(10), p.102515_1 - 102515_17, 2017/10
Times Cited Count:10 Percentile:44.54(Physics, Fluids & Plasmas)Two full-F global gyrokinetic codes are benchmarked to compute flux-driven ion temperature gradient turbulence in tokamak plasmas. For this purpose, the Semi-Lagrangian code GYSELA and the Eulerian code GT5D are employed, which solve the full-F gyrokinetic equation with a realistic fixed flux condition. Using the appropriate settings for the boundary and initial conditions, flux-driven ITG turbulence simulations are carried out. The avalanche-like transport is assessed with a focus on spatio-temporal properties. A statistical analysis is performed to discuss this self-organized criticality (SOC) like behaviors, where we found spectra and a transition to
spectra at high-frequency side in both codes. Based on these benchmarks, it is verified that the SOC-like behavior is robust and not dependent on numerics.
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
Times Cited Count:7 Percentile:51.62(Computer Science, Theory & Methods)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).
Litaudon, X.*; Sakamoto, Yoshiteru; de Vries, P. C.*; Salmi, A.*; Tala, T.*; Angioni, C.*; Benkadda, S.*; Beurskens, M. N. A.*; Bourdelle, C.*; Brix, M.*; et al.
Nuclear Fusion, 51(7), p.073020_1 - 073020_13, 2011/07
Times Cited Count:11 Percentile:42.10(Physics, Fluids & Plasmas)A variety of triggering mechanisms and structures of internal transport barrier (ITB) has been observed in various devices or depending on operation scenarios. Thus identity experiments on ITB in JT-60U and JET have been performed to shed light on the physics behind ITBs. Because of their similar size, the dimensionless parameters between both devices are the same. These experiments were performed with near identical magnetic configurations, heating waveforms and normalized quantities such as safety factor, magnetic shear, normalized Larmor radius, normalized collision frequency, beta, temperatures ratio. Similarities of the ITB triggering mechanism and the ITB strength have been observed when a proper match is achieved of the most relevant profiles of the normalized quantities. This paper will report on the detail comparison of transport properties of ITBs obtained in these JET/JT-60U identity experiments.
Lesur, M.*; Idomura, Yasuhiro; Shinohara, Koji; Garbet, X.*
Proceedings of 23rd IAEA Fusion Energy Conference (FEC 2010) (CD-ROM), 8 Pages, 2011/03
Alfvn Eigenmodes (AEs) with frequency sweeping (chirping AEs) are analyzed based on the so-called Berk-Breizman (BB) model. Kinetic parameter regimes with chirping nonlinear solutions are delimited, and a new quasi-periodic chirping regime is found, where the quasi-period of chirping events depends on the linear growth rate. Fundamental kinetic parameters such as the linear drive and the external damping rate are estimated by fitting nonlinear chirping characteristics between the experiment and the BB model. This approach is applied to Toroidicity-induced AEs on JT-60U, which suggests the existence of modes far from marginal stability. Two collision models are considered, and it is shown that dynamical friction and velocity-space diffusion are essential to reproduce nonlinear features observed in experiments. The results are validated by recovering measured growth and decay of perturbation amplitude, and by estimating collision frequencies from experimental equilibrium data.
Lesur, M.*; Idomura, Yasuhiro; Shinohara, Koji; Garbet, X.*; JT-60 Team
Physics of Plasmas, 17(12), p.122311_1 - 122311_9, 2010/12
Times Cited Count:35 Percentile:76.24(Physics, Fluids & Plasmas)Villard, L.*; Bottino, A.*; Brunner, S.*; Casati, A.*; Chowdhury, J.*; Dannert, T.*; Ganesh, R.*; Garbet, X.*; Grler, T.*; Grandgirard, V.*; et al.
Plasma Physics and Controlled Fusion, 52(12), p.124038_1 - 124038_18, 2010/11
Times Cited Count:21 Percentile:60.19(Physics, Fluids & Plasmas)Garbet, X.*; Idomura, Yasuhiro; Villard, L.*; Watanabe, Tomohiko*
Nuclear Fusion, 50(4), p.043002_1 - 043002_30, 2010/04
Times Cited Count:290 Percentile:74.19(Physics, Fluids & Plasmas)This overview is an assessment of the gyrokinetic framework and simulations to compute the turbulent transport in fusion plasmas. It covers an introduction to the gyrokinetic theory, the principal numerical techniques which are being used to solve the gyrokinetic equations, fundamentals in gyrokinetic turbulence and the main results which have been brought by simulations with regard to transport in fusion devices and actuation measurements.
Lesur, M.*; Idomura, Yasuhiro; Garbet, X.*
Physics of Plasmas, 16(9), p.092305_1 - 092305_12, 2009/09
Times Cited Count:40 Percentile:79.58(Physics, Fluids & Plasmas)Doyle, E. J.*; Houlberg, W. A.*; Kamada, Yutaka; Mukhovatov, V.*; Osborne, T. H.*; Polevoi, A.*; Bateman, G.*; Connor, J. W.*; Cordey, J. G.*; Fujita, Takaaki; et al.
Nuclear Fusion, 47(6), p.S18 - S127, 2007/06
no abstracts in English
Connor, J. W.*; Fukuda, Takeshi*; Garbet, X.*; Gormezano, C.*; Mukhovatov, V.*; Wakatani, Masahiro*; ITB Database Group; ITPA Topical Group on Transport and Internal Barrier Physics*
Nuclear Fusion, 44(4), p.R1 - R49, 2004/04
Times Cited Count:322 Percentile:74.42(Physics, Fluids & Plasmas)This paper first reviews the present state of theoretical and experimental knowledge regarding the formation and characteristics of ITBs in tokamaks. Specifically, the current status of theoretical modeling of ITBs is presented; then, an international ITB database based on experimental information extracted from some nine tokamaks is described and used to draw some general conclusions concerning the necessary conditions for ITBs to appear, comparing these with the theoretical models. The experimental situation regarding the steady-state, or at least quasi-steady-state, operation of tokamaks is reviewed and finally the issues and prospects for achieving such operational modes in ITER are discussed.
Bcoulet, M.*; Huysmans, G.*; Sarazin, Y.*; Garbet, X.*; Ghendrih, P.*; Rimini, F.*; Joffrin, E.*; Litaudon, X.*; Monier-Garbet, P.*; An
, J.-M.*; et al.
Plasma Physics and Controlled Fusion, 45(12A), p.A93 - A113, 2003/12
Times Cited Count:86 Percentile:90.58(Physics, Fluids & Plasmas)no abstracts in English
Asahi, Yuichi*; Grandgirard, V.*; Idomura, Yasuhiro; Sarazin, Y.*; Latu, G.*; Garbet, X.*
no journal, ,
This talk reviews outcomes from BMTFF projects, which was conducted for FY2015-2016. In this project, in order to establish a firm basis of full-f gyrokinetic models, two major full-f gyrokinetic codes in EU and Japan, GYSELA and GT5D, were benchmarked. In FY2015, all the numerical implementations were examined, and boundary conditions were fixed to be the same. With this correction, collisional transport, linear zonal flow damping, and linear stability of the ion temperature gradient driven (ITG) mode were successfully benchmarked. In FY2016, the same source and sink models were implemented in both codes, and nonlinear turbulence simulations were benchmarked. Decaying ITG turbulence simulations without heat sources showed similar profile relaxation processes, and nonlinear critical temperature gradients agreed quantitatively with each other. On the other hand, driven ITG turbulence simulations with heat sources showed intermittent bursts of avalanche like transport, which indicate similar 1/f type frequency spectra.
Asahi, Yuichi; Maeyama, Shinya*; Latu, G.*; Garbet, X.*; Watanabe, Tomohiko*; Aoki, Takayuki*; Ogino, Masao*
no journal, ,
Under the JHPCN international collaboration with French researchers, we have integrated a kinetic electron model into the kinetic plasma turbulence code GYSELA developed in France. We also explore a performance portable for a kinetic fusion plasma turbulence code like GYSELA. For this purpose, we extract the key features of GYSELA such as the high dimensionality and the semi-Lagrangian scheme, and encapsulate them into a mini-application which solves the simpler version of Vlasov-Poisson system as GYSELA. We implement the mini-app with Open ACC and Kokkos to evaluate the advantages and disadvantages of each approach, from the view point of performance portability, readability and productivity.
Asahi, Yuichi; Maeyama, Shinya*; Bigot, J.*; Garbet, X.*; Grandgirard, V.*; Fujii, Keisuke*; Shimokawabe, Takashi*; Watanabe, Tomohiko*; Idomura, Yasuhiro; Onodera, Naoyuki; et al.
no journal, ,
We have established an in-situ data analysis method for large scale fluid simulation data and developed deep learning based surrogate models to predict fluid simulation results. Firstly, we have developed an in-situ data processing approach, which loosely couples the MPI application and python scripts. It has been shown that this approach is simple and efficient which offers the speedup of 2.7 compared to post hoc data processing. Secondly, we have developed a deep learning model for predicting multiresolution steady flow fields. The deep learning model can give reasonably accurate predictions of simulation results with orders of magnitude faster compared to simulations.
Lesur, M.*; Idomura, Yasuhiro; Shinohara, Koji; Garbet, X.*
no journal, ,
Asahi, Yuichi; Idomura, Yasuhiro; Ina, Takuya; Garbet, X.*; Grandgirard, V.*; Latu, G.*
no journal, ,
In the so-called delta-f gyrokinetic simulations, the scale separation between the equilibrium and fluctuation plasmas is assumed, and the time evolution is solved only for the fluctuation part. In contrast, in the full-f gyrokinetic simulations, both of the equilibrium and fluctuation plasmas are solved on the basis of the same first principle, where the self-consistent simulations for the equilibrium and fluctuation plasmas are possible. So far, there are a plenty number of cross-code benchmarks for delta-f gyrokinetic simulations, which helps to improve the robustness of the simulations. However, this is not the case for the full-f simulations since the complicated full-f physics makes benchmarks more difficult. In the presentation, we will show the progress of the full-f benchmarks and discuss the confronting issues.