High-density carbon ablator ignition path with low-density gas-filled rugby hohlraum

Transcription

1 High-density carbon ablator ignition path with low-density gas-filled rugby hohlraum Peter Amendt, Darwin D. Ho and Ogden S. Jones Lawrence Livermore National Laboratory, Livermore CA A recent low gas-fill density (0.6 mg/cc 4 He) cylindrical hohlraum experiment on the National Ignition Facility has shown high laser-coupling efficiency (>96%), reduced phenomenological laser drive corrections, and improved high-density carbon capsule implosion symmetry [Jones et al., BAPS 59 (15), 66 (2014)]. In this Letter an ignition design using a large rugby-shaped hohlraum [Amendt et al., Phys. Plasmas (2014)] for high energetics efficiency and symmetry control with the same low gas-fill density (0.6 mg/cc 4 He) is developed as a potentially robust platform for demonstrating thermonuclear burn. The companion high-density carbon capsule for this hohlraum design is driven by an adiabat-shaped [Betti et al., Phys. Plasmas 9, 2277 (2002)] 4-shock drive profile for robust high gain (>10) 1-D ignition performance and large margin to 2-D perturbation growth. Efforts to demonstrate thermonuclear energy breakeven with indirect drive on the National Ignition Facility (NIF) continue in earnest [1]. Indirect drive uses a laser-heated hohlraum (or high-z enclosure) for generating soft x rays that irradiate and implode a centrally located capsule filled with deuterium (D) and tritium (T) fuel. The first ignition attempts used a 4- shock low-adiabat drive in cylinder hohlraums at moderate gas-fill density (ρ=0.96 mg/cc) and large outer-to-inner laser cone wavelength separation to control time-integrated x-ray drive symmetry, leading to good compression [2] but lower neutron yield than expected [3]. Significant fuel self-heating has been recently demonstrated for a higher adiabat 3-shock high-foot laser power drive history in cylinder hohlraums at higher gas-fill density (ρ=1.6 mg/cc) and still large outer-to-inner laser cone wavelength separation Δλ>5Å [4]. The high gas fill is used to constrain the inward motion of the ablating Au or U hohlraum wall plasma, which can lead to large x-ray drive symmetry swings in time. These symmetry variations result in residual (non-radial) kinetic energy at peak compression [5] and degraded implosion yield. A large wavelength separation (~8Å) is required to compensate for inhibited inner-cone laser propagation at late time and to achieve time-integrated 2 nd Legendre mode (P2) drive symmetry. The combination of high gas-fill density and large cross-beam energy transfer (CBET)[6] from the imposed large Δλ provides for a computationally challenging hohlraum environment, involving 3-D effects and significant plasma-induced laser backscatter (~40-50%) on the inner cones. An additional feature is the prevalence of 20-25% unaccounted drive energy according to high-flux model [7] comparisons with the database that collectively encompasses (1) shock-tuning experiments [8], (2) 1-D streaked radiography experiments for direct measurement of the capsule implosion trajectory [9], and (3) observed x-ray and neutron peak emission bang times [10]. The origin of this missing energy is not known, possibly involving mixing of the high-z hohlraum material with the low-z gas fill [11] and the role of numerical mesh management [12] in mainline radiation-hydrodynamic (RH) simulations [13]. At the other end of hohlraum gas-fill studies is the near-vacuum hohlraum (NVH) platform (ρ=0.032 mg/cc of 4 He) that was successfully used to drive exploding pusher or thin plastic (CH) ablators [14]. Due to the short laser pulse length used (~4 ns), good agreement with RH simulations was achieved without the need for phenomenological drive multipliers [10,15], i.e., no missing energy was inferred. These results have served as the impetus for recent 2- and 3- shock high-density carbon (HDC) ablator studies on the NIF [16, 17]. The strategy is to combine the intrinsically short laser power histories of HDC with the NVH platform for acceptable drive symmetry control without depending on large CBET at high Δλ. A HDC capsule requires an early-time high foot in laser drive to ensure melting of the ablator after first shock passage, thereby avoiding the crystalline diamond phase that can seed hydrodynamic instability growth [16]. To date HDC implosion experiments in NVHs have shown a strong bias towards equator high implosion drive asymmetry, in contrast to RH simulation results [18]. The origin of this difference is

2 not firmly established, but may involve ion diffusive effects from interpenetrating low-z HDC ablator and high-z hohlraum ions that are not captured in mainline fluid-based modeling tools. According to RH simulations of these NVHs the stopping length for a fully ionized ablated carbon ion (~ 350 µm/ns) into the expanding gold wall (~ -250 µm/ns) is on the order of 500 µm, based on a Au mass density of ~1 mg/cc and temperature of ~3 kev just before collision onset. Such long mean free paths challenge the use of RH simulations in a physical regime where kinetic effects are important [19]. For example, ion mixing from kinetic processes may modify the high-z density and temperature distribution above the capsule compared with RH simulations and allow for more efficient inner-cone laser propagation after the stagnation time of C/Au because of reduced inverse bremsstrahlung absorption or enhanced inner-beam refraction. A successful HDC ignition campaign based on the NVHs will require a predictive capability for informing studies of high-convergence implosion symmetry. Low gas-fill densities are attracting interest for potentially simplified hohlraum dynamics: high laserhohlraum coupling efficiency through minimal backscatter and unaccounted energy, reduced need for CBET, and calculable implosion symmetry. In this Letter recent experimental results with a 0.6 mg/cc 4 He gas fill that achieved close to optimal cylinder hohlraum performance in driving a 2-shock HDC capsule form the basis for a novel ignition path. The associated design uses a 4-shock drive profile with a comparatively low fuel adiabat (for high ignition margin) and applies the technique of adiabat shaping [20, 21] for ablation front instability control. The capsule is driven by an oversized rugby-shaped gold hohlraum [22] for high drive efficiency and robust symmetry control that is filled with the same 0.6 mg/cc 4 He gas. The 4-shock HDC ablator design has high 1-D simulated ignition margin (for achieving >1 MJ yield) and good tolerance to ablator and ice surface roughness, according to capsule-only high-resolution 2-D RH simulations. Hohlraum simulations that integrate laser propagation, x-ray transport and capsule/hohlraum hydrodynamics give a 2-D yield of nearly 15 MJ for only ~1.5 MJ absorbed laser energy. Based on demonstrated efficient hohlraum performance at low gas-fill densities and RH simulations that suggest strong 1- and 2-D 4-shock HDC implosion performance, an ignition design leveraging this unexplored parameter space is described and proposed as an alternate path for demonstrating ignition on the NIF. The recent low gas-fill density target that was fielded on the NIF combined the benefits of high efficiency of NVHs and the inherent symmetry control of gas-filled hohlraums for the purpose of driving a 2- shock HDC design over a ~7 ns laser pulse length [23]. The recognized challenges with the 2-shock ignition design are suitable control of implosion symmetry and the low 1-D ignition margin due to the intrinsically high fuel adiabat. Late time symmetry control is challenging because of significant wall motion and the ensuing laser absorption close to the hohlraum symmetry axis. To reduce the risk of excessive hohlraum filling with high-z wall material, the low gas-filled cylindrical hohlraum was enlarged to 6.72 mm diameter from the standard 5.75 mm. The larger hohlraum results in greater geometric smoothing of intrinsic hohlraum radiation modes at the capsule ablation front and ensures greater laser clearance at the laser entrance holes (LEHs). The incident laser energy was 1.5 MJ with a peak power of 370 TW, and no wavelength separation (Δλ=0) was chosen. The demonstrated backscatter was <4%, the effective x-ray drive was inferred to lie much closer to high-flux model predictions for typical gas-filled hohlraums (ρ= mg/cc 4 He), and improved imploded core symmetry in good agreement with modeling was observed [23]. Remarkably, increasing the gas fill for HDC implosions did not incur a drive penalty compared with NVHs. The early experimental success with a low gas-fill hohlraum platform opens up a promising ignition path for HDC ablators that is now described. The three main design challenges are ensuring (1) high 1-D performance margin, (2) robustness to instability growth at the ablation front and fuel-ablator and DT gas-ice interfaces, and (3) favorable hohlraum energetics with robust symmetry and preheat control. The demonstrated favorable symmetry control and energetics with a low gas-filled hohlraum density mean that 4-shock capsule designs with their ignitionrelevant low 1-D fuel adiabats α (or ratio of pressure to Fermi degenerate pressure) now become viable. Reference 17 compared the relevant 1-D implosion properties between 2- and 4-shock HDC designs, but without the potential benefit of adiabat shaping from laser pickets [20, 21]. Both capsule designs have comparable peak implosion speeds υ of ~380 µm/ns, but with a significantly higher adiabat for the 2-shock design (3.7 vs 1.5 for the 4-shock design). The ignition threshold factor ITF=M+1, is related to the ignition

3 Yield P stag υ α At GF [MJ] [Gbar] [µm/ns] Table 1: Calculated 1-D capsule yield, peak fuel stagnation pressure (with burn turned off), peak implosion speed υ, maximum fuel adiabat α, Atwood number of DT-ice/HDC interface At, and associated peak linear growth factor (GF) at peak implosion speed (l 80) for HDC 4-shock capsule [see Fig. 1a]. performance margin M and scales as α 2.6 υ 8 [24], giving more than a 10 higher margin for the 4-shock design. The current 2-shock design in a cylinder NVH uses undoped HDC [18], allowing for larger ablationfront stabilization of Rayleigh-Taylor instability growth but low 1-D margin. Table 1 lists the key properties of the 4-shock design used in the integrated hohlraum tune [described below]. However, adiabat shaping of the drive history may significantly reduce the exposure of the 4-shock target to ablation front instability growth. Adiabat shaping refers to the launch of a decaying first shock with a strong laser picket [see Figs. 2a, 3a] to potentially raise the adiabat of the outer ablator for ablation front stability control while keeping the fuel on a low adiabat. The shock pressure is designed to remain above 7 Mbars at the ablator-fuel interface, minimizing the risk of instability seeding from multi-phase HDC. The linear growth factor ~1200 for the 4-shock HDC design with adiabat shaping is 2 lower compared with the 4-shock lowadiabat CH design [17]. Figure 1b shows the simulated perturbation growth for the 4-shock design at the time of maximum shell kinetic energy using the ( 1 ) Fig. 2a-b: (a) Total requested laser power for 4-shock HDC ignition hohlraum design; (b) requested laser cone fraction with 300 ps outer cone delay at t=0. measured surface roughness spectra on all material interfaces from an as-built (undoped) HDC capsule. Up to ~2.5 the measured roughness on all surfaces is predicted to still lead to ignition, according to RH simulation parameter scans. The simulated clean fuel fraction defined as the fraction of DT fuel with less than 5.0 at.% of carbon at the time of maximum shell kinetic energy is a favorable 80% and lies above the threshold for onset of important 3-D effects. The tradeoff between capsule robustness and hohlraum performance is best illustrated by considering the drive duration. Fuel adiabat control benefits from longer pulse lengths [24], but hohlraum symmetry control and coupling efficiency (through reduced laser backscatter and x-ray drive corrections) favor shorter pulses. In order to accommodate the ~9 ns drive history of the optimized 4-shock capsule tune, added margin is sought by using a rugby-shaped Au hohlraum with a 7 mm diameter, 1.09 cm length and a comparatively large LEH fraction of ~60% [See Fig. 4]. The use of a rugby hohlraum over a cylinder allows for a greater volume above the capsule while preserving the hohlraum drive efficiency through a nearly fixed wall surface area. Added benefits of the rugby hohlraum include potentially benign (gold) M- Fig. 1a-b: (a) Diagram of HDC 4-shock design with indicated materials, dopant and dimensions; (b) 2-D highresolution (l = ) multi-mode capsule-only simulation at time of peak implosion speed with left (right) panel showing materials (density). Fig. 3a-b: (a) Uncompensated (for LEH closure) Dante radiation temperature T rad versus time with a viewing angle of 37 to the hohlraum symmetry axis; (b) Dante gold M- band (>1.8 kev) fraction of spectral flux versus Dante T rad for 4-shock ignition HDC design.

4 Fig. 4: Rugby-shaped (quarter) hohlraum geometry and laser ray geometry at peak power (6.5 ns; see Fig. 2a) for HDC 4- shock design (capsule not shown). Ray bundles as shown include only rays with greater than 50% of the peak (central) intensity for clarity. band x-ray preheat [See Fig. 3b] from reduced laser intensity at the curved wall, less hot electron generation, greatly reduced time-dependent drive symmetry variations in time, reduced exposure to specular laser glint and a reduced gold bubble feature in the wall blow-off that is typically seen in cylinder hohlraum experiments [25] and simulations [11]. The peak M-band fraction shown in Fig. 3b is predicted to remain below 20%, which compares favorably with the measured and simulated ~21% level in NVHs. A lower level of M-band preheat translates into smaller dopant concentrations that are needed in the ablator for higher ablation velocities and, in turn, higher performance and stability margins. The gold bubble may lead to impaired inner beam propagation at late time, but the curvature of a rugby-shaped wall distributes the incident laser energy over a larger area for a reduced risk of non-uniform wall ablation. Figure 4 shows the laser propagation at peak laser power, indicating robust inner cone (23 and 30 ) propagation past the incipient (small) gold bubble produced by the outer cones (44 and 50 ). However, the NIF rugby hohlraum database for 4-shock low-adiabat CH implosions (~20 ns laser pulse length) has shown some symmetry anomalies compared with RH simulations using the high-flux model - unless the outer laser cones are significantly displaced inward toward the capsule [11]. The origin of this anomaly is not established, but may involve hydrodynamic mix or ion diffusion across the gold-helium interface near the LEH. To minimize this risk (beyond the benefit of a ~2 shorter laser pulse with HDC), a comparatively larger LEH (~60% vs 50.6% for low-foot rugby hohlraums with 1.2 mg/cc 4 He fill [11]) is used in the 4-shock design [See Fig. 4]. Figure 2a-b depicts the requested laser power history and cone fraction. Use of an early time high (inner) cone fraction is intended to appreciably ablate the Fig. 5a-b: (a) Normalized 2 nd Legendre coefficient of simulated instantaneous ablation pressure asymmetry for HDC 4-shock capsule in rugby hohlraum (green) and HDC 2-shock in cylinder NVH [shot N140702] (blue) versus time; (b) fuel density map before (~9 kev) ignition onset (~12 kev) with rugby hohlraum symmetry axis along vertical. P2/P0 and P4/P0 distortion of maximum fuel density [ρ~900 g/cc] contour ~16% and 1.2%. window before onset of the ~70 TW picket in order to suppress hot electron generation by the two-plasmon decay instability. There is the potential for generating (non-azimuthal) CBET during the picket from the expanding LEH window even for low Δλ. However, use of just the 44 cones to drive the picket pulse can be shown to remove this risk. By removing CBET on the picket the usual hohlraum modeling uncertainties arising from choice of density saturation parameter are bypassed. There is also some CBET near peak laser power at Δλ=0, but a nominally small value of Δλ can largely suppress this contribution to hohlraum symmetry. Figure 5a shows the associated simulated capsule ablation pressure asymmetry history (P2/P0), indicating very low symmetry swings on the order of ±2% compared with the 2-shock HDC capsule in cylinder NVHs (~ ±5%). For low-adiabat implosions in rugby hohlraums, the relatively low level of P2 symmetry swings compared with cylinders has been experimentally confirmed [11]. The key figure-ofmerit is the degree of symmetry control at the end of the laser pulse. Effective inner-cone laser propagation at late time is predicted because of the large volume between the capsule and rugby equator and the large LEH fraction. Figure 5b shows the imploded fuel just before ignition, suggesting keen symmetry control without relying on late time CBET from large Δλ as in higher density gas-filled hohlraums. The energy output is over 15 MJ, corresponding to a gain of ~10. The low gas-fill hohlraum design also benefits from remaining in the fluid regime, which facilitates modeling with mainline RH simulation methods. The

5 carbon stopping length in the helium gas fill is on the order of 20 microns, largely preventing the carbon blow off from penetrating the Au plasma. In addition, the low gas fill density delays deceleration onset of the classical Au/ 4 He interface compared with high gas-fill platforms [11], lessening the amount of Rayleigh- Taylor instability growth after (wall/gas) stagnation onset. The alternate ignition path described herein is the culmination of promising early data with low-density gas-filled (cylindrical) hohlraums [23], recent capsule designs that both exploit adiabat shaping for improved stability control and optimize performance margin [21], and rugby-shaped hohlraum simulations for improved efficiency and symmetry control. However, the path to ignition is justifiably fraught with uncertainty over the fidelity of the mainline simulation tools, thereby requiring a careful and long-term experimental plan to benchmark the calculations and iterate on the design as necessary. In summary, a low density gas-filled hohlraum regime shows promise towards a novel ignition path that combines the energetics benefits of a NVH and the symmetry control advantages of a gas-filled hohlraum. A 4-shock HDC capsule design uses an oversized rugby-shaped hohlraum filled with the same 4 He density for greater margin to impaired laser propagation and potential symmetry anomalies. Further energetics margins can be obtained with depleted uranium hohlraums and the associated ~6% improvement in energy coupling [26]. Useful input from O.L. Landen and J. Milovich is gratefully acknowledged. The suggestions of L. Berzak-Hopkins, S. Khan, S. Nagel, R.J. Rygg and D. Turnbull are appreciated. Work performed under the auspices of Lawrence Livermore National Security, LLC (LLNS) under Contract DE-AC52-07NA2734. [1] J.D. Lindl, P. Amendt, R.L. Berger, S.G. Glendinning, S.H. Glenzer, S.W. Haan, R.L. Kauffman, O.L. Landen, and L.J. Suter, Phys. Plasmas 21, (2004). [2]H.F. Robey, B.J.MacGowan, O.L. Landen, K.N. LaFortune, C. Widmayer, P.M. Celliers, J.D. Moody, J.S. Ross, J. Ralph, S. LePape et al., Phys. Plasmas 20, (2013). [3] M.J. Edwards, P.K. Patel, J.D. Lindl, B.K. Spears, S.V. Weber, L.J. Atherton, D.L. Bleuel, D.K. Bradley, D.A. Callahan, C.J. Cerjan et al., Phys. Plasmas 20, (2013). [4] O.A. Hurricane, D.A. Callahan, D.T. Casey, P.M. Celliers, C. Cerjan, E.L. Dewald, T.R. Dittrich, T. Doppner, D.E. Hinkel, L.F. Berzak Hopkins et al., Nature 506, 343 (2014). [5] R.P.J. Town, D.K. Bradley, A. Kritcher, O.S. Jones, J.R. Rygg, R. Tommasini, M. Barrios, L.R. Benedetti, L.F. Berzak Hopkins, P.M. Celliers et al., Phys. Plasmas 21, (2014); A.L. Kritcher, R. Town, D. Bradley D. Clark, B. Spears, O. Jones, S. Haan, P.T. Springer, J. Lindl, R.H.H. Scott, et al., Phys. Plasmas 21, (2014). [6] P. Michel, W. Rozmus, E.A. Williams, L. Divol, R.L. Berger, S.H. Glenzer and D.A. Callahan, Phys. Plasmas 20, (2013). [7] M.D. Rosen, H.A. Scott, D.E. Hinkel, E.A. Williams, D.A. Callahan, R.P.J. Town, L. Divol, P.A. Michel, W.L. Kruer, L.J. Suter et al., High Energy Density Physics 7, 180 (2011). [8]H.F. Robey, T.R. Boehly, P.M. Celliers, J.H. Eggert, D. Hicks, R.F. Smith, R. Collins, M.W. Bowers, K.G. Krauter, P.S. Datte et al., Phys. Plasmas 19, (2012). [9] D.G. Hicks, N.B. Meezan, E.L. Dewald, A.J. MacKinnon, R.E. Olson, D.A. Callahan, T. Doppner, L.R. Benedetti, D.K. Bradley, P.M. Celliers et al., Phys. Plasmas 19, (2012). [10] D.S. Clark, D.E. Hinkel, D.C. Eder, O.S. Jones, S.W. Haan, B.A. Hammel, M.M. Marinak, J.L. Milovich, H.F. Robey, L.J. Suter and R.P.J. Town, Phys. Plasmas 20, (2013). [11] P. Amendt, J.S. Ross, J.L. Milovich, M. Schneider, E. Storm, D.A. Callahan, D. Hinkel, B. Lasinski, D. Meeker, P. Michel, J. Moody and D. Strozzi, Phys. Plasmas 21, (2014). [12] C.A. Thomas (private communication, 2014). [13] G.B. Zimmerman and W.L. Kruer, Comments Plasma Phys. Controlled Fusion 2, 51 (1975); M.M. Marinak, R.E. Tipton, O.L. Landen, T.J. Murphy, P. Amendt, S.W. Haan, S.P. Hatchett, C.J. Keane, R. McEachern and R. Wallace, Phys. Plasmas 3, 2070 (1996). [14] S. Le Pape, L. Divol, L. Berzak Hopkins, A. MacKinnon, N.B. Meezan, D. Casey, J. Frenje, H. Herrmann, J. McNaney, T. Ma et al., Phys. Rev. Lett. 112, (2014). [15] O.S. Jones, C.J. Cerjan, M.M. Marinak, J.L. Milovich, H.F. Robey, P.T. Springer, L.R. Benedetti, D.L. Bleuel, E.J. Bond, D.K. Bradley et al., Phys. Plasmas 19, (2012). [16] D. Ho, Bull. Am. Phys. Soc. 52(16), 273 (2007). [17] A.J. MacKinnon, N.B. Meezan, J.S. Ross, S. Le Pape, L. Berzak Hopkins, L. Divol, D. Ho, J. Milovich, A. Pak, J. Ralph et al., Phys. Plasmas 21, (2014). [18] L. Berzak-Hopkins, S. Le Pape, L. Divol, N. Meezan, A. MacKinnon, D. Ho, O. Jones, S. Khan, J. Milovich, J. Ross, P. Amendt, D. Casey, A. Pak, J. Peterson, J. Ralph and J. Rygg, Near-vacuum hohlraums for driving fusion

6 implosions with high density carbon ablators, Phys. Plasmas (submitted). [19] M.J. Rosenberg, H.G. Rinderknecht, N.M. Hoffman, P.A. Amendt, S. Atzeni, A.B. Zylstra, C.K. Li, F.H. Seguin, H. Sio, M. Gatu Johnson et al., Phys. Rev. Lett. 112, (2014). [20] R. Betti, K. Anderson, V.N. Goncharov, R.L. McCrory, D.D. Meyerhofer, S. Skupsky and R.P.J. Town, Phys. Plasmas 9, 2277 (2002). [21] D. Ho, P. Amendt, D. Clark, S. Haan, J. Milovich, J. Salmonson, G. Zimmerman, L. Berzak Hopkins, J. Biener, N. Meezan et al., Bull. Am. Phys. Soc. 59(15), 199 (2014). [22] P. Amendt, C. Cerjan, A. Hamza, D.E. Hinkel, J.L. Milovich and H.F. Robey, Phys. Plasmas 14, (2007); M. Vandenboomgaerde, J. Bastian, A. Casner, D. Galmiche, J.-P. Jadaud, S. Lafitte, S. Liberatore, G. Malinie and F. Philippe, Phys. Rev. Lett. 99, (2007). [23] O. Jones, N. Izumi, L.Berzak-Hopkins, D.J. Strozzi, P.A. Amendt, G.N. Hall, D.D. Ho, S.F. Khan, N.B. Meezan, J.D. Moody et al., Hohlraum fill gas density scaling of x-ray drive, symmetry and laser coupling backscatter in 6.72-mm NIF hohlraums, Bull. Am. Phys. Soc. 59(15), 66 (2014). [24] J. Lindl, O. Landen, J. Edwards and E. Moses, Phys. Plasmas 21, (2014). [25] S.A. MacLaren, M.B. Schneider, K. Widmann, J.H. Hammer, B.E. Yoxall, J.D. Moody, P.M. Bell, L.R. Benedetti, D.K. Bradley, M.J. Edwards et al., Phys. Rev. Lett. 112, (2014). [26] D.A. Callahan, N.B. Meezan, S.H. Glenzer, A.J. MacKinnon, L.R. Benedetti, D.K. Bradley, J.R. Celeste, P.M. Celliers, S.N. Dixit, T. Doppner et al., Phys. Plasmas 19, (2012).

7

8

9

10

11