Calculation of the transport coefficients and slip velocities for molecules in the form of solid spheres

A solution of the Boltzmann equation is carried out by the Monte Carlo method for problems of rarefied gasdynamics in a linear formulation. The problems are solved by calculating the transport coefficients and slip velocities on a solid wall for molecules in the form of solid spheres. The accuracy of the method due to various parameters of the computational scheme in the solution of the problem is investigated by calculating the transport coefficients for pseudo-Maxwellian molecules.The Boltzmann kinetic equation is a complex integro-differential equation which is very difficult to solve and analyze. Hence, the solution of even one-dimensional problems and for the linearized Boltzmann equation turns out to be quite difficult, and such problems are solved by approximate methods (the expansion in Knudsen numbers, the method of moments, the expansion in series, etc. [1]). A method of solving the linearized Boltzmann equation by the Monte Carlo method is proposed in [2]. An exact solution of a number of problems of rarefied gas dynamics has been obtained by this method [3, 4]. However, the method was applied for pseudo-Maxwellian molecules, for which the collision cross section is inversely proportional to the relative velocity of the colliding particles =₀/g.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 155–158, March–April, 1971.In conclusion, the author is grateful to M. N. Kogan for formulating the problem and for great assistance provided during the research, and also to V. I. Vlasov, S. L. Gorelov and V. A. Perepukhov for assistance in compiling the program.     

Investigation of the hypersonic viscous shock layer around blunt bodies in a nonuniform flow

Unseparated viscous gas flow past a body is numerically investigated within the framework of the theory of a thin viscous shock layer [13–15]. The equations of the hypersonic viscous shock layer with generalized Rankine-Hugoniot conditions at the shock wave are solved by a finite-difference method [16] over a broad interval of Reynolds numbers and values of the temperature factor and nonuniformity parameters. Calculation results characterizing the effect of free-stream nonuniformity on the velocity and temperature profiles across the shock layer, the friction and heat transfer coefficients and the shock wave standoff distance are presented. The unseparated flow conditions are investigated and the critical values of the nonuniformity parameter a_k [10] at which reverse-circulatory zones develop on the front of the body are obtained as a function of the Reynolds number. The calculations are compared with the asymptotic solutions [10, 12].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 154–159, May–June, 1987.     

EHD equations and transport coefficients in a strong electric field

The methods of the kinetics theory are used to obtain a closed system of equations describing the behavior of a multicomponent, partially ionized gaseous mixture in an electromagnetic field in which the space charge is significant. A criterion is presented that makes it possible to separate the problem of finding the magnetic field from that of finding the other defining parameters. Expressions are obtained for the viscosity stress tensors and the thermal and diffusive fluxes; the transport coefficients are calculated in the strong electric field. The relations for the friction force and the energy exchange between components during particle collisions are analyzed. The equations for a mixture consisting of neutral particles and charged particles of a single sign are discussed in detail. The dimensionless EHD criteria are written and analyzed. Possible simplifications of the system of equations are examined, and various forms of Ohm's law are discussed. Weak discontinuities in EHD are analyzed. The equations of EHD under various assumptions have also been considered in several studies^* and in [1–3].     

Effect of inelastic collisions on the first-order slip coefficients for a diatomic gas with rotational degrees of freedom

When solving problems of inhomogeneous gas dynamics in the slip regime, it is necessary to know the boundary conditions for the velocity, temperature, heat fluxes, etc., that is, the boundary conditions for the gas macroparameters. In particular, such problems arise in developing the theory of thermophoresis of moderately large aerosol particles [1].The problem of monatomic and molecular (di- and polyatomic) gas slip along a boundary surface is considered in many publications (see, for example, [2–8]). The first-order effects include the isothermal and thermal gas slips characterized by the coefficients C_m and K_TS, respectively.In contrast to a monatomic gas, the molecules of diatomic and polyatomic gases have internal degrees of freedom, which considerably complicates the kinetic equation [9]. The lack of reliable models for the intermolecular interaction potential predetermines the need to construct model kinetic equations [10].In this study, for a diatomic gas whose molecules have rotational degrees of freedom, we propose a model kinetic equation obtained by developing the approach described in [6]. With the use of this model equation, the problem of diatomic gas slip along a plane surface is solved. As a result, for diatomic gases the coefficients C_m and K_TS, which depend on the thermophysical gas parameters and the intensity of inelastic collisions, are obtained.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, 2004, pp. 176–182. Original Russian Text Copyright © 2004 by Poddoskin.     

Flows in the initial sections of channels with permeable walls

Calculations of two types of flows in the initial sections of channels with permeable walls are carried out on the basis of semiempirical turbulence theories during fluid injection only through the walls and during interaction of the external flow with the injected fluid. Experimental studies of the first type [1–3] show that at least within the limits of the lengths L/h<30 and L/a< 50 (2h is the distance between permeable walls of a flat channel anda is the tube radius) the velocity distributions in the laminar and turbulent flow regimes differ little and are nearly self-similar for solutions obtained in [4]. For sufficiently large Reynolds numbers, Re₀>100 (Re₀=v₀h/ or Re₀=v₀ a/, where v₀ is the injection velocity), and small fluid compressibility, the axial velocity component is described by the relations for ideal eddying motion: u=(/2)x× cos (y/2) in a flat channel and u=x cos (y²/2) in atube (the characteristic values for the coordinates are, respectively, h anda). Measurements indicate the existence of a segment of laminar flow; its length depends on the Reynolds number of the injection [3]. In the turbulent regime the maximum generation of turbulent energy occurs significantly farther from the wall than in parallel flow. Flows of the second type in tubes were studied in [5–7]. These studies disclosed that for Reynolds numbers of the flow at the entrance to the porous part of the tube Re=u₀ a/<3.10³ fluid injection with v₀/u₀>0.01 leads to suppression of turbu lence in the initial section of the tube. An analogous phenomenon was observed in the boundary layer with v₀/u₀>0.023 [8, 9]. Laminar-turbulent transition in flows with injection was explained in [10, 11] on the basis of hydrodynamic instability theory, taking into account the non-parallel character of these flows. The mechanisms for the development of turbulence and reverse transition in channels with permeable walls are not theoretically explained. Simple semiempirical turbulence theories apparently are insufficient for this purpose. In the present work results are given of calculations with two-parameter turbulence models proposed in [12, 13] for describing complex flows. Due to the sharp changes of turbulent energy along the channel length, a numerical solution of the complete system of equations of motion was carried out by the finite-difference method [14].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 43–48, September–October, 1976.     

On some new solutions obtainable by means of invariant transformations

With reference to the example of the equations of monoenergetic nonrelativistic beam of particles of like charge, it is shown how new noninvariant solutions can be obtained by means of invariant transformations (§ 1). The conditions under which Lorentz forces can be ignored and the electric field considered a potential field are obtained for nonstationary flows. Solutions that describe the passage through a plane diode of high-frequency current from the emitter in a high-frequency electric field for an arbitrary relationship between the constant component of the collector potential and the amplitude of the ac voltage across it are derived (§2). Multivelocity (the velocity vector is a multivalued function) beams, and also electrostatic beams that can be described by Vlasov's equations are examined (§3).Given a system of differential equations (S) for m 1 unknown functions u^k (k=1,.,m) of n – m 1 independent variables xⁱ (i=1, ., n – m). The set of values (x, u) is considered as the set of coordinates of a point in n-dimensional space E_n. Any solution of this system u=u(x) defines some manifold in E_n. All possible solutions of (SI specify in E_n some set M. Any invariant transformation of system (S) has the property that it does not lead out of M. In a number of cases, this makes it possible to obtain new solutions by means of invariant transformations, no limitations being imposed on the solutions transformed. For a given system (S), all transformations that preserve (Si and form a continuous group, can be obtained by the method developed by L. V. Ovsyannikov [1–3]. Note that new solutions arise only when the principal group G of system (S) allows other than merely elementary transformations: magnifications, rotations, and translations are, as a rule, useless.Below, solutions of the equations of a monoenergetic nonrelativistic beam of particles of like charge are examined as an example [6–8].     

Nonlinear interactions of langmuir waves in a weakly inhomogeneous plasma

Kinetic equations for the scattering of the waves of the one-dimensional spectrum by plasma particles are obtained for a weakly inhomogeneous plasma. The equation for the evolution of the spectrum of the short waves [k² > (m_e/m_i) D_e ^–2] trapped in the inhomogeneities of the plasma density differs significantly from the kinetic equation for the waves in a homogeneous plasma. The problem of localization on the spectrum of the Langmuir waves in regions near the minima of the plasma density is also considered. A solution of the kinetic equation for the waves, which describes this process, is obtained.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 6–13, November–December, 1972.In conclusion, the author thanks A. S. Kingsep for suggesting the problem and for directing the work.     

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).     

The process of charging aerosol particles in an electric field,with allowance for diffusion of the ions

A study is made of the charging of highly conducting particles in disperse media. To this end, in the case of a small volume concentration of particles, there is a study, within the framework of electrohydrodynamics [1, 2], of the charging of an ideally conducting spherical particle in a unipolar charged gas. It is assumed that the particle is at rest relative to the gas, and its size and ion concentration are fairly small, so that the electric self-field of the ions may be neglected. A computer was used to calculate the dependences of the charging current on the charge of the particle for various values of the external electric field intensity. A comparison is made with the experimental data on the charging of particles [3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 177–180, January–February, 1986.The author wishes to express his gratitude to L. T. Chernyi for formulating the problem, and to V. V. Gogosov for a useful discussion of the results.     

The boundary layers of a fully ionized two-temperature plasma with given component temperatures at an electrode

The flow of a plasma with different component temperatures in the boundary layers at the electrodes of an MHD channel is investigated without any assumptions as to self-similarity. For the calculation of the electron temperature, the full energy equation for an electron gas [1] is solved with allowance for the estimates given in [2]. In contrast to [3, 4], the calculation includes the change in temperature of electrons and ions along the channel caused by the collective transport of energy, the work done by the partial pressure forces, and the Joule heating and the energy exchange between the components. The problem of the boundary layers in the flow of a two-temperature, partially ionized plasma past an electrode is solved in simplified form by the local similarity method in [5–7]. In these papers, either the Kerrebrock equation is used [5, 6] or the collective terms are omitted from the electron energy equation [7].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 3–10, September–October, 1972.The author thanks V. V. Gogosov and A. E. Yakubenko for interest in this work.     

Nonlinear waves in a liquid film entrained by a turbulent gas stream

In the long-wavelength approximation and on the basis of a simplified system of equations analogous to the one considered by Shkadov and Nabil' [1, 2], an investigation is made into waves of finite amplitude in thin films of a viscous liquid on the walls of a channel in the presence of a turbulent gas stream. A bibliography on the linear stability of such plane-parallel flows can be found in [3–5]. The nonlinear stability is considered in [6]. A stationary periodic solution is sought in the form of a Fourier expansion whose coefficients are found near the upper curve of neutral stability by Newton's method and near the lower branch of the stability curve by the method of Petviashvili and Tsvelodub [7, 8].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No, 2, pp. 37–42, March–April, 1981.I thank V. Ya. Shkadov for supervising the work and all the participants of G. I. Petrov's seminar for a helpful discussion.     

Nonlinear interaction of convective waves and the onset of turbulence in a rotating horizontal layer

Nonlinear development and interaction of disturbances in the classical Bernard-Rayleigh problem has been considered in a number of works [1, 2]. As a rule, the investigation leads to the construction and study of model equations of Landau's type [3–6]. In the present work the convection equation in the Boussinesq approximation is solved with the aid of direct methods [7–10]. We investigate the evolution of the disturbances which are a superposition of two waves with different wave numbers ₁ and ₂. We consider the appearance of harmonics of the form n₁ ± m₂, where n and m are integers. The main attention is given to problems of the onset of turbulence [11]. The numerical experments carried out showed that the consideration of the interaction of a limited number of harmonics (from 23 to 500) allows one to reproduce some characteristic features of the motion during the onset of turbulence.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 9–15, March–April, 1977.     

Determination of the constants of electrochemical reactions and the ion mobility coefficients in weakly conducting liquids from an analysis of the current—voltage characteristics

The problem is considered of the passage of a direct current through a solution of a weak electrolyte in a two-dimensional cell. Allowance is made for the electrochemical reactions of dissociation and reconbination which take place in the electrolyte when the rate of dissociation of the molecules is regarded as dependent on the electric field intensity [1–3]. For electrolytes whose recombination coefficient is of the order of magnitude of the Langevin coefficient, theoretical current—voltage characteristics are given for the limiting cases of large and small values of the characteristic times for the ion concentrations to be changed by electrochemical reactions and the transport of ions by the electric field. A method of determining the dissociation rate, the recombination coefficient, and the ion mobility coefficients is proposed on the basis of comparison of the theoretical and experimental current—voltage characteristics.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6. pp. 113–120, November–December, 1984.     

Messung der Wärmeleitfähigkeit von Stickstoff und Kohlenmonoxid bei hohen Temperaturen im Stoßwellenrohr

Zusammenfassung Die Meßergebnisse für die Wärmeleitfähigkeit von Stickstoff bei Temperaturen zwischen 1230 und 6000 K und Drückenzwischen 1 und 10 bar und von Kohlenmonoxid zwischen 1150 und 5000 K bei 1 bar werden mitgeteilt. Diese mit dem Stoßwellenrohr gemessenen Werte werden mit jenen verglichen, die sich aus der strengen kinetischen Gastheorie ergeben. Auch verfügbare Daten anderer Autoren werden zum Vergleich herangezogen.

---

Measurement of thermal conductivity of nitrogen and carbon monoxide at high temperatures in a shock tube

The paper presents results of shock-tube measurements of thermal conductivity of nitrogen at temperatures between 1230 and 6000 K and at pressures between 1 and 10 bar and of carbon monoxide at temperatures between 1150 and 5000 K at 1 bar. Experimental results are compared with several variants of theoretical values, computed from rigorous kinetic theory, and with available data of other authors.

---

Bezeichnungen (Einheiten in Klammern) a [m² s^–1] Temperaturleitzahl - C _p[J mol^–1 K^–1] molare Wärmekapazität - k [J K^–1] Boltzmann-Konstante - M [kg mol^–1] molare Masse - p bar Gesamtdruck - R [J mol^–1 K^–1] Gaskonstante - T [K] thermodynamische Temperatur - t [s] Zeit - U [J mol^–1] innere Energie - w [m s^–1] Geschwindigkeit - x [m] Ortskoordinate - x _i [1] Molanteil der Komponentei im Gasgemisch - [Wm^–1 K^–1] Wärmeleitfähigkeit - [mol m^–3] molare Konzentration Indizes i die Komponentei im Gasgemisch - g bezieht sich auf das (kalte) Gas bei der Wandtemperatur - w bezieht sich auf die feste Wand - p bei konstantem Druck Dieser Beitrag wurde auf dem Thermodynamik-Kolloquium des VDI im Oktober 1969 in Zürich vorgetragen.     

Transport coefficients of a dissociating gas

The calculation of the transport coefficients of a dissociating gas involves fundamental difficulties which arise when the internal degrees of freedom of the molecules are taken strictly into account. In practical calculations extensive use is made of the approximation proposed in [1], in the context of which the dependence of the diffusion velocity of the molecule on its internal state is totally neglected. In this case the expressions for the stress tensor and the diffusion velocities coincide with the corresponding expressions for a mixture of structureless particles; in the expression for the heat flux the diffusion transport of internal energy is taken only approximately into account. Here, analytic expressions for the diffusion velocities, heat flux and stress tensor are obtained without introducing simplifying assumptions. The calculation method is based on the results of [2], in which an approximate method of calculating the transport coefficients of a multicomponent mixture of structureless particles was proposed, and [3], in which the transport coefficients of a rotationally excited gas were calculated. The relations obtained are analyzed and compared with the existing results; their accuracy is estimated. A closed system of equations of gas dynamics is presented for a number of cases of practical importance.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 158–165, January–February, 1987.     

Plane one-dimensional stationary flow of an ideal charged gas in its own electric field

The plane one-dimensional flow of an incompressible gas consisting of a neutral and a charged component in its own electric field has been investigated by Stuetzer [1]. Stuetzer's results are valid when the electrostatic pressure is small compared with the hydraulic pressure. In the present paper an analogous problem is considered for a compressible gas under the more general assumption that the pressures are comparable. Three cases are analyzed: a) the velocity of the relative motion of the charged and the neutral particles is equal to zero; b) it is nonvanishing but the flow can be assumed to be approximately isentropic; c) a nonisentropic flow, i.e., one cannot ignore irreversible losses due to the relative motion of the charged and neutral particles. In the first two cases, closed solutions are obtained.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 24–32, January–February, 1971.I should like to thank I. V. Bespalov and Yu. M. Trushin for their interest and helpful comments.     

Convection in a closed cavity heated from below when equilibrium conditions are disrupted

The equilibrium of a fluid is possible in a closed cavity in the presence of a strictly vertical temperature gradient (heating from below) [1]. There is a distinct sequence of critical Rayleigh numbers R_i at which this equilibrium loses its stability relative to low characteristic perturbations. The presence of different finite perturbations, unavoidable in an experiment, leads to the absence of a strict equilibrium when R < R₁. The problem of the influence of the perturbation on the convection conditions near the critical points arises in this context [2, 3]. The case in which the cavity is heated not strictly from below is investigated in [2] and the case in which the perturbation of the equilibrium is due to the slow movement of the upper boundary of the region is investigated in [3]. In [2, 3] the perturbation has the structure of a first critical motion and thus the results of these papers coincide qualitatively. The perturbation of the temperature in the horizontal sections of the boundary, which creates a perturbation with a two-vortex structure corresponding to the second critical point R₂, is examined in this paper. A similar type of perturbation is characteristic for experiments in which the thermal conductivity properties of the fluid and the cavity walls are different. The nonlinear convection conditions are investigated numerically by the net-point method.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 203–207, March–April, 1977.The author wishes to thank D. B. Lyubimova, V. I. Chernatynskii, and A. A, Nepomnyashchii for their helpful comments.     

ON THE STRUCTURE OF SOLUTIONS OF LINEAR PARTIAL DIFFERENTIAL EQUATION sum form i+j≤n aijpiqjφ=0 WITH TWO INDEPENDENT VARIABLES AND CONSTANT COEFFICIENTS

This paper is a continuation of [1]. In this paper, the solutions of the more general linear partial differential equation form i+j≤n a_ijpⁱq^jφ=0 with two independent variables and constant coefficients are discussed in detail. The general solution which can be used in the approximation to the conditions of the definite solution of the practical problems is presented. To illustrate the use of the results obtained in this paper, some practical examples in mechanics are given.     

Investigation of supersonic isobaric submerged turbulent jets

Turbulent supersonic submerged air jets have been investigated on the Mach number interval M_a = 1.5–3.4 and on the interval of ratios of the total enthalpies in the external medium and the jet i₀ = 0.01 – 1. Oxyhydrogen jets with oxidizer ratios = 0.3–5 flowing from a nozzle at Mach numbers M_a = 1 and 2.4 have also been investigated. When < 1 the excess hydrogen in the jet burns up on mixing with the air. Special attention has been paid to obtaining experimental data free of the influence on the level of turbulence in the jet of the initial turbulence in the nozzle forechamber, shock waves occurring in the nozzle or in the jet at the nozzle exit, and the external acoustic field. The jet can be divided into two parts: an initial part and a main part. The initial part extends from the nozzle exit from the section x, in which the dimensionless velocity on the jet axis u_m = u_x/u_d = 0.75. Here, u_x is the velocity on the jet axis at distance x from the nozzle exit, and u_a is the nozzle exit velocity. The main part of the jet extends downstream from the section x. For the dimensionless length of the initial part x_m = x/d_a, where d_a is the diameter of the nozzle outlet section, empirical dependences on M_a and i₀ are obtained. It is shown, that in the main part of the jet the parameters on the flow axis — the dimensionless velocity and temperature — vary in inverse proportion to the distance, measured in units of length x, and do not depend on the flow characteristics which determine the length of the initial part of the jet. The angles of expansion of the viscous turbulent mixing layer in the submerged heated or burning jet increase with decrease in i₀ and M_a and are practically independent of the afterburning process.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza. No. 4, pp. 56–62, July–August, 1988.