• High Power Laser Science and Engineering
  • Vol. 7, Issue 1, 010000e3 (2019)
S. Kawata1、†,*, T. Karino1, and Y. J. Gu2、3
Author Affiliations
  • 1Graduate School of Engineering, Utsunomiya University, Utsunomiya 321-8585, Japan
  • 2Institute of Physics of the ASCR, ELI-Beamlines, Na Slovance 2, 18221 Prague, Czech Republic
  • 3Institute of Plasma Physics of the CAS, Za Slovankou 1782/3, 18200 Prague, Czech Republic
  • show less
    DOI: 10.1017/hpl.2018.61 Cite this Article Set citation alerts
    S. Kawata, T. Karino, Y. J. Gu. Dynamic stabilization of plasma instability[J]. High Power Laser Science and Engineering, 2019, 7(1): 010000e3 Copy Citation Text show less
    An example concept of feedback control. (a) At $t=0$ a perturbation is imposed. The initial perturbation may grow at instability onset. (b) After $\unicode[STIX]{x0394}t$, if the feedback control works on the system, another perturbation, which has an inverse phase with the detected amplitude at $t=0$, is actively imposed, so that (c) the actual perturbation amplitude is mitigated very well after the superposition of the initial and additional perturbations.
    Fig. 1. An example concept of feedback control. (a) At $t=0$ a perturbation is imposed. The initial perturbation may grow at instability onset. (b) After $\unicode[STIX]{x0394}t$, if the feedback control works on the system, another perturbation, which has an inverse phase with the detected amplitude at $t=0$, is actively imposed, so that (c) the actual perturbation amplitude is mitigated very well after the superposition of the initial and additional perturbations.
    Kapitza’s pendulum, which can be stabilized by applying an additional strong and rapid acceleration $A\text{sin}\,\unicode[STIX]{x1D714}t$.
    Fig. 2. Kapitza’s pendulum, which can be stabilized by applying an additional strong and rapid acceleration $A\text{sin}\,\unicode[STIX]{x1D714}t$.
    Superposition of perturbations defined by the wobbling driver beam. At each time the wobbler provides a perturbation, whose amplitude and phase are defined by the wobbler itself. If the system is unstable, each perturbation is a source of instability. At a certain time the overall perturbation is the superposition of the growing perturbations. The superimposed perturbation growth is mitigated by the beam wobbling motion.
    Fig. 3. Superposition of perturbations defined by the wobbling driver beam. At each time the wobbler provides a perturbation, whose amplitude and phase are defined by the wobbler itself. If the system is unstable, each perturbation is a source of instability. At a certain time the overall perturbation is the superposition of the growing perturbations. The superimposed perturbation growth is mitigated by the beam wobbling motion.
    Example simulation results for the Rayleigh–Taylor instability (RTI) mitigation. $\unicode[STIX]{x1D6FF}g$ is 10% of the acceleration $g_{0}$ and oscillates with the frequency of $\unicode[STIX]{x1D6FA}=\unicode[STIX]{x1D6FE}$. As shown above and in Equation (3), the dynamic instability mitigation mechanism works well to mitigate the instability growth.
    Fig. 4. Example simulation results for the Rayleigh–Taylor instability (RTI) mitigation. $\unicode[STIX]{x1D6FF}g$ is 10% of the acceleration $g_{0}$ and oscillates with the frequency of $\unicode[STIX]{x1D6FA}=\unicode[STIX]{x1D6FE}$. As shown above and in Equation (3), the dynamic instability mitigation mechanism works well to mitigate the instability growth.
    Fluid simulation results for the RTI mitigation for the time-dependent $\unicode[STIX]{x1D6FF}g(t)=\unicode[STIX]{x1D6FF}g-\unicode[STIX]{x1D6E5}\sin \unicode[STIX]{x1D6FA}^{\prime }t$ at $t=5/\unicode[STIX]{x1D6FE}$. In the simulations $\unicode[STIX]{x1D6E5}=0.3$, and (a) $\unicode[STIX]{x1D6FA}^{\prime }=\unicode[STIX]{x1D6FA}/3$, (b) $\unicode[STIX]{x1D6FA}^{\prime }=\unicode[STIX]{x1D6FA}$ and (c) $\unicode[STIX]{x1D6FA}^{\prime }=3\unicode[STIX]{x1D6FA}$. The dynamic mitigation mechanism is robust against the time change of the perturbation amplitude $\unicode[STIX]{x1D6FF}g(t)$.
    Fig. 5. Fluid simulation results for the RTI mitigation for the time-dependent $\unicode[STIX]{x1D6FF}g(t)=\unicode[STIX]{x1D6FF}g-\unicode[STIX]{x1D6E5}\sin \unicode[STIX]{x1D6FA}^{\prime }t$ at $t=5/\unicode[STIX]{x1D6FE}$. In the simulations $\unicode[STIX]{x1D6E5}=0.3$, and (a) $\unicode[STIX]{x1D6FA}^{\prime }=\unicode[STIX]{x1D6FA}/3$, (b) $\unicode[STIX]{x1D6FA}^{\prime }=\unicode[STIX]{x1D6FA}$ and (c) $\unicode[STIX]{x1D6FA}^{\prime }=3\unicode[STIX]{x1D6FA}$. The dynamic mitigation mechanism is robust against the time change of the perturbation amplitude $\unicode[STIX]{x1D6FF}g(t)$.
    Fluid simulation results for the RTI mitigation for the time-dependent wobbling frequency $\unicode[STIX]{x1D6FA}(t)=\unicode[STIX]{x1D6FA}(1+\unicode[STIX]{x1D6E5}\sin \unicode[STIX]{x1D6FA}^{\prime }t)$ at $t=5/\unicode[STIX]{x1D6FE}$. In the simulations $\unicode[STIX]{x1D6E5}=0.3$, and (a) $\unicode[STIX]{x1D6FA}^{\prime }=\unicode[STIX]{x1D6FA}/3$, (b) $\unicode[STIX]{x1D6FA}^{\prime }=\unicode[STIX]{x1D6FA}$ and (c) $\unicode[STIX]{x1D6FA}^{\prime }=3\unicode[STIX]{x1D6FA}$. The dynamic mitigation mechanism is also robust against the time change of the perturbation frequency $\unicode[STIX]{x1D6FA}(t)$.
    Fig. 6. Fluid simulation results for the RTI mitigation for the time-dependent wobbling frequency $\unicode[STIX]{x1D6FA}(t)=\unicode[STIX]{x1D6FA}(1+\unicode[STIX]{x1D6E5}\sin \unicode[STIX]{x1D6FA}^{\prime }t)$ at $t=5/\unicode[STIX]{x1D6FE}$. In the simulations $\unicode[STIX]{x1D6E5}=0.3$, and (a) $\unicode[STIX]{x1D6FA}^{\prime }=\unicode[STIX]{x1D6FA}/3$, (b) $\unicode[STIX]{x1D6FA}^{\prime }=\unicode[STIX]{x1D6FA}$ and (c) $\unicode[STIX]{x1D6FA}^{\prime }=3\unicode[STIX]{x1D6FA}$. The dynamic mitigation mechanism is also robust against the time change of the perturbation frequency $\unicode[STIX]{x1D6FA}(t)$.
    Fluid simulation results for the RTI mitigation for the time-dependent wobbling wavelength $k(t)=k_{0}+\unicode[STIX]{x0394}ke^{i\unicode[STIX]{x1D6FA}_{k}^{\prime }t}$, at $t=5/\unicode[STIX]{x1D6FE}$. In the simulations $\unicode[STIX]{x0394}k/k_{0}=0.3$, and (a) $\unicode[STIX]{x1D6FA}_{k}^{\prime }=\unicode[STIX]{x1D6FA}/3$, (b) $\unicode[STIX]{x1D6FA}_{k}^{\prime }=\unicode[STIX]{x1D6FA}$ and (c) $\unicode[STIX]{x1D6FA}_{k}^{\prime }=3\unicode[STIX]{x1D6FA}$. The dynamic mitigation mechanism is also robust against the time change of the perturbation wavelength $k(t)$.
    Fig. 7. Fluid simulation results for the RTI mitigation for the time-dependent wobbling wavelength $k(t)=k_{0}+\unicode[STIX]{x0394}ke^{i\unicode[STIX]{x1D6FA}_{k}^{\prime }t}$, at $t=5/\unicode[STIX]{x1D6FE}$. In the simulations $\unicode[STIX]{x0394}k/k_{0}=0.3$, and (a) $\unicode[STIX]{x1D6FA}_{k}^{\prime }=\unicode[STIX]{x1D6FA}/3$, (b) $\unicode[STIX]{x1D6FA}_{k}^{\prime }=\unicode[STIX]{x1D6FA}$ and (c) $\unicode[STIX]{x1D6FA}_{k}^{\prime }=3\unicode[STIX]{x1D6FA}$. The dynamic mitigation mechanism is also robust against the time change of the perturbation wavelength $k(t)$.
    Filamentation instability. In this case an electron beam has a density perturbation in the transverse direction, and is injected into a plasma. In the plasma return current is induced to compensate for the electron beam current. The perturbed electron beam itself defines the filamentation instability phase, and the e-beam axis oscillates in the $y$ direction in this example case. Therefore, the filamentation instability is mitigated by the dynamic stabilization mechanism.
    Fig. 8. Filamentation instability. In this case an electron beam has a density perturbation in the transverse direction, and is injected into a plasma. In the plasma return current is induced to compensate for the electron beam current. The perturbed electron beam itself defines the filamentation instability phase, and the e-beam axis oscillates in the $y$ direction in this example case. Therefore, the filamentation instability is mitigated by the dynamic stabilization mechanism.
    Dynamic stabilization mechanism for the filamentation instability. (a) A modulated electron beam is imposed to induce the filamentation instability. The electron beam axis is wobbled or oscillates transversally with its frequency of $\unicode[STIX]{x1D6FA}$. (b) At a later time its phase-shifted perturbation is additionally imposed by the electron beam itself. The overall perturbation is the superimposition of all the perturbations, and the filamentation instability is dynamically stabilized.
    Fig. 9. Dynamic stabilization mechanism for the filamentation instability. (a) A modulated electron beam is imposed to induce the filamentation instability. The electron beam axis is wobbled or oscillates transversally with its frequency of $\unicode[STIX]{x1D6FA}$. (b) At a later time its phase-shifted perturbation is additionally imposed by the electron beam itself. The overall perturbation is the superimposition of all the perturbations, and the filamentation instability is dynamically stabilized.
    Filamentation instability simulation results without and with the electron beam oscillation. The current density $J_{x}$ is shown at each time step. When the electron beam axis oscillates in the $y$ direction ((d)–(f) and (g)–(i)), the filamentation instability growth is clearly mitigated.
    Fig. 10. Filamentation instability simulation results without and with the electron beam oscillation. The current density $J_{x}$ is shown at each time step. When the electron beam axis oscillates in the $y$ direction ((d)–(f) and (g)–(i)), the filamentation instability growth is clearly mitigated.
    Magnetic field $B_{z}$ for the filamentation instability without and with the electron beam oscillation.
    Fig. 11. Magnetic field $B_{z}$ for the filamentation instability without and with the electron beam oscillation.
    Histories of the normalized magnetic field energy $U_{Bz}\propto |B_{z}|^{2}$. When the electron beam transverse oscillation frequency $\unicode[STIX]{x1D6FA}$ in $y$ becomes larger than or comparable to $\unicode[STIX]{x1D6FE}_{F}$, the dynamic stabilization effect is remarkable.
    Fig. 12. Histories of the normalized magnetic field energy $U_{Bz}\propto |B_{z}|^{2}$. When the electron beam transverse oscillation frequency $\unicode[STIX]{x1D6FA}$ in $y$ becomes larger than or comparable to $\unicode[STIX]{x1D6FE}_{F}$, the dynamic stabilization effect is remarkable.
    3D PIC simulation results for the filamentation instability growth at (a) $t=0$ and (b) $t=32\unicode[STIX]{x1D714}_{pe}$ without the electron beam wobbling motion, and (c) $t=0$ and (d) $t=32\unicode[STIX]{x1D714}_{pe}$ with the wobbling motion. The histories of the normalized magnetic field energy $U_{Bz}\propto |B_{z}|^{2}$ is shown in (e). When the electron beam transverse oscillation frequency $\unicode[STIX]{x1D6FA}$ in $y$ becomes larger than or comparable to $\unicode[STIX]{x1D6FE}_{F}$, the dynamic stabilization effect is also remarkable. The 3D simulation results also confirm the instability mitigation mechanism in the plasma.
    Fig. 13. 3D PIC simulation results for the filamentation instability growth at (a) $t=0$ and (b) $t=32\unicode[STIX]{x1D714}_{pe}$ without the electron beam wobbling motion, and (c) $t=0$ and (d) $t=32\unicode[STIX]{x1D714}_{pe}$ with the wobbling motion. The histories of the normalized magnetic field energy $U_{Bz}\propto |B_{z}|^{2}$ is shown in (e). When the electron beam transverse oscillation frequency $\unicode[STIX]{x1D6FA}$ in $y$ becomes larger than or comparable to $\unicode[STIX]{x1D6FE}_{F}$, the dynamic stabilization effect is also remarkable. The 3D simulation results also confirm the instability mitigation mechanism in the plasma.
    S. Kawata, T. Karino, Y. J. Gu. Dynamic stabilization of plasma instability[J]. High Power Laser Science and Engineering, 2019, 7(1): 010000e3
    Download Citation