T. Chapman

Electron and ion kinetic effects on non-linearly driven electron plasma and ion acoustic waves

Fully non-linear kinetic simulations of electron plasma and ion acoustic waves (IAWs) have been carried out with a new multi-species, parallelized Vlasov code. The numerical implementation of the Vlasov model and the methods used to compute the wave frequency are described in detail. For the first time, the nonlinear frequency of IAWs, combining the contributions from electron and ion kinetic effects and from harmonic generation, has been calculated and compared to Vlasov results. Excellent agreement of theory with simulation results is shown at all amplitudes, harmonic generation being an essential component at large amplitudes. For IAWs, the positive frequency shift from trapped electrons is confirmed and is dominant for the effective electron-to-ion temperature ratio, Z Te/Ti ≳ 10 with Z as the charge state. Furthermore, numerical results demonstrate unambiguously the dependence [R. L. Dewar, Phys. Fluids 15, 712 (1972)] of the kinetic shifts on details of the distribution of the trapped particles, whi...

High Power Dynamic Polarization Control Using Plasma Photonics

We report the first experimental demonstration of a plasma wave plate based on laser-induced birefringence. An elliptically polarized input was converted into a nearly ideal circularly polarized beam using an optical system composed of a second laser beam and a plasma. The results are in excellent agreement with linear theory and three-dimensional simulations up to phase delays exceeding π/4, thus establishing the feasibility of laser-plasma photonic devices that are ultrafast, damage-resistant, and easily tunable.

Exploring the limits of case-to-capsule ratio, pulse length, and picket energy for symmetric hohlraum drive on the National Ignition Facility Laser

We present a data-based model for low mode asymmetry in low gas-fill hohlraum experiments on the National Ignition Facility {NIF [Moses et al., Fusion Sci. Technol. 69, 1 (2016)]} laser. This model is based on the hypothesis that the asymmetry in these low fill hohlraums is dominated by the hydrodynamics of the expanding, low density, high-Z (gold or uranium) “bubble,” which occurs where the intense outer cone laser beams hit the high-Z hohlraum wall. We developed a simple model which states that the implosion symmetry becomes more oblate as the high-Z bubble size becomes large compared to the hohlraum radius or the capsule size becomes large compared to the hohlraum radius. This simple model captures the trends that we see in data that span much of the parameter space of interest for NIF ignition experiments. We are now using this model as a constraint on new designs for experiments on the NIF.

New insights into the decay of ion waves to turbulence, ion heating, and soliton generation

The decay of a single-frequency, propagating ion acoustic wave (IAW) via two-ion wave decay to a continuum of IAW modes is found to result in a highly turbulent plasma, ion soliton production, and rapid ion heating. Instability growth rates, thresholds, and sensitivities to plasma conditions are studied via fully kinetic Vlasov simulations. The decay rate of IAWs is found to scale linearly with the fundamental IAW potential amplitude ϕ1 for ZTe/Ti≲20, beyond which the instability is shown to scale with a higher power of ϕ1, where Z is the ion charge number and Te (Ti) is the electron (ion) thermal temperature. The threshold for instability is found to be smaller by an order of magnitude than linear theory estimates. Achieving a better understanding of the saturation of stimulated Brillouin scatter levels observed in laser-plasma interaction experiments is part of the motivation for this study.

Quantitative study of the trapped particle bunching instability in Langmuir waves

The bunching instability of particles trapped in Langmuir waves is studied using Vlasov simulations. A measure of particle bunching is defined and used to extract the growth rate from numerical simulations, which are compared with theory [Dodin et al., Phys. Rev. Lett. 110, 215006 (2013)]. In addition, the general theory of trapped particle instability in 1D is revisited and a more accurate description of the dispersion relation is obtained. Excellent agreement between numerical and theoretical predictions of growth rates of the bunching instability is shown over a range of parameters.

Resonance between heat-carrying electrons and Langmuir waves in inertial confinement fusion plasmas

In ignition scale hot plasmas, temperature gradients and thermal transport modify electron distributions in a velocity range resonant with Langmuir waves typical of those produced by stimulated Raman scattering. We examine the resultant changes to the Landau damping experienced by these Langmuir waves and the levels of thermal plasma fluctuations. The form factor and Thomson scattering cross-section in such plasmas display unique characteristics of the background conditions. A theoretical model and high-order Vlasov-Fokker-Planck simulations are used in our analysis. An experiment to measure changes in thermal plasma fluctuation levels due to a thermal gradient is proposed.

Demonstrating the saturation of stimulated Brillouin scattering by ion acoustic decay using fully kinetic simulations

The saturation of stimulated Brillouin scattering (SBS) by the decay to turbulence of the ion acoustic wave (IAW) that participates in the three-wave SBS interaction is demonstrated using a quasi-noiseless one-dimensional numerical solution to the Vlasov-Maxwell system of equations. This simulation technique permits careful examination of the decay process and its role in the complex evolution of SBS. The IAW decay process is shown to be an effective SBS saturation mechanism. In our example, the instantaneous plasma reflectivity saturates at ∼30% and drops to ∼0% as a direct consequence of IAW decay. A contrasting example where the reflectivity is controlled by dephasing due to the nonlinear frequency of the IAW is also discussed.

Longitudinal and Transverse Instability of Ion Acoustic Waves

Ion acoustic waves are found to be susceptible to at least two distinct decay processes. Which process dominates depends on the parameters. In the cases examined, the decay channel where daughter modes propagate parallel to the mother mode is found to dominate at larger amplitudes, while the decay channel where the daughter modes propagate at angles to the mother mode dominates at smaller amplitudes. Both decay processes may occur simultaneously and with onset thresholds below those suggested by fluid theory, resulting in the eventual multidimensional collapse of the mother mode to a turbulent state.

A plasma amplifier to combine multiple beams at NIF

Combining laser beams in a plasma is enabled by seeded stimulated Brillouin scattering which allows cross-beam energy transfer (CBET) to occur and re-distributes the energy between beams that cross with different incident angles and small differences in wavelength [Kirkwood et al. Phys. Plasmas 4, 1800 (1997)]. Indirect-drive implosions at the National Ignition Facility (NIF) [Haynam et al. Appl. Opt. 46, 3276–3303 (2007)] have controlled drive symmetry by using plasma amplifiers to transfer energy between beams [Kirkwood et al., Plasma Phys. Controlled Fusion 55, 103001 (2013); Lindl et al., Phys. Plasmas 21, 020501 (2014); and Hurricane et al. Nature 506, 343–348 (2014)]. In this work, we show that the existing models are well enough validated by experiments to allow a design of a plasma beam combiner that, once optimized, is expected to produce a pulse of light in a single beam with the energy greatly enhanced over existing sources. The scheme combines up to 61 NIF beams with 120 kJ of available energy...

A tesselation-based model for intensity estimation and laser plasma interactions calculations in three dimensions

A three-dimensional laser propagation model for computation of laser-plasma interactions is presented. It is focused on indirect drive geometries in inertial confinement fusion and formulated for use at large temporal and spatial scales. A modified tesselation-based estimator and a relaxation scheme are used to estimate the intensity distribution in plasma from geometrical optics rays. Comparisons with reference solutions show that this approach is well-suited to reproduce realistic 3D intensity field distributions of beams smoothed by phase plates. It is shown that the method requires a reduced number of rays compared to traditional rigid-scale intensity estimation. Using this field estimator, we have implemented laser refraction, inverse-bremsstrahlung absorption, and steady-state crossed-beam energy transfer with a linear kinetic model in the numerical code Vampire. Probe beam amplification and laser spot shapes are compared with experimental results and pf3d paraxial simulations. These results are promising for the efficient and accurate computation of laser intensity distributions in holhraums, which is of importance for determining the capsule implosion shape and risks of laser-plasma instabilities such as hot electron generation and backscatter in multi-beam configurations.A three-dimensional laser propagation model for computation of laser-plasma interactions is presented. It is focused on indirect drive geometries in inertial confinement fusion and formulated for use at large temporal and spatial scales. A modified tesselation-based estimator and a relaxation scheme are used to estimate the intensity distribution in plasma from geometrical optics rays. Comparisons with reference solutions show that this approach is well-suited to reproduce realistic 3D intensity field distributions of beams smoothed by phase plates. It is shown that the method requires a reduced number of rays compared to traditional rigid-scale intensity estimation. Using this field estimator, we have implemented laser refraction, inverse-bremsstrahlung absorption, and steady-state crossed-beam energy transfer with a linear kinetic model in the numerical code Vampire. Probe beam amplification and laser spot shapes are compared with experimental results and pf3d paraxial simulations. These results are pro...