Local and nonlocal parallel heat transport in general magnetic fields

A novel approach that enables the study of parallel transport in magnetized plasmas is presented. The method applies to general magnetic fields with local or nonlocal parallel closures. Temperature flattening in magnetic islands is accurately computed. For a wave number $k$, the fattening time scales as $\chi_{\parallel} \tau \sim k^{-\alpha}$ where $\chi$ is the parallel diffusivity, and $\alpha=1$ ($\alpha=2$) for non-local (local) transport. The fractal structure of the devil staircase temperature radial profile in weakly chaotic fields is resolved. In fully chaotic fields, the temperature exhibits self-similar evolution of the form $T=(\chi_{\parallel} t)^{-\gamma/2} L \left[ (\chi_{\parallel} t)^{-\gamma/2} \delta \psi \right]$, where $\delta \psi$ is a radial coordinate. In the local case, $f$ is Gaussian and the scaling is sub-diffusive, $\gamma=1/2$. In the non-local case, $f$ decays algebraically, $L (\eta) \sim \eta^{-3}$, and the scaling is diffusive, $\gamma=1$

Diffusive transport and self-consistent dynamics in coupled maps

The study of diffusion in Hamiltonian systems has been a problem of interest for a number of years. In this paper we explore the influence of self-consistency on the diffusion properties of systems described by coupled symplectic maps. Self-consistency, i.e. the back-influence of the transported quantity on the velocity field of the driving flow, despite of its critical importance, is usually overlooked in the description of realistic systems, for example in plasma physics. We propose a class of self-consistent models consisting of an ensemble of maps globally coupled through a mean field. Depending on the kind of coupling, two different general types of self-consistent maps are considered: maps coupled to the field only through the phase, and fully coupled maps, i.e. through the phase and the amplitude of the external field. The analogies and differences of the diffusion properties of these two kinds of maps are discussed in detail.Comment: 13 pages, 14 figure

Truncation effects in superdiffusive front propagation with L\'evy flights

A numerical and analytical study of the role of exponentially truncated L\'evy flights in the superdiffusive propagation of fronts in reaction-diffusion systems is presented. The study is based on a variation of the Fisher-Kolmogorov equation where the diffusion operator is replaced by a $\lambda$-truncated fractional derivative of order $\alpha$ where $1/\lambda$ is the characteristic truncation length scale. For $\lambda=0$ there is no truncation and fronts exhibit exponential acceleration and algebraic decaying tails. It is shown that for $\lambda \neq 0$ this phenomenology prevails in the intermediate asymptotic regime $(\chi t)^{1/\alpha} \ll x \ll 1/\lambda$ where $\chi$ is the diffusion constant. Outside the intermediate asymptotic regime, i.e. for $x > 1/\lambda$, the tail of the front exhibits the tempered decay $\phi \sim e^{-\lambda x}/x^{(1+\alpha)}$, the acceleration is transient, and the front velocity, $v_L$, approaches the terminal speed $v_* = (\gamma - \lambda^\alpha \chi)/\lambda$ as $t\to \infty$, where it is assumed that $\gamma > \lambda^\alpha \chi$ with $\gamma$ denoting the growth rate of the reaction kinetics. However, the convergence of this process is algebraic, $v_L \sim v_* - \alpha /(\lambda t)$, which is very slow compared to the exponential convergence observed in the diffusive (Gaussian) case. An over-truncated regime in which the characteristic truncation length scale is shorter than the length scale of the decay of the initial condition, $1/\nu$, is also identified. In this extreme regime, fronts exhibit exponential tails, $\phi \sim e^{-\nu x}$, and move at the constant velocity, $v=(\gamma - \lambda^\alpha \chi)/\nu$.Comment: Accepted for publication in Phys. Rev. E (Feb. 2009

Non-diffusive transport in plasma turbulence: a fractional diffusion approach

Numerical evidence of non-diffusive transport in three-dimensional, resistive pressure-gradient-driven plasma turbulence is presented. It is shown that the probability density function (pdf) of test particles' radial displacements is strongly non-Gaussian and exhibits algebraic decaying tails. To model these results we propose a macroscopic transport model for the pdf based on the use of fractional derivatives in space and time, that incorporate in a unified way space-time non-locality (non-Fickian transport), non-Gaussianity, and non-diffusive scaling. The fractional diffusion model reproduces the shape, and space-time scaling of the non-Gaussian pdf of turbulent transport calculations. The model also reproduces the observed super-diffusive scaling