Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MAGNETIC CONFINEMENT”

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 1,729 records · Page 96

Asymptotic quasisymmetric high-beta three-dimensional magnetohydrodynamic equilibria near axisymmetry

Quasisymmetry (QS), a hidden symmetry of the magnetic field strength, is known to support nested flux surfaces and provide superior particle confinement in stellarators. In this work, we study the ideal magnetohydrodynamic (MHD) equilibrium and stability of high-beta plasma in a large-aspect-ratio stellarator. In particular, we show that the lowest-order description of a near-axisymmetric equilibrium vastly simplifies the problem of three-dimensional quasisymmetric MHD equilibria, which can be reduced to a standard elliptic Grad–Shafranov equation for the flux function. We show that any large-aspect-ratio tokamak, deformed periodically in the vertical direction, is a stellarator with approximate volumetric QS. We discuss exact analytical solutions and numerical benchmarks. Finally, we discuss the ideal ballooning and interchange stability of some of our equilibrium configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the toroidal current density flowing across a poloidal-magnetic-field null in an axisymmetric plasma

The axisymmetry condition and two of Maxwell's equations are used to show that, in general, there are no nested magnetic surfaces around a poloidal-magnetic-field null for a sufficiently small value of the toroidal current density flowing there. Hence, the toroidal current density at the axis of a magnetic configuration with extreme shear reversal cannot continuously approach zero unless nested surfaces are first broken or particular values are assigned to boundary conditions and other plasma parameters. The threshold of the toroidal current-density at which the topology changes is shown to be set by such parameters, and some examples of the predicted topology transition are presented using analytical solutions of the Grad-Shafranov equation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stability of compressible reduced magnetohydrodynamic equilibria-Analogy with magnetorotational instability

Stability analyses for equilibria of the compressible reduced magnetohydrodynamics (CRMHD) model are carried out by means of the Energy-Casimir (EC) method. Stability results are compared with those obtained for ideal magnetohydrodynamics (MHD) from the classical {delta}W criterion. An identification of the terms in the second variation of the free energy functional for CRMHD with those of {delta}W is made: two destabilizing effects present for CRMHD turn out to correspond to the kink and interchange instabilities in usual MHD, while the stabilizing roles of field line bending and compressibility are also identified in the reduced model. Also, using the EC method, stability conditions in the presence of toroidal flow are obtained. A formal analogy between CRMHD and a reduced incompressible model for magnetized rotating disks, due to Julien and Knobloch [EAS Pub. Series, 21, 81 (2006)], is discovered. In light of this analogy, energy stability analysis shows that the condition for magnetorotational instability (MRI) for the latter model corresponds to the condition for interchange instability in CRMHD, with the Coriolis term and shear velocity playing the roles of the curvature term and pressure gradient, respectively. Using the EC method, stability conditions for the rotating disk model, for a large class of equilibria with possible non-uniform magnetic fields, are obtained. In particular, this shows it is possible for the MRI system to undergo, in addition to the MRI, another instability that is analogous to the kink instability. For vanishing magnetic field, the Rayleigh hydrodynamical stability condition is recovered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-consistent simulations of nonlinear magnetohydrodynamics and profile evolution in stellarator configurations

Self-consistent extended MHD framework is used to investigate nonlinear macroscopic dynamics of stellarator configurations. In these calculations, initial conditions are given by analytical 3-D vacuum solutions. Finite beta discharges in a straight stellarator are simulated. Vacuum magnetic fields are applied to produce stellarator-like rotational transform profiles with iota(0) ≤ 0.5 and iota(0) ≥ 0.5. The vacuum magnetic fields are either helically symmetric or spoiled by the presence of magnetic harmonics of incommensurate helicity. As heat is added to the system, pressure-driven instabilities are excited when a critical β is exceeded. These instabilities may grow to large amplitude and effectively terminate the discharge, or they may saturate nonlinearly as the configuration evolves. In all of these studies, anisotropic heat conduction is allowed with κ{sub ∥}/κ{sub ⊥}=10{sup 4}−10{sup 7}.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparisons of linear and nonlinear plasma response models for non-axisymmetric perturbations

With the installation of non-axisymmetric coil systems on major tokamaks for the purpose of studying the prospects of ELM-free operation, understanding the plasma response to the applied fields is a crucial issue. Application of different response models, using standard tools, to DIII-D discharges with applied non-axisymmetric fields from internal coils, is shown to yield qualitatively different results. The plasma response can be treated as an initial value problem, following the system dynamically from an initial unperturbed state, or from a nearby perturbed equilibrium approach, and using both linear and nonlinear models [A. D. Turnbull, Nucl. Fusion 52, 054016 (2012)]. Criteria are discussed under which each of the approaches can yield a valid response. In the DIII-D cases studied, these criteria show a breakdown in the linear theory despite the small 10{sup −3} relative magnitude of the applied magnetic field perturbations in this case. For nonlinear dynamical evolution simulations to reach a saturated nonlinear steady state, appropriate damping mechanisms need to be provided for each normal mode comprising the response. Other issues arise in the technical construction of perturbed flux surfaces from a displacement and from the presence of near nullspace normal modes. For the nearby equilibrium approach, in the absence of a full 3D equilibrium reconstruction with a controlled comparison, constraints relating the 2D system profiles to the final profiles in the 3D system also need to be imposed to assure accessibility. The magnetic helicity profile has been proposed as an appropriate input to a 3D equilibrium calculation and tests of this show the anticipated qualitative behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of magnetohydrodynamic modes during sawteeth in tokamak plasmas

A dynamical analysis applied to a reduced resistive magnetohydrodynamics model is shown to explain the chronology of the nonlinear destabilization of modes observed in tokamak sawteeth. A special emphasis is put on the nonlinear self-consistent perturbation of the axisymmetric m = n = 0 mode that manifests through the q-profile evolution. For the very low fusion-relevant resistivity values, the q-profile is shown to remain almost unchanged on the early nonlinear timescale within the central tokamak region, which supports a partial reconnection scenario. Within the resistive region, indications for a local flattening or even a local reversed-shear of the q-profile are given. The impact of this ingredient in the occurrence of the sawtooth crash is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-bandwidth polarimeter for a high density, accelerated spheromak

A three-beam heterodyne polarimeter has been built to diagnose the PI-1 plasma injector at General Fusion, Inc. The polarimeter measures plasma density and Faraday rotation, which can provide estimates of magnetic field magnitude and direction. Two important calibration steps are required for the polarimeter to produce reasonable Faraday rotation signals. Beam polarization ellipticity must be measured to ensure there is a linear Faraday rotation response. In addition, the two beams travelling through the plasma must be collinear to reduce error due to differences in plasma density. Once calibrated, the Faraday rotation signals are in much better agreement with other diagnostics. For a null signal, the Faraday rotation measurement noise floor is 0.1° at 0.5 MHz bandwidth. Comparing preliminary spheromak Faraday rotation measurements to a model, the maximum error is about 0.3° at 0.5 MHz bandwidth, which is primarily due to electrical noise during the injector's capacitor discharge and limitations of the model. At a bandwidth of 0.5 MHz, the polarimeter has an axial resolution between 6 cm and 30 cm depending on the speed of the spheromak, which varies between 30 km/s and 150 km/s. The spheromak length ranges from 0.75 m to 2 m. Additional polarimeter chords will be added in future upgrades.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Synchrotron radiation from a runaway electron distribution in tokamaks

The synchrotron radiation emitted by runaway electrons in a fusion plasma provides information regarding the particle momenta and pitch-angles of the runaway electron population through the strong dependence of the synchrotron spectrum on these parameters. Information about the runaway density and its spatial distribution, as well as the time evolution of the above quantities, can also be deduced. In this paper, we present the synchrotron radiation spectra for typical avalanching runaway electron distributions. Spectra obtained for a distribution of electrons are compared with the emission of mono-energetic electrons with a prescribed pitch-angle. We also examine the effects of magnetic field curvature and analyse the sensitivity of the resulting spectrum to perturbations to the runaway distribution. The implications for the deduced runaway electron parameters are discussed. We compare our calculations to experimental data from DIII-D and estimate the maximum observed runaway energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resistive magnetohydrodynamic simulations of helicity-injected startup plasmas in National Spherical Torus eXperiment

The generation of helicity-injected startup plasmas in National Spherical Torus eXperiment (NSTX), including flux surface closure, is studied using resistive-magnetohydrodynamic simulations with plasma flows, currents, ohmic heating and anisotropic thermal conduction. An injection-voltage pulse shape is used that separates the injection and closure phases allowing elucidation of the physics. The formation of an X-point near the helicity-injection gap is triggered as the injector voltage drops to zero. Near the forming X-point, magnetic pressure due to toroidal field entrained in the E × B plasma flow from the helicity-injection gap drops, allowing resistive magnetic reconnection even though the total injected current is almost constant. Where appropriate, the simulations are compared with Transient Coaxial Helicity Injection experiments in the NSTX spherical tokamak, which have demonstrated the formation of a promising candidate for non-inductive startup plasmas [Raman et al., Phys. Rev. Lett. 90, 075005 (2003)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Axisymmetric curvature-driven instability in a model divertor geometry

A model problem is presented which qualitatively describes a pressure-driven instability which can occur near the null-point in the divertor region of a tokamak where the poloidal field becomes small. The model problem is described by a horizontal slot with a vertical magnetic field which plays the role of the poloidal field. Line-tying boundary conditions are applied at the planes defining the slot. A toroidal field lying parallel to the planes is assumed to be very strong, thereby constraining the possible structure of the perturbations. Axisymmetric perturbations which leave the toroidal field unperturbed are analyzed. Ideal magnetohydrodynamics is used, and the instability threshold is determined by the energy principle. Because of the boundary conditions, the Euler equation is, in general, non-separable except at marginal stability. This problem may be useful in understanding the source of heat transport into the private flux region in a snowflake divertor which possesses a large region of small poloidal field, and for code benchmarking as it yields simple analytic results in an interesting geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum confinement of a hydrogenic donor in a double quantum well: Through diamagnetic susceptibility

Diamagnetic susceptibility of a randomly distributed donor in a GaAs/Al{sub 0.3}Ga{sub 0.7}As Double Quantum Well has been calculated in its ground state as a function of barrier and well width. It is shown that the modification in the barrier and well dimension significantly influences the dimensional character of the donor through modulating the subband distribution and in turn the localization of the donor. The effect of barrier and well thickness on the interparticle distance has also been observed. Interestingly it opens up the possibility of tuning the susceptibility and monitoring the tunnel coupling among the wells.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effects of dust contamination on the transverse dynamics of a magnetized electron plasma

Complex (dusty) plasmas are characterized by the presence of a fraction of micrometric or sub-micrometric particles which may collect a surface charge up to the order of a few thousand electron charges. The dusty plasmas studied in the experiments generally satisfy a global neutrality condition. By contrast, we present here the investigation of a magnetized nonneutral plasma, i.e., a plasma with a single sign of charge (e.g. electrons) confined in a Penning-Malmberg trap, contaminated by a dust population. We simulate the two-dimensional transverse dynamics of this multi-component plasma with a particle-in-cell code implementing a mass-less fluid (drift-Poisson) approximation for the electrons and a kinetic description for the dust component (including gravity). Simulations with different initial dust distributions and densities have been performed in order to investigate the influence of the dust on the development of the diocotron instability in the electron plasma. In particular, the early stage of the growth of the diocotron modes has been analyzed by Fourier decomposition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Using motional Stark splitting of D$_α$ emission to constrain MHD equilibrium analysis in DIII-D plasmas

Herein we report tests of an alternate technique for constraining MHD equilibrium analysis in tokamak plasmas using internal magnetic field measurements based on |$B$| measurements from the motional Stark splitting of D$_α$ spectral lines emitted by a neutral heating beam (MSE-LS). We compare results using MSE-LS with those of the standard equilibrium analysis technique based on line polarization of the D$_α$ emission (MSE-LP). An alternative to MSE-LP is needed in future devices such as ITER where MSE-LP will be difficult due to a plasma-induced coating of the first optical element. The tests utilized data from 10 DIII-D shots with 7 MSE-LS and 14 MSE-LP views covering a range of radii along the outer midplane of the plasma. Seven MSE-LS measurements can contribute significantly to the equilibrium reconstruction of pressure and q profiles using both synthetic and experimental DIII-D MSE-LS data. For example, 7 MSE-LS plus seven MSE-LP measurements give a fit quality that is as good as the same cases with 14 MSE-LP measurements. Analyzing synthetic data for 14 MSE-LS measurements shows significant improvement in fitting quality over the case with 7 MSE-LS locations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of plasmoid formation on sawtooth process in a tokamak

For realistic values of the Lundquist number in tokamak plasmas, 1/1 magnetic islands lead to the formation of a secondary thin current sheet, which breaks up into a chain of small magnetic islands, called plasmoids. The role of plasmoid dynamics during the sawtooth reconnection process in fusion plasmas remains an unresolved issue. In this study, systematic simulations are performed to investigate the resistive internal kink mode using the full resistive magnetohydrodynamics equations implemented in the NIMROD code in a simplified tokamak geometry. For Lundquist number S≥1.6×10 7 , the secondary current sheet is found to be unstable to plasmoids during the nonlinear resistive kink mode evolution with a critical aspect ratio of the current sheet of ~70. The merging of small plasmoids leads to the formation of a monster plasmoid that can significantly affect the primary island evolution. This may provide an explanation for the partial reconnection observed in sawtooth experiments.

Computer simulation↗

Artificially structured boundary plasma trap

A computer simulation is presented of single-species non-neutral plasma confinement using an artificially structured boundary. The artificially structured boundary produces a spatially periodic static electromagnetic field along the plasma periphery such that the spatial period of the applied field is much smaller than the dimensions of the confined plasma. The simulated non-neutral plasma self-consistently produces an electrostatic potential energy well for oppositely signed charged particles. The results support the prospect of developing plasma space-charge based confinement, with an unmagnetized plasma of one species of charged particles confined by an electric field produced by an edge-confined plasma of a second species of charged particles. The Warp particle-in-cell code is used for the simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Anti-symmetric representation of the extended magnetohydrodynamic equations

We introduce the anti-symmetric representation of the extended magnetohydrodynamic (MHD) equations. In this representation, the use of the anti-symmetric flux operator (∇·ν+ν·∇) results in conservation theorems with discrete analogs. Inherently robust numerical applications are achieved with little effort, and conservation to machine precision is possible with simple numerical schemes. Starting from the two-fluid equations, we construct a single-fluid MHD model based on generalized center-of-mass variables for the mass (√ρ), momentum (√ρν), and pressure (√ p ). This model is shown to possess identical conservation properties to the two-fluid system, with the only restriction being the use of a single temperature. Common approximations to the Braginskii heat fluxes and to the gyroviscous stress tensor are cast into our representation for convenience. The discrete conservation properties are verified using the classic Orszag–Tang vortex problem. In addition to the favorable mass, momentum, and energy conservation properties, the time reversibility of the simulations is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron heating in 2D: Combining Fermi–Ulam acceleration and magnetic-moment non-adiabaticity in a mirror-configuration plasma

We analyze a new mechanism for the creation and confinement of energetic electrons in a mirror-configuration plasma. A Fermi–Ulam-type process, driven by end-localized coherent electrostatic oscillations, provides axial acceleration, while a natural non-adiabaticity of μ provides phase decorrelation and energy isotropization. This novel 2D combination causes the electron energy distribution function, calculated with a diffusive-loss model, to assume a Maxwellian shape with the μ non-adiabaticity, reducing loss-cone escape and annulling the absolute-barrier energy-limiting Chirikov criterion of lower dimensional models. Here, the theoretical predictions are compared with data from an experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of synthetic electron cyclotron emission from the high field side of HL-2M tokamak plasmas

A synthetic electron cyclotron emission (ECE) diagnostic is used to interpret ECE signals from preset plasma equilibrium profiles, including magnetic field, electron density, and electron temperature. According to the simulation results, the electron temperature (T e ) profile covering the harmonic overlap region can be obtained by receiving ECE signals at the high field side (HFS) of the HL-2M plasma. The third harmonic ECE at the low field side (LFS) cannot pass through the second harmonic resonance layer at the HFS unless the optical thickness (τ) of the second harmonic becomes gray (τ ≤ 2). In addition, the impact of the relativistic frequency down-shift has been evaluated and corrected. The measurable range of the HFS ECE has been calculated by scanning different parameters (electron density, temperature, and magnetic field). Higher plasma parameters allow a wider radial range of electron temperature measurements. The minimum inner measurable position can reach R = 120 cm (r/a = -0.89) when the product of core temperature (T e0 ) and density (n e0 ) is greater than 35 × 10 19 keV m -3 , which is extended by more than 30 cm inward compared with that of the LFS measurement. The HFS ECE will greatly improve the diagnostic ability of ECE systems on the HL-2M tokamak.

47 OTHER INSTRUMENTATION↗