A multi-frame shadowgraphy diagnostic has been developed and applied to laser preheat experiments relevant to the Magnetized Liner Inertial Fusion (MagLIF) concept. The diagnostic views the plasma created by laser preheat in MagLIF-relevant gas cells immediately after the laser deposits energy as well as the resulting blast wave evolution later in time. The expansion of the blast wave is modeled with 1D radiation-hydrodynamic simulations that relate the boundary of the blast wave at a given time to the energy deposited into the fuel. This technique is applied to four different preheat protocols that have been used in integrated MagLIF experiments to infer the amount of energy deposited by the laser into the fuel. The results of the integrated MagLIF experiments are compared with those of two-dimensional LASNEX simulations. The best performing shots returned neutron yields ∼40-55% of the simulated predictions for three different preheat protocols.
A series of Magnetized Liner Inertial Fusion (MagLIF) experiments have been conducted in order to investigate the mix introduced from various target surfaces during the laser preheat stage. The material mixing was measured spectroscopically for a variety of preheat protocols by employing mid-atomic number surface coatings applied to different regions of the MagLIF target. The data show that the material from the top cushion region of the target can be mixed into the fuel during preheat. For some preheat protocols, our experiments show that the laser-entrance-hole (LEH) foil used to contain the fuel can be transported into the fuel a significant fraction of the stagnation length and degrade the target performance. Preheat protocols using pulse shapes of a few-ns duration result in the observable LEH foil mix both with and without phase-plate beam smoothing. In order to reduce this material mixing, a new capability was developed to allow for a low energy (∼20 J) laser pre-pulse to be delivered early in time (-20 ns) before the main laser pulse (∼1.5 kJ). In experiments, this preheat protocol showed no indications of the LEH foil mix. The experimental results are broadly in agreement with pre-shot two-dimensional HYDRA simulations that helped motivate the development of the early pre-pulse capability.
The Magnetized Liner Inertial Fusion concept (MagLIF) [Slutz et al., Phys. Plasmas 17, 056303 (2010)] is being studied on the Z facility at Sandia National Laboratories. Neutron yields greater than 1012 have been achieved with a drive current in the range of 17-18 MA and pure deuterium fuel [Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)]. We show that 2D simulated yields are about twice the best yields obtained on Z and that a likely cause of this difference is the mix of material into the fuel. Mitigation strategies are presented. Previous numerical studies indicate that much larger yields (10-1000 MJ) should be possible with pulsed power machines producing larger drive currents (45-60 MA) than can be produced by the Z machine [Slutz et al., Phys. Plasmas 23, 022702 (2016)]. To test the accuracy of these 2D simulations, we present modifications to MagLIF experiments using the existing Z facility, for which 2D simulations predict a 100-fold enhancement of MagLIF fusion yields and considerable increases in burn temperatures. Experimental verification of these predictions would increase the credibility of predictions at higher drive currents.
Material equation-of-state (EOS) models, generally providing the pressure and internal energy for a given density and temperature, are required to close the equations of hydrodynamics. As a result they are an essential piece of physics used to simulate inertial confinement fusion (ICF) implosions. Historically, EOS models based on different physical/chemical pictures of matter have been developed for ICF relevant materials such as the deuterium (D2) or deuterium-tritium (DT) fuel, as well as candidate ablator materials such as polystyrene (CH), glow-discharge polymer (GDP), beryllium (Be), carbon (C), and boron carbide (B4C). The accuracy of these EOS models can directly affect the reliability of ICF target design and understanding, as shock timing and material compressibility are essentially determined by what EOS models are used in ICF simulations. Systematic comparisons of current EOS models, benchmarking with experiments, not only help us to understand what the model differences are and why they occur, but also to identify the state-of-the-art EOS models for ICF target designers to use. For this purpose, the first Equation-of-State Workshop, supported by the US Department of Energy's ICF program, was held at the Laboratory for Laser Energetics (LLE), University of Rochester on 31 May–2nd June, 2017. This paper presents a detailed review on the findings from this workshop: (1) 5–10% model-model variations exist throughout the relevant parameter space, and can be much larger in regions where ionization and dissociation are occurring, (2) the D2 EOS is particularly uncertain, with no single model able to match the available experimental data, and this drives similar uncertainties in the CH EOS, and (3) new experimental capabilities such as Hugoniot measurements around 100 Mbar and high-quality temperature measurements are essential to reducing EOS uncertainty.
Atomic structure of N-electron atoms is often determined by solving the Hartree-Fock equations, which are a set of integro-differential equations. The integral part of the Hartree-Fock equations treats electron exchange, but the Hartree-Fock equations are not often treated as an integro-differential equation. The exchange term is often approximated as an inhomogeneous or an effective potential so that the Hartree-Fock equations become a set of ordinary differential equations (which can be solved using the usual shooting methods). Because the Hartree-Fock equations are an iterative-refinement method, the inhomogeneous term relies on the previous guess of the wavefunction. In addition, there are numerical complications associated with solving inhomogeneous differential equations. This work uses matrix methods to solve the Hartree-Fock equations as an integro-differential equation. It is well known that a derivative operator can be expressed as a matrix made of finite-difference coefficients; energy eigenvalues and eigenvectors can be obtained by using linear-algebra packages. The integral (exchange) part of the Hartree-Fock equation can be approximated as a sum and written as a matrix. The Hartree-Fock equations can be solved as a matrix that is the sum of the differential and integral matrices. We compare calculations using this method against experiment and standard atomic structure calculations. This matrix method can also be used to solve for free-electron wavefunctions, thus improving how the atoms and free electrons interact. This technique is important for spectral line broadening in two ways: it improves the atomic structure calculations, and it improves the motion of the plasma electrons that collide with the atom.
