トカマクプラズマにおける二次元輸送モデリング

京都大学0048新制・課程博士博士(工学)甲第18275号工博第3867号新制||工||1593(附属図書館)31133京都大学大学院工学研究科原子核工学専攻(主査)教授 福山 淳, 教授 功刀 資彰, 准教授 村上 定義学位規則第4条第1項該当Doctor of Philosophy (Engineering)Kyoto UniversityDFA

Drift reduced Landau fluid model for magnetized plasma turbulence simulations in BOUT++ framework

Recently the drift-reduced Landau fluid six-field turbulence model within the BOUT++ framework has been upgraded. In particular, this new model employs a new normalization, adds a volumetric flux-driven source option, the Landau fluid closure for parallel heat flux and a Laplacian inversion solver which is able to capture n=0 axisymmetric mode evolution in realistic tokamak configurations. These improvements substantially extended model's capability to study a wider range of tokamak edge phenomena, and are essential to build a fully self-consistent edge turbulence model capable of both transient (e.g., ELM, disruption) and transport time-scale simulations.Comment: 26 pages, 14 figure

Impact of equilibrium radial electric field on energy loss process after pedestal collapse

An impact of the equilibrium radial electric field on energy loss processes after pedestal collapse is numerically investigated using the BOUT++ framework. Using linear stability analysis, the resistive ballooning mode is shown to be stabilized by the radial shear of the equilibrium radial electric field. On the other hand, the energy loss level after the pedestal collapse increases if the equilibrium radial electric field is taken into account. The spatio‐temporal and phase diagram analyses show that the equilibrium radial electric field partially cancels the fluctuation‐driven toroidally axisymmetric radial electric field and weakens the E  × B shearing rate after pedestal collapse, weakening the turbulence suppression by vortex shearing. The equilibrium radial electric field therefore increases turbulence intensity in nonlinear cyclic oscillations among pressure gradient, E  × B shearing rate, and turbulence intensity, which gives rise to subsequent bursts of turbulent transport and increases the energy loss level

Formulation of Two-Dimensional Transport Modeling in Tokamak Plasmas

A two-dimensional transport modeling applicable to a whole tokamak plasma is proposed. The model is derived from the multi-fluid equations and Maxwell's equations and the moment approach of neoclassical transport is employed as fluid closures. The multi-fluid equations consist of the equations for particle density, momentum, energy and total heat flux transport for each plasma species. The expressions of the parallel viscosity and heat viscosity are extended in order to be applicable to both inside and outside of the last closed flux surface. It is confirmed that our neoclassical transport model is consistent with the ordinary flux-surface-averaged one-dimensional neoclassical transport model. Our transport equations are coupled with the electromagnetic equations in order to describe the time evolution of tokamak plasmas. The procedure for coupling a transport solver based on our transport model with an equilibrium solver is also briefly described