Search NASA⌕ Search

SEARCH · Search NASA

Results for “toroidal”

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 217 records · Page 12

Improved Method of Design for Folding Inflatable Shells

An improved method of designing complexly shaped inflatable shells to be assembled from gores was conceived for original application to the inflatable outer shell of a developmental habitable spacecraft module having a cylindrical mid-length section with toroidal end caps. The method is also applicable to inflatable shells of various shapes for terrestrial use. The method addresses problems associated with the assembly, folding, transport, and deployment of inflatable shells that may comprise multiple layers and have complex shapes that can include such doubly curved surfaces as toroids and spheres. One particularly difficult problem is that of mathematically defining fold lines on a gore pattern in a double- curvature region. Moreover, because the fold lines in a double-curvature region tend to be curved, there is a practical problem of how to implement the folds. Another problem is that of modifying the basic gore shapes and sizes for the various layers so that when they are folded as part of the integral structure, they do not mechanically interfere with each other at the fold lines. Heretofore, it has been a common practice to design an inflatable shell to be assembled in the deployed configuration, without regard for the need to fold it into compact form. Typically, the result has been that folding has been a difficult, time-consuming process resulting in a An improved method of designing complexly shaped inflatable shells to be assembled from gores was conceived for original application to the inflatable outer shell of a developmental habitable spacecraft module having a cylindrical mid-length section with toroidal end caps. The method is also applicable to inflatable shells of various shapes for terrestrial use. The method addresses problems associated with the assembly, folding, transport, and deployment of inflatable shells that may comprise multiple layers and have complex shapes that can include such doubly curved surfaces as toroids and spheres. One particularly difficult problem is that of mathematically defining fold lines on a gore pattern in a double- curvature region. Moreover, because the fold lines in a double-curvature region tend to be curved, there is a practical problem of how to implement the folds. Another problem is that of modifying the basic gore shapes and sizes for the various layers so that when they are folded as part of the integral structure, they do not mechanically interfere with each other at the fold lines. Heretofore, it has been a common practice to design an inflatable shell to be assembled in the deployed configuration, without regard for the need to fold it into compact form. Typically, the result has been that folding has been a difficult, time-consuming process resulting in a

Johnson, Christopher J.↗

Determining the Mode, Frequency, and Azimuthal Wave Number of ULF Waves During a HSS and Moderate Geomagnetic Storm

Ultralow frequency (ULF) waves play a fundamental role in the dynamics of the inner magnetosphere and outer radiation belt during geomagnetic storms. Broadband ULF wave power can transport energetic electrons via radial diffusion, and discrete ULF wave power can energize electrons through a resonant interaction. Using observations from the Magnetospheric Multiscale mission, we characterize the evolution of ULF waves during a highspeed solar wind stream (HSS) and moderate geomagnetic storm while there is an enhancement of the outer radiation belt. The Automated Flare Inference of Oscillations code is used to distinguish discrete ULF wave power from broadband wave power during the HSS. During periods of discrete wave power and utilizing the close separation of the Magnetospheric Multiscale spacecraft, we estimate the toroidal mode ULF azimuthal wave number throughout the geomagnetic storm. We concentrate on the toroidal mode as the HSS compresses the dayside magnetosphere resulting in an asymmetric magnetic field topology where toroidal mode waves can interact with energetic electrons. Analysis of the mode structure and wave numbers demonstrates that the generation of the observed ULF waves is a combination of externally driven waves, via the KelvinHelmholtz instability, and internally driven waves, via unstable ion distributions. Further analysis of the periods and toroidal azimuthal wave numbers suggests that these waves can couple with the core electron radiation belt population via the drift resonance during the storm. The azimuthal wave number and structure of ULF wave power (broadband or discrete) have important implications for the inner magnetospheric and radiation belt dynamics.

Murphy, Kyle R.↗

Application of Synoptic Magnetograms for Prediction of Solar Activity Using Ensemble Kalman Filter

Solar activity predictions using the data assimilation approach have demonstrated great potential to build reliable long-term forecasts of solar activity. In particular, it has been shown that the Ensemble Kalman Filter (EnKF) method applied to a non-linear dynamo model is capable of predicting solar activity up to one sunspot cycle ahead in time, as well as estimating the properties of the next cycle a few years before it begins. These developments assume an empirical relationship between the mean toroidal magnetic field flux and the sunspot number. Estimated from the sunspot number series, variations of the toroidal field have been used to assimilate the data into the Parker-Kleeorin-Ruzmakin (PKR) dynamo model by applying the EnKF method. The dynamo model describes the evolution of the toroidal and poloidal components of the magnetic field and the magnetic helicity. Full-disk magnetograms provide more accurate and complete input data by constraining both the toroidal and poloidal global field components, but these data are available only for the last four solar cycles. In this presentation, using the available magnetogram data, we discuss development of the methodology and forecast quality criteria (including forecast uncertainties and sources of errors). We demonstrate the influence of limited time series observations on the accuracy of solar activity predictions. We present EnKF predictions of the upcoming Solar Cycle 25 based on both the sunspot number series and observed magnetic fields and discuss the uncertainties and potential of the data assimilation approach.

Kitiashvili, Irina N.↗

Charge-exchange losses of beam ions in MAST Upgrade: experiments and modelling

Dedicated experiments were performed on MAST Upgrade to study beam-ion losses caused by charge exchange (CX) with edge neutrals. The fuelling was switched from the high-field side to the low-field side mid-discharge. Direct measurements suggest a strong increase in the neutral density around the plasma and a decrease in the beam-ion density, which is qualitatively explained by CX losses. Measurements by a resistive bolometer have suggested particle bombardment during neutral beam injection, providing a unique opportunity to separate CX from other loss mechanisms. To verify and quantify CX losses, the orbit-following code ASCOT, which accounts both for CX neutralization and reionization, was used to simulate beam-particle power loads on the bolometer. Simulations reproduce measured bolometer power loads during high-field-side fuelling, verifying CX losses of approximately 10% of the off-axis beam power. Toroidally symmetric simulations overestimate power loads on the bolometer during low-field-side fuelling, which is explained by toroidal asymmetry in the neutral density distribution, as is demonstrated by toroidally asymmetric simulations. Results suggest significantly higher CX losses during low-field-side fuelling, up to about 50% of off-axis beam power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of Injected Current Streams on MHD Equilibrium Reconstruction of Local Helicity Injection Plasmas in a Spherical Tokamak

Open field line currents are intrinsic to DC helicity injection plasma startup and pose a challenge for inferring the plasma equilibrium with standard reconstruction analysis. Local helicity injection (LHI) is a type of DC helicity injection which uses small, modular current sources to drive force-free current along helical field lines to produce tokamak plasmas. MHD modeling and magnetic measurements during LHI indicate the injected current streams remain coherent as helical structures on the outboard edge of a core toroidal plasma that is tokamak-like in a toroidally averaged sense. To extract core plasma equilibrium properties, external magnetic diagnostics corrected for contributions from the injected current streams are fitted by a standard Grad-Shafranov equilibrium code. An iterative approach for estimating and subtracting the stream contributions from the diagnostic signals is described and applied to a model equilibrium database to reduce systematic errors introduced by the streams. Convergence is usually attained with 2 to 4 iterations, with derived equilibrium parameters matching the prescribed axisymmetric core values to within estimated experimental uncertainties. Accurate recovery of core parameters occurs when the ratio of the net toroidal windup current from the streams to the core plasma current is less than 0.2, which is typically satisfied in most experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fast solvers for tokamak fluid models with PETSc

Multigrid (MG) is widely recognized as a highly effective solver for the model problem, the Laplacian, but textbook MG fails on most problems of interest. MG methods have been applied to complex, real-world applications with careful consideration of the physical model and discretization. In this work we develop the first step in applying MG methods to science and engineering relevant magnetohydrodynamics (MHD) tokamak models in the M3D-C1 (https://m3dc1.pppl.gov) fusion energy science code. The semi-implicit time integrator in M3D-C1 is composed of many linear solves. The implicit advance of the momentum equation is the most challenging and is the focus of this work. The current production solver in M3D-C1 is a block Jacobi (BJ) preconditioner within a Krylov solver, where blocks group degrees of freedom on planes of constant toroidal coordinate. BJ convergence degrades as the number of planes increases due to the spectral properties of the matrix preconditioned with BJ. The partially magnetic field-aligned, regular toroidal grid structure in M3D-C1 is amenable to semi-coarsening geometric MG in the toroidal direction. This paper develops such a solver and demonstrates competitive performance on a runaway electron model of a SPARC (https://cfs.energy/technology/sparc) disruption, and superior robustness on a stellarator model on which the BJ solver fails to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inferring the scrape-off layer heat flux width in a divertor with a low degree of axisymmetry

Plasma facing components (PFCs) in the next generation of tokamak devices will operate in challenging environments, with heat loads predicted to exceed 10 MWm -2 . The magnitude of these heat loads is set by the width of the channel, the ‘scrape-off layer’ (SOL), into which heat is exhausted, and can be characterised by an e-folding length scale for the decay of heat flux across the channel. It is expected this channel will narrow as tokamaks move towards reactor relevant conditions. Understanding the processes involved in setting the SOL heat flux width is imperative to be able to predict the heat loads PFCs must handle in future devices. Measurements of the SOL width are performed on the high-field spherical tokamak, ST40, using a newly commissioned infrared thermography system. With its high on-axis toroidal magnetic field (≥1.5 T) ST40 is uniquely positioned to investigate the influence of toroidal field on the heat flux width in spherical tokamaks, whilst also extending measurements of the SOL width in spherical tokamaks to increased poloidal field (≥0.3 T). Due to the divertor on ST40 having a low degree of axisymmetry, it is necessary for a set of radial measurements of the heat flux to be taken across the divertor, made possible using an automated toolchain that fully incorporates its 3D geometry. These radial profiles are combined with the magnetic geometry of the plasma to infer the width of the SOL, with both Eich and double exponential profiles of heat flux observed. A reduction in the heat flux is observed toroidally across part of the divertor, along with increased heat loads observed locally around the edges of the tiles. Future work in characterising the impact of tile misalignment and uncertainties in the reconstructed divertor magnetic geometry is required in order to further understand the observed heat flux patterns, as are additional investigations into the role potentially being played by an inhomogeneous sheath electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Full-torus impurity transport simulation in boron powder injection experiments in the Large Helical Device

The toroidal distribution of boron deposition on plasma-facing components (PFCs) in boron powder injection using an impurity power dropper (IPD) was investigated by full-torus simulation and observations in a systematic plasma density-scan experiment. The images of the ablation of dropped boron powders observed with a visible CCD camera were consistently explained by the simulations of the ablation positions of the boron powders considering the size distribution. Simulations assuming full-torus boron deposition on the PFCs did not reproduce the observed intensity profile of boron emission lines for higher plasma densities. It indicated that the density of boron deposited on PFCs installed toroidally far from the IPD was low for higher plasma densities due to the change in the ablation positions of the boron powders toward the outboard side. The experimental results verified the previous full-torus simulation of the toroidal distribution of the boron deposition in both lower and higher plasma densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tungsten erosion during L-mode discharges in the DIII-D SAS-VW divertor

A series of L-mode plasma discharges was performed in the DIII-D tokamak to assess the impact of outer strike point (OSP) position and toroidal magnetic field direction on erosion and core contamination potential of the recently-installed, tungsten-coated Small Angle Slot (SAS-VW) divertor. In one discharge, in-slot emission spectroscopy measured an<48 % increase in the W gross erosion rate when the OSP was moved 3 cm outwards, away from the V-shaped vertex of the slot divertor. However, the effective W yield (erosion rate divided by the incident D flux) was, overall, insensitive to changes in OSP location. Consistently low estimates of the effective W yield based on measurements taken a few cm outwards from the vertex suggest potentially significant C surface contamination. No W emission signal was detected when orienting the toroidal magnetic field such that the ion B×∇B drift direction is pointed away from the X-point. However, measurements of W content in the plasma core for both toroidal magnetic field directions suggest the presence of additional, unmeasured sources of erosion. The difference in the measured core W density with OSP position is much greater than the difference in the measured erosion rates, which may suggest that the leakage of eroded impurities out of the divertor is governed primarily through the parallel ion temperature gradient and friction forces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear coupling of tearing modes in reversed field pinch plasmas with stepped pressure profiles

A theory of the three-wave coupling of triplets of tearing modes in toroidal pinches [i.e., either reversed field pinches (RFPs) or tokamaks] was proposed by R. Fitzpatrick [Phys. Plasmas 6, 1168 (1999)]. However, this theory only applies to toroidal pinches with negligible equilibrium plasma pressure gradients. Such a limitation is particularly inappropriate to RFPs. This paper generalizes the analysis of R. Fitzpatrick [Phys. Plasmas 6, 1168 (1999)] in order to take the equilibrium pressure gradient into account. However, for the sake of simplicity, a stepped pressure profile, rather than a continuous profile, is employed. In the limit in which the number of steps becomes very large, the results obtained from the generalized theory are presumably equivalent to those that would have been achieved using a continuous pressure profile. The generalized theory is used to investigate the formation of the characteristic toroidally localized pattern of phase-locked m = 1 and m = 0 tearing modes in RFP plasmas that is known as the “slinky” pattern. The incorporation of the equilibrium plasma pressure into the analysis is found to be of crucial importance when determining the properties of the pattern. This is the case because the plasma pressure controls the number of unstable m = 1 and m = 0 tearing modes, and also significantly affects the strength of three-wave coupling, as well as the phase relation between the phase-locked m = 1 and m = 0 modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resistive hose modes in tokamak runaway electron beams II

Resistive hose instabilities of runaway electron (RE) beams immersed in resistive background plasmas are examined with analytic and numerical calculations. The RE beam-plasma equilibria considered are characteristic of the situation observed post-thermal quench in a tokamak disruption. An analytic linear dispersion relation is presented for the case of a uniform RE current density profile with a sharp boundary in cylindrical geometry. Initial value linear calculations for a more general profile in toroidal geometry find that reducing aspect ratio increases the resistive hose mode growth rates with fixed safety factor profile. Nonlinear calculations in cylindrical and toroidal geometry find that the resistive hose instability-driven fluctuations relax the gradient of the current density profile. In toroidal geometry, changes to the magnetic topology are observed as a result of the resistive hose activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On resistive interchange, double tearing, and resonant and non-resonant infernal modes in spherical tokamak negative central shear discharges

We examine the linear and nonlinear stability of a sequence of reversed shear NSTX equilibria with the same toroidal current density and pressure profile but with different toroidal field strengths. All equilibria have two q = 2 surfaces and have q > 1 everywhere. For this sequence, which all have β P ∼ 0.5, as the minimum value of q decreases below about 1.4, the localized resistive interchange criterion is strongly violated, and the unstable mode with toroidal mode number n=1 changes its dominant poloidal mode number m from a (m,n)=(2,1) double tearing mode to a non-resonant (1,1) infernal mode. This (1,1) mode nonlinearly flattens both the current density and the pressure near the magnetic axis. The higher-n unstable modes are all resistive infernal modes localized around the surface where q has a minimum and the magnetic shear vanishes. They generally saturate nonlinearly at low amplitude but can cause magnetic surface breakup near the axis.

Spitzer resistivity↗

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗

Validation study of RWM stability in DIII-D high- β N plasmas

The n = 1 (n is the toroidal mode number) resistive wall mode (RWM) stability is numerically investigated for two DIII-D high-β N discharges 176440 and 172461, utilizing the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) codes. Systematic validation efforts are attempted, for the first time, for discharges with very slow or vanishing toroidal flow for a large fraction of the plasma volume. While gaining physics insights in accessing stable operation regime at β N exceeding the Troyon no-wall limit in these slow-rotation experiments, the predictive capability of fluid and non-perturbative magnetohydrodynamic-kinetic hybrid models for the RWM is further confirmed. The MARS-F fluid model, with a strong but numerically tunable viscosity mimicking ion Landau damping of parallel sound waves, finds complete stabilization of the n = 1 RWM in the considered DIII-D plasmas under the experimental flow conditions. Similarly, either full stabilization (for discharge 176440) or marginal stability (for discharge 172461) of the mode is computed by the MARS-K hybrid model, which is first-principle based without free model parameters. In particular, all drift kinetic resonances, including those of thermal and energetic particles, are found to synergistically act to marginally stabilize the RWM in discharge 172461. These MARS-F/K modeling results explain the experimentally observed stable operational regime in DIII-D, as far as the RWM stability is concerned. Extensive numerical sensitivity studies, with respect to the plasma toroidal flow speed as well as the radial location of the resistive wall, are also carried out to further support the validation study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extended-MHD simulations of disruption mitigation via massive gas injection in SPARC

Recent developments to the M3D-C1 code enable higher fidelity modeling of disruptions, and can be applied in the design verification of reactor-scale tokamaks. Among these new capabilities is a method to mesh conducting vessel structures such as coils and passive plates, packing of the toroidal mesh around gas injectors, as well as anisotropic resistivity inside the vessel structures. We present extended-magnetohydrodynamic (MHD) simulations of disruption mitigation via massive gas injection (MGI) in SPARC. The goal of this study is to inform the disruption mitigation layout of SPARC and aid in the design of an effective gas injector configuration. Fully three-dimensional simulations with M3D-C1 are carried out for various injector configurations with the primary goal of determining the effect of different MGI parameters on heat loads and vessel forces. The simulations include a model for impurity ionization, recombination, advection and radiation, as well as spatially resolved conducting structures around the plasma. A localized mixture of deuterium and neon with a small toroidal and poloidal width is injected in up to six locations. We demonstrate that M3D-C1 can model a rapid shutdown via MGI using narrow and more realistic gas plumes than in previous simulations. As a result of the q = 1 surface in the SPARC baseline case a sawtooth is observed early in the simulations. Despite the sawtooth and the onset of edge MHD instabilities, the impurity distribution remains localized around the injector locations, but enables a radiative shutdown of the plasma. We find that using the maximum of six gas injectors results in a lower peaking factor and leads to a more even distribution of radiation toroidally than using two injectors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling resistive-inductive evolution of currents in Wendelstein 7-X

This research investigates the temporal evolution of the toroidal plasma current in the Wendelstein 7-X (W7-X) stellarator under different heating, fueling, and current drive scenarios. The THRIFT code has been modernized and its predictions of the evolution of the toroidal current have been compared against experimentally measured currents in W7-X. Good agreement is found with respect to the characteristic timescale between experimentally measured and simulated toroidal currents. The total bootstrap current is under-predicted owing to the applicability of the BOOTSJ model for the plasma collisionalities in question. Edge plasma resistivity is found to play an important role in the asymptotic behavior of the evolution of the current, indicating a possible limitation of the minimum plasma temperature when applying this model. Simulations of ECCD and heating power steps show THRIFT is capable of capturing the dynamical evolution of the current in response to changes in current sources. Future uses of THRIFT include validating and benchmarking other non-inductive current models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

3D modeling of n = 1 RMP driven heat fluxes on the SPARC tokamak PFCs using HEAT

3D heat flux calculations at the lower outer divertor plate of SPARC using the HEAT code show that 3D fields generated from error field correction coils can lead to enhanced peak heat fluxes up to 15 times larger compared to the axisymmetric case. Previously employed to simulate axisymmetric heat flux on 3D plasma facing components, the HEAT code can now predict 3D heat flux generated by non-axisymmetric plasmas. This is achieved via a new HEAT module which leverages the 3D field line tracing capabilities of MAFOT starting from an M3D-C1 (MHD resistive code) perturbed equilibria. The resulting heat flux is assigned using the magnetic footprint and the heat flux layer model, an extension of the 2D heat flux model also known as the Eich, to 3D non-axisymmetric plasmas. For SPARC, the new capabilities of HEAT are used to calculate the 3D heat loads resulting from n = 1 perturbation fields (with n indicating the toroidal periodicity) applied through a toroidal array of six picture frame coils with different amplitude. The comparison with the unperturbed case shows significant changes in shape and intensity of the heat flux profile. The results show that the application of n = 1 3D field leads to a localized enhancement of the heat flux peak, influenced by the wetted area impacted by the magnetic footprint, and the appearance of a secondary heat flux peak, whose intensity depends on amplitude of the applied 3D field and toroidal location.

3D heat flux↗

MARS-Q modeling of low- n resistive kink-peeling modes and edge harmonic oscillations in DIII-D

Linear and quasilinear magnetohydrodynamic (MHD) modeling is carried out for two DIII-D discharges that both featured a transition from quiescent H-mode (QH) to wide-pedestal QH (WPQH). The MHD perturbations, associated with the edge harmonic oscillations (EHOs) observed during the QH-phase in both discharges, are identified as low-n (n is the toroidal mode number) resistive kink-peeling instabilities, with the ideal MHD counterpart remaining stable. The quasilinear model successfully simulates EHO-like perturbations during the QH phase in both discharges, confirming the experimental observations. A less intuitive finding is the EHO-like behavior involving the n = 2 perturbation, simulated for one of the discharges during the WPQH-phase. This result, while consistent with experimental observations, is obtained despite the fact that the initial perturbation is linearly stable. The eventual growth of the perturbation is solely due to nonlinear interaction between the MHD perturbation and the plasma toroidal flow. The occurrence of the EHO-like perturbation during the WPQH-phase is found to be sensitive to the initial profile of the plasma edge rotation. For DIII-D plasmas considered, the neoclassical toroidal viscosity, generated by three-dimensional low-n perturbations, is found to play a dominant role in modifying the edge flow during EHOs.

EHO↗