Source author record

F. Jenko

F. 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

27works
6topics
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

27 published item(s)

preprint2026arXiv

Plasmoid formation via competing lower-hybrid drift and Kelvin-Helmholtz instabilities: A hybrid kinetic-gyrokinetic simulation study

We investigate the nonlinear formation of plasmoids in 2D low-beta current sheets through the interplay between the Kelvin-Helmholtz instability (KHI) and the lower-hybrid drift instability (LHDI). Using a hybrid kinetic-gyrokinetic model-based Super Simple Vlasov (ssV) code with fully kinetic ions and drift-kinetic electrons, we simulate Harris-type current sheets and velocity shear layers with strong cross-field density gradients. Our central hypothesis is that steep density gradients drive LHDI, which can grow faster than KHI and initiate an inverse cascade from kinetic to fluid scales, potentially suppressing KHI. Our simulations confirm that, in thin current sheets, LHDI develops rapidly at the sheet edges and nonlinearly merges into larger-scale magnetic islands before KHI can evolve. These LHDI-driven structures distort the velocity shear and suppress classical KH vortices. In contrast, for thicker current sheets or weaker density gradients, KHI dominates and produces the expected rolled-up vortices and associated plasmoids. These findings demonstrate that LHDI-induced turbulence can act as both a seed and a regulator of plasmoid-generating instabilities, mediating cross-scale energy transfer. This mechanism is relevant to thin boundary layers in space plasmas, such as the solar wind magnetosphere interface, and suggests that microturbulence can govern large-scale magnetic topology during collisionless reconnection.

preprint2025arXiv

Verification of a hybrid gyrokinetic model using the advanced semi-Lagrange code ssV

The super simple Vlasov (ssV) code was developed to study instabilities, turbulence, and reconnection in weakly magnetized plasmas, such as the solar wind in the dissipation range and the edge of fusion plasmas. The ssV code overcomes the limitations of standard gyrokinetic theory by using a hybrid model that incorporates fully kinetic ions and gyrokinetic electrons. This hybrid gyrokinetic model enables accurate modeling in regimes characterized by steep gradients and high-frequency dynamics. To achieve this, ssV implements a set of semi-Lagrangian numerical schemes, including Positive Flux Conservative (PFC), Flux Conservative fifth-order (FCV), FCV with Umeda limiters, and a Semi-Lagrangian Monotonicity-Preserving fifth-order scheme (SLMP5). Benchmark problems such as Landau damping, ion-acoustic waves, ion Bernstein waves, and kinetic Alfven waves were employed to evaluate the schemes. The SLMP5 scheme consistently delivered the best overall accuracy and numerical stability performance. The code also addresses well-known electromagnetic gyrokinetic simulation issues, such as the Ampere cancellation problem, using carefully chosen velocity-space resolutions and accurate integral evaluation.

preprint2021arXiv

Core transport barriers induced by fast ions in global gyrokinetic GENE simulations

A novel type of internal transport barrier (ITB) called F-ATB (fast ion-induced anomalous transport barrier) has been recently observed in state-of-the-art global gyrokinetic simulations on a properly optimized ASDEX Upgrade experiment [A. Di Siena et al. Phys. Rev. Lett. {\bf 127} 025002 (2021)]. Unlike the transport barriers previously reported in literature, the trigger mechanism for the F-ATB is a basically electrostatic wave-particle resonant interaction between supra-thermal particles - generated via ion cyclotron resonance heating (ICRH) - and ion scale plasma turbulence. This resonant effect strongly depends on the particular shape of the fast ion temperature and density profiles. Therefore, to further improve our theoretical understanding of this transport barrier, we present results exploring the parameter space and physical conditions for the F-ATB generation by performing a systematic study with global GENE simulations. Particular emphasis is given to the transport barrier width and its localization by scanning over different energetic particle temperature profiles. The latter are varied in amplitude, half-width, and radial localization of an ad-hoc Gaussian-like energetic particle logarithmic temperature gradient profile. For the reference parameters at hand, a threshold in the amplitude of the fast ion logarithmic temperature gradient is identified to trigger the transport barrier effectively.

preprint2021arXiv

Nonlocal effects in negative triangularity TCV plasmas

Global gradient driven GENE gyrokinetic simulations are used to investigate TCV plasmas with negative triangularity. Considering a limited L-mode plasma, corresponding to an experimental triangularity scan, numerical results are able to reproduce the actual transport level over a major fraction of the plasma minor radius for a plasma with $δ_{\rm LCFS}=-0.3$ and its equivalent with standard positive triangularity $δ$. For the same heat flux, a larger electron temperature gradient is sustained by $δ<0$, in turn resulting in an improved electron energy confinement. Consistently with the experiments, a reduction of the electron density fluctuations is also seen. Local flux-tube simulations are used to gauge the magnitude of nonlocal effects. Surprisingly, very little differences are found between local and global approaches for $δ>0$, while local results yield a strong overestimation of the heat fluxes when $δ<0$. Despite the high sensitivity of the turbulence level with respect to the input parameters, global effects appear to play a crucial role in the negative triangularity plasma and must be retained to reconcile simulations and experiments. Finally, a general stabilizing effect of negative triangularity, reducing fluxes and fluctuations by a factor dependent on the actual profiles, is recovered.

preprint2021arXiv

Validation of edge turbulence codes against the TCV-X21 diverted L-mode reference case

Self-consistent full-size turbulent-transport simulations of the divertor and SOL of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this work, we perform a series of diverted Ohmic L-mode discharges on the TCV tokamak, building a first-of-a-kind dataset for the validation of edge turbulence models. This dataset, referred to as TCV-X21, contains measurements from 5 diagnostic systems -- giving a total of 45 1- and 2-D comparison observables in two toroidal magnetic field directions. The dataset is used to validate three flux-driven 3D fluid-turbulence models: GBS, GRILLIX and TOKAM3X. With each model, we perform simulations of the TCV-X21 scenario, tuning the particle and power source rates to achieve a reasonable match of the upstream separatrix value of density and electron temperature. We find that the simulations match the experimental profiles for most observables at the OMP -- both in terms of profile shape and absolute magnitude -- while a poorer agreement is found towards the divertor targets. The match between simulation and experiment is seen to be sensitive to the value of the resistivity, the heat conductivities, the power injection rate and the choice of sheath boundary conditions. Additionally, despite targeting a sheath-limited regime, the discrepancy between simulations and experiment also suggests that the neutral dynamics should be included. The results of this validation show that turbulence models are able to perform simulations of existing devices and achieve reasonable agreement with experimental measurements. Where disagreement is found, the validation helps to identify how the models can be improved. By publicly releasing the experimental dataset, this work should help to guide and accelerate the development of predictive turbulence simulations of the edge and SOL.

preprint2020arXiv

Self-consistent modeling of runaway electron generation in massive gas injection scenarios in ASDEX Upgrade

We present the first successful simulation of a induced disruption in ASDEX Upgrade from massive material injection (MMI) up to established runaway electron (RE) beam, thus covering pre-thermal quench, thermal quench and current quench (CQ) of the discharge. For future high-current fusion devices such as ITER, the successful suppression of REs through MMI is of critical importance to ensure the structural integrity of the vessel. To computationally study the interplay between MMI, background plasma response, and RE generation, a toolkit based on the 1.5D transport code coupling ASTRA-STRAHL is developed. Electron runaway is described by state-of-the-art reduced kinetic models in the presence of partially ionized impurities. Applied to argon MMI in ASDEX Upgrade discharge #33108, key plasma parameters measured experimentally, such as temporal evolution of the line averaged electron density, plasma current decay rate and post-CQ RE current, are well reproduced by the simulation presented. Impurity ions are transported into the central plasma by the combined effect of neoclassical processes and additional effects prescribed inside the $q = 2$ rational surface to explain experimental time scales. Thus, a thermal collapse is induced through strong impurity radiation, giving rise to a substantial RE population as observed experimentally.

preprint2019arXiv

Classification of tokamak plasma confinement states with convolutional recurrent neural networks

During a tokamak discharge, the plasma can vary between different confinement regimes: Low (L), High (H) and, in some cases, a temporary (intermediate state), called Dithering (D). In addition, while the plasma is in H mode, Edge Localized Modes (ELMs) can occur. The automatic detection of changes between these states, and of ELMs, is important for tokamak operation. Motivated by this, and by recent developments in Deep Learning (DL), we developed and compared two methods for automatic detection of the occurrence of L-D-H transitions and ELMs, applied on data from the TCV tokamak. These methods consist in a Convolutional Neural Network (CNN) and a Convolutional Long Short Term Memory Neural Network (Conv-LSTM). We measured our results with regards to ELMs using ROC curves and Youden's score index, and regarding state detection using Cohen's Kappa Index.

preprint2018arXiv

Electromagnetic turbulence suppression by energetic particle driven modes

In recent years, a strong reduction of plasma turbulence in the presence of energetic particles has been reported in a number of magnetic confinement experiments and corresponding gyrokinetic simulations. While highly relevant to performance predictions for burning plasmas, an explanation for this primarily nonlinear effect has remained elusive so far. A thorough analysis finds that linearly marginally stable energetic particle driven modes are excited nonlinearly, depleting the energy content of the turbulence and acting as an additional catalyst for energy transfer to zonal modes (the dominant turbulence saturation channel). Respective signatures are found in a number of simulations for different JET and ASDEX Upgrade discharges with reduced transport levels attributed to energetic ion effects.

preprint2016arXiv

Hybrid simulation of a parallel collisionless shock in the Large Plasma Device

We present two-dimensional hybrid kinetic/magnetohydrodynamic simulations of planned laser-ablation experiments in the Large Plasma Device (LAPD). Our results, based on parameters which have been validated in previous experiments, show that a parallel collisionless shock can begin forming within the available space. Carbon-debris ions that stream along the magnetic-field direction with a blow-off speed of four times the Alfven velocity excite strong magnetic fluctuations, eventually transfering part of their kinetic energy to the surrounding hydrogen ions. This acceleration and compression of the background plasma creates a shock front, which satisfies the Rankine-Hugoniot conditions and can therefore propagate on its own. Furthermore, we analyze the upstream turbulence and show that it is dominated by the right-hand resonant instability.

preprint2016arXiv

Improved energy confinement with nonlinear isotope effects in magnetically confined plasmas

The efficient production of electricity from nuclear fusion in magnetically confined plasmas relies on a good confinement of the thermal energy. For more than thirty years, the observation that such confinement depends on the mass of the plasma isotope and its interaction with apparently unrelated plasma conditions has remained largely unexplained and it has become one of the main unsolved issues. By means of numerical studies based on the gyrokinetic theory, we quantitatively show how the plasma microturbulence depends on the isotope mass through nonlinear multiscale microturbulence effects involving the interplay between zonal flows, electromagnetic effects and the torque applied. This finding has crucial consequences for the design of future reactors since, in spite of the fact that they will be composed by multiple ion species, their extrapolation from present day experiments heavily relies on the knowledge obtained from a long experimental tradition based in single isotope plasmas.

preprint2016arXiv

Subproton-scale cascades in solar wind turbulence: driven hybrid-kinetic simulations

A long-lasting debate in space plasma physics concerns the nature of subproton-scale fluctuations in solar wind (SW) turbulence. Over the past decade, a series of theoretical and observational studies were presented in favor of either kinetic Alfvén wave (KAW) or whistler turbulence. Here, we investigate numerically the nature of the subproton-scale turbulent cascade for typical SW parameters by means of unprecedented high-resolution simulations of forced hybrid-kinetic turbulence in two real-space and three velocity-space dimensions. Our analysis suggests that small-scale turbulence in this model is dominated by KAWs at $β\gtrsim1$ and by magnetosonic/whistler fluctuations at lower $β$. The spectral properties of the turbulence appear to be in good agreement with theoretical predictions. A tentative interpretation of this result in terms of relative changes in the damping rates of the different waves is also presented. Overall, the results raise interesting new questions about the properties and variability of subproton-scale turbulence in the SW, including its possible dependence on the plasma $β$, and call for detailed and extensive parametric explorations of driven kinetic turbulence in three dimensions.

preprint2016arXiv

The Structure of Plasma Heating in Gyrokinetic Alfvénic Turbulence

We analyze plasma heating in weakly collisional kinetic Alfvén wave (KAW) turbulence using high resolution gyrokinetic simulations spanning the range of scales between the ion and the electron gyroradii. Real space structures that have a higher than average heating rate are shown not to be confined to current sheets. This novel result is at odds with previous studies, which use the electromagnetic work in the local electron fluid frame, i.e. $\mathbf{J} \!\cdot\! (\mathbf{E} + \mathbf{v}_e\times\mathbf{B})$, as a proxy for turbulent dissipation to argue that heating follows the intermittent spatial structure of the electric current. Furthermore, we show that electrons are dominated by parallel heating while the ions prefer the perpendicular heating route. We comment on the implications of the results presented here.

preprint2015arXiv

Gyrokinetic and kinetic particle-in-cell simulations of guide-field reconnection. I: Macroscopic effects of the electron flows

In this work, we compare gyrokinetic (GK) and fully kinetic Particle-in-Cell (PIC) simulations of magnetic reconnection in the limit of strong guide field. In particular, we analyze the limits of applicability of the GK plasma model compared to a fully kinetic description of force free current sheets for finite guide fields ($b_g$). Here we report the first part of an extended comparison, focusing on the macroscopic effects of the electron flows. For a low beta plasma ($β_i=0.01$), it is shown that both plasma models develop magnetic reconnection with similar features in the secondary magnetic islands if a sufficiently high guide field ($b_g\gtrsim 30$) is imposed in the kinetic PIC simulations. Outside of these regions, in the separatrices close to the X points, the convergence between both plasma descriptions is less restrictive ($b_g\gtrsim 5$). Kinetic PIC simulations using guide fields $b_g \lesssim 30$ reveal secondary magnetic islands with a core magnetic field and less energetic flows inside of them in comparison to the GK or kinetic PIC runs with stronger guide fields. We find that these processes are mostly due to an initial shear flow absent in the GK initialization and negligible in the kinetic PIC high guide field regime, in addition to fast outflows on the order of the ion thermal speed that violate the GK ordering. Since secondary magnetic islands appear after the reconnection peak time, a kinetic PIC/GK comparison is more accurate in the linear phase of magnetic reconnection. For a high beta plasma ($β_i=1.0$) where reconnection rates and fluctuations levels are reduced, similar processes happen in the secondary magnetic islands in the fully kinetic description, but requiring much lower guide fields ($b_g\lesssim 3$).

preprint2015arXiv

Gyrokinetic studies of core turbulence features in ASDEX Upgrade H-mode plasmas

Gyrokinetic validation studies are crucial in developing confidence in the model incorporated in numerical simulations and thus improving their predictive capabilities. As one step in this direction, we simulate an ASDEX Upgrade discharge with the GENE code, and analyze various fluctuating quantities and compare them to experimental measurements. The approach taken is the following. First, linear simulations are performed in order to determine the turbulence regime. Second, the heat fluxes in nonlinear simulations are matched to experimental fluxes by varying the logarithmic ion temperature gradient within the expected experimental error bars. Finally, the dependence of various quantities with respect to the ion temperature gradient is analyzed in detail. It is found that density and temperature fluctuations can vary significantly with small changes in this parameter, thus making comparisons with experiments very sensitive to uncertainties in the experimental profiles. However, cross-phases are more robust, indicating that they are better observables for comparisons between gyrokinetic simulations and experimental measurements.

preprint2015arXiv

Multiscale nature of the dissipation range in gyrokinetic simulations of Alfvénic turbulence

Nonlinear energy transfer and dissipation in Alfvén wave turbulence are analyzed in the first gyrokinetic simulation spanning all scales from the tail of the MHD range to the electron gyroradius scale. For typical solar wind parameters at 1 AU, about 30% of the nonlinear energy transfer close to the electron gyroradius scale is mediated by modes in the tail of the MHD cascade. Collisional dissipation occurs across the entire kinetic range $k_\perpρ_i\gtrsim 1$. Both mechanisms thus act on multiple coupled scales, which have to be retained for a comprehensive picture of the dissipation range in Alfvénic turbulence.

preprint2015arXiv

Understanding nonlinear saturation in zonal-flow-dominated ion temperature gradient turbulence

We propose a quantitative model of ion temperature gradient driven turbulence in toroidal magnetized plasmas. In this model, the turbulence is regulated by zonal flows, i.e. mode saturation occurs by a zonal-flow-mediated energy cascade ("shearing"), and zonal flow amplitude is controlled by nonlinear decay. Our model is tested in detail against numerical simulations to confirm that both its assumptions and predictions are satisfied. Key results include (1) a sensitivity of the nonlinear zonal flow response to the energy content of the linear instability, (2) a persistence of zonal-flow-regulated saturation at high temperature gradients, (3) a physical explanation of the nonlinear saturation process in terms of secondary and tertiary instabilities, and (4) dependence of heat flux in terms of dimensionless parameters.

preprint2014arXiv

Collision-dependent power law scalings in 2D gyrokinetic turbulence

Nonlinear gyrokinetics provides a suitable framework to describe short-wavelength turbulence in magnetized laboratory and astrophysical plasmas. In the electrostatic limit, this system is known to exhibit a free energy cascade towards small scales in (perpendicular) real and/or velocity space. The dissipation of free energy is always due to collisions (no matter how weak the collisionality), but may be spread out across a wide range of scales. Here, we focus on freely-decaying 2D electrostatic turbulence on sub-ion-gyroradius scales. An existing scaling theory for the turbulent cascade in the weakly collisional limit is generalized to the moderately collisional regime. In this context, non-universal power law scalings due to multiscale dissipation are predicted, and this prediction is confirmed by means of direct numerical simulations.

preprint2014arXiv

Electromagnetic stabilization of tokamak microturbulence in a high-$β$ regime

The impact of electromagnetic stabilization and flow shear stabilization on ITG turbulence is investigated. Analysis of a low-$β$ JET L-mode discharge illustrates the relation between ITG stabilization, and proximity to the electromagnetic instability threshold. This threshold is reduced by suprathermal pressure gradients, highlighting the effectiveness of fast ions in ITG stabilization. Extensive linear and nonlinear gyrokinetic simulations are then carried out for the high-$β$ JET hybrid discharge 75225, at two separate locations at inner and outer radii. It is found that at the inner radius, nonlinear electromagnetic stabilization is dominant, and is critical for achieving simulated heat fluxes in agreement with the experiment. The enhancement of this effect by suprathermal pressure also remains significant. It is also found that flow shear stabilization is not effective at the inner radii. However, at outer radii the situation is reversed. Electromagnetic stabilization is negligible while the flow shear stabilization is significant. These results constitute the high-$β$ generalization of comparable observations found at low-$β$ at JET. This is encouraging for the extrapolation of electromagnetic ITG stabilization to future devices. An estimation of the impact of this effect on the ITER hybrid scenario leads to a 20% fusion power improvement.

preprint2014arXiv

Pressure tensor in the presence of velocity shear: stationary solutions and self-consistent equilibria

Observations and numerical simulations of laboratory and space plasmas in almost collisionless regimes reveal anisotropic and non-gyrotropic particle distribution functions. We investigate how such states can persist in the presence of a sheared flow. We focus our attention on the pressure tensor equation in a magnetized plasma and derive analytical self-consistent plasma equilibria which exhibit a novel asymmetry with respect to the magnetic field direction. These results are relevant for investigating, within fluid models that retain the full pressure tensor dynamics, plasma configurations where a background shear flow is present.

preprint2014arXiv

Towards optimal explicit time-stepping schemes for the gyrokinetic equations

The nonlinear gyrokinetic equations describe plasma turbulence in laboratory and astrophysical plasmas. To solve these equations, massively parallel codes have been developed and run on present-day supercomputers. This paper describes measures to improve the efficiency of such computations, thereby making them more realistic. Explicit Runge-Kutta schemes are considered to be well suited for time-stepping. Although the numerical algorithms are often highly optimized, performance can still be improved by a suitable choice of the time-stepping scheme, based on spectral analysis of the underlying operator. Here, an operator splitting technique is introduced to combine first-order Runge-Kutta-Chebychev schemes for the collision term with fourth-order schemes for the remaining terms. In the nonlinear regime, based on the observation of eigenvalue shifts due to the (generalized) $E\times B$ advection term, an accurate and robust estimate for the nonlinear timestep is developed. The presented techniques can reduce simulation times by factors of up to three in realistic cases. This substantial speedup encourages the use of similar timestep optimized explicit schemes not only for the gyrokinetic equation, but also for other applications with comparable properties.

preprint2013arXiv

Applications of Large Eddy Simulation methods to gyrokinetic turbulence

The Large Eddy Simulation (LES) approach - solving numerically the large scales of a turbulent system and accounting for the small-scale influence through a model - is applied to nonlinear gyrokinetic systems that are driven by a number of different microinstabilities. Comparisons between modeled, lower resolution, and higher resolution simulations are performed for an experimental measurable quantity, the electron density fluctuation spectrum. Moreover, the validation and applicability of LES is demonstrated through a series of diagnostics based on the free energetics of the system.

preprint2013arXiv

Ion temperature profile stiffness: non-linear gyrokinetic simulations and comparison with experiment

Recent experimental observations at JET show evidence of reduced ion temperature profile stiffness, hypothesised to be due to concomitant low magnetic shear (s) and significant toroidal rotational flow shear. Non-linear gyrokinetic simulations are performed, aiming to investigate the physical mechanism behind the observations. A comprehensive set of simulations are carried out, comparing the impact on the ion heat flux of various parameters that differ within the data-set. These parameters include q, s, rotation, effect of rotation on the magnetohydrodynamic (MHD) equilibrium, R/L_n, beta_e, Z_eff, and the fast particle content. The effect of toroidal flow shear itself is not predicted by the simulations to lead to a significant reduction in ion heat flux, due both to an insufficient magnitude of flow shear and significant parallel velocity gradient destabilisation. It is however found that non-linear electromagnetic effects due to both thermal and fast-particle pressure gradients, even at low beta_e, can significantly reduce the profile stiffness. A total of five discharges are examined, at both inner and outer radii. For all cases studied, the simulated and experimental ion heat flux values agree within reasonable variations of input parameters around the experimental uncertainties.

preprint2013arXiv

Nonlinear stabilization of tokamak microturbulence by fast ions

Nonlinear electromagnetic stabilization by suprathermal pressure gradients found in specific regimes is shown to be a key factor in reducing tokamak microturbulence, augmenting significantly the thermal pressure electromagnetic stabilization. Based on nonlinear gyrokinetic simulations investigating a set of ion heat transport experiments on the JET tokamak, described by Mantica et al. [Phys. Rev. Lett. 107 135004 (2011)], this result explains the experimentally observed ion heat flux and stiffness reduction. These findings are expected to improve the extrapolation of advanced tokamak scenarios to reactor relevant regimes.

preprint2011arXiv

Gyrokinetic Large Eddy Simulations

The Large Eddy Simulation (LES) approach is adapted to the study of plasma microturbulence in a fully three-dimensional gyrokinetic system. Ion temperature gradient driven turbulence is studied with the {\sc GENE} code for both a standard resolution and a reduced resolution with a model for the sub-grid scale turbulence. A simple dissipative model for representing the effect of the sub-grid scales on the resolved scales is proposed and tested. Once calibrated, the model appears to be able to reproduce most of the features of the free energy spectra for various values of the ion temperature gradient.

preprint2010arXiv

Free energy cascade in gyrokinetic turbulence

In gyrokinetic theory, the quadratic nonlinearity is known to play an important role in the dynamics by redistributing (in a conservative fashion) the free energy between the various active scales. In the present study, the free energy transfer is analyzed for the case of ion temperature gradient driven turbulence. It is shown that it shares many properties with the energy transfer in fluid turbulence. In particular, one finds a forward (from large to small scales), extremely local, and self-similar cascade of free energy in the plane perpendicular to the background magnetic field. These findings shed light on some fundamental properties of plasma turbulence, and encourage the development of large eddy simulation techniques for gyrokinetics.

preprint2009arXiv

Direct multiscale coupling of a transport code to gyrokinetic turbulence codes

Direct coupling between a transport solver and local, nonlinear gyrokinetic calculations using the multiscale gyrokinetic code TRINITY [M. Barnes, Ph.D. thesis, arxiv:0901.2868] is described. The coupling of the microscopic and macroscopic physics is done within the framework of multiscale gyrokinetic theory, of which we present the assumptions and key results. An assumption of scale separation in space and time allows for the simulation of turbulence in small regions of the space-time grid, which are embedded in a coarse grid on which the transport equations are implicitly evolved. This leads to a reduction in computational expense of several orders of magnitude, making first-principles simulations of the full fusion device volume over the confinement time feasible on current computing resources. Numerical results from TRINITY simulations are presented and compared with experimental data from JET and ASDEX Upgrade plasmas.

preprint2009arXiv

Wave-driven dynamo action in spherical MHD systems

Hydrodynamic and magnetohydrodynamic numerical studies of a mechanically forced two-vortex flow inside a sphere are reported. The simulations are performed in the intermediate regime between the laminar flow and developed turbulence where a hydrodynamic instability is found to generate internal waves with a characteristic m=2 zonal wave number. It is shown that this time-periodic flow acts as a dynamo although snapshots of the flow as well as the mean flow are not dynamos. The magnetic fields' growth rate exhibits resonance effects depending on the wave frequency. Furthermore, a cyclic self-killing and self-recovering dynamo based on the relative alignment of the velocity and magnetic fields is presented. The phenomena are explained in terms of a mixing of non-orthogonal eigenstates of the time dependent linear operator of the magnetic induction equation. The potential relevance of this mechanism to dynamo experiments is discussed.