Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Toroid”

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.

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

At least 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A fast, accurate and easy to implement Kapur–Rokhlin quadrature scheme for singular integrals in axisymmetric geometries

Many applications in magnetic confinement fusion require the efficient calculation of surface integrals with singular integrands. The singularity subtraction approaches typically used to handle such singularities are complicated to implement and low-order accurate. In contrast, we demonstrate that the Kapur–Rokhlin quadrature scheme is well-suited for the logarithmically singular integrals encountered for a toroidally axisymmetric confinement system, is easy to implement and is high-order accurate. As an illustration, we show how to apply this quadrature scheme for the efficient and accurate calculation of the normal component of the magnetic field due to the plasma current on the plasma boundary, via the virtual-casing principle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Energetic electron transport in magnetic fields with island chains and stochastic regions

This paper investigates energetic electron transport in magnetized toroidal plasmas with magnetic fields characterized by island chains and regions of stochastic field lines produced by coil perturbations. We report on experiments performed in the DIII-D tokamak, which utilize electron cyclotron heating and current drive pulses to ‘tag’ electron populations within different locations across the discharge. The cross-field transport of these populations is then inferred from electron cyclotron emission measurements and gamma emission signals from scintillator detectors. Two types of energetic particles are distinguished and discussed: non-relativistic suprathermal electrons and relativistic runaway electrons. The magnetic field topology in each discharge is reconstructed with field-line tracing codes, which are also used to determine the location and scale of magnetic islands and stochastic regions. Comparison of simulations and experiments suggests that suprathermal transport is suppressed when the tagging is performed at a smaller radial location than the location of the $q = 1$ island chain and enhanced otherwise. Here q is the safety factor. We further demonstrate that increasing the width of the stochastic region within the edge plasma yields enhancement of the suprathermal electron transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Turbulence spreading effects on the ELM size and SOL width

BOUT++ turbulence simulations were performed to investigate the impact of turbulence spreading on the edge localized mode (ELM) size and divertor heat flux width (λ q ) broadening in small ELM regimes. Here, this study is motivated by EAST experiments. BOUT++ linear simulations of a pedestal radial electric field (E r ) scan show that the dominant toroidal number mode (n) shifts from high-n to low-n, with a narrow mode spectrum, and the maximum linear growth rate increases as the pedestal E r well deepens. The nonlinear simulations show that as the net E × B pedestal flow increases, the pressure fluctuation level and its inward penetration beyond the top of the pedestal both increase. This leads to a transition from small ELMs to large ELMs. Both inward and outward turbulence spreading are sensitive to the scrape-off-layer (SOL) plasma profiles. The inward turbulence spreading increases for the steep SOL profiles, leading to increasing pedestal energy loss in the small ELM regime. The SOL width (λ q ) is significantly broadened progressing from the ELM-free to small ELM regime, due to the onset of strong radial turbulent transport. The extent of the SOL width (λ q ) broadening depends strongly on outward turbulence spreading. The fluctuation energy intensity flux Γ ε at the separatrix can be enhanced by increasing either pedestal E r flow shear or local SOL pressure gradient. The λ q is broadened as the fluctuation energy intensity flux Γ ε at the last close flux surface (LCFS) increases. Local SOL E × B flow shear will restrain outward turbulence spreading and the associated heat flux width broadening. Operating in H-mode with small ELMs has the potential to solve two critical problems: reducing the ELM size and broadening the SOL width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Global fluid turbulence simulations in the scrape-off layer of a stellarator island divertor

Isothermal fluid turbulence simulations have been performed in the edge and scrape-off layer (SOL) of an analytic stellarator configuration with an island divertor, thereby providing numerical insight into edge turbulence in regions around islands in a stellarator. The steady-state transport follows the a curvature drive that is inverse to the major radius ( $1/R$ ) toward the outboard side, but large fluctuations are present throughout the island divertor region, with the average wavelength of similar size to the island width. The system exhibits a prominent $m=2$ , $n=5$ mode, where m is the poloidal mode number and n is the toroidal mode number, although other modes are present. The amplitude and radial extent of the density fluctuations are similar throughout the edge and SOL, but can decrease near island O-points. The fluctuations exhibit a predominantly positive skewness on the outboard midplane, indicating blob-like perturbations for the transport into the outer SOL. It is determined that a point on the separatrix is generally more correlated with regions outside of the SOL than a nearby reference point which does not lie on the separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Direct, simple and efficient computation of all components of the virtual-casing magnetic field in axisymmetric geometries with Kapur–Rokhlin quadrature

In a recent publication (Toler et al., J. Plasma Phys., vol. 89, issue 2, 2023, p. 905890210), we demonstrated that for axisymmetric geometries, the Kapur–Rokhlin quadrature rule provided an efficient and high-order accurate method for computing the normal component, on the plasma surface, of the magnetic field due to the toroidal current flowing in the plasma, via the virtual-casing principle. The calculation was indirect, as it required the prior computation of the magnetic vector potential from the virtual-casing principle, followed by the computation of its tangential derivative by Fourier differentiation, to obtain the normal component of the magnetic field. Our approach did not provide the other components of the virtual-casing magnetic field. In this letter, we show that a more direct and more general approach is available for the computation of the virtual-casing magnetic field. The Kapur–Rokhlin quadrature rule accurately calculates the principal value integrals in the expression for all the components of the magnetic field on the plasma boundary, and the numerical error converges at a rate nearly as high as the indirect method we presented previously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Coil optimization for quasi-helically symmetric stellarator configurations

Filament-based coil optimizations are performed for several quasi-helical stellarator configurations, beginning with the one from Landreman & Paul ( Phys. Rev. Lett. , vol. 128, 2022, 035001), demonstrating that precise quasi-helical symmetry can be achieved with realistic coils. Several constraints are placed on the shape and spacing of the coils, such as low curvature and sufficient plasma–coil distance for neutron shielding. The coils resulting from this optimization have a maximum curvature 0.8 times that of the coils of the Helically Symmetric eXperiment (HSX) and a mean squared curvature 0.4 times that of the HSX coils when scaled to the same plasma minor radius. When scaled up to reactor size and magnetic field strength, no fast particle losses were found in the free-boundary configuration when simulating 5000 alpha particles launched at $3.5\,\mathrm {MeV}$ on the flux surface with a normalized toroidal flux of $s=0.5$ . An analysis of the tolerance of the coils to manufacturing errors is performed using a Gaussian process model, and the coils are found to maintain low particle losses for smooth, large-scale errors up to amplitudes of approximately $0.15\,\mathrm {m}$ . Another coil optimization is performed for the Landreman–Paul configuration with the additional constraint that the coils are purely planar. Visual inspection of the Poincaré plot of the resulting magnetic field-lines reveal that the planar modular coils alone do a poor job of reproducing the target equilibrium. Additional non-planar coil optimizations are performed for the quasi-helical configuration with $5\,\%$ volume-averaged plasma beta from Landreman et al. ( Phys. Plasma , vol. 29, issue 8, 2022, 082501), and a similar configuration also optimized to satisfy the Mercier criterion. The finite beta configurations had larger fast-particle losses, with the free-boundary Mercier-optimized configuration performing the worst, losing approximately $5.5\,\%$ of alpha particles launched at $s=0.5$ .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2024

Radiative and continuum dampings of reversed shear Alfvén eigenmodes and perturbative analysis limitations for tokamaks

A careful theoretical analysis of the excitation of Alfvén eigenmodes (AEs), such as TAE (toroidicity-induced AE) and RSAE (reversed shear AE), by superalfvenic energetic particles is required for reliable predictions of energetic ion relaxation in present day fusion experiments. This includes the evaluation of different AE damping mechanisms including radiative and continuum dampings which are the focus of this study. A recent comprehensive benchmark of different eigenmode solvers including gyrokinetic, gyrofluid and hybrid magenetohydrodynamics (MHD) has shown that employed models may have deficiencies when addressing some of them (Taimourzadeh et al., Nucl. Fusion, vol. 59, 2019, 066006). Here, in this paper, we are studying the radiative and continuum dampings of RSAEs in details which were missing in hybrid NOVA/NOVA-C calculations to prepare a NOVA-C package with a substantial upgrade. Both dampings require the finite Larmor radius (FLR) corrections to AE mode structures to be accounted for. Accurately calculating different damping rates and understanding their parametric dependencies, we resolve the limitation coming out of the perturbative approach. In particular, here, the radiative damping is included perturbatively, whereas the continuum damping is computed non-perturbatively. Our comparison leads to the conclusion that the non-perturbative treatment of the unstable RSAE modes is needed to find the agreement with the gyrokinetic calculations. We expect that the RSAE mode structure modification plays a dominant role in determining the RSAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Turbulence and transport in mirror geometries in the Large Plasma Device

Thanks to advances in plasma science and enabling technology, mirror machines are being reconsidered for fusion power plants and as possible fusion volumetric neutron sources. However, cross-field transport and turbulence in mirrors remains relatively understudied compared with toroidal devices. Turbulence and transport in mirror configurations were studied utilizing the flexible magnetic geometry of the Large Plasma Device (LAPD). Multiple mirror ratios from $M=1$ to $M=2.68$ and three mirror-cell lengths from $L=3.51$ to $L=10.86$ m were examined. Langmuir and magnetic probes were used to measure profiles of density, temperature, potential and magnetic field. The electric field-fluctuation-driven ${\tilde {\boldsymbol{E}}} \times {\boldsymbol{B}}$ particle flux, where $\boldsymbol{B}$ is the background field, was calculated from these quantities. Two probe correlation techniques were used to infer wavenumbers and two-dimensional structure. Cross-field particle flux and density fluctuation power decreased with increased mirror ratio. Core density and temperatures remain similar with mirror ratio, but radial line-integrated density increased. The physical expansion of the plasma in the mirror cell by using a higher field in the source region may have led to reduced density fluctuation power through the increased gradient scale length. This increased scale length reduced the growth rate and saturation level of rotational interchange and drift-like instabilities. Despite the introduction of magnetic curvature, no evidence of mirror-driven instabilities – interchange, velocity space or otherwise – were observed. For curvature-induced interchange, many possible stabilization mechanisms were present, suppressing the visibility of the instability.

Travis, Phil (ORCID:0000000343170432)↗
DOE OSTI2025

Bootstrap current modeling in M3D-C1

Bootstrap current plays a crucial role in the equilibrium of magnetically confined plasmas, particularly in quasi-symmetric stellarators and in tokamaks, where it can represent bulk of the electric current density. Accurate modeling of this current is essential for understanding the magnetohydrodynamic (MHD) equilibrium and stability of these configurations. This study expands the modeling capabilities of M3D-C1, an extended-MHD code, by implementing self-consistent physics models for bootstrap current. It employs two analytical frameworks: a generalized Sauter model (Sauter et al. 1999 Phys. Plasmas vol. 6, no. 7, pp. 2834–2839), and a revised Sauter-like model (Redl et al. 2021 Phys. Plasmas vol. 28, no. 2, pp. 022502). The isomorphism described by Landreman et al. (2022 Phys. Rev. Lett. vol. 128, pp. 035001) is employed to apply these models to quasi-symmetric stellarators. The implementation in M3D-C1 is benchmarked against neoclassical codes, including NEO, XGCa and SFINCS, showing excellent agreement. These improvements allow M3D-C1 to self-consistently calculate the neoclassical contributions to plasma current in axisymmetric and quasi-symmetric configurations, providing a more accurate representation of the plasma behavior in these configurations. A workflow for evaluating the neoclassical transport using SFINCS with arbitrary toroidal equilibria calculated using M3D-C1 is also presented. This workflow enables a quantitative evaluation of the error in the Sauter-like model in cases that deviate from axi- or quasi-symmetry (e.g. through the development of an MHD instability).

fusion plasma↗
DOE OSTI2026

Omnigenous umbilic stellarators

To better understand the dependence of the magnetic field structure in the plasma edge on the plasma boundary shape, in the context of X-point and island divertor designs, we define and develop a class of stellarators called umbilic stellarators. These equilibria are characterised by a single continuous high-curvature edge on the plasma boundary that goes around multiple times toroidally before meeting itself. We develop a technique that allows us to simultaneously optimise the plasma boundary along with a curve lying on the boundary on which we impose a high curvature while imposing omnigenity – a property of the magnetic field that ensures trapped particle confinement throughout the plasma volume. We find that umbilic stellarators naturally tend to favour piecewise omnigenity instead of omnigenity with a specific helicity. After generating omnigenous umbilic stellarators, we design coil sets for some of them and explore the fieldline structure in the edge and its sensitivity to small fluctuations in the plasma. Finally, using single-stage optimisation, we simultaneously modify the plasma and coil shape and propose an experiment to modify an existing tokamak to a finite-𝛽 stellarator using this technique and explore a potentially simpler way to convert a limited tokamak into a diverted stellarator.

fusion plasma↗
DOE OSTI2020

Structural changes in liquid lithium under high-pressure

In this paper, we have experimentally studied the effect of compression on the structure of liquid lithium (Li) by multi-angle energy dispersive X-ray diffraction in a large-volume cupped-Drickamer-Toroidal cell. The structure factors, s(q), of liquid Li have been successfully determined under an isothermal compression at 600 ± 30 K and at pressures up to 11.5 GPa. The first peak position in s(q) is found to increase with increasing pressure and showing an obvious slope change starting at ~7.5 GPa. The slope change is interpreted as a structural change from bcc-like to fcc-like local ordering in liquid Li. At pressures above 8.7 GPa, the liquid Li becomes predominantly fcc-like up to the highest pressure of 11.5 GPa in this study. The observed structural changes in liquid Li are consistent with the recently determined melting curve of Li.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Mechanistic Study of the Li–Air Battery with a Co 3 O 4 Cathode and Dimethyl Sulfoxide Electrolyte

The lithium–air battery, a powerful competitor to replace the traditional lithium-ion battery, has attracted increasing attention due to its extremely high theoretical energy density. However, its development is limited by the cathode and electrolyte properties, which should include high stability, conductivity, and electrocatalytic properties in oxygen-rich environments. Here, we employ a systematic first-principles study of Li–O 2 discharge and charge reactions on the Co 3 O 4 -based cathode with the assistance of dimethyl sulfoxide (DMSO) electrolyte. The structure, stability, and electronic properties of different surface reconstructions of the Co 3 O 4 (100) facet are investigated. In addition, the mechanisms and thermodynamic overpotentials of multi-step reactions between Li + /e – and O 2 are provided, where lithium suboxide products (Li 2 O 2 or Li 3 O 2 ) are formed on the different Co 3 O 4 (100) terminations. The solvation shell of Li + components in explicit DMSO solvent is investigated through ab initio molecular dynamics simulations. In general, we find that the Co 3 O 4 (100)-O (oxidized) surface is the most stable one under standard conditions, and the stable Li + solvation structure is found in a tetrahedral Li(DMSO) 4 + shell in the DMSO-based electrolyte. Moreover, in the system of the Co 3 O 4 (100)-O cathode and DMSO electrolyte, the solution model pathway is energetically favorable for the Li–O 2 discharge reaction. It provides a low constant overpotential of 0.17 V during a long-term discharging process, thus causing the final toroid Li 2 O 2 formation on the cathode. During the charging process, an overpotential of 0.36 V is required to rapidly decompose Li 2 O 2 .

25 ENERGY STORAGE↗
DOE OSTI2022

Observed and Bin Model Simulated Evolution of Drop Size Distributions in High-Based Cumulus Congestus Over the United Arab Emirates

This study examines microphysical processes in developing high-based cumulus congestus over the United Arab Emirates using aircraft observations and a large-eddy-simulation model with bin microphysics. A notable feature of this case is the lack of mm-sized drops despite having a liquid cloud layer >3.5 km deep, contrasting with copious large drops observed in maritime tropical cumulus congestus having a similar vertical extent. Modeled drop size distributions are similar to observations at various temperatures between 9.5 and –12°C, including the lack of mm-sized drops. Cloud dilution leads to low-to-moderate liquid water contents (~0.5–1.5 g m –3 ) in most of the cloud core, several times smaller than adiabatic values. Dilution is enhanced in the inflowing branch of the toroidal circulations associated with individual cloud thermals, which are favored regions for secondary droplet activation (activation above cloud base). Secondary activation in general contributes substantially to the droplet population. Turning it off leads to a sharp decrease in droplet concentration and increase in mean size aloft, but does little to increase rain drop production. Warm rain generation (or lack thereof) in this case is therefore determined more by the sub-cloud aerosol and the cloud base droplet size distribution (DSD) than DSD evolution aloft from secondary droplet activation. Decreasing the aerosol concentration by a factor of 10 greatly increases production of large drops via collision-coalescence. Thus, despite its high base (low temperatures) and substantial dilution, the simulated cloud is thermodynamically and dynamically capable of rapidly producing copious mm-sized drops from collision-coalescence under pristine aerosol conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Europa modifies Jupiter’s 1 plasma sheet

Jupiter's plasma sheet has been understood to be primarily composed of Io-genic sulfur and oxygen, along with protons at lower mass density. These ions move radially away from Jupiter, filling its magnetosphere. The material in the plasma sheet interacts with Europa, which is also a source of magnetospheric pickup ions, primarily hydrogen and oxygen. Juno's thermal plasma instrument JADE, the Jovian Auroral Distributions Experiment, has provided comprehensive in situ observations of the composition of Jupiter's plasma sheet ions with its Time-of-Flight mass-spectrometry capabilities. Here, we present observations of the magnetospheric composition in the Europa-Ganymede region of Jupiter's magnetosphere. We find material from Europa is intermittently present at comparable densities to Io-genic plasma. The intermittency of Europa-genic signatures suggests Europa's neutral oxygen toroidal cloud is more localized to Europa's vicinity than its hydrogen cloud. These observations reveal a more complex and compositionally diverse magnetosphere than previously thought.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Revealing the three-dimensional arrangement of polar topology in nanoparticles

In the early 2000s, low dimensional ferroelectric systems were predicted to have topologically nontrivial polar structures, such as vortices or skyrmions, depending on mechanical or electrical boundary conditions. A few variants of these structures have been experimentally observed in thin film model systems, where they are engineered by balancing electrostatic charge and elastic distortion energies. However, the measurement and classification of topological textures for general ferroelectric nanostructures have remained elusive, as it requires mapping the local polarization at the atomic scale in three dimensions. Here we unveil topological polar structures in ferroelectric BaTiO3 nanoparticles via atomic electron tomography, which enables us to reconstruct the full three-dimensional arrangement of cation atoms at an individual atom level. Our three-dimensional polarization maps reveal clear topological orderings, along with evidence of size-dependent topological transitions from a single vortex structure to multiple vortices, consistent with theoretical predictions. The discovery of the predicted topological polar ordering in nanoscale ferroelectrics, independent of epitaxial strain, widens the research perspective and offers potential for practical applications utilizing contact-free switchable toroidal moments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Demonstration of reduced neoclassical energy transport in Wendelstein 7-X

Research on magnetic confinement of high-temperature plasmas has the ultimate goal of harnessing nuclear fusion for the production of electricity. Although the tokamak is the leading toroidal magnetic-confinement concept, it is not without shortcomings and the fusion community has therefore also pursued alternative concepts such as the stellarator. Unlike axisymmetric tokamaks, stellarators possess a three-dimensional (3D) magnetic field geometry. The availability of this additional dimension opens up an extensive configuration space for computational optimization of both the field geometry itself and the current-carrying coils that produce it. Such an optimization was undertaken in designing Wendelstein 7-X (W7-X), a large helical-axis advanced stellarator (HELIAS), which began operation in 2015 at Greifswald, Germany. A major drawback of 3D magnetic field geometry, however, is that it introduces a strong temperature dependence into the stellarator’s non-turbulent ‘neoclassical’ energy transport. Indeed, such energy losses will become prohibitive in high-temperature reactor plasmas unless a strong reduction of the geometrical factor associated with this transport can be achieved; such a reduction was therefore a principal goal of the design of W7-X. In spite of the modest heating power currently available, W7-X has already been able to achieve high-temperature plasma conditions during its 2017 and 2018 experimental campaigns, producing record values of the fusion triple product for such stellarator plasmas. The triple product of plasma density, ion temperature and energy confinement time is used in fusion research as a figure of merit, as it must attain a certain threshold value before net-energy-producing operation of a reactor becomes possible. Here we demonstrate that such record values provide evidence for reduced neoclassical energy transport in W7-X, as the plasma profiles that produced these results could not have been obtained in stellarators lacking a comparably high level of neoclassical optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Ferroelectric incommensurate spin crystals

Ferroics, especially ferromagnets, can form complextopological spin structures such as vortices and skyrmions when subjected to particular electrical and mechanical boundary conditions. Simple vortex-like, electric-dipole-based topological structures have been observed in dedicated ferroelectric systems, especially ferroelectric-insulator superlattices such as PbTiO 3 /SrTiO 3 , which was later shown to be a model system owing to its high depolarizing field. To date, the electric dipole equivalent of ordered magnetic spin lattices driven by the Dzyaloshinskii-Moriya interaction (DMi) has not been experimentally observed. Here we examine a domain structure in a single PbTiO 3 epitaxial layer sandwiched between SrRuO 3 electrodes. In this work, we observe periodic clockwise and anticlockwise ferroelectric vortices that are modulated by a second ordering along their toroidal core. The resulting topology, supported by calculations, is a labyrinth-like pattern with two orthogonal periodic modulations that form an incommensurate polar crystal that provides a ferroelectric analogue to the recently discovered incommensurate spin crystals in ferromagnetic materials. These findings further blur the border between emergent ferromagnetic and ferroelectrictopologies, clearing the way for experimental realization of further electric counterparts of magnetic DMi-driven phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗