Research Article  Open Access
G. Nath, A. K. Sinha, "A SelfSimilar Flow behind a Magnetogasdynamic Shock Wave Generated by a Moving Piston in a Gravitating Gas with Variable Density: Isothermal Flow", Physics Research International, vol. 2011, Article ID 782172, 8 pages, 2011. https://doi.org/10.1155/2011/782172
A SelfSimilar Flow behind a Magnetogasdynamic Shock Wave Generated by a Moving Piston in a Gravitating Gas with Variable Density: Isothermal Flow
Abstract
The propagation of a cylindrical (or spherical) shock wave in an ideal gas with azimuthal magnetic field and with or without selfgravitational effects is investigated. The shock wave is driven out by a piston moving with time according to power law. The initial density and the initial magnetic field of the ambient medium are assumed to be varying and obeying power laws. Solutions are obtained, when the flow between the shock and the piston is isothermal. The gas is assumed to have infinite electrical conductivity. The shock wave moves with variable velocity, and the total energy of the wave is nonconstant. The effects of variation of the piston velocity exponent (i.e., variation of the initial density exponent), the initial magnetic field exponent, the gravitational parameter, and the AlfvenMach number on the flow field are obtained. It is investigated that the selfgravitation reduces the effects of the magnetic field. A comparison is also made between gravitating and nongravitating cases.
1. Introduction
The explanation and analysis for the internal motion in stars is one of the basic problems in astrophysics. According to the observational data, the unsteady motion of large mass of the gas was followed by sudden release of energy results flareups in novae and supernovae. A qualitative behaviour of the gaseous mass may be discussed with the help of the equations of motion and equilibrium taking gravitational forces into account. Numerical solutions for selfsimilar adiabatic flows in selfgravitating gas were obtained by Sedov [1] and Carrus et al. [2], independently. Purohit [3] and Singh and Vishwakarma [4] have discussed homothermal flows behind a spherical shock wave in a selfgravitating gas using similarity method. Nath et al. [5] have studied the above problem assuming the flow to be adiabatic and selfsimilar and obtained the effects of the presence of a magnetic field. Shock wave through a variable density medium have been treated by Sedov [1], Sakurai [6], Rogers [7], Rosenau and Frankenthal [8], Nath et al. [5], Vishwakarma and Yadav [9], Nath [10], and others. Their results are more applicable to the shock formed in the deep interior of stars. Also, Vishwakarma and Singh [11] obtained the similarity solution for the flow behind a shock wave in a gravitating or nongravitating nonuniform gas with heat conduction and radiation heat flux in the case of adiabatic flow.
In aerodynamics, the analogy between the steady hypersonic flow past slenderblunted (planar or axisymmetric) powerlaw bodies and the onedimensional unsteady selfsimilar flow behind a shock driven out by a piston moving with time according to a powerlaw is well known (see, e.g., [1, 8, 10, 12–19]).
Since at high temperatures that prevail in the problems associated with the shock waves a gas is ionized, electromagnetic effects may also be significant. A complete analysis of such a problem should, therefore, consist of the study of the gasdynamic flow and the electromagnetic field simultaneously. A detailed study towards gaining a better understanding of the interaction between gasdynamic motion of an electrically conducting medium and magnetic field within the context of hyperbolic system has been carried out by many investigators such as Korobeinikov [20], Shang [21], and Lock and Mestel [22]. A detailed review in the field of magnetogasdynamic flows can be seen in the paper of Shang [21]. Lock and Mestel [22] analyzed the annular selfsimilar solutions in ideal magnetogasdynamics by casting the ideal magnetogasdynamic equations to a threedimensional autonomous system in which either the magnetic pressure or the fluid pressure vanishes.
In all of the works, mentioned above, the effect of selfgravitation is not taken into account by any of the authors in the case of isothermal flow with magnetic field.
The purpose of this study is, therefore, to obtain the selfsimilar solutions for the propagation of magnetogasdynamic cylindrical (or spherical) shock wave generated by a moving piston in a nonuniform gas with or without selfgravitational effects, in the presence of an azimuthal magnetic field, under isothermal flow condition. The media ahead and behind the shock front are assumed to be inviscid and to behave as thermally perfect gases. The counter pressure (the pressure ahead of the shock) is taken into account. The azimuthal magnetic field and the density in the ambient medium are assumed to vary and obey the power laws. The assumption of isothermal flow is physically realistic, when radiation heat transfer effects are implicitly present. As the shock propagates, the temperature behind it increases and becomes very large so that there is intense transfer of energy by radiation. This causes the temperature gradient to approach zero; that is, the dependent temperature tends to become uniform behind the shock front, and the flow becomes isothermal [10, 17, 20, 23–27]. A detailed mathematical theory of onedimensional isothermal blast waves in a magnetic field was developed by Lerche [28, 29]. Numerical solutions for the flow field between the shock and the piston are obtained in the case of isothermal flow in Section 3. Effects of viscosity and rotation are not taken into account.
2. Equations of Motion and Boundary Conditions
The fundamental equations governing the unsteady and cylindrically (or spherically) symmetric isothermal flow of an electrically conducting and selfgravitating ideal gas, in the presence of an azimuthal magnetic field, may be written as (c.f. [2, 10, 17, 20, 23–26, 30]) where and are independent space and time coordinates, is the fluid velocity, the density, the pressure, the azimuthal magnetic field, the temperature, the magnetic permeability, the mass contained in a unit cylinder of radius or in a sphere of radius , and the dimension of is taken as , and takes the values 2 and 1 for the respective cases of spherical and cylindrical symmetries, and the gravitational constant. In nongravitating case, (4) and the term in (2) do not occur.
The electrical conductivity of the gas is assumed to be infinite. Therefore, the diffusion term from the magnetic field equation is omitted, and the electrical resistivity is ignored. Also, the effect of viscosity on the flow of the gas is assumed to be negligible.
The above system of equations should be supplemented with an equation of state. An ideal gas behaviour of the medium is assumed, so that where is the gas constant, and is the ratio of specific heats.
A strong cylindrical (or spherical) shock wave is supposed to be propagating in the undisturbed ideal gas with variable density in the presence of an azimuthal magnetic field , where , and are constant.
The flow variables immediately ahead of the shock front are where , , is the shock radius, and the subscript “1” refers to the conditions immediately ahead of the shock.
The jump conditions at the magnetogasdynamic shock wave are given by the principle of conservation of mass, momentum, and energy across the shock; namely, where the subscript “2” denotes the conditions immediately behind the shock front, denotes the velocity of the shock front, and “” is the radiation heat flux.
From (8), we obtain where is the shockMach number referred to the frozen speed of sound , and is the AlfvenMach number. The quantity is obtained by the relation As the shock is strong, we assume to be negligible in comparison with the product of and [10, 17, 18, 23, 27, 31, 32]. Therefore, (10) reduces to Equation (5) together with (6) gives
3. SelfSimilarity Transformations
The inner boundary of the flow field behind the shock is assumed to be an expanding surface (piston). In the frame work of selfsimilarity [1], the velocity of the piston is assumed to follow a power law which results in [10, 16–19] where is the radius of the piston, and denotes the time at a reference state, is the piston velocity at , and is a constant. The consideration of the ambient pressure and the ambient magnetic field imposes a restriction on “”. Also, cylindrical and spherical geometries i.e., for do not permit for physical reasons. Thus, using (17), we obtain . For , the piston velocity suddenly rises at from zero to infinite velocity leading to the formation of a strong shock in the initial phase. The piston is then decelerated. Concerning the shock boundary conditions, selfsimilarity requires that the velocity of the shock is proportional to the velocity of the piston; that is, where “” is a dimensionless constant. Using (14), the time and space coordinate can be changed into a dimensionless selfsimilarity variable as Evidently, at the piston, and at the shock. To obtain the similarity solutions, we write the unknown variables in the following form [10, 16–18]: where , , , , and are functions of only.
For the existence of similarity solutions, “” and should be constants; therefore, Thus, where . Equation (18) shows that the solutions of the present problem can be reduced to the case in which the gas is nongravitating (i.e., the case in which ; the solution obtained by Nath [10] in the case of nonrotating medium for ).
Also, the total energy of the disturbance is given by Applying the similarity transformations (16) in (19), we obtain where , and is the gravitational parameter.
Hence, the total energy of the shock wave is nonconstant and varies as .
Equation (12) with the aid of (16) and (9) yields a relation between and in the form where .
By use of (16) and (21), (1) to (4) can be transformed and simplified to
From (22)–(25), we obtain Using the selfsimilarity transformations (16), the shock conditions (9) transform into The piston path coincides at with a particle path. Using (13) and (16), the relation can be derived. In addition to shock conditions (30), the kinematic condition (31) at the piston surface must be satisfied.
Normalizing the variables , , , , and with their respective values at the shock, we obtain
4. Results and Discussion
Distributions of the flow variables in the flow field behind the shock front are obtained by the numerical integration of (26)–(29) in the selfgravitating case and from (26), (27), and (29) in the nongravitating case with the boundary conditions (30) by the RungeKutta method of the fourth order. The expression for the shockMach number “” in the selfgravitating case with magnetic field is and in the nongravitating case with magnetic field , where “” is the exponent in the law of the variation of initial magnetic field. Also, the expression for the shockMach number “” in the nonmagnetic selfgravitating case (i.e., is and in the nongravitating and nonmagnetic case For the purpose of numerical integration, the values of the constant parameters are taken to be [8, 10, 33], , , . For fully ionized gas , and; therefore, it is applicable to stellar medium. Rosenau and Frankenthal [8] have shown that the effects of magnetic field on the flow field behind the shock are significant when ; therefore, the above values of are taken for calculations in the present problem. The value corresponds to the nonmagnetic case. The value corresponds to the solution in nongravitating case, the solution obtained by Nath [10] in nonrotating case when .
Table 1 shows the variation of density ratio across the shock front and the position of the piston for different values with , , in both the gravitating and nongravitating cases. Also, Table 1 shows that the distance of the piston from the shock front is less in the case of gravitating medium in comparison with that in the case of nongravitating medium. Physically, it means that the gas behind the shock is compressed in gravitating medium; that is, the shock strength is increased in gravitating medium.



Figures 1(a) to 1(c) and 2 show the variation of the flow variables , , , and , with at various values of the parameters , , and .
(a)
(b)
(c)
Figure 1(a) shows that the reduced fluid velocity increases from the shock front to the piston whereas it decreases when , and or 0.02, or in the case of nongravitating medium and or 0.02, or in the case of gravitating medium for .
Figure 1(c) shows that the reduced mass decreases from the shock front to the piston, and Figure 2 shows that the reduced azimuthal magnetic field increases from the shock front to the piston.
Figure 1(b) shows that the reduced density (pressure) decreases; but it increases in the cases when (i) , and (ii) , and starts to decrease after attaining a maximum in the cases (i) and the formation of maxima is absent in the cases (ii). As can be seen from (26) for “,” there is a singularity at the piston where because this equation becomes singular there. This singularity is removable, and a finite solution for “” is obtained as shown in Figure 1(b) (see curves 1–6) in nonmagnetic case . In the magnetic case in Figure 1(b), the curves 7–18 show that the singularity is nonremovable, and the derivative of density tends to negative infinity.
This singularity can be physically interpreted as follows [10, 16, 18]: the path of the piston diverges from the path of the particle immediately ahead rarifying the gas.
The piston position at which , is obtained after numerical integration of (26) to (29) in the selfgravitating case and from (26), (27), and (29) in the nongravitating case with the boundary conditions (30) and is tabulated in Table 1 for different values of , , and . The piston position is related with the velocity ratio of shock and piston, from (22), (23), and (24) as follows:
From Table 1 and Figures 1(a)–1(c) and 2, it is found that the effects of an increase in the value of (i.e., the effect of an increase in the strength of ambient magnetic field) are the following.(i)To decrease when or 10 (i.e., to increase the distance of the piston from the shock front), whereas to increase it when (i.e., to decrease the distance of the piston from the shock front). Physically it means that the gas behind the shock is less compressed when or 10, whereas it is more compressed when , that is, to decrease the shock strength in the case when or 10 and to increase it when (see Table 1). Thus, the effect of an increase in the strength of initial magnetic field is reduced by an increase in the strength of selfgravitation.(ii)To increase the value of (i.e., to decrease the shock strength) (see Table 1).(iii)To increase the reduced fluid velocity and the reduced azimuthal magnetic field in the case of selfgravitating medium but to decrease these flow variables in the case of nongravitating medium (see Figures 1(a) and 2).(iv)To increase the reduced mass but to decrease it when and or 100 (see Figure 1(c)).(v)To increase the reduced density (pressure) (see Figure 1(b)).
The effects of an increase in the value of the piston velocity exponent or the exponent for the variation of the initial density (i.e., the effects of the decrease in the value of the exponent for the variation of the initial magnetic field) of the ambient medium are the following.(i)To increase , that is, to decrease the distance of the piston from the shock front. Physically, it means that the gas behind the shock is more compressed; that is, the shock strength is increased (see Table 1).(ii)To increase the value of , when (i.e., to decrease the shock strength), which is same as in (i) (see Table 1).(iii)To increase the reduced fluid velocity and the reduced azimuthal magnetic field (see Figures 1(a) and 2).(iv)To decrease the reduced mass (see Figure 1(c)). (v)To increase the reduced density (pressure) in general (see Figure 1(b)).
Effects of an increase in the value of the gravitational parameter are(i)to increase , that is, to decrease the distance of the piston from the shock front (see Table 1),(ii)to increase the flow variables the reduced fluid velocity , the reduced density (pressure) , and the reduced azimuthal magnetic field (see Figures 1(a), 1(b) and 2), (iii)to decrease the reduced mass (see Figure 1(c)).
Acknowledgment
The authors are grateful to Dr. J. P. Vishwakarma, Professor of mathematics, DDU Gorakhpur University Gorakhpur273009, India, for many useful discussions.
References
 L. I. Sedov, Similarity and Dimensional Methods in Mechanics, Academic Press, New York, NY, USA, 1959.
 P. Carrus, P. Fox, F. Hass, and Z. Kopal, “The propagation of shock waves in a steller model with continuous density distribution,” The Astrophysical Journal, vol. 113, pp. 496–518, 1951. View at: Google Scholar
 S. C. Purohit, “Selfsimilar homothermal flow of selfgravitating gas behind shock wave,” Journal of the Physical Society of Japan, vol. 36, no. 1, pp. 288–292, 1974. View at: Google Scholar
 J. B. Singh and P. R. Vishwakarma, “Selfsimilar solutions in the theory of flareups in novae, I,” Astrophysics and Space Science, vol. 95, no. 1, pp. 99–104, 1983. View at: Publisher Site  Google Scholar
 O. Nath, S. Ojha, and H. S. Takhar, “A study of stellar point explosion in a selfgravitating radiative magnetohydrodynamic medium,” Astrophysics and Space Science, vol. 183, no. 1, pp. 135–145, 1991. View at: Publisher Site  Google Scholar
 A. Sakurai, “Propagation of spherical shock waves in stars,” Journal of Fluid Mechanics, vol. 1, pp. 436–453, 1956. View at: Google Scholar
 M. H. Rogers, “Analytic solutions for blast wave problem with an atmosphere of varying density,” Astrophysical Journal, vol. 125, pp. 478–493, 1957. View at: Google Scholar
 P. Rosenau and S. Frankenthal, “Equatorial propagation of axisymmetric magnetohydrodynamic shocks,” Physics of Fluids, vol. 19, no. 12, pp. 1889–1899, 1976. View at: Google Scholar
 J. P. Vishwakarma and A. K. Yadav, “Selfsimilar analytical solutions for blast waves in inhomogeneous atmospheres with frozeninmagnetic field,” European Physical Journal B, vol. 34, no. 2, pp. 247–253, 2003. View at: Publisher Site  Google Scholar
 G. Nath, “Magnetogasdynamic shock wave generated by a moving piston in a rotational axisymmetric isothermal flow of perfect gas with variable density,” Advances in Space Research, vol. 37, no. 9, pp. 1463–1471, 2011. View at: Google Scholar
 J. P. Vishwakarma and A. K. Singh, “A selfsimilar flow behind a shock wave in a gravitating or nongravitating gas with heat conduction and radiation heatflux,” Journal of Astrophysics and Astronomy, vol. 30, no. 1, pp. 53–69, 2009. View at: Publisher Site  Google Scholar
 S. S. Grigobian, “The Cauchy problem and piston problem for onedimensional unsteady motions of a gas: selfsimilar motions,” Journal of Applied Mathematics and Mechanics, vol. 22, no. 2, pp. 179–187, 1958. View at: Google Scholar
 N. N. Kochina and N. S. Mel'nikova, “On unsteady motion of a gas forced out by a piston with counter pressure neglected,” Journal of Applied Mathematics and Mechanics, vol. 22, no. 4, pp. 444–451, 1958. View at: Google Scholar
 K. C. Wang, “The piston problem with thermal radiation,” Journal of Fluid Mechanics, vol. 20, pp. 447–455, 1964. View at: Google Scholar
 J. B. Helliwell, “Self similar piston problem with radiative heat transfer,” Journal of Fluid Mechanics, vol. 37, no. 3, pp. 497–512, 1969. View at: Google Scholar
 H. Steiner and T. Hirschler, “A selfsimilar solution of a shock propagation in a dusty gas,” European Journal of Mechanics B/Fluids, vol. 21, no. 3, pp. 371–380, 2002. View at: Publisher Site  Google Scholar
 G. Nath, “Shock waves generated by a piston moving in a nonideal gas in the presence of a magnetic field: isothermal flow,” South East Asian Journal of Mathematics and Mathematical Sciences, vol. 5, pp. 69–83, 2007. View at: Google Scholar
 J. P. Vishwakarma and G. Nath, “A self similar solution of a shock propagation in a mixture of a nonideal gas and small solid particles,” Meccanica, vol. 44, no. 3, pp. 239–254, 2009. View at: Publisher Site  Google Scholar
 J. P. Vishwakarma and G. Nath, “Propagation of a cylindrical shock wave in a rotating dusty gas with heat conduction and radiation heat flux,” Physica Scripta, vol. 81, no. 4, Article ID 045401, 2010. View at: Publisher Site  Google Scholar
 V. P. Korobeinikov, Problems in the Theory of Point Explosion in Gases, Proceedings of the Steklov Institute of Mathematics Series No. 119, American Mathematical Society, 1976.
 J. S. Shang, “Recent research in magnetoaerodynamics,” Progress in Aerospace Sciences, vol. 37, no. 1, pp. 1–20, 2001. View at: Publisher Site  Google Scholar
 R. M. Lock and A. J. Mestel, “Annular selfsimilar solutions in ideal magnetogasdynamics,” Journal of Plasma Physics, vol. 74, no. 4, pp. 531–554, 2008. View at: Publisher Site  Google Scholar
 D. D. Laumbach and R. F. Probstein, “Selfsimilar strong shocks with radiation in a decreasing exponential atmosphere,” Physics of Fluids, vol. 13, no. 5, pp. 1178–1183, 1970. View at: Google Scholar
 P. L. Sachdev and S. Ashraf, “Converging spherical and cylindrical shocks with zero temperature gradient in the rear flow field,” Journal of Applied Mathematics and Physics, vol. 22, no. 6, pp. 1095–1102, 1971. View at: Publisher Site  Google Scholar
 S. Ashraf and Z. Ahmad, “Approximate analytic solution of a strong shock with radiation near the surface of the star,” Indian Journal of Pure and Applied Mathematics, vol. 6, pp. 1090–1098, 1975. View at: Google Scholar
 T. A. Zhuravskaya and V. A. Levin, “The propagation of converging and diverging shock waves under intense heatexchange conditions,” Journal of Applied Mathematics and Mechanics, vol. 60, no. 5, pp. 745–752, 1996. View at: Publisher Site  Google Scholar  MathSciNet
 G. Nath, “Propagation of a strong cylindrical shock wave in a rotational axisymmetric dusty gas with exponentially varying density,” Research in Astronomy and Astrophysics, vol. 10, no. 5, pp. 445–460, 2010. View at: Publisher Site  Google Scholar
 I. Lerche, “Mathematical theory of onedimensional isothermal blast waves in a magnetic field,” Australian Journal of Physics, vol. 32, pp. 491–502, 1979. View at: Google Scholar
 I. Lerche, “Mathematical theory of onedimensional cylindrical isothermal blast waves in a magnetic field,” Australian Journal of Physics, vol. 34, pp. 279–301, 1981. View at: Google Scholar
 G. B. Whitham, “On the propagation of shock waves through regions of nonuniform area or flow,” Journal of Fluid Mechanics, vol. 4, pp. 337–360, 1958. View at: Google Scholar
 J. P. Vishwakarma and G. Nath, “Similarity solutions for unsteady flow behind an exponential shock in a dusty gas,” Physica Scripta, vol. 74, no. 4, article 015, pp. 493–498, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 J. P. Vishwakarma and G. Nath, “Similarity solutions for the flow behind an exponential shock in a nonideal gas,” Meccanica, vol. 42, no. 4, pp. 331–339, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 J. P. Vishwakarma, A. K. Maurya, and K. K. Singh, “Selfsimilar adiabatic flow headed by a magnetogasdynamic cylindrical shock wave in a rotating nonideal gas,” Geophysical and Astrophysical Fluid Dynamics, vol. 101, no. 2, pp. 155–168, 2007. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2011 G. Nath and A. K. Sinha. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.