Low-mode Perturbations Research Articles

Mitigation of mode-one asymmetry in laser-direct-drive inertial confinement fusion implosions

Nonuniformities present in the laser illumination and target in laser-driven inertial confinement fusion experiments lead to an asymmetric compression of the target, resulting in an inefficient conversion of shell kinetic energy to thermal energy of the hot-spot plasma. In this paper, the effects of asymmetric compression of cryogenic deuterium tritium laser-direct-drive implosions are examined using a suite of nuclear and x-ray diagnostics on the OMEGA laser. The neutron-averaged hot-spot velocity (u→hs) and apparent ion temperature (Ti) asymmetry are determined from neutron time-of-flight measurements of the primary deuterium tritium fusion neutron energy spectrum, while the areal density (ρR) of the compressed fuel surrounding the hot spot is inferred from measurements of the scattered neutron energy spectrum. The low-mode perturbations of the hot-spot shape are characterized from x-ray self-emission images recorded along three quasi-orthogonal lines of sight. Implosions with significant mode-1 laser-drive asymmetries show large hot-spot velocities (>100 km/s) in a direction consistent with the hot-spot elongation observed in x-ray images, measured Ti asymmetry, and ρR asymmetry. Laser-drive corrections have been applied through shifting the initial target location in order to mitigate the observed asymmetry. With the asymmetry corrected, a more-symmetric hot spot is observed with reduced u→hs, Ti asymmetry, ρR asymmetry, and a 30% increase in the fusion yield.

Spin depolarization induced by self-generated magnetic fields during cylindrical implosions

Spin-polarized fuels are promising for inertial confinement fusion due to the enhanced fusion cross section. One significant concern of spin-polarized inertial confinement fusion is whether the nuclei polarization could survive in the implosions and contribute to ignitions. Here we present numerical simulation methods and results of spin dynamics of polarized deuterium-tritium fuels in strong self-generated magnetic fields during the implosions of dense cylindrical shells. The magnetic field generation and evolution is modeled with generalized Ohm's laws combined with hydrodynamic equations. The spin dynamics is investigated with a particle-tracking method, by solving the spin precession equations of tracked particles. Rayleigh-Taylor instabilities and Richtmyer-Meshkov instabilities are found to be the main cause of depolarization. Hydrodynamic instabilities lead to depolarization of nuclei near the hot-spot shell interface, and an asymmetric shock front leads to depolarization of nuclei inside a hot spot. Deuterium polarization is more stable than tritium polarization due to its smaller gyromagnetic ratio. Low-mode perturbations can lead to higher depolarization inside a hot spot than high-mode perturbations. In the multimode simulations, the modes around 16-32 are significant for hot-spot depolarization.

Two-dimensional thin shell model for the nonlinear Rayleigh-Taylor instability in spherical geometry

A two-dimensional thin shell model has been developed for the Rayleigh-Taylor instability in spherical geometry. The present model generalizes the work of Ott [Phys. Rev. Lett. 29, 1429 (1972)] to the case of spherical geometry. The governing equations for the shell motion and the deformation are obtained analytically and solved numerically. The linear growth rates from our model agree well with those from Mikaelian’s work but are slightly larger than those from the classical prediction for low-mode perturbations. The perturbation amplitudes and bubble velocities from our model are in general agreement with those from the weakly nonlinear model of Zhang [Phys. Plasmas 24, 062703 (2017)] and Layzer’s model [Astrophys. J. 122, 1 (1955)], respectively. Moreover, the present model can investigate the shell evolution initialized by multimode axisymmetric perturbations in spherical geometry.

Synthetic nuclear diagnostics for inferring plasma properties of inertial confinement fusion implosions

A suite of synthetic nuclear diagnostics has been developed to post-process radiation hydrodynamics simulations performed with the code Chimera. These provide experimental observables based on simulated capsule properties and are used to assess alternative experimental and data analysis techniques. These diagnostics include neutron spectroscopy, primary and scattered neutron imaging, neutron activation, γ-ray time histories and carbon γ-ray imaging. Novel features of the neutron spectrum have been analysed to infer plasma parameters. The nT and nD backscatter edges have been shown to provide a shell velocity measurement. Areal density asymmetries created by low mode perturbations have been inferred from the slope of the downscatter spectrum down to 10 MeV. Neutron activation diagnostics showed significant aliasing of high mode areal density asymmetries when observing a capsule implosion with 3D multimode perturbations applied. Carbon γ-ray imaging could be used to image the ablator at a high convergence ratio. Time histories of both the fusion and carbon γ signals showed a greater time difference between peak intensities for the perturbed case when compared to a symmetric simulation.

Weakly nonlinear incompressible Rayleigh-Taylor instability in spherical and planar geometries

The relationship between the weakly nonlinear (WN) solutions of the Rayleigh-Taylor instability in spherical geometry [Zhang et al., Phys. Plasmas 24, 062703 (2017)] and those in planar geometry [Wang et al., Phys. Plasmas 19, 112706 (2012)] is analyzed. In the high-mode perturbation limit (Pn(cos θ), n≫1), it is found that at the equator, the contributions of mode P2n along with its neighboring modes, mode P3n along with its neighboring modes, and mode Pn at the third order along with its neighboring modes are equal to those of the second harmonic, the third harmonic, and the third-order feedback to the fundamental mode, respectively, in the planar case with a perturbation of the same wave vector and amplitude as those at the equator. The trends of WN results in spherical geometry towards the corresponding planar counterparts are found, and the convergence behaviors of the neighboring modes of Pn, P2n, and P3n are analyzed. Moreover, the spectra generated from the high-mode perturbations in the WN regime are provided. For low-mode perturbations, it is found that the fundamental modes saturate at larger amplitudes than the planar result. The geometry effect makes the bubbles at or near the equator grow faster than the bubbles in planar geometry in the WN regime.

Rayleigh-Taylor instability at spherical interfaces of incompressible fluids**Project supported by the National Natural Science Foundation of China (Grant Nos. 11275031, 11475034, 11575033, 11574390, and 11274026) and the National Basic Research Program of China (Grant Nos. 2013CB834100 and 2013CBA01504).

Rayleigh-Taylor instability (RTI) of three incompressible fluids with two interfaces in spherical geometry is derived analytically. The growth rate on the two interfaces and the perturbation feedthrough coefficients between two spherical interfaces are derived. For low-mode perturbation, the feedthrough effect from outer interface to inner interface is much more severe than the corresponding planar case, while the feedback from inner interface to the outer interface is smaller than that in planar geometry. The low-mode perturbations lead to the pronounced RTI growth on the inner interface of a spherical shell that are larger than the cylindrical and planar results. It is the low-mode perturbation that results in the difference between the RTI growth in spherical and cylindrical geometry. When the mode number of the perturbation is large enough, the results in cylindrical geometry are recovered.

Self-Generated Magnetic Fields in the Stagnation Phase of Indirect-Drive Implosions on the National Ignition Facility.

Three-dimensional extended-magnetohydrodynamic simulations of the stagnation phase of inertial confinement fusion implosion experiments at the National Ignition Facility are presented, showing self-generated magnetic fields over 10^{4} T. Angular high mode-number perturbations develop large magnetic fields, but are localized to the cold, dense hot-spot surface, which is hard to magnetize. When low-mode perturbations are also present, the magnetic fields are injected into the hot core, reaching significant magnetizations, with peak local thermal conductivity reductions greater than 90%. However, Righi-Leduc heat transport effectively cools the hot spot and lowers the neutron spectra-inferred ion temperatures compared to the unmagnetized case. The Nernst effect qualitatively changes the results by demagnetizing the hot-spot core, while increasing magnetizations at the edge and near regions of large heat loss.

Weakly nonlinear Bell-Plesset effects for a uniformly converging cylinder

In this research, a weakly nonlinear (WN) model has been developed considering the growth of a small perturbation on a cylindrical interface between two incompressible fluids which is subject to arbitrary radial motion. We derive evolution equations for the perturbation amplitude up to third order, which can depict the linear growth of the fundamental mode, the generation of the second and third harmonics, and the third-order (second-order) feedback to the fundamental mode (zero-order). WN solutions are obtained for a special uniformly convergent case. WN analyses are performed to address the dependence of interface profiles, amplitudes of inward-going and outward-going parts, and saturation amplitudes of linear growth of the fundamental mode on the Atwood number, the mode number (m), and the initial perturbation. The difference of WN evolution in cylindrical geometry from that in planar geometry is discussed in some detail. It is shown that interface profiles are determined mainly by the inward and outward motions rather than bubbles and spikes. The amplitudes of inward-going and outward-going parts are strongly dependent on the Atwood number and the initial perturbation. For low-mode perturbations, the linear growth of fundamental mode cannot be saturated by the third-order feedback. For fixed Atwood numbers and initial perturbations, the linear growth of fundamental mode can be saturated with increasing m. The saturation amplitude of linear growth of the fundamental mode is typically 0.2λ–0.6λ for m < 100, with λ being the perturbation wavelength. Thus, it should be included in applications where Bell-Plesset [G. I. Bell, Los Alamos Scientific Laboratory Report No. LA-1321, 1951; M. S. Plesset, J. Appl. Phys. 25, 96 (1954)] converging geometry effects play a pivotal role, such as inertial confinement fusion implosions.

Cross-code comparisons of mixing during the implosion of dense cylindrical and spherical shells

We present simulations of the implosion of a dense shell in two-dimensional (2D) spherical and cylindrical geometry performed with four different compressible, Eulerian codes: RAGE, FLASH, CASTRO, and PPM. We follow the growth of instabilities on the inner face of the dense shell. Three codes employed Cartesian grid geometry, and one (FLASH) employed polar grid geometry. While the codes are similar, they employ different advection algorithms, limiters, adaptive mesh refinement (AMR) schemes, and interface-preservation techniques. We find that the growth rate of the instability is largely insensitive to the choice of grid geometry or other implementation details specific to an individual code, provided the grid resolution is sufficiently fine. Overall, all simulations from different codes compare very well on the fine grids for which we tested them, though they show slight differences in small-scale mixing. Simulations produced by codes that explicitly limit numerical diffusion show a smaller amount of small-scale mixing than codes that do not. This difference is most prominent for low-mode perturbations where little instability finger interaction takes place, and less prominent for high- or multi-mode simulations where a great deal of interaction takes place, though it is still present. We present RAGE and FLASH simulations to quantify the initial perturbation amplitude to wavelength ratio at which metrics of mixing agree across codes, and find that bubble/spike amplitudes are converged for low-mode and high-mode simulations in which the perturbation amplitude is more than 1% and 5% of the wavelength of the perturbation, respectively. Other metrics of small-scale mixing depend on details of multi-fluid advection and do not converge between codes for the resolutions that were accessible.

Laser-beam zooming to mitigate crossed-beam energy losses in direct-drive implosions.

Spherically symmetric direct-drive-ignition designs driven by laser beams with a focal-spot size nearly equal to the target diameter suffer from energy losses due to crossed-beam energy transfer (CBET). Significant reduction of CBET and improvements in implosion hydrodynamic efficiency can be achieved by reducing the beam diameter. Narrow beams increase low-mode perturbations of the targets because of decreased illumination uniformity that degrades implosion performance. Initiating an implosion with nominal beams (equal in size to the target diameter) and reducing the beam diameter by ∼ 30%-40% after developing a sufficiently thick target corona, which smooths the perturbations, mitigate CBET while maintaining low-mode target uniformity in ignition designs with a fusion gain ≫ 1.

Shock ignition target design for inertial fusion energy

Continuing work in the design of shock ignition targets is described. Because of reduced implosion velocity requirements, low target adiabats, and efficient drive by short wavelength lasers, these targets produce high gain (>100) at laser energies well below 1 MJ. Effects of hydrodynamic instabilities such as Rayleigh–Taylor or Richtmyer–Meshkov are greatly reduced in these low-aspect ratio targets. Of particular interest is the optimum ratio of ignitor to compression pulse energy. A simple pellet model and simulation-derived coupling coefficients are used to analyze optimal fuel assembly, and determine that shock ignition allows enough control to create theoretically optimum assemblies. The effects on target design due to constraints on the compression and ignitor pulse intensities are also considered and addressed. Significant sensitivity is observed from low-mode perturbations because of large convergence ratios, but a more powerful ignitor can mitigate this.