# Interaction of Dust-Acoustic Shock Waves in a Magnetized Dusty Plasma under the Influence of Polarization Force

Authors: N. S. Saini; Kuldeep Singh; Papihra Sethi
Journal: Laser and Particle Beams
 (2021)

---

## 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.


---

## Body

## 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 [3–6]. 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 [7–12] 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 Fe=qdE (qd=−Zde is dust charge, Zd 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 [17–25]. 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 [26–32]. 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 ofDA waves has been reported, yet the combined effects of polarization force and nonextensivity of ions on the interaction of DA 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]:(1)ni=ni01−q−1eψKBTiq+1/2q−1.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 beqd=−Zde, and ne0, ni0, and nd0 are the unperturbed electron, ion, and dust number densities, respectively. Zd is the dust grain charge. The charge neutrality condition at equilibrium yields ne0+Zdnd0=ni0 which can be expressed as μe+1=μi, where μe=ne0/Zdnd0=1/ρ−1 and μi=ni0/Zdnd0=ρ/ρ−1 and ρ=ni0/ne0. The electrons are considered Maxwellian with density ne=ne0expeψ/KBTe. After Taylor expansion of the RHS of equation (1) for eψ/KBTi≪1, we have(2)ni=ni01−h1eψKBTi+h2eψKBTi2,where h1=q+1/2 and h2=q+13−q/8. The expression for polarization force Fpol is given by [15](3)Fpol=qd22∇λDλD2.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λD for negative dust (i.e., neTi≪niTe) is λD≃KBTi/4πe2nih1 in the presence of nonextensive ions. The expression for nonextensively modified polarization force is obtained as(4)Fpol=−ZdeRh0nini01/2∇ψ,where R=Zde2/4KBTiλDi0 measures the polarization force effect, i.e., interaction of highly negatively charged dust with nonextensive ions with λDi0=KBTi/4πe2ni0h1 and h0=h1−2h2−h12/2eψ/KBTi. Thus, nonextensivity of ions (via q) notably alters the polarization force Fpol and also polarization parameter R. When q⟶1, 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 ofDA shocks waves in polarized dusty plasma is characterized by the following set of normalized fluid equations (continuity, momentum, and Poisson):(5)∂nd∂t+∇.ndud=0,(6)∂ud∂t+ud.∇ud=χp∇ψ+ηl∇2ud−Ωud×z^,(7)∇2ψ=nd+ne−ni,where ∇=x^∂/∂x+y^∂/∂y+z^∂/∂z and ud=x^udx+y^udy+z^udz. nd, ud, and ψ are the number density, velocity, and potential of dust fluid, σ=Ti/Te, and χp=1−Rh1−2h2+h12/2. Density njj=d,e,i is normalized by unperturbed density of dust nd0, dust velocity ud by ZdKBTi/md, ψ by KBTi/e, time t by ωpd−1, where ωpd=4πZde2nd0/md1/2, and x by λDd, where λDd=KBTi/4πnd0Zde21/2. ηl is the dust viscosity normalized by ωpdλDd2. The magnetic field is assumed to be along the z-axis (i.e., B=B0z^), and the dust cyclotron frequency Ω=ZdeB0/md is normalized by ωpd.


## 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:(8)Y=Y0+∑r=1∞ϵr+1Yr,Γ=∑r=1∞ϵr+2Γr,where Y=nd,udz,ψ, Y0=1,0,0, and Γ=udx,udy. The scaling variables x and t are stretched by the new coordinate system ξ, η, and τ in multiple scale variables:(9)ξ=ϵlxx+lyy+lzz−v0t+ϵ2M0η,τ+…,(10)η=ϵlxx+lyy+lzz+v0t+ϵ2N0ξ,τ+…,(11)τ=ϵ3t,where ξ and η denote the opposite side trajectories of DA shock waves and ϵ is a small parameter determining the strength of nonlinearity. The functions M0η,τ and N0ξ,τ are determined later. lx,ly, and lz represent the direction cosines along the x, y, and z directions, respectively, such that lx2+ly2+lz2=1. The value of dust viscosity is considered small and its value is set as ηl=ϵηl0, where ηl0 is O1. Making use of equations (8)–(11) into (5)–(7), at the lowest orders of ϵ, we obtain the following first-order equations:(12)v0−∂∂ξ+∂∂ηnd1+lz∂∂ξ+∂∂ηudz1=0,(13)−v0∂∂ξ+∂∂ηud1+1−Rh1∂∂ξ+∂∂ηψ1=0,(14)1−Rh1lx∂∂ξ+∂∂ηψ1−Ωudy1=0,(15)1−Rh1ly∂∂ξ+∂∂ηψ1+Ωudx1=0,(16)nd1=−a1ψ1,where a1=μeσ+μih1. On simplifying equations (12)–(16), the relations of different first-order physical quantities are obtained as follows:(17)ψ1=ψ1ξξ,τ+ψ1ηη,τ,(18)nd1=−1−Rh1lz2v02ψ1ξξ,τ+ψ1ηη,τ,(19)udx1=−1−Rh1lyωcdψ1ξξ,τ∂ξ+ψ1ηη,τ∂η,(20)udy1=1−Rh1lxωcdψ1ξξ,τ∂ξ+ψ1ηη,τ∂η,(21)udz1=−1−Rh1lzv0ψ1ξξ,τ−ψ1ηη,τ,(22)nd1=−a1ψ1ξξ,τ+ψ1ηη,τ.From equation (17), one predicts that there are two DA shock waves in which ψ1ξξ,τ is propagating towards the right and ψ1ηη,τ 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 v0=lz1−Rh1/a1.Finally, for the next higher orders ofϵ, we obtain(23)−4lz∂2udz3∂ξ∂η=∫∂ψ1ξξ,τ∂τ+Aψ1ξξ,τ∂ψ1ξξ,τ∂ξ+B∂3ψ1ξξ,τ∂ξ3−C∂2ψ1ξξ,τ∂ξ2dη+∫∂ψ1ηη,τ∂τ−Aψ1ηη,τ∂ψ1ηη,τ∂η−B∂3ψ1ηη,τ∂η3+C∂2ψ1ηη,τ∂η2dξ+∫∫2v0∂M0∂η+Dψ1ηη,τ∂2ψ1ξξ,τ∂ξ2dξdη+∫∫2v0∂N0∂ξ+Dψ1ξξ,τ∂2ψ1ηη,τ∂η2dξdη.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:(24)∂ϕξ∂τ+Aϕξ∂ϕξ∂ξ+B∂3ϕξ∂ξ3−C∂2ϕξ∂ξ2=0,(25)∂ϕη∂τ−Aϕη∂ϕη∂η−B∂3ϕη∂η3+C∂2ϕη∂η2=0.In the above equations, we have consideredψ1ξ=ϕξ and ψ1ξ=ϕη for simplicity. The nonlinear, dispersion, and dissipation coefficients, respectively, are given as(26)A=−3lz21−Rh12v0+v0ba1+v0R2h2+h12/221−Rh1,(27)B=1−lz21−Rh1v02a1Ω2+v02a1,(28)C=ηl02,where b=μeσ2/2−μih2. 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]:(29)ϕξ=ϕm1−tanhwξ−6C225Bτ2,(30)ϕη=ϕm1+tanhwη+6C225Bτ2.The electrostatic potential for weak head-on collision is given as(31)Φ=ϕm1−tanhwξ−6C225Bτ2+ϕm1+tanhwη+6C225Bτ2,where ϕm=3C2/25AB is the maximum amplitude and w=10B/C 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](32)2v0∂M0∂η=−Dϕη,(33)2v0∂N0∂ξ=−Dϕξ,where(34)D=v0ba1+1−Rh1lz22v0+v0R2h2+h12/221−Rh1.Equations (32) and (33) are solved using appropriate initial and boundary conditions to determine the phase shifts as(35)ΔM0η,τ≈−ϵ23ηl050v0CDAB,(36)ΔN0ξ,τ≈ϵ23ηl050v0CDAB.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.


## 3.1. Solution of KdVB Equation and Phase Shifts

The steady-state solutions of equations (24) and (25) are obtained as follows [27]:(29)ϕξ=ϕm1−tanhwξ−6C225Bτ2,(30)ϕη=ϕm1+tanhwη+6C225Bτ2.The electrostatic potential for weak head-on collision is given as(31)Φ=ϕm1−tanhwξ−6C225Bτ2+ϕm1+tanhwη+6C225Bτ2,where ϕm=3C2/25AB is the maximum amplitude and w=10B/C 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](32)2v0∂M0∂η=−Dϕη,(33)2v0∂N0∂ξ=−Dϕξ,where(34)D=v0ba1+1−Rh1lz22v0+v0R2h2+h12/221−Rh1.Equations (32) and (33) are solved using appropriate initial and boundary conditions to determine the phase shifts as(35)ΔM0η,τ≈−ϵ23ηl050v0CDAB,(36)ΔN0ξ,τ≈ϵ23ηl050v0CDAB.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], ni0=5×107cm−3,ne0=4×107cm−3,Zd=3×103,nd=107cm−3,Te=50eV,Ti=0.05eV,andR=0−0.14, 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 A (equation (26)) is negative (i.e., A<0) and the dispersion coefficient B (equation (27)) is positive (i.e., B>0) for the chosen set of parameters. It implies that AB<0; 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 (viaq) on the Polarization Force

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


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

Figure2(a) depicts the variation of the amplitude of negative polarity DA shock waves in R−q plane. It is found that the amplitude of negative polarity DA shock waves is increased with rise in both polarization force parameter (via R) and nonextensive parameter (via q). 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 lz). 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.Figure 2
The variation of rarefactive DA shocks profile with respect toξ (a) for different values of polarization parameter (R) and nonextensivity parameter (q) and (b) for different values of viscosity ηl0 and magnetic field strength whereas other parameters ρ=1.3, lz=0.8, σ=0.25, Ω=0.8, and ηl=0.35 are fixed.
(a)(b)


### 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 Figure3. 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 (ηl0). 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 q for different values of polarization parameter R in the presence of magnetic field. It is seen that magnitude of phase shift escalates with rise in nonextensive parameter q and suppresses with increase in polarization force (via R). 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 R), nonextensivity parameter (via q), strength of magnetic field (via Ω), and other plasma parameters.Figure 3
The collision of two DA shock waves at different time intervals forq=0.33, R = 0.1, ρ=1.3, lz=0.8, σ=0.25, Ω=0.8, and ηl=0.2. (a) τ=50, (b) τ=40, (c) τ=30, (d) τ=20, (e) τ=0, (f) τ=−50, (g) τ=−40, (h) τ=−30, and (i) τ=−20.Figure 4
Collision process of two DA shock waves for (a)ηl0=0.1 and (b) ηl0=0.2 and ρ=1.3, lz=0.8, σ=0.25, Ω=0.8, R=0.11, and q=0.33 are fixed.
(a)(b)Figure 5
Collision process of two DA shock waves for (a)Ω=0.2 and (b) Ω=0.3 and ρ=1.3, lz=0.8, σ=0.25, R=0.11, ηl=0.2, and q=0.3 are fixed.
(a)(b)Figure 6
The variation of phase shift of two DA shock waves vs.q for different values of R (a) with magnetic field and (b) without magnetic field and ρ=1.3, lz=0.8, σ=0.25, and ηl=0.2 are fixed.
(a)(b)


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

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


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

Figure2(a) depicts the variation of the amplitude of negative polarity DA shock waves in R−q plane. It is found that the amplitude of negative polarity DA shock waves is increased with rise in both polarization force parameter (via R) and nonextensive parameter (via q). 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 lz). 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.Figure 2
The variation of rarefactive DA shocks profile with respect toξ (a) for different values of polarization parameter (R) and nonextensivity parameter (q) and (b) for different values of viscosity ηl0 and magnetic field strength whereas other parameters ρ=1.3, lz=0.8, σ=0.25, Ω=0.8, and ηl=0.35 are fixed.
(a)(b)


## 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 Figure3. 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 (ηl0). 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 q for different values of polarization parameter R in the presence of magnetic field. It is seen that magnitude of phase shift escalates with rise in nonextensive parameter q and suppresses with increase in polarization force (via R). 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 R), nonextensivity parameter (via q), strength of magnetic field (via Ω), and other plasma parameters.Figure 3
The collision of two DA shock waves at different time intervals forq=0.33, R = 0.1, ρ=1.3, lz=0.8, σ=0.25, Ω=0.8, and ηl=0.2. (a) τ=50, (b) τ=40, (c) τ=30, (d) τ=20, (e) τ=0, (f) τ=−50, (g) τ=−40, (h) τ=−30, and (i) τ=−20.Figure 4
Collision process of two DA shock waves for (a)ηl0=0.1 and (b) ηl0=0.2 and ρ=1.3, lz=0.8, σ=0.25, Ω=0.8, R=0.11, and q=0.33 are fixed.
(a)(b)Figure 5
Collision process of two DA shock waves for (a)Ω=0.2 and (b) Ω=0.3 and ρ=1.3, lz=0.8, σ=0.25, R=0.11, ηl=0.2, and q=0.3 are fixed.
(a)(b)Figure 6
The variation of phase shift of two DA shock waves vs.q for different values of R (a) with magnetic field and (b) without magnetic field and ρ=1.3, lz=0.8, σ=0.25, and ηl=0.2 are fixed.
(a)(b)


## 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, andq-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 q-nonextensive ions yield to escalate the polarization force parameter (R). 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 q), modified polarization force parameter (via R), viscosity (via ηl0), obliqueness (via lz), 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.



---

*Source: 6679085-2021-01-21.xml*