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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Experiment-modeling studies comparing energy dissipation in the DIII-D SAS and SAS-VW divertors

Recent DIII-D experiments on Small Angle Slot (SAS) divertors have confirmed that a combination of divertor closure and target shaping can enhance cooling across the divertor target and increase energy dissipation, but with significant dependence on B T (toroidal magnetic field) direction. In these novel divertors, the roles of closure, target shaping, drifts, and scale lengths are all interconnected in optimizing dissipation, with the separatrix electron density n eSEP being the key parameter associated with the level of dissipation/detachment. After modifying the original flat-targeted graphite SAS to include a V shape with a tungsten coating on the outer side of the divertor (SAS-VW), matched series of discharges were run to compare to detailed SOLPS-ITER modeling. Experimentally, when run as designed with the outer strike point at the slot vertex, SAS-VW requires nearly identical n eSEP for detachment as the original SAS, with little difference in dissipation for the new geometry. This is in contrast to (1) earlier modeling predictions that a small change of the SAS geometry to a V shape should enhance dissipation at the same n eSEP for magnetic configurations having better H-mode access (ion B × ∇B drift directed into the divertor), and (2) despite the achievement of significantly higher (2-7x) neutral pressures and compression in the SAS-VW slot. Comparisons of experimental density scans to the most recent SOLPS-ITER modeling with ExB drifts show reasonable agreement for dissipation/detachment onset when using separatrix density as the independent parameter. In order to help understand the discrepancy in modeled vs actual performance for the new configuration, additional measurements varying gas injection location and impurity injection were undertaken. In-slot D 2 gas fueling is more effective (5–22 %) in promoting detachment, in accord with modeling. In-slot impurity injection (N 2 or Ne) can yield 30 % lower core Z eff and 15 % less confinement degradation after detachment compared to main chamber puffing, as well as relatively lower tungsten leakage from the divertor. Modeling can also reproduce the improved detachment seen as the strike point moves inboard of the slot vertex. While we can explain the effects of the most important parameters causing energy dissipation in these slot divertors, it remains that many aspects of their behavior cannot be accurately modeled using state-of-art codes such as SOLPS-ITER. This is of concern for future model-driven designs utilizing similar V-shaped geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Searching for solitons in heavy-ion reactions near the fermi energy

Solitons are special shape-conserving hydrodynamical solutions that appear in many areas of physics. Here, we explore the existence of such solutions in microscopic descriptions of the heavy ion reaction 12 C + 28 Si → 12 C*+ 7α in the range E/A = 10 – 65 MeV/u. After recognizing the centrality of the collision and time-reversibility as fundamental requirements for the presence of solitonic 12 C, we utilize the Hybrid α-Cluster model for our analysis with a novel methodology. Our results suggest soliton production for E/A = 25 – 45 MeV/u in the forward direction, with a total cross section being at least in the order of a few μb. This, apart from being encouraging for possible experimental studies, maybe connected to possible toroid states of 28 Si at low angular momenta.

Heavy ion reactions↗
DOE OSTI2025

Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory

Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.

Drops and Bubbles: Bubble dynamics↗
DOE OSTI2020

Fast-ion physics in SPARC

Potential loss of energetic ions including alphas and radio-frequency tail ions due to classical orbit effects and magnetohydrodynamic instabilities (MHD) are central physics issues in the design and experimental physics programme of the SPARC tokamak. The expected loss of fusion alpha power due to ripple-induced transport is computed for the SPARC tokamak design by the ASCOT and SPIRAL orbit-simulation codes, to assess the expected surface heating of plasma-facing components. We find good agreement between the ASCOT and SPIRAL simulation results not only in integrated quantities (fraction of alpha power loss) but also in the spatial, temporal and pitch-angle dependence of the losses. If the toroidal field (TF) coils are well-aligned, the SPARC edge ripple is small (0.15–0.30 %), the computed ripple-induced alpha power loss is small ( ~0.25% ) and the corresponding peak surface power density is acceptable ( 244 kW m -2 ). However, the ripple and ripple-induced losses increase strongly if the TF coils are assumed to suffer increasing magnitudes of misalignment. Surface heat loads may become problematic if the TF coil misalignment approaches the centimetre level. Ripple-induced losses of the energetic ion tail driven by ion cyclotron range of frequency (ICRF) heating are not expected to generate significant wall or limiter heating in the nominal SPARC plasma scenario. Because the expected classical fast-ion losses are small, SPARC will be able to observe and study fast-ion redistribution due to MHD including sawteeth and Alfvén eigenmodes (AEs). SPARC's parameter space for AE physics even at moderate Q is shown to reasonably overlap that of the demonstration power plant ARC (Sorbom et al., Fusion Engng Des., vol. 100, 2015, p. 378), and thus measurements of AE mode amplitude, spectrum and associated fast-ion transport in SPARC would provide relevant guidance about AE behaviour expected in ARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Drift kinetic theory of alpha transport by tokamak perturbations

We report that upcoming tokamak experiments fuelled with deuterium and tritium are expected to have large alpha particle populations. Such experiments motivate new attention to the theory of alpha particle confinement and transport. A key topic is the interaction of alpha particles with perturbations to the tokamak fields, including those from ripple and magnetohydrodynamic modes like Alfvén eigenmodes. These perturbations can transport alphas, leading to changed localization of alpha heating, loss of alpha power and damage to device walls. Alpha interaction with these perturbations is often studied with single-particle theory. In contrast, we derive a drift kinetic theory to calculate the alpha heat flux resulting from arbitrary perturbation frequency and periodicity (provided these can be studied drift kinetically). Novel features of the theory include the retention of a large effective collision frequency resulting from the resonant alpha collisional boundary layer, correlated interactions over many poloidal transits and finite orbit effects. Heat fluxes are considered for the example cases of ripple and the toroidal Alfvén eigenmode (TAE). The ripple heat flux is small. The TAE heat flux is significant and scales with the square of the perturbation amplitude, allowing the derivation of constraints on mode amplitude for avoidance of significant alpha depletion. A simple saturation condition suggests that TAEs in one upcoming experiment will not cause significant alpha transport via the mechanisms in this theory. However, saturation above the level suggested by the simple condition, but within numerical and experimental experience, which could be accompanied by the onset of stochasticity, could cause significant transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Vacuum magnetic fields with exact quasisymmetry near a flux surface. Part 1. Solutions near an axisymmetric surface

While several results have pointed to the existence of exactly quasisymmetric fields on a surface (Garren & Boozer, Phys. Fluids B, vol. 3, 1991, pp. 2805–2821; 2822–2834; Plunk & Helander, J. Plasma Phys. , vol. 84, 2018, 905840205), we have obtained the first such solutions using a vacuum surface expansion formalism. We obtain a single nonlinear parabolic partial differential equation 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 quasihelical 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. The local solutions thus obtained can be continued globally only for special initial surfaces.

Physics↗
DOE OSTI2021

Simple, general, realistic, robust, analytic tokamak equilibria. Part 1. Limiter and divertor tokamaks

Tokamak equilibria have been derived that are analytic solutions to the Grad–Shafranov equation. Here in this paper, Part 1, describes a wide range of such equilibria including smooth limiter surfaces, double- and single-null divertor surfaces, arbitrary aspect ratio, elongation, triangularity and beta. Part 2 generalizes the analysis to include edge pedestals and toroidal flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The on-axis magnetic well and Mercier's criterion for arbitrary stellarator geometries

A simplified analytical form of the on-axis magnetic well and Mercier's criterion for interchange instabilities for arbitrary three-dimensional magnetic field geometries is derived. For this purpose, a near-axis expansion based on a direct coordinate approach is used by expressing the toroidal magnetic flux in terms of powers of the radial distance to the magnetic axis. For the first time, the magnetic well and Mercier's criterion are then written as a one-dimensional integral with respect to the axis arclength. When compared with the original work of Mercier, the derivation here is presented using modern notation and in a more streamlined manner that highlights essential steps. Finally, these expressions are verified numerically using several quasisymmetric and non-quasisymmetric stellarator configurations including Wendelstein 7-X.

Physics↗
DOE OSTI2021

Bispectral analysis of broadband turbulence and geodesic acoustic modes in the T-10 tokamak

Local fluctuations of electrostatic potential, poloidal electric field, magnetic potential and electron density are simultaneously measured in the T-10 tokamak by a heavy ion beam probe (HIBP) having a five-slit energy analyser, which allows an estimate of the turbulent particle flux and $\boldsymbol {E}\times \boldsymbol {B}$ rotation velocity in the off-minor-axis gradient zone of the toroidal plasma column. The high spatial and temporal resolution of the modern multichannel HIBP makes it an effective tool to study plasma oscillations. Motivated by previous work that has documented time-resolved interactions between measured plasma parameters using correlation analysis (coherence of $E_{\textrm {pol}}$ and density $n_e$ , and cross-phase), a new result from bicorrelation analysis (bicoherence of magnetic potential $A_\zeta$ and density $n_e$ , and biphase) is reported for documenting the evidence of wave–wave coupling and energy transfer associated with the interaction between geodesic acoustic modes (GAM) and broadband, quasi-coherent modes.

Physics↗
DOE OSTI2021

Phase contrast imaging measurements and numerical simulations of turbulent density fluctuations in gas-fuelled ECRH discharges in Wendelstein 7-X

The fundamental nature of turbulent density fluctuations in standard Wendelstein 7-X (W7-X) stellarator discharges is investigated experimentally via phase contrast imaging (PCI) in combination with gyrokinetic simulations with the code GENE. We find that density fluctuations are ion-temperature-gradient-driven and radially localised in the outer half of the plasma. It is shown that the line-integrated PCI measurements cover the right range of wavenumbers and a favourable toroidal and poloidal location to capture some of the strongest density fluctuations in W7-X. Due to the radial localisation of fluctuations, measured wavenumber–frequency spectra exhibit a dominant phase velocity, which can be related to the E × B rotation velocity at the radial position of a well in the neoclassical radial electric field. The match is robust against variations of heating power and line-integrated density, which is partly due to the localisation of fluctuations and partly due to effects of the radial gradient in the E × B velocity profile on the wavenumber–frequency spectrum. The latter effect is studied with a newly built synthetic PCI diagnostic and global gyrokinetic simulations with GENE-3D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Phase transition to blob-hole coherent structure in the Hasegawa–Mima model for plasmas

An equilibrium statistical mechanics theory for the Hasegawa–Mima equations of toroidal plasmas, with canonical constraint on energy and microcanonical constraint on potential enstrophy, is solved exactly as a spherical model. The use of a canonical energy constraint instead of a fixed-energy microcanonical approach is justified by the preference for viewing real plasmas as an open system. A significant consequence of the results obtained from the partition function, free energy and critical temperature, is the condensation into a ground state exhibiting a blob-hole-like structure observed in real plasmas.

Physics↗
DOE OSTI2022

Description of global EGAM in the maximum of local frequency during current ramp-up discharges in DIII-D

Energetic-particle-induced geodesic acoustic modes, EGAMs (Fu, Phys. Rev. Let., vol. 101, 2008, pp. 185002), driven by neutral beam injection (NBI), have been observed in many DIII-D tokamak experiments (Nazikian et al., Phys. Rev. Lett., vol. 101, 2008, pp. 185001). This mechanism has been theoretically investigated in (Qiu et al., Plasma Phys. Control. Fusion, vol. 52, 2010, pp. 095003), using a sharp energetic particle distribution function, and in (Qu et al., Plasma Phys. Control. Fusion, vol. 59, 2017, pp. 055018), where the dispersion relation and eigenmode behaviour were obtained for the situation of early beam scenario, that is, for times smaller than the beam slowing down time. In this work, we extend these studies determining the eigenmode for beyond the slowing down time, in a scenario with reverse safety factor q profile, where a small concentration of energetic ions can produce an off-axis maximum in the GAM dispersion relation. Here, the characteristics of EGAM are analytically studied with the drift kinetic equation together with the MHD code NOVA. The toroidal energetic ion transit frequency, coupled with the GAM frequency, produces the maximum in the dispersion relation where the eigenmode can be found. The quantitative correspondence of experimental results with the predictions of the proposed model is analysed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Calculation of collisionless pitch-angle scattering of runaway electrons with synchrotron radiation via high-order guiding-centre equation

Recently, the collisionless pitch-angle scattering for relativistic runaway electrons (REs) in toroidal geometries such as tokamaks was discovered through a full orbit simulation approach (Liu et al., Nucl. Fusion, vol. 56, 2016, p. 064002), and it was then theoretically investigated that a new expression for the magnetic moment, including the second-order corrections, could essentially reproduce the so-called collisionless pitch-angle scattering process (Liu et al., Nucl. Fusion, vol. 58, 2018, p. 106018). In this paper, with synchrotron radiation, extensive numerical verification of the validity of the high-order guiding-centre theory is given for simulations involving REs by incorporating such an expression for the magnetic moment into our particle tracing code. A high-order guiding-centre simulation approach with synchrotron radiation (HGSA) is applied. Synchrotron radiation plays an essential role in the life cycle of REs. The energy of REs first increases and then becomes saturated until the electric field acceleration is balanced by the radiation dissipation. Unfortunately, the process cannot be simulated accurately with the standard guiding-centre model, i.e. the first-order guiding-centre model. Remarkably, it is found that the HGSA can effectively produce the fundamental process of REs. Since the time scale of the energy saturation of REs is close to seconds, the computational cost becomes significant. In order to save costs, it is necessary to estimate the time of energy saturation. An analytical estimate is derived for the time it takes for synchrotron drag to balance an accelerating electric field and the provided formula has been numerically verified. Finally, test calculations reveal that HGSA is favourable for exploiting the dynamics of REs in tokamak plasmas.

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