Calculation of nonequilibrium flows in nozzles

This paper studies in the one-dimensional formulation the flow of a reacting gas with account for the nonequilibrium behavior of the chemical reactions; the pressure distribution along the stream filament is given. Viscosity, heat conduction, diffusion, and ionization are not taken into account. It is assumed that there is equilibrium excitation of the translational, rotational, and vibrational degrees of freedom.Several studies have already been made of nonequilibrium flows in nozzles [1–5]. It is known that in the calculation of nonequilibrium flows considerable difficulty arises in selecting the integration step in those regions where the flow is nearly equilibrium. It is found that with the use for numerical integration of the explicit difference schemes of the type of the Euler, Runge-Kutta, etc., methods the integration step for carrying out a stable calculation must be so small that the calculation becomes practically impossible. The present study proposes a method for calculating nonequilibrium flows using a single implicit difference scheme to calculate with a high degree of accuracy and a quite large step (exceeding the step in the explicit schemes by several orders) both those flow regions which are close to equilibrium and those regions where the flow deviates markedly from equilibrium. A program was compiled using this method for the M-20 electronic digital computer which permitted calculating in the one-dimensional approximation flows in nozzles with account for the nonequilibrium behavior of the chemical reactions for mixtures containing H, O, C, and N atoms.Some qualitative peculiarities of the nonequilibrium flows are demonstrated using as an example nonequilibrium air discharge. A comparison is made with experimental and theoretical results of other authors.The authors wish to thank L. F. Kuz'mina for her assistance in carrying out the present study.     

Investigation of supersonic flows in conical nozzles

The article gives the results of an investigation of flows in supersonic axisymmetric conical nozzles with the presence of shock waves in the flow. The method of straight-through calculation [1] and the method of small perturbations [2] are used. An investigation is made of the effect of various geometric parameters and of the adiabatic index on the flow of a gas in conical nozzles. A comparison is made with experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 101–107, May–June, 1974.In conclusion, the authors thank N. V. Drozdov for his participation in carrying out the calculations.     

Two-layer flows of gas in supersonic axisymmetric nozzles

This article describes two methods for calculating two-layer flows. The first is a generalization of a numerical method for solving the inverse problem [1] for the case of two-layer flows, without taking mixing into account. The second is a method of characteristics, for calculating a two-layer flow in a supersonic nozzle. In this case, the usual method of characteristics is changed in such a way that it is possible to calculate a point on the interface between two layers having different adiabatic indices, and different total pressures and temperatures. This article also gives the results of calculation of two-layer flows in nozzles with different adiabatic indices and different ratios of the mass flow rates of the gas in the layers.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 76–81, July–August, 1970.The calculations were programmed and carried out by G. D. Vladimirova and M. F. Tamarovskii, to whom the author expresses his thanks.     

Theory of slightly perturbed three-dimensional flows in nozzles

The theory of slightly perturbed flows in conical nozzles is used to determine the transverse force and moment generated in the presence of asymmetric perturbations. A system of ordinary differential equations is derived for finding the transverse force and moment. An approximate analytical solution of this system is constructed and its qualitative features are studied. A comparison is made with a numerical solution.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 146–154, January–February, 1977.     

Unsteady nonaxisymmetric flows in the hydrodynamic Czochralski model at high Prandtl numbers

Thermal-gravitational and thermocapillary convection is numerically modeled in the axisymmetric and three-dimensional approximations for the hydrodynamic model configuration corresponding to technological regimes of oxide crystal growth and taken as the basis in an international test. The salient features of the interaction between the convection and the flow driven by crystal rotation are studied at high Prandtl numbers. The flow and temperature fields occurring upon the generation of axisymmetric oscillations and the loss of axial symmetry are studied, analyzed, and compared with the results of other authors.     

Investigation of nonequilibrium two-phase flows in axisymmetric laval nozzles

The singularities of two-phase flows in Laval nozzles were investigated within the framework of the model of a two-fluid continuous medium [1, 2] mainly in a quasi-one-dimensional approximation ([3] and the bibliography therein). Two-dimensional computations of such flows were performed only recently by using the method of buildup [4–7]. However, systematic computations to clarify the influence of the second phase on such fundamental nozzle characteristics as the magnitude of the specific impulse, its losses, and discharge coefficient were performed only in the quasi-one-dimensional approximation [8, 9] and only for the supersonic parts of the nozzle in the two-dimensional approximation under the assumption of uniform flow in the throat [10, 3]. Such an investigation is performed in this paper in the two-dimensional case for the nozzle as a whole, including the sub-, trans-, and supersonic flow domains, and a comparative analysis is given of the magnitudes of the loss of a unit pulse obtained in the quasi-one-dimensional approximation [8].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 86–91, November–December, 1977.     

Theory of vortical helical ideal gas flows in laval nozzles

An asymptotic solution is found for the direct problem of the motion of an arbitrarily vortical helical ideal gas flow in a nozzle. The solution is constructed in the form of double series in powers of parameters characterizing the curvature of the nozzle wall at the critical section and the intensity of stream vorticity. The solution obtained is compared with available theoretical results of other authors. In particular, it is shown that it permits extension of the known Hall result for the untwisted flow in the transonic domain [1]. The behavior of the sonic line as a function of the vorticity distribution and the radius of curvature of the nozzle wall is analyzed. Spiral flows in nozzles have been investigated by analytic methods in [2–5] in a one-dimensional formulation and under the assumption of weak vorticity. Such flows have been studied by numerical methods in a quasi-one-dimensional approximation in [6, 7]. An analogous problem has recently been solved in an exact formulation by the relaxation method [8, 9]. A number of important nonuniform effects for practice have consequently been clarified and the boundedness of the analytical approach used in [2–7] is shown.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 126–137, March–April, 1978.The authors are grateful to A. N. Kraiko for discussing the research and for valuable remarks.     

Some peculiarities of two-phase flow in nozzles

Calculations are conducted for unidimensional two-phase flow in nozzles for a wide range of particle concentrations and dimensions. It is established that there exists a maximum in loss of specific momentum due to a lag in particle velocity and temperature relative to the gas. The results obtained are compared with calculations using linearized theory as well as with experimental data. The agreement between calculation and experiment is noted. Equilibrium flow of a two-phase mixture with solidification of liquid particles is considered. The presence of an anomalous flow region is established, where in the model of an ideal unidimensional equilibrium flow in a nozzle with discharge into a vacuum the presence of two successively located minimum sections is necessary.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 50–57, May–June, 1973.     

Study of spatial mixed flows in nozzles with asymmetrical inlets

The results of some investigations into the behavior of the spatial mixed flows of a nonviscous and nonheat-conducting gas in nozzles deviating from axisymmetrical shape in the subsonic region are presented. The investigation is based on the numerical integration of the transient gasdynamic equations using the Godunov scheme [1, 2] generalized to the spatial case.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 167–169, March–April, 1975.The author wishes to thank M. Ya. Ivanov, A. N., Kraiko, and G. G. Chernyi for help and discussion of the results.     

Inviscid gas flows in nozzles with a steeply converging contour

Inviscid gas flows in nozzles with a uniform exit flow and contours profiled starting from the lower point of a steeply converging region with an angle θ = −90° are analyzed. It is shown that there exists a class of convergent-divergent contours, within which the flow is characterized by the fact that the line θ = 0 of zero angle of the velocity vector inclination to the x axis consists of two oppositely-directed regions located partially or even completely ahead of the minimum section, while near the minimum sections their regions convex inward the gas stream are in decelerated flow. The minimum sections of the nozzles with M _e → 1 approach the center of the nozzle from the right.     

Use of flows with straight sonic line in nozzles with corners

The optimal scheme of a Laval nozzle is discussed. In the case of a profiled nozzle with a corner it is possible to use in the region of mixed flow both flows of general form with curvilinear sonic line as well as the special case when the sonic line is straight. It is shown that the latter alternative is preferable: when the supersonic part of the profile is determined by the simple wave method, the velocity at the wall increases monotonically and the flow does not contain shock waves. In contrast, in nozzles with curvilinear sonic line, a section in which the velocity decreases is formed immediately behind the corner, which can lead to boundary layer separation. In addition, for values of the supersonic velocity at the nozzle exit near the velocity of sound it is proved that the characteristics of the simple wave intersect in the flow region.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 168–170, January–February, 1981.     

Numerical computations of internal flows for axisymmetric and two-dimensional nozzles

MacCormack's explicit time-marching scheme is used to solve the full Navier–Stokes unsteady, compressible equations for internal flows. The requirement of a very fine grid to capture shock as well as separated flows is circumvented by employing grid clustering. The numerical scheme is applied for axisymmetric as well as two-dimensional flows. Numerical predictions are compared with experimental data and the qualitative as well as the quantitative agreement is found to be quite satisfactory. © 1997 John Wiley & Sons, Ltd.     

Shock and wave dynamics in cavitating compressible liquid flows in injection nozzles

Due to the exceptional high inlet pressures up to 2,000 bar flow dynamics and efficiency of modern injection systems are controlled by high frequency wave dynamics of the compressible liquid flow. Corresponding to alternating shock and expansion waves the liquid fluid evaporates and recondenses instantaneously. Here we present CFD simulations of the time accurate evolution of cavitating flows in 2-D plane and in six-hole injection nozzles with focus on the wave dynamics just after initialisation of the flow and within the time scale Δt ≤ 10^?4 s of pilot and multi-point injection. Due to shock reflections at the bottom of the sack hole the instantaneous maximum pressure increases more than three times higher as compared with the prescribed pressure at the nozzle inlet. For instance, in case of an inlet pressure of 600 bar the maximum pressure in the sack and therefore ahead of the nozzle bore holes reaches about 2,100 bar. It is quite reasonable that this amplification of the pressure affects the evolution of the convective flow and therefore the mass flow through the nozzle bore holes.     

Small amplitude pressure wave in subsonic and supersonic bubble flows in convergent-divergent nozzles

This paper makes a theoretical analysis of the propagation phenomena of the small amplitude pressure wave in the subsonic and supersonic bubble flow with a velocity slip between bubble and liquid in the convergent-divergent nozzle. From an analysis of the time-mean flow, the nondimensional parameter $\text{m = {u}\text{2}\text{G}\text{·α(1 ? α)ρ}\text{l}\text{β(2 ? 1/S)/P·[αβS + (1 ? α)βS}{}^2\text{+ α(1 ? α)]}}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$ corresponds to Mach number is gasdynamics where u_G is the gas velocity, α: the void fraction, ρ_L: the liquid density, P: the pressure, S: the velocity ratio of the gas and liquid flows and β: the proportional constant for the virtual mass. From a theoretical analysis of the small disturbance field, it is clarified that the parameter m also plays an essential and important role as Mach number, although the propagation performance of the disturbance is very complicated compared with that in gasdynamics. It is also shown that the pressure waves are divided into four groups depending on the velocity ratio S. Two of them are rather realistic, but the other two are required of a further investigation in future.     

Investigation of flows in hypersonic nozzles in the framework of simplified Navier-Stokes equations

Simplified Navier-Stokes equations have found application as an alternative to the complete Navier-Stokes equations for the simulation of viscous gas flows in regions of large dimensions, when there is a predominant direction of the flow [1–4]. In the present paper, flows in wind tunnel nozzles are investigated on the basis of this model. Flows in conical and profiled axisymmetric hypersonic nozzles are calculated in a wide range of Mach and Reynolds numbers. Good agreement with the experiment is obtained. The important role of viscous-inviscid interaction in nozzles for large hypersonic Mach numbers is shown. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 20–26, November–December, 1986. The representation of the nature of flows in hypersonic nozzles given here arose as a result of repeated discussions of the results of the calculations with my colleagues: V. N. Vetlutskii, V. L. Ganimedov, G. P. Klemenkov, Yu. G. Korobeinikov, and V. I. Pinchukov, for which the author is deeply grateful to them.     

Visualization and two-color DPIV measurements of flows in circular and square coaxial nozzles

 High-resolution, reactive Mie scattering laser-sheet visualizations, two-color digital particle image velocimetry (DPIV) and thermal anemometry measurements in flows generated by equivalent coaxial circular and square jets are presented. Visualization results were obtained for three square, coaxial configurations, and a reference circular coaxial nozzle, at two Reynolds numbers of the outer jet (19,000 and 29,000) and for inner-to-outer jet velocity ratios of 0.15, 0.22, and 0.3. These indicated that the internal unmixed region diminished with decreasing velocity ratio. Strong evidence of unsteady recirculation and back-flow was observed at the end of the core of the inner jet, for the low velocity ratios. Comparisons between circular and square jet configurations indicated considerable mixing enhancement when square nozzles were used. Low-coherence, organized large-scale structure was evident from the visualizations and DPIV measurements near the origin of the inner mixing-region shear layers, and more so in the core region of the near field. These observations were confirmed by velocity spectra, which displayed peaks corresponding to a free shear-layer instability mode in the inner mixing-region shear layers, and a wake-type mode in the core region where the mean flow has a wake-like character. Although some large-scale structure was observed in the outer mixing layer during the visualizations, this was found to be incoherent on the basis of the DPIV measurements and the velocity spectra. It is noted that no axis-switching phenomena were observed in the square nozzle flows examined here. This is attributed to the absence of an organized structure in the outer shear layer, which was initially highly turbulent, and the weakly coherent nature of the organized structure observed in the inner mixing-region near field. Received: 2 November 1998/Accepted: 8 September 2000