Source author record

Shinsuke Satake

Shinsuke Satake appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

3works
1topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

3 published item(s)

preprint2020arXiv

Benchmark of a new multi-ion-species collision operator for $δf$ Monte Carlo neoclassical simulation

A numerical method to implement a linearized Coulomb collision operator in the two-weight $δf$ Monte Carlo method for multi-ion-species neoclassical transport simulation is developed. The conservation properties and the adjointness property of the operator in the collisions between two particle species with different temperatures are verified. The linearized operator in a $δf$ Monte Carlo code is benchmarked with other two kinetic simulations, a $δf$ continuum gyrokinetic code with the same linearized collision operator and a full-f PIC code with Nanbu collision operator. The benchmark simulations of the equilibration process of plasma flow and temperature fluctuation among several particle species show very good agreement between $δf$ Monte Carlo code and the other two codes. An error in the H-theorem in the two-weight $δf$ Monte Carlo method is found, which is caused by the weight spreading phenomenon inherent in the two-weight $δf$ method. It is demonstrated that the weight averaging method serves to restoring the H-theorem without causing side effect.

preprint2019arXiv

Improved linearized model collision operator for the highly collisional regime

The linearized model collision operator for multiple species plasmas given by H. Sugama, T.-H. Watanabe, and M. Nunami [Phys.\ Plasmas {\bf 16}, 112503 (2009)] is improved to be properly applicable up to the highly collisional regime. The improved linearized model operator retains conservation laws of particles, momentum, and energy as well as it reproduces the same friction-flow relations as derived by the linearized Landau operator so that this model can be used to correctly evaluate neoclassical transport fluxes in all collisionality regimes. The adjointness relations and Boltzmann's H-theorem are exactly satisfied by the improved operator except in the case of collisions between unlike particle species with unequal temperatures where these relations and H-theorem still holds approximately because there is a large difference between the masses of the two species with significantly different temperatures. Even in the unequal-temperature case, the improved operator can also be modified so as to exactly satisfy the adjointness relations while it causes the values of the friction coefficients to deviate from those given by the Landau operator. In addition, for application to gyrokinetic simulations of turbulent transport, the improved operator is transformed into the gyrophase-averaged form with keeping the finite gyroradius effect.

preprint2015arXiv

Effects of magnetic drift tangential to magnetic surfaces on neoclassical transport in non-axisymmetric plasmas

In evaluating neoclassical transport by radially-local simulations, the magnetic drift tangential to a flux surface is usually ignored in order to keep the phase-space volume conservation. In this paper, effect of the tangential magnetic drift on the local neoclassical transport are investigated. To retain the effect of the tangential magnetic drift in the local treatment of neoclassical transport, a new local formulation for the drift kinetic simulation is developed. The compressibility of the phase-space volume caused by the tangential magnetic drift is regarded as a source term for the drift kinetic equation, which is solved by using a two-weight $δf$ Monte Carlo method for non-Hamiltonian system [G.~Hu and J.~A.~Krommes, Phys. Plasmas $\rm \textbf{1}$, 863 (1994)]. It is demonstrated that the effect of the drift is negligible for the neoclassical transport in tokamaks. In non-axisymmetric systems, however, the tangential magnetic drift substantially changes the dependence of the neoclassical transport on the radial electric field $E_{\rm r}$. The peaked behavior of the neoclassical radial fluxes around $E_{\rm r} = 0$ observed in conventional local neoclassical transport simulations is removed by taking the tangential magnetic drift into account.