SEARCH · Engineering Papers

Results for “Near-axis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Extending near-axis equilibria in $\mathtt{DESC}$

The near-axis description of optimised stellarator fields has proven to be a powerful tool both for the design and understanding of this magnetic confinement concept. The description consists of an asymptotic model of the equilibrium in the distance from its centremost axis, and is thus only approximate. Any practical application therefore requires the eventual construction of a global equilibrium. This paper presents a novel way of constructing global equilibria using the DESC code that guarantees the correct asymptotic behaviour imposed by a given near-axis construction. The theoretical underpinnings of this construction are carefully presented, and benchmarking examples provided. This opens the door to an efficient coupling of the near-axis framework and that of global equilibria for future optimisation efforts.

fusion plasma↗

DOE OSTI2025

Inputs, Outputs and Plotting Scripts for paper Extending near-axis equilibria in DESC

Python scripts and outputs from the DESC and pyQSC/pyQIC codes used to create the results in the paper "Extending near-axis equilibria in DESC". The contained README file has the details of which scripts create which figures, as well as on what versions of the codes were used. The repo also contains python pickle (.p) files and .txt files with the data used to create each figure (which are used by the plotting scripts).

DESC↗

DOE OSTI2026

Crystal Orientation and Defect Mapping in Electron-Beam-Sensitive Zeolites with Near-Axis Transmission Kikuchi Diffraction

Porous materials are vital in catalysis, energy conversion, and environmental remediation. Understanding structural heterogeneity in zeolites is key to linking synthesis, framework intergrowths, and catalytic performance, yet current methods for phase identification and spatial mapping lack sufficient resolution or throughput. We present a high-throughput approach using near-axis transmission Kikuchi diffraction in a scanning electron microscope, achieving high phase and spatial resolution for electron-beam-sensitive zeolites, including ZSM-5 and, for the first time, Zeolite A. Here, this method enables direct visualization of intergrowth features that critically affect catalytic and adsorption behavior, bridging the gap between ensemble-averaged X-ray diffraction and high-resolution but low-throughput transmission electron microscopy. Combining nanoscale mapping with statistical sampling is highly suited for machine-learning pipelines guiding new structure function understanding and could be extended to other beam-sensitive porous materials such as metal–organic or covalent–organic frameworks.

crystallography↗

DOE OSTI2022

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 & Boozer [ Phys. Fluids B, vol. 3, 1991, p. 2805]. Due to the reduction of the equations by the expansion, the computational cost is significantly reduced, to times of 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 $\| \boldsymbol {\nabla }\boldsymbol B \|$ , the norm of the magnetic field gradient, 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.

Landreman, Matt (ORCID:000000027233577X)↗

DOE OSTI2023

Constructing the space of quasisymmetric stellarators through near-axis expansion

Abstract 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

An asymptotic Grad–Shafranov equation for quasisymmetric stellarators

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 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 is the set of 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.

Nikulsin, Nikita (ORCID:0000000318611777)↗

DOE OSTI2023

Direct stellarator coil optimization for nested magnetic surfaces with precise quasi-symmetry

We present a robust optimization algorithm for the design of electromagnetic coils that generate vacuum magnetic fields with nested flux surfaces and precise quasi-symmetry. The method is based on a bilevel optimization problem, where the outer coil optimization is constrained by a set of inner least squares optimization problems whose solutions describe magnetic surfaces. The outer optimization objective targets coils that generate a field with nested magnetic surfaces and good quasi-symmetry. The inner optimization problems identify magnetic surfaces when they exist, and approximate surfaces in the presence of magnetic islands or chaos. We show that this formulation can be used to heal islands and chaos, thus producing coils that result in magnetic fields with precise quasi-symmetry. We show that the method can be initialized with coils from the traditional two-stage coil design process, as well as coils from a near-axis expansion optimization. We present a numerical example where island chains are healed and quasi-symmetry is optimized up to surfaces with aspect ratio 6. Another numerical example illustrates that the aspect ratio of nested flux surfaces with optimized quasi-symmetry can be decreased from 6 to approximately 4. In conclusion, the last example shows that our approach is robust and a cold-start using coils from a near-axis expansion optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Application of Lagrangian techniques for calculating the on-axis rotational transform

The Floquet exponents of periodic field lines are studied through the variations of the magnetic action on the magnetic axis, which is assumed to be elliptical. The near-axis formalism developed by Mercier, Solov'ev and Shafranov is combined with a Lagrangian approach. The on-axis Floquet exponent is shown to coincide with the on-axis rotational transform. A discrete solution suitable for numerical implementation is introduced, which gives the Floquet exponents as solutions to an eigenvalue problem. This discrete formalism expresses the exponents as the eigenvalues of a $6$ X $6$ matrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Minimizing separatrix crossings through isoprominence

Abstract A simple property of magnetic fields that minimizes bouncing to passing type transitions of guiding center orbits is defined and discussed. This property, called isoprominence, is explored through the framework of a near-axis expansion. It is shown that isoprominent magnetic fields for a toroidal configuration exist to all orders in a formal expansion about a magnetic axis. Some key geometric features of these fields are described.

confinement↗

DOE OSTI2024

Stellarator equilibrium axis-expansion to all orders in distance from the axis for arbitrary plasma beta

A systematic theory of the asymptotic expansion of the magnetohydrostatics (MHS) equilibrium in the distance from the magnetic axis is developed to include arbitrary smooth currents near the magnetic axis. Compared with the vacuum and the force-free system, an additional magnetic differential equation must be solved to obtain the pressure-driven currents. It is shown that there exist variables in which the rest of the MHS system closely mimics the vacuum system. Thus, a unified treatment of MHS fields is possible. The mathematical structure of the near-axis expansions to arbitrary order is examined carefully to show that the double-periodicity of physical quantities in a toroidal domain can be satisfied order by order. The essential role played by the leading-order Birkhoff–Gustavson normal form in solving the magnetic differential equations is highlighted. Several explicit examples of vacuum, force-free and MHS equilibrium in different geometries are presented.

Physics↗

DOE OSTI2025

The shear Alfvén continuum of quasisymmetric stellarators

The shear Alfvén wave (SAW) continuum plays a critical role in the stability of energetic particle-driven Alfvén eigenmodes. We develop a theoretical framework to analyze the SAW continuum in three-dimensional (3-D) quasisymmetric magnetic fields, focusing on its implications for stellarator design. By employing a near-axis model and degenerate perturbation theory, the continuum equation is solved, highlighting unique features in 3-D configurations, such as the interactions between spectral gaps. Numerical examples validate the theory, demonstrating the impact of flux-surface shaping and quasisymmetric field properties on continuum structure. The results provide insights into optimizing stellarator configurations to minimize resonance-driven losses of energetic particles. This work establishes a basis for incorporating Alfvénic stability considerations into the stellarator design process, demonstrated through optimization of a quasihelical configuration to avoid high-frequency spectral gaps.

fusion plasma↗

DOE OSTI2023

Nonlinear magnetohydrodynamic modeling of current-drive-induced sawtooth-like crashes in the W7-X stellarator

Sawtooth-like core electron temperature crashes have been observed in W7-X experiments with electron cyclotron current drive. Here we present nonlinear single-fluid magnetohydrodynamic simulations of this phenomenon using the newly developed stellarator modeling capability of the M3D-C 1 code. The near-axis current drive gives rise to two ι = 1 resonances in the equilibrium rotational transform profile so that two consecutive (1, 1) internal kink modes are seen in the simulations. A small-amplitude crash at the inner resonance occurs first, which may correspond to the sawtooth precursors observed in the experiments. A bigger crash at the outer resonance then flattens the core temperature profile, which shows semi-quantitative agreements with experimental measurements on certain metrics such as the crash amplitude and the inversion radius of the temperature change. These results illustrate a likely mechanism of the current-drive-induced sawtooth-like crashes in W7-X and, to some extent, validate the stellarator modeling capability of M3D-C 1 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Orbit-averaging and deposition accuracy for runaway electron beams in hybrid kinetic-MHD simulations of the runaway plateau

We develop a new procedure that combines the kinetic orbit runaway electron code (KORC) and the NIMROD extended-magnetohydrodynamic code to simulate runaway electrons (REs) in the post-disruption plateau. KORC integrates guiding-center orbits, with a barycentric-based binary search strategy providing initial guesses for the Newton–Raphson logical-to-physical coordinate inversion, ensuring reliable particle-to-mesh mapping in NIMROD, whose fields remain static for the present study. Samples are drawn in accord with experimental parallel current profiles of RE beams during the plateau phase. Deposition in NIMROD is verified through comparison with a Python-based finite-element code that ensures periodicity in the poloidal direction and continuity at the magnetic axis. Accurate representation of near-axis fields requires finer mesh resolution to prevent under- and overshoots in current density from orbit inaccuracies. Yet, at a fixed particle count, increasing mesh resolution amplifies statistical noise in the deposited fields. An orbit-averaging method accumulates partial current deposits over multiple kinetic steps and reduces the statistical noise with little added computational cost. By coupling kinetic routines from KORC directly into the NIMROD codebase, these developments lay essential groundwork for future self-consistent KORC–NIMROD coupling.

Algorithms and data structure↗

DOE OSTI2025

Alfvén mode spectrum in the advanced steady-state tokamak scenario with negative triangularity

In the earlier studies [Zheng et al. Fundam. Plasma Phys. 10, 100051 (2024)], the negative triangularity tokamak was shown to have a unique and outstanding feature for the steady-state confinement. In addition to relieving the divertor heat load concern, it can reach a much higher beta limit in the nearly steady-state confinement, about twice the Troyon limit in the positive triangularity case. In this work, we extend the investigation to the Alfvén spectrum and toroidal Alfvén eigenmodes (TAEs). This helps to evaluate the property of energetic particle confinement in this advanced scenario of a negative triangularity tokamak. We found that, because of the large safety factor and reversed magnetic shear in the core plasma region near the magnetic axis, the TAE gap tends to be closed or have only a short radial extension between the continua. in the on-axis core region and open primarily in the off-axis region, where TAEs are shown to be present. This indicates that energetic particle confinement in the core region near the magnetic axis, where nuclear fusion reactions prevail, is not seriously affected by the excitation of TAEs in this scenario. The configuration with off-axis TAEs helps the confinement of energetic particles as compared to the configuration with near-axis TAEs. The results show that the advanced steady-state tokamak scenario with negative triangularity, because of the rapid variation of the safety factor, also has a favorable feature for energetic particle confinement.

Fourier analysis↗

DOE OSTI2025

Periodic Korteweg–de Vries soliton potentials generate quasisymmetric magnetic field strength in a finite plasma- β equilibrium

Quasisymmetry (QS) is a hidden symmetry of the magnetic field strength, B , that enables effective confinement of charged particles in a fully three-dimensional (3D) toroidal plasma equilibrium. Such equilibria are typically modeled by the ideal magnetohydrostatic (MHS) equations. The nonlinear, overdetermined nature of the QS MHS equations severely complicates our understanding of the interplay between 3D shaping, equilibrium properties such as pressure and rotational transform, and B . Progress has been made through expansions near the magnetic axis; however, a more comprehensive theory is desirable. Using a combination of analysis and regression on a large dataset of numerically optimized quasisymmetric stellarators, we demonstrate that there is a hidden lower dimensionality of B on a magnetic flux surface with connections to the theory of periodic solitons. We show that B on a flux surface is determined by three or at most four flux functions, each of which determines a critical value of the derivative of B along the field line. While being consistent with the near-axis models, our results are global and hold even on the last closed flux surface.

Differential geometry↗

DOE OSTI2025

2D photofragmentation LIF imaging of H 2 O 2 and HO 2 in the effluent of an atmospheric-pressure plasma jet: effects of solid and liquid interfaces

Two-dimensional (2D) absolute measurements of hydrogen peroxide (H 2 O 2 ) and approximations of the hydroperoxyl radical (HO 2 ) in the effluent of a COST Reference Microplasma Jet operated with a He/H 2 O feed gas are presented. Gas-phase densities are mapped using photofragmentation laser-induced fluorescence (PF-LIF) under three boundary conditions: open effluent, a solid target, and a liquid target. A novel method is presented for separating PF-LIF signals from H 2 O 2 and HO 2 using comparative measurements in oxygen-rich and oxygen-free environments to exploit the preferential formation of HO 2 in the presence of molecular oxygen. This separation strategy is supported by results from a plug-flow plasma chemistry model. Measured densities agree closely with model predictions in both magnitude and trend, while the 2D experimental distributions provide additional insight into the spatial dependencies of these species. In particular, the results show distinct differences in species transport depending on the target type: solid surfaces induce lateral deflection and reduced centerline densities, whereas liquid interfaces promote axial accumulation and higher near-axis concentrations.

atmospheric-pressure plasma jet (APPJ)↗

DOE OSTI2024

Enhancement of harmonic generation by an intense driving laser with high-order waveguide modes in a high-pressure gas-filled hollow waveguide

High-order harmonics have been widely used as reliable tabletop coherent radiation sources recently, but their applications have often been limited by the available pulse energy. Here, we report that by using an overdriven intense laser in a long waveguide with high-pressure gas, phase matching can be achieved in three distinct “regimes”. In the third regime, favorable phase matching is achieved at near-axis positions to enhance harmonic yields. Our results are supported by a full theoretical analysis, and we demonstrate that coupling of the driving laser with the high-order waveguide modes (instead of the fundamental mode used in most prior experiments) is responsible for achieving phase matching. Furthermore, we establish that this phase matching (and harmonic enhancement) is robust, and a scaling relation is derived for the necessary waveguide and gas parameters, allowing our predictions to be tested immediately in any laboratory today.

Li, Baochang [Nanjing University of Science and Te↗

NASA NTRS1994

Speckle imaging with the MAMA detector: Preliminary results

We report on the first successful speckle imaging studies using the Stanford University speckle interferometry system, an instrument that uses a multianode microchannel array (MAMA) detector as the imaging device. The method of producing high-resolution images is based on the analysis of so-called 'near-axis' bispectral subplanes and follows the work of Lohmann et al. (1983). In order to improve the signal-to-noise ratio in the bispectrum, the frame-oversampling technique of Nakajima et al. (1989) is also employed. We present speckle imaging results of binary stars and other objects from V magnitude 5.5 to 11, and the quality of these images is studied. While the Stanford system is capable of good speckle imaging results, it is limited by the overall quantum efficiency of the current MAMA detector (which is due to the response of the photocathode at visible wavelengths and other detector properties) and by channel saturation of the microchannel plate. Both affect the signal-to-noise ratio of the power spectrum and bispectrum.

Horch, E.↗