BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Matt Landreman

Matt Landreman contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.plasm-phmath.OCphysics.app-phphysics.comp-phmath.NANumerical Analysis

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
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

11 published item(s)

preprint2026arXiv

Efficient calculation of magnetic fields from ferromagnetic materials near strong electromagnets, and application to stellarator coil optimization

In fusion reactor design, steels under consideration for the blanket are ferromagnetic, so the steel's effect on the plasma physics must be examined. For efficient calculation of these fields, we can exploit the fact that the magnetic material gives a small perturbation relative to the fields from the electromagnetic coils and plasma. Moreover the magnetization is saturated due to the strong fields in typical fusion systems. These approximations significantly reduce the nonlinearity of the problem, so the magnetic materials can be described by an array of point dipoles of known magnitude, oriented in the direction of the coil and plasma field. The approach is verified by comparison to finite-element calculations with commercial software and shown to be accurate. As no linear or nonlinear solve is required, only evaluation of Biot-Savart-type integrals, the method here is significantly simpler to implement than other methods, and extremely fast. The method is compatible with arbitrary CAD geometry, and also allows rapid computation of the magnetic forces. We demonstrate adding the ferromagnetic effects to free-boundary MHD equilibrium calculations, assessing the effect on plasma properties such as confinement and stability. Moreover, it is straightforward to differentiate through the model to get the derivative of the field with respect to the electromagnet parameters. We thereby demonstrate gradient-based coil optimization for a quasi-isodynamic stellarator in which the field contribution from a ferromagnetic blanket is included. Even a significant steel volume is found to have little impact on the plasma physics properties, with the main effects being a slight destabilization of ballooning modes and a radial shift of the edge islands due to decrease in rotational transform. Both issues are corrected by minor reoptimization of the coil shapes to account for the field from the steel.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Direct computation of magnetic surfaces in Boozer coordinates and coil optimization for quasi-symmetry

We propose a new method to compute magnetic surfaces that are parametrized in Boozer coordinates for vacuum magnetic fields. We also propose a measure for quasi-symmetry on the computed surfaces and use it to design coils that generate a magnetic field that is quasi-symmetric on those surfaces. The rotational transform of the field and complexity measures for the coils are also controlled in the design problem. Using an adjoint approach, we are able to obtain analytic derivatives for this optimization problem, yielding an efficient gradient-based algorithm. Starting from an initial coil set that presents nested magnetic surfaces for a large fraction of the volume, our method converges rapidly to coil systems generating fields with excellent quasi-symmetry and low particle losses. In particular for low complexity coils, we are able to significantly improve the performance compared to coils obtained from the standard two-stage approach, e.g.~reduce losses of fusion-produced alpha particles born at half-radius from $17.7\%$ to $6.6\%$. We also demonstrate 16-coil configurations with alpha loss < $1\%$ and neoclassical transport magnitude $ε_{\mathrm{eff}}^{3/2}$ less than approximately $5\times 10^{-9}.$

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mapping the space of quasisymmetric stellarators using optimized near-axis expansion

A method is demonstrated to rapidly calculate the shapes and properties of quasi-axisymmetric and quasi-helically symmetric stellarators. In this approach, optimization is applied to the equations of magnetohydrodynamic equilibrium and quasisymmetry, expanded in the small distance from the magnetic axis, as formulated by Garren and Boozer [Phys. Fluids B, 3, 2805 (1991)]. Due to the reduction of the equations by the expansion, the computational cost is significantly reduced, to times on the order of 1 cpu-second, enabling wide and high-resolution scans over parameter space. In contrast to traditional stellarator optimization, here the cost function serves to maximize the volume in which the expansion is accurate. A key term in the cost function is $|| \nabla\mathbf B ||$, to maximize scale lengths in the field. Using this method, a database of $5\times 10^5$ optimized configurations is calculated and presented. Quasisymmetric configurations are observed to exist in continuous bands, varying in the ratio of the magnetic axis length to average major radius. Several qualitatively new types of configuration are found, including quasi-helically symmetric fields in which the number of field periods is two or more than six.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimization of quasisymmetric stellarators with self-consistent bootstrap current and energetic particle confinement

Quasisymmetry can greatly improve the confinement of energetic particles and thermal plasma in a stellarator. The magnetic field of a quasisymmetric stellarator at high plasma pressure is significantly affected by the bootstrap current, but the computational cost of accurate stellarator bootstrap calculations has precluded use inside optimization. Here, a new efficient method is demonstrated for optimization of quasisymmetric stellarator configurations such that the bootstrap current profile is consistent with the geometry. The approach is based on the fact that all neoclassical phenomena in quasisymmetry are isomorphic to those in axisymmetry. Therefore accurate formulae for the bootstrap current in tokamaks, which can be evaluated rapidly, can be applied also in stellarators. The deviation between this predicted parallel current and the actual parallel current in the magnetohydrodynamic equilibrium is penalized in the objective function, and the current profile of the equilibrium is included in the parameter space. Quasisymmetric configurations with significant pressure are thereby obtained with self-consistent bootstrap current and excellent confinement. In a comparison of fusion-produced alpha particle confinement across many stellarators, the new configurations have significantly lower alpha energy losses than many previous designs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Single-stage gradient-based stellarator coil design: Optimization for near-axis quasi-symmetry

We present a new coil design paradigm for magnetic confinement in stellarators. Our approach directly optimizes coil shapes and coil currents to produce a vacuum quasi-symmetric magnetic field with a target rotational transform on the magnetic axis. This approach differs from the traditional two-stage approach in which first a magnetic configuration with desirable physics properties is found, and then coils to approximately realize this magnetic configuration are designed. The proposed single-stage approach allows us to find a compromise between confinement and engineering requirements, i.e., find easy-to-build coils with good confinement properties. Using forward and adjoint sensitivities, we derive derivatives of the physical quantities in the objective, which is constrained by a nonlinear periodic differential equation. In two numerical examples, we compare different gradient-based descent algorithms and find that incorporating approximate second-order derivative information through a quasi-Newton method is crucial for convergence. We also explore the optimization landscape in the neighborhood of a minimizer and find many directions in which the objective is mostly flat, indicating ample freedom to find simple and thus easy-to-build coils.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Figures of merit for stellarators near the magnetic axis

A new paradigm for rapid stellarator configuration design has been recently demonstrated, in which the shapes of quasisymmetric or omnigenous flux surfaces are computed directly using an expansion in small distance from the magnetic axis. To further develop this approach, here we derive several other quantities of interest that can be rapidly computed from this near-axis expansion. First, the $\nabla\vec{B}$ and $\nabla\nabla\vec{B}$ tensors are computed, which can be used for direct derivative-based optimization of electromagnetic coil shapes to achieve the desired magnetic configuration. Moreover, if the norm of these tensors is large compared to the field strength for a given magnetic field, the field must have a short length scale, suggesting it may be hard to produce with coils that are suitably far away. Second, we evaluate the minor radius at which the flux surface shapes would become singular, providing a lower bound on the achievable aspect ratio. This bound is also shown to be related to an equilibrium beta limit. Finally, for configurations that are constructed to achieve a desired magnetic field strength to first order in the expansion, we compute the error field that arises due to second order terms.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Magnetic fields with precise quasisymmetry for plasma confinement

Quasisymmetry is an unusual symmetry that can be present in toroidal magnetic fields, enabling confinement of charged particles and plasma. Here it is shown that both quasi-axisymmetry and quasi-helical symmetry can be achieved to a much higher precision than previously thought over a significant volume, resulting in exceptional confinement. For a 1 Tesla mean field far from axisymmetry (vacuum rotational transform $>$ 0.4), symmetry-breaking mode amplitudes throughout a volume of aspect ratio 6 can be made as small as the typical $\sim 50$ $μ$T geomagnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Single-stage gradient-based stellarator coil design: stochastic optimization

We extend the single-stage stellarator coil design approach for quasi-symmetry on axis from [Giuliani et al, 2020] to additionally take into account coil manufacturing errors. By modeling coil errors independently from the coil discretization, we have the flexibility to consider realistic forms of coil errors. The corresponding stochastic optimization problems are formulated using a risk-neutral approach and risk-averse approaches. We present an efficient, gradient-based descent algorithm which relies on analytical derivatives to solve these problems. In a comprehensive numerical study, we compare the coil designs resulting from deterministic and risk-neutral stochastic optimization and find that the risk-neutral formulation results in more robust configurations and reduces the number of local minima of the optimization problem. We also compare deterministic and risk-neutral approaches in terms of quasi-symmetry on and away from the magnetic axis, and in terms of the confinement of particles released close to the axis. Finally, we show that for the optimization problems we consider, a risk-averse objective using the Conditional Value-at-Risk leads to results which are similar to the risk-neutral objective.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Calculation of permanent magnet arrangements for stellarators: A linear least-squares method

A problem arising in several engineering areas is to design magnets outside a volume that produce a desired magnetic field inside it. One instance of this problem is stellarator design, where it has recently been shown that permanent magnets can provide the required shaping of the magnetic field. Here we demonstrate a robust and efficient algorithm REGCOIL_PM to calculate the spatial distribution of these permanent magnets. The procedure involves a small number of fixed-point iterations, with a linear least-squares problem solved at each step. The method exploits the Biot-Savart Law&#39;s exact linearity in magnetization density and approximate linearity in magnet size, for magnets far from the target region. No constraint is placed on the direction of magnetization, so Halbach solutions are found naturally, and the magnitude of the magnetization can be made uniformly equal to a target value.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Impurity temperature screening in stellarators close to quasisymmetry

Impurity temperature screening is a favorable neoclassical phenomenon involving an outward radial flux of impurity ions from the core of fusion devices. Quasisymmetric magnetic fields are intrinsically ambipolar and give rise to temperature screening for low enough $η^{-1}\equiv d\ln n/d\ln T$. In contrast, neoclassical fluxes in generic stellarators will depend on the radial electric field, which is predicted to be inward for ion-root plasmas, potentially leading to impurity accumulation. Here we examine the impurity particle flux in a number of approximately quasisymmetric stellarator configurations and parameter regimes while varying the amount of symmetry-breaking in the magnetic field. The amount of symmetry-breaking (from perfect quasisymmetry) required to lose temperature screening is dependent on the particular configuration and flux surface. Neoclassical fluxes have been obtained using the SFINCS drift-kinetic equation solver with the electrostatic potential $Φ=Φ(ψ)$, where $ψ$ is a flux surface label. Results indicate that achieving temperature screening is possible, but unlikely, at reactor-relevant conditions in the core. Thus, the small departures from symmetry in nominally quasisymmetric stellarators are large enough to significantly alter the neoclassical impurity transport. A further look at these fluxes when compared to a gyro-Bohm turbulence estimate suggests that neoclassical fluxes are small in optimized configurations compared to turbulent fluxes. Therefore, although neoclassical impurity accumulation is expected in most situations, the strength of the turbulence may render this irrelevant.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Adjoint approach to calculating shape gradients for three-dimensional magnetic confinement equilibria. Part II: Applications

The shape gradient is a local sensitivity function that provides the change in a figure of merit associated with a perturbation to the shape of the object. The shape gradient can be used for gradient-based optimization, sensitivity analysis, and tolerance calculations. However, it is generally expensive to compute from finite-difference derivatives for shapes which are described by many parameters, as is the case for stellarator geometry. In an accompanying work (Antonsen et al. 2019), generalized self-adjointness relations are obtained for MHD equilibria. These describe the relation between perturbed equilibria due to changes in the rotational transform or toroidal current profiles, displacements of the plasma boundary, modifications of currents in the vacuum region, or the addition of bulk forces. These are applied to efficiently compute the shape gradient of functions of magnetohydrodynamic (MHD) equilibria with an adjoint approach. In this way, the shape derivative with respect to any perturbation applied to the plasma boundary or coil shapes can be computed with only one additional MHD equilibrium solution. We demonstrate that this approach is applicable for several figures of merit of interest for stellarator configuration optimization: the magnetic well, the magnetic ripple on axis, the departure from quasisymmetry, the effective ripple in the low-collisionality $1/ν$ regime ($ε_{\text{eff}}^{3/2})$ (Nemov et al. 1999), and several finite collisionality neoclassical quantities. Numerical verification is demonstrated for the magnetic well figure of merit with the VMEC code (Hirshman & Whitson 1983) and for the magnetic ripple with modification of the ANIMEC code (Cooper et al. 1992). Comparisons with the direct approach demonstrate that to obtain agreement within several percent, the adjoint approach provides a factor of $\mathcal{O}(10^3)$ in computational savings.

0 citations0 reviewsNo code linkedOpen access