scispace - formally typeset
Search or ask a question

Showing papers in "Physics of Plasmas in 2019"


Journal ArticleDOI
TL;DR: In this paper, the authors present a pedagogical review of hydrodynamic instability-induced turbulent flows, including Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz instabilities.
Abstract: In diverse areas of science and technology, including inertial confinement fusion (ICF), astrophysics, geophysics, and engineering processes, turbulent mixing induced by hydrodynamic instabilities is of scientific interest as well as practical significance. Because of the fundamental roles they often play in ICF and other applications, three classes of hydrodynamic instability-induced turbulent flows—those arising from the Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz instabilities—have attracted much attention. ICF implosions, supernova explosions, and other applications illustrate that these phases of instability growth do not occur in isolation, but instead are connected so that growth in one phase feeds through to initiate growth in a later phase. Essentially, a description of these flows must encompass both the temporal and spatial evolution of the flows from their inception. Hydrodynamic instability will usually start from potentially infinitesimal spatial perturbations, will eventually transition to a turbulent flow, and then will reach a final state of a true multiscale problem. Indeed, this change in the spatial scales can be vast, with hydrodynamic instability evolving from just a few microns to thousands of kilometers in geophysical or astrophysical problems. These instabilities will evolve through different stages before transitioning to turbulence, experiencing linear, weakly, and highly nonlinear states. The challenges confronted by researchers are enormous. The inherent difficulties include characterizing the initial conditions of such flows and accurately predicting the transitional flows. Of course, fully developed turbulence, a focus of many studies because of its major impact on the mixing process, is a notoriously difficult problem in its own right. In this pedagogical review, we will survey challenges and progress, and also discuss outstanding issues and future directions.

138 citations


Journal ArticleDOI
TL;DR: Wendelstein 7-X as discussed by the authors is the first comprehensively optimized stellarator aiming at good confinement with plasma parameters relevant to a future stellarator power plant, which achieved the highest triple product (6.5× 1019 keV m−3 s) achieved in a stellarator until now.
Abstract: Wendelstein 7-X is the first comprehensively optimized stellarator aiming at good confinement with plasma parameters relevant to a future stellarator power plant. Plasma operation started in 2015 using a limiter configuration. After installing an uncooled magnetic island divertor, extending the energy limit from 4 to 80 MJ, operation continued in 2017. For this phase, the electron cyclotron resonance heating (ECRH) capability was extended to 7 MW, and hydrogen pellet injection was implemented. The enhancements resulted in the highest triple product (6.5 × 1019 keV m−3 s) achieved in a stellarator until now. Plasma conditions [Te(0) ≈ Ti(0) ≈ 3.8 keV, τE > 200 ms] already were in the stellarator reactor-relevant ion-root plasma transport regime. Stable operation above the 2nd harmonic ECRH X-mode cutoff was demonstrated, which is instrumental for achieving high plasma densities in Wendelstein 7-X. Further important developments include the confirmation of low intrinsic error fields, the observation of current-drive induced instabilities, and first fast ion heating and confinement experiments. The efficacy of the magnetic island divertor was instrumental in achieving high performance in Wendelstein 7-X. Symmetrization of the heat loads between the ten divertor modules could be achieved by external resonant magnetic fields. Full divertor power detachment facilitated the extension of high power plasmas significantly beyond the energy limit of 80 MJ.

88 citations


Journal ArticleDOI
TL;DR: The quantum hydrodynamics (QHD) results have not found application in astrophysics or in experiments in condensed matter physics as mentioned in this paper, and these results practically did not have any impact on the former quantum plasma theory community.
Abstract: Quantum plasmas are an important topic in astrophysics and high pressure laboratory physics for more than 50 years. In addition, many condensed matter systems, including the electron gas in metals, metallic nanoparticles, or electron-hole systems in semiconductors and heterostructures, exhibit—to some extent—plasmalike behavior. Among the key theoretical approaches that have been applied to these systems are quantum kinetic theory, Green function theory, quantum Monte Carlo, semiclassical and quantum molecular dynamics, and more recently, density functional theory simulations. These activities are in close contact with the experiments and have firmly established themselves in the fields of plasma physics, astrophysics, and condensed matter physics. About two decades ago, a second branch of quantum plasma theory emerged that is based on a quantum fluid description and has attracted a substantial number of researchers. The focus of these studies has been on collective oscillations and linear and nonlinear waves in quantum plasmas. Even though these papers pretend to address the same physical systems as the more traditional papers mentioned above, the former appear to form a rather closed community that is largely isolated from the rest of the field. The quantum hydrodynamics (QHD) results have—with a few exceptions—not found application in astrophysics or in experiments in condensed matter physics. Moreover, these results practically did not have any impact on the former quantum plasma theory community. One reason is the unknown accuracy of the QHD for dense plasmas. In this paper, we present a novel derivation, starting from reduced density operators that clearly point to the deficiencies of QHD, and we outline possible improvements. It is also to be noted that some of the QHD results have attracted negative attention being criticized as unphysical. Examples include the prediction of “novel attractive forces” between protons in an equilibrium quantum plasma, the notion of “spinning quantum plasmas,” or the new field of “quantum dusty plasmas.” In the present article, we discuss the latter system in some detail because it is a particularly disturbing case of formal theoretical investigations that are detached from physical reality despite bold and unproven claims of importance for, e.g., dense astrophysical plasmas or microelectronics. We stress that these deficiencies are not a problem of QHD itself, which is a powerful and efficient method, but rather are due to ignorance of its properties and limitations. We analyze the common flaws of these works and come up with suggestions to improve the situation of QHD applications to quantum plasmas.

74 citations


Journal ArticleDOI
TL;DR: Spaeth et al. as discussed by the authors reviewed the current state of detailed modeling of NIF implosions, the scaling to ignition from recent experiments that that modeling implies, and areas for future improvements in modeling technique that could increase understanding and further enhance predictive capabilities.
Abstract: The goal of an inertially confined, igniting plasma on the National Ignition Facility (NIF) [M. L. Spaeth, Fusion Sci. Technol. 69, 25 (2016)] remains elusive. However, there is a growing understanding of the factors that appear to be limiting current implosion performance. And with this understanding, the question naturally arises: What conditions will ultimately be required to achieve ignition, either by continuing to improve the quality of current implosions, or by hydrodynamically scaling those implosions to larger driver energies on some future facility? Given the complexity of NIF implosions, answering this question must rely heavily on sophisticated numerical simulations. In particular, those simulations must respect the three-dimensionality of real NIF implosions and also resolve the wide range of scales for the many perturbation sources that degrade them. This prospectus article reviews the current state of detailed modeling of NIF implosions, the scaling to ignition from recent experiments that that modeling implies, and areas for future improvements in modeling technique that could increase understanding and further enhance predictive capabilities. Given the uncertainties inherent in any extrapolation, particularly for a process as nonlinear as ignition, there will be no definitive answer on the requirements for ignition until it is actually demonstrated experimentally. However, with continuing improvements in modeling technique and a growing experience base from NIF, the requirements for ignition are becoming clearer.

72 citations


Journal ArticleDOI
TL;DR: In this article, the locus of National Ignition Facility (NIF) inertial confinement fusion (ICF) implosion data, in hot-spot burn-average areal density (ρR) and Brysk temperature (T) space, is shown and illustrates that several implosions are nearing a burning plasma state, where α-heating is the dominant source of plasma heating.
Abstract: The locus of National Ignition Facility (NIF) inertial confinement fusion (ICF) implosion data, in hot-spot burn-average areal density (ρR) and Brysk temperature (T) space, is shown and illustrates that several implosions are nearing a burning plasma state, where α-heating is the dominant source of plasma heating. A formula for diagnosing a burning plasma using measured/inferred data from ICF implosion experiments is given with the underlying derivation. Plotting ICF implosion performance against inferred hot-spot energy illustrates the key need to maximize the delivery of energy to an implosion hot-spot. A very compact analytical equation for α-heating is given, which shows that fundamentally, α-heating provides a discrete “boost” to pVγ of an implosion hot-spot. It is then shown numerically that the analytical expression for the amplification of pVγ is simply related to yield amplification and to other more famous metrics used for inferring yield amplification in experiments. Interestingly, the argument of the pVγ boost equation appears to provide a fundamental ignition condition that implicitly includes the effects of asymmetry, and it also exposes the origin of why there is uncertainty in defining single ignition criteria. The resulting analysis of NIF implosion data indicates that an increase in burn-average hot-spot temperature will be needed in order to ignite, and the strategy being pursued to achieve this goal is outlined.

70 citations


Journal ArticleDOI
TL;DR: In this paper, the authors used 2.5D kinetic particle-in-cell simulations to simulate reconnection conditions appropriate for the magnetosheath and solar wind, i.e., plasma beta greater than 1 and low magnetic shear (strong guide field).
Abstract: Using 2.5 dimensional kinetic particle-in-cell simulations, we simulate reconnection conditions appropriate for the magnetosheath and solar wind, i.e., plasma beta (ratio of gas pressure to magnetic pressure) greater than 1 and low magnetic shear (strong guide field). Changing the simulation domain size, we find that the ion response varies greatly. For reconnecting regions with scales comparable to the ion inertial length, the ions do not respond to the reconnection dynamics leading to “electron-only” reconnection with very large quasisteady reconnection rates. Note that in these simulations, the ion Larmor radius is comparable to the ion inertial length. The transition to a more traditional “ion-coupled” reconnection is gradual as the reconnection domain size increases, with the ions becoming frozen-in in the exhaust when the magnetic island width in the normal direction reaches many ion inertial lengths. During this transition, the quasisteady reconnection rate decreases until the ions are fully coupled, ultimately reaching an asymptotic value. The scaling of the ion outflow velocity with the exhaust width during this electron-only to ion-coupled transition is found to be consistent with a theoretical model of a newly reconnected field line. In order to have a fully frozen-in ion exhaust with ion flows comparable to the reconnection Alfven speed, an exhaust width of at least several ion inertial lengths is needed. In turbulent systems with reconnection occurring between magnetic bubbles associated with fluctuations, using geometric arguments, we estimate that fully ion-coupled reconnection requires magnetic bubble length scales of at least several tens of ion inertial lengths.

53 citations


Journal ArticleDOI
TL;DR: In this article, the authors analyzed linear slow magnetoacoustic waves in a plasma in thermal equilibrium formed by a balance of optically thin radiative losses, field-align thermal conduction, and an unspecified heating process.
Abstract: Slow magnetoacoustic waves are omnipresent in both natural and laboratory plasma systems. The wave-induced misbalance between plasma cooling and heating processes causes the amplification or attenuation, and also dispersion, of slow magnetoacoustic waves. The wave dispersion could be attributed to the presence of characteristic time scales in the system, connected with the plasma heating or cooling due to the competition of the heating and cooling processes in the vicinity of thermal equilibrium. We analyzed linear slow magnetoacoustic waves in a plasma in thermal equilibrium formed by a balance of optically thin radiative losses, field-align thermal conduction, and an unspecified heating process. The dispersion is manifested by the dependence of the effective adiabatic index of the wave on the wave frequency, making the phase and group speeds frequency-dependent. The mutual effect of the wave amplification and dispersion is shown to result in the occurrence of an oscillatory pattern in an initially broadband slow wave, with the characteristic period determined by the thermal misbalance time scales, i.e., by the derivatives of the combined radiation loss and heating function with respect to the density and temperature, evaluated at the equilibrium. This effect is illustrated by estimating the characteristic period of the oscillatory pattern, appearing because of thermal misbalance in the plasma of the solar corona. It is found that by an order of magnitude, the period is about the typical periods of slow magnetoacoustic oscillations detected in the corona.

53 citations


Journal ArticleDOI
TL;DR: In this article, the authors show that at low collisionality (ν*e < 0.5), low pedestal density is required for resonant field penetration at the pedestal top (ne,ped ≈ 2.5
Abstract: The density dependence of edge-localized-mode (ELM) suppression and density pump-out (density reduction) by n = 2 resonant magnetic perturbations (RMPs) is consistent with the effects of narrow well-separated magnetic islands at the top and bottom of the H-mode pedestal in DIII-D low-collisionality plasmas. Nonlinear two-fluid MHD simulations for DIII-D ITER similar shape discharges show that, at low collisionality (ν*e < 0.5), low pedestal density is required for resonant field penetration at the pedestal top (ne,ped ≈ 2.5 × 1019 m−3 at ψN ≈ 0.93), consistent with the ubiquitous low density requirement for ELM suppression in these DIII-D plasmas. The simulations predict a drop in the pedestal pressure due to parallel transport across these narrow width (ΔψN ≈ 0.02) magnetic islands at the top of the pedestal that is stabilizing to Peeling-Ballooning-Modes and comparable to the pedestal pressure reduction observed in experiment at the onset of ELM suppression. The simulations predict density pump-out at experimentally relevant levels (Δne/ne ≈ −20%) at low pedestal collisionality (ν*e ≈ 0.1) due to very narrow (ΔψN ≈ 0.01–0.02) RMP driven magnetic islands at the pedestal foot at ψN ≈ 0.99. The simulations show decreasing pump-out with increasing density, consistent with experiment, resulting from the inverse dependence of parallel particle transport on collisionality at the foot of the pedestal. The robust screening of resonant fields is predicted between the top and bottom of the pedestal during density pump-out and ELM suppression, consistent with the preservation of strong temperature gradients in the edge transport barrier as seen in experiment.The density dependence of edge-localized-mode (ELM) suppression and density pump-out (density reduction) by n = 2 resonant magnetic perturbations (RMPs) is consistent with the effects of narrow well-separated magnetic islands at the top and bottom of the H-mode pedestal in DIII-D low-collisionality plasmas. Nonlinear two-fluid MHD simulations for DIII-D ITER similar shape discharges show that, at low collisionality (ν*e < 0.5), low pedestal density is required for resonant field penetration at the pedestal top (ne,ped ≈ 2.5 × 1019 m−3 at ψN ≈ 0.93), consistent with the ubiquitous low density requirement for ELM suppression in these DIII-D plasmas. The simulations predict a drop in the pedestal pressure due to parallel transport across these narrow width (ΔψN ≈ 0.02) magnetic islands at the top of the pedestal that is stabilizing to Peeling-Ballooning-Modes and comparable to the pedestal ...

49 citations


Journal ArticleDOI
TL;DR: In this paper, a statistical framework for the calibration of ICF simulations using data collected at the National Ignition Facility (NIF) is described, and Bayesian inferences for a series of laser shots using an approach that is designed to respect the physics simulation as much as possible and then build a second model that links the individual-shot inferences together.
Abstract: Computer models of inertial confinement fusion (ICF) implosions play an essential role in experimental design and interpretation as well as our understanding of fundamental physics under the most extreme conditions that can be reached in the laboratory. Building truly predictive models is a significant challenge, with the potential to greatly accelerate progress to high yield and ignition. One path to more predictive models is to use experimental data to update the underlying physics in a way that can be extrapolated to new experiments and regimes. We describe a statistical framework for the calibration of ICF simulations using data collected at the National Ignition Facility (NIF). We perform Bayesian inferences for a series of laser shots using an approach that is designed to respect the physics simulation as much as possible and then build a second model that links the individual-shot inferences together. We show that this approach is able to match multiple X-ray and neutron diagnostics for a whole series of NIF “BigFoot” shots. Within the context of 2D radiation hydrodynamic simulations, our inference strongly favors a significant reduction in fuel compression over other known degradation mechanisms (namely, hohlraum issues and engineering perturbations). This analysis is expanded using a multifidelity technique to pick fuel-ablator mix from several candidate causes of the degraded fuel compression (including X-ray preheat and shock timing errors). Finally, we use our globally calibrated model to investigate the extra laser drive energy that would be required to overcome the inferred fuel compression issues in NIF BigFoot implosions.

49 citations


Journal ArticleDOI
TL;DR: In this paper, a new computational model suitable for exploring the self-consistent production of energetic electrons during magnetic reconnection in macroscale systems is presented, based on the recent discovery that parallel electric fields are ineffective drivers of energetic particles during reconnection so that the kinetic scales which control the development of such fields can be ordered out of the equations.
Abstract: A new computational model suitable for exploring the self-consistent production of energetic electrons during magnetic reconnection in macroscale systems is presented. The equations are based on the recent discovery that parallel electric fields are ineffective drivers of energetic particles during reconnection so that the kinetic scales which control the development of such fields can be ordered out of the equations. The resulting equations consist of a magnetohydrodynamic (MHD) backbone with the energetic component represented by macro-particles described by the guiding center equations. Crucially, the energetic component feeds back on the MHD equations so that the total energy of the MHD fluid and the energetic particles is conserved. The equations correctly describe the firehose instability, whose dynamics plays a key role in throttling reconnection and in controlling the spectra of energetic particles. The results of early tests of the model, including the propagation of Alfven waves in a system with pressure anisotropy and the growth of firehose modes, establish that the basic algorithm is stable and produces reliable physics results in preparation for further benchmarking with particle-in-cell models of reconnection.

46 citations


Journal ArticleDOI
TL;DR: An overview of recent advances in the field of electron kinetics in low-temperature plasmas (LTPs) can be found in this article, where the authors provide author's views on where the field is headed and suggest promising strategies for further development.
Abstract: This article presents an overview of recent advances in the field of electron kinetics in low-temperature plasmas (LTPs). It also provides author's views on where the field is headed and suggests promising strategies for further development. The authors have selected several problems to illustrate multidisciplinary nature of the subject (space and laboratory plasma, collisionless and collisional plasmas, and low-pressure and high-pressure discharges) and to illustrate how cross-disciplinary research efforts could enable further progress. Nonlocal electron kinetics and nonlocal electrodynamics in low-pressure rf plasmas resemble collisionless effects in space plasma and hot plasma effects in fusion science, terahertz technology, and plasmonics. The formation of electron groups in dc and rf discharges has much in common with three groups of electrons (core, strahl, and halo) in solar wind. Runaway electrons in LTPs are responsible for a wide range of physical phenomena from nano- and picoscale breakdown of dielectrics to lightning initiation. Understanding electron kinetics of LTPs could promote scientific advances in a number of topics in plasma physics and accelerate modern plasma technologies.

Journal ArticleDOI
TL;DR: In this article, the absolute instability thresholds for stimulated Raman scattering (SRS) and two-plasmon decay (TPD) driven by a broadband laser pulse are evaluated numerically.
Abstract: Absolute instability thresholds for stimulated Raman scattering (SRS) and two-plasmon decay (TPD) driven by a broadband laser pulse are evaluated numerically. The scalings of the calculated thresholds with the density scale length, temperature, and central wavelength are qualitatively similar to the existing analytical results. The threshold values, however, exhibit significant quantitative differences. Comparisons between the thresholds calculated for various broadband power spectra indicate a universal scaling of the threshold intensity with laser coherence time (τc). For SRS, Ithr∝τc−1/3, and for TPD, Ithr∝τc−1/2.

Journal ArticleDOI
TL;DR: Falessi and Zonca as discussed by the authors used gyrokinetic theory to extract the phase space zonal structure from the flux surface averaged particle response, that is, the nonlinear response that is undamped by collisionless processes.
Abstract: We adopt gyrokinetic theory to extract the phase space zonal structure from the flux surface averaged particle response, that is, the nonlinear response that is undamped by collisionless processes. We argue that phase space zonal structures are a proper definition for the nonlinear distortion of the plasma reference state and, thus, of the generally non-Maxwellian neighboring nonlinear equilibria consistent with toroidal symmetry breaking fluctuations. Evolution equations for phase space zonal structures are derived and discussed, along with the corresponding density and energy transport equations. It is shown that this approach is consistent with the usual evolution of macroscopic plasma profiles under the action of fluctuation induced fluxes, when the deviation of the reference state from local Maxwellian response is small. In particular, the present results recover those of a previous article [M. V. Falessi and F. Zonca, Phys. Plasmas 25, 032306 (2018)], where transport equations holding on the reference state length scale have been derived using the moment approach introduced in the classical review work by Hinton and Hazeltine.

Journal ArticleDOI
TL;DR: In this article, a generalized analytical approach is developed which shows good agreement with self-consistent quantum mechanical calculations, empirically discovered that anomalous strong scaling in the analytical model provides agreement with data obtained in a regime where the lattice structure still prevails.
Abstract: The radial dependence of the free electron density within the ion sphere radius in finite temperature dense plasmas shows characteristic scaling laws that permit us to derive analytical plasma screening potentials. A generalized analytical approach is developed which shows good agreement with self-consistent quantum mechanical calculations. It is empirically discovered that anomalous strong scaling in the analytical model provides agreement with data obtained in a regime where the lattice structure still prevails.

Journal ArticleDOI
TL;DR: In this paper, a micro-channel was aligned to be collinear with the incident laser pulse, confining the majority of the laser energy within the channel, and the measured electron spectrum showed a large increase in the cutoff energy and slope temperature when compared to that from a 2.5
Abstract: We present an experimental demonstration of the efficient acceleration of electrons beyond 60 MeV using micro-channel plasma targets. We employed a high-contrast, 2.5 J, 32 fs short pulse laser interacting with a 5 μm inner diameter, 300 μm long micro-channel plasma target. The micro-channel was aligned to be collinear with the incident laser pulse, confining the majority of the laser energy within the channel. The measured electron spectrum showed a large increase in the cut-off energy and slope temperature when compared to that from a 2 μm flat Copper target, with the cutoff energy more than doubled and the total energy in electrons >5 MeV enhanced by over 10 times. Three-dimensional particle-in-cell simulations confirm efficient direct laser acceleration enabled by the novel structure as the dominant acceleration mechanism for the high energy electrons. The simulations further reveal the guiding effect of the channel that successfully explains preferential acceleration on the laser/channel axis observed in experiments. Finally, systematic simulations provide scalings for the energy and charge of the electron pulses. Our results show that the micro-channel plasma target is a promising electron source for applications such as ion acceleration, Bremsstrahlung X-ray radiation, and THZ generation.

Journal ArticleDOI
TL;DR: In this paper, a negative triangularity shape has been created on the DIII-D tokamak that, despite maintaining standard L-mode edge radial profiles, reach volume averaged pressure levels typical of H-mode scenarios.
Abstract: Plasmas with a negative triangularity shape have been created on the DIII-D tokamak that, despite maintaining standard L-mode edge radial profiles, reach volume averaged pressure levels typical of H-mode scenarios. Within the auxiliary power available for these experiments, plasmas exhibit near-zero power degradation while sustaining βN = 2.7 and H98,y2 = 1.2 for several energy confinement times. Detailed comparison with matched discharges at positive triangularity indicates that Trapped Electron Modes are weakened at negative triangularity, consistent with increased confinement and reduced intensity of fluctuations in electron density, electron temperature, and ion density. These results indicate that a negative triangularity plasma operating without an edge pedestal might provide an attractive scenario for operations in future reactors.

Journal ArticleDOI
TL;DR: In this article, the authors developed a general quasi-optical theory for mode-converting electromagnetic beams in a weakly inhomogeneous linear medium with no sources and implemented it in a numerical algorithm.
Abstract: This work opens a series of papers where we develop a general quasi-optical theory for mode-converting electromagnetic beams in plasma and implement it in a numerical algorithm. Here, the basic theory is introduced. We consider a general quasimonochromatic multicomponent wave in a weakly inhomogeneous linear medium with no sources. For any given dispersion operator that governs the wave field, we explicitly calculate the approximate operator that governs the wave envelope ψ to the second order in the geometrical-optics parameter. Then, we further simplify this envelope operator by assuming that the gradient of ψ transverse to the local group velocity is much larger than the corresponding parallel gradient. This leads to a parabolic differential equation for ψ (“quasioptical equation”) on the basis of the geometrical-optics polarization vectors. Scalar and mode-converting vector beams are described on the same footing. We also explain how to apply this model to electromagnetic waves in general. In the next papers of this series, we report successful quasioptical modeling of radio frequency wave beams in magnetized plasma based on this theory.This work opens a series of papers where we develop a general quasi-optical theory for mode-converting electromagnetic beams in plasma and implement it in a numerical algorithm. Here, the basic theory is introduced. We consider a general quasimonochromatic multicomponent wave in a weakly inhomogeneous linear medium with no sources. For any given dispersion operator that governs the wave field, we explicitly calculate the approximate operator that governs the wave envelope ψ to the second order in the geometrical-optics parameter. Then, we further simplify this envelope operator by assuming that the gradient of ψ transverse to the local group velocity is much larger than the corresponding parallel gradient. This leads to a parabolic differential equation for ψ (“quasioptical equation”) on the basis of the geometrical-optics polarization vectors. Scalar and mode-converting vector beams are described on the same footing. We also explain how to apply this model to electromagnetic waves in general. In the next...

Journal ArticleDOI
TL;DR: In this article, the authors report on a high convergence ratio liquid layer capsule implosion performed on the National Ignition Facility and contrast it to two previously reported layered implosions, in order to better understand how the capsule design impacts the hydrodynamic stability properties.
Abstract: We report on a high convergence ratio liquid layer capsule implosion performed on the National Ignition Facility and contrast it to two previously reported layered implosions, in order to better understand how the capsule design impacts the hydrodynamic stability properties of implosions. Three implosions were performed with similar convergence ratios, fuel entropy, in-flight aspect ratios, and unablated shell mass; these qualities are important for determining hydrodynamic stability. Nevertheless, while two of these implosions exhibited robustness to asymmetries, including our recent experiment that had abnormally large amplitude long-wavelength capsule asymmetries, and produced more than 80% or the yield predicted by one-dimensional (1D) simulations, which do not account for the impacts of hydrodynamic instabilities, the third implosion produced only 14% of the yield from a 1D simulation. We perform a detailed computational analysis of these three shots, which suggests that the combination of several large asymmetry seeds result in the significantly degraded performance: a large 30 μm fill tube, the presence of a microstructure in the high density carbon ablator, and a higher level of drive asymmetry. This indicates that while it is possible to stabilize a high convergence ratio implosion through various means, the factors that determine stability cannot be considered independently. Furthermore, when these asymmetries are combined in 2D simulations, they can exhibit destructive interference and underpredict the yield degradation compared to experiment and three-dimensional simulations.

Journal ArticleDOI
TL;DR: In this paper, the effects of mix on target performance in magnetized liners inertial fusion (MagLIF) experiments at Sandia National Laboratories were analyzed, showing that mix is likely produced from a variety of sources, approximately half of which originates during the laser heating phase and the remainder near stagnation.
Abstract: In magneto-inertial-fusion experiments, energy losses such as a radiation need to be well controlled in order to maximize the compressional work done on the fuel and achieve thermonuclear conditions. One possible cause for high radiation losses is high-Z material mixing from the target components into the fuel. In this work, we analyze the effects of mix on target performance in Magnetized Liner Inertial Fusion (MagLIF) experiments at Sandia National Laboratories. Our results show that mix is likely produced from a variety of sources, approximately half of which originates during the laser heating phase and the remainder near stagnation, likely from the liner deceleration. By changing the “cushion” component of MagLIF targets from Al to Be, we achieved a 10× increase in neutron yield, a 60% increase in ion temperature, and an ∼50% increase in fuel energy at stagnation.

Journal ArticleDOI
TL;DR: In this paper, 3D particle-in-cell simulations are used to explore a more extensive laser-target parameter space than previously reported in the literature as well as to study the laser pulse coupling to the target, the structure of the fields, and the properties of the accelerated ion beam in the MVA scheme.
Abstract: Magnetic Vortex Acceleration (MVA) from near critical density targets is one of the promising schemes of laser-driven ion acceleration. 3D particle-in-cell simulations are used to explore a more extensive laser-target parameter space than previously reported in the literature as well as to study the laser pulse coupling to the target, the structure of the fields, and the properties of the accelerated ion beam in the MVA scheme. The efficiency of acceleration depends on the coupling of the laser energy to the self-generated channel in the target. The accelerated proton beams demonstrate a high level of collimation with achromatic angular divergence, and carry a significant amount of charge. For petawatt-class lasers, this acceleration regime provides a favorable scaling of the maximum ion energy with the laser power for the optimized interaction parameters. The megatesla-level magnetic fields generated by the laser-driven coaxial plasma structure in the target are a prerequisite for accelerating protons to the energy of several hundred mega-electron-volts.

Journal ArticleDOI
TL;DR: In this article, the authors present simulation and experimental results on momentum transfer to different layers in a double shell and also present the details of the NIF cylindrical hohlraum double shell platform including an imaging shell design with a mid-Z inner shell necessary for imaging the inner shell shape and the trajectory.
Abstract: Advances in target fabrication have made double shell capsule implosions a viable platform to study burning fusion plasmas. Central to the double shell capsule is a high-Z (e.g., Au) metal pusher that accesses the volume-burn regime by reducing radiative losses through radiation trapping and compressing a uniform fuel volume at reduced velocities. A double shell implosion relies on a series of energy transfer processes starting from x-ray absorption by the outer shell, followed by transfer of kinetic energy to an inner shell, and finally conversion of kinetic energy to fuel internal energy. We present simulation and experimental results on momentum transfer to different layers in a double shell. We also present the details of the development of the NIF cylindrical hohlraum double shell platform including an imaging shell design with a mid-Z inner shell necessary for imaging the inner shell shape and the trajectory with the current 2DConA platform capability. We examine 1D energy transfer between shell layers using trajectory measurements from a series of surrogate targets; the series builds to a complete double shell layer by layer, isolating the physics of each step of the energy transfer process. The measured energy transfer to the foam cushion and the inner shell suggests that our radiation-hydrodynamics simulations capture most of the relevant collision physics. With a 1 MJ laser drive, the experimental data indicate that 22% ± 3% of the ablator kinetic energy couples into inner shell KE, compared to a 27% ± 2% coupling in our xRAGE simulations. Thus, our xRAGE simulations match experimental energy transfer to ∼5%, without inclusion of higher order 2D and 3D effects.

Journal ArticleDOI
TL;DR: In this paper, the authors used the Advanced Radiographic Capability (ARC) short-pulse laser at the National Ignition Facility in the Lawrence Livermore National Laboratory (LLNL).
Abstract: New short-pulse kilojoule, Petawatt-class lasers, which have recently come online and are coupled to large-scale, many-beam long-pulse facilities, undoubtedly serve as very exciting tools to capture transformational science opportunities in high energy density physics. These short-pulse lasers also happen to reside in a unique laser regime: very high-energy (kilojoule), relatively long (multi-picosecond) pulse-lengths, and large (10s of micron) focal spots, where their use in driving energetic particle beams is largely unexplored. Proton acceleration via Target Normal Sheath Acceleration (TNSA) using the Advanced Radiographic Capability (ARC) short-pulse laser at the National Ignition Facility in the Lawrence Livermore National Laboratory is demonstrated for the first time, and protons of up to 18 MeV are measured using laser irradiation of >1 ps pulse-lengths and quasi-relativistic (∼1018 W/cm2) intensities. This is indicative of a super-ponderomotive electron acceleration mechanism that sustains acceleration over long (multi-picosecond) time-scales and allows for proton energies to be achieved far beyond what the well-established scalings of proton acceleration via TNSA would predict at these modest intensities. Furthermore, the characteristics of the ARC laser (large ∼100 μm diameter focal spot, flat spatial profile, multi-picosecond, relatively low prepulse) provide acceleration conditions that allow for the investigation of 1D-like particle acceleration. A high flux ∼ 50 J of laser-accelerated protons is experimentally demonstrated. A new capability in multi-picosecond particle-in-cell simulation is applied to model the data, corroborating the high proton energies and elucidating the physics of multi-picosecond particle acceleration.

Journal ArticleDOI
TL;DR: In this article, the authors used 3D extended-magnetohydrodynamics Gorgon simulations to investigate how thermal conduction suppression, the Lorentz force, and α-particle magnetisation affect three hot-spot perturbation scenarios: a cold fuel spike, a time-dependent radiation drive asymmetry, and a multi-mode perturbations.
Abstract: Pre-magnetisation of inertial confinement fusion implosions on the National Ignition Facility has the potential to raise current high-performing targets into the ignition regime [Perkins et al. “The potential of imposed magnetic fields for enhancing ignition probability and fusion energy yield in indirect-drive inertial confinement fusion,” Phys. Plasmas 24, 062708 (2017)]. A key concern with this method is that the application of a magnetic field inherently increases asymmetry. This paper uses 3-D extended-magnetohydrodynamics Gorgon simulations to investigate how thermal conduction suppression, the Lorentz force, and α-particle magnetisation affect three hot-spot perturbation scenarios: a cold fuel spike, a time-dependent radiation drive asymmetry, and a multi-mode perturbation. For moderate magnetisations (B0 = 5 T), the single spike penetrates deeper into the hot-spot, as thermal ablative stabilisation is reduced. However, at higher magnetisations (B0 = 50 T), magnetic tension acts to stabilise the spike. While magnetisation of α-particle orbits increases the peak hot-spot temperature, no impact on the perturbation penetration depth is observed. The P4-dominated radiation drive asymmetry demonstrates the anisotropic nature of the thermal ablative stabilisation modifications, with perturbations perpendicular to the magnetic field penetrating deeper and perturbations parallel to the field being preferentially stabilised by increased heat-flows. Moderate magnetisations also increase the prevalence of high modes, while magnetic tension reduces vorticity at the hot-spot edge for larger magnetisations. For a simulated high-foot experiment, the yield doubles through the application of a 50 T magnetic field-an amplification which is expected to be larger for higher-performing configurations.

Journal ArticleDOI
TL;DR: In this article, the time-resolved cross-field electron anomalous collision frequency in a Hall thruster was inferred from minimally invasive laser-based measurements, which was employed to characterize the collision frequency.
Abstract: The time-resolved cross-field electron anomalous collision frequency in a Hall thruster is inferred from minimally invasive laser-based measurements. This diagnostic is employed to characterize the...

Journal ArticleDOI
TL;DR: Hariri et al. as mentioned in this paper extended the fluid-turbulence code GRILLIX to a global model with drift reduced Braginskii equations, which consistently accounts for large variation in plasma parameters.
Abstract: Turbulent dynamics in the scrape-off layer of magnetic fusion devices is intermittent with large fluctuations in density and pressure. Therefore, a model is required that allows perturbations of similar or even larger magnitude to the time-averaged background value. The fluid-turbulence code GRILLIX is extended to such a global model, which consistently accounts for large variation in plasma parameters. Derived from the drift reduced Braginskii equations, the new GRILLIX model includes electromagnetic and electron-thermal dynamics, retains global parametric dependencies and the Boussinesq approximation is not applied. The penalization technique is combined with the flux-coordinate independent approach [F. Hariri and M. Ottaviani, Comput. Phys. Commun. 184, 2419 (2013) and A. Stegmeir et al., Comput. Phys. Commun. 198, 139 (2016)], which allows to study realistic diverted geometries with X-point(s) and general boundary contours. We characterize results from turbulence simulations and investigate the effect of geometry by comparing simulations in circular geometry with toroidal limiter against realistic diverted geometry at otherwise comparable parameters. Turbulence is found to be intermittent with relative fluctuation levels of up to 40% showing that a global description is indeed important. At the same time via direct comparison, we find that the Boussinesq approximation has only a small quantitative impact in a turbulent environment. In comparison to circular geometry, the fluctuations are reduced in diverted geometry, which is related to a different zonal flow structure. Moreover, the fluctuation level has a more complex spatial distribution in diverted geometry. Due to local magnetic shear, which differs fundamentally in circular and diverted geometries, turbulent structures become strongly distorted in the perpendicular direction and are eventually damped away toward the X-point.

Journal ArticleDOI
TL;DR: In this article, the crossover between high energy density plasmas (HEDP) and Ultracold Neutral Plasmas was discussed and theoretical and computational models were developed to understand a subset of HEDP experiments.
Abstract: In this paper, we present ideas that were part of the miniconference on the crossover between High Energy Density Plasmas (HEDP) and Ultracold Neutral Plasmas (UNPs) at the 60th Annual Meeting of the American Physical Society Division of Plasma Physics, November 2018. We give an overview of UNP experiments with an emphasis on measurements of the time-evolving ion density and velocity distributions, the electron-ion thermalization rate, and plasma self-assembly—all just inside the strongly coupled plasma regime. We also present theoretical and computational models that were developed to understand a subset of HEDP experiments. However, because HEDP experiments display similar degrees of strong coupling, many aspects of these models can be vetted using precision studies of UNPs. This comparison is important because some statistical assumptions used for ideal plasmas are of questionable validity in the strongly coupled plasma regime. We summarize two theoretical approaches that extend kinetic theories into the strong-coupling regime and show good agreement for momentum transfer and self-diffusion. As capabilities improve, both computationally and experimentally, UNP measurements may help guide the ongoing development of HEDP-appropriate plasma models. Future opportunities in viscosity, energy relaxation, and magnetized plasmas are discussed.In this paper, we present ideas that were part of the miniconference on the crossover between High Energy Density Plasmas (HEDP) and Ultracold Neutral Plasmas (UNPs) at the 60th Annual Meeting of the American Physical Society Division of Plasma Physics, November 2018. We give an overview of UNP experiments with an emphasis on measurements of the time-evolving ion density and velocity distributions, the electron-ion thermalization rate, and plasma self-assembly—all just inside the strongly coupled plasma regime. We also present theoretical and computational models that were developed to understand a subset of HEDP experiments. However, because HEDP experiments display similar degrees of strong coupling, many aspects of these models can be vetted using precision studies of UNPs. This comparison is important because some statistical assumptions used for ideal plasmas are of questionable validity in the strongly coupled plasma regime. We summarize two theoretical approaches that extend kinetic theories into t...

Journal ArticleDOI
TL;DR: In this paper, the authors review some of the mass filter concepts proposed to date and underline how the practicality of these concepts is conditioned upon the ability to sustain a suitable perpendicular electric field in a plasma for parameters compatible with high-throughput operation.
Abstract: High-throughput plasma separation based on atomic mass holds promise for offering unique solutions to a variety of high-impact societal applications. Through the mass differential effects they exhibit, crossed-field configurations can in principle be exploited in various ways to separate ions based on atomic mass. Here, we review some of the E × B mass filter concepts proposed to date and underline how the practicality of these concepts is conditioned upon the ability to sustain a suitable perpendicular electric field in a plasma for parameters compatible with high-throughput operation. We show that while the limited present predictive capabilities do not make it possible to confirm this possibility, past experimental results suggest that end-electrode biasing may be effective, at least for certain electric field values. We conclude that a better understanding of cross-field conductivity is needed to confirm these results and confirm the potential of crossed-field configurations for high-throughput separation.

Journal ArticleDOI
TL;DR: In this article, the authors describe a systematic development of kinetic entropy as a diagnostic in fully kinetic particle-in-cell (PIC) simulations and use it to interpret plasma physics processes in heliospheric, planetary, and astrophysical systems.
Abstract: We describe a systematic development of kinetic entropy as a diagnostic in fully kinetic particle-in-cell (PIC) simulations and use it to interpret plasma physics processes in heliospheric, planetary, and astrophysical systems. First, we calculate kinetic entropy in two forms—the “combinatorial” form related to the logarithm of the number of microstates per macrostate and the “continuous” form related to flnf, where f is the particle distribution function. We discuss the advantages and disadvantages of each and discuss subtleties about implementing them in PIC codes. Using collisionless PIC simulations that are two-dimensional in position space and three-dimensional in velocity space, we verify the implementation of the kinetic entropy diagnostics and discuss how to optimize numerical parameters to ensure accurate results. We show the total kinetic entropy is conserved to three percent in an optimized simulation of antiparallel magnetic reconnection. Kinetic entropy can be decomposed into a sum of a position space entropy and a velocity space entropy, and we use this to investigate the nature of kinetic entropy transport during collisionless reconnection. We find the velocity space entropy of both electrons and ions increases in time due to plasma heating during magnetic reconnection, while the position space entropy decreases due to plasma compression. This project uses collisionless simulations, so it cannot address physical dissipation mechanisms; nonetheless, the infrastructure developed here should be useful for studies of collisional or weakly collisional heliospheric, planetary, and astrophysical systems. Beyond reconnection, the diagnostic is expected to be applicable to plasma turbulence and collisionless shocks.We describe a systematic development of kinetic entropy as a diagnostic in fully kinetic particle-in-cell (PIC) simulations and use it to interpret plasma physics processes in heliospheric, planetary, and astrophysical systems. First, we calculate kinetic entropy in two forms—the “combinatorial” form related to the logarithm of the number of microstates per macrostate and the “continuous” form related to flnf, where f is the particle distribution function. We discuss the advantages and disadvantages of each and discuss subtleties about implementing them in PIC codes. Using collisionless PIC simulations that are two-dimensional in position space and three-dimensional in velocity space, we verify the implementation of the kinetic entropy diagnostics and discuss how to optimize numerical parameters to ensure accurate results. We show the total kinetic entropy is conserved to three percent in an optimized simulation of antiparallel magnetic reconnection. Kinetic entropy can be decomposed into a sum of a posit...

Journal ArticleDOI
Charlson C. Kim, Yueqiang Liu1, Paul Parks1, Lang L. Lao1, Michael Lehnen2, A. Loarte2 
TL;DR: Izzo et al. as mentioned in this paper developed a particle-based SPI model in the NIMROD code coupled to its modified single fluid equations with impurity and radiation, and validated simulations of the thermal quench and comparisons to DIII-D impurity scan experiments.
Abstract: Optimal strategies for disruption mitigation benefit from the understanding of details both spatially and temporally. Beyond the assessment of the efficacy of a particular proposed Disruption Mitigation System (DMS), ITER's longevity will require accounting of both mitigated and unmitigated disruptions. Accurate models and validated simulations that detail multiple ITER scenarios with mitigated and unmitigated disruptions are essential for accurate estimates of load damage. The primary candidate for ITER's DMS is Shattered Pellet Injection (SPI); its efficacy must be evaluated within the next several years. To perform critical time dependent 3-D nonlinear simulations, we have developed a particle based SPI model in the NIMROD code coupled to its modified single fluid equations with impurity and radiation [Izzo, Nucl. Fusion 46(5), 541 (2006)]. SPI validation simulations of the thermal quench and comparisons to DIII-D impurity scan experiments [Shiraki et al., Phys. Plasmas 23(6), 062516 (2016)] are presented. We also present an initial ITER Q = 10 pure neon SPI simulation and compare it with the DIII-D SPI simulations. NIMROD SPI simulations demonstrate that the ablating fragment drives strong parallel flows that transport the impurities and governs the thermal quench. Analysis of SPI simulations shows that the mixed deuterium/neon SPI results in a more benign thermal quench due to the enhanced transport caused by the additional deuterium. These results suggest that an optimal pellet mixture exists for the SPI system.

Journal ArticleDOI
TL;DR: In this paper, the authors presented the first systematic experimental study of indirectly driven beryllium capsules with a cryogenic deuterium-tritium fuel layer, and they showed that the capsule scale is predicted to reduce the impact of fuel-ablator mix and achieve high performance.
Abstract: Many inertial fusion designs use capsules made of beryllium, as its high mass ablation rate is advantageous. We present the first systematic experimental study of indirectly driven beryllium capsules with a cryogenic deuterium-tritium fuel layer. “Subscale” capsules, 80% of the nominal National Ignition Facility point design radius, show optimal performance with the remaining mass of ∼6–7%. A buoyancy-drag mix model explains the implosion performance, suggesting that fuel-ablator mix is the dominant degradation mechanism. Increasing the capsule scale is predicted to reduce the impact of fuel-ablator mix and achieve high performance.