Source author record

M. Landreman

M. Landreman 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
1topics
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)

preprint2021arXiv

Stellarator optimization for nested magnetic surfaces at finite $β$ and toroidal current

Good magnetic surfaces, as opposed to magnetic islands and chaotic field lines, are generally desirable for stellarators. In previous work, M. Landreman et al. [Phys. of Plasmas 28, 092505 (2021)] showed that equilibria computed by the Stepped-Pressure Equilibrium Code (SPEC) [S. P. Hudson et al., Phys. Plasmas 19, 112502 (2012)] could be optimized for good magnetic surfaces in vacuum. In this paper, we build upon their work to show the first finite-$β$, fixed- and free-boundary optimization of SPEC equilibria for good magnetic surfaces. The objective function is constructed with the Greene's residue of selected rational surfaces and the optimization is driven by the SIMSOPT framework [M. Landreman et al., J. Open Source Software 6, 3525 (2021)]. We show that the size of magnetic islands and the consequent regions occupied by chaotic field lines can be minimized in a classical stellarator geometry by optimizing either the injected toroidal current profile, the shape of a perfectly conducting wall surrounding the plasma (fixed-boundary case), or the coils (free-boundary case), in a reasonable amount of computational time. This work shows that SPEC can be used as an equilibrium code both in a two-step or single-step stellarator optimization loop.

preprint2020arXiv

Construction of Quasisymmetric Stellarators Using a Direct Coordinate Approach

Optimized stellarator configurations and their analytical properties are obtained using a near-axis expansion approach. Such configurations are associated with good confinement as the guiding center particle trajectories and neoclassical transport are isomorphic to those in a tokamak. This makes them appealing as fusion reactor candidates. Using a direct coordinate approach, where the magnetic field and flux surface functions are found explicitly in terms of the position vector at successive orders in the distance to the axis, the set of ordinary differential equations for first and second order quasisymmetry is derived. Examples of quasi-axisymmetric shapes are constructed using a pseudospectral numerical method. Finally, the direct coordinate approach is used to independently verify two hypotheses commonly associated with quasisymmetric magnetic fields, namely that the number of equations exceeds the number of parameters at third order in the expansion and that the near-axis expansion does not prohibit exact quasisymmetry from being achieved on a single flux surface.

preprint2020arXiv

Near-Axis Expansion of Stellarator Equilibrium at Arbitrary Order in the Distance to the Axis

A direct construction of equilibrium magnetic fields with toroidal topology at arbitrary order in the distance from the magnetic axis is carried out, yielding an analytical framework able to explore the landscape of possible magnetic flux surfaces in the vicinity of the axis. This framework can provide meaningful analytical insight on the character of high-aspect-ratio stellarator shapes, such as the dependence of the rotational transform and the plasma beta-limit on geometrical properties of the resulting flux surfaces. The approach developed here is based on an asymptotic expansion on the inverse aspect-ratio of the ideal MHD equation. The analysis is simplified by using an orthogonal coordinate system relative to the Frenet-Serret frame at the magnetic axis. The magnetic field vector, the toroidal magnetic flux, the current density, the field line label and the rotational transform are derived at arbitrary order in the expansion parameter. Moreover, a comparison with a near-axis expansion formalism employing an inverse coordinate method based on Boozer coordinates (the so-called Garren-Boozer construction) is made, where both methods are shown to agree at lowest order. Finally, as a practical example, a numerical solution using a W7-X equilibrium is presented, and a comparison between the lowest order solution and the W7-X magnetic field is performed.

preprint2019arXiv

Direct construction of optimized stellarator shapes. III. Omnigenity near the magnetic axis

The condition of omnigenity is investigated, and applied to the near-axis expansion of Garren and Boozer (1991a). Due in part to the particular analyticity requirements of the near-axis expansion, we find that, excluding quasi-symmetric solutions, only one type of omnigenity, namely quasi-isodynamicity, can be satisfied at first order in the distance from the magnetic axis. Our construction provides a parameterization of the space of such solutions, and the cylindrical reformulation and numerical method of Landreman and Sengupta (2018); Landreman et al. (2019), enables their efficient numerical construction.

preprint2016arXiv

Global effects on neoclassical transport in the pedestal with impurities

We present a numerical study of collisional transport in a tokamak pedestal in the presence of non-trace impurities, using the radially global $δf$ neoclassical solver PERFECT [M. Landreman et al. 2014 Plasma Phys. Control. Fusion 56 045005]. It is known that in a tokamak core with non-trace impurities present the radial impurity flux opposes the bulk ion flux to provide an ambipolar particle transport, with the electron transport being negligibly small. However, in a sharp density pedestal with sub-sonic ion flows the electron transport can be comparable to the ion and impurity flows. Furthermore, the neoclassical particle transport is not intrinsically ambipolar, and the non-ambipolarity of the fluxes extends outside the pedestal region by the radial coupling of the perturbations. The neoclassical momentum transport, which is finite in the presence of ion orbit-width scale profile variations, is significantly enhanced when impurities are present in non-trace quantities, even if the total parallel mass flow is dominated by the bulk ions.

preprint2016arXiv

Kinetic modelling of runaway electrons in dynamic scenarios

Improved understanding of runaway-electron formation and decay processes are of prime interest for the safe operation of large tokamaks, and the dynamics of the runaway electrons during dynamical scenarios such as disruptions are of particular concern. In this paper, we present kinetic modelling of scenarios with time-dependent plasma parameters; in particular, we investigate hot-tail runaway generation during a rapid drop in plasma temperature. With the goal of studying runaway-electron generation with a self-consistent electric-field evolution, we also discuss the implementation of a collision operator that conserves momentum and energy and demonstrate its properties. An operator for avalanche runaway-electron generation, which takes the energy dependence of the scattering cross section and the runaway distribution into account, is investigated. We show that the simplified avalanche model of Rosenbluth & Putvinskii [Nucl. Fusion 1997 37 1355] can give inaccurate results for the avalanche growth rate (either lower or higher) for many parameters, especially when the average runaway energy is modest, such as during the initial phase of the avalanche multiplication. The developments presented pave the way for improved modelling of runaway-electron dynamics during disruptions or other dynamic events.

preprint2016arXiv

NORSE: A solver for the relativistic non-linear Fokker-Planck equation for electrons in a homogeneous plasma

Energetic electrons are of interest in many types of plasmas, however previous modeling of their properties has been restricted to the use of linear Fokker-Planck collision operators or non-relativistic formulations. Here, we describe a fully non-linear kinetic-equation solver, capable of handling large electric-field strengths (compared to the Dreicer field) and relativistic temperatures. This tool allows modeling of the momentum-space dynamics of the electrons in cases where strong departures from Maxwellian distributions may arise. As an example, we consider electron runaway in magnetic-confinement fusion plasmas and describe a transition to electron slide-away at field strengths significantly lower than previously predicted.

preprint2016arXiv

Runaway-electron formation and electron slide-away in an ITER post-disruption scenario

Mitigation of runaway electrons is one of the outstanding issues for the reliable operation of ITER and other large tokamaks, and accurate estimates for the expected runaway-electron energies and current are needed. Previously, linearized tools (which assume the runaway population to be small) have been used to study the runaway dynamics, but these tools are not valid in the cases of most interest, i.e. when the runaway population becomes substantial. We study runaway-electron formation in a post-disruption ITER plasma using the newly developed non-linear code NORSE, and describe a feedback mechanism by which a transition to electron slide-away can be induced at field strengths significantly lower than previously expected. If the electric field is actively imposed using the control system, the entire electron population is quickly converted to runaways in the scenario considered. We find the time until the feedback mechanism sets in to be highly dependent on the details of the mechanisms removing heat from the thermal electron population.

preprint2013arXiv

Synchrotron radiation from a runaway electron distribution in tokamaks

The synchrotron radiation emitted by runaway electrons in a fusion plasma provides information regarding the particle momenta and pitch-angles of the runaway electron population through the strong dependence of the synchrotron spectrum on these parameters. Information about the runaway density and its spatial distribution, as well as the time evolution of the above quantities, can also be deduced. In this paper we present the synchrotron radiation spectra for typical avalanching runaway electron distributions. Spectra obtained for a distribution of electrons are compared to the emission of mono-energetic electrons with a prescribed pitch-angle. We also examine the effects of magnetic field curvature and analyse the sensitivity of the resulting spectrum to perturbations to the runaway distribution. The implications for the deduced runaway electron parameters are discussed. We compare our calculations to experimental data from DIII-D and estimate the maximum observed runaway energy.

preprint2011arXiv

Impurity flows and plateau-regime poloidal density variation in a tokamak pedestal

In the pedestal of a tokamak, the sharp radial gradients of density and temperature can give rise to poloidal variation in the density of impurities. At the same time, the flow of the impurity species is modified relative to the conventional neoclassical result. In this paper, these changes to the density and flow of a collisional impurity species are calculated for the case when the main ions are in the plateau regime. In this regime it is found that the impurity density can be higher at either the inboard or outboard side. This finding differs from earlier results for banana- or Pfirsch-Schlüter-regime main ions, in which case the impurity density is always higher at the inboard side in the absence of rotation. Finally, the modifications to the impurity flow are also given for the other regimes of main-ion collisionality.