Interaction physics for the shock ignition scheme of inertial ... - CiteSeerX

Nov 14, 2011 - the time of maximum implosion to produce a lateral hot spot where fusion ... absorption (inverse bremsstrahlung) of the laser electromagnetic ...
491KB taille 0 téléchargements 374 vues
IOP PUBLISHING

PLASMA PHYSICS AND CONTROLLED FUSION

Plasma Phys. Control. Fusion 53 (2011) 124034 (10pp)

doi:10.1088/0741-3335/53/12/124034

Interaction physics for the shock ignition scheme of inertial confinement fusion targets S Depierreux1 , C Goyon1 , K Lewis2 , H Bandulet2 , D T Michel 1,2 , G Loisel2 , V Yahia1,2 , V Tassin1 , C Stenz3 , N G Borisenko4 , W Nazarov5 , J Limpouch6 , P E Masson Laborde1 , P Loiseau1 , M Casanova1 , 7 ¨ Ph Nicola¨ı3 , S Huller , D Pesme7 , C Riconda8 , V T Tikhonchuk3 and C Labaune2 1

CEA, DAM, DIF, F-91297 Arpajon, France Laboratoire pour l’Utilisation des Lasers Intenses, Ecole Polytechnique, Palaiseau, France 3 Centre Lasers Intenses et Applications, Universit´ e Bordeaux 1, CEA, CNRS, Talence, France 4 P. N. Lebedev Physical Institute, Russian Academy of Sciences, Moscow, Russia 5 University of St Andrews, School of Chemistry, Fife KY16 9ST, Scotland, UK 6 FNSPE, Czech Technical University in Prague, 115 19 Prague 1, Czech Republic 7 Centre de Physique Th´ eorique, Ecole Polytechnique, 91128 Palaiseau cedex, France 8 Laboratoire pour l’Utilisation des Lasers Intenses, (TIPS) Universit´ e Paris 6, CNRS, CEA, Ecole Polytechnique, 94200 Ivry-sur-Seine, France 2

Received 24 June 2011, in final form 21 September 2011 Published 14 November 2011 Online at stacks.iop.org/PPCF/53/124034 Abstract This paper presents an analysis of laser–plasma interaction risks of the shock ignition (SI) scheme and experimental results under conditions relevant to the corona of a compressed target. Experiments are performed on the LIL facility at the 10 kJ level, on the LULI 2000 facility with two beams at the kJ level and on the LULI 6-beam facility with 100 J in each beam. Different aspects of the interaction of the SI pulse are studied exploiting either the flexibility of the LULI 6-beam facility to produce a very high intensity pulse or the high energy of the LIL to produce long and hot plasmas. A continuity is found allowing us to draw some conclusions regarding the coupling quality and efficiency of the SI spike pulse. It is shown that the propagation of the SI beams in the underdense plasma present in the corona of inertial confinement fusion targets could strongly modify the initial spot size of the beam through filamentation. Detailed experimental studies of the growth and saturation of backscattering instabilities in these plasmas indicate that significant levels of stimulated scattering reflectivities (larger than 40%) may be reached at least for some time during the SI pulse. (Some figures may appear in colour only in the online journal)

0741-3335/11/124034+10$33.00

© 2011 IOP Publishing Ltd

Printed in the UK & the USA

1

Plasma Phys. Control. Fusion 53 (2011) 124034

S Depierreux et al

1. Introduction Laser-driven inertial confinement fusion (ICF) relies on the use of laser beams to compress and ignite a thermonuclear fuel with the ultimate goal of producing energy [1]. Two main irradiation schemes have been proposed for ICF. The simplest one, called direct drive, consists in direct focusing of the laser beams on the fuel capsule. Due to the high illumination uniformity required to spherically compress the capsule, which is difficult to obtain with high-energy laser beams, a second scheme has been proposed. It uses a high-Z hohlraum to first convert the laser energy into x-rays, which then implode the fuel capsule. In the prospect of reducing the laser driver energy, alternative approaches have been proposed that separate the compression and ignition phases using two independent sets of laser pulses. The fast ignition (FI) scheme is based on short and intense laser pulses that are fired at the time of maximum implosion to produce a lateral hot spot where fusion reactions start. The analysis of the hot spot heating during the fuel stagnation time shows that very high intensities (∼1020 –1021 W cm−2 ) and very short (∼10–20 ps) pulses are needed. Another approach, which was proposed in the early times of ICF [2], is based on central ignition triggered by a strong shock generated at the end of the compression pulse. A dedicated laser pulse separated from the compression pulse is used to generate a strong convergent shock which is timed to collide with the return shock. The collision of the two shocks results in a region of high pressure which finally compresses a small part of the fuel up to the high pressure required for ignition. The initial thermonuclear fusion reactions then propagate to the surrounding cold high-density fuel. Recent analyses [3, 4] have demonstrated the efficiency of this scheme leading to significant fusion gains with laser energies lower than those required in conventional approaches. As the shock ignition (SI) scheme requires laser beams that directly hit the fuel capsule, the numerical simulations have been conducted in the direct drive geometry. The spike pulse intensity must be in the range of a few 1015 to a few 1016 W cm−2 and its duration is of the order of a few hundreds of picoseconds [4]. This pulse is fired when the corona surrounding the compressed fuel extends already on a millimetre distance containing a hot (∼2–5 keV) underdense plasma. The first analyses of the SI scheme [3, 4] have considered only collisional absorption of the SI pulse in this plasma ignoring excitation of parametric instabilities as there is no reliable model to completely describe their non-linear behaviour. The first experiments were conducted at the OMEGA Laser Facility [5] showing that backscattering is dominated by stimulated Raman scattering (SRS). Understanding the development of parametric instabilities of the SI pulse is of interest not only because it determines the efficiency of the coupling of laser energy to the plasma but also because it determines the way this energy is absorbed. In this paper, we present experimental data collected on different laser systems that are relevant to the interaction physics of the shock pulse with the corona of a compressed ICF target. A long-term program on laser–plasma interaction physics is pursued, associating theoretical, numerical and experimental studies, to understand the propagation and the coupling of laser pulses in the plasma for a wide range of parameters. This physics is very important for the success of implosion and ignition experiments where high efficiency and high quality of the coupling between the laser pulses and the target are needed. Apart from collisional absorption (inverse bremsstrahlung) of the laser electromagnetic wave, the laser beam can couple with the plasma modes to produce other plasma or electromagnetic waves [6]. These mechanisms develop as soon as the laser intensity exceeds a threshold that depends on the plasma and laser parameters. The main objective of SI interaction control is to ensure a highefficiency coupling with the plasma while avoiding all possible losses through parametric instabilities. 2

Plasma Phys. Control. Fusion 53 (2011) 124034

S Depierreux et al

The concern for the SI pulse mainly comes from the high intensity required to drive the high-amplitude shock at the end of the compression phase when the plasma has expanded for more than 10 ns. At this time, the SI pulse interacts with the underdense plasma heated to electron temperatures of 2–4 keV with an exponential density profile (∼ exp(−z/Ln )) along the propagation axis z of characteristic scale length Ln larger than 300 µm. The production of such plasma conditions requires multi-kilojoule lasers to maintain the high temperature during the plasma expansion. This constraint is in conflict with the necessity of producing a highintensity beam, which is usually achieved on small-scale facilities through careful wavefront corrections. Complementary measurements are thus needed using different systems. In the following, we first discuss the main issues of laser–plasma interaction of the SI pulse under the plasma conditions expected at the end of the compression pulse, as given in numerical simulations of SI of an ICF target [3, 4]: absorption by inverse bremsstrahlung, selffocusing (SF) and linear gain calculations of stimulated Brillouin (SBS) and Raman (SRS) backscattering in an exponential density profile. We then present in section 3 the results relevant to the SI scheme from experiments performed with the LULI 6-beam, the LULI 2000 and the LIL facilities. We conclude in section 4 with the implications of these experimental results for the SI scheme. 2. Laser–plasma coupling in the SI context Target designs for direct drive have been considered at different laser driver energies with optimized laser pulse durations and capsule size. The corresponding plasma parameters slightly differ but the density profile always has an approximately exponential shape with a characteristic scale length that depends on the laser energy. Direct drive designs for the NIF at the megajoule level, with 15 ns laser pulse duration, result in a 750 µm scalelength at the time of maximum intensity [7]. Designs for HiPER, which correspond to a compression energy ∼5 times lower than NIF, give a 300 µm scalelength at the end of the compression pulse when the SI pulse is fired. In all cases, the electron temperature varies in the 2–4 keV range along the laser beam path. We first present a calculation of the inverse bremsstrahlung absorption of the SI pulse in the hot plasma produced at the end of the compression phase and then analyse the propagation of the SI pulse in the underdense plasma in front of the absorption region. A rapid evaluation of the linear SRS and SBS gains is also provided. 2.1. Inverse bremsstrahlung absorption The absorption of the laser beam at a wavelength of 351 nm is considered along the exponential density gradient with 300 and 750 µm scalelengths for two constant electron temperatures (2 and 4 keV). Collisional absorption is calculated analytically using the well-known formula for electron–ion collision frequency [6]. The corresponding results are shown in figure 1. It is found that for beams incident along the density gradient (along the initial target normal), the energy is fully absorbed before reaching the critical density. These first estimates provide the location of the main absorption region: significant absorption occurs above 0.5nc for Ln = 300 µm and above 0.35nc for Ln = 750 µm. 2.2. Critical power for SF SF/filamentation is an important process that can modify the propagation of the high-intensity beam in the lower density part of the plasma located in front of the absorption region. SF of the smoothed laser beam takes place in the laser speckles and is characterized by the 3

Plasma Phys. Control. Fusion 53 (2011) 124034

S Depierreux et al

ne/nc 300 µm I/I0 2keV I/I0 4 keV

(a) 1.2

1.2

1

1

0.8

0.8

0.6

0.6

0.4

0.4

0.2

0.2

0

0

500

1000

ne/nc 750µm I/I0 2keV I/I0 4keV

(b)

0

1500

0

500

z (µm)

1000

1500

2000

z (µm)

Figure 1. Exponential density profiles along the propagation axis z and associated intensity profiles resulting from inverse bremsstrahlung absorption for Ln = 300 µm (a) and Ln = 750 µm (b).

/Pc8 (a) Ln=300 µm

/Pc5 (b) Ln=750 µm

7

4

6 5

3

4

2

3 2

1

1 0

0

500

1000

z (µm)

1500

2 keV 5×1014 W/cm² 2 keV 5×1015 W/cm² 2 keV 1016 W/cm²

0

0

500

1000

z (µm)

1500

2000

4 keV 5×1014W/cm² 4 keV 5×1015 W/cm² 4 keV 1016 W/cm²

Figure 2. Dependence of the SF parameter, #P $/Pc , for a speckle with a radius of 2.5 µm at different laser intensities in plasmas at 2 and 4 keV with Ln = 300 µm (a) and Ln = 750 µm (b).

quantity #P $/Pc , which is the ratio of the average speckle power to the critical power with Pc = 34Te (keV)(1 − ne /nc )1/2 /(ne /nc ) [6]. The average speckle power is calculated for a speckle radius of 2.5 µm, which corresponds to an f /10 0.351 µm laser beam typical of the high-energy laser facilities. We calculate #P $ = 982 MW for I = 5 × 1015 W cm−2 and #P $ = 1963 MW for I = 1016 W cm−2 . The quantity #P $/Pc is given in figure 2 along the propagation axis under the same conditions as in figure 1. We find that for laser intensities considered for the SI beam, #P $/Pc is well above unity in front of the absorption region. Moreover, the forward stimulated Brillouin scattering may develop even below the SF threshold leading to strong beam spraying [8] and resulting in a reduced local intensity of the laser beam inside the plasma. Consequently, a smaller shock pressure will be generated compared to what is expected. These results indicate that some strategy must be identified to control SF and beam spraying, such as using a shock laser pulse with a smaller f -number or a smaller wavelength. 4

Plasma Phys. Control. Fusion 53 (2011) 124034

S Depierreux et al

The use of laser smoothing through spectral dispersion and multiple polarizations has been shown to be efficient and should be taken into consideration. 2.3. Backscattering instabilities The SI pulse is fired at the very end of the compression pulse when stagnation has begun so that the expansion velocity is almost constant over the underdense plasma of interest. The exponential density profile can be linearized for SRS and SBS gain calculations with an associated scalelength equal to Ln . The electron temperature, as given by hydrodynamics simulations, is inhomogeneous. When taking everything into account, it is found that the SBS amplification is limited either by the density or temperature gradients and that the corresponding Rosenbluth gain reaches a very high value (ARos > 50) for laser intensities above 5 × 1015 W cm−2 . These linear gain calculations indicate that non-linear saturation mechanisms will finally determine the level of SBS reflectivities of the SI pulse. The situation may be different for SRS. Indeed, the SRS amplification is limited by the density gradient and the corresponding Rosenbluth gain is equal to 7 for 0.351 µm laser wavelength and an intensity of 1016 W cm−2 below the quarter critical density in the exponential density profile with Ln = 300 µm. For this scalelength, it is so expected that the SRS amplification will remain limited by the density inhomogeneity even under these conditions of high laser intensities. With the longer density scalelength (Ln > 750 µm), significant SRS gains (ARos > 20) may be reached below quarter critical for the highest intensity beam (1016 W cm−2 ). These results are in agreement with those discussed in [9]. The question of the SRS growth close to quarter critical is discussed in detail in this reference. 3. Experimental results Three series of experiments were conducted with three different installations, taking advantage of their particular characteristics to make complementary studies. Single hot spot experiments were designed with the LULI 6-beam facility to study the propagation of the laser pulse at a high intensity (I ∼ 1016 W cm−2 ). The LULI 2000 and the LIL installations were used to measure the stimulated Brillouin and Raman scattering in long and hot plasmas that can be produced with kilojoule systems. The results concerning filamentation, SF, SBS and SRS are presented in this section. 3.1. SF and filamentation Experiments were performed with the LULI 6-beam facility in a preformed plasma with an interaction beam fired at 1.053 µm with a pulse duration of 600 ps [10]. The laser interaction beam was diffraction limited, using a deformable mirror and apodizing the beam, with a corresponding spot size measured at half-maximum of 21.5 µm. The maximum intensity of this beam was I λ2 = 1016 W cm−2 µm2 . The plasma was preformed and heated using three beams at 0.53 µm. During the interaction pulse, the plasma density was below critical and decreased as a function of time. The density had an approximately inverse parabolic shape along the laser axis with a scale length of ∼700 µm. The electron temperature was 0.6 keV. A sophisticated diagnostic was set up to analyse the intensity distribution inside the plasma with high spatial and temporal resolutions. A transverse plane to the beam propagation axis across the plasma was imaged with high magnification in transmission [11] both on a streak camera and on two optical gated cameras. 5

Plasma Phys. Control. Fusion 53 (2011) 124034

S Depierreux et al

Figure 3. Time- and space-resolved images of the transmitted light for P /Pc = 0.23, 2.3 and 23.

(a)

RSBS (%)

3ω 3⊃ 2⊃ 2ω 1⊃ 1ω

RSBS (%) 25

(b)

20 10

15 10

1

5 0.1 0

15

2×10 4×10

15

6×10

15

8×10

I (W/cm²)

15

1×10

16

0

1

1.5

2

2.5

3

3.5

time (ns)

Figure 4. (a) Stimulated Brillouin backscattering reflectivities measured at 1ω in the LULI 6-beam experiment and at 2 and 3ω in the LULI 2000 experiment; (b) time-resolved SBS reflectivity measured at 2ω in the LULI 2000 experiment.

Experimental images of the intensity distribution in the plasma as a function of time are shown for three laser intensities corresponding to P /Pc = 0.23 (I = 1013 W cm−2 ), 2.3 (I = 1014 W cm−2 ) and 22.8 (I = 1015 W cm−2 ) in figure 3. At low laser intensity, the beam is little affected by its propagation through the plasma. At P /Pc ∼ 2, some modification of the trajectory, such as a dancing filament, of the beam can be observed. For P /Pc & 1, the initial beam structure is completely lost: the beam has split in many smaller speckles extending over a size at least four times larger than the initial beam diameter. The observed transition between linear and non-linear propagation of the laser pulse happens for a laser intensity as predicted by the critical power law. This effect is important for the SI scheme as the average intensity inside the plasma may be lowered by at least one order of magnitude. 3.2. Stimulated Brillouin scattering (SBS) SBS reflectivities were measured in preformed plasmas for three laser wavelengths (ω, 2ω and 3ω) in LULI experiments. At 1.053 µm, in the 6-beam experiment described in section 3.1, the 6

Plasma Phys. Control. Fusion 53 (2011) 124034

S Depierreux et al

SBS reflectivities, represented by the red squares in figure 4(a), were measured as a function of the laser intensity. The SBS threshold is observed at 5 × 1014 W cm−2 , it is followed by a rapid rise over one and a half order of magnitude and a saturation of the SBS reflectivity at the level of 10% for intensities higher than 1015 W cm−2 . Experiments were performed with the LULI 2000 facility to study the interaction of 526 and 351 nm beams (with maximum I λ2 = 1015 W cm−2 µm2 ) with a plasma preformed by irradiating a thick plastic foil. At the time of interaction, the density profile was exponential with the scalelength Ln = 100 µm. The electron temperature was ∼1 keV. The results are also plotted in figure 4(a). The growth and saturation of SBS under these conditions of very high SBS gains (ARos max > 100) have been reported in [13]. A first interesting result is that, in all the cases, we observe the saturation of the SBS reflectivity with increasing laser intensity and that the saturation level is the same for the three laser wavelengths. The growth and saturation of SBS under these conditions have been the subjects of numerical and theoretical investigations [12, 13]. For the three experiments, it appears that the SBS saturation arises mainly from the interplay with SF and filamentation. This results either from density depletion in self-focused speckles which saturates the SBS gain or from spraying of the laser beam in the front part of the plasma that reduces the growth of SBS deeper inside the plasma. These scenarios for the SBS saturation should also apply to the interaction of the SI beams which takes place in an exponential density profile [13] with a significant value of #P $/Pc . Our experimental results indicate that, despite high Rosenbluth gains, limited levels of SBS reflectivities ( 750 µm are used. Furthermore, another absorption mechanism operating at quarter critical density and related to SRS has recently been identified in this high-intensity regime that could provide an alternative interesting way of efficient absorption of the SI beams [9]. Stimulated Brillouin backscattering reflectivities are measured in different experiments performed in CH plasmas. Either experiments performed at 1 µm in exploding foils or experiments performed at 526 and 351 nm in an exponential density profile with very high SBS gains relevant to those expected for the SI beams have shown time-integrated backscattering SBS reflectivities limited at the 10% level. In these experiments, self-focusing and beam spray are the principal saturation mechanisms of SBS reflectivity. This scenario is certainly relevant also to the SI interaction conditions [13] where the quantity #P $/Pc is larger than unity in a large part of the underdense plasma. Under these conditions, significant level of near backscattering SBS is expected. When including this contribution, we find a maximum, time-resolved SBS reflectivity that reaches the 30% level. Another risk may result from the stimulated Raman backscattering instability. Our experiments have shown that significant SRS reflectivities in the few per cent range can be observed in exponential density profiles with scalelengths ∼300 µm at moderate laser intensities even if the Rosenbluth gain is as small as 2. In the SI highest intensity case, even higher SRS reflectivities could be expected. However, as already mentioned, SRS could also provide an efficient way of converting laser energy into hot electrons with energies of a few tens of kilovolts that may be beneficial for the SI scheme. Acknowledgments The experiment on LIL was coordinated under the auspice of the Institute Lasers and Plasmas. The authors acknowledge the support of the operation team of the LIL and LULI 2000 facilities who made these experiments possible. The authors acknowledge the support of the ANR contract (CORPARIN). Participation of JL was supported by the Czech Ministry of Education (projects LC528 and MSM6840770022). References [1] Lindl J D 1998 Inertial Confinement Fusion (New York: Springer) Atzeni S and Meyer-ter-Vehn J 2004 The Physics of Inertial Fusion (Oxford: Oxford University Press) [2] Shcherbakov V A 1983 Sov. J. Plasma Phys. 9 240 [3] Betti R et al 2007 Phys. Rev. Lett. 98 155001 Ribeyre X et al 2009 Plasma Phys. Control. Fusion 51 015013 Perkins L J et al 2009 Phys. Rev. Lett. 103 045004 [4] Schmitt A J et al 2010 Phys. Plasmas 17 042701 [5] Theobald W et al 2009 Plasma Phys. Control. Fusion 51 124052 [6] Pesme D et al Interaction Collisionnelle et Collective, La Fusion Thermonucl´eaire par Laser vol 1, ed R Dautray and J P Watteau (Paris: Eyrolles) 9

Plasma Phys. Control. Fusion 53 (2011) 124034 [7] Regan S et al 1999 Phys. Plasmas 6 2072 [8] Schmitt A J and Afeyan B B 1998 Phys. Plasmas 5 503 Lushnikov P M and Rose H A 2006 Plasma Phys. Control. Fusion 48 1501 Grech M et al 2009 Phys. Rev. Lett. 102 155001 [9] Klimo O et al 2010 Plasma Phys. Control. Fusion 52 055013 Klimo O et al 2011 Phys. Plasmas 18 082709 [10] Labaune C et al 2004 Plasma Phys. Control. Fusion 46 B301 Bandulet H C et al 2004 Phys. Rev. Lett. 93 035002 [11] Lewis K et al 2005 Rev. Sci. Inst. 76 093502 [12] Masson-Laborde P E et al 2000 J. Phys. IV 133 247 [13] Depierreux S et al 2009 Phys. Rev. Lett. 103 115001 Depierreux S et al 2011 Phys. Plasmas submitted [14] Depierreux S et al 2009 Phys. Rev. Lett. 102 195005

10

S Depierreux et al