Showing papers by "Princeton Plasma Physics Laboratory published in 2018"

Space micropropulsion systems for Cubesats and small satellites: from proximate targets to furthermost frontiers

Abstract: Rapid evolution of miniaturized, automatic, robotized, function-centered devices has redefined space technology, bringing closer the realization of most ambitious interplanetary missions and intense near-Earth space exploration. Small unmanned satellites and probes are now being launched in hundreds at a time, resurrecting a dream of satellite constellations, i.e., wide, all-covering networks of small satellites capable of forming universal multifunctional, intelligent platforms for global communication, navigation, ubiquitous data mining, Earth observation, and many other functions, which was once doomed by the extraordinary cost of such systems. The ingression of novel nanostructured materials provided a solid base that enabled the advancement of these affordable systems in aspects of power, instrumentation, and communication. However, absence of efficient and reliable thrust systems with the capacity to support precise maneuvering of small satellites and CubeSats over long periods of deployment remains a real stumbling block both for the deployment of large satellite systems and for further exploration of deep space using a new generation of spacecraft. The last few years have seen tremendous global efforts to develop various miniaturized space thrusters, with great success stories. Yet, there are critical challenges that still face the space technology. These have been outlined at an inaugural International Workshop on Micropropulsion and Cubesats, MPCS-2017, a joint effort between Plasma Sources and Application Centre/Space Propulsion Centre (Singapore) and the Micropropulsion and Nanotechnology Lab, the G. Washington University (USA) devoted to miniaturized space propulsion systems, and hosted by CNR-Nanotec—P.Las.M.I. lab in Bari, Italy. This focused review aims to highlight the most promising developments reported at MPCS-2017 by leading world-reputed experts in miniaturized space propulsion systems. Recent advances in several major types of small thrusters including Hall thrusters, ion engines, helicon, and vacuum arc devices are presented, and trends and perspectives are outlined.

Magnetic configuration effects on the Wendelstein 7-X stellarator

Abstract: The two leading concepts for confining high-temperature fusion plasmas are the tokamak and the stellarator. Tokamaks are rotationally symmetric and use a large plasma current to achieve confinement, whereas stellarators are non-axisymmetric and employ three-dimensionally shaped magnetic field coils to twist the field and confine the plasma. As a result, the magnetic field of a stellarator needs to be carefully designed to minimize the collisional transport arising from poorly confined particle orbits, which would otherwise cause excessive power losses at high plasma temperatures. In addition, this type of transport leads to the appearance of a net toroidal plasma current, the so-called bootstrap current. Here, we analyse results from the first experimental campaign of the Wendelstein 7-X stellarator, showing that its magnetic-field design allows good control of bootstrap currents and collisional transport. The energy confinement time is among the best ever achieved in stellarators, both in absolute figures (τE > 100 ms) and relative to the stellarator confinement scaling. The bootstrap current responds as predicted to changes in the magnetic mirror ratio. These initial experiments confirm several theoretically predicted properties of Wendelstein 7-X plasmas, and already indicate consistency with optimization measures. Results from the first experimental campaign of the Wendelstein 7-X stellarator demonstrate that its magnetic-field design grants good control of parasitic plasma currents, leading to long energy confinement times.

A fast low-to-high confinement mode bifurcation dynamics in the boundary-plasma gyrokinetic code XGC1

Abstract: A fast edge turbulence suppression event has been simulated in the electrostatic version of the gyrokinetic particle-in-cell code XGC1 in a realistic diverted tokamak edge geometry under neutral particle recycling. The results show that the sequence of turbulent Reynolds stress followed by neoclassical ion orbit-loss driven together conspire to form the sustaining radial electric field shear and to quench turbulent transport just inside the last closed magnetic flux surface. The main suppression action is located in a thin radial layer around ψN≃0.96–0.98, where ψN is the normalized poloidal flux, with the time scale ∼0.1 ms.

3D field phase-space control in tokamak plasmas

Abstract: A small relaxation of the axisymmetric magnetic field of a tokamak into a non-axisymmetric three-dimensional (3D) configuration can be effective to control magnetohydrodynamic instabilities, such as edge-localized modes. However, a major challenge to the concept of 3D tokamaks is that there are virtually unlimited possible choices for a 3D magnetic field, and most of them will only destabilize or degrade plasmas by symmetry breaking. Here, we demonstrate the phase-space visualization of the full 3D field-operating windows of a tokamak, which allows us to predict which configurations will maintain high confinement without magnetohydrodynamic instabilities in an entire region of plasmas. We test our approach at the Korean Superconducting Tokamak Advanced Research (KSTAR) facility, whose 3D coils with many degrees of freedom in the coil space make it unique for this purpose. Our experiments show that only a small subset of coil configurations can accomplish edge-localized mode suppression without terminating the discharge with core magnetohydrodynamic instabilities, as predicted by the perturbative 3D expansion of plasma equilibrium and the optimizing principle of local resonance. The prediction provided excellent guidance, implying that our method can substantially improve the efficiency and fidelity of the 3D optimization process in tokamaks. A theoretical and numerical approach, validated by experiments at the KSTAR facility, shows how magnetohydrodynamic instabilities in tokamak plasmas can be efficiently controlled by a small relaxation of the confining field into a 3D configuration.

The Twisted Configuration of the Martian Magnetotail: MAVEN Observations

Abstract: Measurements provided by the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft are analyzed to investigate the Martian magnetotail configuration as a function of interplanetary magnetic field (IMF) BY. We find that the magnetotail lobes exhibit a ~45deg twist, either clockwise or counterclockwise from the ecliptic plane, up to a few Mars radii downstream. Moreover, the associated cross-tail current sheet is rotated away from the expected location for a Venus-like induced magnetotail based on nominal IMF draping. Data-model comparisons using magnetohydrodynamic simulations are in good agreement with the observed tail twist. Model field line tracings indicate that a majority of the twisted tail lobes are composed of open field lines, surrounded by draped IMF. We infer that dayside magnetic reconnection between the crustal fields and draped IMF creates these open fields and may be responsible for the twisted tail configuration, similar to what is observed at Earth.

First Direct Observation of Runaway-Electron-Driven Whistler Waves in Tokamaks.

Abstract: DIII-D experiments at low density (${n}_{e}\ensuremath{\sim}{10}^{19}\text{ }\text{ }{\mathrm{m}}^{\ensuremath{-}3}$) have directly measured whistler waves in the 100--200 MHz range excited by multi-MeV runaway electrons. Whistler activity is correlated with runaway intensity (hard x-ray emission level), occurs in novel discrete frequency bands, and exhibits nonlinear limit-cycle-like behavior. The measured frequencies scale with the magnetic field strength and electron density as expected from the whistler dispersion relation. The modes are stabilized with increasing magnetic field, which is consistent with wave-particle resonance mechanisms. The mode amplitudes show intermittent time variations correlated with changes in the electron cyclotron emission that follow predator-prey cycles. These can be interpreted as wave-induced pitch angle scattering of moderate energy runaways. The tokamak runaway-whistler mechanisms have parallels to whistler phenomena in ionospheric plasmas. The observations also open new directions for the modeling and active control of runaway electrons in tokamaks.

Liquid metals as a divertor plasma-facing material explored using the Pilot-PSI and Magnum-PSI linear devices

Abstract: For DEMO and beyond, liquid metal plasma-facing components are considered due to their resilience to erosion through flowed replacement, potential for cooling beyond conduction and inherent immunity to many of the issues of neutron loading compared to solid materials. The development curve of liquid metals is behind that of e.g. tungsten however, and tokamak-based research is currently somewhat limited in scope. Therefore, investigation into linear plasma devices can provide faster progress under controlled and well-diagnosed conditions in assessing many of the issues surrounding the use of liquid metals. The linear plasma devices Magnum-PSI and Pilot-PSI are capable of producing DEMO-relevant plasma fluxes, which well replicate expected divertor conditions, and the exploration of physics issues for tin (Sn) and lithium (Li) such as vapour shielding, erosion under high particle flux loading and overall power handling are reviewed here. A deeper understanding of erosion and deposition through this work indicates that stannane formation may play an important role in enhancing Sn erosion, while on the other hand the strong hydrogen isotope affinity reduces the evaporation rate and sputtering yields for Li. In combination with the strong redeposition rates, which have been observed under this type of high-density plasma, this implies that an increase in the operational temperature range, implying a power handling range of 20-25 MW m-2 for Sn and up to 12.5 MW m-2 for Li could be achieved. Vapour shielding may be expected to act as a self-protection mechanism in reducing the heat load to the substrate for off-normal events in the case of Sn, but may potentially be a continual mode of operation for Li.

Role of the Plasmoid Instability in Magnetohydrodynamic Turbulence

Abstract: The plasmoid instability in evolving current sheets has been widely studied due to its effects on the disruption of current sheets, the formation of plasmoids, and the resultant fast magnetic reconnection. In this Letter, we study the role of the plasmoid instability in two-dimensional magnetohydrodynamic (MHD) turbulence by means of high-resolution direct numerical simulations. At a sufficiently large magnetic Reynolds number (${R}_{m}={10}^{6}$), the combined effects of dynamic alignment and turbulent intermittency lead to a copious formation of plasmoids in a multitude of intense current sheets. The disruption of current sheet structures facilitates the energy cascade towards small scales, leading to the breaking and steepening of the energy spectrum. In the plasmoid-mediated regime, the energy spectrum displays a scaling that is close to the spectral index $\ensuremath{-}2.2$ as proposed by recent analytic theories. We also demonstrate that the scale-dependent dynamic alignment exists in 2D MHD turbulence and the corresponding slope of the alignment angle is close to 0.25.

Dual Phase-space Cascades in 3D Hybrid-Vlasov-Maxwell Turbulence

Abstract: To explain energy dissipation via turbulence in collisionless, magnetized plasmas, the existence of a dual real- and velocity-space cascade of ion-entropy fluctuations below the ion gyroradius has been proposed. Such a dual cascade, predicted by the gyrokinetic theory, has previously been observed in gyrokinetic simulations of two-dimensional, electrostatic turbulence. For the first time we show evidence for a dual phase-space cascade of ion-entropy fluctuations in a three-dimensional simulation of hybrid-kinetic, electromagnetic turbulence. Some of the scalings observed in the energy spectra are consistent with a generalized theory for the cascade that accounts for the spectral anisotropy of critically balanced, intermittent, sub-ion-Larmor-scale fluctuations. The observed velocity-space cascade is also anisotropic with respect to the magnetic-field direction, with linear phase mixing along magnetic-field lines proceeding mainly at spatial scales above the ion gyroradius and nonlinear phase mixing across magnetic-field lines proceeding at perpendicular scales below the ion gyroradius. Such phase-space anisotropy could be sought in heliospheric and magnetospheric data of solar-wind turbulence and has far-reaching implications for the dissipation of turbulence in weakly collisional astrophysical plasmas.

Modeling Martian Atmospheric Losses over Time: Implications for Exoplanetary Climate Evolution and Habitability

Abstract: In this Letter, we make use of sophisticated 3D numerical simulations to assess the extent of atmospheric ion and photochemical losses from Mars over time. We demonstrate that the atmospheric ion escape rates were significantly higher (by more than two orders of magnitude) in the past at $\sim 4$ Ga compared to the present-day value owing to the stronger solar wind and higher ultraviolet fluxes from the young Sun. We found that the photochemical loss of atomic hot oxygen dominates over the total ion loss at the current epoch whilst the atmospheric ion loss is likely much more important at ancient times. We briefly discuss the ensuing implications of high atmospheric ion escape rates in the context of ancient Mars, and exoplanets with similar atmospheric compositions around young solar-type stars and M-dwarfs.

Experimental conditions to suppress edge localised modes by magnetic perturbations in the ASDEX Upgrade tokamak

Abstract: Access conditions for full suppression of Edge Localised Modes (ELMs) by Magnetic Perturbations (MP) in low density high confinement mode (H-mode) plasmas are studied in the ASDEX Upgrade tokamak. The main empirical requirements for full ELM suppression in our experiments are: 1. The poloidal spectrum of the MP must be aligned for best plasma response from weakly stable kink-modes, which amplify the perturbation, 2. The plasma edge density must be below a critical value, $3.3 \\times 10^{19}$~m$^{-3}$. The edge collisionality is in the range $\ u^*_i = 0.15-0.42$ (ions) and $\ u^*_e = 0.15-0.25$ (electrons). However, our data does not show that the edge collisionality is the critical parameter that governs access to ELM suppression. 3. The pedestal pressure must be kept sufficiently low to avoid destabilisation of small ELMs. This requirement implies a systematic reduction of pedestal pressure of typically 30\\% compared to unmitigated ELMy H-mode in otherwise similar plasmas. 4. The edge safety factor $q_{95}$ lies within a certain window. Within the range probed so far, $q_{95}=3.5-4.2$, one such window, $q_{95}=3.57-3.95$ has been identified. Within the range of plasma rotation encountered so far, no apparent threshold of plasma rotation for ELM suppression is found. This includes cases with large cross field electron flow in the entire pedestal region, for which two-fluid MHD models predict that the resistive plasma response to the applied MP is shielded.

Benchmark of gyrokinetic, kinetic MHD and gyrofluid codes for the linear calculation of fast particle driven TAE dynamics

Abstract: Fast particles in fusion plasmas may drive Alfven modes unstable leading to fluctuations of the internal electromagnetic fields and potential loss of particles. Such instabilities can have an impact on the performance and the wall-load of machines with burning plasmas such as ITER. A linear benchmark for a toroidal Alfven eigenmode (TAE) is done with 11 participating codes with a broad variation in the physical as well as the numerical models. A reasonable agreement of around 20% has been found for the growth rates. Also, the agreement of the eigenfunctions and mode frequencies is satisfying. However, they are found to depend strongly on the complexity of the used model.

Correlation of the tokamak H-mode density limit with ballooning stability at the separatrix

Abstract: We show for JET and ASDEX Upgrade, based on Thomson-scattering measurements, a clear correlation of the density limit of the tokamak H-mode high-confinement regime with the approach to the ideal ballooning instability threshold at the periphery of the plasma. It is shown that the MHD ballooning parameter at the separatrix position increases about linearly with the separatrix density normalized to Greenwald density, for a wide range of discharge parameters in both devices. The observed operational space is found to reach at maximum –0.5 at values for –2.5, in the range of theoretical predictions for ballooning instability. This work supports the hypothesis that the H-mode density limit may be set by ballooning stability at the separatrix.

Plasma mass separation

Abstract: This tutorial describes mechanisms for separating ions in a plasma device with respect to their atomic or molecular mass for practical applications. The focus here is not on separating isotopes of a single atomic species but rather on systems with a much lower mass resolution and a higher throughput. These separation mechanisms include ion gyro-orbit separation, drift-orbit separation, vacuum arc centrifugation, steady-state rotating plasmas, and several other geometries. Generic physics issues are discussed such as the ion charge state, neutrals and molecules, collisions, radiation loss, and electric fields and fluctuations. Generic technology issues are also discussed such as plasma sources and ion heating, and suggestions are made for future research.

Role of Kinetic Instability in Runaway-Electron Avalanches and Elevated Critical Electric Fields

Abstract: The effects of kinetic whistler wave instabilities on the runaway-electron (RE) avalanche is investigated. With parameters from experiments at the DIII-D National Fusion Facility, we show that RE scattering from excited whistler waves can explain several poorly understood experimental results. We find an enhancement of the RE avalanche for low density and high electric field, but for high density and low electric field the scattering can suppress the avalanche and raise the threshold electric field, bringing the present model much closer to observations. The excitation of kinetic instabilities and the scattering of resonant electrons are calculated self-consistently using a quasilinear model and local approximation. We also explain the observed fast growth of electron cyclotron emission signals and excitation of very low-frequency whistler modes observed in the quiescent RE experiments at DIII-D tokamak. Simulations using ITER parameters show that by controlling the background thermal plasma density and temperature, the plasma waves can also be excited spontaneously in tokamak disruptions and the avalanche generation of runaway electrons may be suppressed.

Statistical Analysis of Torus and Kink Instabilities in Solar Eruptions

Abstract: A recent laboratory experiment of ideal magnetohydrodynamic (MHD) instabilities reveals four distinct eruption regimes readily distinguished by the torus instability (TI) and helical kink instability (KI) parameters \citep{Myers2015}. To establish its observational counterpart, we collect 38 solar flares (stronger than GOES class M5 in general) that took place within 45$^{\circ}$ of disk center during 2011$-$2017, 26 of which are associated with a halo or partial halo coronal mass ejection (CME) (i.e., ejective events), while the others are CMEless (i.e., confined events). This is a complete sample of solar events satisfying our selection criteria detailed in the paper. For each event, we calculate decay index $n$ of the potential strapping field above the magnetic flux rope (MFR) in and around the flaring magnetic polarity inversion line (a TI parameter), and the unsigned twist number $T_w$ of the non-linear force-free (NLFF) field lines forming the same MFR (a KI parameter). We then construct a $n-T_w$ diagram to investigate how the eruptiveness depends on these parameters. We find: (1) $T_w$ appears to play little role in discriminating between confined and ejective events; (2) the events with $n\gtrsim0.8$ are all ejective and all confined events have $n\lesssim0.8$. However, $n\gtrsim0.8$ is not a necessary condition for eruption, because some events with $n\lesssim0.8$ also erupted. In addition, we investigate the MFR's geometrical parameters, apex height and distance between footpoints, as a possible factor for the eruptiveness. We briefly discuss the difference of the present result for solar eruptions with that of the laboratory result in terms of the role played by magnetic reconnection.

New method to design stellarator coils without the winding surface

Abstract: Finding an easy-to-build coils set has been a critical issue for stellarator design for decades. Conventional approaches assume a toroidal 'winding' surface, but a poorly chosen winding surface can unnecessarily constrain the coil optimization algorithm, This article presents a new method to design coils for stellarators. Each discrete coil is represented as an arbitrary, closed, one-dimensional curve embedded in three-dimensional space. A target function to be minimized that includes both physical requirements and engineering constraints is constructed. The derivatives of the target function with respect to the parameters describing the coil geometries and currents are calculated analytically. A numerical code, named flexible optimized coils using space curves (FOCUS), has been developed. Applications to a simple stellarator configuration, W7-X and LHD vacuum fields are presented.

ELM elimination with Li powder injection in EAST discharges using the tungsten upper divertor

Abstract: We report the first successful use of lithium (Li) to eliminate edge-localized modes (ELMs) with tungsten divertor plasma-facing components in the EAST device. Li powder injected into the scrape-off layer of the tungsten upper divertor successfully eliminated ELMs for 3–5 s in EAST. The ELM elimination became progressively more effective in consecutive discharges at constant lithium delivery rates, and the divertor D α baseline emission was reduced, both signatures of improved wall conditioning. A modest decrease in stored energy and normalized energy confinement was also observed, but the confinement relative to H98 remained well above 1, extending the previous ELM elimination results via Li injection into the lower carbon divertor in EAST (Hu et al 2015 Phys. Rev. Lett. 114 055001). These results can be compared with recent observations with lithium pellets in ASDEX-Upgrade that failed to mitigate ELMs (Lang et al 2017 Nucl. Fusion 57 016030), highlighting one comparative advantage of continuous powder injection for real-time ELM elimination.

Core radial electric field and transport in Wendelstein 7-X plasmas

Abstract: The results from the investigation of neoclassical core transport and the role of the radial electric field profile (Er) in the first operational phase of the Wendelstein 7-X (W7-X) stellarator are presented In stellarator plasmas, the details of the Er profile are expected to have a strong effect on both the particle and heat fluxes Investigation of the radial electric field is important in understanding neoclassical transport and in validation of neoclassical calculations The radial electric field is closely related to the perpendicular plasma flow (u⊥) through the force balance equation This allows the radial electric field to be inferred from measurements of the perpendicular flow velocity, which can be measured using the x-ray imaging crystal spectrometer and correlation reflectometry diagnostics Large changes in the perpendicular rotation, on the order of Δu⊥∼ 5 km/s (ΔEr ∼ 12 kV/m), have been observed within a set of experiments where the heating power was stepped down from 2 MW to 06 MW The

Orchestrating TRANSP Simulations for Interpretative and Predictive Tokamak Modeling with OMFIT

Abstract: TRANSP simulations are being used in the OMFIT workflow manager to enable a machine-independent means of experimental analysis, postdictive validation, and predictive time-dependent simulations on ...

OMFIT Tokamak Profile Data Fitting and Physics Analysis

Abstract: One Modeling Framework for Integrated Tasks (OMFIT) has been used to develop a consistent tool for interfacing with, mapping, visualizing, and fitting tokamak profile measurements. OMFIT is used to...

Electron Physics in 3‐D Two‐Fluid 10‐Moment Modeling of Ganymede's Magnetosphere

Abstract: We studied the role of electron physics in 3D two-fluid 10-moment simulation of the Ganymede's magnetosphere. The model captures non-ideal physics like the Hall effect, the electron inertia, and anisotropic, non-gyrotropic pressure effects. A series of analyses were carried out: 1) The resulting magnetic field topology and electron and ion convection patterns were investigated. The magnetic fields were shown to agree reasonably well with in-situ measurements by the Galileo satellite. 2) The physics of collisionless magnetic reconnection were carefully examined in terms of the current sheet formation and decomposition of generalized Ohm's law. The importance of pressure anisotropy and non-gyrotropy in supporting the reconnection electric field is confirmed. 3) We compared surface \"brightness\" morphology, represented by surface electron and ion pressure contours, with oxygen emission observed by the Hubble Space Telescope (HST). The correlation between the observed emission morphology and spatial variability in electron/ion pressure was demonstrated. Potential extension to multi-ion species in the context of Ganymede and other magnetospheric systems is also discussed.

Grassy-ELM regime with edge resonant magnetic perturbations in fully noninductive plasmas in the DIII-D tokamak

Abstract: Resonant magnetic perturbations (n = 3 RMPs) are used to suppress large amplitude ELMs and mitigate naturally occurring 'grassy'-ELMs in DIII-D plasmas relevant to the ITER steady-state mission. Fully non-inductive discharges in the ITER shape and pedestal collisionality (nu(e)* approximate to 0.05-0.15) are routinely achieved in DIII-D with RMP suppression of the Type-I ELMs. The residual grassy-ELMs deliver a low peak heat flux to the divertor as low as 1.2 x the inter-ELM heat flux in plasmas with sustained high H-factor (H-98y2 approximate to 1.2). The operating window for the RMP grassy-ELM regime is q(95) = 5.3-7.1 and external torque in the range 9-0.7 Nm in the co-Ip direction, which is in the range required for a steady-state tokamak reactor. The RMP grassy-ELM regime is associated with a two-step pedestal, with strong flattening of the density around the zero crossing in the E x B shear. The edge magnetic response of the plasma to the n = 3 RMP is found to be approximate to 2-3x larger than for comparable ITER baseline plasmas (beta(N) approximate to 1.8, q(95) approximate to 3.1). The amplification of the RMP is consistent with the weak magnetic perturbation level (delta B/B approximate to 1. x 10(-4)) required for effective Type-I ELM suppression. Cyclic variations in the pedestal pressure, width, and toroidal rotation are observed in these plasmas, correlated with cyclic variations in the strength and frequency of the grassy-ELMs. Extended MHD analysis and magnetic measurements indicate that these pedestal pulsations are driven by cyclic variations in the resonant field strength at the top of the pedestal. These pedestal pulsations reveal that the grassy-ELMs is correlated with the proximity of the pedestal to the low-n peeling-ballooning mode stability boundary. The use of low amplitude magnetic fields to access grassy-ELM conditions free of Type-I ELMs in high beta poloidal plasmas (beta(p) approximate to 1.5-2.0) opens the possibility for the further optimization of the steady-state tokamak by use of edge resonant magnetic perturbations.