Source author record

Frank Jenko

Frank Jenko 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

21works
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

21 published item(s)

preprint2023arXiv

Zonal Flow Excitation in Electron-Scale Tokamak Turbulence

The derivation of an intermediate-scale gyrokinetic-electron theory in nonuniform tokamak plasmas [Chen H. et al 2021 Nucl. Fusion 61 066017] has shown that a Navier-Stokes type nonlinearity couples electron-temperature-gradient (ETG) modes and zonal flow (ZF) modes with wavelengths much shorter than the ion gyroradius but much longer than the electron gyroradius. This intermediate-scale ETG-ZF coupling is typically stronger than the Hasegawa-Mima type nonlinearity characteristic of the fluid approximation and is predicted to lead to relevant zonal flow generation and ETG mode regulation. Electron-scale, continuum, gyrokinetic simulation results are presented here which include both single-mode ETG and full-spectrum ETG turbulence. The zonal flow generation due to single ETG modes is investigated and the single-mode intermediate-scale results are found to be in agreement with theory. The full-spectrum results are then presented and explained qualitatively in terms of the single-mode results. It is found that the ETG-driven zonal flows regulate intermediate-scale electron heat flux transport to levels in the predicted range.

preprint2022arXiv

Gyrokinetic modelling of anisotropic energetic particle driven instabilities in tokamak plasmas

Energetic particles produced by neutral beams are observed to excite energetic-particle-driven geodesic acoustic modes (EGAMs) in tokamaks. We study the effects of anisotropy of distribution function of the energetic particles on the excitation of such instabilities with ORB5, a gyrokinetic particle-in-cell code. Numerical results are shown for linear electrostatic simulations with ORB5. The growth rate is found to be sensitively dependent on the phase-space shape of the distribution function. The behavior of the instability is qualitatively compared to the theoretical analysis of dispersion relations. Realistic neutral beam energetic particle anisotropic distributions are obtained from the heating solver RABBIT and are introduced into ORB5 as input distribution function. Results show a dependence of the growth rate on the injection angle. A qualitative comparison to experimental measurements is presented and few disagreements between them are found, being the growth rate in the simulations much lower than that from experiments. An explanation for the difference is advanced.

preprint2020arXiv

Gyrokinetic full-f particle-in-cell simulations on open field lines with PICLS

While in recent years, gyrokinetic simulations have become the workhorse for theoretical turbulence and transport studies in the plasma core, their application to the edge and scrape-off layer (SOL) region presents significant challenges. In particular, steep density and temperature gradients as well as large fluctuation amplitudes call for a "full-f" treatment. To specifically study problems in the SOL region, the gyrokinetic particle-in-cell (PIC) code PICLS has been developed. The code is based on an electrostatic full-f model with linearised field equations and uses kinetic electrons. Here, the well-studied parallel transport problem during an edge-localized mode (ELM) in the SOL shall be investigated for one spatial dimension. The results are compared to previous gyrokinetic continuum and fully kinetic PIC simulations and show good agreement.

preprint2020arXiv

Nonlinear symmetry breaking in electron temperature gradient driven turbulence

Nonlinear symmetry breaking may occur in systems with two or more states whose linear dynamics displays certain symmetries, one of which is preferred nonlinearly. We have identified a regime of electron temperature gradient (ETG) instabilities in a tokamak plasma with circular concentric flux surfaces that has its largest growth rate at a finite ballooning angle, establishing a symmetry that is nonlinearly broken to favor one sign for the ballooning angle. This is the first example of nonlinear symmetry breaking in simulations of a drift instability in the absence of externally imposed flow shear or asymmetry in the plasma column.

preprint2020arXiv

Turbulence suppression by energetic particle effects in modern optimized stellarators

Turbulent transport is known to limit the plasma confinement of present-day optimized stellarators. To address this issue, a novel method to strongly suppress turbulence in such devices is proposed, namely the resonant wave-particle interaction of supra-thermal particles - e.g., from ion-cyclotron-resonance-frequency (ICRF) heating - with turbulence-driving microinstabilities like ion-temperature-gradient (ITG) modes. The effectiveness of this mechanism is demonstrated via analytic theory and large-scale gyrokinetic simulations, revealing an overall turbulence reduction by up to $80\%$. These results hold the promise of new and still unexplored stellarator scenarios with enhanced confinement and improved performance, essential for achieving burning plasmas in future devices.

preprint2019arXiv

Collisional gyrokinetic full-f particle-in-cell simulations on open field lines with PICLS

Applying gyrokinetic simulations for theoretical turbulence and transport studies to the plasma edge and scrape-off layer (SOL) presents significant challenges. To in particular account for steep density and temperature gradients in the SOL, the "full-f" code PICLS was developed. PICLS is a gyrokinetic particle-in-cell (PIC) code and is based on an electrostatic model with a linearized field equation and uses kinetic electrons. In previously published results we were applying PICLS to the well-studied 1D parallel transport problem during an edge-localized mode (ELM) in the SOL without collisions. As an extension to this collision-less case and in preparation for 3D simulations, in this work a collisional model will be introduced. The implemented Lenard-Bernstein collision operator and its Langevin discretization will be shown. Conservation properties of the collision operator as well as a comparison of the collisional and non-collisional case will be discussed.

preprint2019arXiv

Treatment of Advanced Divertor Configurations in the Flux-Coordinate Independent turbulence code GRILLIX

Advanced divertor configurations modify the magnetic geometry of the diverter to achieve a combination of strong magnetic flux expansion, increased connection length and higher divertor volume - to improve detachment stability, neutral/impurity confinement and heat-channel broadening. In this paper, we discuss the modification of the Flux-Coordinate Independent (FCI) turbulence code GRILLIX to treat generalised magnetic geometry, to allow for the investigation of the effect of magnetic geometry on turbulent structures in the edge and SOL. The development of grids and parallel operators from numerically-defined magnetic equilibria is discussed, as is the application of boundary conditions via penalisation, with the finite-width method generalised to treat complex non-conformal boundaries. Initial testing of hyperbolic (advection) and parabolic (diffusion) test cases is presented for the Snowflake scenario.

preprint2016arXiv

A linear dispersion relation for the hybrid kinetic-ion/fluid-electron model of plasma physics

A dispersion relation for a commonly used hybrid model of plasma physics is developed, which combines fully kinetic ions and a massless-electron fluid description. Although this model and variations of it have been used to describe plasma phenomena for about 40 years, to date there exists no general dispersion relation to describe the linear wave physics contained in the model. Previous efforts along these lines are extended here to retain arbitrary wave propagation angles, temperature anisotropy effects, as well as additional terms in the generalized Ohm's law which determines the electric field. A numerical solver for the dispersion relation is developed, and linear wave physics is benchmarked against solutions of a full Vlasov-Maxwell dispersion relation solver. This work opens the door to a more accurate interpretation of existing and future wave and turbulence simulations using this type of hybrid model.

preprint2016arXiv

Comparative study of gyrokinetic, hybrid-kinetic and fully kinetic wave physics for space plasmas

A set of numerical solvers for the linear dispersion relations of the gyrokinetic, the hybrid-kinetic, and the fully kinetic model is employed to study the physics of the kinetic Alfvén wave and the fast magnetosonic mode in these models. In particular, we focus on parameters that are relevant for solar wind oriented applications (using a homogeneous, isotropic background), which are characterized by wave propagation angles averaging close to 90°. It is found that the gyrokinetic model, while lacking high-frequency solutions and cyclotron effects, faithfully reproduces the fully kinetic Alfvén wave physics close to, and sometimes significantly beyond, the boundaries of its range of validity. The hybrid-kinetic model, on the other hand, is much more complete in terms of high-frequency waves, but owing to its simple electron model it is found to severely underpredict wave damping rates even on ion spatial scales across a large range of parameters, despite containing full kinetic ion physics.

preprint2016arXiv

Numerical simulations of the Princeton magneto-rotational instability experiment with conducting axial boundaries

We investigate numerically the Princeton magneto-rotational instability (MRI) experiment and the effect of conducting axial boundaries or endcaps. MRI is identified and found to reach a much higher saturation than for insulating endcaps. This is probably due to stronger driving of the base flow by the magnetically rather than viscously coupled boundaries. Although the computations are necessarily limited to lower Reynolds numbers ($\Re$) than their experimental counterparts, it appears that the saturation level becomes independent of $\Re$ when $\Re$ is sufficiently large, whereas it has been found previously to decrease roughly as $\Re^{-1/4}$ with insulating endcaps. The much higher saturation levels will allow for the first positive detection of MRI beyond its theoretical and numerical predictions.

preprint2015arXiv

Cosmic-ray pitch-angle scattering in imbalanced MHD turbulence simulations

Pitch-angle scattering rates for cosmic-ray particles in magnetohydrodynamic (MHD) simulations with imbalanced turbulence are calculated for fully evolving electromagnetic turbulence. We compare with theoretical predictions derived from the quasilinear theory of cosmic-ray diffusion for an idealized slab spectrum and demonstrate how cross helicity affects the shape of the pitch-angle diffusion coefficient. Additional simulations in evolving magnetic fields or static field configurations provide evidence that the scattering anisotropy in imbalanced turbulence is not primarily due to coherence with propagating Alfven waves, but an effect of the spatial structure of electric fields in cross-helical MHD turbulence.

preprint2015arXiv

How turbulence regulates biodiversity in systems with cyclic competition

Cyclic, nonhierarchical interactions among biological species represent a general mechanism by which ecosystems are able to maintain high levels of biodiversity. However, species coexistence is often possible only in spatially extended systems with a limited range of dispersal, whereas in well-mixed environments models for cyclic competition often lead to a loss of biodiversity. Here we consider the dispersal of biological species in a fluid environment, where mixing is achieved by a combination of advection and diffusion. In particular, we perform a detailed numerical analysis of a model composed of turbulent advection, diffusive transport, and cyclic interactions among biological species in two spatial dimensions and discuss the circumstances under which biodiversity is maintained when external environmental conditions, such as resource supply, are uniform in space. Cyclic interactions are represented by a model with three competitors, resembling the children's game of rock-paper-scissors, whereas the flow field is obtained from a direct numerical simulation of two-dimensional turbulence with hyperviscosity. It is shown that the space-averaged dynamics undergoes bifurcations as the relative strengths of advection and diffusion compared to biological interactions are varied.

preprint2015arXiv

Large-Eddy Simulations of Magnetohydrodynamic Turbulence in Heliophysics and Astrophysics

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions.

preprint2015arXiv

The anisotropic redistribution of free energy for gyrokinetic plasma turbulence in a Z-pinch

For a Z-pinch geometry, we report on the nonlinear redistribution of free energy across scales perpendicular to the magnetic guide field, for a turbulent plasma described in the framework of gyrokinetics. The analysis is performed using a local flux-surface approximation, in a regime dominated by electrostatic fluctuations driven by the entropy mode, with both ion and electron species being treated kinetically. To explore the anisotropic nature of the free energy redistribution caused by the emergence of zonal flows, we use a polar coordinate representation for the field-perpendicular directions and define an angular density for the scale flux. Positive values for the classically defined (angle integrated) scale flux, which denote a direct energy cascade, are shown to be also composed of negative angular sections, a fact that impacts our understanding of the backscatter of energy and the way in which it enters the modeling of sub-grid scales for turbulence. A definition for the flux of free energy across each perpendicular direction is introduced as well, which shows that the redistribution of energy in the presence of zonal flows is highly anisotropic.

preprint2014arXiv

Acceleration of Particles in Imbalanced Magnetohydrodynamic Turbulence

The present work investigates the acceleration of test particles in balanced and imbalanced Alfvénic turbulence, relevant to the solar-wind problem. These turbulent states, obtained numerically by prescribing the injection rates for the ideal invariants, are evolved dynamically with the particles. While the energy spectrum for balanced and imbalanced states is known, the impact made on particle heating is a matter of debate, with different considerations giving different results. By performing direct numerical simulations, resonant and non-resonant particle accelerations are automatically considered and the correct turbulent phases are taken into account. For imbalanced turbulence, it is found that the acceleration rate of charged particles is reduced and the heating rate diminished. This behaviour is independent of the particle gyroradius, although particles that have a stronger adiabatic motion (smaller gyroradius) tend to experience a larger heating.

preprint2014arXiv

The energetic coupling of scales in gyrokinetic plasma turbulence

In magnetized plasma turbulence, the couplings of perpendicular spatial scales that arise due to the nonlinear interactions are analyzed from the perspective of the free-energy exchanges. The plasmas considered here, with appropriate ion or electron adiabatic electro-neutrality responses, are described by the gyrokinetic formalism in a toroidal magnetic geometry. Turbulence develops due to the electrostatic fluctuations driven by temperature gradient instabilities, either ion temperature gradient (ITG) or electron temperature gradient (ETG). The analysis consists in decomposing the system into a series of scale structures, while accounting separately for contributions made by modes possessing special symmetries (e.g., the zonal flow modes). The interaction of these scales is analyzed using the energy transfer functions, including a forward and backward decomposition, scale fluxes and locality functions. The comparison between the ITG and ETG cases shows that ETG turbulence has a more pronounced classical turbulent behavior, exhibiting a stronger energy cascade, with implications for gyrokinetic turbulence modeling.

preprint2013arXiv

Identification of vortexes obstructing the dynamo mechanism in laboratory experiments

The magnetohydrodynamic dynamo effect explains the generation of self-sustained magnetic fields in electrically conducting flows, especially in geo- and astrophysical environments. Yet the details of this mechanism are still unknown, e.g., how and to which extent the geometry, the fluid topology, the forcing mechanism and the turbulence can have a negative effect on this process. We report on numerical simulations carried out in spherical geometry, analyzing the predicted velocity flow with the so-called Singular Value Decomposition, a powerful technique that allows us to precisely identify vortexes in the flow which would be difficult to characterize with conventional spectral methods. We then quantify the contribution of these vortexes to the growth rate of the magnetic energy in the system. We identify an axisymmetric vortex, whose rotational direction changes periodically in time, and whose dynamics are decoupled from those of the large scale background flow, is detrimental for the dynamo effect. A comparison with experiments is carried out, showing that similar dynamics were observed in cylindrical geometry. These previously unexpected eddies, which impede the dynamo effect, offer an explanation for the experimental difficulties in attaining a dynamo in spherical geometry.

preprint2013arXiv

Nonuniversal power-law spectra in turbulent systems

Turbulence is generally associated with universal power-law spectra in scale ranges without significant drive or damping. Although many examples of turbulent systems do not exhibit such an inertial range, power-law spectra may still be observed. As a simple model for such situations, a modified version of the Kuramoto-Sivashinsky equation is studied. By means of semi-analytical and numerical studies, one finds power laws with nonuniversal exponents in the spectral range for which the ratio of nonlinear and linear time scales is (roughly) scale-independent.

preprint2012arXiv

Locality properties of the free energy fluxes in gyrokinetic turbulence

The nature of the nonlinear interactions in gyrokinetic (GK) turbulence, driven by an ion-temperature gradient instability, is investigated using numerical simulations of single ion species plasma in three-dimensional flux tube geometry. To account for the level of separation existing between scales involved in an energetic interaction, the degree of locality of the free energy scale flux is analyzed employing Kraichnan's infrared (IR) and ultraviolet locality functions. Due to the nontrivial dissipative nature of GK turbulence, an asymptotic level for the locality exponents, indicative of a universal dynamical regime for GK's, is not recovered and an accentuated non-local behavior of the IR interactions is found instead, in spite of the local energy cascade observed.

preprint2011arXiv

Dynamic Procedure for Filtered Gyrokinetic Simulations

Large Eddy Simulations (LES) of gyrokinetic plasma turbulence are investigated as interesting candidates to decrease the computational cost. A dynamic procedure is implemented in the GENE code, allowing for dynamic optimization of the free parameters of the LES models (setting the amplitudes of dissipative terms). Employing such LES methods, one recovers the free energy and heat flux spectra obtained from highly resolved Direct Numerical Simulations (DNS). Systematic comparisons are performed for different values of the temperature gradient and magnetic shear, parameters which are of prime importance in Ion Temperature Gradient (ITG) driven turbulence. Moreover, the degree of anisotropy of the problem, that can vary with parameters, can be adapted dynamically by the method that shows Gyrokinetic Large Eddy Simulation (GyroLES) to be a serious candidate to reduce numerical cost of gyrokinetic solvers.

preprint2010arXiv

Impurity transport in ITG and TE mode dominated turbulence

In the present study, turbulent impurity transport in Deuterium tokamak plasmas, driven by Ion Temperature Gradient (ITG) and Trapped Electron (TE) modes, has been investigated using fluid and gyrokinetic models. The impurity diffusivity ($D_Z$) and convective velocity ($V_Z$) are calculated, and from these the zero-flux peaking factor ($PF_0$) is derived. This quantity expresses the impurity density gradient at which the convective and diffusive transport of impurities are exactly balanced. The sign of $PF_0$ is of special interest, as it determines whether the impurities are subject to an inward ($PF_0 > 0$) or outward ($PF_0 < 0$) pinch.