Source author record

M. Hoelzl

M. Hoelzl 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

10works
3topics
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

10 published item(s)

preprint2025arXiv

Runaway electron avalanche and macroscopic beam formation: simulations of the DTT full power scenario

The transition of the Divertor Tokamak Test (DTT) facility from its initial commissioning phase (Day-0, plasma current $I_{p}=2$ MA) to the full power scenario ($I_{p}=5.5$ MA) introduces a critical shift in the dynamics of runaway electrons (REs) generation. While previous predictive studies of the low-current scenario indicated a robust safety margin against RE beam formation, this work reveals that the exponential scaling of the RE avalanche gain with plasma current severely narrows the safe operational window in the full power scenario. Using the non-linear magnetohydrodynamic code JOREK, we perform comprehensive 2D simulations of the current quench (CQ) phase of several disruption scenarios, systematically scanning initial RE seed currents and injected impurity levels. The results demonstrate that in the full power scenario, the avalanche multiplication factor is sufficiently high ($G_\text{av} \approx 1.3 \cdot 10^5$) to convert a mere 5.5 A seed current into macroscopic RE beams of $\approx 0.7$ MA when large amounts of impurities are present. For even higher RE seeds, the RE current can peak at $ \approx 3.2$ MA, constituting up to $\approx$ 80% of the total plasma current during the CQ. These findings suggest that, unlike the Day-0 phase, the disruption mitigation strategy for the full power scenario involves a careful balance between thermal load mitigation and RE avoidance, necessitating a well-chosen quantity of injected impurities. This work provides the baseline needed for future estimations of RE loads on the plasma-facing components of DTT, which will be essential for designing and positioning mitigation components like sacrificial limiters.

preprint2021arXiv

MHD simulations of small ELMs at low triangularity in ASDEX Upgrade

The development of small- and no-ELM regimes for ITER is a high priority topic due to the risks associated to type-I ELMs. By considering non-linear extended MHD simulations of the ASDEX Upgrade tokamak with the JOREK code, we probe a regime that avoids type-I ELMs completely provided that the separatrix density is high enough. The dynamics of the pedestal in this regime are observed to be qualitatively similar to the so-called quasi-continuous exhaust (QCE) regime in several ways. Repetitive type-I ELMs are substituted by roughly constant levels of outwards transport caused by peeling-ballooning modes (with dominant ballooning characteristics) which are localised in the last 5\% of the confined region (in normalised poloidal flux). The simulated low triangularity plasma transitions to a type-I ELMy H-mode if the separatrix density is sufficiently reduced or if the input heating power is sufficiently increased. The stabilising factors that play a role in the suppression of the small ELMs are also investigated by analysing the simulations, and the importance of including diamagnetic effects in the simulations is highlighted. By considering a scan in the pedestal resistivity and by measuring the poloidal velocity of the modes (and comparing to theoretical estimates for ideal and resistive modes), we identify the underlying instabilities as resistive peeling-ballooning modes. Decreasing the resistivity below experimentally-relevant conditions (i.e., going towards ideal MHD), the peeling-ballooning modes that constrain the pedestal below the type-I ELM stability boundary display sharply decreasing growth rates.

preprint2021arXiv

The development of an implicit full f method for electromagnetic particle simulations of Alfvén waves and energetic particle physics

In this work, an implicit scheme for particle-in-cell/Fourier electromagnetic simulations is developed and applied to studies of Alfvén waves in one dimension and three-dimensional tokamak plasmas. An analytical treatment is introduced to achieve efficient convergence of the iterative solution of the implicit field-particle system. First, its application to the one-dimensional uniform plasma demonstrates its applicability in a broad range of $β/m_e$ values. Second, toroidicity induced Alfvén eigenmodes (TAE) are simulated in a three-dimensional axisymmetric tokamak plasma, using the widely studied case defined by the International Tokamak Physics Activity (ITPA) Energetic Particle (EP) Topical Group. The real frequency and the growth (or damping) rate of the TAE with (or without) EPs agree with previous results reasonably well. The full f electromagnetic particle scheme established in this work provides a possible natural choice for EP transport studies where large profile variation and arbitrary particle distribution functions need to be treated in kinetic simulations.

preprint2020arXiv

Development and testing of an unstructured mesh method for whole plasma gyrokinetic simulations in realistic tokamak geometry

In this work, we have formulated and implemented a mixed unstructured mesh-based finite element (FE)-Fourier decomposition scheme for gyrokinetic simulations in realistic tokamak geometry. An efficient particle positioning (particle-triangle mapping) scheme for the charge deposition and field scattering using an intermediate grid as the search index for triangles has been implemented and a significant speed-up by a factor of $\sim30$ is observed as compared with the brute force scheme for a medium-size simulation. The TRIMEG (TRIangular MEsh based Gyrokinetic) code has been developed. As an application, the ion temperature gradient (ITG) mode is simulated using the simplified gyrokinetic Vlasov-Poisson model. Our simulation and that using the ORB5 code for the DIII-D Cyclone case show reasonable agreement. As an additional application, ITG simulations using an ASDEX Upgrade equilibrium have been performed with density and temperature gradient profiles similar to the Cyclone case. Capabilities of the TRIMEG code for simulations with realistic experimental equilibria in the plasma core and in the whole plasma volume with open field lines are demonstrated.

preprint2020arXiv

Understanding the reduction of the edge safety factor during hot VDEs and fast edge cooling events

In the present work a simple analytical approach is presented in order to clarify the physics behind the edge current density behaviour of a hot plasma entering in contact with a resistive conductor. When a plasma enters in contact with a highly resistive wall, large current densities appear at the edge of the plasma. The model shows that this edge current originates from the plasma response, which attempts to conserve the poloidal magnetic flux ($Ψ$) when the outer current is being lost. The loss of outer current is caused by the high resistance of the outer current path compared to the plasma core resistance. The resistance of the outer path may be given by plasma contact with a very resistive structure or by a sudden decrease of the outer plasma temperature (e.g. due to a partial thermal quench or due to a cold front penetration caused by massive gas injection). For general plasma geometries and current density profiles the model shows that given a small change of minor radius ($δa$) the plasma current is conserved to first order ($δI_p = 0 + \mathcal{O}(δa^2)$). This conservation comes from the fact that total inductance remains constant ($δL = 0$) due to an exact compensation of the change of external inductance with the change of internal inductance ($δL_\text{ext}+δL_\text{int} = 0$). As the total current is conserved and the plasma volume is reduced, the edge safety factor drops according to $q_a \propto a^2/I_p$. Finally the consistency of the resulting analytical predictions is checked with the help of free-boundary MHD simulations.

preprint2019arXiv

Axisymmetric simulations of vertical displacement events in tokamaks: A benchmark of M3D-C$^1$, NIMROD and JOREK

A benchmark exercise for the modeling of vertical displacement events (VDEs) is presented and applied to the 3D nonlinear magneto-hydrodynamic codes M3D-C$^1$, JOREK and NIMROD. The simulations are based on a vertically unstable NSTX equilibrium enclosed by an axisymmetric resistive wall with rectangular cross section. A linear dependence of the linear VDE growth rates on the resistivity of the wall is recovered for sufficiently large wall conductivity and small temperatures in the open field line region. The benchmark results show good agreement between the VDE growth rates obtained from linear NIMROD and M3D-C$^1$ simulations as well as from the linear phase of axisymmetric nonlinear JOREK, NIMROD and M3D-C$^1$ simulations. Axisymmetric nonlinear simulations of a full VDE performed with the three codes are compared and excellent agreement is found regarding plasma location and plasma currents as well as eddy and halo currents in the wall.

preprint2019arXiv

First predictive simulations for deuterium shattered pellet injection in ASDEX Upgrade

First simulations of deuterium shattered pellet injection (SPI) into an ASDEX Upgrade H-Mode plasma with the JOREK MHD code are presented. Resistivity is increased by one order of magnitude in most simulations to reduce computational costs and allow for extensive parameter scans. The effect of various physical parameters onto MHD activity and thermal quench (TQ) dynamics is studied and the influence of MHD onto ablation is shown. TQs are obtained quickly after injection in most simulations with a typical duration of 100 microseconds, which slows down at lower resistivity. Although the n=1 magnetic perturbation dominates in the simulations, toroidal harmonics up to n=10 contribute to stochastization and stochastic transport in the plasma core. The post-TQ density profile remains hollow for a few hundred microseconds. However, when flux surfaces re-form around the magnetic axis, the density has become monotonic again suggesting a beneficial behaviour for runaway electron avoidance/mitigation. With $10^{21}$ atoms injected, the TQ is typically incomplete and triggered when the shards reach the q=2 rational surface. At a larger number of injected atoms, the TQ can set in even before the shards reach this surface. For low field side injection considered here, repeated formation of outward convection cells is observed in the ablation region reducing material assimilation. Responsible is a sudden rise of pressure in the high density cloud when the stochastic region expands further releasing heat from the hot core. After the TQ, strong sheared poloidal rotation is created by Maxwell stress, which contributes to re-formation of flux surfaces.

preprint2016arXiv

Non-linear modeling of the plasma response to RMPs in ASDEX Upgrade

The plasma response to Resonant Magnetic Perturbations (RMPs) in ASDEX Upgrade is modeled with the non-linear resistive MHD code JOREK, using input profiles that match those of the experiments as closely as possible. The RMP configuration for which Edge Localized Modes are best mitigated in experiments is related to the largest edge kink response observed near the X-point in modeling. On the edge resonant surfaces $q = m/n$, the coupling between the m + 2 kink component and the m resonant component is found to induce the amplification of the resonant magnetic perturbation. The ergodicity and the 3D-displacement near the X-point induced by the resonant amplification can only partly explain the density pumpout observed in experiments.

preprint2014arXiv

Non-linear Simulations of MHD Instabilities in Tokamaks Including Eddy Current Effects and Perspectives for the Extension to Halo Currents

The dynamics of large scale plasma instabilities can strongly be influenced by the mutual interaction with currents flowing in conducting vessel structures. Especially eddy currents caused by time-varying magnetic perturbations and halo currents flowing directly from the plasma into the walls are important. The relevance of a resistive wall model is directly evident for Resistive Wall Modes (RWMs) or Vertical Displacement Events (VDEs). However, also the linear and non-linear properties of most other large-scale instabilities may be influenced significantly by the interaction with currents in conducting structures near the plasma. The understanding of halo currents arising during disruptions and VDEs, which are a serious concern for ITER as they may lead to strong asymmetric forces on vessel structures, could also benefit strongly from these non-linear modeling capabilities. Modeling the plasma dynamics and its interaction with wall currents requires solving the magneto-hydrodynamic (MHD) equations in realistic toroidal X-point geometry consistently coupled with a model for the vacuum region and the resistive conducting structures. With this in mind, the non-linear finite element MHD code JOREK has been coupled with the resistive wall code STARWALL, which allows to include the effects of eddy currents in 3D conducting structures in non-linear MHD simulations. This article summarizes the capabilities of the coupled JOREK-STARWALL system and presents benchmark results as well as first applications to non-linear simulations of RWMs, VDEs, disruptions triggered by massive gas injection, and Quiescent H-Mode. As an outlook, the perspectives for extending the model to halo currents are described.

preprint2012arXiv

Coupling JOREK and STARWALL for Non-linear Resistive-wall Simulations

The implementation of a resistive-wall extension to the non-linear MHD-code JOREK via a coupling to the vacuum-field code STARWALL is presented along with first applications and benchmark results. Also, non-linear saturation in the presence of a resistive wall is demonstrated. After completion of the ongoing verification process, this code extension will allow to perform non-linear simulations of MHD instabilities in the presence of three-dimensional resistive walls with holes for limited and X-point plasmas.