Researcher profile

Matthias Hoelzl

Matthias Hoelzl contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2022arXiv

Confinement of passing and trapped runaway electrons in the simulation of an ITER current quench

Runaway electrons (REs) present a high-priority issue for ITER but little is known about the extent to which RE generation is affected by the stochastic field intrinsic to disrupting plasmas. RE generation can be modelled with reduced kinetic models and there has been recent progress in involving losses due to field stochasticity, either via a loss-time parameter or radial transport coefficients which can be estimated by tracing test electrons in 3D fields. We evaluate these terms in ITER using a recent JOREK 3D MHD simulation of plasma disruption to provide the stochastic magnetic fields where RE markers are traced with the built-in particle tracing module. While the MHD simulation modelled only the current quench phase, the case is MHD unstable and exhibits similar relaxation as would be expected during the thermal quench. Therefore, the RE simulations can be considered beginning right after the thermal quench but before the MHD relaxation is complete. The plasma is found to become fully stochastic for 8 ms and the resulting transport is sufficient to overcome RE avalanche before flux surfaces are reformed. We also study transport mechanisms for trapped REs and find those to be deconfined as well during this phase. While the results presented here are not sufficient to assess the magnitude of the formed RE beam, we show that significant RE losses could be expected to arise due to field stochasticity.

preprint2022arXiv

Modelling of saturated external MHD instabilities in tokamaks: a comparison of 3D free boundary equilibria and nonlinear stability calculations

3D free boundary equilibrium computations have recently been used to model external kinks and edge harmonic oscillations (EHOs), comparing with linear MHD stability codes, and nonlinear analytic theory [Kleiner et al, PPCF 61 084005 (2019)]. In this study, results of the VMEC equilibrium code are compared further with nonlinear reduced MHD simulations, using the JOREK code. The purpose of this investigation is to understand the extent to which the modelling approaches agree, and identify the important physical effects which can modify the dynamics. For the simulated external kink, which is dominated by a single toroidal harmonic, good agreement is found when a large Lundquist number is used in the JOREK simulation, such that resistive effects are sub-dominant. Modelling EHOs where multiple toroidal harmonics are linearly unstable, the saturated perturbation observed can differ in the dominant toroidal harmonic. On the ideal timescale, a n=2 EHO is observed in JOREK, while the saturated perturbation predicted by VMEC is a n=1 mode. Extending simulations into timescales where resistive effects can play a role, similar n=1 perturbations can be found. The coupling of different linearly unstable toroidal harmonics in the JOREK simulation broadens the magnetic energy spectrum and ergodises the plasma edge region, resulting in a more localised pressure perturbation. These effects are not observed in VMEC, because closed magnetic flux surfaces are enforced. Despite the sensitivity of JOREK results on the assumed resistivity, saturated states can be found using both approaches that are in reasonable agreement, even for this more advanced case.

preprint2021arXiv

An E & B Gyrokinetic Simulation Model for Kinetic Alfvén Waves inTokamak Plasmas

The gyrokinetic particle simulation serves as a powerful tool for the studies of transport, nonlinear phenomenon, and energetic particle physics in tokamak plasmas. While most gyrokinetic simulations make use of the scalar and vector potentials, a new model (GK-E&B) has been developed by using the E and B field in a general and comprehensive form and has been implemented in simulating kinetic Alfvén waves in uniform plasma [Chen et al, Science China Phys. Mechanics & Astronomy 64 (2021)]. In our work, the Chen et al. GK-E&B model has been expressed in general tokamak geometry explicitly using specific coordinates. Its reduction to the uniform plasma is verified and the numerical results show good agreement with the work by Chen et al. The theoretical dispersion relation and numerical results in the local screw pinch model are in excellent agreement. Numerical results show excellent performance in a realistic parameter regime of burning plasmas in terms of high values of $β/( M_e k_\perp^2 ρ_i^2)$, which is challenging for traditional methods due to the ``cancellation' problem. As one application, the GK-E&B model is implemented with kinetic electrons in the local limit. With the matched ITPA-TAE parameters, numerical results show the capability of the GK-E&B in treating the parallel electron Landau damping for realistic tokamak plasma parameters. As another application, the global GK-E&B model is implemented with the dominant electron contribution to $E_\|$ in the cold electron limit. Its capability in simulating the finite $E_\|$ due to the finite electron mass is demonstrated.

preprint2018arXiv

Simulation of the electromagnetic wall response during Vertical Displacement Events (VDE) in ITER tokamak

The key basis for tokamak plasma disruption modeling is to understand how currents flow to the plasma facing surfaces during plasma disruption events. In ITER tokamak, the occurrence of a limited number of major disruptions will definitively damage the chamber with no possibility to restore the device. In the current exchange plasma-wall-plasma, according to the Helmholtz decomposition theorem, our surface current density in the conducting shell - the unknown of our problem - being a vector field twice continuously differentiable in 3D, has been splited into two components: an irrotational (curl-free) vector field and a solenoidal (divergence-free) vector field. Developing a weak formulation form and minimizing the correspondent energy functionals in a Finite Element approach, we have obtained the space and time distribution of the surface currents. We verified successfully our numerical simulation with an analytical solution with pure homogeneous Neumann B.C. and satisfying the necessary existence condition. By considering the iron core presence in JET tokamak, we have split the magnetization currents - the unknowns in some integral equations - into two components, the first producing a magnetic field in the iron region only and the second producing a magnetic field in the vacuum, obtaining thus a better evaluation of the influence of the iron core on the plasma equilibrium. To reduce the influence of the singularities appearing during the surface currents determination in multiply connected domains (L-shaped domains) we have used a conformal transformation method.

preprint2015arXiv

Energy conservation and numerical stability for the reduced MHD models of the non-linear JOREK code

In this paper we present a rigorous derivation of the reduced MHD models with and without parallel velocity that are implemented in the non-linear MHD code JOREK. The model we obtain contains some terms that have been neglected in the implementation but might be relevant in the non-linear phase. These are necessary to guarantee exact conservation with respect to the full MHD energy. For the second part of this work, we have replaced the linearized time stepping of JOREK by a non-linear solver based on the Inexact Newton method including adaptive time stepping. We demonstrate that this approach is more robust especially with respect to numerical errors in the saturation phase of an instability and allows to use larger time steps in the non-linear phase.

preprint2013arXiv

Nonlinear excitation of low-n harmonics in reduced magnetohydrodynamic simulations of edge-localized modes

Nonlinear simulations of the early ELM phase based on a typical type-I ELMy ASDEX Upgrade discharge have been carried out using the reduced MHD code JOREK. The analysis is focused on the evolution of the toroidal Fourier spectrum. It is found that during the nonlinear evolution, linearly subdominant low-n Fourier components, in particular the n = 1, grow to energies comparable with linearly dominant harmonics. A simple model is developed, based on the idea that energy is transferred among the toroidal harmonics via second order nonlinear interaction. The simple model reproduces and explains very well the early nonlinear evolution of the toroidal spectrum in the JOREK simulations. Furthermore, it is shown for the n = 1 harmonic, that its spatial structure changes significantly during the transition from linear to nonlinearly driven growth. The rigidly growing structure of the linearly barely unstable n = 1 reaches far into the plasma core. In contrast, the nonlinearly driven n = 1 has a rigidly growing structure localized at the plasma edge, where the dominant toroidal harmonics driving the n = 1 are maximal and in phase. The presented quadratic coupling model might explain the recent experimental observation of strong low-n components in magnetic measurements [Wenninger et al., Non-linear magnetic perturbations during edge localized modes in TCV dominated by low n mode components, submitted to Nuclear Fusion].

preprint2012arXiv

Non regression testing for the JOREK code

Non Regression Testing (NRT) aims to check if software modifications result in undesired behaviour. Suppose the behaviour of the application previously known, this kind of test makes it possible to identify an eventual regression, a bug. Improving and tuning a parallel code can be a time-consuming and difficult task, especially whenever people from different scientific fields interact closely. The JOREK code aims at investing Magnetohydrodynamic (MHD) instabilities in a Tokamak plasma. This paper describes the NRT procedure that has been tuned for this simulation code. Automation of the NRT is one keypoint to keeping the code healthy in a source code repository.

preprint2012arXiv

Reduced-MHD Simulations of Toroidally and Poloidally Localized ELMs

We use the non-linear reduced-MHD code JOREK to study ELMs in the geometry of the ASDEX Upgrade tokamak. Toroidal mode numbers, poloidal filament sizes, and radial propagation speeds of filaments into the scrape-off layer are in good agreement with observations for type-I ELMs in ASDEX Upgrade. The observed instabilities exhibit a localization of perturbations which is compatible with the "solitary magnetic perturbations" recently discovered in ASDEX Upgrade [R.Wenninger et.al., Solitary Magnetic Perturbations at the ELM Onset, Nucl.Fusion, submitted]. This localization can only be described in numerical simulations with high toroidal resolution.