BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Matt Landreman

Matt Landreman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.plasm-phphysics.comp-phmath.NAmath.OCphysics.app-phNumerical Analysis

Catalog footprint

What is connected

28works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

28 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'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
preprint2016arXiv

Efficient magnetic fields for supporting toroidal plasmas

The magnetic field that supports tokamak and stellarator plasmas must be produced by coils well separated from the plasma. However the larger the separation, the more difficult it is to produce a given magnetic field in the plasma region, so plasma configurations should be chosen that can be supported as efficiently as possible by distant coils. The efficiency of an externally-generated magnetic field is a measure of the field's shaping component magnitude at the plasma compared to the magnitude near the coils; the efficiency of a plasma equilibrium can be measured using the efficiency of the required external shaping field. Counterintuitively, plasma shapes with low curvature and spectral width may have low efficiency, whereas plasma shapes with sharp edges may have high efficiency. Two precise measures of magnetic field efficiency, which correctly identify such differences in difficulty, will be examined. These measures, which can be expressed as matrices, relate the externally-produced normal magnetic field on the plasma surface to the either the normal field or current on a distant control surface. A singular value decomposition (SVD) of either matrix yields an efficiency ordered basis for the magnetic field distributions. Calculations are carried out for both tokamak and stellarator cases. For axisymmetric surfaces with circular cross-section, the SVD is calculated analytically, and the range of poloidal and toroidal mode numbers that can be controlled to a given desired level is determined. If formulated properly, these efficiency measures are independent of the coordinates used to parameterize the surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Transport and deceleration of fusion products in microturbulence

The velocity-space distribution of alpha particles born in fusion devices is subject to modification at moderate energies due to turbulent transport. Therefore, one must calculate the evolution of an equilibrium distribution whose functional form is not known \emph{a priori}. Using a novel technique, applicable to any trace impurity, we have made this calculation not only possible, but particularly efficient. We demonstrate a microturbulence-induced departure from the local slowing-down distribution, an inversion of the energy distribution, and associated modifications to the alpha heating and pressure profiles in an ITER-like scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accurate Spectral Numerical Schemes for Kinetic Equations with Energy Diffusion

We examine the merits of using a family of polynomials that are orthogonal with respect to a non-classical weight function to discretize the speed variable in continuum kinetic calculations. We consider a model one-dimensional partial differential equation describing energy diffusion in velocity space due to Fokker-Planck collisions. This relatively simple case allows us to compare the results of the projected dynamics with an expensive but highly accurate spectral transform approach. It also allows us to integrate in time exactly, and to focus entirely on the effectiveness of the discretization of the speed variable. We show that for a fixed number of modes or grid points, the non-classical polynomials can be many orders of magnitude more accurate than classical Hermite polynomials or finite-difference solvers for kinetic equations in plasma physics. We provide a detailed analysis of the difference in behavior and accuracy of the two families of polynomials. For the non-classical polynomials, if the initial condition is not smooth at the origin when interpreted as a three-dimensional radial function, the exact solution leaves the polynomial subspace for a time, but returns (up to roundoff accuracy) to the same point evolved to by the projected dynamics in that time. By contrast, using classical polynomials, the exact solution differs significantly from the projected dynamics solution when it returns to the subspace. We also explore the connection between eigenfunctions of the projected evolution operator and (non-normalizable) eigenfunctions of the full evolution operator, as well as the effect of truncating the computational domain.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Generalized universal instability: Transient linear amplification and subcritical turbulence

In this work we numerically demonstrate both significant transient (i.e. non-modal) linear amplification and sustained nonlinear turbulence in a kinetic plasma system with no unstable eigenmodes. The particular system considered is an electrostatic plasma slab with magnetic shear, kinetic electrons and ions, weak collisions, and a density gradient, but with no temperature gradient. In contrast to hydrodynamic examples of non-modal growth and subcritical turbulence, here there is no sheared flow in the equilibrium. Significant transient linear amplification is found when the magnetic shear and collisionality are weak. It is also demonstrated that nonlinear turbulence can be sustained if initialized at sufficient amplitude. We prove these two phenomena are related: when sustained turbulence occurs without unstable eigenmodes, states that are typical of the turbulence must yield transient linear amplification of the gyrokinetic free energy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Impurities in a non-axisymmetric plasma: transport and effect on bootstrap current

Impurities cause radiation losses and plasma dilution, and in stellarator plasmas the neoclassical ambipolar radial electric field is often unfavorable for avoiding strong impurity peaking. In this work we use a new continuum drift-kinetic solver, the SFINCS code (the Stellarator Fokker-Planck Iterative Neoclassical Conservative Solver) [M. Landreman et al., Phys. Plasmas 21 (2014) 042503] which employs the full linearized Fokker-Planck-Landau operator, to calculate neoclassical impurity transport coefficients for a Wendelstein 7-X (W7-X) magnetic configuration. We compare SFINCS calculations with theoretical asymptotes in the high collisionality limit. We observe and explain a 1/nu-scaling of the inter-species radial transport coefficient at low collisionality, arising due to the field term in the inter-species collision operator, and which is not found with simplified collision models even when momentum correction is applied. However, this type of scaling disappears if a radial electric field is present. We also use SFINCS to analyze how the impurity content affects the neoclassical impurity dynamics and the bootstrap current. We show that a change in plasma effective charge Zeff of order unity can affect the bootstrap current enough to cause a deviation in the divertor strike point locations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Less constrained omnigeneous stellarators

A stellarator is said to be omnigeneous if all particles have vanishing average radial drifts. In omnigeneous stellarators, particles are perfectly confined in the absence of turbulence and collisions, whereas in non-omnigeneous configurations, particle can drift large radial distances. One of the consequences of omnigeneity is that the unfavorable inverse scaling with collisionality of the stellarator neoclassical fluxes disappears. In the pioneering and influential article [Cary~J~R and Shasharina~S~G 1997 {\it Phys. Plasmas} {\bf 4} 3323], the conditions that the magnetic field of a stellarator must satisfy to be omnigeneous are derived. However, reference [Cary~J~R and Shasharina~S~G 1997 {\it Phys. Plasmas} {\bf 4} 3323] only considered omnigeneous stellarators in which all the minima of the magnetic field strength on a flux surface must have the same value. The same is assumed for the maxima. We show that omnigenenous magnetic fields can have local minima and maxima with different values. Thus, the parameter space in which omnigeneous stellarators are possible is larger than previously expected. The analysis presented in this article is only valid for orbits with vanishing radial width, and in principle it is not applicable to energetic particles. However, one would expect that improving neoclassical confinement would improve energetic particle confinement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Universal instability for wavelengths below the ion Larmor scale

We demonstrate that the universal mode driven by the density gradient in a plasma slab can be absolutely unstable even in the presence of reasonable magnetic shear. Previous studies from the 1970s that reached the opposite conclusion used an eigenmode equation limited to $L_x \gg ρ_i$, where $L_x$ is the scale length of the mode in the radial direction, and $ρ_i$ is the ion Larmor radius. Here we instead use a gyrokinetic approach which does not have this same limitation. Instability is found for perpendicular wavenumbers $k_y$ in the range $0.7 \lesssim k_y ρ_i \lesssim 100$, and for sufficiently weak magnetic shear: $L_s / L_n \ge 17$, where $L_s$ and $L_n$ are the scale lengths of magnetic shear and density. Thus, the plasma drift wave in a sheared magnetic field may be unstable even with no temperature gradients, no trapped particles, and no magnetic curvature.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Comparison of particle trajectories and collision operators for collisional transport in nonaxisymmetric plasmas

In this work, we examine the validity of several common simplifying assumptions used in numerical neoclassical calculations for nonaxisymmetric plasmas, both by using a new continuum drift-kinetic code and by considering analytic properties of the kinetic equation. First, neoclassical phenomena are computed for the LHD and W7-X stellarators using several versions of the drift-kinetic equation, including the commonly used incompressible-ExB-drift approximation and two other variants, corresponding to different effective particle trajectories. It is found that for electric fields below roughly one third of the resonant value, the different formulations give nearly identical results, demonstrating the incompressible ExB-drift approximation is quite accurate in this regime. However, near the electric field resonance, the models yield substantially different results. We also compare results for various collision operators, including the full linearized Fokker-Planck operator. At low collisionality, the radial transport driven by radial gradients is nearly identical for the different operators, while in other cases it is found to be important that collisions conserve momentum.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical calculation of the runaway electron distribution function and associated synchrotron emission

Synchrotron emission from runaway electrons may be used to diagnose plasma conditions during a tokamak disruption, but solving this inverse problem requires rapid simulation of the electron distribution function and associated synchrotron emission as a function of plasma parameters. Here we detail a framework for this forward calculation, beginning with an efficient numerical method for solving the Fokker-Planck equation in the presence of an electric field of arbitrary strength. The approach is continuum (Eulerian), and we employ a relativistic collision operator, valid for arbitrary energies. Both primary and secondary runaway electron generation are included. For cases in which primary generation dominates, a time-independent formulation of the problem is described, requiring only the solution of a single sparse linear system. In the limit of dominant secondary generation, we present the first numerical verification of an analytic model for the distribution function. The numerical electron distribution function in the presence of both primary and secondary generation is then used for calculating the synchrotron emission spectrum of the runaways. It is found that the average synchrotron spectra emitted from realistic distribution functions are not well approximated by the emission of a single electron at the maximum energy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Radially global $δf$ computation of neoclassical phenomena in a tokamak pedestal

Conventional radially-local neoclassical calculations become inadequate if the radial gradient scale lengths of the H-mode pedestal become as small as the poloidal ion gyroradius. Here, we describe a radially global $δf$ continuum code that generalizes neoclassical calculations to allow stronger gradients. As with conventional neoclassical calculations, the formulation is time-independent and requires only the solution of a single sparse linear system. We demonstrate precise agreement with an asymptotic analytic solution of the radially global kinetic equation in the appropriate limits of aspect ratio and collisionality. This agreement depends crucially on accurate treatment of finite orbit width effects.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

New velocity-space discretization for continuum kinetic calculations and Fokker-Planck collisions

Numerical techniques for discretization of velocity space in continuum kinetic calculations are described. An efficient spectral collocation method is developed for the speed coordinate - the radius in velocity space - employing a novel set of non-classical orthogonal polynomials. For problems in which Fokker-Planck collisions are included, a common situation in plasma physics, a procedure is detailed to accurately and efficiently treat the field term in the collision operator (in the absence of gyrokinetic corrections). When species with disparate masses are included simultaneously, a careful extrapolation of the Rosenbluth potentials is performed. The techniques are demonstrated in neoclassical calculations of the bootstrap current and plasma flows in a tokamak.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radio frequency induced and neoclassical asymmetries and their effects on turbulent impurity transport in a tokamak

Poloidal asymmetries in the impurity density can be generated by radio frequency heating in the core and by neoclassical effects in the edge of tokamak plasmas. In a pedestal case study, using global neoclassical simulations we find that finite orbit width effects can generate significant poloidal variation in the electrostatic potential, which varies on a small radial scale. Gyrokinetic modeling shows that these poloidal asymmetries can be strong enough to significantly modify turbulent impurity peaking. In the pedestal the ExB drift in the radial electric field can give a larger contribution to the poloidal motion of impurities than that of their parallel streaming. Under such circumstances we find that up-down asymmetries can also affect impurity peaking.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Changes to neoclassical flow and bootstrap current in a tokamak pedestal

In a tokamak pedestal, radial scale lengths can become comparable to the ion orbit width, invalidating conventional neoclassical calculations of flow and bootstrap current. In this work we illustrate a non-local approach that allows strong radial density variation while maintaining small departures from a Maxwellian distribution. Non-local effects alter the magnitude and poloidal variation of the flow and current. The approach is implemented in a new global delta-f continuum code using the full linearized Fokker-Planck collision operator. Arbitrary collisionality and aspect ratio are allowed as long as the poloidal magnetic field is small compared to the total magnetic field. Strong radial electric fields, sufficient to electrostatically confine the ions, are also included. These effects may be important to consider in any comparison between experimental pedestal flow measurements and theory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Local and global Fokker-Planck neoclassical calculations showing flow and bootstrap current modification in a pedestal

In transport barriers, particularly H-mode edge pedestals, radial scale lengths can become comparable to the ion orbit width, causing neoclassical physics to become radially nonlocal. In this work, the resulting changes to neoclassical flow and current are examined both analytically and numerically. Steep density gradients are considered, with scale lengths comparable to the poloidal ion gyroradius, together with strong radial electric fields sufficient to electrostatically confine the ions. Attention is restricted to relatively weak ion temperature gradients (but permitting arbitrary electron temperature gradients), since in this limit a delta-f (small departures from a Maxwellian distribution) rather than full-f approach is justified. This assumption is in fact consistent with measured inter-ELM H-Mode edge pedestal density and ion temperature profiles in many present experiments, and is expected to be increasingly valid in future lower collisionality experiments. In the numerical analysis, the distribution function and Rosenbluth potentials are solved for simultaneously, allowing use of the exact field term in the linearized Fokker-Planck collision operator. In the pedestal, the parallel and poloidal flows are found to deviate strongly from the best available conventional neoclassical prediction, with large poloidal variation of a different form than in the local theory. These predicted effects may be observable experimentally. In the local limit, the Sauter bootstrap current formulae appear accurate at low collisionality, but they can overestimate the bootstrap current near the plateau regime. In the pedestal ordering, ion contributions to the bootstrap and Pfirsch-Schluter currents are also modified.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Omnigenity as generalized quasisymmetry

Any viable stellarator reactor will need to be nearly omnigenous, meaning the radial guiding-center drift velocity averages to zero over time for all particles. While omnigenity is easier to achieve than quasisymmetry, we show here that several properties of quasisymmetric plasmas also apply directly or with only minor modification to the larger class of omnigenous plasmas. For example, concise expressions exist for the flow and current, closely resembling those for a tokamak, and these expressions are explicit in that no magnetic differential equations remain. A helicity (M,N) can be defined for any omnigenous field, based on the topology by which |B| contours close on a flux surface, generalizing the helicity associated with quasisymmetric fields. For generalized quasi-poloidal symmetry (M=0), the bootstrap current vanishes, which may yield desirable equilibrium and stability properties. A concise expression is derived for the radial electric field in any omnigenous plasma that is not quasisymmetric. The fact that tokamak-like analytical calculations are possible in omnigenous plasmas despite their fully-3D magnetic spectrum makes these configurations useful for gaining insight and benchmarking codes. A construction is given to produce omnigenous B(theta, zeta) patterns with stellarator symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The monoenergetic approximation in stellarator neoclassical calculations

In the standard "monoenergetic" approach to numerical calculation of stellarator neoclassical transport, to expedite computation, ad-hoc changes are made to the kinetic equation so speed enters only as a parameter. Here we examine the validity of this approach by considering the effective particle trajectories in a model magnetic field. We find monoenergetic codes systematically under-predict the true trapped particle fraction, with the error in the trapped ion fraction being of order unity when the electric field is large, suggesting some results of these codes may be unreliable in this regime. This inaccuracy is independent of any errors introduced by approximation of the collision operator.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Geometry and neoclassical theory in a quasi-isodynamic stellarator

We show that in perfectly quasi-isodynamic magnetic fields, which are generally non-quasisymmetric and which can approximate fields of experimental interest, neoclassical calculations can be carried out analytically more completely than in a general stellarator. Here, we define a quasi-isodynamic field to be one in which the longitudinal adiabatic invariant is a flux function and in which the constant-B contours close poloidally. We first derive several geometric relations among the magnetic field components and the field strength. Using these relations, the forms of the flow and current are obtained for arbitrary collisionality. The flow, radial electric field, and bootstrap current are also determined explicitly for the long-mean-free-path regime.

0 citations0 reviewsNo code linkedOpen access