-
Direct prediction of saturated neoclassical tearing modes in slab using an equilibrium approach
Authors:
Erol Balkovic,
Joaquim Loizu,
Jonathan P. Graves,
Yi-Min Huang,
Christopher B. Smiet
Abstract:
We demonstrate for the first time that the nonlinear saturation of neoclassical tearing modes (NTMs) can be found directly using a variational principle based on Taylor relaxation, without needing to simulate the intermediate, resistivity-dependent dynamics. As in previous investigations of classical tearing mode saturation (Loizu et al. 2020; Loizu & Bonfiglio 2023), we make use of SPEC (Hudson e…
▽ More
We demonstrate for the first time that the nonlinear saturation of neoclassical tearing modes (NTMs) can be found directly using a variational principle based on Taylor relaxation, without needing to simulate the intermediate, resistivity-dependent dynamics. As in previous investigations of classical tearing mode saturation (Loizu et al. 2020; Loizu & Bonfiglio 2023), we make use of SPEC (Hudson et al. 2012), an equilibrium solver based on the variational principle of the Multi-Region relaxed MHD, featuring stepped pressure profiles and arbitrary magnetic topology. We work in slab geometry and employ a simple bootstrap current model $J_\textrm{bs} = C \nabla p$ to study the bootstrap-driven tearing modes, scanning over the asymptotic matching parameter $Δ'$ and the bootstrap current strength. Saturated island widths produced by SPEC agree well with the predictions of an initial value resistive MHD code (Huang & Bhattacharjee 2016) while being orders of magnitude faster to calculate. Additionally, we observe good agreement with a simple analytical Modified Rutherford Equation, without requiring any fitting coefficients. The match is obtained for both linearly unstable classical tearing modes in the presence of bootstrap current, and neoclassical tearing modes, which are linearly stable but nonlinear-unstable due to the effects of the bootstrap current
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Efficient single-stage optimization of islands in finite-$β$ stellarator equilibria
Authors:
Christopher Berg Smiet,
Joaquim Loizu,
Erol Balkovic,
Antoine Baillod
Abstract:
We present the first single-stage optimization of islands in finite-$β$ stellarator equilibria. Stellarator optimization is traditionally performed as a two-stage process; in the first stage, an optimal equilibrium is calculated which balances a set of competing constraints, and in the second stage a set of coils is found that supports said equilibrium. Stage one is generally performed using a rep…
▽ More
We present the first single-stage optimization of islands in finite-$β$ stellarator equilibria. Stellarator optimization is traditionally performed as a two-stage process; in the first stage, an optimal equilibrium is calculated which balances a set of competing constraints, and in the second stage a set of coils is found that supports said equilibrium. Stage one is generally performed using a representation for the equilibrium that assumes nestedness of flux surfaces, even though this is not warranted and occasionally undesired. The second stage optimization of coils is never perfect, and the mismatch leads to worse performing equilibria, and further deteriorates if additional constraints such as force minimization, coil torsion or port access are included. The higher fidelity of single-stage optimization is especially important for the optimization of islands as these are incredibly sensitive to changes in the field. In this paper we demonstrate an optimization scheme capable of optimizing islands in finite $β$ stellarator equilibria directly from coils. We furthermore develop and demonstrate a method to reduce the dimensionality of the single-stage optimization problem to that of the first stage in the two-stage approach.
△ Less
Submitted 2 July, 2024;
originally announced July 2024.
-
A topological approach to magnetic nulls
Authors:
Ben Y. Israeli,
Christopher Berg Smiet
Abstract:
Magnetic nulls are locations where the magnetic field vanishes. Nulls are the location of magnetic reconnection, and they determine to a large degree the magnetic connectivity in a system. We describe a novel approach to understanding movement, appearance, and disappearance of nulls in magnetic fields. This approach is based on the novel concept of isotropes, or lines where the field direction is…
▽ More
Magnetic nulls are locations where the magnetic field vanishes. Nulls are the location of magnetic reconnection, and they determine to a large degree the magnetic connectivity in a system. We describe a novel approach to understanding movement, appearance, and disappearance of nulls in magnetic fields. This approach is based on the novel concept of isotropes, or lines where the field direction is constant. These lines are streamlines of a vector field whose flux is sourced by the topological indices of nulls, and can be conceptualized as corresponding "lines of force" between nulls. We show how this topological approach can be used to generate analytical expressions for the location of nulls in the presence of external fields for dipoles and for a field defined by the Hopf fibration.
△ Less
Submitted 8 August, 2023;
originally announced August 2023.
-
EUV Debris Mitigation using Magnetic Nulls
Authors:
Ben Y. Israeli,
Christopher Berg Smiet,
Marien Simeni Simeni,
Ahmed Diallo
Abstract:
Next generation EUV sources for photolithography use light produced by laser-produced plasmas (LPP) from ablated tin droplets. A major challenge for extending the lifetime of these devices is mitigating damage caused by deposition of tin debris on the sensitive collection mirror. Especially difficult to stop are high energy (up to 10 keV) highly charged tin ions created in the plasma. Existing sol…
▽ More
Next generation EUV sources for photolithography use light produced by laser-produced plasmas (LPP) from ablated tin droplets. A major challenge for extending the lifetime of these devices is mitigating damage caused by deposition of tin debris on the sensitive collection mirror. Especially difficult to stop are high energy (up to 10 keV) highly charged tin ions created in the plasma. Existing solutions include the use of stopping gas, electric fields, and magnetic fields. One common configuration consists of a magnetic field perpendicular to the EUV emission direction, but such a system can result in ion populations that are trapped rather than removed. We investigate a previously unconsidered mitigation geometry consisting of a magnetic null by performing full-orbit integration of the ion trajectories in an EUV system with realistic dimensions, and optimize the coil locations for the null configuration. The magnetic null prevents a fraction of ions from hitting the mirror comparable to that of the perpendicular field, but does not trap any ions due to the chaotic nature of ion trajectories that pass close to the null. This technology can potentially improve LPP-based EUV photolithography system efficiency and lifetime, and may allow for a different, more efficient formulation of buffer gas.
△ Less
Submitted 20 July, 2023;
originally announced July 2023.
-
Bifurcations of the magnetic axis and the alternating-hyperbolic sawtooth
Authors:
C. B. Smiet,
G. J. Kramer,
S. R. Hudson
Abstract:
We present a sawtooth model that explains observations where the central safety factor, $q_0$, stays well below one, which is irreconcilable with current models that predict a reset to $q_0=1$ after the crash. We identify the structure of the field around the magnetic axis with elements of the Lie group $\mathrm{SL}(2,\mathbb{R})$ and find a transition to an alternating-hyperbolic geometry when…
▽ More
We present a sawtooth model that explains observations where the central safety factor, $q_0$, stays well below one, which is irreconcilable with current models that predict a reset to $q_0=1$ after the crash. We identify the structure of the field around the magnetic axis with elements of the Lie group $\mathrm{SL}(2,\mathbb{R})$ and find a transition to an alternating-hyperbolic geometry when $q_0=2/3$. This transition is driven by an ideal MHD instability and leads to a chaotic magnetic field near the axis.
△ Less
Submitted 8 July, 2020;
originally announced July 2020.
-
Resistive evolution of toroidal field configurations and their relation to magnetic clouds
Authors:
Christopher B. Smiet,
Hugo J. de Blank,
Tobias A. de Jong,
David N. L. Kok,
Dirk Bouwmeester
Abstract:
We study the resistive evolution of a localized self-organizing magnetohydrodynamic equilibrium. In this configuration the magnetic forces are balanced by a pressure force caused by a toroidal depression in the pressure. Equilibrium is attained when this low pressure region prevents further expansion into the higher-pressure external plasma. We find that, for the parameters investigated, the resis…
▽ More
We study the resistive evolution of a localized self-organizing magnetohydrodynamic equilibrium. In this configuration the magnetic forces are balanced by a pressure force caused by a toroidal depression in the pressure. Equilibrium is attained when this low pressure region prevents further expansion into the higher-pressure external plasma. We find that, for the parameters investigated, the resistive evolution of the structures follows a universal pattern when rescaled to resistive time. The finite resistivity causes both a decrease in the magnetic field strength and a finite slip of the plasma fluid against the static equilibrium. This slip is caused by a Pfirsch-Schlüter type diffusion, similar to what is seen in tokamak equilibria. The net effect is that the configuration remains in Magnetostatic equilibrium whilst it slowly grows in size. The rotational transform of the structure becomes nearly constant throughout the entire structure, and decreases according to a power law. In simulations this equilibrium is observed when highly tangled field lines relax in a high-pressure (relative to the magnetic field strength) environment, a situation that occurs when the twisted field of a coronal loop is ejected into the interplanetary solar wind. In this paper we relate this localized MHD equilibrium to magnetic clouds in the solar wind.
△ Less
Submitted 3 December, 2018;
originally announced December 2018.
-
Ideal Relaxation of the Hopf Fibration
Authors:
Christopher Berg Smiet,
Simon Candelaresi,
Dirk Bouwmeester
Abstract:
Ideal MHD relaxation is the topology-conserving reconfiguration of a magnetic field into a lower energy state where the net force is zero. This is achieved by modeling the plasma as perfectly conducting viscous fluid. It is an important tool for investigating plasma equilibria and is often used to study the magnetic configurations in fusion devices and astrophysical plasmas. We study the equilibri…
▽ More
Ideal MHD relaxation is the topology-conserving reconfiguration of a magnetic field into a lower energy state where the net force is zero. This is achieved by modeling the plasma as perfectly conducting viscous fluid. It is an important tool for investigating plasma equilibria and is often used to study the magnetic configurations in fusion devices and astrophysical plasmas. We study the equilibrium reached by a localized magnetic field through the topology conserving relaxation of a magnetic field based on the Hopf fibration in which magnetic field lines are closed circles that are all linked with one another. Magnetic fields with this topology have recently been shown to occur in non-ideal numerical simulations. Our results show that any localized field can only attain equilibrium if there is a finite external pressure, and that for such a field a Taylor state is unattainable. We find an equilibrium plasma configuration that is characterized by a lowered pressure in a toroidal region, with field lines lying on surfaces of constant pressure. Therefore, the field is in a Grad-Shafranov equilibrium. Localized helical magnetic fields are found when plasma is ejected from astrophysical bodies and subsequently relaxes against the background plasma, as well as on earth in plasmoids generated by e.g.\ a Marshall gun. This work shows under which conditions an equilibrium can be reached and identifies a toroidal depression as the characteristic feature of such a configuration.
△ Less
Submitted 25 October, 2017; v1 submitted 15 October, 2016;
originally announced October 2016.
-
Self-organizing Knotted Magnetic Structures in Plasma
Authors:
C. B. Smiet,
S. Candelaresi,
A. Thompson,
J. Swearngin,
J. W. Dalhuizen,
D. Bouwmeester
Abstract:
We perform full-MHD simulations on various initially helical configurations and show that they reconfigure into a state where the magnetic field lines span nested toroidal surfaces. This relaxed configuration is not a Taylor state, as is often assumed for relaxing plasma, but a state where the Lorentz force is balanced by the hydrostatic pressure, which is lowest on the central ring of the nested…
▽ More
We perform full-MHD simulations on various initially helical configurations and show that they reconfigure into a state where the magnetic field lines span nested toroidal surfaces. This relaxed configuration is not a Taylor state, as is often assumed for relaxing plasma, but a state where the Lorentz force is balanced by the hydrostatic pressure, which is lowest on the central ring of the nested tori. Furthermore, the structure is characterized by a spatially slowly varying rotational transform, which leads to the formation of a few magnetic islands at rational surfaces. We then obtain analytic expressions that approximate the global structure of the quasi-stable linked and knotted plasma configurations that emerge, using maps from $S^3$ to $S^2$ of which the Hopf fibration is a special case. The knotted plasma configurations have a highly localized magnetic energy density and retain their structure on time scales much longer than the Alfvenic time scale.
△ Less
Submitted 2 September, 2015; v1 submitted 31 July, 2015;
originally announced July 2015.