[4] We study one-dimensional stationary disturbance in weakly ionized nonisothermal plasma with Te Ti Tn generated by a strong shock wave of neutral component. We take this wave in a form
(1) |
(2) |
(3) |
(4) |
(5) |
[5] In undisturbed conditions we will consider the plasma as a weakly ionized one: d ni0nn0-1 1. This fact makes it possible to simplify the description of the process of interaction of the shock wave with plasma components. Solving the problem at the first stage, a step-by-step consideration of the process is possible for obtaining estimates. Taking into account weak ionization of the plasma, we first neglect the retroaction of charged components on the plasma neutral components.
[7] Taking into account that the process is stationary (the fields depend on one variable x ), the equations in partial derivatives (1)-(5) in the case hi=0 (the finite viscosity coefficient h will be taken into account below) are reduced to the system of two connected ordinary differential equations for the functions Vi and Y.
(6) |
(7) |
[8] There are two limiting situation. The first situation corresponds to the case De2 L2 ( De2 is a small parameter in the approximation of long waves, L is a special scale of the fields variation). The order of magnitude of L can be estimated later after receiving equation (12) by comparison of the first two terms in this equation:
[9] The second limiting situation corresponds to neglecting collisions ni=0 and is well studied in publications. It should be mentioned that for ni=0, equation (6) takes the form
(8) |
(9) |
[10] Under Mi >M* 1.6, the solution in the form of a soliton is absent. This property is typical only for the case of a strong nonlinearity. At the presence of a weak nonlinearity, there is no such critical parameter Mi =M* (see, for example, stationary solution of the Korteweg-de Vries equations). At a transition from the value Mi =1 to the limiting value Mi =M*, there occurs an increase of the amplitude of soliton Y from zero up to Y* = Mi2/2. In this case a limiting value of the ion component velocity Vi = M* V s is reached.
[11] We take into account simultaneous influence of collisions ( ni 0 ) and dispersion ( De 0 ) under the condition Te Ti. To do that, we use the relation between Y and Vi:
(10) |
(11) |
[12] Using (10) and (11) we obtain the following formulae:
(12) |
[13] Taking into account equation (12), we obtain a description of the process in the problem of the Stage 1 for the dimensionless functions V=Vi/c
(13) |
[14] Equation (13) is of the third degree. To solve this equation, one has to formulate three boundary conditions. At analytical solution, one condition is put at the infinity behind the shock wavefront: V() =0 which corresponds to the limiting transition to the undisturbed condition. Two more conditions should provide a joining of the fields at the front. They may be obtained from the solution of the closed problem in the region behind the front (xi<0). Not doing that, we formulate the second condition approximately. Assuming an idealization of the complete carrying away of charges by the front a strong shock wave, we obtain the condition V(0) 2 (g+1)-1, g 1.4 being the ratio of the specific heat capacities. At this approximated approach, the formulation of the third condition presents some difficulties and its choice may be ambiguous. For example, the rule of selection of the solution by some sign may be such condition.
[15] Integrating equation (13) numerically (for example, by the Runge-Kutta method), one has to give all three boundary conditions in the initial point xi=0, but in such a way that at xi there occurred the transition V 0. Thus, two conditions on the derivatives V'(0) and V''(0) should simultaneously provide the limiting transition to the undisturbed state and the rule of the solution selection.
[16] It should be noted that at d=0 and A=0, (13) is transformed into an equation of the first order, but the solution of the physical problem should fulfill two boundary conditions: V(0) and V(). In the case Mi<1, there exists a continuous solution, whereas at Mi>1, the condition at the infinity V() may be fulfilled only due to introduction of a break [Pavlov, 1996]. The influence of the dispersion d 0 and viscosity A>0 leads to a smoothing of the break in the plasma precursor.
Figure 1 |
Figure 2 |
Figure 3 |
Figure 4 |
Figure 5 |
Figure 6 |
(14) |
and for the Mach ion numbers in the range Mi = 1.7 1.94. We note the most interesting properties of the plasma precursor of a shock wave. For the chosen version of the undisturbed plasma parameters in the Mach ion numbers Mi = 1.7 < 1.92, in the precursor there are two maxima of the ion concentration. They are two solitons. The soliton most distant from the point xi=0 provides a "smoothing" of the brake at the front of the ion-sonic precursor. Under fulfillment of the condition Mi = 1.7 < 1.85 (see Figures 1 and 2), the largest of the solitons is located closer to the front of the shock wave xi=0. Under Mi>1.86 (see Figures 3 and 4), the structure of the disturbance changes: the smaller of the solitons is located closer to the shock wavefront (the largest of the maxima moves to the front edge of the precursor). Then there happens a further sharp change in the structure: at Mi 1.92 there is formed one maximum (see Figure 5). Further on at Mi 1.94 (see Figure 6), a "destruction" of this plasma coagulum occurs and a ion concentration disturbance changing in a monotonous way is formed. The dimensionless spatial scale of the plasma precursor is of the order of the unity (this spatial scale has a value of the order of x1 nin V s ). An increase of the shock front velocity leads to a decrease of the dimension of the precursor (see Figures 1-5).
[18] Thus, the shock wave precursor has an oscillating structure (it can have several maxima). The precursor has nonlinear resonance properties relative to the motion velocity c. Depending on the plasma parameters in the undisturbed condition, the maximum values of the disturbance are reached at different values of the shock wave velocity in the range c (1.6 2)V s c*. The only maximally dense local coagulum of charged particles is formed in the precursor (in the considered version it occurs under c 1.93 V s, see Figure 5). If the shock wave exceeds the critical velocity c*, a collapse of this plasma coagulum occurs (see Figure 6). It is interesting that even the presence of a strong dissipation does not destroy the sharp resonance nonlinear effect. The nature of this effect is due to the rivaling of strong nonlinearity, strong dispersion, and dissipation.
[19] The gas in this coagulum is no more a weakly ionized one, and charged particles may influence on the neutral component and shock wave.
(15) |
(16) |
Powered by TeXWeb (Win32, v.2.0).