-
Runaway electron beam formation, vertical motion, termination and wall loads in EU-DEMO
Authors:
F. Vannini,
V. Bandaru,
H. Bergstroem,
N. Schwarz,
F. J. Artola,
M. Hoelzl,
G. Pautasso,
E. Nardon,
F. Maviglia,
M. L. Richiusa,
E. Emanuelli,
the JOREK team
Abstract:
Runaway electron loads onto material structures are a major concern for future large tokamaks due to the efficient avalanching at high plasma currents. Here, we perform predictive studies using the JOREK code for a plausible plasma configuration in the European DEMO fusion power plant with focus on a pessimistic scenario in which a multi mega-ampere runaway electron beam is formed. The work first…
▽ More
Runaway electron loads onto material structures are a major concern for future large tokamaks due to the efficient avalanching at high plasma currents. Here, we perform predictive studies using the JOREK code for a plausible plasma configuration in the European DEMO fusion power plant with focus on a pessimistic scenario in which a multi mega-ampere runaway electron beam is formed. The work first comprises axisymmetric predictions of runaway electron beam formation in a mitigated scenario and of the simultaneous vertical motion of the beam due to loss of position control. The subsequent runaway electron beam termination triggered by a burst of MHD activity during the course of the vertical motion is then simulated in 3D with the runaway electron fluid self-consistently coupled to the MHD modes. Finally, the resulting deposition pattern of the runaway electrons onto wall structures is calculated with a relativistic test particle approach. This way, the suitability of a possible sacrificial limiter concept for the protection of first wall components is assessed.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Axisymmetric predictions for mitigated and vertically unstable disruptions in ITER with runaway electrons
Authors:
V. Bandaru,
M. Hoelzl,
F. J. Artola,
M. Lehnen,
JOREK team
Abstract:
We present two-dimensional global simulations of mitigated and vertically unstable disruptions in ITER in the presence of runaway electrons. An elongated plasma in free-boundary equilibrium is subjected to an artificial thermal quench and current-profile flattening, followed by one or more massive material injections and RE avalanche. Scenarios of major disruptions as well as upward and downward v…
▽ More
We present two-dimensional global simulations of mitigated and vertically unstable disruptions in ITER in the presence of runaway electrons. An elongated plasma in free-boundary equilibrium is subjected to an artificial thermal quench and current-profile flattening, followed by one or more massive material injections and RE avalanche. Scenarios of major disruptions as well as upward and downward vertical displacement events are considered. Results provide important insights into the effects of runaway electron formation, post thermal quench current-profile, injection quantities and timings, and impurity flushout on the overall evolution of disruption and the plasma vertical motion thereof. Interplay between the various effects offers scope for potentially beneficial runaway electron mitigation scenarios.
△ Less
Submitted 14 June, 2024;
originally announced June 2024.
-
The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas
Authors:
M Hoelzl,
GTA Huijsmans,
SJP Pamela,
M Becoulet,
E Nardon,
FJ Artola,
B Nkonga,
CV Atanasiu,
V Bandaru,
A Bhole,
D Bonfiglio,
A Cathey,
O Czarny,
A Dvornova,
T Feher,
A Fil,
E Franck,
S Futatani,
M Gruca,
H Guillard,
JW Haverkort,
I Holod,
D Hu,
SK Kim,
SQ Korving
, et al. (28 additional authors not shown)
Abstract:
JOREK is a massively parallel fully implicit non-linear extended MHD code for realistic tokamak X-point plasmas. It has become a widely used versatile code for studying large-scale plasma instabilities and their control developed in an international community. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and phys…
▽ More
JOREK is a massively parallel fully implicit non-linear extended MHD code for realistic tokamak X-point plasmas. It has become a widely used versatile code for studying large-scale plasma instabilities and their control developed in an international community. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench and current quench triggered by massive gas injection (MGI) and shattered pellet injection (SPI), runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events (VDEs). Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring thermal quenches triggered by locked modes, and radiative collapses are being studied.
△ Less
Submitted 21 April, 2021; v1 submitted 18 November, 2020;
originally announced November 2020.
-
Simulating the nonlinear interaction of relativistic electrons and tokamak plasma instabilities: Implementation and validation of a fluid model
Authors:
V. Bandaru,
M. Hoelzl,
F. J. Artola,
G. Papp,
G. T. A. Huijsmans
Abstract:
For the simulation of disruptions in tokamak fusion plasmas, a fluid model describing the evolution of relativistic runaway electrons and their interaction with the background plasma is presented. The overall aim of the model is to self-consistently describe the nonlinear coupled evolution of runaway electrons (REs) and plasma instabilities during disruptions. In this model, the runaway electrons…
▽ More
For the simulation of disruptions in tokamak fusion plasmas, a fluid model describing the evolution of relativistic runaway electrons and their interaction with the background plasma is presented. The overall aim of the model is to self-consistently describe the nonlinear coupled evolution of runaway electrons (REs) and plasma instabilities during disruptions. In this model, the runaway electrons are considered as a separate fluid species in which the initial seed is generated through the Dreicer source, which eventually grows by the avalanche mechanism (further relevant source mechanisms can easily be added). Advection of the runaway electrons is considered primarily along field lines, but also taking into account the ExB drift. The model is implemented in the nonlinear magnetohydrodynamic (MHD) code JOREK based on Bezier finite elements, with current coupling to the thermal plasma. Benchmarking of the code with the one-dimensional runaway electron code GO is done using an artificial thermal quench on a circular plasma. As a first demonstration, the code is applied to the problem of an axisymmetric cold vertical displacement event in an ITER plasma, revealing significantly different dynamics between cases computed with and without runaway electrons. Though it is not yet feasible to achieve fully realistic runaway electron velocities close to the speed of light in complete simulations of slowly evolving plasma instabilities, the code is demonstrated to be suitable to study various kinds of MHD-RE interactions in MHD-active and disruption relevant plasmas.
△ Less
Submitted 28 June, 2019;
originally announced June 2019.
-
Transitional boundary layers in low-Prandtl-number convection
Authors:
Jörg Schumacher,
Vinodh Bandaru,
Ambrish Pandey,
Janet D. Scheel
Abstract:
The boundary layer structure of the velocity and temperature fields in turbulent Rayleigh-Benard flows in closed cylindrical cells of unit aspect ratio is revisited from a transitional and turbulent viscous boundary layer perspective. When the Rayleigh number is large enough, the dynamics at the bottom and top plates can be separated into an impact region of downwelling plumes, an ejection region…
▽ More
The boundary layer structure of the velocity and temperature fields in turbulent Rayleigh-Benard flows in closed cylindrical cells of unit aspect ratio is revisited from a transitional and turbulent viscous boundary layer perspective. When the Rayleigh number is large enough, the dynamics at the bottom and top plates can be separated into an impact region of downwelling plumes, an ejection region of upwelling plumes and an interior region away from the side walls. The latter is dominated by the shear of the large-scale circulation (LSC) roll which fills the whole cell and continuously varies its orientation. The working fluid is liquid mercury or gallium at a Prandtl number Pr=0.021 for Rayleigh numbers between Ra=3e+5 and 4e+8. The generated turbulent momentum transfer corresponds to macroscopic flow Reynolds numbers with values between 1800 and 46000. It is shown that the viscous boundary layers for the largest Rayleigh numbers are highly transitional and obey properties that are directly comparable to transitional channel flows at friction Reynolds numbers Re_tau slightly below 100. The transitional character of the viscous boundary layer is also underlined by the strong enhancement of the fluctuations of the wall stress components with increasing Rayleigh number. An extrapolation of our analysis data suggests that the friction Reynolds number Re_tau in the velocity boundary layer can reach values of 200 for Ra beyond 1e+11. Thus the viscous boundary layer in a liquid metal flow would become turbulent at a much lower Rayleigh number than for turbulent convection in gases and gas mixtures.
△ Less
Submitted 9 December, 2016;
originally announced December 2016.
-
A hybrid finite-difference/boundary element procedure for the simulation of turbulent MHD duct flow at finite magnetic Reynolds numbers
Authors:
Vinodh Bandaru,
Thomas Boeck,
Dmitry Krasnov,
Jörg Schumacher
Abstract:
A conservative coupled finite difference-boundary element computational procedure for the simulation of turbulent magnetohydrodynamic flow in a straight rectangular duct at finite magnetic Reynolds number is presented. The flow is assumed to be periodic in the streamwise direction and is driven by a mean pressure gradient. The duct walls are considered to be electrically insulating. The co-evoluti…
▽ More
A conservative coupled finite difference-boundary element computational procedure for the simulation of turbulent magnetohydrodynamic flow in a straight rectangular duct at finite magnetic Reynolds number is presented. The flow is assumed to be periodic in the streamwise direction and is driven by a mean pressure gradient. The duct walls are considered to be electrically insulating. The co-evolution of the velocity and magnetic fields as described respectively by the Navier-Stokes and the magnetic induction equations, together with the coupling of the magnetic field between the conducting domain and the non-conducting exterior is solved using the magnetic field formulation. The aim is to simulate localized magnetic fields interacting with turbulent duct flow. Detailed verification of the implementation of the numerical scheme is conducted in the limiting case of low magnetic Reynolds number by comparing with the results obtained using a quasistatic approach that has no coupling with the exterior. The rigorous procedure with non-local magnetic boundary conditions is compared versus simplified pseudo-vacuum boundary conditions and the differences are quantified. Our first direct numerical simulations of turbulent Hartmann duct flow at moderate magnetic Reynolds numbers and a low flow Reynolds number show significant differences in the duct flow turbulence, even at low interaction level between the flow and magnetic field
△ Less
Submitted 5 October, 2015;
originally announced October 2015.
-
Role of critical points of the skin friction field in formation of plumes in thermal convection
Authors:
Vinodh Bandaru,
Anastasiya Kolchinskaya,
Kathrin Padberg-Gehle,
Jörg Schumacher
Abstract:
The dynamics in the thin boundary layers of temperature and velocity is the key to a deeper understanding of turbulent transport of heat and momentum in thermal convection. The velocity gradient at the hot and cold plates of a Rayleigh-Bénard convection cell forms the two-dimensional skin friction field and is related to the formation of thermal plumes in the respective boundary layers. Our analys…
▽ More
The dynamics in the thin boundary layers of temperature and velocity is the key to a deeper understanding of turbulent transport of heat and momentum in thermal convection. The velocity gradient at the hot and cold plates of a Rayleigh-Bénard convection cell forms the two-dimensional skin friction field and is related to the formation of thermal plumes in the respective boundary layers. Our analysis is based on a direct numerical simulation of Rayleigh-Bénard convection in a closed cylindrical cell of aspect ratio $Γ=1$ and focused on the critical points of the skin friction field. We identify triplets of critical points, which are composed of two unstable nodes and a saddle between them, as the characteristic building block of the skin friction field. Isolated triplets as well as networks of triplets are detected. The majority of the ridges of line-like thermal plumes coincide with the unstable manifolds of the saddles. From a dynamical Lagrangian perspective, thermal plumes are formed together with an attractive hyperbolic Lagrangian Coherent Structure of the skin friction field. We also discuss the differences from the skin friction field in turbulent channel flows from the perspective of the Poincaré-Hopf index theorem for two-dimensional vector fields.
△ Less
Submitted 16 September, 2015;
originally announced September 2015.
-
Simulation of flux expulsion and associated dynamics in a two-dimensional magnetohydrodynamic channel flow
Authors:
Vinodh Bandaru,
Julian Pracht,
Thomas Boeck,
Jörg Schumacher
Abstract:
We consider a plane channel flow of an electrically conducting fluid which is driven by a mean pressure gradient in the presence of an applied magnetic field that is streamwise periodic with zero mean. Magnetic flux expulsion and the associated bifurcation in such a configuration is explored using direct numerical simulations (DNS). The structure of the flow and magnetic fields in the Hartmann reg…
▽ More
We consider a plane channel flow of an electrically conducting fluid which is driven by a mean pressure gradient in the presence of an applied magnetic field that is streamwise periodic with zero mean. Magnetic flux expulsion and the associated bifurcation in such a configuration is explored using direct numerical simulations (DNS). The structure of the flow and magnetic fields in the Hartmann regime (where the dominant balance is through Lorentz forces) and the Poiseuille regime (where viscous effects play a significant role) are studied and detailed comparisons to the existing one-dimensional model of Kamkar and Moffatt (J. Fluid. Mech., Vol.90, pp 107-122, 1982) are drawn to evaluate the validity of the model. Comparisons show good agreement of the model with DNS in the Hartmann regime, but significant diferences arising in the Poiseuille regime when non-linear effects become important. The effects of various parameters like the magnetic Reynolds number, imposed field wavenumber etc. on the bifurcation of the flow are studied. Magnetic field line reconnections occuring during the dynamic runaway reveal a specific two-step pattern that leads to the gradual expulsion of flux in the core region.
△ Less
Submitted 9 May, 2015;
originally announced May 2015.