On the symmetrical deformation of an orthotropic toroidal shell
Solution for symmetrical deformation problem of orthotropic toroidal shell
TOROIDAL SHELL↗
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.
Stressed-strained state of a toroidal shell at bending, applicable to calculating bellows
Stress-strain state of toroidal shell at bending, applicable to calculating bellows
BELLOWS↗
Calculation of a ribbed toroidal shell under a symmetrical load
Model of structural anisotropy for determining basic stresses in circumference and external fibers of symmetrically loaded ribbed toroidal shell
ANISOTROPY↗
Equilibrium of electron clouds in toroidal magnetic fields.
Electron cloud equilibrium in toroidal magnetic and electric fields due to space and image charges
PLASMA CONFINEMENT↗
ITER Toroidal Interferometer and Polarimeter (TIP) beam refraction in 3D density profiles
Calculations of the expected refraction of the five ITER Toroidal Interferometer and Polarimeter (TIP) chords are presented for a range of conditions including high-density axisymmetric, ELMing H-mode cases and shattered pellet disruption mitigation cases. The calculations are carried out with a newly developed ray tracing code capable of following TIP’s 10.59 μm laser beams through arbitrary 3D electron density profiles. Using JOREK simulations of ELM density perturbations in a 15MA ITER baseline plasma, it is shown that refraction from ELMs is expected to be negligible. It is found, however, that TIP interferometers will be able to clearly resolve the line-integrated density perturbation from ELMs and contribute to the ITER measurement requirement “14. H-mode, ELMs and L-H mode transition indicator”. Calculations of the expected refraction in a NIMROD simulated shattered pellet disruption mitigation scenario with peak local electron densities of n e = 3.7 x 10 21 m -3 (max line-averaged densities of n avg = 1.3 x 10 21 m -3 ) also show tolerable refraction and it is likely that most, if not all chords, would avoid signal loss. In conclusion, for both axisymmetric and structured plasmas with line-averaged densities in the mid to upper 10 21 m -3 range, values likely present only during disruption mitigation, refraction becomes significant and could limit the ability of TIP to make reliable density measurements.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Runaway electron transport in stochastic toroidal magnetic fields
Here, we study the transport and confinement properties of runaway electrons (RE) in the presence of magnetic fields with perturbations producing different levels of stochasticity. We use Kinetic Orbit Runaway Electron Code (KORC) [Carbajal et al. , Phys. Plasmas 24, 042512 (2017) and del-Castillo-Negrete et al. , Phys. Plasmas 25, 056104 (2018)] for simulating the full-orbit (FO) and guiding-center (GC) dynamics of RE in perturbed magnetic fields that exhibit magnetic islands. We extend previous works on this problem [Wingen et al. , Nucl. Fusion 46, 941 (2006); Izzo et al. , Nucl. Fusion 51, 063032 (2011); Papp et al. , Nucl. Fusion 51, 043004 (2011); V. Izzo and P. Parks, Phys. Plasmas 24, 060705 (2017); and Sommariva et al. , Nucl. Fusion 58, 016043 (2018)] by studying in detail full-orbit effects on the RE dynamics. We quantify FO effects on RE transport by performing one-to-one comparisons between FO and GC simulations. It is found that, for the magnetic field configurations considered, GC simulations predict twice the RE losses of FO simulations for 1 MeV and four times the RE losses of FO simulations for 25 MeV. Similarly, we show how different GC and FO dynamics of RE moving around magnetic islands can be, especially in the scenario where the RE Larmor radius is on the order of the size of the magnetic island. We also study the role of rotation of the magnetic islands on RE confinement, and we find that low-frequency toroidal rotation has no observable effect on RE transport in the cases considered. These results shed some light into the potential of avoidance or mitigation mechanisms based on magnetic perturbations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Toroidal modeling of runaway electron loss due to 3-D fields in DIII-D and COMPASS
The 3-D field induced relativistic runaway electron (RE) loss has been simulated for DIII-D, COMPASS and ITER plasmas, utilizing the MARS-F code incorporated with the recently developed and updated RE orbit module (REORBIT). Modeling shows effectively 100% loss of a post-disruption, high-current runaway beam in DIII-D, due to the 1 kG level of magnetic field perturbation produced by a fast growing n = 1 resistive kink instability. The RE loss is shown to be independent of the particle energy or the initial location of particles in the configuration space. Applied resonant magnetic perturbation (RMP) fields from in-vessel coils are not effective for RE beam mitigation in DIII-D, but do produce finite (>10%) RE loss in COMPASS, consistent with experimental observations in above two devices. The major reasons for this difference in RE control by RMP between these two devices are (i) the coil proximity to the RE beam and (ii) the effective coil current scaling versus the machine size and the toroidal magnetic field. About 10% RE loss fraction is also predicted for an ITER 15 MA scenario with pre-disruption plasma, highlighting the role of the plasma response. Up to 30% loss is computed, however, by artificially scaling the equilibrium pressure to zero. This is due to the more resistive plasma response and stronger resulting field line stochasticity. Distributions of the lost REs to the limiting surface show poloidally peaked profile near the high-field-side in both DIII-D and COMPASS, covering about 100o poloidal angle. Higher perturbation field level and/or higher particle energy also result in REs being lost to the low-field-side of the limiting surface of these two devices, increasing the effective wetted area. Finally, for the modeled ITER plasmas, the applied RMP field, with optimal poloidal spectrum that maximizes the edge localized mode control, leads to REs being lost to the lower divertor region of the limiting surface within a narrow poloidal band.
3-D field↗
Hybrid-gyrokinetic simulations of low- n toroidal Alfvén eigenmodes using gKPSP
Here, we report a benchmark study of toroidal Alfvén eigenmode (TAE) simulation using the hybrid-gyrokinetic code GyroKinetic Plasma Simulation Program (gKPSP). A simulation capability for energetic particles based on the gyrokinetic δf method has been newly implemented in the gKPSP code. Benchmark simulations have been performed in both circular and realistic tokamak geometries. Good agreement has been found with previously reported results, demonstrating the new capability of the gKPSP code. We have investigated the effects of the distribution function on TAE stability by examining both isotropic and anisotropic slowing-down distributions of energetic particles. The slowing-down distribution produces a higher linear growth rate than a Maxwellian distribution, while a growth rate scan with the anisotropy shows an opposite trend. This can be attributed to competition between Landau damping and the linear drive, which are correlated with the fraction of resonant passing particles and their distribution in phase space.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Verification of a fully kinetic ion model for electromagnetic simulations of high-frequency waves in toroidal geometry
For the research of high-frequency electromagnetic waves in tokamaks, an electromagnetic simulation model, in which the ion dynamics is described by a six-dimensional Vlasov equation and the electron dynamics is described by a drift kinetic equation, is formulated and implemented in the global gyrokinetic toroidal code (GTC). Analytic dispersion relations are derived in reduced systems and compared with various theories to verify the model. Linear simulations of a generalized ion Bernstein wave and ion cyclotron emission are verified by comparing the GTC simulation results with analytic dispersion relation theory and magnetoacoustic cyclotron instability theory, respectively, in cylindrical geometry.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
ATEQ: Adaptive toroidal equilibrium code
A radially adaptive numerical scheme is developed to solve the Grad–Shafranov equation for axisymmetric magnetohydrodynamic equilibrium. Here, a decomposition with independent solutions is employed in the radial direction, and Fourier decomposition is used in the poloidal direction. The independent solutions are then obtained using an adaptive shooting scheme together with the multi-region matching technique in the radial direction. Accordingly, the adaptive toroidal equilibrium (ATEQ) code is constructed for axisymmetric equilibrium studies. The adaptive numerical scheme in the radial direction improves considerably the accuracy of the equilibrium solution. The decomposition with independent solutions effectively reduces the matrix size in solving the magnetohydrodynamic equilibrium problem. The reduction of the matrix size is about an order of magnitude as compared with the conventional radially grid-based numerical schemes. Also, in this ATEQ numerical scheme, no matter how accuracy in the radial direction is imposed, the size of matrices basically does not change. The small matrix size scheme gives ATEQ more flexibility to address the requirement of the number of Fourier components in the poloidal direction in tough equilibrium problems. These two unique features, the adaptive shooting and small matrix size, make ATEQ useful to improve tokamak equilibrium solutions.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
All-in-one probe for exploring self-organized two-fluid equilibria in toroidal plasmas
This paper presents the development of an all-in-one probe to simultaneously measure all components of the generalized Ohm’s law in reversed-field pinch plasmas and tokamaks. The polyhedral configuration of the Mach probe is achieved through the specific arrangement, angle, and depth of the collimator channel apertures drilled into the surface of a hollow boron nitride cylinder encasing it. This probe includes a central Mach probe to assess the ion velocity field in three dimensions. Initial tests at the RELAX and Madison Symmetric Torus machines have confirmed the probe’s effectiveness, revealing an octahedron form similar to a tetrahedron. The probe seems to function correctly and is expected to facilitate the empirical validation of two-fluid equilibria at the periphery of toroidal plasmas.
Instruments & Instrumentation↗
Toroidal Alfvén mode instability driven by plasma current in low-density Ohmic plasmas of the spherical tori
Due to intrinsically low magnetic fields, at low density the drift speed of the current-carrying electrons in spherical tokamaks can exceed fraction of the Alfvén speed sufficient for the excitation of the Alfvén gap modes. A particular case of the toroidal mode, observed during minor disruptions in Ohmic shots on SUNIST [Liu et al., Phys. Plasmas 23, 120706 (2016)], is considered in the present communication. Due to the negligible effect of the electron pressure gradient, the growth rate scales linearly with the drift speed, with slope inversely proportional to electron thermal velocity. In the absence of continuum damping, the threshold value of the drift speed for TAE excitation is independent of electron temperature.
Marchenko, V. S. (ORCID:0000000261480533)↗
An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry
A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ ). Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Efficient analysis of magnetic field line behavior in toroidal plasmas
The confinement of plasmas in tokamaks and stellarators depends on magnetic field lines lying in nested toroidal surfaces. The transition near the plasma edge away from the lines lying in magnetic surfaces defines properties of divertors. The transition in time defines properties of disruptions. Divertor design and disruption analyses require a detailed understanding of these transitions. The use of a Fourier transform coupled with a Gaussian window function allows far more information to be extracted about these transitions using far shorter field line integrations than can be obtained using traditional methods based on Poincaré plots. The physics of divertors and disruptions is reviewed to clarify why the type of information that can be gained from more efficient methods of analysis is of central importance to the fusion program based on magnetic confinement.
Boozer, Allen H. [Columbia Univ., New York, NY (Un↗
Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor
Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx , and inter-null orientations, θ, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm , d xx , and θ. For small values of θ, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.
Computer simulation↗
Simulations of plasma current induced by toroidal field ramping down on tokamaks
In tokamak experiments, the BT ramping down can drive off-axis and parallel inductive current. This approach leads to a decrease in q95 and an increase in normalized beta, βN. The off-axis inductive current also broadens the current profile. Typically, the 1.5D transport code is used to simulate the time evolution of the plasma current profile, which is based on the flux-surface-averaged Faraday’s law. The ONETWO code is one of those transport codes. However, this code cannot simulate the situation of the evolving toroidal field BT. This article proposes a new time-dependent model to take into account the BT ramping down situation. A modified formula of flux-surface-averaged Faraday’s law was derived to consider the effect of the BT ramping rate on the current evolution, and it was implemented in the ONETWO code. Then, the modified ONETWO code was used to simulate the BT ramping down experiment on DIII-D. The simulation result of the plasma current evolution with BT ramping down shows a broader current profile with smaller ohmic current induced by the poloidal field, compared to that without BT ramping down.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Magnetic toroidicity
Abstract Directional non-reciprocity refers to the phenomenon where the motion in one direction differs from the motion in the opposite direction. This behavior is observed across various systems, such as one-way traffic and materials displaying electronic/optical directional dichroism, characterized by the symmetry of velocity vectors. Magnetic toroidal moments (MTMs), which typically arise from rotational spin arrangements, also possess the symmetry of velocity vectors, making them inherently directionally non-reciprocal. In this paper, we examine magnetic point groups (MPGs) that exhibit MTMs, subsequently leading to off-diagonal linear magnetoelectricity. Our focus is on the induction of MTMs through electric fields, magnetic fields, or shear stress, while enumerating the relevant MPGs. The findings of our study will serve as valuable guidance for future investigations on directional non-reciprocity, MTMs, and off-diagonal linear magnetoelectric effects.
Physics↗
Fast neural Poincaré maps for toroidal magnetic fields
Poincaré maps for toroidal magnetic fields are routinely employed to study gross confinement properties in devices built to contain hot plasmas. In most practical applications, evaluating a Poincaré map requires numerical integration of a magnetic field line, a process that can be slow and that cannot be easily accelerated using parallel computations. We propose a novel neural network architecture, the HénonNet, and show that it is capable of accurately learning realistic Poincaré maps from observations of a conventional field-line-following algorithm. After training, such learned Poincaré maps evaluate much faster than the field-line integration method. Moreover, the HénonNet architecture exactly reproduces the primary physics constraint imposed on field-line Poincaré maps: flux preservation. Furthermore, this structure-preserving property is the consequence of each layer in a HénonNet being a symplectic map. We demonstrate empirically that a HénonNet can learn to mock the confinement properties of a large magnetic island by using coiled hyperbolic invariant manifolds to produce a sticky chaotic region at the desired island location. This suggests a novel approach to designing magnetic fields with good confinement properties that may be more flexible than ensuring confinement using KAM tori.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗