## Abstract

The chapter discusses the evolution of Rayleigh-Taylor instability (RTI) in ordinary fluids and in a plasma fluid. RT instability exits in many situations from overturn of the outer portion of the collapsed core of a massive star to laser implosion of deuterium-tritium fusion targets. In the mixture of fluids, the instability is triggered by the gravitational force acting on an inverted density gradient. The motivation behind the study of the instability has been explored by discussing the applications of RT instability. The basic magnetohydrodynamics equations are used to derive the dispersion relation (for an ordinary fluid and plasmas) for two fluids of unequal densities. The conditions of the growth rate of the instability and the propagating modes are obtained by linearizing the fluid equations. The perturbed potential is found to increase with the plasma parameters in a Hall thruster.

### Keywords

- instabilities
- plasma
- Navier-Stokes
- growth rate
- Hall thruster

## 1. Introduction

Flow instabilities are used to increase the heat and mass transfer rates as well as to fuse the fluids of dissimilar properties (viscosity, elasticity, density, etc.). In other technological applications, these instabilities are accountable to unstable the multilayer and free-surface flows. Multilayer flows are used in coating processes and lubricated pipeline transport. The presence of the instabilities in the system leads to nonuniform film thickness and defects, where good optical finishing and smooth edges are required by the industry, which further leads to poor product quality. Suppression of these instabilities has been a major task from a long time by the researchers to improve the product quality [1, 2]. Rayleigh-Taylor (RT) instability takes place when a lighter fluid supports a heavy fluid, then any perturbation of the interface grows and leads to spikes of the heavier fluid penetrating into the lighter one and the interface becomes unstable. The contact discontinuity between the two fluids is unstable to perturbations that grow by converting potential energy to kinetic energy, causing bubbles of the low-density fluid to rise, and spikes of the high-density fluid to sink. If the light fluid is above the heavy fluid, the interface is stable. In a magnetized plasma, the Rayleigh-Taylor instability can occur because the magnetic field acts as a light fluid supporting a heavy fluid (the plasma).

In curved magnetic fields, the centrifugal force on the plasma due to the charged particle motion along the curved field lines acts as an equivalent gravity force. When forces associated with the density gradient and gravity oppose each other, the RT instability sets in [3, 4]. The box of fluid shown in Figure 1 is now filled with two incompressible fluids of differing densities, separated by an interface with a perturbation imposed as shown in Figure 1 . Here, RTI is seen to play a wider role in many branches of science from astrophysical systems to industries.

## 2. Review of status of research

This instability occurs in many interesting physical situations, such as implosion of inertial confinement fusion capsules, core collapse of supernovae, or electromagnetic implosions of metal liners. The Rayleigh-Taylor problem was first studied by Lord Rayleigh in 1883 and Sir G.I. Taylor in 1950 [3]. Taylor used the theory of linearization for the small oscillations at the interface and obtained an exponential growth rate. Chandrasekhar, in 1961, studied the magnetic field case analytically for the fluids that are incompressible, inviscid, and have zero resistivity. Qin et al. [5] reported the synthesis of chains of metal nanoparticles with well-controlled particle sizes and spacing induced by the Rayleigh instability. Bychkov et al. [6] derived the dispersion relation for the internal waves and the RT instability in a nonuniform unmagnetized quantum plasma with a constant gravitational field. They have shown that the quantum effects always play a stabilizing role for the RT wave instability. Cao et al. [7] studied the RT instability incorporating the quantum magnetohydrodynamic equations and solved the second-order differential equation under different boundary conditions with quantum effects. Khomenko et al. [8] modeled the growth rate of the instability and the evolution of velocity and magnetic field vector in the prominence plasma (closer to Sun’s surface) under the presence of neutral atoms. Diaz et al. [9] derived the criterion for the growth rate of the RT instability in partially ionized plasma using single fluid theory. Ibrahim and Marshall theoretically investigated the impact of velocity profile on RTI within the jet to examine the effects of its relaxation on intact length [10]. Carlyle and Hillier experimentally verified that stronger magnetic fields can suppress the growth of the rising bubbles of the RTI [11]. Litvak and Fisch derived the necessary instability conditions of azimuthally propagating perturbations in a Hall thruster plasma [12]. Recently, investigators derived the dispersion for the Rayleigh-Taylor instabilities in a Hall thruster using the two - fluid theory [13, 14]. Shorbagy and Shukla investigated the RT instability in a nonuniform multi-ion plasma in a Hall thruster to obtain the growth rate of the instability [15]. Ali et al. [16] derived the modified dispersion relation for the Rayleigh-Taylor instability under the quantum corrections incorporating the terms of Fermi pressure and the Bohm potential force.

## 3. Basic fluid equations and Bernoulli’s theorem

First, we consider the two simple fluids separated by a smooth interface to derive the dispersion relation. Let us assume that in each separate region, the density is constant. The coordinate *x* is in the horizontal, *z* in the vertical, and *y* is going into the page. We consider a flow in the *x*-direction, which in the lower half-space

Here, we have used total time derivative. Partial time derivative keeps an eye on a point and represents the rate of velocity change at that point. Total time derivative keeps an eye on fluid element and measures its velocities at

Let us consider the fluid is inviscid, so that we take viscosity

If the density remains constant in one region, we can write Eq. (3) as

Now integrating the above equation in horizontal and vertical directions, we get unsteady equation for the Bernoulli theorem.

That is, the total mechanical energy of the moving fluid comprising the gravitational potential energy of elevation, the energy associated with the fluid pressure, and the kinetic energy of the fluid motion remains constant.

Let

For the incompressible fluid, Eq. (1) yields that

The Bernoulli theorem state that quintiles

To understand the interface, we must impose boundary conditions. First of all the vertical velocities of the fluids must match with the interface, so we impose the kinematic boundary condition. Now we need to introduce the location of the interface by assigning variable *z*-direction. In addition, at the interface, the velocity of both fluids must be continuous.

Let us say that the pressure is continuous along the interface, that is

## 4. Asymptotic boundary conditions at far field

We are looking for changes only on the interface at

## 5. Linear analysis

Eq. (9) contains a nonlinear term

## 6. Eigenvalue solution

Let us consider all the perturbed variables

The general solution of Eq. (16) is written as

The above two solutions must satisfy the boundary conditions such that

This solution further leads to the following form in Fourier mode,

Here

Then Eq. (20) changes into the form

Eq. (15) gives

Using Eq. (21) in Eq. (23) results in

Since the perturbed quantity

Eq. (25) contains complete information about the linear stability of the two superposed fluid layers of different densities. The Atwood number

### 6.1 First case: capillary-gravity waves (ρ 2 = 0 )

It is classical dispersion relation for gravity-capillary waves in deep water [17, 18]. These are also called short gravity waves. In this category the longer waves travel faster. Any initial disturbance may be regarded as the superposition of waves of a broad spectrum of lengths. The above relation then says that waves of different lengths will eventually separate, that is, disperse. This phenomenon is called dispersion, hence above relations are also known as the dispersion relation.

### 6.2 Second case: propagating modes (ρ 1 > ρ 2 )

If the lighter fluid is supported by heavier fluid, that is,

Figure 2 shows the variations of phase velocity of RT instability with (a) density ratio and (b) wave number respectively.

### 6.3 Third case: Rayleigh-Taylor instability (ρ 2 > ρ 1 )

The frequency of oscillations will be negative imaginary and unstable if

Substituting the value of

The term

## 7. Rayleigh-Taylor instability in plasma thruster

In the previous section, the general idea of RT instability has been explored. Here we have derived the RT equation for a plasma fluid using two fluid theory. In a Hall thruster, the propellant (plasma) is ionized and then accelerated by electrostatic forces. It has high thrust resolution, so it is best suited for the adjustment of the location of the satellite onboard [19, 20, 21, 22, 23, 24, 25, 26, 27]. Let us consider a plasma with nonuniform density confined under the crossed electric and magnetic fields.

Figure 4 shows the typical diagram of a Hall plasma thruster [26]. RT instability is common in Hall thrusters. Studies show that Rayleigh instability is driven by the presence of gradients in axial density, magnetic field, and velocity of the plasma species. Here we deduce a Rayleigh equation under the presence of ion temperature and check the variations of perturbed potential with plasma parameters.

### 7.1 Theoretical model for RTI in plasma

We consider plasma comprising of ions and electrons immersed in a magnetic field *x-*axis (axis of the thruster) and the magnetic field *z*-axis (along the radius of the thruster). Hence, the azimuthal dimension is along the *y-*axis. We use

We use the linearized form of the above equations for small perturbations of the ion and electron densities, their velocities, and electric field. We write perturbed densities of ions (electrons) by *y*-direction. The unperturbed density (electric field) is taken as

The unperturbed ions’ velocity

Eq. (38) provides the velocity components of electron

In the above equations, the coordinate

The electron cyclotron frequency is almost

The electron continuity equation gives the perturbed electron density

The plasma frequency of oscillations for ion (electron) is defined as

Using Eqs. (39) and (46) in Eq. (48) gives the perturbed potential in the following form:

In the case of high frequency of oscillations and in the absence of ion thermal pressure, Eq. (49) turns into Rayleigh’s equation of fluid dynamics as below

Here

**Resonance condition for the RT instability**

From Eq. (49), it is clear that propagating mode may lead to instability if parameter

### 7.2 Variations of perturbed potential

The RT Eq. (49) is solved numerically for the perturbed potential

Figure 5 shows the variation of the perturbed potential with the magnetic field and it has been observed that the potential increases with the increasing magnetic field. These results are consistent with Keidar and Boyd model [28] and that other investigators [13, 14] for the potential of plasma plume. This situation is correspond to the plasma jet enters a transverse magnetic field with a high velocity under the condition that the magnetic field is relatively weak so that only the electrons are magnetized whereas the ions move out of the effect of magnetic field. However, ambipolar (both electrons and ions moving in opposite directions) plasma flow across the magnetic field may require an electric field to appear under the above conditions. Therefore, we can expect the potential to increase across the magnetic field.

The perturbed potential gets increased with the higher value of electron’s initial drift velocity (shown in Figure 6 ). Similar behavior of the potential was reported experimentally by King et al. [29] for the potential of plasma plume. Similar results are also reported in Refs. [13, 14]. The enhanced perturbed potential

## 8. Discussions and summary

In conclusion, we can say that short-wavelength perturbations blow up exponentially much more quickly in RTI. The primary source by which this instability is triggered is the gravitational force acting on an inverted density gradient (e.g., a heavy fluid supported by a light fluid). Stable and steady flows may become unstable depending on the ranges of the flow parameters. The instability takes free energy from the mean flow or externally supplied heat and the amplitude of waves grows exponentially. The instabilities exist in all natural and artificial phenomena (in smoke from chimneys, in rivers, in flickering flames) and their effects result in turbulence or random waves. The presence of plasma density and magnetic field gradients is one of the main sources for plasma instabilities in Hall thrusters. It is found that perturbed potential increases with the higher value of electrons’ drift velocity, magnetic field, and ion temperature.

## Acknowledgments

The University Grants Commission (UGC), New Delhi, India is acknowledged for providing the startup Grant (No. F. 30-356/2017/BSR).