Showing papers in "Physics of Plasmas in 2014"

Review of the National Ignition Campaign 2009-2012

Abstract: The National Ignition Campaign (NIC) was a multi-institution effort established under the National Nuclear Security Administration of DOE in 2005, prior to the completion of the National Ignition Facility (NIF) in 2009. The scope of the NIC was the planning and preparation for and the execution of the first 3 yr of ignition experiments (through the end of September 2012) as well as the development, fielding, qualification, and integration of the wide range of capabilities required for ignition. Besides the operation and optimization of the use of NIF, these capabilities included over 50 optical, x-ray, and nuclear diagnostic systems, target fabrication facilities, experimental platforms, and a wide range of NIF facility infrastructure. The goal of ignition experiments on the NIF is to achieve, for the first time, ignition and thermonuclear burn in the laboratory via inertial confinement fusion and to develop a platform for ignition and high energy density applications on the NIF. The goal of the NIC was to develop and integrate all of the capabilities required for a precision ignition campaign and, if possible, to demonstrate ignition and gain by the end of FY12. The goal of achieving ignition can be divided into three main challenges. The first challenge is defining specifications for the target, laser, and diagnostics with the understanding that not all ignition physics is fully understood and not all material properties are known. The second challenge is designing experiments to systematically remove these uncertainties. The third challenge is translating these experimental results into metrics designed to determine how well the experimental implosions have performed relative to expectations and requirements and to advance those metrics toward the conditions required for ignition. This paper summarizes the approach taken to address these challenges, along with the progress achieved to date and the challenges that remain. At project completion in 2009, NIF lacked almost all the diagnostics and infrastructure required for ignition experiments. About half of the 3 yr period covered in this review was taken up by the effort required to install and performance qualify the equipment and experimental platforms needed for ignition experiments. Ignition on the NIF is a grand challenge undertaking and the results presented here represent a snapshot in time on the path toward that goal. The path forward presented at the end of this review summarizes plans for the Ignition Campaign on the NIF, which were adopted at the end of 2012, as well as some of the key results obtained since the end of the NIC.

Low temperature plasma biomedicine: A tutorial reviewa)

Abstract: Gas discharge plasmas formed at atmospheric pressure and near room temperature have recently been shown to be potentially useful for surface and wound sterilization, antisepsis, bleeding cessation, wound healing, and cancer treatment, among other biomedical applications. This tutorial review summarizes the field, stressing the likely role of reactive oxygen and nitrogen species created in these plasmas as the biologically and therapeutically active agents. Reactive species, including radicals and non-radical compounds, are generated naturally within the body and are now understood to be essential for normal biological functions. These species are known to be active agents in existing therapies for wound healing, infection control, and cancer treatment. But they are also observed at elevated levels in persons with many diseases and are associated with aging. The physical and chemical complexity of plasma medical devices and their associated biochemical effects makes the development of safe, effective plasma medical devices and procedures a challenge, but encouragingly rapid progress has been reported around the world in the last several years.

The link between shocks, turbulence, and magnetic reconnection in collisionless plasmas

Abstract: Global hybrid (electron fluid, kinetic ions) and fully kinetic simulations of the magnetosphere have been used to show surprising interconnection between shocks, turbulence, and magnetic reconnection. In particular, collisionless shocks with their reflected ions that can get upstream before retransmission can generate previously unforeseen phenomena in the post shocked flows: (i) formation of reconnecting current sheets and magnetic islands with sizes up to tens of ion inertial length. (ii) Generation of large scale low frequency electromagnetic waves that are compressed and amplified as they cross the shock. These “wavefronts” maintain their integrity for tens of ion cyclotron times but eventually disrupt and dissipate their energy. (iii) Rippling of the shock front, which can in turn lead to formation of fast collimated jets extending to hundreds of ion inertial lengths downstream of the shock. The jets, which have high dynamical pressure, “stir” the downstream region, creating large scale disturbances such as vortices, sunward flows, and can trigger flux ropes along the magnetopause. This phenomenology closes the loop between shocks, turbulence, and magnetic reconnection in ways previously unrealized. These interconnections appear generic for the collisionless plasmas typical of space and are expected even at planar shocks, although they will also occur at curved shocks as occur at planets or around ejecta.

Edge-localized-modes in tokamaksa)

Abstract: Edge-localized-modes (ELMs) are a ubiquitous feature of H-mode in tokamaks. When gradients in the H-mode transport barrier grow to exceed the MHD stability limit the ELM instability grows explosively, rapidly transporting energy and particles onto open field lines and material surfaces. Though ELMs provide additional particle and impurity transport through the H-mode transport barrier, enabling steady operation, the resulting heat flux transients to plasma facing surfaces project to large amplitude in future low collisionality burning plasma tokamaks. Measurements of the ELM heat flux deposition onto material surfaces in the divertor and main chamber indicate significant broadening compared to inter-ELM heat flux, with a timescale for energy deposition that is consistent with sonic ion flow and numerical simulation. Comprehensive ELM simulation is highlighting the important physics processes of ELM transport including parallel transport due to magnetic reconnection and turbulence resulting from collapse of the H-mode transport barrier. Encouraging prospects for ELM control and/or suppression in future tokamaks include intrinsic modes of ELM free operation, ELM triggering with frequent small pellet injection and the application of 3D magnetic fields.

The mechanisms of electron heating and acceleration during magnetic reconnection

Abstract: The heating of electrons in collisionless magnetic reconnection is explored in particle-in-cell simulations with non-zero guide fields so that electrons remain magnetized. In this regime, electric fields parallel to B accelerate particles directly, while those perpendicular to B do so through gradient-B and curvature drifts. The curvature drift drives parallel heating through Fermi reflection, while the gradient B drift changes the perpendicular energy through betatron acceleration. We present simulations in which we evaluate each of these mechanisms in space and time in order to quantify their role in electron heating. For a case with a small guide field (20% of the magnitude of the reconnecting component), the curvature drift is the dominant source of electron heating. However, for a larger guide field (equal to the magnitude of the reconnecting component) electron acceleration by the curvature drift is comparable to that of the parallel electric field. In both cases, the heating by the gradient B drift is negligible in magnitude. It produces net cooling because the conservation of the magnetic moment and the drop of B during reconnection produce a decrease in the perpendicular electron energy. Heating by the curvature drift dominates in the outflow exhausts where bent field lines expand to relax their tension and is therefore distributed over a large area. In contrast, the parallel electric field is localized near X-lines. This suggests that acceleration by parallel electric fields may play a smaller role in large systems where the X-line occupies a vanishing fraction of the system. The curvature drift and the parallel electric field dominate the dynamics and drive parallel heating. A consequence is that the electron energy spectrum becomes extremely anisotropic at late time, which has important implications for quantifying the limits of electron acceleration due to synchrotron emission. An upper limit on electron energy gain that is substantially higher than earlier estimates is obtained by balancing reconnection drive with radiative loss.

The high-foot implosion campaign on the National Ignition Facilitya)

Abstract: The “High-Foot” platform manipulates the laser pulse-shape coming from the National Ignition Facility laser to create an indirect drive 3-shock implosion that is significantly more robust against instability growth involving the ablator and also modestly reduces implosion convergence ratio. This strategy gives up on theoretical high-gain in an inertial confinement fusion implosion in order to obtain better control of the implosion and bring experimental performance in-line with calculated performance, yet keeps the absolute capsule performance relatively high. In this paper, we will cover the various experimental and theoretical motivations for the high-foot drive as well as cover the experimental results that have come out of the high-foot experimental campaign. At the time of this writing, the high-foot implosion has demonstrated record total deuterium-tritium yields (9.3×1015) with low levels of inferred mix, excellent agreement with implosion simulations, fuel energy gains exceeding unity, and evidence for the “bootstrapping” associated with alpha-particle self-heating.

Stability analysis for two-dimensional ion-acoustic waves in quantum plasmas

Abstract: The quantum hydrodynamic model is applied to two-dimensional ion-acoustic waves in quantum plasmas. The two-dimensional quantum hydrodynamic model is used to obtain a deformed Kortewegde Vries (dKdV) equation by reductive perturbation method. By using the solution of auxiliary ordinary equations, a extended direct algebraic method is described to construct the exact solutions for nonlinear quantum dKdV equation. The present results are describing the generation and evolution of such waves, their interactions, and their stability.

Metrics for long wavelength asymmetries in inertial confinement fusion implosions on the National Ignition Facility

Abstract: We investigate yield degradation due to applied low mode P2 and P4 asymmetries in layered inertial confinement fusion implosions. This study has been performed with a large database of >600 2D simulations. We show that low mode radiation induced drive asymmetries can result in significant deviation between the core hot spot shape and the fuel ρR shape at peak compression. In addition, we show that significant residual kinetic energy at peak compression can be induced by these low mode asymmetries. We have developed a metric, which is a function of the hot spot shape, fuel ρR shape, and residual kinetic energy at peak compression, that is well correlated to yield degradation due to low mode shape perturbations. It is shown that the ρR shape and residual kinetic energy cannot, in general, be recovered by inducing counter asymmetries to make the hot core emission symmetric. In addition, we show that the yield degradation due to low mode asymmetries is well correlated to measurements of time dependent shape throughout the entire implosion, including early time shock symmetry and inflight fuel symmetry.

Computing the reconnection rate in turbulent kinetic layers by using electron mixing to identify topology

Abstract: Three-dimensional kinetic simulations of magnetic reconnection for parameter regimes relevant to the magnetopause current layer feature the development of turbulence, driven by the magnetic and velocity shear, and dominated by coherent structures including flux ropes, current sheets, and flow vortices. Here, we propose a new approach for computing the global reconnection rate in the presence of this complexity. The mixing of electrons originating from separate sides of the magnetopause layer is used as a proxy to rapidly identify the magnetic topology and track the evolution of magnetic flux. The details of this method are illustrated for an asymmetric current layer relevant to the subsolar magnetopause and for a flow shear dominated layer relevant to the lower latitude magnetopause. While the three-dimensional reconnection rates show a number of interesting differences relative to the corresponding two-dimensional simulations, the time scale for the energy conversion remains very similar. These results suggest that the mixing of field lines between topologies is more easily influenced by kinetic turbulence than the physics responsible for the energy conversion.

Improving the hot-spot pressure and demonstrating ignition hydrodynamic equivalence in cryogenic deuterium–tritium implosions on OMEGAa)

Abstract: Reaching ignition in direct-drive (DD) inertial confinement fusion implosions requires achieving central pressures in excess of 100 Gbar. The OMEGA laser system [T. R. Boehly et al., Opt. Commun. 133, 495 (1997)] is used to study the physics of implosions that are hydrodynamically equivalent to the ignition designs on the National Ignition Facility (NIF) [J. A. Paisner et al., Laser Focus World 30, 75 (1994)]. It is shown that the highest hot-spot pressures (up to 40 Gbar) are achieved in target designs with a fuel adiabat of α ≃ 4, an implosion velocity of 3.8 × 107 cm/s, and a laser intensity of ∼1015 W/cm2. These moderate-adiabat implosions are well understood using two-dimensional hydrocode simulations. The performance of lower-adiabat implosions is significantly degraded relative to code predictions, a common feature between DD implosions on OMEGA and indirect-drive cryogenic implosions on the NIF. Simplified theoretical models are developed to gain physical understanding of the implosion dynamics that dictate the target performance. These models indicate that degradations in the shell density and integrity (caused by hydrodynamic instabilities during the target acceleration) coupled with hydrodynamics at stagnation are the main failure mechanisms in low-adiabat designs. To demonstrate ignition hydrodynamic equivalence in cryogenic implosions on OMEGA, the target-design robustness to hydrodynamic instability growth must be improved by reducing laser-coupling losses caused by cross beam energy transfer.

Design of magnetized liner inertial fusion experiments using the Z facilitya)

Abstract: The magnetized liner inertial fusion concept has been presented as a path toward obtaining substantial thermonuclear fusion yields using the Z accelerator [S. A. Slutz et al., Phys. Plasmas 17, 056303 (2010)]. We present the first integrated magnetohydrodynamic simulations of the inertial fusion targets, which self-consistently include laser preheating of the fuel, the presence of electrodes, and end loss effects. These numerical simulations provided the design for the first thermonuclear fusion neutron-producing experiments on Z using capabilities that presently exist: peak currents of Imax = 18–20 MA, pre-seeded axial magnetic fields of Bz0=10 T, laser preheat energies of about Elas = 2 kJ delivered in 2 ns, DD fuel, and an aspect ratio 6 solid Be liner imploded to 70 km/s. Specific design details and observables for both near-term and future experiments are discussed, including sensitivity to laser timing and absorbed preheat energy. The initial experiments measured stagnation radii rstag<75 μm, temper...

Fluid description of multi-component solar partially ionized plasma

Abstract: We derive self-consistent formalism for the description of multi-component partially ionized solar plasma, by means of the coupled equations for the charged and neutral components for an arbitrary number of chemical species, and the radiation field. All approximations and assumptions are carefully considered. Generalized Ohm's law is derived for the single-fluid and two-fluid formalism. Our approach is analytical with some order-of-magnitude support calculations. After general equations are developed, we particularize to some frequently considered cases as for the interaction of matter and radiation.

20 years of research on the Alcator C-Mod tokamak

Abstract: The object of this review is to summarize the achievements of research on the Alcator C-Mod tokamak [Hutchinson et al., Phys. Plasmas 1, 1511 (1994) and Marmar, Fusion Sci. Technol. 51, 261 (2007)] and to place that research in the context of the quest for practical fusion energy. C-Mod is a compact, high-field tokamak, whose unique design and operating parameters have produced a wealth of new and important results since it began operation in 1993, contributing data that extends tests of critical physical models into new parameter ranges and into new regimes. Using only high-power radio frequency (RF) waves for heating and current drive with innovative launching structures, C-Mod operates routinely at reactor level power densities and achieves plasma pressures higher than any other toroidal confinement device. C-Mod spearheaded the development of the vertical-target divertor and has always operated with high-Z metal plasma facing components—approaches subsequently adopted for ITER. C-Mod has made ground-breaking discoveries in divertor physics and plasma-material interactions at reactor-like power and particle fluxes and elucidated the critical role of cross-field transport in divertor operation, edge flows and the tokamak density limit. C-Mod developed the I-mode and the Enhanced Dα H-mode regimes, which have high performance without large edge localized modes and with pedestal transport self-regulated by short-wavelength electromagnetic waves. C-Mod has carried out pioneering studies of intrinsic rotation and demonstrated that self-generated flow shear can be strong enough in some cases to significantly modify transport. C-Mod made the first quantitative link between the pedestal temperature and the H-mode's performance, showing that the observed self-similar temperature profiles were consistent with critical-gradient-length theories and followed up with quantitative tests of nonlinear gyrokinetic models. RF research highlights include direct experimental observation of ion cyclotron range of frequency (ICRF) mode-conversion, ICRF flow drive, demonstration of lower-hybrid current drive at ITER-like densities and fields and, using a set of novel diagnostics, extensive validation of advanced RF codes. Disruption studies on C-Mod provided the first observation of non-axisymmetric halo currents and non-axisymmetric radiation in mitigated disruptions. A summary of important achievements and discoveries are included.

High-density carbon ablator experiments on the National Ignition Facilitya)

Abstract: High Density Carbon (HDC) is a leading candidate as an ablator material for Inertial Confinement Fusion (ICF) capsules in x-ray (indirect) drive implosions. HDC has a higher density (3.5 g/cc) than plastic (CH, 1 g/cc), which results in a thinner ablator with a larger inner radius for a given capsule scale. This leads to higher x-ray absorption and shorter laser pulses compared to equivalent CH designs. This paper will describe a series of experiments carried out to examine the feasibility of using HDC as an ablator using both gas filled hohlraums and lower density, near vacuum hohlraums. These experiments have shown that deuterium (DD) and deuterium-tritium gas filled HDC capsules driven by a hohlraum filled with 1.2 mg/cc He gas, produce neutron yields a factor of 2× higher than equivalent CH implosions, representing better than 50% Yield-over-Clean (YoC). In a near vacuum hohlraum (He = 0.03 mg/cc) with 98% laser-to-hohlraum coupling, such a DD gas-filled capsule performed near 1D expectations. A cryogenic layered implosion version was consistent with a fuel velocity = 410 ± 20 km/s with no observed ablator mixing into the hot spot.

An in-flight radiography platform to measure hydrodynamic instability growth in inertial confinement fusion capsules at the National Ignition Facility

Abstract: A new in-flight radiography platform has been established at the National Ignition Facility (NIF) to measure Rayleigh–Taylor and Richtmyer–Meshkov instability growth in inertial confinement fusion capsules. The platform has been tested up to a convergence ratio of 4. An experimental campaign is underway to measure the growth of pre-imposed sinusoidal modulations of the capsule surface, as a function of wavelength, for a pair of ignition-relevant laser drives: a “low-foot” drive representative of what was fielded during the National Ignition Campaign (NIC) [Edwards et al., Phys. Plasmas 20, 070501 (2013)] and the new high-foot [Dittrich et al., Phys. Rev. Lett. 112, 055002 (2014); Park et al., Phys. Rev. Lett. 112, 055001 (2014)] pulse shape, for which the predicted instability growth is much lower. We present measurements of Legendre modes 30, 60, and 90 for the NIC-type, low-foot, drive, and modes 60 and 90 for the high-foot drive. The measured growth is consistent with model predictions, including much less growth for the high-foot drive, demonstrating the instability mitigation aspect of this new pulse shape. We present the design of the platform in detail and discuss the implications of the data it generates for the on-going ignition effort at NIF.

The effect of turbulent kinetic energy on inferred ion temperature from neutron spectra

Abstract: Measuring the width of the energy spectrum of fusion-produced neutrons from deuterium (DD) or deuterium-tritium (DT) plasmas is a commonly used method for determining the ion temperature in inertial confinement fusion (ICF) implosions. In a plasma with a Maxwellian distribution of ion energies, the spread in neutron energy arises from the thermal spread in the center-of-mass velocities of reacting pairs of ions. Fluid velocities in ICF are of a similar magnitude as the center-of-mass velocities and can lead to further broadening of the neutron spectrum, leading to erroneous inference of ion temperature. Motion of the reacting plasma will affect DD and DT neutrons differently, leading to disagreement between ion temperatures inferred from the two reactions. This effect may be a contributor to observations over the past decades of ion temperatures higher than expected from simulations, ion temperatures in disagreement with observed yields, and different temperatures measured in the same implosion from DD and DT neutrons. This difference in broadening of DD and DT neutrons also provides a measure of turbulent motion in a fusion plasma.

Theory on excitations of drift Alfvén waves by energetic particles. II. The general fishbone-like dispersion relation

Abstract: The theoretical framework of the general fishbone-like dispersion relation (GFLDR), presented and discussed in the Companion Paper [Phys. Plasmas 21, 072120 (2014)], is applied to cases of practical interest of shear/drift Alfven waves (SAWs/DAWs) excited by energetic particles (EPs) in toroidal fusion plasmas. These applications demonstrate that the GFLDR provides a unified approach that allows analytical and numerical calculations of stability properties, as well as mode structures and, in general, nonlinear evolutions, based on different models and with different levels of approximation. They also show the crucial importance of kinetic descriptions, accurate geometries and boundary conditions for predicting linear as well as nonlinear SAW/DAW and EP behaviors in burning plasmas. Thus, the GFLDR unified theoretical framework elevates the interpretative capability for both experimental and numerical simulation results.

Three dimensional simulations of plasma filaments in the scrape off layer: A comparison with models of reduced dimensionality

Abstract: This paper presents simulations of isolated 3D filaments in a slab geometry obtained using a newly developed 3D reduced fluid code, written using the BOUT++ framework. First, systematic scans were performed to investigate how the dynamics of a filament are affected by its amplitude, perpendicular size, and parallel extent. The perpendicular size of the filament was found to have a strong influence on its motions, as it determined the relative importance of parallel currents to polarization and viscous currents, whilst drift-wave instabilities were observed if the initial amplitude of the blob was increased sufficiently. Next, the 3D simulations were compared to 2D simulations using different parallel closures; namely, the sheath dissipation closure, which neglects parallel gradients, and the vorticity advection closure, which neglects the influence of parallel currents. The vorticity advection closure was found to not replicate the 3D perpendicular dynamics and overestimated the initial radial acceleration of all the filaments studied. In contrast, a more satisfactory comparison with the sheath dissipation closure was obtained, even in the presence of significant parallel gradients, where the closure is no longer valid. Specifically, it captured the contrasting dynamics of filaments with different perpendicular sizes that were observed in the 3D simulations which the vorticity advection closure failed to replicate. However, neither closure successfully replicated the Boltzmann spinning effects nor the associated poloidal drift of the blob that was observed in the 3D simulations. Although the sheath dissipation closure was concluded to be more successful in replicating the 3D dynamics, it is emphasized that the vorticity advection closure may still be relevant for situations where the parallel current is inhibited from closing through the sheath due to effects such as strong magnetic shear around X points or increased resistivity near the targets.

Electron heating during magnetic reconnection: A simulation scaling study

Abstract: Electron bulk heating during magnetic reconnection with symmetric inflow conditions is examined using kinetic particle-in-cell simulations. Inflowing plasma parameters are varied over a wide range of conditions, and the increase in electron temperature is measured in the exhaust well downstream of the x-line. The degree of electron heating is well correlated with the inflowing Alfven speed cAr based on the reconnecting magnetic field through the relation ΔTe=0.033 mi cAr2, where ΔTe is the increase in electron temperature. For the range of simulations performed, the heating shows almost no correlation with inflow total temperature Ttot=Ti+Te or plasma β. An out-of-plane (guide) magnetic field of similar magnitude to the reconnecting field does not affect the total heating, but it does quench perpendicular heating, with almost all heating being in the parallel direction. These results are qualitatively consistent with a recent statistical survey of electron heating in the dayside magnetopause (Phan et al.,...

Variational formulation of macro-particle plasma simulation algorithmsa)

Abstract: A variation formulation of macro-particle kinetic plasma models is discussed. In the electrostatic case, the use of symplectic integrators is investigated and found to offer advantages over typical generic methods. For the electromagnetic case, gauge invariance and momentum conservation are considered in detail. It is shown that, while the symmetries responsible for these conservation laws are broken in the presence of a spatial grid, the conservation laws hold in an average sense. The requirements for exact invariance are explored and it is shown that one viable option is to represent the potentials with a truncated Fourier basis.

An ITPA joint experiment to study runaway electron generation and suppression

Abstract: Recent results from an ITPA joint experiment to study the onset, growth, and decay of relativistic electrons (REs) indicate that loss mechanisms other than collisional damping may play a dominant role in the dynamics of the RE population, even during the quiescent Ip flattop. Understanding the physics of RE growth and mitigation is motivated by the theoretical prediction that disruptions of full-current (15 MA) ITER discharges could generate up to 10 MA of REs with 10–20 MeV energies. The ITPA MHD group is conducting a joint experiment to measure the RE detection threshold conditions on a number of tokamaks under quasi-steady-state conditions in which Vloop, ne, and REs can be well-diagnosed and compared to collisional theory. Data from DIII-D, C-Mod, FTU, KSTAR, and TEXTOR have been obtained so far, and the consensus to date is that the threshold E-field is significantly higher than predicted by relativistic collisional theory, or conversely, the density required to damp REs is significantly less than predicted, which could have significant implications for RE mitigation on ITER.

Mode 1 drive asymmetry in inertial confinement fusion implosions on the National Ignition Facility

Abstract: Mode 1 radiation drive asymmetry (pole-to-pole imbalance) at significant levels can have a large impact on inertial confinement fusion implosions at the National Ignition Facility (NIF). This asymmetry distorts the cold confining shell and drives a high-speed jet through the hot spot. The perturbed hot spot shows increased residual kinetic energy and reduced internal energy, and it achieves reduced pressure and neutron yield. The altered implosion physics manifests itself in observable diagnostic signatures, especially the neutron spectrum which can be used to measure the neutron-weighted flow velocity, apparent ion temperature, and neutron downscattering. Numerical simulations of implosions with mode 1 asymmetry show that the resultant simulated diagnostic signatures are moved toward the values observed in many NIF experiments. The diagnostic output can also be used to build a set of integrated implosion performance metrics. The metrics indicate that P1 has a significant impact on implosion performance a...

Dynamic symmetry of indirectly driven inertial confinement fusion capsules on the National Ignition Facilitya)

Abstract: In order to achieve ignition using inertial confinement fusion it is important to control the growth of low-mode asymmetries as the capsule is compressed. Understanding the time-dependent evolution of the shape of the hot spot and surrounding fuel layer is crucial to optimizing implosion performance. A design and experimental campaign to examine sources of asymmetry and to quantify symmetry throughout the implosion has been developed and executed on the National Ignition Facility (NIF) [E. I. Moses et al., Phys. Plasmas 16, 041006 (2009)]. We have constructed a large simulation database of asymmetries applied during different time intervals. Analysis of the database has shown the need to measure and control the hot-spot shape, areal density distribution, and symmetry swings during the implosion. The shape of the hot spot during final stagnation is measured using time-resolved imaging of the self-emission, and information on the shape of the fuel at stagnation can be obtained from Compton radiography [R. Tommasini et al., Phys. Plasmas 18, 056309 (2011)]. For the first time on NIF, two-dimensional inflight radiographs of gas-filled and cryogenic fuel layered capsules have been measured to infer the symmetry of the radiation drive on the capsule. These results have been used to modify the hohlraum geometry and the wavelength tuning to improve the inflight implosion symmetry. We have also expanded our shock timing capabilities by the addition of extra mirrors inside the re-entrant cone to allow the simultaneous measurement of shock symmetry in three locations on a single shot, providing asymmetry information up to Legendre mode 4. By diagnosing the shape at nearly every step of the implosion, we estimate that shape has typically reduced fusion yield by about 50% in ignition experiments.

Comparison of particle trajectories and collision operators for collisional transport in nonaxisymmetric plasmas

Abstract: In this work, we examine the validity of several common simplifying assumptions used in numerical neoclassical calculations for nonaxisymmetric plasmas, both by using a new continuum drift-kinetic code and by considering analytic properties of the kinetic equation. First, neoclassical phenomena are computed for the LHD and W7-X stellarators using several versions of the drift-kinetic equation, including the commonly used incompressible-E × B-drift approximation and two other variants, corresponding to different effective particle trajectories. It is found that for electric fields below roughly one third of the resonant value, the different formulations give nearly identical results, demonstrating the incompressible E × B-drift approximation is quite accurate in this regime. However, near the electric field resonance, the models yield substantially different results. We also compare results for various collision operators, including the full linearized Fokker-Planck operator. At low collisionality, the radial transport driven by radial gradients is nearly identical for the different operators; while in other cases, it is found to be important that collisions conserve momentum.

Gamma-ray flares in the Crab Nebula: A case of relativistic reconnection?

Abstract: The Crab Nebula was formed after the collapse of a massive star about a thousand years ago, leaving behind a pulsar that inflates a bubble of ultra-relativistic electron-positron pairs permeated with magnetic field. The observation of brief but bright flares of energetic gamma rays suggests that pairs are accelerated to PeV energies within a few days; such rapid acceleration cannot be driven by shocks. Here, it is argued that the flares may be the smoking gun of magnetic dissipation in the Nebula. Using 2D and 3D particle-in-cell simulations, it is shown that the observations are consistent with relativistic magnetic reconnection, where pairs are subject to strong radiative cooling. The Crab flares may highlight the importance of relativistic magnetic reconnection in astrophysical sources.

Energy partition, γ-ray emission, and radiation reaction in the near-quantum electrodynamical regime of laser-plasma interaction

Abstract: When extremely intense lasers (I ≥ 1022 W/cm2) interact with plasmas, a significant fraction of the pulse energy is converted into photon emission in the multi-MeV energy range. This emission results in a radiation reaction (RR) force on electrons, which becomes important at ultrahigh intensities. Using three-dimensional particle-in-cell simulations which include a quantum electrodynamics model for the γ–photons emission, the corresponding RR force and electron-positron pair creation, the energy partition in the laser-plasma system is investigated. At sufficiently high laser amplitudes, the fraction of laser energy coupled to electrons decreases, while the energy converted to γ-photons increases. The interaction becomes an efficient source of γ-rays when I > 1024 W/cm2, with up to 40% of the laser energy converted to high-energy photons. A systematic study of energy partition and γ-photon emission angle shows the influence of laser intensity and polarization for two plasma conditions: high-density carbon ...

Hybrid Vlasov-Maxwell simulations of two-dimensional turbulence in plasmas

Abstract: Turbulence in plasmas is a very challenging problem since it involves wave-particle interactions, which are responsible for phenomena such as plasma dissipation, acceleration mechanisms, heating, temperature anisotropy, and so on. In this work, a hybrid Vlasov-Maxwell numerical code is employed to study local kinetic processes in a two-dimensional turbulent regime. In the present model, ions are treated as a kinetic species, while electrons are considered as a fluid. As recently reported in [S. Servidio, Phys. Rev. Lett. 108, 045001 (2012)], nearby regions of strong magnetic activity, kinetic effects manifest through a deformation of the ion velocity distribution function that consequently departs from the equilibrium Maxwellian configuration. Here, the structure of turbulence is investigated in detail in phase space, by evaluating the high-order moments of the particle velocity distribution, i.e., temperature, skewness, and kurtosis. This analysis provides quantitative information about the non-Maxwellian character of the system dynamics. This departure from local thermodynamic equilibrium triggers several processes commonly observed in many astrophysical and laboratory plasmas.

Multiple-beam laser–plasma interactions in inertial confinement fusiona)

Abstract: The experimental evidence for multiple-beam laser-plasma instabilities of relevance to laser driven inertial confinement fusion at the ignition scale is reviewed, in both the indirect and direct-drive approaches. The instabilities described are cross-beam energy transfer (in both indirectly driven targets on the NIF and in direct-drive targets), multiple-beam stimulated Raman scattering (for indirect-drive), and multiple-beam two-plasmon decay instability (in direct drive). Advances in theoretical understanding and in the numerical modeling of these multiple beam instabilities are presented.

A two-dimensional (azimuthal-axial) particle-in-cell model of a Hall thruster

Abstract: We have developed a two-dimensional Particle-In-Cell model in the azimuthal and axial directions of the Hall thruster. A scaling method that consists to work at a lower plasma density to overcome constraints on time-step and grid-spacing is used. Calculations are able to reproduce the breathing mode due to a periodic depletion of neutral atoms without the introduction of a supplementary anomalous mechanism, as in fluid and hybrid models. Results show that during the increase of the discharge current, an electron-cyclotron drift instability (frequency in the range of MHz and wave number on the order of 3000 rad s−1) is formed in the region of the negative gradient of magnetic field. During the current decrease, an axial electric wave propagates from the channel toward the exhaust (whose frequency is on the order of 400 kHz) leading to a broadening of the ion energy distribution function. A discussion about the influence of the scaling method on the calculation results is also proposed.