Showing papers in "Physics of Plasmas in 2016"

Three-dimensional simulations of low foot and high foot implosion experiments on the National Ignition Facility

Abstract: In order to achieve the several hundred Gbar stagnation pressures necessary for inertial confinement fusion ignition, implosion experiments on the National Ignition Facility (NIF) [E. I. Moses et al., Phys. Plasmas 16, 041006 (2009)] require the compression of deuterium-tritium fuel layers by a convergence ratio as high as forty. Such high convergence implosions are subject to degradation by a range of perturbations, including the growth of small-scale defects due to hydrodynamic instabilities, as well as longer scale modulations due to radiation flux asymmetries in the enclosing hohlraum. Due to the broad range of scales involved, and also the genuinely three-dimensional (3D) character of the flow, accurately modeling NIF implosions remains at the edge of current simulation capabilities. This paper describes the current state of progress of 3D capsule-only simulations of NIF implosions aimed at accurately describing the performance of specific NIF experiments. Current simulations include the effects of hohlraum radiation asymmetries, capsule surface defects, the capsule support tent and fill tube, and use a grid resolution shown to be converged in companion two-dimensional simulations. The results of detailed simulations of low foot implosions from the National Ignition Campaign are contrasted against results for more recent high foot implosions. While the simulations suggest that low foot performance was dominated by ablation front instability growth, especially the defect seeded by the capsule support tent, high foot implosions appear to be dominated by hohlraum flux asymmetries, although the support tent still plays a significant role. For both implosion types, the simulations show reasonable, though not perfect, agreement with the data and suggest that a reliable predictive capability is developing to guide future implosions toward ignition.

Radiation pressure acceleration of protons to 93 MeV with circularly polarized petawatt laser pulses

Abstract: The radiation pressure acceleration (RPA) of charged particles has been a challenging task in laser-driven proton/ion acceleration due to its stringent requirements in laser and target conditions. The realization of radiation-pressure-driven proton acceleration requires irradiating ultrathin targets with an ultrahigh contrast and ultraintense laser pulses. We report the generation of 93-MeV proton beams achieved by applying 800-nm 30-fs circularly polarized laser pulses with an intensity of 6.1×1020 W/cm2 to 15-nm-thick polymer targets. The radiation pressure acceleration was confirmed from the obtained optimal target thickness, quadratic energy scaling, polarization dependence, and three-dimensional particle-in-cell simulations. We expect this clear demonstration of RPA to facilitate the realization of laser-driven proton/ion sources delivering energetic and short-pulse particle beams for novel applications.

Laser absorption via quantum electrodynamics cascades in counter propagating laser pulses

Abstract: A model for laser light absorption in electron–positron plasmas self-consistently created via QED cascades is described. The laser energy is mainly absorbed due to hard photon emission via nonlinear Compton scattering. The degree of absorption depends on the laser intensity and the pulse duration. The QED cascades are studied with multi-dimensional particle-in-cell simulations complemented by a QED module and a macro-particle merging algorithm that allows to handle the exponential growth of the number of particles. Results range from moderate-intensity regimes ( ∼10 PW) where the laser absorption is negligible to extreme intensities ( >100 PW) where the degree of absorption reaches 80%. Our study demonstrates good agreement between the analytical model and simulations. The expected properties of the hard photon emission and the generated pair-plasma are investigated, and the experimental signatures for near-future laser facilities are discussed.

General theory of the plasmoid instability

Abstract: A general theory of the onset and development of the plasmoid instability is formulated by means of a principle of least time. The scaling relations for the final aspect ratio, transition time to rapid onset, growth rate, and number of plasmoids are derived and shown to depend on the initial perturbation amplitude (w0), the characteristic rate of current sheet evolution (1/τ), and the Lundquist number (S). They are not simple power laws, and are proportional to Sατβ[ln f(S,τ,w0)]σ. The detailed dynamics of the instability is also elucidated, and shown to comprise of a period of quiescence followed by sudden growth over a short time scale.

Divertor plasma detachment

Abstract: Regime with the plasma detached from the divertor targets (detached divertor regime) is a natural continuation of the high recycling conditions to higher density and stronger impurity radiation loss. Both the theoretical considerations and experimental data show clearly that the increase of the impurity radiation loss and volumetric plasma recombination causes the rollover of the plasma flux to the target when the density increases, which is the manifestation of detachment. Plasma-neutral friction (neutral viscosity effects), although important for the sustainment of high density/pressure plasma upstream and providing the conditions for efficient recombination and power loss, is not directly involved in the reduction of the plasma flux to the targets. The stability of detachment is also discussed.

Beyond the ponderomotive limit:Direct laser acceleration of relativistic electrons in sub-critical plasmas

Abstract: We examine a regime in which a linearly polarized laser pulse with relativistic intensity irradiates a sub-critical plasma for much longer than the characteristic electron response time. A steady-state channel is formed in the plasma in this case with quasi-static transverse and longitudinal electric fields. These relatively weak fields significantly alter the electron dynamics. The longitudinal electric field reduces the longitudinal dephasing between the electron and the wave, leading to an enhancement of the electron energy gain from the pulse. The energy gain in this regime is ultimately limited by the superluminosity of the wave fronts induced by the plasma in the channel. The transverse electric field alters the oscillations of the transverse electron velocity, allowing it to remain anti-parallel to laser electric field and leading to a significant energy gain. The energy enhancement is accompanied by the development of significant oscillations perpendicular to the plane of the driven motion, making trajectories of energetic electrons three-dimensional. Proper electron injection into the laser beam can further boost the electron energy gain.

Theory for the anomalous electron transport in Hall effect thrusters. I. Insights from particle-in-cell simulations

Abstract: Using a 1D particle-in-cell simulation with perpendicular electric, E0, and magnetic, B0, fields, and modelling the azimuthal direction (i.e., the E0 × B0 direction), we study the cross-field electron transport in Hall effect thrusters (HETs). For low plasma densities, the electron transport is found to be well described by classical electron-neutral collision theory, but at sufficiently high densities (representative of typical HETs), a strong instability is observed to significantly enhance the electron mobility, even in the absence of electron-neutral collisions. This instability is associated with correlated high-frequency (of the order of MHz) and short-wavelength (of the order of mm) fluctuations in both the electric field and the plasma density, which are shown to be the cause of the anomalous transport. Saturation of the instability is observed to occur due to a combination of ion-wave trapping in the E0 × B0 direction, and convection in the E0 direction.

Theory for the anomalous electron transport in Hall effect thrusters. II. Kinetic model

Abstract: In Paper I [T. Lafleur et al., Phys. Plasmas 23, 053502 (2016)], we demonstrated (using particle-in-cell simulations) the definite correlation between an anomalously high cross-field electron transport in Hall effect thrusters (HETs), and the presence of azimuthal electrostatic instabilities leading to enhanced electron scattering. Here, we present a kinetic theory that predicts the enhanced scattering rate and provides an electron cross-field mobility that is in good agreement with experiment. The large azimuthal electron drift velocity in HETs drives a strong instability that quickly saturates due to a combination of ion-wave trapping and wave-convection, leading to an enhanced mobility many orders of magnitude larger than that expected from classical diffusion theory. In addition to the magnetic field strength, B0, this enhanced mobility is a strong function of the plasma properties (such as the plasma density) and therefore does not, in general, follow simple 1/B02 or 1/B0 scaling laws.

Two decades of progress in understanding and control of laser plasma instabilities in indirect drive inertial fusion

Abstract: Our understanding of laser-plasma instability (LPI) physics has improved dramatically over the past two decades through advancements in experimental techniques, diagnostics, and theoretical and modeling approaches. We have progressed from single-beam experiments—ns pulses with ∼kJ energy incident on hundred-micron-scale target plasmas with ∼keV electron temperatures—to ones involving nearly 2 MJ energy in 192 beams onto multi-mm-scale plasmas with temperatures ∼4 keV. At the same time, we have also been able to use smaller-scale laser facilities to substantially improve our understanding of LPI physics and evaluate novel approaches to their control. These efforts have led to a change in paradigm for LPI research, ushering in an era of engineering LPI to accomplish specific objectives, from tuning capsule implosion symmetry to fixing nonlinear saturation of LPI processes at acceptable levels to enable the exploration of high energy density physics in novel plasma regimes. A tutorial is provided that reviews the progress in the field from the vantage of the foundational LPI experimental results. The pedagogical framework of the simplest models of LPI will be employed, but attention will also be paid to settings where more sophisticated models are needed to understand the observations. Prospects for the application of our improved understanding for inertial fusion (both indirect- and direct-drive) and other applications will also be discussed.

The role of zonal flows in the saturation of multi-scale gyrokinetic turbulence

Abstract: The 2D spectrum of the saturated electric potential from gyrokinetic turbulence simulations that include both ion and electron scales (multi-scale) in axisymmetric tokamak geometry is analyzed. The paradigm that the turbulence is saturated when the zonal (axisymmetic) ExB flow shearing rate competes with linear growth is shown to not apply to the electron scale turbulence. Instead, it is the mixing rate by the zonal ExB velocity spectrum with the turbulent distribution function that competes with linear growth. A model of this mechanism is shown to be able to capture the suppression of electron-scale turbulence by ion-scale turbulence and the threshold for the increase in electron scale turbulence when the ion-scale turbulence is reduced. The model computes the strength of the zonal flow velocity and the saturated potential spectrum from the linear growth rate spectrum. The model for the saturated electric potential spectrum is applied to a quasilinear transport model and shown to accurately reproduce the electron and ion energy fluxes of the non-linear gyrokinetic multi-scale simulations. The zonal flow mixing saturation model is also shown to reproduce the non-linear upshift in the critical temperature gradient caused by zonal flows in ion-scale gyrokinetic simulations.

A hybrid-drive nonisobaric-ignition scheme for inertial confinement fusion

Abstract: A new hybrid-drive (HD) nonisobaric ignition scheme of inertial confinement fusion (ICF) is proposed, in which a HD pressure to drive implosion dynamics increases via increasing density rather than temperature in the conventional indirect drive (ID) and direct drive (DD) approaches. In this HD (combination of ID and DD) scheme, an assembled target of a spherical hohlraum and a layered deuterium-tritium capsule inside is used. The ID lasers first drive the shock to perform a spherical symmetry implosion and produce a large-scale corona plasma. Then, the DD lasers, whose critical surface in ID corona plasma is far from the radiation ablation front, drive a supersonic electron thermal wave, which slows down to a high-pressure electron compression wave, like a snowplow, piling up the corona plasma into high density and forming a HD pressurized plateau with a large width. The HD pressure is several times the conventional ID and DD ablation pressure and launches an enhanced precursor shock and a continuous compression wave, which give rise to the HD capsule implosion dynamics in a large implosion velocity. The hydrodynamic instabilities at imploding capsule interfaces are suppressed, and the continuous HD compression wave provides main pdV work large enough to hotspot, resulting in the HD nonisobaric ignition. The ignition condition and target design based on this scheme are given theoretically and by numerical simulations. It shows that the novel scheme can significantly suppress implosion asymmetry and hydrodynamic instabilities of current isobaric hotspot ignition design, and a high-gain ICF is promising.

Parallel electric fields are inefficient drivers of energetic electrons in magnetic reconnection

Abstract: We present two-dimensional kinetic simulations, with a broad range of initial guide fields, which isolate the role of parallel electric fields ( E∥) in energetic electron production during collisionless magnetic reconnection. In the strong guide field regime, E∥ drives essentially all of the electron energy gains, yet fails to generate an energetic component. We suggest that this is due to the weak energy scaling of particle acceleration from E∥ compared to that of a Fermi-type mechanism responsible for energetic electron production in the weak guide-field regime. This result has important implications for energetic electron production in astrophysical systems and reconnection-driven dissipation in turbulence.

Intermittency, coherent structures and dissipation in plasma turbulence

Abstract: Collisionless dissipation in turbulentplasmas such as the solar wind and the solar corona has been an intensively studied subject recently, with new insights often emerging from numerical simulation. Here we report results from high resolution, fully kinetic simulations of plasma turbulence in both two (2D) and three (3D) dimensions, studying the relationship between intermittency and dissipation. The simulations show development of turbulent coherent structures, characterized by sheet-like current density structures spanning a range of scales. An approximate dissipation measure is employed, based on work done by the electromagnetic field in the local electron fluid frame. This surrogate dissipation measure is highly concentrated in small subvolumes in both 2D and 3D simulations. Fully kinetic simulations are also compared with magnetohydrodynamics(MHD) simulations in terms of coherent structures and dissipation. The interesting result emerges that the conditional averages of dissipation measure scale very similarly with normalized current densityJ in 2D and 3D particle-in-cell and in MHD. To the extent that the surrogate dissipation measure is accurate, this result implies that on average dissipation scales as ∼J2 in turbulent kinetic plasma. Multifractal intermittency is seen in the inertial range in both 2D and 3D, but at scales ∼ion inertial length, the scaling is closer to monofractal.

Three-dimensional modeling of direct-drive cryogenic implosions on OMEGA

Abstract: The effects of large-scale (with Legendre modes ≲10) laser-imposed nonuniformities in direct-drive cryogenic implosions on the OMEGA Laser System are investigated using three-dimensional hydrodynamic simulations performed using the newly developed code ASTER. Sources of these nonuniformities include an illumination pattern produced by 60 OMEGA laser beams, capsule offsets (∼10–20 μm), and imperfect pointing, power balance, and timing of the beams (with typical σrms∼10 μm, 10%, and 5 ps, respectively). Two implosion designs using 26-kJ triple-picket laser pulses were studied: a nominal design, in which an 874-μm-diameter capsule is illuminated by about the same-diameter beams, and a more hydrodynamically efficient “R75” design using a 900-μm-diameter capsule and beams of 75% of this diameter. Simulations show that nonuniformities caused by capsule offsets and beam imbalance have the largest effect on implosion performance. These nonuniformities lead to significant distortions of implosion cores, resulting ...

Initial measurements of two- and three-dimensional ordering, waves, and plasma filamentation in the Magnetized Dusty Plasma Experiment

Abstract: The Magnetized Dusty Plasma Experiment at Auburn University has been operational for over one year. In that time, a number of experiments have been performed at magnetic fields up to B = 2.5 T to explore the interaction between magnetized plasmas and charged, micron-sized dust particles. This paper reports on the initial results from studies of: (a) the formation of imposed, ordered structures, (b) the properties of dust wave waves in a rotating frame, and (c) the generation of plasma filaments.

Gyrokinetic neoclassical study of the bootstrap current in the tokamak edge pedestal with fully non-linear Coulomb collisions

Abstract: As a follow-up on the drift-kinetic study of the non-local bootstrap current in the steep edge pedestal of tokamak plasma by Koh et al. [Phys. Plasmas 19, 072505 (2012)], a gyrokinetic neoclassical study is performed with gyrokinetic ions and drift-kinetic electrons. Besides the gyrokinetic improvement of ion physics from the drift-kinetic treatment, a fully non-linear Fokker-Planck collision operator—that conserves mass, momentum, and energy—is used instead of Koh et al.'s linearized collision operator in consideration of the possibility that the ion distribution function is non-Maxwellian in the steep pedestal. An inaccuracy in Koh et al.'s result is found in the steep edge pedestal that originated from a small error in the collisional momentum conservation. The present study concludes that (1) the bootstrap current in the steep edge pedestal is generally smaller than what has been predicted from the small banana-width (local) approximation [e.g., Sauter et al., Phys. Plasmas 6, 2834 (1999) and Belli et al., Plasma Phys. Controlled Fusion 50, 095010 (2008)], (2) the plasma flow evaluated from the local approximation can significantly deviate from the non-local results, and (3) the bootstrap current in the edge pedestal, where the passing particle region is small, can be dominantly carried by the trapped particles in a broad trapped boundary layer. A new analytic formula based on numerous gyrokinetic simulations using various magnetic equilibria and plasma profiles with self-consistent Grad-Shafranov solutions is constructed.

Kinetic interpretation of resonance phenomena in low pressure capacitively coupled radio frequency plasmas

Abstract: Low pressure capacitive radio frequency (RF) plasmas are often described by equivalent circuit models based on fluid approaches that predict the self-excitation of resonances, e.g., high frequency oscillations of the total current in asymmetric discharges, but do not provide a kinetic interpretation of these effects. In fact, they leave important questions open: How is current continuity ensured in the presence of energetic electron beams generated by the expanding sheaths that lead to a local enhancement of the conduction current propagating through the bulk? How do the beam electrons interact with cold bulk electrons? What is the kinetic origin of resonance phenomena? Based on kinetic simulations, we find that the energetic beam electrons interact with cold bulk electrons (modulated on a timescale of the inverse local electron plasma frequency) via a time dependent electric field outside the sheaths. This electric field is caused by the electron beam itself, which leaves behind a positive space charge, that attracts cold bulk electrons towards the expanding sheath. The resulting displacement current ensures current continuity by locally compensating the enhancement of the conduction current. The backflow of cold electrons and their interaction with the nonlinear plasma sheath cause the generation of multiple electron beams during one phase of sheath expansion and contribute to a strongly non-sinusoidal RF current. These kinetic mechanisms are the basis for a fundamental understanding of the electron power absorption dynamics and resonance phenomena in such plasmas, which are found to occur in discharges of different symmetries including perfectly symmetric plasmas.

Multi-scale gyrokinetic simulations: Comparison with experiment and implications for predicting turbulence and transport

Abstract: To better understand the role of cross-scale coupling in experimental conditions, a series of multi-scale gyrokinetic simulations were performed on Alcator C-Mod, L-mode plasmas. These simulations, performed using all experimental inputs and realistic ion to electron mass ratio ((mi/me)1∕2 = 60.0), simultaneously capture turbulence at the ion ( kθρs∼O(1.0)) and electron-scales ( kθρe∼O(1.0)). Direct comparison with experimental heat fluxes and electron profile stiffness indicates that Electron Temperature Gradient (ETG) streamers and strong cross-scale turbulence coupling likely exist in both of the experimental conditions studied. The coupling between ion and electron-scales exists in the form of energy cascades, modification of zonal flow dynamics, and the effective shearing of ETG turbulence by long wavelength, Ion Temperature Gradient (ITG) turbulence. The tightly coupled nature of ITG and ETG turbulence in these realistic plasma conditions is shown to have significant implications for the interpretat...

Validation metrics for turbulent plasma transport

Abstract: Developing accurate models of plasma dynamics is essential for confident predictive modeling of current and future fusion devices. In modern computer science and engineering, formal verification and validation processes are used to assess model accuracy and establish confidence in the predictive capabilities of a given model. This paper provides an overview of the key guiding principles and best practices for the development of validation metrics, illustrated using examples from investigations of turbulent transport in magnetically confined plasmas. Particular emphasis is given to the importance of uncertainty quantification and its inclusion within the metrics, and the need for utilizing synthetic diagnostics to enable quantitatively meaningful comparisons between simulation and experiment. As a starting point, the structure of commonly used global transport model metrics and their limitations is reviewed. An alternate approach is then presented, which focuses upon comparisons of predicted local fluxes, ...

Theory of anomalous backscattering in second harmonic X-mode ECRH experiments

Abstract: A quantitative model explaining generation of the anomalous backscattering signal in the second harmonic X-mode electron cyclotron resonance heating (ECRH) experiments at TEXTOR tokamak as a secondary nonlinear process which accompanies a primary low-threshold parametric decay instability (PDI) leading to excitation of two—upper hybrid (UH)—plasmons trapped in plasma is developed. The primary absolute PDI enhancing the UH wave fluctuations from the thermal noise level is supposed to be saturated due to a cascade of secondary low-threshold decays of the daughter UH wave leading to excitation of the secondary UH waves down-shifted in frequency and the ion Bernstein wave. A set of equations describing the cascade is derived and solved numerically. The results of numerical modelling are shown to be in agreement with the analytical estimations of the growth rate of the initial and secondary parametric decays and the saturation level. The generation of backscattering signal is explained by coupling of the daughter UH waves. The fine details of the frequency spectrum of the anomalously reflected extraordinary wave and the absolute value of the observed backscattering signal in the second harmonic X-mode ECRH experiments at TEXTOR are reproduced.

Integrated modeling of cryogenic layered highfoot experiments at the NIF

Abstract: Integrated radiation hydrodynamic modeling in two dimensions, including the hohlraum and capsule, of layered cryogenic HighFoot Deuterium-Tritium (DT) implosions on the NIF successfully predicts important data trends. The model consists of a semi-empirical fit to low mode asymmetries and radiation drive multipliers to match shock trajectories, one dimensional inflight radiography, and time of peak neutron production. Application of the model across the HighFoot shot series, over a range of powers, laser energies, laser wavelengths, and target thicknesses predicts the neutron yield to within a factor of two for most shots. The Deuterium-Deuterium ion temperatures and the DT down scattered ratios, ratio of (10–12)/(13–15) MeV neutrons, roughly agree with data at peak fuel velocities <340 km/s and deviate at higher peak velocities, potentially due to flows and neutron scattering differences stemming from 3D or capsule support tent effects. These calculations show a significant amount alpha heating, 1–2.5× for shots where the experimental yield is within a factor of two, which has been achieved by increasing the fuel kinetic energy. This level of alpha heating is consistent with a dynamic hot spot model that is matched to experimental data and as determined from scaling of the yield with peak fuel velocity. These calculations also show that low mode asymmetries become more important as the fuel velocity is increased, and that improving these low mode asymmetries can result in an increase in the yield by a factor of several.

Scaling magnetized liner inertial fusion on Z and future pulsed-power accelerators

Abstract: The MagLIF (Magnetized Liner Inertial Fusion) concept [S. A. Slutz et al., Phys. Plasmas 17, 056303 (2010)] has demonstrated fusion–relevant plasma conditions [M. R. Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)] on the Z accelerator with a peak drive current of about 18 MA. We present 2D numerical simulations of the scaling of MagLIF on Z as a function of drive current, preheat energy, and applied magnetic field. The results indicate that deuterium-tritium (DT) fusion yields greater than 100 kJ could be possible on Z when all of these parameters are at the optimum values: i.e., peak current = 25 MA, deposited preheat energy = 5 kJ, and Bz = 30 T. Much higher yields have been predicted [S. A. Slutz and R. A. Vesey, Phys. Rev. Lett. 108, 025003 (2012)] for MagLIF driven with larger peak currents. Two high performance pulsed-power accelerators (Z300 and Z800) based on linear-transformer-driver technology have been designed [W. A. Stygar et al., Phys. Rev. ST Accel. Beams 18, 110401 (2015)]. The Z300 design would provide 48 MA to a MagLIF load, while Z800 would provide 65 MA. Parameterized Thevenin-equivalent circuits were used to drive a series of 1D and 2D numerical MagLIF simulations with currents ranging from what Z can deliver now to what could be achieved by these conceptual future pulsed-power accelerators. 2D simulations of simple MagLIF targets containing just gaseous DT have yields of 18 MJ for Z300 and 440 MJ for Z800. The 2D simulated yield for Z800 is increased to 7 GJ by adding a layer of frozen DT ice to the inside of the liner.

Integrated fusion simulation with self-consistent core-pedestal coupling

Abstract: Accurate prediction of fusion performance in present and future tokamaks requires taking into account the strong interplay between core transport, pedestal structure, current profile, and plasma equilibrium. An integrated modeling workflow capable of calculating the steady-state self-consistent solution to this strongly coupled problem has been developed. The workflow leverages state-of-the-art components for collisional and turbulent core transport, equilibrium and pedestal stability. Testing against a DIII-D discharge shows that the workflow is capable of robustly predicting the kinetic profiles (electron and ion temperature and electron density) from the axis to the separatrix in a good agreement with the experiments. An example application is presented, showing self-consistent optimization for the fusion performance of the 15 MA D-T ITER baseline scenario as functions of the pedestal density and ion effective charge Zeff.

Numerical study of the interaction of a helium atmospheric pressure plasma jet with a dielectric material

Abstract: This is a computational modeling study of a cold atmospheric pressure helium plasma jet impinging on a dielectric surface placed normal to the jet axis. This study provides insights into the propagation mechanism of the plasma jet, the electrical properties, and the total accumulated charge density at the dielectric surface. For the radial streamer propagation along the dielectric surface, Penning ionization and the electron impact ionization of helium atoms are the major ionization reactions in the streamer head, while Penning ionization is the only dominant contributor along the streamer body. In addition, the plasma bullet velocity along the dielectric surface is 10–100 times lower than that in the plasma column. Increasing tube radius or helium flow rate lowers air entrainment in the plasma jet, leading to a decrease of the radial electric field and the accumulated charge density at the dielectric surface. Furthermore, the tube radius has weaker influence on the plasma properties as tube radius increases. For a target dielectric with lower relative permittivity, a higher radial electric field penetrates into the material, and the surface ionization wave along the dielectric surface extends farther. Higher relative permittivity of the treated dielectric results in more charging at the dielectric surface and more electron density in the plasma column.

Discovery of stationary operation of quiescent H-mode plasmas with net-zero neutral beam injection torque and high energy confinement on DIII-D

Abstract: Recent experiments in DIII-D [J. L. Luxon et al., in Plasma Physics and Controlled Nuclear Fusion Research 1996 (International Atomic Energy Agency, Vienna, 1987), Vol. I, p. 159] have led to the discovery of a means of modifying edge turbulence to achieve stationary, high confinement operation without Edge Localized Mode (ELM) instabilities and with no net external torque input. Eliminating the ELM-induced heat bursts and controlling plasma stability at low rotation represent two of the great challenges for fusion energy. By exploiting edge turbulence in a novel manner, we achieved excellent tokamak performance, well above the H98y2 international tokamak energy confinement scaling (H98y2 = 1.25), thus meeting an additional confinement challenge that is usually difficult at low torque. The new regime is triggered in double null plasmas by ramping the injected torque to zero and then maintaining it there. This lowers E × B rotation shear in the plasma edge, allowing low-k, broadband, electromagnetic turbul...

Ultrafast proton radiography of the magnetic fields generated by a laser-driven coil current

Abstract: Magnetic fields generated by a current flowing through a U-shaped coil connecting two copper foils were measured using ultrafast proton radiography. Two ∼1.25 kJ, 1-ns laser pulses propagated through laser entrance holes in the front foil and were focused to the back foil with an intensity of ∼3 × 1016 W/cm2. The intense laser-solid interaction induced a high voltage between the copper foils and generated a large current in the connecting coil. The proton data show ∼40–50 T magnetic fields at the center of the coil ∼3–4 ns after laser irradiation. The experiments provide significant insight for future target designs that aim to develop a powerful source of external magnetic fields for various applications in high-energy-density science.

Fast ignition realization experiment with high-contrast kilo-joule peta-watt LFEX laser and strong external magnetic field

Abstract: A petawatt laser for fast ignition experiments (LFEX) laser system [N. Miyanaga et al., J. Phys. IV France 133, 81 (2006)], which is currently capable of delivering 2 kJ in a 1.5 ps pulse using 4 laser beams, has been constructed beside the GEKKO-XII laser facility for demonstrating efficient fast heating of a dense plasma up to the ignition temperature under the auspices of the Fast Ignition Realization EXperiment (FIREX) project [H. Azechi et al., Nucl. Fusion 49, 104024 (2009)]. In the FIREX experiment, a cone is attached to a spherical target containing a fuel to prevent a corona plasma from entering the path of the intense heating LFEX laser beams. The LFEX laser beams are focused at the tip of the cone to generate a relativistic electron beam (REB), which heats a dense fuel core generated by compression of a spherical deuterized plastic target induced by the GEKKO-XII laser beams. Recent studies indicate that the current heating efficiency is only 0.4%, and three requirements to achieve higher efficiency of the fast ignition (FI) scheme with the current GEKKO and LFEX systems have been identified: (i) reduction of the high energy tail of the REB; (ii) formation of a fuel core with high areal density using a limited number (twelve) of GEKKO-XII laser beams as well as a limited energy (4 kJ of 0.53-μm light in a 1.3 ns pulse); (iii) guiding and focusing of the REB to the fuel core. Laser–plasma interactions in a long-scale plasma generate electrons that are too energetic to efficiently heat the fuel core. Three actions were taken to meet the first requirement. First, the intensity contrast of the foot pulses to the main pulses of the LFEX was improved to >109. Second, a 5.5-mm-long cone was introduced to reduce pre-heating of the inner cone wall caused by illumination of the unconverted 1.053-μm light of implosion beam (GEKKO-XII). Third, the outside of the cone wall was coated with a 40-μm plastic layer to protect it from the pressure caused by imploding plasma. Following the above improvements, conversion of 13% of the LFEX laser energy to a low energy portion of the REB, whose slope temperature is 0.7 MeV, which is close to the ponderomotive scaling value, was achieved. To meet the second requirement, the compression of a solid spherical ball with a diameter of 200-μm to form a dense core with an areal density of ∼0.07 g/cm2 was induced by a laser-driven spherically converging shock wave. Converging shock compression is more hydrodynamically stable compared to shell implosion, while a hot spot cannot be generated with a solid ball target. Solid ball compression is preferable also for compressing an external magnetic field to collimate the REB to the fuel core, due to the relatively small magnetic Reynolds number of the shock compressed region. To meet the third requirement, we have generated a strong kilo-tesla magnetic field using a laser-driven capacitor-coil target. The strength and time history of the magnetic field were characterized with proton deflectometry and a B-dot probe. Guidance of the REB using a 0.6-kT field in a planar geometry has been demonstrated at the LULI 2000 laser facility. In a realistic FI scenario, a magnetic mirror is formed between the REB generation point and the fuel core. The effects of the strong magnetic field on not only REB transport but also plasma compression were studied using numerical simulations. According to the transport calculations, the heating efficiency can be improved from 0.4% to 4% by the GEKKO and LFEX laser system by meeting the three requirements described above. This efficiency is scalable to 10% of the heating efficiency by increasing the areal density of the fuel core.

Second order gyrokinetic theory for particle-in-cell codes

Abstract: The main idea of the gyrokinetic dynamical reduction consists in a systematical removal of the fast scale motion (the gyromotion) from the dynamics of the plasma, resulting in a considerable simplification and a significant gain of computational time. The gyrokinetic Maxwell–Vlasov equations are nowadays implemented in for modeling (both laboratory and astrophysical) strongly magnetized plasmas. Different versions of the reduced set of equations exist, depending on the construction of the gyrokinetic reduction procedure and the approximations performed in the derivation. The purpose of this article is to explicitly show the connection between the general second order gyrokinetic Maxwell–Vlasov system issued from the modern gyrokinetic theory and the model currently implemented in the global electromagnetic Particle-in-Cell code ORB5. Necessary information about the modern gyrokinetic formalism is given together with the consistent derivation of the gyrokinetic Maxwell–Vlasov equations from first principles. The variational formulation of the dynamics is used to obtain the corresponding energy conservation law, which in turn is used for the verification of energy conservation diagnostics currently implemented in ORB5. This work fits within the context of the code verification project VeriGyro currently run at IPP Max-Planck Institut in collaboration with others European institutions.

Stochastic modelling of intermittent fluctuations in the scrape-off layer: Correlations, distributions, level crossings, and moment estimation

Abstract: A stochastic model is presented for intermittent fluctuations in the scrape-off layer of magnetically confined plasmas. The fluctuations in the plasma density are modeled by a super-position of uncorrelated pulses with fixed shape and duration, describing radial motion of blob-like structures. In the case of an exponential pulse shape and exponentially distributed pulse amplitudes, predictions are given for the lowest order moments, probability density function, auto-correlation function, level crossings, and average times for periods spent above and below a given threshold level. Also, the mean squared errors on estimators of sample mean and variance for realizations of the process by finite time series are obtained. These results are discussed in the context of single-point measurements of fluctuations in the scrape-off layer, broad density profiles, and implications for plasma–wall interactions due to the transient transport events in fusion grade plasmas. The results may also have wide applications for modelling fluctuations in other magnetized plasmas such as basic laboratory experiments and ionospheric irregularities.

Extended application of Kohn-Sham first-principles molecular dynamics method with plane wave approximation at high energy—From cold materials to hot dense plasmas

Abstract: An extended first-principles molecular dynamics (FPMD) method based on Kohn-Sham scheme is proposed to elevate the temperature limit of the FPMD method in the calculation of dense plasmas. The extended method treats the wave functions of high energy electrons as plane waves analytically and thus expands the application of the FPMD method to the region of hot dense plasmas without suffering from the formidable computational costs. In addition, the extended method inherits the high accuracy of the Kohn-Sham scheme and keeps the information of electronic structures. This gives an edge to the extended method in the calculation of mixtures of plasmas composed of heterogeneous ions, high-Z dense plasmas, lowering of ionization potentials, X-ray absorption/emission spectra, and opacities, which are of particular interest to astrophysics, inertial confinement fusion engineering, and laboratory astrophysics.