
Citation: | Jiachen TONG, Haiying LI, Bin XU, Songyang WU, Lu BAI. Excitation and power spectrum analysis of electromagnetic radiation for the plasma wake of reentry vehicles[J]. Plasma Science and Technology, 2023, 25(5): 055301. DOI: 10.1088/2058-6272/aca7ad |
The plasma wake of reentry vehicles has the advantages of extensive space range and long traceability, which provides new possibilities for the detection and monitoring of reentry vehicles. Based on the Zakharov model, this work investigates the excitation and power spectrum characteristics of electromagnetic radiation for the plasma wake of a typical reentry vehicle. With the aid of parametric decay instability, the excitation condition of electromagnetic radiation for a typical plasma wake is evaluated first. The power spectrum characteristics of electromagnetic radiation, including the effects of both the flight parameters and incident wave parameters are analyzed in detail. The results show that when the phenomenon of excited electromagnetic radiation occurs, plasma wakes closer to the bottom of the vehicle and with faster speeds require higher incident frequencies and thresholds of the electric field. As the frequency of the incident wave increases, peaks appear in the power spectra of plasma wakes, and their magnitudes increase gradually. The frequency shifts of the secondary peaks are equal, whereas, the peaks of the downshifted spectral lines are generally larger than those of the upshifted spectral lines. The work in this paper provides a new idea and method for the tracking of reentry vehicles, which has potential application value in the field of reentry vehicle detection.
Atmospheric pressure plasma jets (APPJs) have attracted lots of interest in the past decades due to their urgent applications [1–18]. To better understand the plasma characteristic, novel plasma diagnostic methods are needed. The main plasma parameters include the electron density, electric field, electron temperature, reactive oxygen species density, and so on [19–21]. Among these parameters, electric field is one of the most important parameters, which is related to the dynamic process of electron and high energy electron distribution [22].
On the other hand, it is extremely difficult to measure the electric field of APPJs. Recently, several methods for atmospheric pressure plasma electric field diagnosis have been proposed, such as the electro-optic probe method [18] and optical emission spectroscopy method [23–25]. However, the electro-optic probe usually needs to invade into plasma, which interferes with the gas flow. Surface charge accumulation on the probe also affects the plasma characteristic. Regarding the optical emission spectroscopy method, it cannot be applied when the optical emission of the plasma is weak. For example, it cannot be used to diagnose the electric field of the dark channel of a plasma jet. In addition, it requires that the excited states are generated by electron collision excitation rather than excited by metastable state species.
Recently, an electric field diagnosis method based on the electric field induced second harmonic (E-FISH) method was proposed [26–30], which has been used to measure the nonlinear polarization coefficients of different gas components, static electric field of DC corona, transient electric field of dielectric barrier discharge and blade electrode discharge [27–30]. The main advantages of E-FISH are its high spatial and temporal resolution. Additionally, it can be used to measure the electric field for most situations; the working gas of the plasma can be molecular gas or atomic gas. In addition, the electric field components in different directions can be obtained by using a laser with different polarizations. Furthermore, the measurement does not depend on the optical emission of the plasma; thus, it can be used to measure the electric field of nonluminous regions such as the dark channel of APPJs.
In this work, an E-FISH method based on a nanosecond pulse laser is used to obtain the time-resolved electric field of APPJs driven by positive or negative pulse voltage. Additionally, the effects of the percentages of nitrogen and oxygen mixture of a helium plasma jet on the temporal behavior of the electric field are investigated.
Atmospheric pressure plasma jets (APPJs) have attracted lots of interest in the past decades due to their urgent applications [1–18]. To better understand the plasma characteristic, novel plasma diagnostic methods are needed. The main plasma parameters include the electron density, electric field, electron temperature, reactive oxygen species density, and so on [19–21]. Among these parameters, electric field is one of the most important parameters, which is related to the dynamic process of electron and high energy electron distribution [22].
On the other hand, it is extremely difficult to measure the electric field of APPJs. Recently, several methods for atmospheric pressure plasma electric field diagnosis have been proposed, such as the electro-optic probe method [18] and optical emission spectroscopy method [23–25]. However, the electro-optic probe usually needs to invade into plasma, which interferes with the gas flow. Surface charge accumulation on the probe also affects the plasma characteristic. Regarding the optical emission spectroscopy method, it cannot be applied when the optical emission of the plasma is weak. For example, it cannot be used to diagnose the electric field of the dark channel of a plasma jet. In addition, it requires that the excited states are generated by electron collision excitation rather than excited by metastable state species.
Recently, an electric field diagnosis method based on the electric field induced second harmonic (E-FISH) method was proposed [26–30], which has been used to measure the nonlinear polarization coefficients of different gas components, static electric field of DC corona, transient electric field of dielectric barrier discharge and blade electrode discharge [27–30]. The main advantages of E-FISH are its high spatial and temporal resolution. Additionally, it can be used to measure the electric field for most situations; the working gas of the plasma can be molecular gas or atomic gas. In addition, the electric field components in different directions can be obtained by using a laser with different polarizations. Furthermore, the measurement does not depend on the optical emission of the plasma; thus, it can be used to measure the electric field of nonluminous regions such as the dark channel of APPJs.
In this work, an E-FISH method based on a nanosecond pulse laser is used to obtain the time-resolved electric field of APPJs driven by positive or negative pulse voltage. Additionally, the effects of the percentages of nitrogen and oxygen mixture of a helium plasma jet on the temporal behavior of the electric field are investigated.
The E-FISH method is used to obtain the time-resolved electric field of a helium plasma jet. When a laser passes through a nonlinear medium, an external photoelectric field interacts with molecules, atoms and electrons in the plasma to produce a polarized photoelectric field. This field is expressed by the electric polarization intensity P in the strong field environment. When the laser is in a strong electric field, there is a nonlinear optical effect and the field induced second harmonic corresponds to the third-order nonlinearity [31, 32]. E-FISH generation can be described as a third-order nonlinear process involving both a light field and an external electric field [26]:
P(ω2)i=32Nχ(3)i,j,k,l(−ω2,0,ω1,ω1)EFjEω1kEω1l | (1) |
where
In the process of light propagation, the relationship between light field intensity and polarization intensity can be described in the form of a Maxwell wave equation, which for the second harmonic intensity
∇2Eω2(x)+ω22c2Eω2(x)=−μ0ω22P(ω2)i(x) | (2) |
where c is the speed of light in a vacuum and
By integrating the region along the x-axis and replacing
Eω2=3Nμ0ω2c4n2χ(3)i,j,k,l·EFj·(Eω1)2·∫eiΔkxdx | (3) |
where
According to the relationship between the intensity of the laser and the intensity of the light field, replacing
Iω2=(3μ0ω2c4n2)2N2·(EFj)2(Iω1)2L2·[sinΔk·L2Δk·L2]2 | (4) |
where L = 2ZR = π/Δk is the coherence length and ZR is the Rayleigh length. There is a quadratic relationship between the measured second harmonic signal and the applied electric field:
Iω2∝N2·(EFj)2·(Iω1)2 | (5) |
Hence, the measured electric field can be rewritten as:
EFj=A√Iω2 | (6) |
where A is the calibration constant, which is calibrated with a known electric field, and all unknown parameters are included in the calibration constant.
During the experiment, the Rayleigh length ZR determines the spatial resolution along the laser direction and the resolution in the vertical laser direction is determined by the waist diameter of the laser beam. The beam waist diameter is affected by the laser itself and the focusing ability of the focusing lens. The strong focusing ability makes the beam waist diameter smaller and the spatial resolution higher. However, the strong focusing of the laser increases the energy density, thus increasing the risk of photo-breakdown. Finally, a focusing lens with focal length of 10 cm is selected, and the corresponding Rayleigh length is 2 mm.
In this experiment, a Nd: YAG nanosecond pulse laser is used as the pump laser. It has a frequency of 10 Hz, pulse width of 8 ns, wavelength of 1064 nm, and energy per shot of 60 mJ during the measurements. A long pass filter (LP, through above 950 nm) is used to filter out the mixed 532 nm light and other stray light. A half wave plate (HWP) is used to change the laser polarization direction from horizontal polarization to vertical polarization, so the electric field along the vertical direction is measured. Then the laser focuses through a focus lens (FL) with a focal length of 10 cm. The focus should be located in the center of position to be measured, so as to ensure the Rayleigh range covers the region of the electric field. After the laser passes through the center area, a collimator lens (CL) is used to change the laser back to parallel light, and a 532 nm reflector mirror (HF) is used to separate the E-FISH signal from the source laser. The laser emission is collected by a beam dump. The E-FISH signal is further filtered by a spectrometer and finally collected by an ICCD camera.
Because the E-FISH signal is weak, especially measured in helium, it is essential to reduce stray noise and increase the signal-to-noise ratio of E-FISH. Firstly, strict time synchronization should be carried out. The flash lamp and Q-switch of the laser are synchronized with the nanosecond pulse power supply and the ICCD camera through a digital delay generator (DG645). The exposure time of the ICCD camera is set to be the same as the laser pulse width to minimize stray light. Secondly, the refractive index of the collimator lens is slightly different between 1064 nm and 532 nm, and the position of the lens needs to be adjusted to ensure E-FISH light is accurately collimated. Then, the E-FISH signal is separated by a high-reflectivity mirror and focused into spectrometer by a focus lens. Finally, the signal is accumulated 600 times to enhance the intensity of E-FISH and reduce the measurement error. Figure 1 shows a schematic of the experimental setup.
Helium or helium mixed with a small amount of oxygen or nitrogen is used as the working gas. The total gas flow rate is 2 l min−1. A schematic of the plasma jet device is shown in figure 2. The high-voltage electrode is made of a tungsten needle, which is placed in a glass tube. The tip of the tungsten needle is 5 mm away from the nozzle of the glass tube. The inner diameter of the glass tube is 5 mm, the inner diameter of the nozzle is 1 mm, and there is no extra electrical ground during the measurements. The position of the nozzle is labelled as z = 0 mm. The high-voltage electrode is connected to a nanosecond pulsed high-voltage power supply with peak voltage up to 10 kV, frequency up to 10 kHz, and pulse width adjustable from 200 ns to DC. The voltage is measured by a voltage probe (Tektronix P6015) and recorded by an oscilloscope (Tektronix DPO7104). The dynamics of the plasma are captured by an ICCD camera (Andor DH-334). The plasma jet device is placed on a movable station adjusted by a micrometer to change the laser focusing position. The time-resolved electric field of the plasma can be obtained by changing the delay time of the power supply, the laser and the ICCD camera through the digital delay generator (DG645).
The E-FISH signal intensity is related to several factors, including the electric field, fundamental frequency laser intensity, interaction length between laser and the electric field, particle number density, etc. Therefore, the system needs to be calibrated to establish the relationship between the E-FISH signal and the electric field. During the calibration, the factors affecting the E-FISH signal intensity should be kept the same as the measurement conditions.
The E-FISH method is used to obtain the time-resolved electric field of a helium plasma jet. When a laser passes through a nonlinear medium, an external photoelectric field interacts with molecules, atoms and electrons in the plasma to produce a polarized photoelectric field. This field is expressed by the electric polarization intensity P in the strong field environment. When the laser is in a strong electric field, there is a nonlinear optical effect and the field induced second harmonic corresponds to the third-order nonlinearity [31, 32]. E-FISH generation can be described as a third-order nonlinear process involving both a light field and an external electric field [26]:
P(ω2)i=32Nχ(3)i,j,k,l(−ω2,0,ω1,ω1)EFjEω1kEω1l | (1) |
where
In the process of light propagation, the relationship between light field intensity and polarization intensity can be described in the form of a Maxwell wave equation, which for the second harmonic intensity
∇2Eω2(x)+ω22c2Eω2(x)=−μ0ω22P(ω2)i(x) | (2) |
where c is the speed of light in a vacuum and
By integrating the region along the x-axis and replacing
Eω2=3Nμ0ω2c4n2χ(3)i,j,k,l·EFj·(Eω1)2·∫eiΔkxdx | (3) |
where
According to the relationship between the intensity of the laser and the intensity of the light field, replacing
Iω2=(3μ0ω2c4n2)2N2·(EFj)2(Iω1)2L2·[sinΔk·L2Δk·L2]2 | (4) |
where L = 2ZR = π/Δk is the coherence length and ZR is the Rayleigh length. There is a quadratic relationship between the measured second harmonic signal and the applied electric field:
Iω2∝N2·(EFj)2·(Iω1)2 | (5) |
Hence, the measured electric field can be rewritten as:
EFj=A√Iω2 | (6) |
where A is the calibration constant, which is calibrated with a known electric field, and all unknown parameters are included in the calibration constant.
During the experiment, the Rayleigh length ZR determines the spatial resolution along the laser direction and the resolution in the vertical laser direction is determined by the waist diameter of the laser beam. The beam waist diameter is affected by the laser itself and the focusing ability of the focusing lens. The strong focusing ability makes the beam waist diameter smaller and the spatial resolution higher. However, the strong focusing of the laser increases the energy density, thus increasing the risk of photo-breakdown. Finally, a focusing lens with focal length of 10 cm is selected, and the corresponding Rayleigh length is 2 mm.
In this experiment, a Nd: YAG nanosecond pulse laser is used as the pump laser. It has a frequency of 10 Hz, pulse width of 8 ns, wavelength of 1064 nm, and energy per shot of 60 mJ during the measurements. A long pass filter (LP, through above 950 nm) is used to filter out the mixed 532 nm light and other stray light. A half wave plate (HWP) is used to change the laser polarization direction from horizontal polarization to vertical polarization, so the electric field along the vertical direction is measured. Then the laser focuses through a focus lens (FL) with a focal length of 10 cm. The focus should be located in the center of position to be measured, so as to ensure the Rayleigh range covers the region of the electric field. After the laser passes through the center area, a collimator lens (CL) is used to change the laser back to parallel light, and a 532 nm reflector mirror (HF) is used to separate the E-FISH signal from the source laser. The laser emission is collected by a beam dump. The E-FISH signal is further filtered by a spectrometer and finally collected by an ICCD camera.
Because the E-FISH signal is weak, especially measured in helium, it is essential to reduce stray noise and increase the signal-to-noise ratio of E-FISH. Firstly, strict time synchronization should be carried out. The flash lamp and Q-switch of the laser are synchronized with the nanosecond pulse power supply and the ICCD camera through a digital delay generator (DG645). The exposure time of the ICCD camera is set to be the same as the laser pulse width to minimize stray light. Secondly, the refractive index of the collimator lens is slightly different between 1064 nm and 532 nm, and the position of the lens needs to be adjusted to ensure E-FISH light is accurately collimated. Then, the E-FISH signal is separated by a high-reflectivity mirror and focused into spectrometer by a focus lens. Finally, the signal is accumulated 600 times to enhance the intensity of E-FISH and reduce the measurement error. Figure 1 shows a schematic of the experimental setup.
Helium or helium mixed with a small amount of oxygen or nitrogen is used as the working gas. The total gas flow rate is 2 l min−1. A schematic of the plasma jet device is shown in figure 2. The high-voltage electrode is made of a tungsten needle, which is placed in a glass tube. The tip of the tungsten needle is 5 mm away from the nozzle of the glass tube. The inner diameter of the glass tube is 5 mm, the inner diameter of the nozzle is 1 mm, and there is no extra electrical ground during the measurements. The position of the nozzle is labelled as z = 0 mm. The high-voltage electrode is connected to a nanosecond pulsed high-voltage power supply with peak voltage up to 10 kV, frequency up to 10 kHz, and pulse width adjustable from 200 ns to DC. The voltage is measured by a voltage probe (Tektronix P6015) and recorded by an oscilloscope (Tektronix DPO7104). The dynamics of the plasma are captured by an ICCD camera (Andor DH-334). The plasma jet device is placed on a movable station adjusted by a micrometer to change the laser focusing position. The time-resolved electric field of the plasma can be obtained by changing the delay time of the power supply, the laser and the ICCD camera through the digital delay generator (DG645).
The E-FISH signal intensity is related to several factors, including the electric field, fundamental frequency laser intensity, interaction length between laser and the electric field, particle number density, etc. Therefore, the system needs to be calibrated to establish the relationship between the E-FISH signal and the electric field. During the calibration, the factors affecting the E-FISH signal intensity should be kept the same as the measurement conditions.
To calibrate the E-FISH system, two parallel plate electrodes with a diameter of 3 cm are used. In order to avoid the influence of different gases on the third-order nonlinear polarizability
EFj=A√Iω2,A=5.624 |
However, the different working gases or electrode structures have an effect on generation of E-FISH; further calibration is needed [33]. In figure 4(a), when the applied voltage is 8 kV and no plasma is generated, the electric field at 5 mm away from the nozzle (z = 5 mm) of the plasma jet device is simulated by COMSOL. Then, we obtain the electric field gradient by setting different voltages and the corresponding E-FISH signal is measured. The calibration result is shown in figure 4(b); it can be found that mixing with oxygen has some influence on the intensity of E-FISH, but it still has a good linearity between the electric field and square of E-FISH, and the slope of the line indicates the calibration coefficient A in different situations. Note that when pure He is used for calibration, the calibration coefficient of jet A′ = 5.903 is very close to that of parallel plate A = 5.624, which indicates that only negligible air is mixed into the interaction length.
In the next part, the temporal electric field at 5 mm away from the nozzle (z = 5 mm) of the plasma jet device is measured when no plasma is generated. It can be seen from figure 5 that the waveform of the electric field has a similar shape to the applied voltage, which indicates that the measurement of the electric field has a good time resolution. Meanwhile, according to figure 4(a), when the applied voltage is 8 kV, the electric field is about 6 kV cm−1 at 5 mm away from the nozzle (z = 5 mm), and the calibration results are consistent with the simulation results. Note that the E-FISH signal after calibration obtains a line-integrated electric field along the laser path, and we use EL to indicate the line-integrated electric field.
Then, helium (99.999%) flows through the plasma jet device with a gas flow rate of 2 l min−1. It can be seen from figure 6 that the plasma propagates out of the nozzle at 60 ns, reaches the laser position at 90 ns, and the plasma bullet keeps propagating, leaving a dark channel behind after 130 ns. After 1000 ns, the falling edge of the applied voltage appears and a secondary discharge is ignited. At about 1085 ns, the plasma reaches the laser position. Then, the luminous region becomes weaker and eventually disappears.
Furthermore, the electric field at the position of 1 cm away from the nozzle (z = 1 cm) is measured. The exposure time of the ICCD is set at 10 ns, which is slightly larger than the laser pulse width of 8 ns. Figure 7 shows the electric field obtained through the E-FISH method. It is shown that the electric field reaches its peak of about 25 kV cm−1 at t1 = 85 ns. Then the plasma keeps propagating and the electric field at z = 1 cm decreases gradually to 5 kV cm−1 at about 450 ns and stays at about 5 kV cm−1 until the falling edge of the voltage pulse at 1000 ns. Then the electric field drops quickly and increases reversely to a peak of −14 kV cm−1 at t2 = 1085 ns. The results of electric field measurement are consistent with the dynamics of the plasma bullet. For a given position, the electric field is the highest at the moment when the head of the plasma bullet reaches the position. With the propagation of the plasma bullet, a dark channel appears; the electric field in the dark channel is about 5 kV cm−1. Note that the E-FISH method only obtains the intensity of the electric field; the polarity of the electric field needs additional judgment. According to the measured results and the dynamics of the plasma, there is no sudden change in electric field. Therefore, we define the positive direction of the z-axis as the positive polarity of the electric field; when the electric field drops to zero, the polarity of electric field changes to negative polarity.
As we know, the polarity of the applied voltage affects the dynamics of the APPJs, which is related to the difference in the electric field of the plasma channel [34–37]. In the next part, a negative polarity of the nanosecond pulsed power supply is used. The applied voltage is adjusted to −8 kV. All of the other parameters are the same as above, i.e. the pulse repetition rate is fixed at 5 kHz, and the pulse width is 1 μs. The dynamics of the plasma bullet behavior are captured by the ICCD camera, as shown in figure 8. When the applied voltage changes from 0 to −8 kV, the plasma propagates with a long tail. The maximum propagation distance of the plasma during the increase in the voltage is shorter than the case of 8 kV. On the other hand, the propagation of the plasma is more like a plasma bullet at the rising edge of the applied voltage.
For the applied voltage of −8 kV, the time-resolved electric field at 1 cm away from the nozzle (z = 1 cm) obtained by E-FISH method is shown in figure 9. The electric field increases quickly when the applied voltage changed from 0 to −8 kV. The electric field reaches its peak of −16 kV cm−1 at t1 = 85 ns. Then it decreases slowly to about −6 kV cm−1 and stays as a constant until 1000 ns. Afterwards, with the drop of the applied voltage from −8 kV to zero, the electric field quickly increases to a peak of about 20 kV cm−1. Subsequently it decreases to zero after about 500 ns. It is worthy to point out that the peak electric field at t2 during the secondary discharge is significantly higher than that at t1 during the first discharge.
According to figures 7 and 9, it is found that the peak electric field is about 20–25 kV cm−1 at the edge when the applied voltage changes from 0 to 8 kV or −8 to 0 kV. In addition, the plasma appears like a bullet during these periods. On the other hand, at the edge when the applied voltage changes from 8 to 0 kV or 0 to −8 kV, the peak electric field is about −14 to −16 kV cm−1. The plasma propagates much shorter distance and it has a relatively bright long tail.
In the next part, the effects of the pulse width and the pulse rising time on the electric field of the plasma bullet are measured. Firstly, the pulse width is reduced to 500 ns. All of the other parameters remain unchanged. The applied voltage waveform along with the measured electric field is shown in figure 10, which is similar to the results shown in figure 7. It reaches a peak value of about 21 kV cm−1 at 85 ns, and there is also a secondary discharge at the falling edge with a peak value of about −15 kV cm−1.
Secondly, the pulse rising time is increased from 60 to 200 ns. The applied voltage waveform and the obtained electric field at 1 cm away from the nozzle (z = 1 cm) by E-FISH method are shown in figure 11. It can be seen that the electric field increases slowly and reaches a peak of about 10.5 kV cm−1 at about 350 ns, which is much weaker than the case of 60 ns pulse rising time. Accordingly, the electric field induced by the secondary discharge at the falling edge is also much weaker, which is about −8 kV cm−1. Compared to the rapid falling edge of the voltage, the longer falling time makes discharge occur in the process of the voltage drop. This phenomenon also happens at the rising edge of the voltage. However, in figure 11, due to the imperfect rising edge of the voltage, the voltage still slowly rises after 200 ns, which causes the time of the electric field reaching its peak to be delayed to about 350 ns.
When APPJs are used for various applications, a trace amount of nitrogen or oxygen is often added to the working gas, which can increase the concentration of reactive nitrogen species (RNS) or reactive oxygen species (ROS) [38–41]. However, how the trace amount of nitrogen or oxygen affects the electric field and thus the propagation of the plasma jets has not been investigated. Note that the hyperpolarizability between helium and oxygen (or nitrogen) is significantly different; only a tiny amount of oxygen or nitrogen can be mixed into helium in this experiment. Moreover, in order to reduce the mixing of ambient air, sufficient helium flow rate is required and the measurement position should be close to the nozzle. Corresponding calibration is carried out under the same experimental conditions, which can be found in figure 4. Although the hyperpolarizability of argon is much closer to that of air, the helium jets are more widely used in applications, and the propagation of the helium plasma is like a bullet but that of argon plasma is like filaments; a more regular shape is convenient for simulation and correction [42].
Firstly, 0%, 0.5% and 1.0% oxygen are added into helium with a total gas flow rate of 2 l min−1. All the other parameters are kept the same. The electric field at 5 mm away from the nozzle (z = 5 mm) is measured by the E-FISH method as shown in figure 12. It is clearly shown that, when 0.5% of oxygen is added, the peak electric field at the rising edge reduced from 22 to 15 kV cm−1. Moreover, the electric field decreases much faster and it stays at a value of about 2 kV cm−1 lower than in the case of pure helium before the falling edge of the applied voltage. When 1% of oxygen is added, the peak electric field becomes even weaker.
Secondly, 0%, 0.5%, 1.0% and 1.5% nitrogen are added to the working gas and the corresponding electric field is measured by the E-FISH method as shown in figure 13. It can be seen that, when 0.5% nitrogen is added, the peak electric field actually slightly increases, especially at the falling edge of the voltage pulse, and it increases reversely from −12 to −16 kV cm−1. However, when 1.5% nitrogen is added, the peak electric field decreases from 24 to 19 kV cm−1 at the rising edge of the voltage pulse. Moreover, it is found that the electric field decreases much faster when a trace amount of nitrogen is added compared with the pure helium case.
We know that the higher the electric field is, the faster the plasma bullet propagates. To see if this is true, high-speed photographs of the plasma for He only, He + 0.5%N2, and He + 1.5%O2 are captured by an ICCD camera with exposure time of 5 ns. The obtained results are shown in figure 14. For the case of He only, the plasma already reaches 1.5 cm (red dotted line) out of the nozzle at 90 ns. For the case of He + 0.5%N2, the plasma propagates faster and reaches 1.5 cm at about 70 ns. When the working gas is changed to He + 1.5%O2, it does not reach 1.5 cm at 90 ns and it propagates even more slowly later.
To calibrate the E-FISH system, two parallel plate electrodes with a diameter of 3 cm are used. In order to avoid the influence of different gases on the third-order nonlinear polarizability
EFj=A√Iω2,A=5.624 |
However, the different working gases or electrode structures have an effect on generation of E-FISH; further calibration is needed [33]. In figure 4(a), when the applied voltage is 8 kV and no plasma is generated, the electric field at 5 mm away from the nozzle (z = 5 mm) of the plasma jet device is simulated by COMSOL. Then, we obtain the electric field gradient by setting different voltages and the corresponding E-FISH signal is measured. The calibration result is shown in figure 4(b); it can be found that mixing with oxygen has some influence on the intensity of E-FISH, but it still has a good linearity between the electric field and square of E-FISH, and the slope of the line indicates the calibration coefficient A in different situations. Note that when pure He is used for calibration, the calibration coefficient of jet A′ = 5.903 is very close to that of parallel plate A = 5.624, which indicates that only negligible air is mixed into the interaction length.
In the next part, the temporal electric field at 5 mm away from the nozzle (z = 5 mm) of the plasma jet device is measured when no plasma is generated. It can be seen from figure 5 that the waveform of the electric field has a similar shape to the applied voltage, which indicates that the measurement of the electric field has a good time resolution. Meanwhile, according to figure 4(a), when the applied voltage is 8 kV, the electric field is about 6 kV cm−1 at 5 mm away from the nozzle (z = 5 mm), and the calibration results are consistent with the simulation results. Note that the E-FISH signal after calibration obtains a line-integrated electric field along the laser path, and we use EL to indicate the line-integrated electric field.
Then, helium (99.999%) flows through the plasma jet device with a gas flow rate of 2 l min−1. It can be seen from figure 6 that the plasma propagates out of the nozzle at 60 ns, reaches the laser position at 90 ns, and the plasma bullet keeps propagating, leaving a dark channel behind after 130 ns. After 1000 ns, the falling edge of the applied voltage appears and a secondary discharge is ignited. At about 1085 ns, the plasma reaches the laser position. Then, the luminous region becomes weaker and eventually disappears.
Furthermore, the electric field at the position of 1 cm away from the nozzle (z = 1 cm) is measured. The exposure time of the ICCD is set at 10 ns, which is slightly larger than the laser pulse width of 8 ns. Figure 7 shows the electric field obtained through the E-FISH method. It is shown that the electric field reaches its peak of about 25 kV cm−1 at t1 = 85 ns. Then the plasma keeps propagating and the electric field at z = 1 cm decreases gradually to 5 kV cm−1 at about 450 ns and stays at about 5 kV cm−1 until the falling edge of the voltage pulse at 1000 ns. Then the electric field drops quickly and increases reversely to a peak of −14 kV cm−1 at t2 = 1085 ns. The results of electric field measurement are consistent with the dynamics of the plasma bullet. For a given position, the electric field is the highest at the moment when the head of the plasma bullet reaches the position. With the propagation of the plasma bullet, a dark channel appears; the electric field in the dark channel is about 5 kV cm−1. Note that the E-FISH method only obtains the intensity of the electric field; the polarity of the electric field needs additional judgment. According to the measured results and the dynamics of the plasma, there is no sudden change in electric field. Therefore, we define the positive direction of the z-axis as the positive polarity of the electric field; when the electric field drops to zero, the polarity of electric field changes to negative polarity.
As we know, the polarity of the applied voltage affects the dynamics of the APPJs, which is related to the difference in the electric field of the plasma channel [34–37]. In the next part, a negative polarity of the nanosecond pulsed power supply is used. The applied voltage is adjusted to −8 kV. All of the other parameters are the same as above, i.e. the pulse repetition rate is fixed at 5 kHz, and the pulse width is 1 μs. The dynamics of the plasma bullet behavior are captured by the ICCD camera, as shown in figure 8. When the applied voltage changes from 0 to −8 kV, the plasma propagates with a long tail. The maximum propagation distance of the plasma during the increase in the voltage is shorter than the case of 8 kV. On the other hand, the propagation of the plasma is more like a plasma bullet at the rising edge of the applied voltage.
For the applied voltage of −8 kV, the time-resolved electric field at 1 cm away from the nozzle (z = 1 cm) obtained by E-FISH method is shown in figure 9. The electric field increases quickly when the applied voltage changed from 0 to −8 kV. The electric field reaches its peak of −16 kV cm−1 at t1 = 85 ns. Then it decreases slowly to about −6 kV cm−1 and stays as a constant until 1000 ns. Afterwards, with the drop of the applied voltage from −8 kV to zero, the electric field quickly increases to a peak of about 20 kV cm−1. Subsequently it decreases to zero after about 500 ns. It is worthy to point out that the peak electric field at t2 during the secondary discharge is significantly higher than that at t1 during the first discharge.
According to figures 7 and 9, it is found that the peak electric field is about 20–25 kV cm−1 at the edge when the applied voltage changes from 0 to 8 kV or −8 to 0 kV. In addition, the plasma appears like a bullet during these periods. On the other hand, at the edge when the applied voltage changes from 8 to 0 kV or 0 to −8 kV, the peak electric field is about −14 to −16 kV cm−1. The plasma propagates much shorter distance and it has a relatively bright long tail.
In the next part, the effects of the pulse width and the pulse rising time on the electric field of the plasma bullet are measured. Firstly, the pulse width is reduced to 500 ns. All of the other parameters remain unchanged. The applied voltage waveform along with the measured electric field is shown in figure 10, which is similar to the results shown in figure 7. It reaches a peak value of about 21 kV cm−1 at 85 ns, and there is also a secondary discharge at the falling edge with a peak value of about −15 kV cm−1.
Secondly, the pulse rising time is increased from 60 to 200 ns. The applied voltage waveform and the obtained electric field at 1 cm away from the nozzle (z = 1 cm) by E-FISH method are shown in figure 11. It can be seen that the electric field increases slowly and reaches a peak of about 10.5 kV cm−1 at about 350 ns, which is much weaker than the case of 60 ns pulse rising time. Accordingly, the electric field induced by the secondary discharge at the falling edge is also much weaker, which is about −8 kV cm−1. Compared to the rapid falling edge of the voltage, the longer falling time makes discharge occur in the process of the voltage drop. This phenomenon also happens at the rising edge of the voltage. However, in figure 11, due to the imperfect rising edge of the voltage, the voltage still slowly rises after 200 ns, which causes the time of the electric field reaching its peak to be delayed to about 350 ns.
When APPJs are used for various applications, a trace amount of nitrogen or oxygen is often added to the working gas, which can increase the concentration of reactive nitrogen species (RNS) or reactive oxygen species (ROS) [38–41]. However, how the trace amount of nitrogen or oxygen affects the electric field and thus the propagation of the plasma jets has not been investigated. Note that the hyperpolarizability between helium and oxygen (or nitrogen) is significantly different; only a tiny amount of oxygen or nitrogen can be mixed into helium in this experiment. Moreover, in order to reduce the mixing of ambient air, sufficient helium flow rate is required and the measurement position should be close to the nozzle. Corresponding calibration is carried out under the same experimental conditions, which can be found in figure 4. Although the hyperpolarizability of argon is much closer to that of air, the helium jets are more widely used in applications, and the propagation of the helium plasma is like a bullet but that of argon plasma is like filaments; a more regular shape is convenient for simulation and correction [42].
Firstly, 0%, 0.5% and 1.0% oxygen are added into helium with a total gas flow rate of 2 l min−1. All the other parameters are kept the same. The electric field at 5 mm away from the nozzle (z = 5 mm) is measured by the E-FISH method as shown in figure 12. It is clearly shown that, when 0.5% of oxygen is added, the peak electric field at the rising edge reduced from 22 to 15 kV cm−1. Moreover, the electric field decreases much faster and it stays at a value of about 2 kV cm−1 lower than in the case of pure helium before the falling edge of the applied voltage. When 1% of oxygen is added, the peak electric field becomes even weaker.
Secondly, 0%, 0.5%, 1.0% and 1.5% nitrogen are added to the working gas and the corresponding electric field is measured by the E-FISH method as shown in figure 13. It can be seen that, when 0.5% nitrogen is added, the peak electric field actually slightly increases, especially at the falling edge of the voltage pulse, and it increases reversely from −12 to −16 kV cm−1. However, when 1.5% nitrogen is added, the peak electric field decreases from 24 to 19 kV cm−1 at the rising edge of the voltage pulse. Moreover, it is found that the electric field decreases much faster when a trace amount of nitrogen is added compared with the pure helium case.
We know that the higher the electric field is, the faster the plasma bullet propagates. To see if this is true, high-speed photographs of the plasma for He only, He + 0.5%N2, and He + 1.5%O2 are captured by an ICCD camera with exposure time of 5 ns. The obtained results are shown in figure 14. For the case of He only, the plasma already reaches 1.5 cm (red dotted line) out of the nozzle at 90 ns. For the case of He + 0.5%N2, the plasma propagates faster and reaches 1.5 cm at about 70 ns. When the working gas is changed to He + 1.5%O2, it does not reach 1.5 cm at 90 ns and it propagates even more slowly later.
In this work, the electric field induced second harmonic (E-FISH) method is applied to diagnose the electric field of a helium plasma jet. Time-resolved electric fields under different conditions are investigated, including the effects of the applied voltage polarity, voltage pulse rising time and working gases percentage on the electric field of the plasma jet and its dark channel.
When positive pulse voltage is applied, the electric field has a peak of about 25 kV cm−1 at the rising edge of the voltage pulse. A dark channel is left behind the plasma bullet and the electric field in the dark channel is about 5 kV cm−1. At the falling edge of the voltage pulse, the electric field increases reversely to a peak of −14 kV cm−1. On the other hand, when negative pulse voltage is applied, the electric field has a peak of −16 kV cm−1 when the voltage is increased to −8 kV. A relatively bright channel is left behind the plasma head and the electric field in this relatively bright channel is about −6 kV cm−1. When the applied voltage changed from −8 kV to zero, the electric field reaches a peak of about 20 kV cm−1, which is higher than that of the primary discharge. It is worth emphasizing that the electric field in the dark channel when positive voltage pulse is used is lower than that in the relatively bright channel when negative voltage pulse is used, which is the reason why one is dark and the other one is relatively bright. Because low electric field indicates low electron temperature, a high-energy electron is not enough to excite adequate species, which leads to a relatively dark channel.
When the pulse width is reduced from 1 μs to 500 ns, it does not affect the peak electric field of the plasma bullet. On the other hand, when the pulse rising time increases from 60 to 200 ns, the peak electric fields at both the rising edge and the falling edge of the voltage decrease significantly to about 10.5 kV cm−1 and −8 kV cm−1, respectively.
The percentage of gas mixture also affects the electric field of the plasma bullet. When 0.5% of oxygen is added to the main working gas helium, the peak electric field at the rising edge is only about 15 kV cm−1. Moreover, the electric field decreases much faster and stays at a value of about 2 kV cm−1 lower than the case of pure helium before the falling edge of the applied voltage. On the other hand, when 0.5% nitrogen is added, the peak electric field actually increases, especially at the falling edge of the voltage pulse, where it increases reversely from −12 to −16 kV cm−1. However, when 1.5% nitrogen is added, the peak electric field decreases from 24 to 19 kV cm−1 at the rising edge of the voltage pulse. Note that the influence of the temporal electric field profile is not considered in this work; more modification can be carried out for future diagnostics [42].
In this work, the electric field induced second harmonic (E-FISH) method is applied to diagnose the electric field of a helium plasma jet. Time-resolved electric fields under different conditions are investigated, including the effects of the applied voltage polarity, voltage pulse rising time and working gases percentage on the electric field of the plasma jet and its dark channel.
When positive pulse voltage is applied, the electric field has a peak of about 25 kV cm−1 at the rising edge of the voltage pulse. A dark channel is left behind the plasma bullet and the electric field in the dark channel is about 5 kV cm−1. At the falling edge of the voltage pulse, the electric field increases reversely to a peak of −14 kV cm−1. On the other hand, when negative pulse voltage is applied, the electric field has a peak of −16 kV cm−1 when the voltage is increased to −8 kV. A relatively bright channel is left behind the plasma head and the electric field in this relatively bright channel is about −6 kV cm−1. When the applied voltage changed from −8 kV to zero, the electric field reaches a peak of about 20 kV cm−1, which is higher than that of the primary discharge. It is worth emphasizing that the electric field in the dark channel when positive voltage pulse is used is lower than that in the relatively bright channel when negative voltage pulse is used, which is the reason why one is dark and the other one is relatively bright. Because low electric field indicates low electron temperature, a high-energy electron is not enough to excite adequate species, which leads to a relatively dark channel.
When the pulse width is reduced from 1 μs to 500 ns, it does not affect the peak electric field of the plasma bullet. On the other hand, when the pulse rising time increases from 60 to 200 ns, the peak electric fields at both the rising edge and the falling edge of the voltage decrease significantly to about 10.5 kV cm−1 and −8 kV cm−1, respectively.
The percentage of gas mixture also affects the electric field of the plasma bullet. When 0.5% of oxygen is added to the main working gas helium, the peak electric field at the rising edge is only about 15 kV cm−1. Moreover, the electric field decreases much faster and stays at a value of about 2 kV cm−1 lower than the case of pure helium before the falling edge of the applied voltage. On the other hand, when 0.5% nitrogen is added, the peak electric field actually increases, especially at the falling edge of the voltage pulse, where it increases reversely from −12 to −16 kV cm−1. However, when 1.5% nitrogen is added, the peak electric field decreases from 24 to 19 kV cm−1 at the rising edge of the voltage pulse. Note that the influence of the temporal electric field profile is not considered in this work; more modification can be carried out for future diagnostics [42].
This work was supported by National Natural Science Foundation of China (Nos. 62171355, 61875156), the 111 Project (No. B17035), the Natural Science Basic Research Program of Shaanxi Province (No. 2021JM-135), and the Stable Support Project of Basic Scientific Research Institutes (Nos. A131901W14, A132001W12).
[1] |
Chen W et al 2016 IEEE Trans. Plasma Sci. 44 3235 doi: 10.1109/TPS.2016.2617680
|
[2] |
Guo L X et al 2018 Phys. Plasmas 25053707 doi: 10.1063/1.5020628
|
[3] |
Yao B et al 2017 Phys. Plasmas 24102104 doi: 10.1063/1.4998302
|
[4] |
Yao B et al 2019 IEEE IEEE Trans. Veh. Technol. 68 4328 doi: 10.1109/TVT.2019.2902962
|
[5] |
Zhang D M et al 2022 Infrared Phys. Technol. 122104054 doi: 10.1016/j.infrared.2022.104054
|
[6] |
Liu S H and Guo L X 2016 IEEE Trans. Plasma Sci. 44 2838 doi: 10.1109/TPS.2016.2612003
|
[7] |
Niu Z Y, Su T and Zhang J B 2019 IEEE Access 7 7678 doi: 10.1109/ACCESS.2018.2886243
|
[8] |
Liu S H 2018 Phys. Plasmas 25062105 doi: 10.1063/1.5030194
|
[9] |
Zhou R 2017 IET Radar Sonar Navig. 11 1235 doi: 10.1049/iet-rsn.2016.0464
|
[10] |
Bian Z et al 2021 IEEE Trans. Geosci. Remote Sens. 602002413 doi: 10.1109/TGRS.2021.3112801
|
[11] |
Zhang J et al 2022 IET Radar Sonar Navig. 16 1162 doi: 10.1049/rsn2.12251
|
[12] |
Schwartz R et al 2012 Remote infrared imaging of the space shuttle during hypersonic flight: HYTHIRM mission operations and coordination In: Proc. of the 42nd AIAA Thermophysics Conf. Honolulu: AIAA
|
[13] |
Niu Q L et al 2019 Infrared Phys. Technol. 99 28 doi: 10.1016/j.infrared.2019.04.005
|
[14] |
Niu Q L, He Z H and Dong S K 2017 Chin. J. Aeronaut. 30 1101 doi: 10.1016/j.cja.2017.04.003
|
[15] |
Zupanc F J and Degroot W A 1993 Propul. Eng. Res. Center Pennsylvania State Univ., NASA
|
[16] |
Thidé B et al 1983 Radio Sci. 18 851 doi: 10.1029/RS018i006p00851
|
[17] |
Leyser T B 2001 Space Sci. Rev. 98 223 doi: 10.1023/A:1013875603938
|
[18] |
Stenflo L 2004 Phys. Scr. 4 262 doi: 10.1238/Physica.Topical.107a00262
|
[19] |
Norin L et al 2009 Phys. Rev. Lett. 102065003 doi: 10.1103/PhysRevLett.102.065003
|
[20] |
Leyser T B et al 2009 Phys. Rev. Lett. 102065004 doi: 10.1103/PhysRevLett.102.065004
|
[21] |
Sergeev E N et al 2016 Radio Sci. 51 1118 doi: 10.1002/2015RS005936
|
[22] |
Chen F F 1984 Introduction to Plasma Physics and Controlled Fusion 2nd edn (New York: Plenum Press)
|
[23] |
Sun X et al 2018 Calculation of characteristic parameters of solid rocket exhaust plume plasma In: Proc. of the 2018 9th Int. Conf. on Mechanical and Aerospace Engineering Budapest: IEEE
|
[24] |
Takahashi Y 2016 J. Phys. D: Appl. Phys. 49015201 doi: 10.1088/0022-3727/49/1/015201
|
[25] |
Fu H Y and Scales W A 2018 Kinetic modeling of stimulated electromagnetic emissions during ionospheric heating experiment In: Proc. of the 2018 12th Int. Symp. on Antennas, Propagation and EM Theory Hangzhou: IEEE
|
[26] |
Barati H and Habibi M 2020 IEEE Trans. Plasma Sci. 48 3862 doi: 10.1109/TPS.2020.3026451
|
[27] |
Dehghan M and Abbaszadeh M 2017 Eng. Comput. 33 961 doi: 10.1007/s00366-017-0509-y
|
[28] |
Strozzi D J et al 2006 J. Plasma Phys. 72 1299 doi: 10.1017/S0022377806005599
|
[29] |
Scales W A, Cheng K T and Srivastava S 1997 J. Atmos. Solar Terr. Phys. 59 2373 doi: 10.1016/S1364-6826(96)00131-9
|
[30] |
Kolber T, Rozmus W and Tikhonchuk V T 1995 Phys. Plasmas 2 256 doi: 10.1063/1.871446
|
[31] |
Almutairi A et al 2022 Results Phys. 37105521 doi: 10.1016/j.rinp.2022.105521
|
[32] |
Sary G and Gremillet L 2022 Phys. Plasmas 29072103 doi: 10.1063/5.0090211
|
[33] |
Çınar M, Onder I and Secer A 2022 Electronic Research Archive 30 335 doi: 10.3934/era.2022018
|
[34] |
Jie T and Jie J 2011 Stress limiter consideration for k-omega turbulence models in shock-wave/turbulent boundary-layer interactions in supersonic and hypersonic flows In: Proc. of the 20th AIAA Computational Fluid Dynamics Conf. Honolulu: AIAA
|
[35] |
Ichimaru S 1973 Basic Principles of Plasma Physics (Boca Raton: CRC Press)
|
[36] |
Hanssen A et al 1992 J. Geophys. Res. 97 12073 doi: 10.1029/92JA00874
|
[37] |
Zakharov V E and Rubenchik A M 1972 J. Appl. Mech. Tech. Phys. 13 669 doi: 10.1007/BF00850881
|
[38] |
Sprague R A 1995 A numerical study of plasma turbulence with application to ionospheric heating PhD Thesis University of California (San Diego, USA)
|
[39] |
Mjølhus E, Hanssen A and DuBois D F 1995 J. Geophys. Res. 100 17527 doi: 10.1029/95JA01158
|
[1] | Yan YANG, Peiyu JI, Maoyang LI, Yaowei YU, Jianjun HUANG, Bin YU, Xuemei WU, Tianyuan HUANG. The application of a helicon plasma source in reactive sputter deposition of tungsten nitride thin films[J]. Plasma Science and Technology, 2022, 24(6): 065503. DOI: 10.1088/2058-6272/ac5c27 |
[2] | Min ZHU (朱敏), Chao YE (叶超), Xiangying WANG (王响英), Amin JIANG (蒋阿敏), Su ZHANG (张苏). Effect of radio-frequency substrate bias on ion properties and sputtering behavior of 2 MHz magnetron sputtering[J]. Plasma Science and Technology, 2019, 21(1): 15507-015507. DOI: 10.1088/2058-6272/aae7dd |
[3] | Gao ZHAO (赵高), Wanying ZHU (朱婉莹), Huihui WANG (王慧慧), Qiang CHEN (陈强), Chang TAN (谭畅), Jiting OUYANG (欧阳吉庭). Study of axial double layer in helicon plasma by optical emission spectroscopy and simple probe[J]. Plasma Science and Technology, 2018, 20(7): 75402-075402. DOI: 10.1088/2058-6272/aab4f1 |
[4] | Ling ZHANG (张玲), Guo CHEN (陈果), Zhibing HE (何智兵), Xing AI (艾星), Jinglin HUANG (黄景林), Lei LIU (刘磊), Yongjian TANG (唐永建), Xiaoshan HE (何小珊). Investigation of working pressure on the surface roughness controlling technology of glow discharge polymer films based on the diagnosed plasma[J]. Plasma Science and Technology, 2017, 19(7): 75505-075505. DOI: 10.1088/2058-6272/aa6618 |
[5] | Liang SONG (宋亮), Xianping WANG (王先平), Le WANG (王乐), Ying ZHANG (张营), Wang LIU (刘旺), Weibing JIANG (蒋卫斌), Tao ZHANG (张涛), Qianfeng FANG (方前锋), Changsong LIU (刘长松). Fabrication and characterization of He-charged ODS-FeCrNi films deposited by a radio-frequency plasma magnetron sputtering technique[J]. Plasma Science and Technology, 2017, 19(4): 45502-045502. DOI: 10.1088/2058-6272/aa57f0 |
[6] | WU Zhonghang (吴忠航), LIANG Rongqing (梁荣庆), Masaaki NAGATSU (永津雅章), CHANG Xijiang (昌锡江). The Characteristics of Columniform Surface Wave Plasma Excited Around a Quartz Rod by 2.45 GHz Microwaves[J]. Plasma Science and Technology, 2016, 18(10): 987-991. DOI: 10.1088/1009-0630/18/10/04 |
[7] | Umm-i-KALSOOM, R. AHMAD, Nisar ALI, I. A. KHAN, Sehrish SALEEM, Uzma IKHLAQ, et al. Effect of Power and Nitrogen Content on the Deposition of CrN Films by Using Pulsed DC Magnetron Sputtering Plasma[J]. Plasma Science and Technology, 2013, 15(7): 666-672. DOI: 10.1088/1009-0630/15/7/12 |
[8] | LU Wenqi (陆文琪), JIANG Xiangzhan (蒋相站), LIU Yongxin (刘永新), YANG Shuo (杨烁), et al. Improved Double-Probe Technique for Spatially Resolved Diagnosis of Dual-Frequency Capacitive Plasmas[J]. Plasma Science and Technology, 2013, 15(6): 511-515. DOI: 10.1088/1009-0630/15/6/05 |
[9] | PANG Jianhua (庞见华), LU Wenqi (陆文琪), XIN Yu (辛煜), WANG Hanghang (王行行), HE Jia (贺佳), XU Jun (徐军). Plasma Diagnosis for Microwave ECR Plasma Enhanced Sputtering Deposition of DLC Films[J]. Plasma Science and Technology, 2012, 14(2): 172-176. DOI: 10.1088/1009-0630/14/2/17 |
[10] | RU Lili (汝丽丽), HUANG Jianjun (黄建军), GAO Liang (高亮), QI Bing (齐冰). Influence of Microwave Power on the Properties of Hydrogenated Diamond-Like Carbon Films Prepared by ECR Plasma Enhanced DC Magnetron Sputtering[J]. Plasma Science and Technology, 2010, 12(5): 551-555. |
1. | Maciel, C.V.T., da Silva, J.F., da Silva, R.F. et al. Intrinsic heating in NdxY1.00-xAl3(BO3)4 particles excited at 808 nm leads to bright multi-band up conversion emission via ladder-thermal excitation. Journal of Alloys and Compounds, 2025. DOI:10.1016/j.jallcom.2024.178119 |
2. | Abdel-Fattah, E.M.. Plasmonic ZnO-Au Nanocomposites: A Synergistic Approach to Enhanced Photocatalytic Activity through Nonthermal Plasma-Assisted Synthesis. Crystals, 2024, 14(10): 890. DOI:10.3390/cryst14100890 |
3. | Heng, Y., Yu, L., Chen, Y. et al. Plasma-Assisted Material Preparation Strategies and Property Optimization. Physica Status Solidi (A) Applications and Materials Science, 2024. DOI:10.1002/pssa.202400702 |