Within the linear theory of stability, the process of isothermal mixing of three-component gas mixtures in a channel of final dimensions in the absence of mass-transfer through its walls is considered. The comparison of experimental data with the results of theoretical calculations for the mixtures He+Ar-N2 and H2+N2-CH4 is shown, that a stable diffusion process as the temperature increased will remain the same and be described by the ordinary diffusion laws, but unstable one lost its intensity and tend to the stable diffusion.
The study of a diffusion process in multicomponent gas mixtures showed that mutual influence of components to each other under certain conditions lead to phenomena no taking place in ordinary binary diffusion. One of the most interesting peculiarities both in scientific and practical application is diffusion instability [
The determination of a transition boundary of mixture from the stable state to the unstable one and back is one of the main factors in the study of regime change of mixing. In the paper [
Experimental data reported in [2-6] in their physical meaning resemble the problem of instability which arises under conditions of thermal convection [8-10]. Application of stability theory [
The aim of this study is to examine the transition from the state of diffusion to the regime of concentration gravitational convection (diffusion instability) in a channel of finite size in the absence of mass-transfer through its wall. In addition, the obtained data are compared with the experiments presented in [
The macroscopic flow of the isothermal ternary gas mixture is described by the general system of the hydrodynamic equations, that includes the Navier-Stokes equations, equations for conservation of the number of particles in the mixture and the components. Taking into account the conditions of independent diffusion, for which the and are valid, the system of equations takes the following form [
where “practical” coefficients of three-component diffusion are defined in accordance with expressions:
where are the mutual diffusion coefficients indicating interference of two components.
The system of Equation (1) is supplemented with the environmental state equation
interrelating the thermodynamic parameters entering into the system of Equation (1).
The method of small perturbations [
where is the Prandtl diffusion number, is the Rayleigh partial number, denotes the parameters, which determine the relationship between the “practical” diffusion coefficients.
Three-dimensional motions are essential in the diffusion cylindrical channel of final length. Therefore, when approximating the velocity, it is necessary to consider all components of the vector differ from zero. Examining the periodical motions along φ in the cylindrical coordinate and satisfying the conditions on the hard boundaries (is the geometry parameter characterizing stability), the velocity approximation can be written as [
The radial functions u, υ, ω should be vanish on the hard lateral surface of cylinder (at). The relation binding these functions follows from the continuity equation:
though
where is n-order Bessel function, but the parameter k can be found from the equation:
Assuming for the first two equations in the set of Equations (2), that, we determine the concentration of components from the following equations:
where
We will consider, that, then
Supposing that the perturbations of concentrations vanish on the end surface, then the additional condition implying that the second derivative vanishes at the ends of diffusion channel results from Equation (6):
that permits to select the following approximation:
whereis the radial function of concentration.
For determination we use Kantorovich method. Substituting Equation (8) into Equation (6), multiplying by depending on z part of the function and integrating between – h and h we obtain the equation:
where Finding components concentrations of Equation (9) were solved with the boundary condition, and then the final solution in the center becomes (see Equation (10)):
where Jn and In are n-order Bessel functions of the first kind.
In order to determine the monotonous stability boundary of the problem under consideration, the third equation of the system (2) can be scalarly multiplied by the vector and integrated all over the volume V of the diffusion channel. This can be done under the conditionsthat. Then we have:
This equation in the coordinates gives a straight line MM, dividing the region of molecular transport and the region of the diffusion instability. Figures 1 and 2 show the location of the neutral line of monotonic instability for the systems He+Ar-N2 and H2+N2-CH4 for. The region that lies below the line MM corresponds to diffusion.
From the condition of zero density gradient of the mixture and with allowance for the determined values of partial Rayleigh numbers (2), we obtain the following equation for the line in the plane:
The mutual position of the line of monotonic instability MM (Equation (11)) and the line (Equation (12)) for is shown in Figures 1 and 2. As follows from figures, there exists a region in the plane where the line MM is situated below the line (12). In this region, the mixture appears to be unstable.
To compare the approach proposed in Paragraph 2 with the experimental data is shown in [
where is the molecular mass of the i-th component,. If conditions of the experiment are known (pressure, temperature, composition of mixtures in each of the flasks, the size of the diffusion channel), then according to Equation (13) we can find R1 and R2 and thus determine the point representing this experiment on the plane. From experiment, we know what the regime (diffusion or convection) occurs under the given conditions. In addition, we assume that the number 1 refers to the component with the lowest density, while the numbers 2 and 3 refer to the components with the highest and intermediate densities respecttively.
Figures 1 and 2 show the experimental data in terms of the Rayleigh numbers for the systems 0.4688Не(1)+0.5312Ar(2)-N2(3) and 0.6179Н2(1)+0.3821N2(2)-CH4(3) and for the systems 0.6192Не(1)+0.3808Ar(2)-N2(3) and 0.6463Н2(1)+0.3537N2(2)-CH4(3) respectively, obtained by varying the temperature. The full circles correspond to the convective mixing process while the open circles conform to the diffusion one.
On Figures 1, all points displaying the experiment are situated in the area of instability (region between the line stability MM and the line of zero density gradient), but with temperature rise tend to the region of stable diffusion (region below the line MM).
As is seen, in
Thus, the results reveal that the unstable diffusion process tends to stable one by raising the temperature, but a decrease in the temperature can cause transition to the unstable state of a diffusion-stable system.
In this paper, we have studied the instability of mechanical equilibrium (diffusion instability) in ternary gas mixtures at different temperatures. Linear stability theory is applied to determine the boundary of transition from the state of diffusion to the region of the concentration convection. The problem has been solved for a cylindrical channel of finite size and boundary conditions, suggesting the absence of component transfer through the walls of the diffusion channel. The comparison of theoretical results with the experimental data for the study of the
temperature dependence of the diffusion mixing of ternary gas mixtures He+Ar-N2 and H2+N2-CH4 indicates qualitative and quantitative agreement. Therefore, it can be concluded that the proposed calculation technique permits estimation of the critical parameters of the transition to the unstable regime.
This work has been performed under the support of the Fund of Science of the Ministry of Education and Science of Republic of Kazakhstan (the project number: No. 1134/GF).