Investigation of the Rayleigh-Taylor instability in charged fluids

Weiping Zhang,Zhongzheng Li and Wenshan Duan

1 College of Physics and Electronics Engineering,Northwest Normal University,Lanzhou 730070,China

2 Department of Physics,Gansu Normal University For Nationalities,Hezuo 747000,China

Abstract The present study shows that the Rayleigh–Taylor (RT) instability and its growth rate are strongly dependent on the charge-mass ratio of charged particles in a charged fluid.A higher charge-mass ratio of the charged fluid appears to result in a stronger effect of the magnetic field to suppress the RT instability.We study the RT instabilities for both dusty plasma(small chargemass ratio of charged particles) and ion-electron plasma (large charge-mass ratio of charged particles).It is found that the impact of the external magnetic field to suppress the RT instability for ion-electron plasma is much greater than that for dusty plasma.It is also shown that,for a dusty plasma,in addition to region parameters such as the external magnetic field,region length,its gradient,as well as dust particle parameters such as number density,mass,and charge of dust particles,the growth rate of the RT instability in a dusty plasma also depends on parameters of both electrons and ions such as the number densities and temperatures of both electrons and ions.

Keywords: Rayleigh–Taylor (RT) instability,charged fluid,magnetized plasmas,dusty plasma

1.Introduction

A dusty plasma is a unique type of plasma that consists of three distinct components: electrons,ions,and charged dust particles.These dust particles are typically found within a size range of nanometers to several hundreds of micrometers[1,2].Dusty plasma,also known as complex plasma,has been the subject of extensive research over the past decades[3–5].Dusty plasma is present in a wide range of cosmic environments [3,6–8],including the solar system [6],planetary rings [3,7],interstellar clouds,circumstellar clouds and various space environments [8].Understanding the properties and behavior of dusty plasmas is crucial for our understanding of many natural phenomena in space.

There are significant differences between dusty plasma and ion-electron plasma due to the relatively large mass of dust particles [9–11] compared to electrons and ions.One of the obvious differences of the dusty plasma is that the influence of gravity on a dust particle cannot be ignored,particularly on Earth.What is more important is that dusty plasma exhibits complex behavior,including dust acoustic waves [4],dust ion acoustic waves [12],and others[5,13–17].These unique properties make dusty plasma an intriguing subject for research in astrophysics,space science,and plasma physics.

Over the past few years,researchers have studied a binary dusty plasma that contains two types of microparticles of different sizes[18,19].This binary dusty plasma can either be mixed [20–22] or form a phase-separated system [23–26].The separation may be caused by spinodal decomposition[25]or an imbalance of external forces[26].In the latter case,an interface emerges between the separated phases.

As is widely recognized,the Rayleigh–Taylor (RT)instability is a hydrodynamic instability that occurs at the interface between two fluids of different densities when the lighter fluid is accelerated to move into the denser fluid due to an external force [27–29].This instability is a result of the growth of small perturbations or waves that occur at the interface,causing the fluids to mix and form complex patterns.The RT instability has significant applications in various scientific and engineering fields [30–36],including astrophysics [37],plasma physics [38–41],and others[42–46].Understanding the RT instability is essential for predicting and controlling the behavior of fluids in these applications.Furthermore,there have been many studies on the RT instability in dusty plasmas [27,28,47–49].For instance,the stability of stratified dust clouds suspended in anode plasma was analyzed,and the RT instability was detected [28].Additionally,the RT instability finds broad applications in inertial confinement fusion[50,51],laboratory plasma [52,53],as well as astrophysical settings like crab nebulae [54] and supernova remnants [55].

Research on RT instability in dusty plasmas has been conducted [56],and while recent studies have investigated this phenomenon,the effects of external magnetic and electric fields on RT instability are still unclear.The present paper aims to study the influence of the magnetic field on RT instability and has found that the charge-mass ratio of the charged particles in the plasma plays a significant role in RT instability and its growth rate.It is worth noting that the larger the charge-mass ratio of the charged particles,the greater the effect of the magnetic field in suppressing the RT instability.Thus,the external magnetic field’s ability to suppress the RT instability in ion-electron plasma is much stronger than in dusty plasma.In addition,it has been discovered that in a dusty plasma,the RT instability and its growth rate are influenced not only by region parameters(such as the external magnetic field,the length and gradient of the region)and dust particle parameters (such as dust particle number density,mass,and charge),but also by electron and ion parameters,including their number densities and temperatures.

2.Theoretical model

2.1.RT instability in a dusty plasma

We will study two typical charged fluids.One is the dusty plasma,the other is ion-electron plasma.Following,we first study RT instability in a dusty plasma,which is composed of dust particles,free electrons,and free ions.

Dusty plasma has been found in space plasma,as well as in the laboratory on earth [57–61].However,it is difficult to form three-dimensional dusty plasma lattice in the laboratory on Earth because of the asymmetric gravity acting on dust particles.Recently,three-dimensional dusty plasma lattices have been achieved under micro-gravity conditions[24]in the Russian–German PK-3 Plus Laboratory aboard the International Space Station(ISS).Furthermore,an interface emerges between two separated phases [18,19].

Fig. 1.The RT model of two regions.The region I is from z=β to z=α,while the region II is from z=α to z=α+Δ.

To describe this type of dusty plasma,we utilize the generalized hydrodynamic model,which takes into account the forces produced by the pressure gradient as well as those produced by the electric and magnetic fields acting on the dust particles.The resulting momentum equation for a dusty plasma can be expressed as follows=(0,k y,kz).Then the hydrodynamic equations to describe

We consider a model shown in figure 1,where α is the interface,the region I is from z=β to z=α,while the region II is from z=α to z=α+Δ.The gravity force is along the negative direction of the z-axis.

For simplicity and generality,we study the waves propagating only in the yoz plane,i.e.the wave number→this kind of dusty plasma are as follows

where nd0is the number density of the dust particle at equilibrium state,nd1is the perturbed number density of the dust particle,vd1and wd1are the perturbed velocities of the dust fluid in the y and z directions,respectively,φ1is the perturbed electric potential,B0is a constant.

In order to understand the RT instability of the system,we assume that the system is at an equilibrium state initially.Then by linearizing equations (2)–(5),neglecting any nonlinear terms and substituting equations (7)–(11) into equations (2)–(5),we have

Simplifying equations (12)–(15),we have

We now consider a dusty plasma in which there are two different regions,Region I: β

Rewriting equations (18) and (19),we have

Fig. 2.(a)Dependence of the growth rate of the RT instability on the wave number k for different parameters of Zd=2700,3600,4500,respectively.(b)Dependence of the growth rate of the RT instability on the wave number k for different parameters of B0x=0.4×10-4 T,0.5×10-4 T,0.6×10-4 T,respectively.(c) Dependence of the growth rate of the RT instability on the wave number k for different parameters of md=3.6×10-15 kg,8.5×10-15 kg,1.3×10-14 kg,respectively.(d)Dependence of the growth rate of the RT instability on the wave number k for different parameters of L=2.0 m,2.5 m,3.0 m,respectively.(e)Dependence of the growth rate of the RT instability on the wave number k for different parameters of λD=6×10-5 m,8×10-5 m,1×10-4 m,respectively.

The dispersion relation is

The growth rate of the RT instability can be expressed as follows

Figure 2 shows the dependence of the growth rate of the RT instability on the wave number k for different Zdin figure 2(a),different B0xin figure 2 (b),different mdin figure 2 (c),different values of L in figure 2(d),different λDin figure 2(e).The plasma parameters are chosen from the previous work[18]md=1.3×10-14kg,B0x=0.4×10-4T,L=2.0 m,λD=1.0×10-4m in figure 2(a).md=1.3×10-14kg,Zd=4500,λD=1.0×10-4m in figure 2 (b).Zd=4500,B0x=0.4×10-4T,L=2.0 m,λD=1.0×10-4m in figure 2 (c).md=1.3×10-14kg,Zd=4500,B0x=0.4×10-4T,λD=1.0×10-4m in figure 2 (d).md=1.3×10-14kg,Zd=4500,B0x=0.4×10-4T,L=2.0 m in figure 2 (e).

It seems that the number density and the temperature of both electrons and ions can affect the RT instability.It indicates that the larger the number densities of both electrons and ions,the smaller the growth rate,while the higher the temperatures of both electrons and ions,the larger the growth rate.It is found that the RT instability and its growth rate not only depend on the region parameters (the external magnetic field,the region length and its gradient)and the dust particles parameters (the number density,the mass and the charge of the dust particle),but also on both electron and ion characters(the number densities,the temperatures of both electrons and ions).

2.2.RT instability in an ion-electron plasma

Using the incompressibility condition,

We assume that

where ni0is the number density of ions at equilibrium state,ni1is the perturbed number density of ions,vi1and wi1are the perturbed velocities of the ion in the y and z directions respectively,φ1is the perturbed potential,B0is a constant.

In order to understand the RT instability of the system,we assume that the system is at an equilibrium state initially.Then by linearizing equations (24)–(27),neglecting any nonlinear terms and substituting equations (28)–(33) into equations (24)–(27),we have

Simplifying equations (34)–(37),we have

Rewriting equations (40) and (41),we have

Moreover,the eigenvalues corresponding to equation (43) are as follows

The dispersion relation is

When B0x=0,the dispersion relation is

The growth rate of the RT instability can be expressed as follows

It seems that the dispersion relation equation (46) is consistent with that of the general fluid [67] and therefore there is RT instability.

2.3.Effects of the charge to mass ratio of charged fluids on the growth rate of the RT instability

The present results show that an external magnetic field is more effective on suppressing the RT instability in ion-electron plasma than that in dusty plasma.The explanation for this difference can be inferred by examining both equations (23)and(45).We can define a refrain factor to RT instability which is given bywhere Q is the electric charge of the charged particle and M is the mass of the charged particle.Put simply,the degree to which the RT instability is suppressed depends on the product of the external magnetic field and the charge-mass ratio of the charged particles.Due to the much larger charge-mass ratio of ions compared to dust particles,the effect of the external magnetic field in suppressing the RT instability is much more pronounced in ion-electron plasma than that in dusty plasma.In other words,the greater the charge-mass ratio of the charged particle,the stronger the effect of the magnetic field in suppressing the RT instability.

3.Summary and conclusion

The present paper provides a theoretical investigation of the RT instability in different kinds of charged fluid.One is a dusty plasma,while the other is ion-electron plasma.It is found that the charge-mass ratio of charged particles in a charged fluid plays a crucial role in the RT instability and its growth rate.It appears that a larger charge-mass ratio leads to a stronger effect of the external magnetic field in restraining the RT instability.Consequently,the effect of the magnetic field in refraining the RT instability is much more significant in ion-electron plasma compared to dusty plasma due to the ions’ larger charge-mass ratio.

It is also found that,for a dusty plasma,the growth rate of the RT instability increases as the perturbed wave number increases.Moreover,the RT instability and its growth rate depend not only on the external parameters of the region(the external magnetic field,the length of the region,and its gradient),the parameters of the dust particles (the number density,the mass,and the charge),but also on the characteristics of both the electrons and ions such as their number densities and their temperatures.

Currently,our research may have potential applications.For example,it is almost impossible to refrain RT instability by an external magnetic field for a dusty plasma.However,it is easy to refrain the RT instability by an external magnetic field for an ion-electron plasma.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Nos.12275223,11965019) and by the Foundation of Gansu Educational Committee (No.2022QB-178).