Skip to main content

Similarity method for the study of strong shock waves in magnetogasdynamics

Abstract

In this paper, a non-dimensional unsteady adiabatic flow of a plane or cylindrical strong shock wave propagating in plasma is studied. The plasma is assumed to be an ideal gas with infinite electrical conductivity permeated by a transverse magnetic field. A self-similar solution of the problem is obtained in terms of density, velocity and pressure in the presence of magnetic field. We use the method of Lie group invariance to determine the class of self-similar solutions. The arbitrary constants, occurring in the expressions of the generators of the local Lie group of transformations, give rise to different cases of possible solutions with a power law, exponential or logarithmic shock paths. A particular case of the collapse of an imploding shock is worked out in detail. Numerical calculations have been performed to obtain the similarity exponents and the profiles of flow variables. Our results are found in good agreement with the known results. All computational work is performed by using software package MATHEMATICA.

Introduction

The spread of shock waves under the control of strong magnetic field is a problem of great interest to researchers in a variety of fields such as nuclear science, geophysics, plasma physics and astrophysics. Hunter [[1]], Guderley [[2]], Greifinger and Cole [[3]] studied the problem of blast wave propagation in homogeneous and inhomogeneous media. Most recently, van Dyke and Guttmann [[4]], Sharma and Radha [[5], [6]], Madhumita and Sharma [[7]], Pandey et al. [[8]], Sharma and Arora [[9]], Arora et al. [[10], [11]] presented high accuracy results and alternative approaches for the investigation of blast wave by using the self-similar solutions method. In the same decade, a number of analytical solutions for the blast wave propagations have been obtained by Sachdev [[12]], Chisnell [[13]] and Singh et al. [[14], [15]]. Chisnell [[13]] provided analytical solutions to the problem of converging shock waves by using the singular points method. Singh et al. [[14], [15]] used the method of Lie group of transformations to obtain an exact solution for unsteady equation of non-ideal gas and magnetogasdynamics.

The magnetic fields have important roles in a variety of astrophysical situations. Complex filamentary structures in molecular clouds, shapes and the shaping of planetary nebulae, synchrotron radiation from supernova remnants, magnetized stellar winds, galactic winds, gamma-ray bursts, dynamo effects in stars, galaxies and galaxy clusters as well as other interesting problems all involve magnetic fields. When the internal disturbances accompanied by an increase in pressure take place in the central region of a star, a shock wave is formed. It travels from the central region to the periphery and emerges at the surface of the star. In the present paper, we consider the problem of propagation of a one-dimensional unsteady flow of an inviscid ideal gas permeated by a transverse magnetic field with infinite electrical conductivity as it approaches the surface of a star. It is assumed that mass density distribution in the medium follows a power law of the radial distance from the point of explosion.

In flows with imploding shocks, conditions of very high temperature and pressure can be produced near the center (axis) of implosion on account of the self-amplifying nature of imploding shocks. As a result of high temperatures attained by gases in motion, the effects of nonequilibrium thermodynamics on the dynamic motion of a converging shock wave can be important.

In this paper, we use the method of Lie group invariance under infinitesimal point transformations [[16]–[18]] to study the problem of propagation of strong shock waves in a radiating and electrically conducting gas permeated by a transverse magnetic field. The arbitrary constants, occurring in the expressions for the generators of the local Lie group of transformations, give rise to different cases of possible solutions with a power law, exponential or logarithmic shock paths. The Lie symmetry approach does not necessarily take into account the boundary and initial conditions unless the same are invariant under the change of variable transformations.

Basic equations

The basic equations describing the one-dimensional unsteady non-planar motion in which the direction of magnetic field is orthogonal to the trajectories of gas particles and electrical conductivity is infinite can be written as follows (Whitham [[19]]):

ρ t + ρ u x + u ρ x + m ρ u / x = 0 , u t + u u x + 1 ρ ( p x + h x ) = 0 , p t + u p x + γ p ( u x + m u x ) = 0 , h t + u h x + 2 h u x + 2 h m u / x = 0 ,
(1)

where p is the gas pressure, ρ is the density, u is the velocity, γ is the constant specific heat ratio, x is the single spatial co-ordinate, t is the time, h is the magnetic pressure defined by h=μ H 2 /2 with μ as magnetic permeability and H is the transverse magnetic field; m=0 and 1 correspond to planar and cylindrical symmetry, respectively, and the non-numeric subscripts denote the partial differentiation with respect to the indicated variable unless stated otherwise. The equation of state is taken to be of the form

p=ρRT.
(2)

Now, we consider the motion of a shock front, x=X(t), with the shock speed V=dX/dt propagating into a medium specified by

u=0,p= p 0 (x),h= h 0 (x),ρ= ρ 0 (x),
(3)

where T is the temperature, R is the gas constant, p 0 (x), h 0 (x) and ρ 0 (x) are some functions of x. The Rankine-Hugoniot jump conditions for the strong shocks are as follows (Whitham [[19]]):

u = 2 γ + 1 V , ρ = γ + 1 γ 1 ρ 0 ( X ( t ) ) , p = 2 γ + 1 ρ 0 ( X ( t ) ) V 2 1 2 ( γ + 1 ) 2 ( γ 1 ) 2 C 0 ρ 0 ( X ( t ) ) V 2 , h = 1 2 ( γ + 1 ) 2 ( γ 1 ) 2 C 0 ρ 0 ( X ( t ) ) V 2 ,
(4)

where C 0 is the shock Cowling number which is defined as C 0 = 2 h 0 ρ 0 V 2 , and the suffix 0 denotes the evaluation of the flow parameters just ahead of the shock.

Similarity analysis by invariance groups

Here, we suppose that there exists a solution of system (1) along a family of curves, called similarity curves, for which system (1) of partial differential equations reduces to a system of ordinary differential equations; this type of solution is called a similarity solution. In order to obtain the similarity solutions of system (1), we derive its symmetry group such that the system is invariant under this group of transformations. The idea of the calculation is to find a one-parameter infinitesimal group of transformations (see Sharma and Arora [[9]]):

x = x + ε χ ( x , t , ρ , u , p , h ) , t = t + ε ψ ( x , t , ρ , u , p , h ) , u = u + ε U ( x , t , ρ , u , p , h ) , ρ = ρ + ε S ( x , t , ρ , u , p , h ) , p = p + ε P ( x , t , ρ , u , p , h ) , h = h + ε E ( x , t , ρ , u , p , h ) ,
(5)

where the generators χ, ψ, U, S, P and E are to be determined in such a way that system (1) of partial differential equations together with conditions (3) and (4) is invariant with respect to transformations (5); the entity ε is so small that its square and higher powers may be neglected. The existence of such a group allows the number of independent variables in the problem to be reduced by one, and thereby allowing system (1) to be replaced by a system of ordinary differential equations.

In continuation, we shall use the summation convention and introduce the notation x 1 =t, x 2 =x, u 1 =ρ, u 2 =u, u 3 =p, u 4 =h and p j i = u i x j , where i=1,2,3,4 and j=1,2.

The system of basic equations (1), which is represented as

F k ( x j , u i , p j i ) =0,k=1,2,3,4,

is said to be constantly conformally invariant under the infinitesimal group (5) if there exist constants α k r (k,r=1,2,3,4) such that for all smooth surfaces, u i = u i ( x j ), we have

L F k = α k r F r ,
(6)

where L is the Lie derivative in the direction of the extended vector field

L= ξ j x j + η i u i + β j i p j i

with

ξ 1 =ψ, ξ 2 =χ, η 1 =S, η 2 =U, η 3 =P, η 4 =E

and

β j i = η i x j + η i u k p j k ξ l x j p l i ξ l u n p l i p j n ,
(7)

where l=1,2, j=1,2, i=1,2,3,4, n=1,2,3,4 and k=1,2,3,4. Here, repeated indices imply summation convention and β j i is the generalization of the derivative transformation.

System (1) implies

ξ j F k x j + η i F k u i + β j i F k p j i = α k r F r ,k=1,2,3,4,r=1,2,3,4.
(8)

Substitution of β j i from (7) into (8) yields a polynomial equations in the P j i . Setting the coefficients of P j i and P j i P l k to zero yields a system of first order, linear partial differential equations in the generators ψ, χ, S, U, P and E. This system, which is called the system of determining equations, is given as follows:

S t + u S x + ρ U x + m x ( ρ U + u S ρ u χ x ) = α 13 γ m p u x + α 11 m ρ u x + α 14 2 h m u x , S ρ ψ t u ψ x = α 11 , S u ρ ψ x = α 12 , S p = α 13 , S h = α 14 , U χ t + u S ρ u χ x + ρ U ρ = α 11 u ,
(9)
S + u S u + ρ ( U u χ x ) = α 11 ρ + α 12 u + α 13 γ p + α 14 2 h , ρ U p + u S p = α 12 ρ 1 + α 13 u , ρ U h + u S h = α 14 u + α 12 ρ 1 ; U t + u U x + ρ 1 ( P x + E x ) = α 21 m ρ u x + α 23 γ p m u x + α 24 2 h m u x , U ρ = α 21 , U u ψ t u ψ x = α 22 , U h ρ 1 ψ x = α 24 , U p ρ 1 ψ x = α 23 , U χ t + u ( U u χ x ) + ρ 1 ( P u + E u ) = α 21 ρ + α 22 u + α 23 γ p + α 24 2 h ,
(10)
u U ρ + ρ 1 ( P ρ + E ρ ) = α 21 u , S ρ 2 + u U p + ρ 1 ( P p χ x + E p ) = α 22 ρ 1 + α 23 u , S ρ 2 + u U h + ρ 1 ( P h χ x ) = α 22 ρ 1 + α 24 u ; P t + u P x + γ p U x + m γ x 1 ( u P + p U p u χ x 1 ) = α 31 m ρ u x 1 + α 33 m u x 1 + α 34 2 h m u / x , P ρ = α 31 , P u γ p ψ x = α 32 , P p ψ t u ψ x = α 33 , P h = α 34 , u P ρ + γ p U ρ = α 31 u ,
(11)
γ P + u P u + γ p ( U u χ x ) = α 31 ρ + α 32 u + α 33 γ p + α 34 2 h , U χ t + γ p U p + u P p u χ x = α 32 ρ 1 + α 33 u , u P h + γ p U h = α 32 ρ 1 + α 34 u ; E t + u E x + 2 h U x 2 h m u x 2 χ + 2 h m x 1 U + 2 u m x 1 E = α 41 m ρ u x 1 + α 43 m u x 1 + α 44 2 h m u x 1 , E ρ = α 41 , E h ψ t u ψ x = α 44 , E u 2 h ψ x = α 42 , E p = α 43 , 2 h U ρ + u E ρ = α 41 u , 2 E + 2 h U u 2 h χ x + u E u = α 41 ρ + α 42 u + α 43 γ p + α 44 2 h , 2 h U p + u E p = α 42 ρ 1 + α 43 u , U χ t + 2 h U h + u E h u χ x = α 44 u + α 42 ρ 1 , ψ = ψ ( x , t ) , χ = χ ( x , t ) .
(12)

Solving the above system of determining equations, we obtain the group of transformations as follows:

S = ( α 11 + a ) ρ , U = { ( α 22 + a ) u + k 1 , m = 0 , ( α 22 + a ) u , m = 1 , P = ( 2 α 22 + α 11 + 3 a ) p , E = ( α 44 + a ) h , χ = { ( α 22 + 2 a ) x + k 1 t + c , m = 0 , ( α 22 + 2 a ) x + c , m = 1 , ψ = a t + b ,
(13)

where α 11 , α 22 , a, b, c, and k 1 are the arbitrary constants.

Self-similar solutions

The arbitrary constants, which appear in the expressions for the infinitesimals of the invariant group of transformations, yield different cases of possible solutions as discussed below.

Case I. When a0 and α 22 +2a0, the change of variables from (x,t) to ( x ¯ , t ¯ ), defined as

x ¯ = { x if  m = 1 , x + c ( α 22 + 2 a ) 1 if  m = 0 , t ¯ =t+ b a ,
(14)

does not change system (1). Thus, rewriting the set of equations (13) in terms of the new variables x ¯ and t ¯ , and then suppressing the bar sign, we obtain

S = ( α 11 + a ) ρ , U = { ( α 22 + a ) u + k 1 , m = 0 , ( α 22 + a ) u , m = 1 , P = ( 2 α 22 + a 11 + 3 a ) p , E = ( α 44 + a ) h , χ = { ( α 22 + 2 a ) x + k 1 t , m = 0 , ( α 22 + 2 a ) x , m = 1 , ψ = a t .
(15)

The similarity variable and the form of similarity solutions for ρ, u, p and h readily follow from the invariant surface condition which yields

ψ ρ t +χ ρ x =S,ψ u t +χ u x =U,ψ p t +χ p x =P,ψ h t +χ h x =E.
(16)

The set of equations (16) together with (15) yields on integration the following forms of the flow variables:

ρ = t ( 1 + α 11 α ) S ( ξ ) , u = { t δ 1 U ( ξ ) if  m = 1 , t δ 1 U ( ξ ) k if  m = 0 , p = t ( 2 δ 1 + α 11 a ) P ( ξ ) , h = t ( 1 + α 44 a ) E ( ξ ) ,
(17)

where

k = k 1 ( δ 1 ) a andδ= α 22 + 2 a a .
(18)

The functions S , U , P and E depend only on the dimensionless form of the similarity variable ξ, which is determined as follows:

ξ= { x A t δ if  m = 1 , x A t δ + k t ( 1 δ ) A if  m = 0 ,
(19)

where A is a dimensional constant whose dimensions are obtained by the similarity exponent δ. Since the shock must be a similarity curve, it may be normalized to be at ξ=1. The shock path X and the shock velocity V are, then, given by

X= { A t δ if  m = 1 , A t [ t ( δ 1 ) k A ] if  m = 0 ,
(20)
V= { δ X t if  m = 1 , A δ t ( δ 1 ) k if  m = 0 .
(21)

At the shock, we have the following conditions on the functions S , U , P and E :

ρ | ξ = 1 = t ( 1 + α 11 a ) S ( 1 ) , u | ξ = 1 = { t δ 1 U ( 1 ) if  m = 1 , t δ 1 U ( 1 ) k if  m = 0 , p | ξ = 1 = t 2 δ 1 + α 11 a P ( 1 ) , h | ξ = 1 = t ( 1 + α 44 a ) E ( 1 ) .
(22)

Equations (22), in view of the invariance of jump conditions, suggest the following forms of ρ 0 (x) and h 0 (x):

ρ 0 ( x ) = ρ c ( x / x 0 ) θ , h 0 ( x ) = h c ( x / x 0 ) μ ,
(23)

and the following conditions on the functions S , U , P and E at the shock:

U ( 1 ) = 2 δ A γ + 1 , E ( 1 ) = 1 2 ( γ + 1 ) 2 ( γ 1 ) 2 C 0 ρ 0 h c A μ h 0 x 0 μ , P ( 1 ) = 2 ρ c δ 2 A 2 + θ ( γ + 1 ) x 0 θ 1 2 ( γ + 1 ) 2 ( γ 1 ) 2 C 0 ρ c δ 2 A 2 + θ x 0 θ , S ( 1 ) = γ + 1 γ 1 ρ c A θ x 0 θ ,
(24)

where

μ= α 44 + a α 22 + 2 a ,
(25)
θ= α 11 + a δ a ,δ= α 22 + 2 a a ,
(26)

where ρ c and h c are some reference constants associated with the medium.

Using (24), we rewrite equations (17) as follows:

ρ = ρ 0 ( X ( t ) ) S ( ξ ) , u = V U ( ξ ) , p = ρ 0 ( X ( t ) ) V 2 P ( ξ ) , h = ρ 0 ( X ( t ) ) V 2 E ( ξ ) .
(27)

Substituting (27) in system (1), we obtain the following system of ordinary differential equations in S , U , P and E , which on suppressing the asterisk sign becomes:

( U ξ ) S + S ( θ + U + m U ξ ) = 0 , ( δ 1 ) δ 1 U + ( U ξ ) U + S 1 ( P + E ) = 0 , 2 ( δ 1 ) δ 1 P + ( U ξ ) P + θ P + γ P ( U + m U ξ ) = 0 , 2 ( δ 1 ) δ 1 E + ( θ + 2 m U ξ 1 + 2 U ) E + ( U ξ ) E = 0 .
(28)

The jump conditions for the strong shock are:

U ( 1 ) = 2 γ + 1 , S ( 1 ) = γ + 1 γ 1 , P ( 1 ) = 2 γ + 1 1 2 ( γ + 1 ) 2 ( γ 1 ) 2 C 0 , E ( 1 ) = 1 2 ( γ + 1 ) 2 ( γ 1 ) 2 C 0 .
(29)

Case II. When a=0 and α 22 0, the change of variables from (x,t) to ( x ¯ , t ¯ ) is defined as

x ¯ =x+c ( α 22 ) 1 , t ¯ =t,

the basic equations in system (1) remain unchanged. The form of similarity solutions for the flow variables readily follows from (13), and can be expressed in the following form on suppressing the bar sign:

ρ = ρ 0 ( X ( t ) ) S ( ξ ) , u = V U ( ξ ) , p = ρ 0 ( X ( t ) ) V 2 P ( ξ ) , h = ρ 0 ( X ( t ) ) V 2 E ( ξ ) ,
(30)

together with the initial density and initial magnetic pressure:

ρ 0 (x)= ρ c ( x x 0 ) θ , h 0 (x)= h c ( x x 0 ) μ ,

where ξ is the dimensionless similarity variable, X(t) is the shock location and V is the shock velocity given by

ξ= x x 0 e δ t A ,X= x 0 e δ t A ,V= δ x 0 A e δ t A ,
(31)

with A as a dimensional constant. Substituting (30) in the equations in system (1) and using (31), we obtain the following system of ordinary differential equations in S , U , P and E , which on suppressing the asterisk sign becomes:

( U ξ ) S + S ( θ + U + m U ξ 1 ) = 0 , U + ( U ξ ) U + S 1 ( P + E ) = 0 , ( 2 + θ ) P + ( U ξ ) P + γ P ( U + m U ξ 1 ) = 0 , ( 2 + θ + 2 U + 2 m U ξ 1 ) E + ( U ξ ) E = 0 .
(32)

The jump conditions are

U ( 1 ) = { 2 γ + 1 if  m = 1 , 2 γ + 1 + k V if  m = 0 , S ( 1 ) = γ + 1 γ 1 , P ( 1 ) = 2 γ + 1 1 2 ( γ + 1 ) 2 ( γ 1 ) 2 C 0 , E ( 1 ) = 1 2 ( γ + 1 ) 2 ( γ 1 ) 2 C 0 .
(33)

Case III. When Q and α 22 +2a=0, the study reveals that this condition cannot be obtained in an axially symmetric (m=1) flow as it does not permit for the existence of a similarity solution in such a flow pattern. However, this condition can arise in a plane (m=0) flow where the change of variables from (x,t) to ( x ¯ , t ¯ ), defined as

x ¯ =x, t ¯ =t+ b a ,

does not change the basic equations in system (1). Accordingly, the similarity variable and the similarity solutions for the flow variables readily follow from (13), and can be expressed in the following form on suppressing the bar sign:

ρ = ρ 0 ( X ( t ) ) S ( ξ ) , u = V U ( ξ ) , p = ρ 0 ( X ( t ) ) V 2 P ( ξ ) , h = ρ 0 ( X ( t ) ) V 2 E ( ξ ) ,
(34)
ξ= [ x x 0 δ ln ( t A ) ] x 0 ,X= x 0 δln ( t A ) ,V= δ x 0 t ,
(35)
ρ 0 (x)= ρ c e θ ( x x 0 ) , h 0 (x)= h c ( x x 0 ) μ x .
(36)

Substituting (34), (35) and (36) in the equations in system (1), we obtain the following system of ordinary differential equations in S , U , P and E , which on suppressing the asterisk sign becomes:

( U 1 ) S + θ S + S U = 0 , ( U 1 ) U U δ + S 1 ( P + E ) = 0 , ( U 1 ) P + ( θ 2 δ ) P + γ P U = 0 , ( U 1 ) E + ( θ 2 δ + 2 U ) E = 0 .
(37)

The jump conditions are:

U ( 0 ) = 2 γ + 1 , S ( 0 ) = γ + 1 γ 1 , P ( 0 ) = 2 γ + 1 1 2 ( γ + 1 ) 2 ( γ 1 ) 2 C 0 , E ( 0 ) = 1 2 ( γ + 1 ) 2 ( γ 1 ) 2 C 0 .
(38)

Case IV. When a=0 and α 22 =0, this situation is similar to the previous case in the sense that it does not permit for the existence of a self-similar solutions in an axially symmetric flow. However, the plane flow involving a shock wave moving at constant speed admits a self-similar solution. Accordingly, the similarity variable and the similarity solutions for the flow variables follow from (13), and can be expressed in the following form:

ρ = ρ 0 ( X ( t ) ) S ( ξ ) , u = V U ( ξ ) , p = ρ 0 ( X ( t ) ) V 2 P ( ξ ) , h = ρ 0 ( X ( t ) ) V 2 E ( ξ ) , ξ = [ x x 0 δ t A ] x 0 , X = x 0 ( 1 + δ t A ) , V = x 0 δ A , ρ 0 ( x ) = ρ c e θ x x 0 x 0 , h 0 ( x ) = h c e μ x x 0 x 0 .
(39)

Substituting (39) in the equations in system (1) for m=0, and using (2), we obtain the following system of ordinary differential equations in S , U , P and E , which on suppressing the asterisk sign becomes:

( U 1 ) S + θ S + S U = 0 , ( U 1 ) U + S 1 ( P + E ) = 0 , ( U 1 ) P + θ P + γ P U = 0 , ( U 1 ) E + ( θ + 2 U ) E = 0 .
(40)

The jump conditions are:

U ( 1 ) = 2 γ + 1 , S ( 1 ) = γ + 1 γ 1 , P ( 1 ) = 2 γ + 1 1 2 ( γ + 1 ) 2 ( γ 1 ) 2 C 0 , E ( 1 ) = 1 2 ( γ + 1 ) 2 ( γ 1 ) 2 C 0 .
(41)

Imploding shocks

Here, we consider Case I of an imploding strong shock in the neighborhood of implosion. For the problem of a converging shock collapsing at the axis, the origin of time t is taken to be the instant at which the shock reaches the axis so that t0 in (28). Therefore, the definition of the similarity variable is slightly modified by setting

X=A ( t ) δ ,ξ=x/A ( t ) δ ,
(42)

so that the intervals of the variables are <t0, Xx< and 1ξ<. At the instant of collapse (t=0), the gas velocity, pressure, density and the sound speed at any finite axial distance x are bounded, but with t=0 and finite x, ξ=. In order for the quantities u, p, ρ and h to be bounded when t=0 and x is finite, we must have the following boundary conditions at ξ=:

U()=0, P ( ) S ( ) =0,E()=0.
(43)

In the matrix notation, system (28) can be written as

C W =B,
(44)

where W= ( U , S , P , E ) t r , and the matrix C and the column vector B can be identified by observing system (28). In system (28), there is an unknown parameter δ, which cannot be obtained from an energy balance or the dimensional considerations; it is computed only by solving a non-linear eigenvalue problem for a system of ordinary differential equations. The range of similarity variable is 1ξ< for the implosion problem, and system (44) can be solved for the derivatives U , S , P and E in the following form:

U = Δ 1 Δ , S = Δ 2 Δ , P = Δ 3 Δ , E = Δ 4 Δ ,
(45)

where Δ, defined as the determinant of the matrix C, is given by

Δ= ( U ξ ) 2 [ ( U ξ ) 2 ( γ P + 2 E ) S ] ,
(46)

and Δ k (k=1,2,3,4) are the determinants obtained from Δ by replacing the k th column by the column vector B. It can be observed that U<ξ in the interval [1,), while Δ is positive at ξ=1 and negative at ξ= indicating that there exists a ξ[1,) at which Δ vanishes, and consequently the solutions become singular. In order to get a non-singular solution of (44) in the interval [1,), we select the value of the exponent δ such that Δ vanishes only at the points where the determinant Δ 1 is zero too. It can be checked that at points where Δ and Δ 1 vanish, the determinants Δ 2 , Δ 3 and Δ 4 also vanish simultaneously. To find the exponent δ in such a manner, we introduce the variable Z as follows:

Z(ξ)= ( U ξ ) 2 γ P ( ξ ) + 2 E ( ξ ) S ( ξ ) .
(47)

This, in view of (45), implies

Z = { 2 ( U ξ ) ( Δ 1 Δ ) γ Δ 3 + 2 Δ 4 S + γ P + 2 E S 2 Δ 2 } /Δ.
(48)

Equations (45), in view of (48), become

d U d Z = Δ 1 Δ 5 , d S d Z = Δ 2 Δ 5 , d P d Z = Δ 3 Δ 5 , d E d Z = Δ 4 Δ 5 ,
(49)

where Δ 5 =2(Uξ)( Δ 1 Δ) γ Δ 3 + 2 Δ 4 S + γ P + 2 E S 2 Δ 2 , with ξ=U+ { Z + γ P + 2 E S } 1 / 2 .

Numerical results and discussion

We integrate equations (45) from the shock Z=Z(1) to the singular point Z=0 by choosing a trial value of δ, and compute the values of U, S, P, E and Δ 1 at Z=0; the value of δ is corrected by successive approximations in such a way that for these values, the determinant Δ 1 vanishes at Z=0. The values of the similarity exponent δ, obtained from the numerical calculations for different values of C 0 , m and θ are given in Table 1.

Table 1 Similarity exponent δ for planar and cylindrically symmetric flows and the density exponent θ with γ=1.4

It may be noticed that the magnetic field effects enter through the parameter C 0 . The values of the flow variables before collapse and at the instant of collapse are depicted in Figures 1 and 2.

Figure 1
figure 1

Flow patterns: (a) density, (b) pressure, (c) temperature, (d) velocity for m=0 (planar flow) and γ=1.4 and C 0 =0 .

Figure 2
figure 2

Flow patterns: (a) density, (b) pressure, (c) temperature, (d) velocity for m=1 (cylindrically symmetric flow) and γ=1.4 and C 0 =0 .

The typical flow profiles show that the density, pressure, temperature and velocity increase behind the shock wave with the increase in the value of θ; this is because a gas particle passing through the shock is subjected to a shock compression. Indeed, this increase in pressure and density behind the shock may also be attributed to the geometrical convergence or area contraction of the shock wave. Figures 1 and 2 show that the growth of the flow variables is slower in cylindrical symmetry as compared with that in planar symmetry. Figures 1 and 2 also confirm the generation of higher pressure near the axis of symmetry, i.e., near ξ=. The difference between flow profiles in cylindrical waves and those in planar waves is attributed to the fact that for planar waves, the flow distribution is relatively less influenced by the interaction between the gasdynamic phenomena as compared to cylindrical waves.

Conclusions

In the present investigation a self-similar method is used to study the flow pattern behind an exponential shock driven by a piston in ideal magnetogasdynamics. The general behavior of density, velocity and pressure profiles remains unaffected due to presence of magnetic field in ideal gas. However, there is a decrease in values of density, velocity and pressure in the case of magnetogasdynamics as compared to non-magnetic case. It may be noted that the effect of magnetic field on the flow pattern is more significant in the case of isothermal flow as compared to that of adiabatic flow.

References

  1. Hunter C: Similarity solutions for the flow into a cavity. J. Fluid Mech. 1963, 15: 289-305. 10.1017/S0022112063000252

    Article  MathSciNet  Google Scholar 

  2. Guderley G: Starke kugelige und zylindrische Verdichtungsstosse in der Nahe des Kugelmittelpunktes bzw. der Zylinderachse. Luftfahrt-Forsch. 1942, 19: 302-312.

    MathSciNet  Google Scholar 

  3. Greifinger C, Cole JD: Similarity solutions for cylindrical magnetohydrodynamic blast waves. Phys. Fluids 1962, 5: 1597-1607. 10.1063/1.1706571

    Article  MathSciNet  Google Scholar 

  4. van Dyke M, Guttmann AJ: The converging shock wave from a spherical or cylindrical piston. J. Fluid Mech. 1982, 120: 451-462. 10.1017/S0022112082002845

    Article  MathSciNet  Google Scholar 

  5. Sharma VD, Radha C: On one-dimensional planar and non-planar shock waves in a relaxing gas. Phys. Fluids 1994, 6: 2177-2190. 10.1063/1.868220

    Article  MathSciNet  Google Scholar 

  6. Sharma VD, Radha C: Similarity solutions for converging shocks in a relaxing gas. Int. J. Eng. Sci. 1995, 33: 535-553. 10.1016/0020-7225(94)00086-7

    Article  MathSciNet  Google Scholar 

  7. Madhumita G, Sharma VD: Propagation of strong converging shock waves in a gas of variable density. J. Eng. Math. 2003, 46(1):55-68. 10.1023/A:1022816118817

    Article  MathSciNet  Google Scholar 

  8. Pandey M, Sharma VD, Radha R: Symmetry analysis and exact solution of magnetogasdynamic equations. Q. J. Mech. Appl. Math. 2008, 61: 291-310. 10.1093/qjmam/hbn011

    Article  MathSciNet  Google Scholar 

  9. Sharma VD, Arora R: Similarity solutions for strong shocks in an ideal gas. Stud. Appl. Math. 2005, 114: 375-394. 10.1111/j.0022-2526.2005.01557.x

    Article  MathSciNet  Google Scholar 

  10. Arora R, Tomar A, Singh VP: Similarity solutions for strong shocks in a non-ideal gas. Math. Model. Anal. 2012, 17(3):351-365. 10.3846/13926292.2012.685957

    Article  MathSciNet  Google Scholar 

  11. Arora R, Siddiqui MJ, Singh VP: Similarity method for imploding strong shocks in a non-ideal relaxing gas. Int. J. Non-Linear Mech. 2013, 57: 1-9. 10.1016/j.ijnonlinmec.2013.06.009

    Article  Google Scholar 

  12. Sachdev PL: Propagation of a blast wave in uniform or non-uniform media: a uniformly valid analytic solution. J. Fluid Mech. 1972, 52: 369-378. 10.1017/S0022112072001478

    Article  MathSciNet  Google Scholar 

  13. Chisnell RF: An analytic description of converging shock waves. J. Fluid Mech. 1998, 354: 357-375. 10.1017/S0022112097007775

    Article  MathSciNet  Google Scholar 

  14. Singh LP, Husain A, Singh M: An analytical study of strong non-planar shock waves in magnetogasdynamics. Adv. Theor. Appl. Mech. 2010, 6: 291-297.

    Google Scholar 

  15. Singh LP, Husain A, Singh M: An approximate analytical solution of imploding strong shocks in a non-ideal gas through Lie group analysis. Chin. Phys. Lett. 2010., 27: 10.1088/0256-307X/27/1/014702

    Google Scholar 

  16. Bluman GW, Kumei S: Symmetries and Differential Equations. Springer, New York; 1989.

    Book  Google Scholar 

  17. Bluman GW, Cole JD: Similarity Methods for Differential Equations. Springer, Berlin; 1974.

    Book  Google Scholar 

  18. Logan JD, Perez JDJ: Similarity solutions for reactive shock hydrodynamics. SIAM J. Appl. Math. 1980, 39: 512-527. 10.1137/0139042

    Article  MathSciNet  Google Scholar 

  19. Whitham GB: Linear and Nonlinear Waves. Wiley-Interscience, New York; 1974.

    Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mohd Junaid Siddiqui.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Arora, R., Yadav, S. & Siddiqui, M.J. Similarity method for the study of strong shock waves in magnetogasdynamics. Bound Value Probl 2014, 142 (2014). https://doi.org/10.1186/s13661-014-0142-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-014-0142-2

Keywords