Next Article in Journal
Rheological Characterization of a Concentrated Phosphate Slurry
Previous Article in Journal
A Computational Fluid Dynamics Model for the Small-Scale Dynamics of Wave, Ice Floe and Interstitial Grease Ice Interaction
Previous Article in Special Issue
Gas–Liquid Two-Phase Flow and Heat Transfer without Phase Change in Microfluidic Heat Exchanger
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diffusion Mechanisms for the Occurrence of the Instability of Mechanical Equilibrium of a Ternary Gas Mixture Containing Carbon Dioxide

1
Department of Physics, Institute of Mathematics, Physics and Informatics, Abai Kazakh National Pedagogical University, 13 Dostyk Avenue, Almaty 050010, Kazakhstan
2
Institute of Experimental and Theoretical Physics, Al-Farabi Kazakh National University, 71 Al-Farabi Avenue, Almaty 050040, Kazakhstan
*
Authors to whom correspondence should be addressed.
Submission received: 22 March 2021 / Revised: 16 April 2021 / Accepted: 20 April 2021 / Published: 1 May 2021
(This article belongs to the Special Issue Flow and Heat Transfer Intensification in Chemical Engineering)

Abstract

:
Mixing of carbon dioxide dissolved in a multicomponent gas mixture at different pressures was researched. It was found that the mechanical equilibrium of the ternary gas mixture 0.4163H2 (1) + 0.5837CO2 (2) − N2 (3) is violated at a pressure of p = 0.7 MPa and structured flows appear in the system. The pressure area (from 0.7 to 1.5 MPa) at which the conditions of priority transfer of components with the highest molecular weight in the mixture are realised in the system is fixed. To analyse the effect of pressure on the process of changing “diffusion–convection” modes, a mathematical model, which takes into account the kinetic features of multicomponent mixing, was applied. It was shown that the change in the modes of mass transfer is associated with a significant difference in the diffusion ability of the components. It is noted that the difference in the diffusion coefficients of components results in the nonlinearity of the concentration distribution, which leads to the inversion of the density gradient of the gas mixture, which is the cause of convective flows.

1. Introduction

Separation of substances during evaporation of solutions from a free surface, separation of vapour-gas mixtures into components with specified properties is a common process in nature and in many practical applications. Despite many researches on the study of pressure, temperature, and concentration distributions of the main mixing elements in specified technological segments, these issues require further study, including those related to the occurrence of concentration gravitational flows caused by diffusion. It should be noted that convective motions in one-component media with a monotonic dependence of density on temperature are well studied and described in [1,2,3]. In these works, a variety of approaches related to the study of natural gravitational convection in channels of various shapes were generalised. These generalisations make it possible to develop methods for determining the main parameters of heat and mass transfer connected with the properties of the substance under study, temperature head (Rayleigh number), hydraulic losses, and other factors, which allow for specific conditions to predict the modes of mixing.
Researches carried out in [4,5] showed that in the absence of an interphase boundary in the system under study, the reason for the occurrence of thermo-concentration gravitational currents is the inhomogeneity of the density of the medium caused by variations in temperature or composition. Development of convective instability depends on the mutual orientation of the density gradient and the direction of gravity. If the direction of the density gradient coincides with the direction of the gravity force, in this case a stable stratification of the medium occurs. The occurrence of convective motion in this situation is impossible. However, the extension of this approach to the case of non-isothermic binary and multicomponent mixtures showed that diffusion could have a destabilising influence on the development of convective flows in mixtures with different diffusion rates of the components [6,7,8]. The arising movement is an important mechanism for intensifying heat and mass transfer and selective separation of components [8,9].
One of the features of the researches carried out in [4,6,7,8,9] was the assumption that heat and solute act as components influencing the overall density distribution of the medium. It should be assumed that the effects of separation associated with the difference in the diffusion coefficients of the components, pressure, and composition can appear in multicomponent mixtures during isothermal mixing. The formation of density inversion, which causes gravitational convection in the isothermal case of multicomponent mixing, is probably associated with the appearance of effects that are not typical for non-isothermal diffusion. Experimental and computational theoretical researches, given in [10,11,12], showed that convective flows possibly occur in the case when the condition for a decrease in density with height is realised in isothermal multicomponent gas mixtures. In addition, the predominant transfer of the component with the highest molecular weight was recorded in [13], which is not typical for diffusion. Therefore, the solution of issues related to the study and refinement of the separation mechanisms in multicomponent gas mixtures is important for the problems of combined mass transfer.
The aim of the work is to experimentally study the change of “diffusion–convection” modes in a three-component gas mixture containing carbon dioxide at different pressures and isothermal conditions, for the situation in which the density of the mixture decreases with height. In this work, a calculation was carried out to determine the boundary of the change in the “diffusion–convection” modes. The density distribution in a vertical cylindrical channel at the transition pressure was analysed. The obtained results of the numerical research were compared with the experimental data.

2. Experimental Research

Separation of a multicomponent gas mixture 0.4163H2 (1) + 0.5837CO2 (2) − N2 (3) was studied at different pressures. Numbers in front of the chemical element correspond to the initial concentration in mole fractions. The numbers in parentheses define the component numbering. Mixing of the components was carried out in a device representing a variation of the two-flask method, which is widely applied in experimental researches at various pressures and temperatures [14]. Flasks 2 and 3, connected by a vertical channel 1, represent a diffusion cell (Figure 1 [15]), in which multicomponent mixing is carried out under the given experimental conditions. The geometric characteristics of the diffusion cell are as follows: volumes of the upper and lower flasks are equal, and the amount of Vu = Vl = (55.0 ± 0.5) × 10−6 m3, the diameter and length of the vertical cylindrical channel are equal to d = (4.00 ± 0.05) × 10−3 m, L = (70.00 ± 0.05) × 10−3 m. The temperature in all experiments was constant and equal to T = 298.0 K. Mixing time for all experiments was 30 min. Experimental conditions assumed that the density of nitrogen in the lower flask 3 would be greater than the density of the mixture of hydrogen and carbon dioxide, which was located in flask 2.
The range of pressures under study is 0.2–1.7 MPa. The design features of the experimental stand determine the minimum pressure limit. The maximum pressure value is associated with the fact that at p ≥ 1.7 MPa due to a significant nonlinear dependence of the density of carbon dioxide on pressure [16] in the mixture under study, the sign of the density gradient changes, which does not correspond to the conditions of the problem. In this regard, the analysis of multicomponent mixing at pressures exceeding 1.7 MPa was not carried out. The concentration of a specified component was measured by the gas chromatography method by comparing the chromatograms by the height of the peaks of the components before and after mixing at a given pressure and temperature. For the mixture under study, the error in determining the component concentration was found by the estimation method of indirect measurement errors and was ±1%.
Experimental methods involved the implementation of certain procedures described, for example, in [17]. After evacuating the upper 2 and lower 3 flasks with the channel 1 connecting them with the valve 4 closed, the gases under study were fed into the diffusion cell. The upper flask 2 should be filled with a gas mixture of hydrogen and carbon dioxide from cylinder 5 through the opening of valves 9 and 11. Valve 13 must be closed. Pressurisation of the upper flask 2 was carried out up to the pressure of the experiment and was controlled by the manometer 7. After reaching the set pressure, valves 9 and 11 were closed. Similarly, the lower flask 3 was filled with nitrogen from the cylinder 6 through the open valves 10, 12 with the closed valve 14. The manometer 8 recorded the pressure of the experiment. Upon completion of the filling procedure, valves 10 and 12 were closed. Mixing of gases occurs when valve 4 is opened with synchronous recording of the mixing time. Valve 4 was closed at the end of the experiment. Then, the mixtures under study were fed through the outlets 15 to the chromatograph for analysis.
In the ternary mixture, the concentrations of carbon dioxide and hydrogen were directly measured. Concentration of nitrogen can be obtained from the condition of conservation number of particles i = 1 n c i = 1 , where ci is the concentration of i-th component. The experiment was repeated several times, and through averaging of the measured values, the specific value of the concentration of the components was determined. One of the main points of the experimental study is the registration of the boundary parameters that determine the transition from one state to another. Obviously, this can be done by comparing experimental data (for example, the concentrations of components) with those calculated by the Stefan-Maxwell equations [18] under the assumption of diffusion in ternary mixtures. If the experimental concentrations cexp and the theoretical values ctheor calculated according to the Stefan-Maxwell equations are close to each other, this indicates a diffusion-mixing mode. The significant difference between the concentrations in experiment and theory shows the presence of other mixing mechanisms. The transition boundary can be fixed if we trace the ratio of the diffused concentrations of the components α = cexp/ctheor. Figure 2 shows the functional dependence of the dimensionless parameter α on pressure for the analysed system.
Analysis of the data, which is depicted in Figure 2, shows that diffusion is realised in the system in the pressure area of 0.1–0.6 MPa. This is evidenced by the equality to unity of the dimensionless parameter αi for all three components of the mixture under research. Starting from a pressure of 0.6 MPa, an increase in αi for all components in the system is fixed. Unusually significant increase was noted for CO2 and to a lesser extent for the most diffusionally mobile H2. Apparently, the pressure p* ≥ 0.6 MPa provides a change in the transfer mode from diffusion to convective. Further increase of pressure leads to an increase in the intensity of convective currents. This mixing is not typical for diffusion, in which the mixing intensity decreases with increasing pressure. The highest values of αi for CO2 at pressures exceeding the transition pressure indicate that under the given conditions there is a predominant transfer of CO2. With a further increase in pressure, the linear dependence is violated and the mixing intensity of the component behaves in a non-monotonic manner.

3. Determination of the Boundary of the Change of Diffusion-Convection Modes in Isothermal Ternary Gas Mixtures

Generalised description of the combined mass transfer for multicomponent mixtures seems to be a difficult problem. However, to determine the boundary of the change of “diffusion–convection” modes, some assumptions are possible. These assumptions allow performing numerical calculations to determine the characteristic features of partial mass transfer. Parameters determining the kinetic transition from diffusion to convective mixing can be obtained within the framework of stability theory [4,5] for the case of an isothermal ternary gas mixture.
In this case, the condition of independent diffusion is taken into account:
i = 1 3 j i = 0 , i = 1 3 c i = 1 .
The system of equations describing convective mixing in three-component mixtures consists of the Navier-Stokes equation, the convective diffusion equation, and the equation for the conservation of the number of particles of the components [12]:
ρ [ u t + ( u u ) ] = p + η 2 u + ( η 3 + ξ ) div u + ρ g n t + div ( n v ) = 0 , c i t + v c i = div j i , j 1 = ( D 11 * c 1 + D 12 * c 2 ) , j 2 = ( D 21 * c 1 + D 22 * c 2 ) ,
where j i is the density of the diffusion flux of the i-th component, ci is the concentration of the i-th component, p is the pressure, g is the acceleration of gravity, u = ρ 1 u 1 + ρ 2 u 2 + ρ 3 u 3 ρ is the mass average velocity, v = n 1 u 1 + n 2 u 2 + n 3 u 3 n is the number average velocity of the ternary mixture, ρ is the density, and η and ξ are the coefficients of shear and bulk viscosity, respectively. Practical diffusion coefficients D i j * are determined as follows:
D 11 * = D 13 [ c 1 D 32 + ( c 2 + c 3 ) D 12 ] D , D 12 * = c 1 D 23 ( D 12 D 13 ) D , D 22 * = D 23 [ c 2 D 13 + ( c 1 + c 3 ) D 12 ] D , D 21 * = c 2 D 13 ( D 12 D 23 ) D , D = c 1 D 23 + c 2 D 13 + c 3 D 12 ,
where Dij are the interdiffusion coefficients calculated within the framework of traditional kinetic concepts.
Equation (1) are completed by the equation of the medium state,
ρ = ρ ( c 1 , c 2 , p ) , T = const .
To solve the system of Equations (1) and (3), we apply the method of small perturbations [4], which allows presenting ci and p in the following way:
c i = c i + c i / , p = p + p / ,
where c i , p are the constant average values taken as the point of origin, ci/, p/ are the perturbed values which do not lead to a significant deviation of the density ρ/ from the average value ρ0. Furthermore, we assume that under the condition L d / 2 the differences between the perturbations of the number average v and mass average u velocities in the Navier-Stokes equation will be insignificant. Then, Equations (1) and (3) can be reduced to a system of equations for the perturbations undimensionalised relative to the given scales: distance d, time d2/ν, velocity D*22/d, concentration of the i-th component Aid, and pressure ρ0νD*22/d2. Finally, the system of perturbed dimensionless equations (the strokes are omitted) takes the form:
P r 22 c 1 t ( u γ ) = τ 11 2 c 1 + A 2 A 1 τ 12 2 c 2 , P r 22 c 2 t ( u γ ) = A 1 A 2 τ 21 2 c 1 + 2 c 2 , u t = p + 2 u + ( R a 1 τ 11 c 1 + R a 2 c 2 ) γ , div u = 0 ,
where P r i i = ν / D i i * is the diffusion Prandtl number, R a i = g β i A i d 4 / ν D i i * is the partial Rayleigh number, ν is the kinematic viscosity coefficient, β i = 1 ρ 0 ( ρ c i ) p , T , τ i j = D i j * / D 22 * are the parameters determining the relationship between the practical diffusion coefficients, γ is the unit vector directed along the upward vertical. Parameters τ i j = D i j * / D 22 * are calculated through D i j * defined by Formula (2). At p = 0.101 MPa and T = 273.0 K, the interdiffusion coefficients of gas vapour have the following values: D H 2 N 2 = 6.883 × 10−5 m2·s−1, D H 2 C O 2 = 5.75 × 10−5 m2·s−1, D N 2 C O 2 = 1.44 × 10−5 m2·s−1 [16]. The calculations took into account the change of IDC (interdiffusion coefficient) as a function of pressure and temperature.
The solution of the system of Equation (4) can be conditionally divided into several stages [19]. At the first stage, taking into account the existence of three-dimensional motion in a cylindrical channel of finite length, an approximation of the velocity is written, considering all components of the vector u to be nonzero. Furthermore, on the basis of the continuity equation, the ratio connecting the radial constituents of the velocity components is written. This ratio is determined through the n-order Bessel functions. At the second stage, from the first two equations of System (4), the concentrations of the components are determined by means of the Kantorovich method under boundary conditions that suggest the vanishing of perturbations of the velocity and flow of matter at the boundaries of the diffusion channel. At the third stage, the boundary of monotonic stability is determined.
On the basis of the described approach to solving the system of equations of concentration convection (4) for a vertical cylindrical channel of finite height with boundary conditions suggesting zero perturbations of the velocity and flow of matter at the boundaries, in [19] the boundary relation in terms of the Rayleigh numbers Rai, defining the change of “diffusion–convection” modes, was obtained in the form:
u 2 u d V + R a 1 τ 11 u z c 1 d V + R a 2 u z c 2 d V = 0 .
Based on the obtained Solution (5), it is possible to determine a stability map in coordinates (Ra1, Ra2) and display the experimental data on it. Experimental data are restored through the partial Rayleigh numbers, for example, according to the formulas given in [19],
R a 1 = g n r 4 Δ m 1 ρ ν D 11 * c 1 z , R a 2 = g n r 4 Δ m 2 ρ ν D 22 * c 2 z ,
where Δ m 1 = m 1 m 2 , Δ m 2 = m 2 m 3 , mi is the molecular mass of i-th component, n = p / k B T , and then represented on the coordinate plane (Ra1, Ra2). Zero density gradient line, the position of which is determined by the following equation
τ 11 R a 1 = R a 2 ,
is also drawn on the cartogram (Ra1, Ra2).
Crossing of the boundary line defined by Equation (5) with axes Ra1 and Ra2 allows obtaining two critical Rayleigh numbers. Results of a numerical study of the system 0.4163H2 + 0.5837CO2 − N2 are shown in Figure 3. It is known from the experiment which mixing mode is realised under given conditions. Let us agree that for diffusion we use the notation in the form of lightened points «o». The mode, where mixing is determined by the convective perturbations, is denoted by the dark symbol «●». A collection of points on the plane (Ra1, Ra2) determines the correspondence of the observed mixing modes to the calculated boundary relations. Critical Rayleigh numbers, which determine the stability boundary, have the following values: Ra1 = 94.4773 and Ra2 = 61.6623. In a range of pressures of 0.2 to 0.6 MPa, the system under consideration is in the diffusion area, which is in agreement with the experimental data shown in Figure 2. Numerical calculations showed that the transition to the convective area occurs at a pressure of 0.775 MPa (point 4). Based on the presented comparison, it can be concluded that concentration convection in the system under study can occur under conditions in which the density of the mixture decreases with height.
The possibility of the system transition from the area, where diffusion mechanisms operate to the area, where the conditions of mechanical equilibrium are violated is associated with the existence of density inversion in the vertical channel due to a significant difference in the diffusion coefficients of the components.
Density distribution in the diffusion channel plays one of the important roles in creating conditions for the gravitational convection, therefore, it is necessary to analyse the behaviour of the total density of the gas mixture and concentrations of the components. For this, it is necessary to solve the system of equations of stable three-component diffusion for a closed system of two flasks connected by a vertical channel:
p = const , t = const , i = 1 3 c i = 1 ; n i = 1 3 c i u i = 0 , div ( n c i u i ) = 0 ,   i = 1 ,   2 ; j i c i c j D i j ( u i u j ) = grad ( c i ) ,   i = 1 ,   2 ;   j = 1 ,   2 ,   3 .
Solution of the system of Equation (7) numerically makes it possible to obtain the distribution of component concentrations along the diffusion channel and the distribution of the mixture density. Figure 4 shows the distribution of the density of the mixture under consideration and the concentration of components along the length of the channel at a pressure of 0.775 MPa, which determines the transition of the system from one kinetic regime to another one. The value of the extremum coordinate in the channel corresponds to zex = 0.0478 m.
According to Figure 4a, for the heaviest in density component of the mixture (carbon dioxide), a nonlinear concentration distribution is observed along the diffusion channel. The appearance of nonlinearity in the concentration profile of the heavy component is due to the fact that a significant difference in the interdiffusion coefficients ( D 12 , D 13 D 23 ) leads to a diffusion “lock” of carbon dioxide. The diffusion “lock” is conditioned by the fact that a more mobile gas (hydrogen) diffusing into the pure component (nitrogen) causes the appearance of a diffusion baroeffect. In this case, the diffusion baroeffect generates a hydrodynamic flow. The intensity of this flow is such that the counter diffusion flow of the heavy component (carbon dioxide) cannot overcome this flow, i.e., a decrease in the flux of the heavy component is observed. This state of the component with the highest molecular weight continues until the concentration of the components in the flasks is redistributed to such values that the magnitude of the baroeffect cannot generate a hydrodynamic flow of the required intensity to keep the heavy component in a locked state. Thus, a significant difference in the interdiffusion coefficients of the components creates conditions for a nonlinear distribution of the concentration of the heaviest component in terms of density. Nonlinearity of the concentration distribution of the component with the highest molecular weight results in the inversion of the density gradient of the gas mixture (Figure 4b), which is the reason for the occurrence of convective flows in the gas mixture under consideration.

4. Conclusions

Experimental data obtained for the system 0.4163H2 (1) + 0.5837CO2 (2) − N2 (3) indicate that in the system under consideration at a pressure of 0.6 MPa, a transition from a stable diffusion process to an unstable one is observed. The increase of pressure in the experiment creates conditions under which the diffusion tendency of the process is disturbed and convection occurs. The behaviour analysis of the dimensionless parameter α on pressure showed that in the pressure range from 0.6 to 1.5 MPa, the largest transfer for the heavy component (carbon dioxide) is observed. This behaviour of the concentration of carbon dioxide on pressure indicates an increase in the intensity of convective currents. The performed stability analysis made it possible to determine the pressure area (p = 0.775 MPa), where the conditions for the occurrence of convective flows and the priority transfer of carbon dioxide are realised. The obtained distributions of component concentrations and the average density of mixture for the system under consideration have shown a significant role of the inversion of density gradient in the occurrence of convective flows. For the mixture under consideration and the given geometry of the diffusion channel, the inversion of the density gradient is observed at the extreme coordinate zex = 0.0478 m. The discovered mechanisms of the predominant transfer of carbon dioxide in gas mixtures can be used in the technological processes associated with the removal of greenhouse gases. Further experimental studies will be aimed at establishing the boundaries of the transition from diffusion to convection in non-isothermal conditions, in vertical and inclined channels of various geometries. Theoretical research involves the study of developed convective mixing and the determination of the boundaries of stability of various modes of motion.

Author Contributions

Conceptualization, V.K. and O.F.; methodology, V.K. and A.K.; software, O.F. and A.Z.; validation, V.K. and A.K.; formal analysis, A.Z.; investigation, A.K.; resources, A.Z.; data curation, A.K.; writing—original draft preparation, V.K. and O.F.; writing—review and editing, V.K. and O.F.; visualization, A.K.; supervision, V.K. and O.F.; project administration, O.F.; funding acquisition, O.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Science Committee of the Ministry of Education and Science of the Republic of Kazakhstan, grant number AP08955418.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bejan, A. Convection Heat Transfer; Wiley: New York, NY, USA, 2004. [Google Scholar]
  2. Jaluria, Y.; Torrance, K.E. Computational Heat Transfer, 2nd ed.; CRC Press: New York, NY, USA, 2002. [Google Scholar]
  3. Ho, C.; Webb, S. Gas Transport in Porous Media; Springer: New York, NY, USA, 2006. [Google Scholar]
  4. Gershuni, G.Z.; Zhukhovitskii, E.M. Convective Stability of Incompressible Fluids; Keter Publishing House: Jerusalem, Israel, 1976. [Google Scholar]
  5. Nield, D.A.; Bejan, A. Convection in Porous Media; Springer: New York, NY, USA, 2006. [Google Scholar]
  6. Matsuura, H.; Nagasaka, Y. Soret Forced Rayleigh Scattering Instrument for Simultaneous Detection of Two-Wavelength Signals to Measure Soret Coefficient and Thermodiffusion Coefficient in Ternary Mixtures. Rev. Sci. Instrum. 2018, 89, 024903. [Google Scholar] [CrossRef] [PubMed]
  7. Hu, Y.-P.; Li, Y.-R.; Li, M.-H.; Zhang, L.; Li, S.-Z. Effects of Enclosure Geometry and Thermal Boundary Condition on Rayleigh-Benard Convection of Cold Water near Its Maximum Density. Int. J. Therm. Sci. 2017, 120, 220–232. [Google Scholar] [CrossRef]
  8. Lyubimova, T.P.; Sadilov, E.S.; Prokopev, S.A. Onset of Soret-induced Convection in a Horizontal Layer of Ternary Fluid with Fixed Vertical Heat Flux at the Boundaries. Eur. Phys. J. E 2017, 40, 15. [Google Scholar] [CrossRef] [PubMed]
  9. Shevtsova, V.; Santos, C.; Sechenyh, V.; Legros, J.C.; Mialdun, A. Diffusion and Soret in Ternary Mixtures. Preparation of the DCMIX2 Experiment on the ISS. Microgravity Sci. Technol. 2014, 25, 275–283. [Google Scholar] [CrossRef]
  10. Dil’man, V.V.; Lipatov, D.A.; Lotkhov, V.A.; Kaminskii, V.A. Instability in Unsteady-state Evaporation of Binary Solutions into an Inert Gas. Theor. Found. Chem. Eng. 2005, 39, 566–572. [Google Scholar] [CrossRef]
  11. Dil’man, V.V.; Lotkhov, V.A. Molecular Turbulent Evaporation in a Gravitational Field. Theor. Found. Chem. Eng. 2015, 49, 102–106. [Google Scholar] [CrossRef]
  12. Kossov, V.; Fedorenko, O.; Zhakebayev, D. Features of Multicomponent Mass Transfer in Gas Mixtures Containing Hydrocarbon Components. Chem. Eng. Technol. 2019, 42, 896–902. [Google Scholar] [CrossRef]
  13. Kosov, V.N.; Seleznev, V.D.; Zhavrin, Y.I. Separation of Components during Isothermal Mixing of Ternary Gas Systems under Free-Convection Conditions. Tech. Phys. 1997, 42, 1236–1237. [Google Scholar] [CrossRef]
  14. Nezovitina, N.A.; Bogatyrev, A.F.; Makeenkova, O.A. Dependence of the Thermodiffusion Constant of a Mixture of Two Gases on the Addition of a Third Gas to It. J. Eng. Phys. Thermophys. 2016, 89, 733–740. [Google Scholar] [CrossRef]
  15. Kossov, V.N.; Fedorenko, O.V.; Krasikov, S.A.; Kalimov, A.B. Separation Features of Multicomponent Gas Mixtures at the Border of Regime Change Diffusion—Concentration Gravitational Convection. J. Phys. Conf. Ser. 2020, 1683, 022045. [Google Scholar] [CrossRef]
  16. Vargaftik, N.B. Handbook of Physical Properties of Liquids and Gases. Pure Substances and Mixtures; Springer: Berlin/Heidelberg, Germany; GmbH & Co. KG.: Berlin, Germany, 2014. [Google Scholar]
  17. Asembaeva, M.K.; Kosov, V.N.; Krasikov, S.A.; Fedorenko, O.V. Effect of the Channel Inclination Angle on Convective Mixing Caused by Instability of Mechanical Equilibrium of Ternary Gas Mixture at Isothermal Diffusion. Tech. Phys. Lett. 2019, 45, 1071–1074. [Google Scholar] [CrossRef]
  18. Bird, R.B.; Stewart, W.E.; Lightfoot, E.N. Transport Phenomena, 2nd ed.; John Wiley & Sons: New York, NY, USA, 2007. [Google Scholar]
  19. Kossov, V.; Krasikov, S.; Fedorenko, O. Diffusion and Convective Instability in Multicomponent Gas Mixtures at Different Pressures. Eur. Phys. J. Spec. Top. 2017, 226, 1177–1187. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the experimental stand [15].
Figure 1. Schematic diagram of the experimental stand [15].
Fluids 06 00177 g001
Figure 2. Dependence of α on pressure. The line is the calculation with the assumption of diffusion. The points correspond to: Δ—hydrogen; ○—carbon dioxide.
Figure 2. Dependence of α on pressure. The line is the calculation with the assumption of diffusion. The points correspond to: Δ—hydrogen; ○—carbon dioxide.
Fluids 06 00177 g002
Figure 3. Map of stability for the system 0.4163H2 + 0.5837CO2 − N2: I is the neutral line of monotonic perturbations; II is the line of zero density gradient; 1–7 are the numerical data at different values of p: 1—0.2; 2—0.4; 3—0.6; 4—0.7; 5—0.8; 6—1.0; 7—1.2 MPa.
Figure 3. Map of stability for the system 0.4163H2 + 0.5837CO2 − N2: I is the neutral line of monotonic perturbations; II is the line of zero density gradient; 1–7 are the numerical data at different values of p: 1—0.2; 2—0.4; 3—0.6; 4—0.7; 5—0.8; 6—1.0; 7—1.2 MPa.
Fluids 06 00177 g003
Figure 4. Changes in the concentrations of components (a) and the total density (b) along the length of the channel for the system 0.4163H2 + 0.5837CO2 − N2: 1—H2; 2—N2; 3—CO2.
Figure 4. Changes in the concentrations of components (a) and the total density (b) along the length of the channel for the system 0.4163H2 + 0.5837CO2 − N2: 1—H2; 2—N2; 3—CO2.
Fluids 06 00177 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kossov, V.; Fedorenko, O.; Kalimov, A.; Zhussanbayeva, A. Diffusion Mechanisms for the Occurrence of the Instability of Mechanical Equilibrium of a Ternary Gas Mixture Containing Carbon Dioxide. Fluids 2021, 6, 177. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids6050177

AMA Style

Kossov V, Fedorenko O, Kalimov A, Zhussanbayeva A. Diffusion Mechanisms for the Occurrence of the Instability of Mechanical Equilibrium of a Ternary Gas Mixture Containing Carbon Dioxide. Fluids. 2021; 6(5):177. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids6050177

Chicago/Turabian Style

Kossov, Vladimir, Olga Fedorenko, Adilet Kalimov, and Aiym Zhussanbayeva. 2021. "Diffusion Mechanisms for the Occurrence of the Instability of Mechanical Equilibrium of a Ternary Gas Mixture Containing Carbon Dioxide" Fluids 6, no. 5: 177. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids6050177

Article Metrics

Back to TopTop