Three-dimensional chemically nonequilibrium viscous shock layer on a catalytic surface with allowance for associated heat transfer

A three-dimensional flow of dissociating air past blunt bodies is investigated in the framework of the thin viscous shock layer theory. Multicomponent diffusion and homogeneous chemical reactions, including dissociation, recombination, and exchange reactions, are taken into account. The generalized Rankine-Kugoniot conditions are specified on the shock wave and the conditions which take into account the heterogeneous catalytic reactions, on the surface of the body. The viscous shock layer equations are solved together with the heat equations inside the coating, which is carbon with a deposited thin film of SiO₂, or quartz. The case of a thermally insulated surface is also considered. The problem for the case of the motion of a body along the re-entry trajectory into Earth's atmosphere is investigated numerically. The temperature of the surface and the heat flux toward it are given as a dependence on the height (tine) of the flight for different cases of the specification of the catalytic reactions. It is shown that the difference between the heat fluxes towards the thermally insulated surface and the fluxes toward the heat-conducting surface in the neighborhood of the stagnation point is of the order of 6–12% for all the cases considered. This makes it possible to decouple the solution of the problem of heat conduction in the body.Translated fron Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 140–146, November–December, 1985.deceased     

Turbulent boundary layer on a catalytic surface in a nonequilibrium dissociating gas

The majority of the studies which consider the flow of a dissociating gas in a turbulent boundary layer are devoted to the investigation of either frozen or equilibrium flows on a flat plate.The frozen turbulent boundary layer has been studied by Dorrance [1], Kutateladze and Leont'ev [2], and Lapin and Sergeev [3]. A study of the effect of catalytic recombination processes at the plate surface on the heat transfer in a frozen turbulent boundary layer was made by Lapin [4].Kosterin and Koshmarov [5], Ginzburg [6], Dorrance [7], and Lapin [8] have studied the turbulent boundary layer on a plate in equilibrium dissociating gas.The calculation of the heat transfer in a turbulent boundary layer on a catalytic plate surface with nonequilibrium dissociation was made by Kulgein [9]. In this study the nonequilibrium nature of the dissociation process was taken into account only in the laminar sublayer, while the flow in the turbulent core was considered frozen. The solution was found numerically using a computer by means of a laborious iteration process.The present paper reports a method for calculating the turbulent boundary layer on a flat catalytic plate with arbitrary dissociation rate. The method, constructed using the assumptions customary for turbulent boundary layer theory, is a successive approximation method. Good convergence of the method is assured by the fact that the effect of the nonequilibrium nature of the dissociation process on the parameter distribution in the boundary layer and, consequently, on the friction and heat transfer may be allowed for merely by finding corrections, usually relatively small, to the distribution of these parameters in the equilibrium or frozen flows. The basis of the study is the two-layer scheme of the turbulent boundary layer. The Prandtl and Schmidt numbers and also their turbulent analogs are taken equal to unity. As the model of the dissociating gas we use the Lighthill model of the ideal dissociating gas [10], extended by Freeman [11] to nonequilibrium flows.     

A hypersonic chemically nonequilibrium viscous shock layer oh wings with a catalytic surface

Within the framework of the theory of a hypersonic viscous shock layer a study is made of flow round wings of infinite span with blunt leading edges at various angles of attack and slip. Account is taken of multicomponent diffusion, and homogeneous chemical reactions, including dissociation-recombination reactions and exchange reactions. On the shock wave the generalized Rankine-Hugoniot conditions are given, and on the surface of the body conditions which allow for heterogeneous catalytic reactions of the first order with reaction rate constants depending [1] or not depending [2] on the temperature. The cases of an ideally catalytic and a noncatalytic surface are also considered. The surface of the body is assumed to be heatinsulated. A numerical study was made of the problem in a broad range of variation in the angles of attack and slip for different cases of prescribed constants representing the rates of the heterogeneous reactions. The conditions of the flow corresponded to the motion of a body which possess a lifting force along the trajectory of entry into the Earth's atmosphere [3]. The dependences are given of the equilibrium temperature of the surface along the stagnation line of the wing on the height of the flight and the distribution of this temperature along the surface of wings with parabolic and hyperbolic contours. It is shown that for flow regimes with a relatively high degree of dissociation in cases when the proportion of atoms recombined on the surface of the body is small, the dependences of the heat flow and the temperature of the surface on the angle of slip are of a nonmonotonic nature.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhldkosti i Gaza., No. 6, pp. 127–135, November–December, 1984.     

Increased heat transfer to a titanium surface with oxygen injection into the nonequilibrium boundary layer

The results of an experimental and numerical investigation of the heat transfer between a subsonic jet of dissociated nitrogen and a titanium surface, through which molecular oxygen is blown into the jet, are presented. It is established that in the nonequilibrium boundary layer regime the dependence of the heat flux on the injected oxygen flow rate is nonmonotonic. At a certain flow rate the heat transfer to the titanium surface reaches a maximum that considerably exceeds (by 20%) the heat transfer to an impermeable wall. The observed increase in heat transfer in the presence of injection is attributed to the interaction of the gas-phase exchange reactions and the recombination of atoms on the titanium surface, which has sharply different catalytic properties with respect to the recombination of nitrogen and oxygen atoms.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 148–155, July–August, 1991.     

Determination of the catalytic activity of materials by solving the equations of a nonequilibrium multicomponent boundary layer on a flat plate

A method is presented for determining the dependence of the probability of heterogeneous recombination γw from results of measurements of the heat flux Q_w to the surface of a catalytic sensor exposed to a pulsed supersonic flow of gas dissociated by an incident shock wave propagating in a shock tube. It is shown that the accuracy of the determination of γw depends not only on the accuracy of the measurements in the experiment, but also on the results of mathematical modeling of the flow of the dissociated gas over the surface of the body. Results from an analysis of an experiment are presented. Institute of Applied Mathematics and Mechanics at Tomsk University, Tomsk 634050. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 4, pp. 110–117, July–August, 1998.     

Multicomponent turbulent boundary layer on a chemically active surface

When a gaseous mixture flows past chemically active surfaces the boundary layer formed on the wetted body may contain a large number of components with different diffusion properties. This leads to the necessity for studying the diffusion of the components in the multicomponent boundary layer.The use of thebinary boundary layer concept in the general case cannot yield satisfactory results, since replacement of the mutual diffusion coefficients D_ij of the various pairs of components by a single diffusion coefficient D in many cases is a rough approximation.In the general case the number of different diffusion coefficients is equal to N(N–1)/2 (N is the number of components). Usually it is possible to identify groups of components with similar molecular weights. Then the number of different diffusion coefficients may be reduced without large error. However, even in the comparatively simple case when it is possible to divide all the components into two groups with similar molecular weights we must take account of three different diffusion coefficients (one diffusion coefficient in each group and also the diffusion coefficient for the components of one group relative to the components of the other group). Only in particular cases when the gaseous mixture consists of only two components with arbitrary molecular weights, or if all the components of the gaseous mixture have similar molecular weights, can we with justification introduce a single diffusion coefficient (if in this case there are no limitations on the direction of the diffusion).Studies have been published covering the laminar multicomponent boundary layer. An analytic method for solving the equations of the laminar multicomponent boundary layer was developed by Tirskii [1]. There are also studies in which concrete results were obtained by numerical methods with the use of computers (for example, [2, 3]).As far as the author knows, for turbulent flow there are studies (for example, [4, 5]) covering flow with chemical reactions only in the case when all the diffusion coefficients are equal (D_ij=D).The present paper presents a method for calculating the turbulent multicomponent boundary layer with account for several different diffusion coefficients.Notation x, y coordinates - u, v velocity components - density - T temperature - h heat content - H enthalpy - c_i mass concentration of the i-th component - c ₁ ⁽¹⁾ element concentrations in solid body - J_i diffusion flux of the i-th component - m molecular weight - dynamic viscosity coefficient - kinematic viscosity coefficient - heat conduction coefficient - c_p specific heat - adiabatic index - D_ij binary diffusion coefficients - P Prandtl number - S_ij Schmidt number - S_t Stanton number - M Mach number - friction - q radiant thermal flux - boundary layer thickness - D rate of displacement of gas-solid interface - degree of gasification - r_ij weight fraction of element i in component j - _ij stoichiometric coefficients - K_i reaction equilibrium constants - l number of components for which I_i0 Indices i, j component number - w quantities for y=0 - * quantities on the edge of the laminar sublayer - ⁽¹⁾ quantities at the solid body - quantities at the outer edge of the boundary layer - molar transport coefficients     

Slip boundary conditions on a catalytic surface in a multicomponent gas flow

The boundary conditions for the velocity slip and temperature and concentration jumps on the surface of a body in a rarefied multicomponent gas flow are obtained. The mathematical treatment is given in detail because of the need to examine critically some previous results which disagree with each other in spite of the fact that the initial premises and the methods of solution were the same. The results of this study, which are given in a convenient form, represent the boundary conditions for both the simplified and the complete Navier-Stokes equations in problems of hypersonic rarefied gas flow past bodies with a catalytically active surface.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 159–168, January–February, 1996.     

Rotation-vibration-dissociation interaction in a multicomponent nonequilibrium viscous shock layer

A new, simple and physically adequate method of calculating vibrationally nonequilibrium dissociation constants is proposed on the basis of a dissociation model which takes into account the equilibrium excitation of the rotational degrees of freedom of the molecules and the nonequilibrium excitation of vibrational quantum states. This rotation-vibration-dissociation interaction model contains only the indeterminacy associated with the indeterminacy of the experimental data on the interaction potentials and the collision cross sections of the components. In the case of thermodynamic equilibrium the model gives values of the dissociation constants close to those generally accepted. The use of this model in multicomponent nonequilibrium total viscous shock layer calculations gives values for the shock detachment distance within 5% of the experimental values. The indeterminacy in the values of the vibrational energy lost by air molecules during dissociation and recovered during recombination does not lead to serious errors in the macrocharacteristics of the flow. The nonequilibrium excitation of vibrational degrees of freedom proves to be not so important in computing the macrocharacteristics of the flow as previously assumed and the existing algorithms for calculating chemically nonequilibrium flows on the assumption of thermodynamic equilibrium can be used with satisfactory accuracy for calculating the values of the heat flux, the position of the shock wave, and the temperature and pressure in the shock layer for partially dissociated and ionized air.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 166–180, November–December, 1994.     

Definition and calculation of effective ambipolar diffusion coefficients for a laminar multicomponent ionized boundary layer

Effective diffusion coefficients substantially facilitate solution of detailed problems for multicomponent boundary layers with frozen-in reactions in the flow and heterogeneous reactions at the wall; they provide physically lucid correlation formulas and final equations for the convective heat flow to the undamaged solid as well as for the mass loss rate or effective erosion enthalpy if the walls are subject to thermochemical attack [1–5].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 60–72, July–August, 1970.     

Calculation of the chemically nonequilibrium flow of a multicomponent gas through a nozzle

We shall describe a method for calculating the flow of a complex gas mixture not in chemical equilibrium, based on selection of the time scale for the problem solution and sectioning of the reaction velocities in the near-equilibrium mode of their flow. The method permits derivation of calculation programs possessing wide applicability for systems with a large number of reactions. The application of the method is illustrated through calculations on a system, the components of which contain H,O,C,N, and Cl elements. Calculation of the parameters of the flow of a multicomponent gas through a nozzle under given conditions requires consideration of the kinetics of the chemical reactions which occur in the system. Mathematically, such a problem leads to the simultaneous numerical solution of the gas-dynamics equations for the flow parameters and the chemical kinetic equations for the concentration of the individual components. Several difficulties exist, the most basic of which are calculation of the region of near equilibrium flow and the transonic flow region with transition across the sonic point, and calculation of a large number of chemical reaction velocities of greatly differing magnitudes, in the case of a complex gas mixture. In order to obtain a stable solution in the near-equilibrium flow region, several methods have recently been proposed, which permit consolidation of the integration step. We note the use of a local linearization of the chemical kinetic equation system, as employed in [1]. This method in practice is useful for relatively slow change in component concentrations. In [2] at each integration step the kinetic equations are transformed into a system of L algebraic equations (where L is the number of reactions), and with an increase in the number of reactions (L20) the laboriousness of such a calculation increases sharply. The implicit differential schemes of integration presented in [3, 4] appear more acceptable, but in fact they too have been tested only for systems with a relatively small number of reactions. The difficulty of calculating the transonic flow region, as is well known, is connected with selection of the unique value of mass flow G, at which the transition to supersonic flow is realized. This may be avoided by defining over the length of the nozzle one of the gas-dynamic functions (for example, pressure distribution [4]), which are not highly sensitive to chemical nonequilibrium, the values being taken from supplementary calculations of nonequilibrium flow through the nozzle. Several investigators have limited their examination to the supersonic flow region (see, for example, [5]), but with this method the results may lack sufficient accuracy, since in some cases (for high gradients of the gas-dynamic magnitudes) the transonic region produces a comparable contribution to the general effect of nonequilibrium. We will describe below a method with which a practically universal system of calculating the nonequilibrium flow of a complex gas mixture through a nozzle can be realized. In practice, up to 60 of the most significant reactions may be considered, out of a practically unlimited number initially present. The method is based on sectioning of the reaction velocities in the near-equilibrium mode. This permits attainment of a stable solution in the near-equilibrium flow region with acceptable machine-time expenditures. The method describes the transition through the sonic point, since the systematic error introduced by its use (within the limits of calculation accuracy) improves the convergence of the iteration process used in finding G. In applying the method it is useful to select the more important of those reactions theoretically possible, and also to conduct calculations for equilibrium flow conditions of the individual chemical reactions; the latter permits evaluation of the maximum possible contribution of reactions for which the velocity constants are unknown.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 159–163, September–October, 1971.The authors are grateful to A. I. Vol'pert and L. N. Stesik for their interest in and evaluation of the study.     

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.     

Translationally nonequilibrium free jet boundary layer in the wake behind a wedge

Macroscopic equations obtained as a thin-layer version of the 13-moment Grad equations derived from kinetic considerations are used for describing the translationally nonequilibrium monatomic gas flow in a hypersonic free jet boundary layer formed in the wake behind a wedge. This model makes it possible to investigate flows with strong violations of equilibrium with respect to the translational degrees of freedom. A method of constructing the solution of this kinetically justified problem based on the solution of an analogous problem in the Navier-Stokes interpretation is proposed. It is established that for the kinetic variant of the problem considered the gas flow velocity distribution along the separating streamline in a plane orthogonal to the wedge generator coincides with the distribution obtained in solving the Navier-Stokes variant. It is found that taking into account the nonequilibrium nature of the flow with respect to the translational degrees of freedom of the gas particles has no effect on the base pressure and the wake angle.     

Nonlocal approach to nonequilibrium thermodynamics and nonlocal heat diffusion processes

We study some aspects of nonequilibrium thermodynamics and heat diffusion processes based on Suykens’s nonlocal-in-time kinetic energy approach recently introduced in the literature. A number of properties and insights are obtained in particular the emergence of oscillating entropy and nonlocal diffusion equations which are relevant to a number of physical and engineering problems. Several features are obtained and discussed in details.     

Effect of surface catalytic activity on nonequilibrium heat transfer in a subsonic jet of dissociated nitrogen

The present paper investigates experimentally and numerically the effect of the heterogeneous recombination of atoms on the heat transfer of models in a subsonic jet of dissociated nitrogen for the conditions of an experiment in the VGU-2 plasma generator and determines the effective probabilities of the heterogeneous recombination of nitrogen atoms for a number of materials at high temperatures.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 166–172, May–June, 1985.