^{1}

^{1}

^{2}

The interaction mechanism between the plasma and liquid is a key problem for the electrothermal chemical launch technology. To investigate this problem, a simulated experiment for the expansion process of a plasma jet in the working fluid is carried on. Based on this experiment, a two-dimensional axisymmetric unsteady theoretical model is established to reveal the plasma-liquid interaction flow field pattern. The results show that a typical Taylor cavity forms as the plasma jet expands in liquid. The induction effect of the stepped-wall structure enhances the radial expansion of the plasma jet. An arc-shaped pressure wave is generated at the front of the plasma jet and then evolves into the plane wave. A high-pressure area forms at the head of the plasma jet and then moves downstream. There is a strong plasma-liquid turbulent mixing at the interface, especially near the steps and the nozzle exit area.

Electrothermal chemical launch technology, which is an advanced high-velocity propulsion technology, is widely interested by scholars [

Electrothermal chemical guns can be classified into solid propellant electrothermal chemical guns [

The concept of using a stepped-wall structure in the liquid chamber was introduced as a potential mean of controlling the combustion of bulk-loaded liquid propellant [

In this work, the evolution process of a plasma jet in the stepped-wall liquid chamber with four steps was observed by the experiment. Based on the experiment, a two-dimensional axisymmetric unsteady mathematical physics model of a plasma jet expanding in water was established to study the expansion characteristics of a plasma jet, as well as the pressure evolution in the plasma-liquid interaction flow field. The results can provide a reference to investigate the interaction mechanism of the plasma jet and the ambient liquid.

To study the interaction mechanism of the plasma jet and the working fluid, an experimental device was designed including a pulse power source, a plasma generator, and liquid chamber. The pulse power source mainly included a pulse-forming network, as shown in

The plasma generator (shown in

The anode was sealed and the cathode was connected to the nozzle bottom of liquid chamber. This liquid chamber has four steps and the size details are shown in

During the experiment process, the pulse power source loaded electric energy between the anode and cathode of the plasma generator. The exploding wire in the plasma generator exploded, and the high-temperature metal plasma was generated. Then the metal plasma ablated the polyethylene on the capillary wall, and the plasma mixture is generated. When the pressure in the nozzle was high enough to break the copper film, the plasma jet was injected into the stepped-wall liquid chamber. The expansion process of a plasma jet in water, shown as

A plasma cavity, called Taylor cavity, forms within the working fluid inside the liquid chamber. This cavity’s axial expansion goes faster than its radial expansion as the plasma cavity expands along the step structure in the liquid chamber. These images show that the cavity boundary is blurred, especially at the later stage of the expansion process, because there is strong turbulent mixing

Structure | Length/mm | Diameter/mm | Diameter Increment/mm |
---|---|---|---|

step1 | 20 | 18 | / |

step2 | 20 | 30 | 12 |

step3 | 20 | 42 | 12 |

step4 | 38 | 54 | 12 |

of plasma and water at the interface. Further, this turbulence is stronger near the step corner than other areas, because the stepped-wall structure causes the reflow liquid which enhancing the plasma-liquid interaction at the step corner. In addition, the plasma cavity is brighter during the period of 2 - 4 milliseconds, indicating the thermal energy of the plasma jet is greater at this period. Certainly, the shading effect caused by many bubbles and vapors should be not ignored in these graphs.

An unsteady 2D theoretical model was established to simulate the expansion process of a plasma jet in liquid. Some physical assumptions are made as follows:

1) The Volume of Fluid model calculates the multiphase flow and the Standard k-ε model simulates the plasma-liquid turbulent mixing.

2) The plasma-liquid interaction process studied in this work is short (not more than 10 milliseconds), so the liquid vaporization is ignored. The effect of the plasma sheath, as well as the effect of gravity, is also ignored to simplify the calculation.

3) The plasma jet is approximated as a high-temperature and high-pressure compressible fluid.

Based on the assumptions, governing equations are described as following:

1) Equation of continuity:

∂ ( α q ) ∂ t + u ⋅ ∇ α q = 0 (1)

where u is velocity, α_{q} is the volume fraction of qth phase. The volume fraction of the primary phase (plasma), α_{1}, will be calculated based on the constraint: α 1 + α 2 + α 3 = 1 .

2) Momentum equation:

∂ ( ρ u ) ∂ t + ∇ ⋅ ( ρ u u ) = − ∇ p + ∇ ⋅ [ μ ( ∇ u + ∇ u T ) ] (2)

where ρ is density, p is static pressure and μ is dynamic viscosity. u^{T} is the transposed matrix of u.

3) Energy equation:

∂ ( ρ E ) ∂ t + ∇ ⋅ ( u ( ρ E + p ) ) = ∇ ⋅ ( k eff ∇ T ) + S h (3)

where k_{eff} is effective thermal conductivity, E is energy and T is temperature. The radiant energy S_{h} is given by

S h = f σ ( T − T 0 ) 4 (4)

where f is emissivity, σ is Stefan-Boltzmann constant and T_{0} is environment temperature.

4) State equation:

p = ρ T R (5)

where R is mixing gas constant.

5) Turbulence equation:

∂ ∂ t ( ρ k ) + ∂ ∂ x ( ρ k u x ) = ∂ ∂ r [ ( μ + μ t σ k ) ∂ k ∂ r ] + G k + G b − ρ ε − Y M (6)

∂ ∂ t ( ρ ε ) + ∂ ∂ x ( ρ ε u x ) = ∂ ∂ r [ ( μ + μ t σ ε ) ∂ ε ∂ r ] + C 1ε ε k ( G k + C 3ε G b ) − C 2ε ρ ε 2 k (7)

where k is the turbulence kinetic energy, and ε is its dissipation rate, σ_{k} and σ_{ε} are the turbulent Prandtl numbers for k and ε. In the equation, μ_{t} = 0.09ρk^{2}/ε, C_{1ε}, C_{2ε}, and C_{3ε} are constants, x is axial displacement, r is radial displacement, u_{x} is axial velocity. G_{k} represents the generation of turbulence kinetic energy due to the mean velocity gradients, G_{b} is the generation of turbulence kinetic energy due to buoyancy, Y_{M} represents the contribution of the fluctuating dilatation during the compressible turbulence to the overall dissipation rate.

The axial symmetry of the theoretical model enables the computational domain adopted half of the plasma-liquid interaction flow field. The computational domain, as shown in

environment. The geometric size of the liquid chamber is the same as the liquid chamber in experimental devices. The nozzle is 10 mm in length and 1 mm in radius, and the atmospheric environment is 100 mm in length and 81 mm in radius.

The computational domain is divided into the rectangular grid, and the gradient mesh is adopted in the atmospheric environment to shorten calculation time. The mesh has 90 thousand cells, and mesh independence has been verified.

In the computational domain, the left exit of the nozzle is defined as the inlet boundary whose pressure is determined by the experiment, the temperature is 5000 K. The pressure outlet parameters are the same as atmospheric parameters. The no-slip heat insulation condition is adopted at the wall, and the standard wall function method is used near the wall to deal with turbulence. At the initial time, the plasma jet has not been ejected into the liquid chamber filled with water, so the pressure of the computational domain is 101,325 Pa and the temperature is 300 K.

Expansion process of a plasma jet in the stepped-wall liquid chamber filled with water was analyzed by the numerical calculation. The phase distribution of the plasma-water interaction flow field is shown as

We defined the maximum axial position where the plasma volume fraction is more than 90% as the axial penetration distance of plasma jet. The experimental and calculated time histories of the axial penetration distance of plasma jet are

shown in

The pressure distribution of plasma-liquid interaction flow field is shown in

The evolution process of a plasma jet in the stepped-wall liquid chamber with four steps was observed by the experiment. Based on the experiment, a two-di- mensional axisymmetric unsteady mathematical physics model of a plasma jet expanding in water was established to study the plasma-liquid interaction mechanism. The calculated results agree well with the experimental results. The conclusions are obtained as follows.

1) As the plasma jet expands in water, a plasma cavity, which forms within the water inside the liquid chamber, expands along the stepped wall. This cavity’s axial expansion goes faster than its radial expansion. The radial induction effect of the stepped-wall structure enhances the radial expansion of the plasma jet. The plasma cavity is brighter during the period of 2 - 4 milliseconds in the sampling period. Besides, there is strong plasma-liquid turbulent mixing at the interface, especially near the steps and nozzle exit.

2) The plasma-liquid interaction has an effect on the pressure distribution in the flow field. An arc-shaped pressure wave is generated at the front of the plasma jet and then evolves into the plane wave. A high-pressure area forms at the head of the plasma jet and then moves downstream. Besides, there is a local low-pressure area at the corner of the step.

This work is greatly supported by the National Natural Science Foundation of China (Grant No. 51506094).

The authors declare no conflicts of interest regarding the publication of this paper.

Liu, Y., Yu, Y.G. and Mang, S.S. (2019) Expansion Characteristics of a Plasma Jet in the Stepped-Wall Chamber Filled with Water. World Journal of Engineering and Technology, 7, 293-301. https://doi.org/10.4236/wjet.2019.72020