Source author record

Johan Anderson

Johan Anderson 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

23works
9topics
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

23 published item(s)

preprint2019arXiv

Elucidating plasma dynamics in Hasegawa-Wakatani turbulence by information geometry

The impact of adiabatic electrons on drift-wave turbulence, modelled by the Hasegawa-Wakatani equations, is studied using information length. Information length is a novel theoretical method for measuring distances between statistical states represented by different probability distribution functions (PDFs) along the path of a system. Specifically, the time-dependent PDFs of turbulent fluctuations for a given adiabatic index $A$ is computed. The changes in fluctuation statistics are then quantified in time by using information length. The numerical results provide time traces exhibiting intermittent plasma dynamics, and such behaviour is identified by a rapid change in the information length. The effects of $A$ are discussed.

preprint2019arXiv

Global scaling of the heat transport in fusion plasmas

A global heat flux model based on a fractional derivative of plasma pressure is proposed for the heat transport in fusion plasmas. The degree of the fractional derivative of the heat flux, $α$, is defined through the power balance analysis of the steady state. The model was used to obtain the experimental values of $α$ for a large database of the JET Carbon-wall as well as ITER Like-wall plasmas. The findings show that the average fractional degree of the heat flux over the database for electrons is $α\sim 0.8$, suggesting a global scaling between the net heating and the pressure profile in the JET plasmas. The model is expected to provide an accurate and a simple description of heat transport that can be used in transport studies of fusion plasmas.

preprint2016arXiv

Charged particle dynamics in the presence of non-Gaussian Lévy electrostatic fluctuations

Full orbit dynamics of charged particles in a $3$-dimensional helical magnetic field in the presence of $α$-stable Lévy electrostatic fluctuations and linear friction modeling collisional Coulomb drag is studied via Monte Carlo numerical simulations. The Lévy fluctuations are introduced to model the effect of non-local transport due to fractional diffusion in velocity space resulting from intermittent electrostatic turbulence. The probability distribution functions of energy, particle displacements, and Larmor radii are computed and showed to exhibit a transition from exponential decay, in the case of Gaussian fluctuations, to power law decay in the case of Lévy fluctuations. The absolute value of the power law decay exponents are linearly proportional to the Lévy index $α$. The observed anomalous non-Gaussian statistics of the particles' Larmor radii (resulting from outlier transport events) indicate that, when electrostatic turbulent fluctuations exhibit non-Gaussian Lévy statistics, gyro-averaging and guiding centre approximations might face limitations and full particle orbit effects should be taken into account.

preprint2016arXiv

Phase dependent advection-diffusion in drift wave - zonal flow turbulence

In plasma turbulence theory, due to the complexity of the system with many non-linearly interacting waves, the dynamics of the phases is often disregarded and the so-called random-phase approximation (RPA) is used assuming the existence of a Chirikov-like criterion for the onset of wave stochasticity. The dynamical amplitudes are represented as complex numbers, $ψ= ψ_r + iψ_i = ae^{iθ}$, with the amplitudes slowly varying whereas the phases are rapidly varying and, in particular, distributed uniformly over the interval $[0;2π)$. However, one could expect that the phase dynamics can play a role in the self-organisation and the formation of coherent structures. In the same manner it is also expected that the RPA falls short to take coherent interaction between phases into account. In this work therefore, we studied the role of phase dynamics and the coupling of phases between different modes on the characteristic time evolution of the turbulent. We assume a simple turbulent system where the so-called stochastic oscillator model can be employed. The idea of interpreting turbulence by stochastic oscillators. The stochastic oscillator models can be derived from radical simplifications of the nonlinear terms in the Navier-Stokes or Gyro-Kinetic equations. In this particular case we adopt the basic equation for the stochastic oscillator model with passive advection and random forcing from Ref.

preprint2016arXiv

Predator-prey model for the self-organisation of stochastic oscillators in dual populations

A predator-prey model of dual populations with stochastic oscillators is presented. A linear cross-coupling between the two populations is introduced following the coupling between the motions of a Wilberforce pendulum in two dimensions: one in the longitudinal and the other in torsional plain. Within each population a Kuramoto type competition between the phases is assumed. Thus, the synchronisation state of the whole system is controlled by these two types of competitions. The results of the numerical simulations show that by adding the linear cross-coupling interactions predator-prey oscillations between the two populations appear which results in self-regulation of the system by a transfer of synchrony between the two populations. The model represents several important features of the dynamical interplay between the drift wave and zonal flow turbulence in magnetically confined plasmas, and a novel interpretation of the coupled dynamics of drift wave-zonal flow turbulence using synchronisation of stochastic oscillator is discussed.

preprint2016arXiv

The probability density function tail of the Kardar-Parisi-Zhang equation in the strongly non-linear regime

An analytical derivation of the probability density function (PDF) tail describing the strongly correlated interface growth governed by the nonlinear Kardar-Parisi-Zhang equation is provided. The PDF tail exactly coincides with a Tracy-Widom distribution i.e. a PDF tail proportional to $\exp( - c w_2^{3/2})$, where $w_2$ is the the width of the interface. The PDF tail is computed by the instanton method in the strongly non-linear regime within the Martin-Siggia-Rose framework using a careful treatment of the non-linear interactions. In addition, the effect of spatial dimensions on the PDF tail scaling is discussed. This gives a novel approach to understand the rightmost PDF tail of the interface width distribution and the analysis suggests that there is no upper critical dimension.

preprint2016arXiv

The role of phase dynamics in a stochastic model of a passively advected scalar

Collective synchronous motion of the phases is introduced in a model for the stochastic passive advection-diffusion of a scalar with external forcing. The model for the phase coupling dynamics follows the well known Kuramoto model paradigm of limit-cycle oscillators. The natural frequencies in the Kuramoto model are assumed to obey a given scale dependence through a dispersion relation of the drift-wave form $-β\frac{k}{1+k^2}$, where $β$ is a constant representing the typical strength of the gradient. The present aim is to study the importance of collective phase dynamics on the characteristic time evolution of the fluctuation energy and the formation of coherent structures. Our results show that the assumption of a fully stochastic phase state of turbulence is more relevant for high values of $β$, where we find that the energy spectrum follows a $k^{-7/2}$ scaling. Whereas for lower $β$ there is a significant difference between a-synchronised and synchronised phase states, and one could expect the formation of coherent modulations in the latter case.

preprint2015arXiv

Statistical properties of Charney-Hasegawa-Mima zonal flows

A theoretical interpretation of numerically generated probability density functions (PDFs) of intermittent plasma transport events in unforced zonal flows is provided within the Charney-Hasegawa-Mima (CHM) model. The governing equation is solved numerically with various prescribed density gradients that are designed to produce different configurations of parallel and anti-parallel streams. Long-lasting vortices form whose flow is governed by the zonal streams. It is found that the numerically generated PDFs can be matched with analytical predictions of PDFs based on the instanton method by removing the autocorrelations from the time series. In many instances the statistics generated by the CHM dynamics relaxes to Gaussian distributions for both the electrostatic and vorticity perturbations, whereas in areas with strong nonlinear interactions it is found that the PDFs are exponentially distributed.

preprint2014arXiv

A Fractional Fokker-Planck Model for Anomalous Diffusion

In this paper we present a study of anomalous diffusion using a Fokker-Planck description with fractional velocity derivatives. The distribution functions are found using numerical means for varying degree of fractionality observing the transition from a Gaussian distribution to a Lévy distribution. The statistical properties of the distribution functions are assessed by a generalized expectation measure and entropy in terms of Tsallis statistical mechanics. We find that the ratio of the generalized entropy and expectation is increasing with decreasing fractionality towards the well known so-called sub-diffusive domain, indicating a self-organising behavior.

preprint2014arXiv

Effects of the Second Harmonic and Plasma Shaping on the Geodesic Acoustic Mode

The effects of second harmonics of the density and temperature perturbations on the linear Geodesic Acoustic Mode (GAM) frequency and non-linear generation of the GAM are investigated, using a fluid model. We show that the second harmonics contribute to the frequency through the density gradient scale length and the wave number of the GAM. In addition, the linear frequency of the GAM is generally increased by coupling to the higher harmonic.

preprint2014arXiv

Effects of the Second Harmonic on the GAM in Electron Scale Turbulence

The effects higher order harmonics have been self-consistently included in the derivation of the electron branch of the electron Geodesic Acoustic Mode (el-GAM) in an Electron-Temperature-Gradient (ETG) turbulence background. The work is based on a two-fluid model including finite $β$-effects while retaining non-adiabatic ions. In solving the linear dispersion relation, it is found that the due to the coupling to the $m=2$ mode the real frequency may be significantly altered and yield higher values.

preprint2013arXiv

Effects of the equilibrium model on impurity transport in tokamaks

Gyrokinetic simulations of ion temperature gradient mode and trapped electron mode driven impurity transport in a realistic tokamak geometry are presented and compared with results using simplified geometries. The gyrokinetic results, obtained with the GENE code in both linear and non-linear modes are compared with data and analysis for a dedicated impurity injection discharge at JET. The impact of several factors on heat and particle transport is discussed, lending special focus to tokamak geometry and rotational shear. To this end, results using s-alpha and concentric circular equilibria are compared with results with magnetic geometry from a JET experiment. To further approach experimental conditions, non-linear gyrokinetic simulations are performed with collisions and a carbon background included. The impurity peaking factors, computed by finding local density gradients corresponding to zero particle flux, are discussed. The impurity peaking factors are seen to be reduced by a factor of ~2 in realistic geometry compared with the simplified geometries, due to a reduction of the convective pinch. It is also seen that collisions reduce the peaking factor for low-Z impurities, while increasing it for high charge numbers, which is attributed to a shift in the transport spectra towards higher wavenumbers with the addition of collisions. With the addition of roto-diffusion, an overall reduction of the peaking factors is observed, but this decrease is not sufficient to explain the flat carbon profiles seen at JET.

preprint2013arXiv

High Frequency Geodesic Acoustic Modes in Electron Scale Turbulence

In this work the finite $β$-effects of an electron branch of the geodesic acoustic mode (el-GAM) driven by electron temperature gradient (ETG) modes is presented. The work is based on a fluid description of the ETG mode retaining non-adiabatic ions and the dispersion relation for el-GAMs driven non-linearly by ETG modes is derived. The ETG growth rate from the fluid model is compared to the results found from gyrokinetic simulations with good agreement. A new saturation mechanism for ETG turbulence through the interaction with el-GAMs is found, resulting in a significantly enhanced ETG turbulence saturation level compared to the mixing length estimate. It is shown that the el-GAM may be stabilized by an increase in finite $β$ as well as by increasing non-adiabaticity. The decreased GAM growth rates is due to the inclusion of the Maxwell stress.

preprint2012arXiv

Derivation and quantitative analysis of the differential self-interrogation Feynman-alpha method

A stochastic theory for a branching process in a neutron population with two energy levels is used to assess the applicability of the differential self-interrogation Feynman-alpha method by numerically estimated reaction intensities from Monte Carlo simulations. More specifically, the variance to mean or Feynman-alpha formula is applied to investigate the appearing exponentials using the numerically obtained reaction intensities.

preprint2012arXiv

High Frequency Geodesic Acoustic Modes in Electron Temperature Gradient Mode Turbulence

In this work the first demonstration of a high frequency branch of the geodesic acoustic mode (GAM) driven by electron temperature gradient (ETG) modes is presented. The work is based on a fluid description of the ETG mode retaining non-adiabatic ions and the dispersion relation for high frequency GAMs driven nonlinearly by ETG modes is derived. A new saturation mechanism for ETG turbulence through the interaction with high frequency GAMs is found, resulting in a significantly enhanced ETG turbulence saturation level compared to the mixing length estimate.

preprint2012arXiv

Non-local gyrokinetic model of linear ion-temperature-gradient modes

A theory of non-local linear ion-temperature-gradient (ITG) drift modes while retaining non-adiabatic electrons is presented, extending the previous work [S. Moradi, et al {\em Phys. Plasmas} {\bf 18}, 062106 (2011)]. A dispersion relation is derived to quantify the effects of the fractional velocity operator in the Fokker-Planck equation modified by temperature gradients and non-adiabatic electrons on the real frequency and growth rate. Solving the dispersion relation, it is shown here that as the plasma becomes more turbulent, it deviates from a Maxwellian distribution and becomes Lévy distributed. The resulting Lévy distribution of the plasma may thus significantly alter the transport. The relative effect of the fractional derivative is larger on the real frequency than on the growth rate of the ITG mode.

preprint2012arXiv

Two-point theory for the differential self-interrogation Feynman-alpha method

A Feynman-alpha formula has been derived in a two region domain pertaining the stochastic differential self-interrogation (DDSI) method and the differential die-away method (DDAA). Monte Carlo simulations have been used to assess the applicability of the variance to mean through determination of the physical reaction intensities of the physical processes in the two domains. More specifically, the branching processes of the neutrons in the two regions are described by the Chapman - Kolmogorov equation, including all reaction intensities for the various processes, that is used to derive a variance to mean relation for the process. The applicability of the Feynman-alpha or variance to mean formulae are assessed in DDSI and DDAA of spent fuel configurations.

preprint2011arXiv

Secondary instability of electromagnetic ion-temperature-gradient modes for zonal flow generation

An analytical model for zonal flow generation by toroidal ion-temperature-gradient (ITG) modes, including finite $β$ electromagnetic effects, is derived. The derivation is based on a fluid model for ions and electrons and takes into account both linear and nonlinear $β$ effects. The influence of finite plasma $β$ on the zonal flow growth rate ($γ_{ZF}$) scaling is investigated for typical tokamak plasma parameters. The results show the importance of the zonal flows close to marginal stability where $γ_{ZF}/γ_{ITG}>>1$ is obtained. In this region the parameter $γ_{ZF}/γ_{ITG}$ increases with $β$, indicating that the ITG turbulence and associated transport would decrease with $β$ at a faster rate than expected from a purely linear or quasi-linear analysis.

preprint2010arXiv

Predicting PDF tails in systems with logarithmic non-linearity

The probability density function (PDF) of flux $R$ is computed in systems with logarithmic non-linearity using a model non-linear dynamical equation. The PDF tails of the first moment flux are analytically predicted to be power law. These PDF tails are shown to be broader than a Gaussian distribution and are a manifestation of intermittency caused by short lived coherent structures (instantons).

preprint2010arXiv

Signature of a universal statistical description for drift-wave plasma turbulence

This Letter provides a theoretical interpretation of numerically generated probability density functions (PDFs) of intermittent plasma transport events. Specifically, nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence produce time series of heat flux which exhibit manifestly non-Gaussian PDFs with enhanced tails. It is demonstrated that, after the removal of autocorrelations, the numerical PDFs can be matched with predictions from a fluid theoretical setup, based on the instanton method. This result points to a universality in the modeling of intermittent stochastic process, offering predictive capability.

preprint2009arXiv

Non-perturbative statistical theory of intermittency in ITG drift wave turbulence with zonal flows

The probability distribution functions (PDFs) of momentum flux and zonal flow formation in ion-temperature-gradient (ITG) turbulence are investigated, including the effect of the shear flow on the PDFs. While ITG turbulence maintains high level of transport, this may be suppressed by shear flow. Zonal flows are also shown here to have an enhanced likelihood of generation further from marginal stability which will then regulate the ITG turbulence, which is more prominent with increased shear flow, leading to a self-regulating system.

preprint2009arXiv

Predicting PDF tails of flux in plasma sheath region

This letter provides the first prediction of the probability density function (PDF) of flux $R$ in plasma sheath sheath region in magnetic fusion devices which is characterized by dynamical equations with exponential non-linearities. By using a non-perturbative statistical theory (instantons), the PDF tails of first moment are shown to be modified Gumbel distribution which represents a frequency distribution of the extreme values of the ensemble. The non-Gaussian PDF tails that are enhanced over Gaussian predictions are the result of intermittency caused by short lived coherent structures (instantons).