On the Effects of Viscosity on the Shock Waves for a Hydrodynamical Case—Part I: Basic Mechanism

The interaction of shock waves with viscosity is one of the central problems in the supersonic regime of compressible fluid flow. In this work, numerical solutions of unmagnetised fluid equations, with the viscous stress tensor, are investigated for a onedimensional shock wave. In the algorithm developed the viscous stress terms are expressed in terms of the relevant Reynolds number. The algorithm concentrated on the compression rate, the entropy change, pressures, and Mach number ratios across the shock wave. The behaviour of solutions is obtained for the Reynolds and Mach numbers defining the medium and shock wave in the supersonic limits.


Introduction
Shock waves are rather common phenomena in the supersonic flows of any fluid.They arise in many areas that are related with hydrodynamics such as fluid mechanics, aerodynamics, astrophysics, solar physics, and space physics.If a medium is shocked, particles behind the shock front experience both compressive and shear forces.They push the particles away from their original equilibrium positions.
Shock waves are studied by many authors.Somow and Spector [1] studied the basic mechanisms of the hydrodynamic shocks in the solar atmosphere during flares.Effects of inhomogeneities in solar wind plasma on the interplanetary shock waves are also studied by Heinemann [2].Their observational studies are given in the work of De Lucas et al. [3].Magara and Shibata [4] worked on the formation of shock waves due to plasma ejections in the solar atmosphere.Khidr and Mahmoud [5] obtained results for the case of an arbitrary Prandtl number in strong shock waves, using a modified power law for viscosity in terms of temperature and Mach number.Kuznetsov [6] studied the stabilities of shock waves in hydrodynamic discontinuities and flows in the relaxation zone.On the other hand, Hamad and El-Fayes [7] studied the entropy change for the structure of inviscid plasma by neglecting viscosity in the gas phase.Compression rate in bow shocks was studied by Kabin [8], and the Euler equations for a one-dimensional hydrodynamic model were considered [9].In a more recent study, Swift et al. [10] derived expressions for shock formation, based on the local curvature of the flow characteristics during dynamic compression.
Motion occurs in a continuous medium because of some external causes like pressure gradients and body forces; this way, the fluid velocity distribution becomes inhomogeneous in general.In consequence, resistance to variations in the distribution of cohesive forces in fluids may such come into play to remove inhomogeneities in velocities.Such resistive effects produce the phenomenon of viscosity in fluid motions [11].Therefore, viscosity, which measures the resistance of a fluid system as a function of velocity, should also be taken into account as one of the most important effects in the equations of motion.For the shock heating of solar corona, Orta et al. [12] obtained that the shock thickness and profile depend on the viscosity and resistivity, and heating ultimately occurs because of them.Dispersive shock waves were studied by Ballai et al. [13].They showed that dispersion will alter the amplitude and propagation speed of a shock wave.
For supersonic flows the coexistence of shock waves with viscous effects produces interesting features.The interaction of viscous effects with a shock wave is a key feature in many fluid dynamics systems [14].Considerable theoretical and Advances in Astronomy numerical interest should be expressed in the strengthening of a shock wave propagating through a viscous medium.The viscous interactions of a solar wind stream were studied by Korzhov et al. [15] and it was found that the Kelvin-Helmholtz instability is excited due to the presence of shear flows.Coronal mass ejections (CMEs) and solar wind are the main results of the solar activity.These events can drive interplanetary shock waves and produce geomagnetic storms.Study of these events is very important for space weather purposes.The shock waves occur where the solar wind changes from being supersonic (with respect to the surrounding interplanetary medium) to being subsonic.In the supersonic regime of compressible gas flow, the interaction of shock waves with viscosity is one of the central problems.To define such type of shock process, the Navier-Stokes equations should be solved.Mathematically this process can be approximated to a hydrodynamical case.
The Reynolds number is another defining feature of fluid motion.Reynolds, a British scientist, showed that the transition from laminar to turbulent flow is directly related to a dimensionless number defined as the Reynolds number [16].This dimensionless number gives the relation between inertial and viscous forces in a fluid flow.If the inertial forces are large (the flow rate is high), then the flow will rather occur turbulently.When the viscous forces are large enough compared to the inertial ones, then a laminar fluid flow character will be observed.
The main purpose of this paper is to describe the basic mechanism of the shock wave problem, in which the viscous terms cannot be neglected in deriving the jump relations.To define this type of shock process, the Navier-Stokes equations defining the momentum conservation should be solved.The viscous behaviour of a fluid can be described by using the Reynolds number.The formalism used for a one-dimensional study is explained in more detail in the following section.Solutions are presented in terms of physical parameters such as compression rate, Mach number and pressure ratios, and entropy change across the shock waves.Results will be given in Section 3. Our results are compared with other similar works in Section 4, together with a discussion and conclusion.

Physical Formulation of the Problem
2.1.Basic Equations.The basic formulation for unmagnetised plasma is also known as the system of hydrodynamical equations in conservative form.For a compressible viscous shock wave for steady flow, the equations, are [17][18][19] These basic equations to be solved consist of the conservation of the mass (1), the momentum (2), and the energy for adiabatic flow (3) together with the ideal gas law (4).Here , , ,   ,   ,   , , and  denote the density, velocity, infinitesimal surface area, unit normal ( = 1, 2, 3), the Kronecker delta, viscous stress tensor, pressure, and the total enthalpy, respectively.

One-Dimensional Problem.
Using these steady viscous flow relations with letting subscript 1 and 2 as the upstream and downstream shock front, respectively, fluid equations can be written for a one-dimensional case (say along -axes).The Rankine-Hugoniot jump relations [20] can be obtained through ( 1)-( 3).We use the viscous tensor in one dimension as where  is the dynamic viscosity coefficient.We define the compression ratio, , in terms of the upstream and downstream velocities: We can define the Reynolds number, Re, through With the use of ( 5)-( 7) in ( 1)-( 4), we can obtain the equation as follows (after a little algebra): Here Re 1 and Re 2 are the upstream and downstream Reynolds number, respectively, and  1 is the value of the upstream Mach number.The last equation is also called the general equation and it can be solved numerically by imposing the appropriate conditions.Since 1 <  < 2, this equation has just one positive root.Moreover the solution reduces to the inviscid hydrodynamic value for the vanishing values of dynamic viscosity.

Value of the Downstream Reynolds Number.
In order to simplify the problem, the downstream Reynolds number, Re 2 , can be expressed as a function of the adiabatic index, Re 1 , and  1 is similar to Bruhn et al. [21].Figure 1 shows the change of Re 2 with respect to Re 1 for different values of  1 for the monatomic with gas  = 5/3.There is a limitation on the upstream Mach number ( 1 ) in order to have a shock wave.Shock waves are common phenomena in supersonic upstream flows ( 1 > 1) of any fluid.From the quadratic structure of ( 8) and the limitation on the value

Model Results
We have attempted to find some special solutions of ( 8) and the Rankine-Hugoniot jump relations using the method given in the previous section with an algorithm developed and adapted to Maple 9.5.Table 1 represents the variation of the parameters describing the physical structure for the inviscid shock waves.These parameters are compression rate (i.e., densities ratio), the Mach numbers ratio ( 2 / 1 ), pressures ratio ( 2 / 1 ), and the entropy difference ( 2 −  1 ) provided by the solutions of inviscid forms of the hydrodynamic equations.These values will be used as a reference model.Finally, one can easily conclude that vanishing values of the viscosity means higher values of Reynolds number from its definition.As explained in the previous section, distributions of physical parameters can be obtained from the solutions of ( 8) and the Rankine-Hugoniot relations.The related quantities, Mach numbers and pressures ratios, compression rate, and entropy change are presented in Figures 3-7 and Table 2 for a monatomic gas with  = 5/3.Figures 3 and 4 show the changes in the compression rate () with respect to the upstream Mach number,  1 , and Re 2 /Re 1 , respectively.Figure 3 was drawn for two different values of Re 1 .It shows that, for weak shocks ( 1 < 2), there is no variation with respect to upstream Mach number in accordance with Khidr and Mahmoud [5].However, for strong shocks ( 1 > 2), there are considerable differences between the values of compression rate for different values of Re 1 .Another important result of strong shocks is the fact that the value of  approaches its inviscid value (3.368) as given in Table 1 for the higher values of Re 1 (e.g., 4500) shown as the upper dashed line in Figure 3.
Figure 4 represents the Re 2 /Re 1 dependence of .It tends to decrease with increasing values of ratio.It also exhibits an interesting feature.For the small values of the Reynolds numbers ratio, it gives a larger compression rate at the smaller value of Re 1 (i.e., 50).It has comparable values in the vicinity of the Re 2 /Re 1 = 1 for both high and low values of Re 1 .However, for the larger values of ratio, the shape is visibly changed, and it gives larger compression rates at the larger value of Re 1 = 4500.
Table 2 summarizes the results obtained for pressures ratio ( 2 / 1 ) with the values of the upstream Mach number in the supersonic region for two different values of Re 1 . 2 / 1 approaches its inviscid value ( 2 / 1 = 19.750)given in Table 1 for the higher values of Re 1 in the strong shock regime.
In Figure 5, variation of  2 / 1 is presented as a function Re 2 /Re 1 for two different values of the upstream Reynolds number.It has the tendency of growing with the   The entropy change depicted in Figures 6 and 7 can be expressed in terms of pressures ratio and compression rate as follows [22]: In Figure 6, variation of  2 −  1 with respect to Re 2 / Re 1 ratio has the tendency to decrease with increasing ratio.For higher values of Re 1 ,  2 −  1 has larger values as expected since higher values of the Reynolds number close the case of inviscid flow value given in Table 1.In Figure 7, dependence of entropy change,  2 −  1 , with respect to  1 is given.It has the tendency to increase with increasing  1 .For the higher values of Re 1 (i.e., vanishing values of viscosity),  2 −  1 has larger values as expected.

Variations with Respect to Reynolds Number.
In this subsection, the distributions of physical parameters are given as a function of the upstream Reynolds number, Re 1 , in Figures 8 and 9 and Table 3 again for a monatomic gas.
In Figure 8, the values of compression rate () are given as a function of Re 1 for different values of  1 .It is increasing with the values of the upstream Mach number and also reaches very high values for strong shock cases of  1 = 2.5 and 4.
From Figure 9, it is seen that the entropy change is negative for smaller values of  1 ≤ 1.2 in the region of very small Re 1 (i.e., <50) across the shock.This case violates the second law of thermodynamics (i.e.,  2 −  1 < 0).Thus the weak shock waves (i.e.,  1 = 1.2) are nearly isentropic.
In Table 3, Re 1 dependencies of the Mach numbers and pressures ratios are given for different values of M 1 .One of the most important features of these results is that the values of  2 / 1 have a decreasing tendency with the increasing values of Re 1 for Re 2 /Re 1 > 1 and a decreasing tendency for Re 2 /Re 1 < 1.The critical value of the Mach number for this turning point was found to be 2.045 which equates the Reynolds numbers ratio to unity as in Figure 2.This point is important not only for the Reynolds numbers ratio but also for the strength of shock waves.The pressures ratio ( 2 / 1 ) has the increasing tendency and reaches a value close to the inviscid flow value given in Table 1 for the higher values of the Reynolds number with increasing values of Mach number (i.e., strong shocks).

Discussion and Conclusion
When a fluid is shocked, particles after the shock front experience both compressive and stress forces.In such a movement viscous forces largely depend on plasma particles motions.In the present work, the detailed analysis of the structure of steady hydrodynamical equations including the viscous stress tensor in conservative forms was attempted in order to investigate the effects of the Reynolds number on the shock waves.
Since the Reynolds number gives the relation between inertial and viscous forces of the fluid flow, the character of flow can be determined from the Reynolds number.If inertial forces are dominant and the Reynolds number Re > 2000, flow is turbulent.If viscous forces surpass, then the value of the Reynolds number Re < 2000 and flow becomes laminar.Using fluid-velocity dependence, viscosity can change the character of flow.As expected, a high Reynolds number can be a result of vanishing viscosity and its lower values are the result of higher viscosity.Since the viscous stress is written as a function of fluid velocity, the viscosity can be determined by the change in the fluid velocity.Such a dissipative mechanism controls the values of flow variables.The viscous terms in the Rankine-Hugoniot jump relations can be expressed in terms of the Reynolds number using its definition.In this present paper, the basic mechanism of the problem is given in order to investigate the effects of viscosity and the Reynolds number.The behaviour of shock waves in supersonic flows with varying Reynolds and Mach numbers was examined using analytical and numerical methods.
We obtained that, in strong shocks, physical parameters have the increasing tendency to reach a value close to the inviscid flow values obtained for the higher values of the upstream Reynolds number making the flow turbulent.Depending on increasing values of viscosity, the values parameters decrease across the shock.
As in Mace and Adamson Jr. [23], for the shock waves in a transonic region, compression rate is independent of the Reynolds number.The results in this regime agree well with the inviscid values.These results are in good concordance with the results obtained in the work of Kabin [8] for weak shocks.Viscosity makes the profiles of compression ratio smooth for very high values of the Mach number.The compression rate is increasing with the values of the upstream Mach number.However, there are considerable differences between the compression rates for different values of Re 1 for strong shocks.That is, compression rate is higher for laminar upstream flow and is lower for turbulent upstream flow.
The Mach numbers ratio has a decreasing tendency for the larger values of Re 1 in the region of Re 2 /Re 1 > 1.

Advances in Astronomy
In other words, a shock wave becomes a Prandtl-Meyer expansion wave for very high values of Re 2 /Re 1 for the transonic values of the upstream Mach number (<2) as it is seen from Figure 1. 2 / 1 decreases for Re 2 /Re 1 < 1 in the region of  1 > 2. When the upstream Reynolds number is dominant, the shock waves become very strong (see Figure 2).Viscosity makes the problem smooth, and the downstream shock waves become more turbulent if Re 1 is turbulent in the transonic limits of  1 in accordance with the results of Jamme et al. [24] As given in the work of Gamba [25], for the turbulent case of Re 1 (i.e., >2000), we obtained that the entropy change increases with increasing values of upstream shock speed.It is also clear that the entropy change decreases with the increasing values of the Reynolds numbers ratio Re 2 /Re 1 and it has close results for both extreme values of the Re 1 .It is also found that the weak shocks are nearly isentropic for smaller values of upstream Re similar to the results of De Sterck et al. [26].
In summary, the present study provides results useful for future studies including dissipative shocks.These types of shock waves arise frequently in the cases where the velocity of fluid is greater than the local sound speed; they find an application in gas dynamics, fluid mechanics, aerodynamics, astrophysics, solar physics, and space physics, for both magnetised and unmagnetised fluid motions.
It is planned to use the present algorithm and its results to be applied to solar coronal shock waves.The effects of viscous flows on the shock wave happened after 13/12/2006-CME is under study in order to apply the present results.The simulations of the effects of the Reynolds number on the shock waves for magnetized plasma in both one and multidimensions are also under study.

2 Figure 1 :
Figure 1: Behaviour of the Re 2 as a function of Re 1 for different values of  1 .

Figure 6 :
Figure 6: Distributions of entropy difference  2 −  1 as a function of Re 2 /Re 1 for two different values of Re 1 .

Figure 7 :
Figure 7: Same as Figure 6, but as a function of  1 .

2 Figure 8 :
Figure 8: Variations of compression rate as a function of Re 1 at four different upstream Mach number values.

1 Figure 9 :
Figure 9: Variations of entropy difference  2 −  1 with Re 1 for four different values of  1 .
3.1.Variations of Parameters with respect to Mach Numbers.

Table 3 :
Reynolds number, Re 1 , dependence of  2 / 1 and  2 / 1 for different values of  1 and Re 2 /Re 1 values.Figure 3: Variations of compression rate with respect to  1 for different values of upstream Re.Figure 4: Variations of compression rate with respect to Re 2 /Re 1 .