Search

Distribution functions of collisionless systems are known to show non-thermal power law tails. Interestingly, collisionless plasmas in various physical scenarios, (e.g., the ion population of the solar wind) feature a $v^{-5}$ tail in the velocity ($v$) distribution, whose origin has been a long-standing mystery. We show this power law tail to be a natural outcome of the self-consistent collisionless relaxation of driven electrostatic plasmas. We perform a quasilinear analysis of the perturbed Vlasov-Poisson equations to show that the coarse-grained mean distribution function (DF), $f_0$, follows a quasilinear diffusion equation with a diffusion coefficient $D(v)$ that depends on $v$ through the plasma dielectric constant. If the plasma is isotropically forced on scales much larger than the Debye length with a white noise-like electric field, then $D(v)\sim v^4$ for $\sigma

A first-order model is derived for quasisymmetric stellarators where the vacuum field due to coils is dominant, but plasma-current-induced terms are not negligible and can contribute to magnetic differential equations, with $\beta$ of the order of the ratio of induced to vacuum fields. Under these assumptions, it is proven that the aspect ratio must be large and a simple expression can be obtained for the lowest-order vacuum field. The first-order correction, which involves both vacuum and current-driven fields, is governed by a Grad-Shafranov equation and the requirement that flux surfaces exist. These two equations are not always consistent, and so this model is generally overconstrained, but special solutions exist that satisfy both equations simultaneously. One family of such solutions are the first-order near-axis solutions. Thus, the first-order near-axis model is a subset of the model presented here. Several other solutions outside the scope of the near-axis model are also found. A case study comparing one such solution to a VMEC-generated solution shows good agreement., Comment: 15 pages, 3 figures. Zenodo dataset for this article available at https://doi.org/10.5281/zenodo.10674506

We consider a nearly collisionless plasma consisting of a species of `test particles' in 1D-1V, stirred by an externally imposed stochastic electric field. The mean effect on the particle distribution function is stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant $C_2 \propto \iint dx dv \, \langle f^2 \rangle$. We find that $C_2$ is transferred from large scales to small scales in both position and velocity space via a phase-space cascade enabled by both particle streaming and nonlinear interactions between particles and the stochastic electric field. We compute the steady-state fluxes and spectrum of $C_2$ in Fourier space, with $k$ and $s$ denoting spatial and velocity wavenumbers, respectively. Whereas even the linear phase mixing alone would lead to a constant flux of $C_2$ to high $s$ (towards the collisional dissipation range) at every $k$, the nonlinearity accelerates this cascade by intertwining velocity and position space so that the flux of $C_2$ is to both high $k$ and high $s$ simultaneously. Integrating over velocity (spatial) wavenumbers, the $k$-space ($s$-space) flux of $C_2$ is constant down to a dissipation length (velocity) scale that tends to zero as the collision frequency does, even though the rate of collisional dissipation remains finite. The resulting spectrum in the inertial range is a self-similar function in the $(k,s)$ plane, with power-law asymptotics at large $k$ and $s$. We argue that stochastic heating is made irreversible by this entropy cascade and that, while collisional dissipation accessed via phase mixing occurs only at small spatial scales rather than at every scale as it would in a linear system, the cascade makes phase mixing even more effective overall in the nonlinear regime than in the linear one., Comment: Minor edits, final version accepted to PRE

A simplified view of the space of optimised stellarators has the potential to guide and aid the design efforts of magnetic confinement configurations suitable for future fusion reactors. We present one such view for the class of quasisymmetric stellarators based on their approximate description near their centre (magnetic axis). The result is a space that captures existing designs and presents new ones, providing a common framework to study them. Such a simplified construction offers a basic topological approach, guided by certain theoretical and physical choices, which this paper presents in detail.

We explore the structure of the space of quasisymmetric configurations identifying them by their magnetic axes, described as 3D closed curves. We demonstrate that this topological perspective divides the space of all configurations into well-separated quasisymmetric phases. Each phase is characterized by the self-linking number (a topological invariant), defining different symmetry configurations (quasi-axisymmetry or quasi-helical symmetry). The phase-transition manifolds correspond to quasi-isodynamic configurations. By considering some models for closed curves (most notably torus unknots), general features associated with these phases are explored. Some general criteria are also built and leveraged to provide a simple way to describe existing quasisymmetric designs. This constitutes the first step in a program to identify quasisymmetric configurations with a reduced set of functions and parameters, to deepen understanding of configuration space, and offer an alternative approach to stellarator optimization that begins with the magnetic axis and builds outward.

We show that the equations satisfied by weakly quasisymmetric magnetic fields can be solved to arbitrarily high order in powers of the distance from the magnetic axis. This demonstration does not consider force balance. The existence of solutions requires an appropriate choice of parameters, most notably the toroidal current or rotational transform profiles. We do not prove that the expansion converges (it is likely divergent but asymptotic), and thus the demonstration here should not be taken as definitive proof of the existence of global solutions. Instead, we provide a systematic construction of solutions to arbitrarily high order.

We prove the existence of a straight-field-line coordinate system we call generalized Boozer coordinates. This coordinate system exists for magnetic fields with nested toroidal flux surfaces} provided $ \oint\mathrm{d}l/B\:(\mathbf{j}\cdot\nabla\psi)=0$, where symbols have their usual meaning, and the integral is taken along closed magnetic field lines. All quasisymmetric fields, regardless of their associated form of equilibria, must satisfy this condition. This coordinate system presents itself as a convenient form in which to describe general quasisymmetric configurations and their properties. Insight can be gained analytically into the difference between strong and weak forms of quasisymmetry, as well as axisymmetry, and the interaction of quasisymmetry with different forms of equilibria.

Steady plasma flows have been studied almost exclusively in systems with continuous symmetry or in open domains. In the absence of continuous symmetry, the lack of a conserved quantity makes the study of flows intrinsically challenging. In a toroidal domain, the requirement of double-periodicity for physical quantities adds to the complications. In particular, the magnetohydrodynamics (MHD) model of plasma steady-state with the flow in a non-symmetric toroidal domain allows the development of singularities when the rotational transform of the magnetic field is rational, much like the equilibrium MHD model. In this work, we show that steady flows can still be maintained provided the rotational transform is close to rational and the magnetic shear is weak. We extend the techniques developed in carrying out perturbation methods to all orders for static MHD equilibrium by Weitzner (Physics of Plasmas 21, 022515 (2014)) to MHD equilibrium with flows. We construct perturbative MHD equilibrium in a doubly-periodic domain with nearly parallel flows by systematically eliminating magnetic resonances order by order. We then utilize an additional symmetry of the flow problem, first discussed by E. Hameiri in (J. Math. Phys. \textbf{22}, 2080 (1981) Sec. III), to obtain a generalized Grad-Shafranov equation for a class of non-symmetric three-dimensional MHD equilibrium with flows both parallel and perpendicular to the magnetic field. For this class of flows, we are able to obtain non-symmetric generalizations of integrals of motion, such as Bernoulli's function and angular momentum. Finally, we obtain the generalized Hamada conditions, which are consistency conditions necessary to suppress singular currents in such a system when the magnetic field lines are closed. We do not attempt to address the question of neoclassical damping of flows., Comment: Updated and revised version

While several results have pointed to the existence of exactly quasisymmetric fields on a surface (Garren & Boozer 1991a,b; Plunk & Helander 2018), we have obtained the first such solutions using a vacuum surface expansion formalism. We obtain a single nonlinear parabolic PDE for a function $\eta$ such the field strength satisfies $B = B(\eta)$. Closed-form solutions are obtained in cylindrical, slab, and isodynamic geometries. Numerical solutions of the full nonlinear equations in general axisymmetric toroidal geometry are obtained, resulting in a class of quasi-helical local vacuum equilibria near an axisymmetric surface. The analytic models provide additional insight into general features of the nonlinear solutions, such as localization of the surface perturbations on the inboard side.

Non-symmetric vacuum magnetic fields with closed magnetic field lines are of interest in the construction of stellarator equilibria. Beyond the result of D.Lortz (ZAMP \textbf{21}, 196 (1970)), few results are available. This work presents a closed-form expression for a class of vacuum magnetic fields in a topological torus with closed field lines. We explicitly obtain the invariants of such a field. We finally show that a three-dimensional low beta magnetohydrodynamic (MHD) equilibrium may be constructed in a topological torus starting with these closed line vacuum magnetic fields.

A method is given to rapidly compute quasisymmetric stellarator magnetic fields for plasma confinement, without the need to call a three-dimensional magnetohydrodynamic equilibrium code inside an optimization iteration. The method is based on direct solution of the equations of magnetohydrodynamic equilibrium and quasisymmetry using Garren and Boozer's expansion about the magnetic axis (Phys Fluids B 3, 2805 (1991)), and it is several orders of magnitude faster than the conventional optimization approach. The work here extends the method of Landreman, Sengupta and Plunk (J Plasma Phys 85, 905850103 (2019)), which was limited to flux surfaces with elliptical cross-section, to higher order in the aspect ratio expansion. As a result, configurations can be generated with strong shaping that achieve quasisymmetry to high accuracy. Using this construction, we give the first numerical demonstrations of Garren and Boozer's ideal scaling of quasisymmetry-breaking with the cube of inverse aspect ratio. We also demonstrate a strongly nonaxisymmetric configuration (vacuum $\iota > 0.4$) in which symmetry-breaking mode amplitudes throughout a finite volume are $< 2\times 10^{-7}$, the smallest ever reported. To generate boundary shapes of finite-minor-radius configurations, a careful analysis is given of the effect of substituting a finite minor radius into the near-axis expansion. The approach here can provide analytic insight into the space of possible quasisymmetric stellarator configurations, and it can be used to generate good initial conditions for conventional stellarator optimization.

We carry out expansions of non-symmetric toroidal ideal magnetohydrodynamic (MHD) equilibria with nested flux surfaces about a periodic cylinder, in which physical quantities are periodic of period $2\pi$ in the cylindrical angle $\theta$ and $z$. The cross-section of a flux surface at a constant toroidal angle is assumed to be approximately circular, and data is given on the cylindrical flux surface $r=1$. Furthermore, we assume that the magnetic field lines are closed on the lowest order flux surface, and the magnetic shear is relatively small. We extend earlier work in a flat-torus by Weitzner [Physics of Plasmas 23, 062512 (2016)]and demonstrate that a power series expansion can be carried out to all orders using magnetic flux as an expansion parameter. The cylindrical metric introduces certain new features to the expansions compared to the flat-torus. However, the basic methodology of dealing with resonance singularities remains the same. The results, even though lacking convergence proofs, once again support the possibility of smooth, low-shear non-symmetric toroidal MHD equilibria., Comment: Part of NYU GSTEM, a program to encourage girls in high school to take up STEM subjects

Quasisymmetric stellarators are appealing intellectually and as fusion reactor candidates since the guiding center particle trajectories and neoclassical transport are isomorphic to those in a tokamak, implying good confinement. Previously, quasisymmetric magnetic fields have been identified by applying black-box optimization algorithms to minimize symmetry-breaking Fourier modes of the field strength $B$. Here instead we directly construct magnetic fields in cylindrical coordinates that are quasisymmetric to leading order in distance from the magnetic axis, without using optimization. The method involves solution of a 1-dimensional nonlinear ordinary differential equation, originally derived by Garren and Boozer [Phys. Fluids B 3, 2805 (1991)]. We demonstrate the usefulness and accuracy of this optimization-free approach by providing the results of this construction as input to the codes VMEC and BOOZ_XFORM, confirming the purity and scaling of the magnetic spectrum. The space of magnetic fields that are quasisymmetric to this order is parameterized by the magnetic axis shape along with three other real numbers, one of which reflects the on-axis toroidal current density, and another one of which is zero for stellarator symmetry. The method here could be used to generate good initial conditions for conventional optimization, and its speed enables exhaustive searches of parameter space.

The confinement of guiding center trajectories in a stellarator is determined by the variation of the magnetic field strength $B$ in Boozer coordinates $(r, \theta, \varphi)$, but $B(r,\theta,\varphi)$ depends on the flux surface shape in a complicated way. Here we derive equations relating $B(r,\theta,\varphi)$ in Boozer coordinates and the rotational transform to the shape of flux surfaces in cylindrical coordinates, using an expansion in distance from the magnetic axis. A related expansion was done by Garren and Boozer [Phys. Fluids B 3, 2805 (1991)] based on the Frenet-Serret frame, which can be discontinuous anywhere the magnetic axis is straight, a situation that occurs in the interesting case of omnigenity with poloidally closed $B$ contours. Our calculation in contrast does not use the Frenet-Serret frame. The transformation between the Garren-Boozer approach and cylindrical coordinates is derived, and the two approaches are shown to be equivalent if the axis curvature does not vanish. The expressions derived here help enable optimized plasma shapes to be constructed that can be provided as input to VMEC and other stellarator codes, or to generate initial configurations for conventional stellarator optimization.

Stellarators are generically small current and low plasma beta ($\beta= p/B^2\ll1$) devices. Often the construction of vacuum magnetic fields with good magnetic surfaces is the starting point for an equilibrium calculation. Although in cases with some continuous spatial symmetry flux functions can always be found for vacuum magnetic fields, an analogous function does not, in general, exist in three dimensions. This work examines several simple equilibria and vacuum magnetic field problems with the intent of demonstrating the possibilities and limitations in the construction of such states. Starting with a simple vacuum magnetic field with closed field lines in a topological torus (toroidal shell with a flat metric), we obtain a self-consistent formal perturbation series using the amplitude of the nonsymmetric vacuum fields as a small parameter. We show that systems possessing stellarator symmetry allow the construction order by order. We further indicate the significance of stellarator symmetry in the amplitude expansion of the full ideal magnetohydrodynamics (MHD) problem as well. We then investigate the conditions that guarantee neighboring flux surfaces given the data on one surface, by expanding in the distance from that surface. We show that it is much more difficult to find low shear vacuum fields with surfaces than force-free fields or ideal MHD equilibrium. Finally, we demonstrate the existence of a class of vacuum magnetic fields, analogous to `snakes' observed in tokamaks, which can be expanded to all orders.

Quasisymmetry and omnigeneity of an equilibrium magnetic field are two distinct properties proposed to ensure radial localization of collisionless trapped particles in any stellarator. These constraints are incompletely explored, but have stringent restrictions on a magnetic geometry. This work employs an analytic approach to understand the implications of the constraints. The particles move in an intrinsically three dimensional equilibrium whose representation is given by earlier work of Weitzner and its extension here. For deeply trapped particles a local equilibrium expansion around a minimum of the magnetic field strength along a magnetic line suffices. This analytical non-symmetric equilibrium solution, enables explicit representation of the constraints. The results show that it is far easier to satisfy the omnigeneity condition than the quasisymmetry requirement. Correspondingly, there exists a large class of equilibrium close to quasisymmetry that remain omnigeneous while allowing inclusion of error fields, which may destroy quasisymmetry.

A new analytical operator method is discussed which solves linear ordinary differential equations with regular singularities. Solutions are obtained in analytic series form and also in Mellin-Barnes-type contour integral form. Exact series solution is obtained without having to calculate series coefficients by recurrence relation.Both homogeneous and inhomogeneous equations are solved identically without having to calculate the Green's function explicitly in the case of inhomogeneous equation.Closed-form solutions are obtained for all the special functions appearing in mathematical physics. For a second-order equation both the independent solutions are obtained without invoking Wronskians, even when the indices differ by an integer., Comment: 14 pages 0 figures