Comparison of different interpolation operators including nonlinear subdivision schemes in the simulation of particle trajectories

In this work, we compare different interpolation operators in the context of particle tracking with an emphasis on situations involving velocity field with steep gradients. Since, in this case, most classical methods give rise to the Gibbs phenomenon (generation of oscillations near discontinuities), we present new methods for particle tracking based on subdivision schemes and especially on the Piecewise Parabolic Harmonic (PPH) scheme which has shown its advantage in image processing in presence of strong contrasts. First an analytic univariate case with a discontinuous velocity field is considered in order to highlight the effect of the Gibbs phenomenon on trajectory calculation. Theoretical results are provided. Then, we show, regardless of the interpolation method, the need to use a conservative approach when integrating a conservative problem with a velocity field deriving from a potential. Finally, the PPH scheme is applied in a more realistic case of a time-dependent potential encountered in the edge turbulence of magnetically confined plasmas, to compare the propagation of density structures (turbulence bursts) with the dynamics of test particles. This study highlights the difference between particle transport and density transport in turbulent fields.

[1]  R. Marchand Test-Particle Simulation of Space Plasmas , 2010 .

[2]  Aleksej F. Filippov,et al.  Differential Equations with Discontinuous Righthand Sides , 1988, Mathematics and Its Applications.

[3]  G. Voth,et al.  Fluid particle accelerations in fully developed turbulence , 2000, Nature.

[4]  M. Crouzeix,et al.  Analyse Numérique des équations Différentielles , 1987 .

[5]  Jacques Liandrat,et al.  On the stability of the PPH nonlinear multiresolution , 2005 .

[6]  F. Toschi,et al.  Particle trapping in three-dimensional fully developed turbulence , 2005 .

[7]  P. Yeung,et al.  Eulerian acceleration statistics as a discriminator between Lagrangian stochastic models in uniform shear flow , 2000 .

[8]  Lagrangian dynamics of drift-wave turbulence , 2010 .

[9]  C. Basdevant,et al.  La décomposition en ondelettes périodiques, un outil pour l'analyse de champs inhomogènes: théorie et algorithmes , 1989 .

[10]  M. Maxey,et al.  Methods for evaluating fluid velocities in spectral simulations of turbulence , 1989 .

[11]  K. Spatschek,et al.  NONLINEAR DRIFT-WAVE STRUCTURES AND THEIR INFLUENCE ON PARTICLE TRANSPORT , 1996 .

[12]  Sergio Plata A review on the piecewise polynomial harmonic interpolation , 2008 .

[13]  Rainer Grauer,et al.  Impact of the floating-point precision and interpolation scheme on the results of DNS of turbulence by pseudo-spectral codes , 2007, Comput. Phys. Commun..

[14]  P. Gibbon,et al.  Introduction to Plasma Physics , 2017, 2007.04783.

[15]  Y. Sarazin,et al.  Intermittent particle transport in two-dimensional edge turbulence , 1998 .

[16]  Peter S. Bernard,et al.  The effect of interpolation errors on the Lagrangian analysis of simulated turbulent channel flow , 1994 .

[17]  Charles K. Chui,et al.  An Introduction to Wavelets , 1992 .

[18]  P. Ghendrih,et al.  Control of Hamiltonian chaos as a possible tool to control anomalous transport in fusion plasmas. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[19]  Marco Pettini,et al.  Geometric approach to Hamiltonian dynamics and statistical mechanics , 2000 .

[20]  R. Grauer,et al.  Lagrangian statistics of Navier–Stokes and MHD turbulence , 2007, Journal of Plasma Physics.

[21]  James Myra,et al.  Convective transport by intermittent blob-filaments: Comparison of theory and experiment , 2011 .

[22]  X. Garbet,et al.  Theoretical analysis of the influence of external biasing on long range turbulent transport in the scrape-off layer , 2003 .

[23]  J. Rasmussen,et al.  Dispersion of ideal particles in a two-dimensional model of electrostatic turbulence , 1999 .

[24]  J. Rasmussen,et al.  Turbulent flux and the diffusion of passive tracers in electrostatic turbulence , 2003 .

[25]  S. Pope,et al.  An algorithm for tracking fluid particles in numerical simulations of homogeneous turbulence , 1988 .

[26]  Manfredi,et al.  Test-particle transport in strong electrostatic drift turbulence with finite Larmor radius effects. , 1996, Physical review letters.

[27]  Test-particle simulations of collisional impurity transport in rotating spherical tokamak plasmas , 2008 .