Abstract

The interaction of dust-acoustic (DA) shock waves in a magnetized dusty plasma under the influence of nonextensively modified polarization force is investigated. The plasma model consists of negatively charged dust, Maxwellian electrons, nonextensive ions, and polarization force. In this investigation, we have derived the expression of polarization force in the presence of nonextensive ions and illustrated the head-on collision between two DA shock waves. The extended Poincaré–Lighthill–Kuo (PLK) method is employed to obtain the two-sided Korteweg–de Vries–Burgers (KdVB) equations and phase shifts of two shock waves. The trajectories and phase shifts of negative potential dust-acoustic shock waves after collision are examined. The combined effects of various physical parameters such as polarization force, nonextensivity of ions, viscosity of dust, and magnetic field strength on the phase shifts of DA shock waves have been studied. The present investigation might be useful to study the process of collision of nonlinear structures in space dusty plasma such as planetary rings where non-Maxwellian particles such as nonextensive ions, negatively charged dust, and electrons are present.

1. Introduction

Dusty plasma has fascinated the attention of many researchers for the study of different kinds of linear as well as nonlinear structures over the past many years. The presence of dust in an electron-ion plasma modifies the behaviour of nonlinear structures and generates new modes such as dust ion-acoustic waves [1] and dust-acoustic waves [2], which have been frequently studied theoretically and experimentally [36]. Dust-acoustic waves are low-frequency waves in which the inertia is provided by the dust mass and the restoring force is provided by the inertialess electrons and ions. Numerous investigations have reported the prevalent role of dust encountered in noctilucent clouds, planetary rings, comets, and plasma crystals [712] as well as in semiconductor manufacturing and ion implantation [13, 14]. Dust particulates in a dusty plasma are affected by various forces due to the large amount of dust charge on them, among which electrostatic force ( is dust charge, is the number of electrons residing on dust grain surface, and E is the electric field) has a significant influence on dust. As the bulk plasma is quasineutral, the electric field there is usually much weaker than the electric field in the sheath region. Hamaguchi Farouki [15, 16] reported that polarization force acts on the dust in the background of nonuniform plasma as a result of deformation of Debye sheath around the particulates. Spatial nonuniform dust density means that the dust charge is nonuniformly distributed in plasmas. The polarization force acts in the direction in which Debye length decreases and also acts in the direction opposite to the electrical force independent of the sign of dust charge [17]. Thus, the polarization force has a significant influence on the propagation characteristics of low-frequency linear DA waves. Various researchers have reported the effect of polarization force on different nonlinear structures (solitons and shocks) in the framework of Maxwellian and non-Maxwellian distribution in various plasma environments [1725]. Sethi et al. [23] studied DA solitary and shock waves under the influence of superthermally modified polarization force in an unmagnetized dusty plasma comprising of negatively charged dust grains, superthermally distributed ions, and Maxwellian electrons. The superthermally modified polarization force was found to be greater for more superthermal particles. Singh et al. [24] illustrated that the superthermality of ions has a significant influence on polarization force which further altered the dynamics of DA cnoidal waves in dusty plasma. In the recent past, Singh and Saini [25] examined the impact of superthermal polarization force on the modulational instability of DA waves and transition of rogue wave triplets to super rogue waves in space dusty plasmas.

Wave-wave interaction is a very fascinating phenomenon apart from the wave propagation in plasmas that can modify the properties of waves. This unique phenomenon of wave-wave interaction can be described by the use of the extended PLK method and has been investigated by numerous authors [2632]. The interaction between waves can take place as overtaking collision or as head-on collision. It is well acknowledged that solitons preserve their asymptotic form while undergoing collision. During the before and after collisions, the solitary waves retain their shape, and also they have a novel impact on phase shifts. In this investigation, we have focused our study on the head-on collision of two DA shock waves under the effect of polarization force and their trajectories after the collision. The investigations of the head-on collision of nonlinear shock waves under the influence of polarization force have been reported by numerous researchers in the framework of Maxwellian and non-Maxwellian distributions in different plasma systems [26, 27, 31]. Head-on collision of modulated dust-acoustic waves in strongly coupled dusty plasma is investigated by El-Labany et al. [26]. They considered negatively charged dust grains and Maxwellian electrons and ions and studied modulational instability and effect of polarization force on the phase shifts of dust-acoustic solitary waves. Head-on collision of dust-acoustic shock waves in strongly coupled dusty plasmas was investigated by El-Shamy and Al-Asbali [27]. In the recent past, Singh et al. [31] investigated the head-on collision in a dusty plasma among DA multisolitons with ions following the hybrid distribution under the effect of the polarization force. It was found that rarefactive DA multisolitons are formed and phase shifts are strongly influenced by the polarization force.

Although prolific literature on the interaction of waves has been reported, yet the combined effects of polarization force and nonextensivity of ions on the interaction of shock waves have not been investigated by researchers in dusty plasmas with application to planetary rings and comet tails. The manuscript is structured as follows. Section 2 presents the fluid model for a given plasma system. Derivations of two different types of KdV–Burgers (KdVB) equations and their corresponding phase shifts are displayed in Section 3. Numerical analysis and discussion are presented in Section 4. The last Section 5 is devoted to the conclusions.

2. The Fluid Equations

In statistical mechanics, there exist generally two types of systems such as extensive and nonextensive systems, depending on the range of interparticle forces. The extensive systems are considered for short-range interparticle forces and need Boltzmann–Gibbs (B-G) statistics. On the other hand, nonextensive systems hold for long-range interparticle forces, need generalized B-G statistics, and also include dissipative systems having nonvanishing thermodynamic currents. Tsallis [33] introduced q-nonextensive distribution which is used to model astrophysical regions (i.e., gravitational systems and stellar polytropes) [34, 35] where Maxwellian distribution might become inadequate. Tsallis q-entropy and the probabilistically independent postulate are known to be the basis of the nonextensive statistics. The principle of maximum q-entropy gives the power-law q-distribution which is called q-exponentials in nonextensive statistics. The nonextensivity is also known as pseudoadditivity or nonadditivity which gives both the q-entropy and the energy (for more details, see [36]). The parameter q measures the degree of nonextensivity of the system. The q-nonextensive number density of ions is given as follows [37]:

We consider a magnetized dusty plasma comprising of inertial negatively charged dust, nonextensive ions, and Maxwellian electrons to study the head-on collision and phase shifts of two DA shock waves under the influence of nonextensively modified polarization force. Consider charge on negatively charged dust grain be , and , , and are the unperturbed electron, ion, and dust number densities, respectively. is the dust grain charge. The charge neutrality condition at equilibrium yields which can be expressed as , where and and . The electrons are considered Maxwellian with density . After Taylor expansion of the RHS of equation (1) for , we havewhere and . The expression for polarization force is given by [15]

We have now derived the polarization force (i.e., deformation of Debye sphere around the dust grains) acting on dust grains modified by nonextensive ions in a nonuniform dusty plasma comprising of nonextensive ions and Maxwellian electrons where the expression for modified Debye length for negative dust (i.e., ) is in the presence of nonextensive ions. The expression for nonextensively modified polarization force is obtained aswhere measures the polarization force effect, i.e., interaction of highly negatively charged dust with nonextensive ions with and . Thus, nonextensivity of ions (via ) notably alters the polarization force and also polarization parameter R. When , the polarization force expression equation (4) agrees well with the expression for the case of Maxwellian ions and electrons derived by Asaduzzaman et al. [18].

The dynamics of shocks waves in polarized dusty plasma is characterized by the following set of normalized fluid equations (continuity, momentum, and Poisson):where and . , , and are the number density, velocity, and potential of dust fluid, , and . Density is normalized by unperturbed density of dust , dust velocity by , by , time t by , where , and by , where . is the dust viscosity normalized by . The magnetic field is assumed to be along the z-axis (i.e., ), and the dust cyclotron frequency is normalized by .

3. Derivation of KdVB Equations and Phase Shifts

We employ an extended PLK perturbation method to investigate the collision between two DA shock waves in a magnetized dusty plasma. The perturbed quantities are described as follows:where , , and . The scaling variables and are stretched by the new coordinate system , , and in multiple scale variables:where and denote the opposite side trajectories of DA shock waves and is a small parameter determining the strength of nonlinearity. The functions and are determined later. , and represent the direction cosines along the , , and directions, respectively, such that . The value of dust viscosity is considered small and its value is set as , where is . Making use of equations (8)–(11) into (5)–(7), at the lowest orders of , we obtain the following first-order equations:where . On simplifying equations (12)–(16), the relations of different first-order physical quantities are obtained as follows:

From equation (17), one predicts that there are two DA shock waves in which is propagating towards the right and is propagating towards the left. On solving equations (17)–(22), the expression for phase velocity in a magnetized dusty plasma under the effect of nonextensively modified polarization force is obtained as .

Finally, for the next higher orders of , we obtain

From equation (23), it is noticed that the first (second) term of the right side is proportional to (), and thus, the integrating function is independent of (). These two secular terms in the above equation must be removed to avoid spurious resonance. Hence, we obtain the following nonlinear equations:

In the above equations, we have considered and for simplicity. The nonlinear, dispersion, and dissipation coefficients, respectively, are given aswhere . Equations (24) and (25) are the KdV–Burgers equations of two shock waves travelling towards each other in the reference frame of and .

3.1. Solution of KdVB Equation and Phase Shifts

The steady-state solutions of equations (24) and (25) are obtained as follows [27]:

The electrostatic potential for weak head-on collision is given aswhere is the maximum amplitude and is the width of DA shock waves.

In the next higher order, the third and fourth terms on the right-hand side of equation (23) become secular and we obtain [27]where

Equations (32) and (33) are solved using appropriate initial and boundary conditions to determine the phase shifts as

These phase shifts are functions of different coefficients that are dependent on various physical parameters. Therefore, any variation in physical parameters makes a significant change in the phase shifts.

4. Numerical Analysis and Discussion

In this investigation, the head-on collision of DA shock waves in a magnetized dusty plasma consisting of nonextensive ions and Maxwellian electrons using the extended PLK perturbation method is presented. The numerical values of the typical physical parameters of planetary rings [20, 24], , are used for analysis. The shock waves are formed as a result of balance of nonlinearity, dispersion, and dissipation. The KdVB equations have been derived in equations (24) and (25). Numerically, it is seen that the nonlinear coefficient (equation (26)) is negative (i.e., ) and the dispersion coefficient (equation (27)) is positive (i.e., ) for the chosen set of parameters. It implies that ; as a result, negative potential shock structures are formed. Since the various physical parameters such as nonextensivity of ions, polarization force, obliqueness, magnetic field strength, and other plasma parameters have a profound influence on the characteristics of DA shock waves, it is of paramount importance to carry out numerical analysis.

4.1. Influence of Nonextensivity of Ions (via ) on the Polarization Force

Figure 1 depicts the variation of the polarization force parameter with nonextensive parameter . Since the parameter is a function of nonextensivity of ions (via q), any variation in will change the value of . It is observed that as the nonextensive parameter is increased, the polarization force parameter increases. Thus, nonextensivity of ions has a profound influence on the polarization force parameter (i.e., polarization force).

4.2. The Variation of Profile of Negative Polarity DA Shock Waves

Figure 2(a) depicts the variation of the amplitude of negative polarity DA shock waves in plane. It is found that the amplitude of negative polarity DA shock waves is increased with rise in both polarization force parameter (via ) and nonextensive parameter (via ). Figure 2(b) shows that the amplitude of negative polarity shock waves is also enhanced with an increase in both magnetic field strength (via ) and obliqueness (via ). It is remarked that the amplitude of DA shock waves is modified with the variation in plasma parameters due to the change in nonlinear as well as dispersive effects under the influence of polarization force in the presence of dissipation effect.

4.3. Evolution and Phase Shift of Two Negative Polarity DA Shock Waves

The collision process of DA shock waves in a magnetized dusty plasma for different time intervals is depicted in Figure 3. The mesial line indicated in the figure separates the before and after collision process. The upper red curves illustrate the profiles of DA shocks before collision and lower blue curves for after collision. This whole process illustrates that two opposite directional propagating negative polarity DA shock waves come closer, interact, and then depart (for more clarity, see animation of collision in the supplementary section (available here)). Figures 4(a) and 4(b) depict the variation of phase shift occurred during the collision process between two negative polarity DA shock waves for different values of dust fluid viscosity (). It is recognised that the magnitude of phase shift escalates with rise in viscosity (clearly seen in color bar). Figures 5(a) and 5(b) depict the variation of phase shift occurred during the collision process between two negative polarity DA shock waves for different values of magnetic field strength (via ). It is found that the magnitude of phase shift enhances as the value of increases (see side color bar). Figure 6(a) illustrates the variation of phase shift with nonextensive parameter for different values of polarization parameter in the presence of magnetic field. It is seen that magnitude of phase shift escalates with rise in nonextensive parameter and suppresses with increase in polarization force (via ). Figure 6(b) shows the same variation in the absence of magnetic field. It is found that magnitude of phase shift is more (almost four times greater) than the case in the absence of magnetic field. It is remarked that magnetic field strength has a significant influence on the profile of DA shocks. It is concluded that the phase shift of DA shock waves is significantly affected under the influence of polarization force (via ), nonextensivity parameter (via ), strength of magnetic field (via ), and other plasma parameters.

5. Conclusions

In the present investigation, we have studied the effects of nonextensively modified polarization force on the head-on collision between two DA shock waves in a magnetized dusty plasma composed of negatively charged dust, Maxwellian electrons, and -nonextensive ions. The extended PLK method is employed to derive the two opposite directional KdVB equations and the expressions for phase shifts after collision. The -nonextensive ions yield to escalate the polarization force parameter (). Only negative polarity DA shock waves are formed under the influence of polarization force. The amplitude of DA shock waves and their phase shifts are significantly modified by the physical parameters such as nonextensivity of ions (via ), modified polarization force parameter (via ), viscosity (via ), obliqueness (via ), and magnetic field strength (via ). The time evolution of rarefactive DA shock waves under the influence of polarization force has been illustrated. The physical parameters have an emphatic effect on the characteristics of DA shock waves and their phase shifts occurred due to head-on collision among DA shock waves. It is important to conclude that the results of our present investigation might be useful for understanding the nonlinear features of disturbances in space/astrophysical dusty plasma systems (e.g., Saturn’s rings, interstellar clouds, and cometary environment) where negatively charged dust and nonextensive ions are present.

Data Availability

The derived data generated in this research will be shared on reasonable request to the corresponding author.

Disclosure

This paper was presented as an invited talk at the 12th International Conference on Plasma Science and Applications, held from 11 to 14 November 2019, Lucknow, India.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflicts of interest.

Acknowledgments

This work was supported by DRS-II (SAP) no. F 530/17/DRS-II/2015(SAP-I) UGC, New Delhi.

Supplementary Materials

Shocks-collision video. (Supplementary Materials)