# Magnetogasdynamic Shock Waves in a Rotating Gas with Exponentially Varying Density

Authors: J. P. Vishwakarma; G. Nath
Journal: ISRN Mathematical Physics
 (2012)

---

## Abstract

Nonsimilar solutions are obtained for one-dimensional adiabatic flow behind a magnetogasdynamic cylindrical shock wave propagating in a rotating or nonrotating
perfect gas in presence of a constant azimuthal magnetic field. The density of the gas is assumed to be varying and obeying an exponential law. In order to obtain the solutions, the angular velocity of the ambient medium is assumed to be decreasing exponentially as the distance from the axis increases. The shock wave moves with variable velocity and the total energy of the wave is nonconstant. The effects of variation of Alfven-Mach number and time are obtained. Also, a comparison between the solutions in the cases of rotating and non-rotating media with or without magnetic field is made.


---

## Body

## 1. Introduction

Hayes [1], Laumbach and Probstein [2], DebRay [3], Verma and Vishwakarma [4, 5], Vishwakarma [6], Vishwakarma and Nath [7], and Vishwakarma et al. [8] have discussed the propagation of shock waves in a medium where density varies exponentially and obtained similarity or nonsimilarity solutions. These authors have not taken into account the effects of rotation of the medium.The formation of self-similar problems and examples describing the adiabatic motion of nonrotating gas models of stars is considered by Sedov [9], Zel'dovich and Raizer [10], Lee and Chen [11], and Summers [12]. Rotation of the stars significantly affects the process taking place in their outer layers. Therefore questions connected with the explosions in rotating gas atmospheres are of definite astrophysical interest. Chaturani [13] studied the propagation of cylindrical shock waves through a gas having solid body rotation and obtained the solution by similarity method adopted by Sakurai [14]. Nath et al. [15] obtained the similarity solutions for the flow behind spherical shock waves propagating in a nonuniform and rotating interplanetary atmosphere with increasing energy. Vishwakarma et al. [16] obtained the similarity solution for the magnetogasdynamics cylindrical shock waves propagation in a rotating medium in the presence of a constant azimuthal magnetic field by assuming that the density of the medium ahead of the shock wave is constant. Levin and Skopina [17] studied the propagation of detonation wave in rotational gas flows by taking into account the variable azimuthal and axial fluid velocities, and the components of the vorticity vector. Nath [18] obtained the similarity solution for the magnetogasdynamic shock wave generated by a moving piston in a rotational axisymmetric isothermal flow of perfect gas with variable density according to power law. Also, Nath [19] obtained the nonsimilarity solutions for the propagation of a strong cylindrical shock waves in a rotational axisymmetric dusty gas with exponentially varying density.In the present work, we investigated the effects of the presence of an ambient azimuthal magnetic field and the rotation of the ambient medium on the flow field behind a magnetogasdynamic cylindrical shock wave. Nonsimilarity solutions for the flow field behind the shock wave are obtained. The density in the medium ahead of the shock is assumed to obey an exponential law. The angular velocity of rotation of the ambient medium is also assumed to be obeying exponential law and to be decreasing as the distance from the axis increases. It is expected that such an angular velocity may occur in the atmospheres of rotating planets and stars. The medium is assumed to be a perfect gas and the initial magnetic field to be constant.Variation of the flow variables behind the shock for different values of Alfven-Mach number and time is obtained. Also, a comparison between the solutions in the cases of rotating and nonrotating media is made for both the magnetic and nonmagnetic cases.


## 2. Fundamental Equations and Boundary Conditions

The fundamental equations governing the unsteady adiabatic cylindrically symmetric flow of an electrically conducting gas, which is rotating about the axis of symmetry in the presence of an azimuthal magnetic field may, in Eulerian coordinates, can be expressed as [13, 16, 20]
(2.1)∂ρ∂t+u∂ρ∂r+ρ∂u∂r+ρur=0,(2.2)∂u∂t+u∂u∂r+1ρ(∂p∂r+μh∂h∂r+μh2r)-v2r=0,(2.3)∂h∂t+u∂h∂r+h∂u∂r=0,(2.4)∂v∂t+u∂v∂r+uvr=0,(2.5)∂U∂t+u∂U∂r-pρ2(∂ρ∂t+u∂ρ∂r)=0,
where r and t are independent space and time coordinates, ρ the density, p the pressure, h the azimuthal magnetic field, U the internal energy per unit mass, u and v are the radial and azimuthal components of the fluid velocity, and μ is the magnetic permeability. The electrical conductivity of the gas is assumed to be infinite and the effects of viscosity and heat conduction are not considered.Also,(2.6)v=Ar,
where "A"  is the angular velocity of the medium at radial distance r from the axis of symmetry. The magnetic field equation (2.3) contains all the relevant information needed from Maxwell's equations and Ohm's law; the diffusion term is omitted from it by virtue of the assumed infinite electrical conductivity. The assumption of infinite electrical conductivity of the gas is physically realistic in the case of astrophysical phenomena where the magnetic Reynolds number is very high (or infinite) due to astrophysical scale. The magnetic Reynolds number is a dimensionless parameter defined by Rm=U¯L¯/ηm, where U¯,L¯ are the characteristic velocity and characteristic length of the flow field, and ηm=1/μσ is the magnetic diffusivity (magnetic viscosity), σ being the electrical conductivity of the medium  [21, page 169]. In this case, the Reynolds number Re of the flow is also very high, as Re is defined by Re=U¯L¯/ν, where ν is the kinematic coefficient of viscosity. It is well known that the effects of viscosity and heat conduction are negligibly small for the high Reynolds number flow except in the boundary layer region near the solid boundary or in any other region of large variations in velocity and temperature such as the inside of a shock wave [22, page 139]. Therefore omission of the effects of viscosity and heat conduction on the flow field may be justified in the present study.The above system of equations should be supplemented with an equation of state. A perfect gas behaviour of the medium is assumed, so that(2.7)p=ΓρT,U=pρ(γ-1),
where Γ is the gas constant and γ is the ratio of specific heats at constant pressure and volume.We assume that a cylindrical shock is propagating outwordally in a perfect gas with infinite electrical conductivity and variable density in presence of a constant azimuthal magnetic field.The ambient density of the medium is assumed to obey the exponential law, namely,(2.8)ρa=ρ0eδR,
where R is the shock radius and ρ0, δ are suitable constants.In order to obtain the solution, it is assumed that the ambient angular velocity of the mediumAa varies as (2.9)Aa=A0eαR,
where A0 and α are constants.The jump conditions at the shock wave are given by the principles of conservation of mass, momentum, magnetic field, and energy across the shock [20, 23], namely, (2.10)ρaV=ρn(V-un),haV=hn(V-un),pa+12μha2+ρaV2=pn+12μhn2+ρn(V-un)2,Ua+paρa+μha2ρa+12V2=Un+pnρn+μhn2ρn+12(V-un)2,va=vn,
where the subscripts “a”   and  “n”   denote the conditions immediately ahead and behind of the shock front, and V(=dR/dt) denotes the velocity of the shock front.If the shock is a strong one, then the jump conditions (2.10) become (2.11)un=(1-β)V,ρn=ρaβ,hn=haβ,pn=[(1-β)+12MA2(1-1β2)]ρaV2,vn=(A0eαR)R,
where MA=(ρaV2/μha2)1/2 is the Alfven-Mach number. The quantity β(0<β<1) is obtained by the relation (2.12)β2-β(γ(MA-2+1)-1(γ+1))+(γ-2)MA-2(γ+1)=0.Let the solution of (2.1)–(2.5) be of the form [4, 6, 19]
(2.13)u=1tU(η),ρ=tΩD(η),p=tΩ-2P(η),v=1tK(η),μh=t(Ω-2)/2H(η),
where
(2.14)η=teλr
while Ω and λ are constants. The variable η assumes constant value η0 at the shock surface. Hence
(2.15)V=-1λt,
which represents an outgoing shock surface, if λ<0.The solution of (2.1)–(2.5) in the form (2.13) to (2.15) are compatible with the shock conditions, if (2.16)Ω=2,λ=α=-δ2.
Since necessarily λ<0, relation (2.16) shows that δ>0, meaning thereby that the shock surface expands outwardly in an exponentially increasing medium.From (2.15) and (2.16), we obtain(2.17)R=2δlog⁡(tt0),
where t0 is the duration of the almost instantaneous explosion.


## 3. Solution to the Equations

The flow variables in the flow field behind the shock front will be obtained by solving the (2.1) to (2.5). From (2.13), (2.15), and (2.16), we obtain(3.1)∂u∂t=λuV-V∂u∂r,(3.2)∂ρ∂t=-2ρλV-V∂ρ∂r,(3.3)∂p∂t=-V∂p∂r,(3.4)∂h∂t=-V∂h∂r,(3.5)∂v∂t=λVv-V∂v∂r.Using (3.1) to (3.5) and the transformation, (3.6)r′=rR,u′=uV,v′=vV,ρ′=ρρn,h′=hhn,p′=ppnin the fundamental equations (2.1) to (2.5), we obtain(3.7)dρ′dr′=ρ′(1-u′)[du′dr′+2log⁡tt0+u′r′],(3.8)dp′dr′=ρ′[(1-β)β+(MA-2/2)(β-(1/β))]×[{(1-u′)-MA-2h′βρ′(1-u′)}du′dr′-u′(log⁡tt0)-MA-2h′2βρ′r′+v′2r′],(3.9)dh′dr′=h′(1-u′)du′dr′,(3.10)dv′dr′=1(1-u′)[(log⁡tt0)v′+u′v′r′],(3.11)du′dr′=(1-u′)[ρ′u′(log⁡(t/t0))+(MA-2h′2/βr′)-(v′2ρ′/r′)][ρ′(1-u′)2-(MA-2h′/β)-γp′β{(1-β)+(MA-2/2)(1-(1/β2))}].Also, the total energy of the disturbance is given by(3.12)E=2π∫r¯R‍ρ[U+12(u2+v2)+μh22ρ]rdr,
where r¯ is the position of inner boundary of the disturbance. Using (2.7), (3.4), and (2.11), (3.12) becomes
(3.13)E=8πμδ2η02R2∫r¯′1‍[{(1-u)+MA-22(1-1β2)}p′(γ-1)+(1/2)ρ′(u′2+v′2)β+MA-2h′22β]r′dr′.
Hence, the total energy of the shock wave is nonconstant and varies as R2. The increase of total energy may be achieved by the pressure exerted on the fluid by the inner expanding surface (a contact surface or a piston). A situation very much of the same kind may prevail during the formation of a cylindrical spark channel from exploding wires. In addition, in the usual cases of spark break down, time-dependent energy input is a more realistic assumption than instantaneous energy input [24, 25].In terms of dimensionless variablesr′, u′, ρ′, p′, h′, and v′ the shock conditions (2.11) take the form(3.14)r′=1,u′=(1-β),ρ′=1,p′=1,h′=1,v′=(A0η0)(log⁡tt0).
Equations (3.7) to (3.12) along with the boundary conditions (3.14) give the solution of our problem. The solution so obtained is a nonsimilar one, since the motion behind the shock can be determined only when a definite value for time is prescribed.


## 4. Results and Discussion

Distribution of the flow variables behind the shock front are obtained from (3.7) to (3.11). For the purpose of numerical integration, the values of the constant parameters are taken to be [16, 22, 23, 26, 27] γ=1.4; MA-2=0,0.05,0.1; t/t0=2,5; the value MA-2=0 corresponds to the nonmagnetic case. Starting from the shock front the numerical integration is carried out until the singularity of the solution (4.1)p′γβ{(1-β)+MA-22(1-1β2)}+MA-2h′β-ρ′(1-u′)2=0,
is reached. This marks the inner boundary of the disturbance and at this surface the value of r′(=r¯′) remains constant. The results are shown in Figures 1(a) to 2. These figures show that the rotation has significant effects on the flow variables. Value of r¯′ (the reduced position of the inner expanding surface) are shown in Table 1 for different cases.Table 1
The position of inner expanding surfacer-′ for MA-2 = 0, 0.05, 0.1 and t/t0 = 2, 5.
M
A
-
2
t
/
t
0
Position of inner expanding surface(r-′)
Rotating gas
Non-rotating gas
0
2
0.13100
0.21000
0
5
0.20000
0.37000
0.05
2
0.33678
0.37041
0.05
5
0.28146
0.39486
0.1
2
0.36408
0.39330
0.1
5
0.31025
0.41210Variation of reduced flow variables in the region behind the shock front forγ=1.4. (1)   MA-2 = 0, t/t0 = 2;   (2)  MA-2 = 0.05, t/t0 = 2;  (3)  MA-2 = 0.1, t/t0 = 2;  (4)  MA-2 = 0, t/t0 = 5;    (5) MA-2 = 0.05, t/t0 = 2;    (6)   MA-2 = 0.1, t/t0 = 5.
(a)
reduced densityρ′
(b)
reduced pressurep′
(c)
reduced radial velocityu′
(d)
reduced azimuthal velocityv′Figure 2
Variation of reduced magnetic fieldh′ in the region behind the shock front for γ=1.4.(1)  MA-2 = 0.05, t/t0 = 2;   (2)  MA-2 = 0.1, t/t0 = 2;  (3)  MA-2 = 0.05, t/t0 = 5;  (4)  MA-2 = 0.1, t/t0 = 5.From Table1, it is obvious that r¯′ increases by an increase in MA-2 in both the rotating and nonrotating cases but it increases with the increase of t/t0 in nonrotating case and decreases in rotating case except for MA-2=0.Figures1(a) and 1(b) show the distributions of reduced density and the reduced pressure, respectively. The density and pressure (except for MA-2=0, in the case of pressure) both decrease from the shock front and approach to zero near the inner expanding surface. In the nonmagnetic case, that is, for MA-2=0, the pressure increases and becomes constant after attaining a maximum value.Figure1(c) shows the distribution of the reduced radial velocity u′ in the flow field behind the shock front. The reduced radial velocity increases as we move from the shock front to the inner expanding surface, in general, whereas it becomes constant after attaining a maximum value in the nonmagnetic case.Figure1(d) shows that the reduced azimuthal velocity v′ has a similar behaviour as the reduced density ρ′.Figure2 shows that the reduced magnetic field h′ increases slowly behind the shock front and tends to infinity near the inner expanding surface in both the rotating and nonrotating cases.From Figures1(a) to 1(d) and 2, it is found that the effects of an increase in the ambient magnetic field (i.e. the effects of increasing MA-2) in both the rotating and nonrotating cases are(i)
a reduction of the radial velocity and a growth in the azimuthal velocity, in general,(ii)
to increase the pressure except forMA-2=0,(iii)
to increase the density,(iv)
to decrease the distance of the inner expanding surface from the shock front (see Table1).Effects of an increase in the time(t/t0) are  (i)
to decrease the density in both the the rotating and nonrotating cases, at a point in the flow field behind the shock,(ii)
to increase the azimuthal velocity,(iii)
to decrease the radial velocity and the azimuthal magnetic field in rotating case and to increase both in nonrotating case,(iv)
to decrease the pressure in rotating case and to increase that in nonrotating case,(v)
to increase the distance of the inner expanding surface from the shock font in rotating case except forMA-2=0 and to decrease that in nonrotating case.From Figures1(a) to 1(c), 2, and Table 1, it is found that the effects of rotation of the gas are   (i)
to decrease the radial velocity and the azimuthal magnetic field and to increase the density at any point in the flow field behind the shock,(ii)
to increase the pressure in the magnetic case and to decrease it in the nonmagnetic case,(i)
to increase the distance between the shock front and the inner expanding surface (see Table1).Present nonsimilar model may be used to describe some of the overall features of a “driven” shock wave produced by a flare energy releaseE (see (3.13)), that is time dependent. The energy E increases with time and the solutions then correspond to a blast wave produced by intense, prolonged flare activity in a rotating star when the wave is driven by fresh erupting plasma for some time, and its energy tends to increase as it propagates from the star into a cold atmosphere whose density varies exponentially with altitude. The atmospheric scale height of a star is generally small compared to its radius so that the solutions still describe a stellar explosion.



---

*Source: 168315-2012-08-05.xml*