Analysis of temperature jump and partial pressure coefficients of a binary evaporating gas mixture

The Maxwell-Loyalka method is used to derive expressions for the boundary temperature and partial pressure jumps of a binary mixture of gases evaporated (condensed) on a wall. The evaporation (condensation) is assumed to be weak. The expressions for the jumps are written in terms of the components normal to the wall of the reduced heat flux of the gas mixture and the average velocities of its components, which makes possible direct generalization to a multicomponent mixture, invariance of the coefficients of the expressions for the jumps with respect to the number of evaporating components, and symmetry of the cross effects. These coefficients are analyzed in the same way as in [1], and it is shown that they can be considerably simplified with an accuracy acceptable for practical purposes.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 150–159, November–December, 1991.     

Transformations of the Burnett components of transport relations in gases

The well-known expressions for the Burnett components of the stress tensor and the heat flux vector in a monatomic gas are transformed to a more effective form. The greatest simplifications are obtained for the stress tensor when using the usually employed approximate values of the Burnett transport coefficients, which are exact for the gas consisting of Maxwellian molecules.     

Requirements to the accuracy of the Burnett transport coefficients

Some effects in gases are completely or partially determined by the Burnett terms of the stress tensor and the heat flux vector. Usually, approximate values of the Burnett transport coefficients for a monatomic gas are used to calculate the above-mentioned quantities. “Exact” expressions for these coefficients correct to values less than one percent are obtained for elastic-sphere molecules. The results of calculations of certain such effects are compared with one another using the approximate and exact Burnett transport coefficients for a monatomic gas consisting of elastic-sphere molecules.     

Calculation of a turbulent two-phase isobaric jet

A numerical calculation is carried out by the finite-difference method based on proposed equations for a turbulent submerged jet containing an admixture of solid particles. The relative longitudinal particle velocity and the influence of particles on the turbulence intensity are taken into account. The calculated results adequately agree with available experimental data. A turbulent two-phase jet is examined in [1] on the basis of the theory for a variable density jet, assuming equal mean velocities for the gas and particles and not considering the influence of particles on the turbulence intensity. Particles are analogously taken into account by a noninertial gas mixture in [2, 3], and a particle Schmidt number of 1.1 is assumed in [4]. A model is proposed in [5] which takes into account the influence of particles on the turbulence intensity of the gas phase. Problems concerning the initial and main sections of a submerged jet were solved in [6] by the integral method on the basis of this model and the assumed equality of the mean velocities of the gas and particles. Turbulent mixing of homogeneous two-phase flows with allowance made for dynamic nonequilibrium of the phases is considered in [7]. However, the neglect of turbulent transfer of particle mass and momentum led to a physically unrealistic solution for the particle concentration in the far field of the mixture. A two-phase jet is considered in the present work on the basis of the theory of a two-velocity continuous medium [8, 9] with allowance made for turbulent transfer of particle mass and momentum. The influence of particles on the turbulence intensity of the gas phase is taken into account with the model of [5].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 57–63, September–October, 1976.The author acknowledges useful comments and discussion.of the work by G. N. Abramovich and participants of his seminar. The author sincerely thanks I. N. Murzinov for scientific supervision of the work.     

Transport properties of nonequilibrium polyatomic gas mixtures

It is shown that the well-known Hirschfelder-Euken correction to the thermal conductivity of a polyatomic gas mixture given by the first approximation in Sonine polynomials can be less than the corresponding exact value (for a Lorentz mixture of light and heavy molecules interacting in accordance with Coulomb's law) by a factor of 3.4. Fairly high accuracy is achieved in the second approximation in Sonine polynomials. Within the framework of the latter, similar corrections to the nonequilibrium heat and diffusion fluxes are found. On the basis of the generalized Chapman-Enskog method a more general case is studied. In this case some of the nonelastic collision integrals is also taken into account in calculating the transport coefficients. The transport coefficients are either represented in terms of the well-known formulas for fast and retarded internal molecular energy exchange or convenient approximate expressions are obtained.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 183–189, March–April, 1995.     

Laminar boundary layer in flow of a nonequilibrium ionized gas

The motion of a hypersonic body is accompanied by an increase in the gas temperature in the boundary layer up to tens of thousands of degrees, which causes the gas to ionize. Under these conditions there are problems in calculating coefficients of viscosity, diffusion, and heat conduction. Investigations have shown that it is invalid to extrapolate the widely used approximations for transport coefficients in the high temperature region [1–3]. This paper considers the laminar boundary layer in the vicinity of the stagnation point of a blunt body in a stream of monatomic nonequilibrium ionized gas. The main thrust is a more accurate calculation of transport coefficients and an investigation of their effect on profiles of the gasdynamic parameters. A specific calculation is performed for argon by way of example.     

Modification of the first approximation of the Chapman-Enskog method for a gas mixture

The authors propose a transformation of the equations of the first approximation of the Chapman-Enskog method for a gas mixture with frozen internal degrees of freedom. As a result, the solution for the perturbation of the distribution functions can be written in terms of the diffusion velocities and temperature gradient, and the derivation of the Stefan-Maxwell relations and the heat flux calculations can be much simplified. The transformations are extended to a mixture of polyatomic gases with nonequilibrium excitation of the internal degrees of freedom of the molecules. The modification of the first approximation of the Chapman-Enskog method does not affect the relations used for calculating the viscosity of the mixture. Accordingly, that part of the solution is not considered.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 178–185, July–August, 1992.     

Determination of heat fluxes in the neighborhood of a double curvature stagnation point in a flow of dissociating gas of arbitrary chemical composition

Numerical solutions are obtained for the equations of a uniform compressible boundary layer with variable physical properties in the vicinity of a stagnation point with different principal curvatures in the presence of an injected gas with the same properties as the incident flow. The results of the numerical solutions are approximated for the heat flux in the form of a relation that depends on the variation of the product of viscosity and density across the boundary layer and on the ratio of the principal radii of curvature.Using the concepts of effective diffusion coefficients in a multicomponent boundary layer, previously introduced by the author in [1], and the generalized analogy between heat and mass transfer in the presence of injection, together with the numerical solutions obtained, it is always possible, even without additional solutions of the boundary-layer equations, to derive final formulas for the heat fluxes in a flow of dissociating gas of arbitrary chemical composition, provided that we make the fundamental assumption that all recombination reactions take place at the surface.By way of example, formulas are given for the heat transfer to the surface of a body from dissociating air, regarded as a five-component mixture of the gases O, N, NO, O₂, N₂, and from a dissociating mixture of carbon dioxide and molecular nitrogen of arbitrary composition, regarded as an eleven-component mixture of the gases O, N, C, NO, C₂, O₂, N₂, CO, CN, C₃, CO₂.In the process of obtaining and analyzing these solutions it was found that, in computing the heat flux, a multicomponent mixture can be replaced with an effective binary mixture with a single diffusion coefficient only when the former can be divided into two groups of components with different (but similar) diffusion properties. In this case the concentrations of one group at the surface must be zero, while the diffusion flows of the second group at the surface are expressible, using the laws of mass conservation of the chemical elements, in terms of the diffusion flows of the first. Then the single effective diffusion coefficient is the binary diffusion coefficient D(A,M), where A relates to one group of components and M to the other.In view of the small amount of NO(c(NO) < 0.05), the diffusion transport of energy in dissociated air maybe described with the aid of a single binary diffusion coefficient D(A, M)(A=O, N, M=O₂, N₂, NO). However even in the case of complete dissociation into O and C atoms at the outer edge of the boundary layer, the diffusion transport of energy in dissociated carbon dioxide can not be described accurately enough by means of a model of a binary mixture with a single diffusion coefficient, since the diffusion properties of the O and C atoms are distinctly different.     

Calculation of convective heat flux in the vicinity of the stagnation point for partially ionized gas flow past a body

From numerical solutions of the boundary layer equations for a four-component gas mixture (E, N⁺, N₂, and N) with gas injection, approximate formulas for the heat flux as a function of the variation of λρ/c_p and h^* across the boundary layer and the magnitude of the objection are obtained (λ is the thermal conductivity of the mixture,ρ is density, c_p is the specific heat, and h^* is the enthalpy of the ideal gas state of the mixture). An effective ambipolar diffusion coefficient D^(a)(i) is introduced, making possible finite formulas for the convective heat fluxes in the “frozen” boundary layer. We study the behavior of these coefficients within the boundary layer. A formula is obtained for convective heat flux to the wall from partially ionized air for a nine-component mixture (E, O⁺, N⁺, NO⁺, O, N, NO, O₂ N₂). Even for simpler four-component gas model three effective ambipolar diffusion coefficients are necessary: $$\begin{gathered} D^{(a)} (A) = D (A, M) D^{(a)} (I) = 2D (A, M), \hfill \\ D^{(a)} (M) = [ 1 + c_e (I)] D(A, M). \hfill \\ \end{gathered} $$ Here D(A, M) is the binary diffusion coefficient of the atoms into molecules, and c_e(I) is the ion concentration at the outer edge of the boundary layer. The assumption of an infinitely large charge-exchange cross section and the other simplifying assumptions used in [1] lead to overestimation of the magnitude of the dimensionless heat flux by 7–15% for the “frozen” boundary layer case.     

Effect of pressure gradient of friction and heat transfer in a dusty boundary layer

Approximate analytic expressions for the local friction and heat transfer coefficients in a dusty laminar boundary layer are obtained and tested in the case of an incompressible carrier phase, power-law variation of the external gas flow velocity and small velocity and temperature phase disequilibrium. These expressions supplement the numerical analysis of the dusty boundary layer on a blunt body [1, 2] and the asymptotic calculation of the friction and heat transfer in a quasiequilibrium dusty gas boundary layer on a plate [3]. The combined effect of dustiness and pressure gradient on the friction and heat transfer coefficients is discussed. The results obtained can be used for the practical calculation of the friction and heat transfer in a quasiequilibrium dusty laminar boundary layer and for interpreting the corresponding experimental data. Tomsk. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 105–108, September–October, 1988.     

Boundary layer on a flat plate with strong longitudinal magnetic field

The influence of a magnetic field on the boundary layer on a flat plate in a sufficiently strongly ionized gas stream is studied. The magnetic field is parallel to the plate and to the velocity of the free stream, and it is so strong that the transport coefficients become anisotropic (the cyclotron rotation frequency of the charged particles is greater than or equal to the order of the frequency of the particle collisions). Using the results of [1–3] it is shown that the effect of the strong longitudinal magnetic field with a sufficiently high degree of gas ionization leads to a reduction in the thermal flux to the plate. For low degrees of ionization this effect is very small, since the viscosity and heat conduction in this case are determined by the neutral component of the gas.Results are presented of numerical calculations of the considered problem with account for the dependence of the transport coefficients on the thermodynamic parameters. It is assumed throughout that the magnetic Reynolds number is small (R_m1).     

Local structure of a binary finely dispersed fluidized bed

Results are given of calculations of the quantities characterizing the random pseudoturbulent motions of the phases in a homogeneous fluidized bed consisting of particles of two sorts, differing in size. The dependence of the coefficients of pseudoturbulent diffusion of the particles, the mean-square velocities of the pulsations, etc., on the partial concentrations of the particles, the ratio of their sizes, and other parameters is evaluated. For granular beds, fluidized by a gas or a drop-type liquid, intense chaotic fluctuations of both phases are characteristic; these determine to a considerable degree the observed macroscopic properties of the bed and affect its effectiveness as a working body in various types of heat exchangers and chemical reactors. Such random (pseudoturbulent) motions are particularly considerable for beds of small particles under homogeneous fluidization conditions, where mixing due to the rise of cavities in the bed, filled only with the fluidizing medium, is practically absent. A similar situation is encountered in reactor and regenerating units for catalytic cracking [1, 2], in beds with a drop-type liquid phase, in rarefied two-phase systems under the conditions of strong fluidization or of the transport of bulk materials in a dilute phase, etc. The characteristics of pseudoturbulence in locally homogeneous flows of monodisperse two-phase systems have been investigated, for example, in [3–5]. However, real fluidized beds are generally polydis-perse; the presence of particles of different sizes in the bed has a very considerable effect, on the intensity of the pulsations, the effective diffusion coefficients of the phases of the bed, the effective viscosities, etc. [1, 6]. In addition, the chaotic mixing in polydisperse beds determines some of the technological characteristics, specifically, the rate of entrainment of small particles by the flow of the fluidizing medium and the settling of large particles, the degree of separation of the fractions of the disperse phase, which is very important in determination of the limits of the existence of the fluidized state, and in the modeling of numerous processes of the separation of particles with respect to size or density [1, 6].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 118–126. March–April, 1976.     

Governing equations in non-isothermal diffusion

Generalizations of Fick's law for the diffusion flux are often considered in the literature by analogy with those for the heat flux. The paper reviews the balance equations for a fluid mixture and provides the equations for the diffusion fluxes. As a consequence, the mass densities are shown to satisfy a system of hyperbolic equations. Moreover, for a binary mixture of ideal gases in stationary conditions, Fick's law is recovered. Next, diffusion fluxes are regarded as constitutive functions and a whole set of thermodynamic restrictions are determined which account for diffusion, heat conduction, viscosity and inhomogeneities. Hyperbolic models for diffusion and heat fluxes are established which involve the co-rotational derivative. The driving term of diffusion turns out to be the gradient of chemical potential rescaled by the temperature.     

Two problems of convective diffusion to the surfaces of poorly streamlined bodies

Problems of diffusion to particles of nonspherical shape at large Peclet numbers have been analyzed in many papers (see [1–7], for example). The solution of the problem of mass exchange of an ellipsoidal bubble at low Reynolds numbers was obtained in [1] while the solution at high Reynolds numbers was obtained in [2, 3]. In [4] an expression is given for the diffusional flux to the surface of a solid ellipsoidal particle over which a uniform Stokes stream flows. Generalization to the case of particles of arbitrary shape was done in [5, 6], while generalization to any number of critical lines on the surface of the body was done in [7, 8]. The two-dimensional problem of steady convective diffusion to the surface of a body of arbitrary shape is analyzed in the approximation of a diffusional boundary layer (ADBL). The simple analytical expressions obtained are more suitable for practical calculations than those in [5-8], since they allow one to determine at once, in the same coordinate system in which the field of flow over the particle was analyzed, the value of the diffusional flux to its surface (from the corresponding hydrodynamic characteristics). The plane problem of the diffusion to an elliptical cylinder in a uniform Stokes stream is solved. The problems of the diffusion to a plate of finite dimensions (in the plane case) and a disk (in the axisymmetric case) whose planes are normal to the direction of the incident stream are considered. It is shown that, in contrast to the results known earlier (see [4, 6-15], for example), where the total diffusional flux was proportional to the cube root of the Peclet number, here it is proportional to the one-fourth power.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 104–109, November–December, 1978.The authors thank Yu. P. Gupalo, Yu. S. Ryazantsev, and Yu. A. Sergeev for a useful discussion.     

Chapman–Enskog solutions to arbitrary order in Sonine polynomials IV: Summational expressions for the diffusion- and thermal conductivity-related bracket integrals

The Chapman–Enskog solutions of the Boltzmann equations provide a basis for the computation of important transport coefficients for both simple gases and gas mixtures. These coefficients include the viscosity, the thermal conductivity, and the diffusion coefficient. In a preceding paper on simple gases, we have shown that the use of higher-order Sonine polynomial expansions enables one to obtain results of arbitrary precision that are free of numerical error. In two subsequent papers, we have extended our original simple gas work to encompass binary gas mixture computations of the viscosity, thermal conductivity, diffusion, and thermal diffusion coefficients to high-order. In all of this previous work we retained the full dependence of our solutions on the molecular masses, the molecular sizes, the mole fractions, and the intermolecular potential model via the omega integrals up to the final point of solution via matrix inversion. The elements of the matrices to be inverted are, in each case, determined by appropriate combinations of bracket integrals which contain, in general form, all of the various dependencies. Since accurate, explicit, general expressions for bracket integrals are not available in the literature beyond order 3, and since such general expressions are necessary for any extensive program of computations of the transport coefficients involving Sonine polynomial expansions to higher orders, we have investigated alternative methods of constructing appropriately general bracket integral expressions that do not rely on the term-by-term, expansion and pattern matching techniques that we developed for our previous work. It is our purpose in this paper to report the results of our efforts to obtain useful, alternative, general expressions for the bracket integrals associated with the diffusion- and thermal conductivity-related Chapman–Enskog solutions for gas mixtures. Specifically, we have obtained such expressions in summational form that are conducive to use in high-order transport coefficient computations for arbitrary gas mixtures and have computed and reported explicit expressions for all of the orders up to 5.     

Motion of a flat piston in a gas at high Reynolds numbers taking account of chemical reactions

One-dimensional flows originating during motion of a heat-conducting piston in a gas at high values of the Reynolds number are studied. The influence of diffusion and chemical reactions is considered in the case of a binary gas mixture. A binomial external expansion taking account of the boundary-layer-displacement thickness formed ahead of the piston is found. A solution is obtained which describes the boundary layer, which includes accommodation effects. An analogous problem about plane shock reflection from a heat-conducting wall has been considered in [1–3], but without taking account of diffusion and chemical reactions. Accomodation effects were taken into account in later work, which improved the agreement between theoretical and experimental results for short times.     

Effect of particles suspended in a fluid flow on the decay of isotropic turbulence

Dynamic equations have been obtained for the two-point double correlations of the fluctuation velocities of a fluid and the particles suspended in it at low volume concentrations of the solid phase. In the case of uniform isotropic turbulence these equations can be considerably simplified. The final period of decay of isotropic turbulence has been studied in detail. At this stage in the case of high-inertia particles the inhomogeneous-fluid turbulence is similar to the turbulence of a homogeneous fluid (without particles) in the sense that the presence of the particles affects only the fluctuation energy but leaves unchanged the spatial scales of turbulence and the spatial energy spectrum function. The suspended particles lead to exponential damping of the turbulent pulsations.Little theoretical information is available on the hydrodynamics of a suspension of fine particles in a turbulent liquid or gas. Research has been mainly confined to the behavior of the individual particles in a given turbulence field [1]. The problem of the turbulent motion of the mixture as a whole has been examined by Barenblatt [2], who derived the equations of motion of the mixture, using Kolmogorov's hypothesis to close them. Hinze [3] has also attempted to derive equations for turbulent pulsations of the mixture. However, as Murray showed [4], Hinze' s equations contradict Newton' s third law.The effect of suspended particles on the turbulence of a two-phase flow is governed by the noncorrespondence of the local velocities of the particles and the medium. The forces of resistance to the motion of the particles relative to the fluid lead to additional dissipation of fluctuation energy and decay of turbulence [2]. On the other hand, if the averaged velocities of particles and medium do not correspond, the suspended particles may also have a destabilizing effect [5, 6], causing energy transfer from the averaged to the pulsating motion. Below we shall consider the case where the averaged velocities of the two phases coincide, i.e., we shall deal only with the first of the two above-mentioned effects.The authors thank G.I. Barenblatt for his useful advice.     

Kinetic model of a gas suspension of particles with internal structure

Model kinetic equations are constructed for different relationships between the characteristic flow parameters of a two-phase medium consisting of diatomic molecules and solid particles with internal structure. The interphase collision integrals are represented as expansions with respect to the parameter of the ratio of the masses and in terms of physical quantities such as the transfer and diffusion coefficients, which depend on the model chosen for the transform of inelastic scattering on the surface. Transport equations are obtained for the light component with additional terms that take into account the interphase interactions. For one definite expression for the inelastic scattering transform simple analytic expressions are obtained for the additional terms in the equations of the gas dynamics.     

Supersonic two-phase flow around a wedge

Supersonic two-phase flow around bodies is encountered in calculating the flow around the last stages of blades of condensing turbines, in studying the motion of airplanes under cloudy conditions, etc. In the latter case, there is, along with erosion of the forward edges of the wing profiles, a change in the wave structure and interference situation in the flow about the airplane, leading to off-design regimes of motion. Supersonic flow of a two-phase mixture around a wedge, without taking account of the influence of the particles on the flow, was investigated in [1–3]. In [4], also in this kind of simplified setting, a study was made of the interaction of particles with the surface of a wedge in which reflection of the particles from the wall was taken into account. Morganthaler [5] made an experimental study of the flow of a mixture of air and aluminum oxide particles around a wedge. In [6] a theoretical study was made of a supersonic two-phase flow around thin flat axially-symmetric bodies. In particular, for the flow around a wedge, closed form solutions were obtained for the form of the shock wave, the gas streamlines and particle paths, and the distribution of all the parameters along the surface of the wedge. On the basis of the equations given in [7] and the method of characteristics, which were developed for flows consisting of a mixture of a gas and heterogeneous particles in nozzles [8,9], we present below a study of a supersonic two-phase flow around a wedge.Moscow. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 2, pp. 83–88, March–April, 1972.