Author Affiliations
1LULI–CNRS, CEA, Sorbonne Université, Ecole Polytechnique, Institut Polytechnique de Paris, F-91128 Palaiseau cedex, France2Sorbonne Université, Observatoire de Paris, Université PSL, CNRS, LERMA, F-75005 Paris, France3Office National d’Etudes et de Recherches Aérospatiales (ONERA), F-91123 Palaiseau, France4IAP, Russian Academy of Sciences, 603950 Nizhny Novgorod, Russia5LNCMI, UPR 3228, CNRS-UGA-UPS-INSA, F-31400 Toulouse, Franceshow less
Fig. 1. Schematic of the simulation setup. The solid target is in the xy plane, the direction of the externally applied magnetic field is along the x axis, and the plasma plume expands along the z axis.
Fig. 2. Plasma dynamics and slab formation. (a)–(d) show global 3D renderings of the decimal logarithm of the mass density at different times (5, 15, 30, and 50 ns) after the laser pulse.
Fig. 3. Decimal logarithm of the mass density ρ (A) and temperatures Te (B) and Ti (C) of the plasma plume sliced in the middle of the xz plane (a) and the yz plane (b) at different times 8, 18, 28, 38, and 48 ns (1–5). The temperatures in (B) and (C) share the same color map. The magnetic field directions are shown in (A1a) and (A1b).
Fig. 4. Plasma expansion and diamagnetic cavity formation: decimal logarithm of the electron number density integrated along or perpendicular to the magnetic field, i.e., the xz plane (a) or the yz plane (b), at 5 ns after the laser ablation. The color map corresponds to log10∫nedy in (a) and log10∫nedx in (b), in cm−2. The black arrows show the direction and magnitude of the velocities, and the light blue lines and arrows show the direction of the magnetic field in the xz plane (they do not appear as continuous lines, because they are taken out of the sliced plane).
Fig. 5. MRTI growth during slab formation. (a) Zoomed view of the plasma/vacuum interface at t = 8 ns. Since the density map here is a slice at the middle plane x = 0, the color map corresponds to log10 ne in cm−3. Cyan lines show the contours of the current density magnitude, while the dashed line contours show the ion temperature Ti. (b) Temperature dependence of the fastest growing mode and the growth time for MRTI.
Fig. 6. Slab formation. (a) Collapse of the cavity at t = 15 ns around 3 mm < z < 5 mm (highlighted by the gray box). (b) Fully grown flutes at t = 30 ns reach the boundary of y = 4 mm. The collapsed region propagates along the z direction, arriving at 5 mm < z < 8 mm, and the slab is formed. The color map corresponds to log10∫nedx in cm−2. The black arrows are the velocity vectors.
Fig. 7. Slab propagation with different initial magnetic field strengths in the simulation with asymmetric perturbation. (a)–(c) are for Bx0 = 10, 20, and 30 T, respectively. The color map corresponds to log10∫nedx in cm−2.
Fig. 8. Typical spatial pattern of the laser intensity deposition on target in related experiments. The pattern is due to the fact that the laser is actually defocused on the plane of the target, so that its intensity is not too large, and hence we have on target a laser pattern with diffraction rings, which is typical of an intermediate plane before focus. The color map represents the normalized intensity.
Fig. 9. Spatial distribution of the Bessel-like perturbation of the plasma velocity: (a) u1,1, for the m = 1, n = 1 mode; (b) u1,2, for the m = 1, n = 2 mode. Other parameters are r0 = 0.5 mm, α1,1 = 3.8317, and α1,2 = 7.0156. The color map corresponds to the perturbation strength.
Fig. 10. Slab propagation with different modes of Bessel-like perturbation for different magnetic field strengths at the end of the simulation at t = 50 ns. The color map corresponds to log10∫nedx in cm−2, i.e., to the decimal logarithm of the electron number density integrated along the external magnetic field. (a) and (b) are for Bx = 10 T, while (c) and (d) are for Bx = 30 T. In (a) and (c), the initial Bessel-like perturbation is in the m = 1, n = 1 mode, while in (b) and (d), it is in the m = 1, n = 2 mode.
Time | 5 ns | 30 ns | 50 ns |
---|
Region | Shock envelope | MRTI flutes | Slab |
---|
Local measurements | Characteristic length L0 (μm) | 100 | 100 | 100 | Averaged atomic number A | 17.3 | 17.3 | 17.3 | Effective charge state Zeff | 6.5 | 8.0 | 8.0 | Electron number density ne (cm−3) | 1.0 × 1019 | 5.0 × 1017 | 1.0 × 1018 | Ion number density ni (cm−3) | 1.5 × 1018 | 7.0 × 1016 | 1.5 × 1017 | Mass density ρ (g/cm3) | 4.0 × 10−5 | 2.0 × 10−6 | 4.0 × 10−6 | Electron temperature Te (eV) | 150.0 | 120.0 | 100.0 | Ion temperature Ti (eV) | 400.0 | 90.0 | 100.0 | Flow velocity V0 (km/s) | 350.0 | 150.0 | 400.0 | Magnetic field strength B (T) | 25.0 | 29.0 | 30.0 | Calculated dimensionless parameters | Electron collisionality λmfp,e/L0 | 0.06 | 0.6 | 0.2 | Ion collisionality λmfp,i/L0 | 0.1 | 0.1 | 0.1 | Electron magnetization (λmfp,e/rL,e) | 5.0 | 680.0 | 280.0 | Ion magnetization (λmfp,i/rL,i) | 0.03 | 0.5 | 0.3 | Mach number M | 3.1 | 1.6 | 4.7 | Alfvénic Mach number MA | 3.3 | 0.3 | 1.0 | Magnetosonic Mach number Mms | 2.3 | 0.3 | 1.0 | Plasma thermal beta βther | 1.4 | 3.0 × 10−2 | 5.0 × 10−2 | Plasma dynamic beta βdyn | 22.0 | 0.1 | 2.0 |
|
Table 1. Measured plasma conditions and calculated dimensionless parameters for the case with initial magnetic field Bx0 = 30 T in different regions as indicated in Fig. 2, i.e., the shock envelope at t = 5 ns, the MRTI flutes at t = 30 ns, and the slab at t = 50 ns. λmfp,s (with s = i and e for ions and electrons, respectively) is the mean free path39 and rL,s is the Larmor radius. The Mach number is the ratio of the flow velocity to the sound velocity, while the Alfvénic Mach number is the ratio of the flow velocity to the Alfvén velocity. The thermal and dynamic beta parameters are the ratios of the plasma thermal and ram pressures, respectively, to the magnetic pressure.