Existence and Stability of Multidimensional Transonic Flows through an Infinite Nozzle of Arbitrary Cross-Sections

We establish the existence and stability of multidimensional steady transonic flows with transonic shocks through an infinite nozzle of arbitrary cross-sections, including a slowly varying de Laval nozzle. The transonic flow is governed by the inviscid potential flow equation with supersonic upstream flow at the entrance, uniform subsonic downstream flow at the exit at infinity, and the slip boundary condition on the nozzle boundary. Our results indicate that, if the supersonic upstream flow at the entrance is sufficiently close to a uniform flow, there exists a solution that consists of a C ^1,α subsonic flow in the unbounded downstream region, converging to a uniform velocity state at infinity, and a C ^1,α multidimensional transonic shock separating the subsonic flow from the supersonic upstream flow; the uniform velocity state at the exit at infinity in the downstream direction is uniquely determined by the supersonic upstream flow; and the shock is orthogonal to the nozzle boundary at every point of their intersection. In order to construct such a transonic flow, we reformulate the multidimensional transonic nozzle problem into a free boundary problem for the subsonic phase, in which the equation is elliptic and the free boundary is a transonic shock. The free boundary conditions are determined by the Rankine–Hugoniot conditions along the shock. We further develop a nonlinear iteration approach and employ its advantages to deal with such a free boundary problem in the unbounded domain. We also prove that the transonic flow with a transonic shock is unique and stable with respect to the nozzle boundary and the smooth supersonic upstream flow at the entrance.     

Stability of Transonic Shock Solutions for One-Dimensional Euler–Poisson Equations

In this paper, both structural and dynamical stabilities of steady transonic shock solutions for one-dimensional Euler–Poisson systems are investigated. First, a steady transonic shock solution with a supersonic background charge is shown to be structurally stable with respect to small perturbations of the background charge, provided that the electric field is positive at the shock location. Second, any steady transonic shock solution with a supersonic background charge is proved to be dynamically and exponentially stable with respect to small perturbations of the initial data, provided the electric field is not too negative at the shock location. The proof of the first stability result relies on a monotonicity argument for the shock position and the downstream density, and on a stability analysis for subsonic and supersonic solutions. The dynamical stability of the steady transonic shock for the Euler–Poisson equations can be transformed to the global well-posedness of a free boundary problem for a quasilinear second order equation with nonlinear boundary conditions. The analysis for the associated linearized problem plays an essential role.     

Transonic Shocks in Compressible Flow Passing a Duct for Three-Dimensional Euler Systems

In this paper we study the transonic shock in steady compressible flow passing a duct. The flow is a given supersonic one at the entrance of the duct and becomes subsonic across a shock front, which passes through a given point on the wall of the duct. The flow is governed by the three-dimensional steady full Euler system, which is purely hyperbolic ahead of the shock and is of elliptic–hyperbolic composed type behind the shock. The upstream flow is a uniform supersonic one with the addition of a three-dimensional perturbation, while the pressure of the downstream flow at the exit of the duct is assigned apart from a constant difference. The problem of determining the transonic shock and the flow behind the shock is reduced to a free-boundary value problem. In order to solve the free-boundary problem of the elliptic–hyperbolic system one crucial point is to decompose the whole system to a canonical form, in which the elliptic part and the hyperbolic part are separated at the level of the principal part. Due to the complexity of the characteristic varieties for the three-dimensional Euler system the calculus of symbols is employed to complete the decomposition. The new ingredient of our analysis also contains the process of determining the shock front governed by a pair of partial differential equations, which are coupled with the three-dimensional Euler system. The paper is partially supported by National Natural Science Foundation of China 10531020, the National Basic Research Program of China 2006CB805902, and the Doctorial Foundation of National Educational Ministry 20050246001.     

Supersonic flow of a combustible gas mixture past a sphere

In recent years considerable interest has developed in the problems of steady-state supersonic flow of a mixture of gases about bodies with the formation of detonation waves and slow combustion fronts. This is due in particular to the problem of fuel combustion in a supersonic air stream.In [1] the problem of supersonic flow past a wedge with a detonation wave attached to the wedge apex is solved. This solution is based on using the equation of the detonation polar obtained in [2]-the analog of the shock polar for the case of an exothermic discontinuity. In [3] a solution is given of the problem of cone flow with an attached detonation wave, and [4] presents solutions of the problems of supersonic flow past the wedge and cone with the formation of attached adiabatic shocks with subsequent combustion of the mixture in slow combustion fronts. In the two latter studies two different solutions were also found for the problem of flow past a point ignition source, one solution with gas combustion in the detonation wave, the other with gas combustion in the slow combustion front following the adiabatic shock. These solutions describe two different asymptotic pictures of flow of a combustible gas mixture past bodies.In an experimental study of the motion of a sphere in a combustible gas mixture [5] it was found that the detonation wave formed ahead of the sphere splits at some distance from the body into an ordinary (adiabatic) shock and a slow combustion front. Arguments are presented in [6] which make it possible to explain this phenomenon and in certain cases to predict its occurrence.The present paper presents examples of the calculation of flow of a combustible gas mixture past a sphere with a detonation wave in the case when the wave does not split. In addition, the flow near the point at which the detonation wave splits is analyzed for the case when splitting occurs where the gas velocity behind the wave is greater than the speed of sound. This analysis shows that in the given case the flow calculation may be carried out without any particular difficulties. On the other hand, the calculation of the flow for the case when the point of splitting is located in the subsonic portion of the flow behind the wave (or in the region of influence of the subsonic portion of the flow) presents difficulties. This flow case is similar to the problem of the supersonic jet of finite width impacting on an obstacle.     

Inverse problem of wing aerodynamics in a supersonic flow

The inverse problem of wing aerodynamics—the determination of the lifting surface shape from a specified load—is solved within the framework of linear theory. Volterra's solution of the wave equation is used. Solutions are found in the class of bounded functions if certain conditions imposed on the governing parameters of the problem are satisfied. Solutions of inverse problems of supersonic flow are presented for an infinite-span wing, a triangular wing with completely subsonic edges, and a rectangular wing. Institute of Theoretical and Applied Mechanics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 3, pp. 86–91, May–June, 1998.     

Laval nozzle flow calculation

The inverse problem of the theory of the Laval nozzle is considered, which leads to the Cauchy problem for the gasdynamic equations; the streamlines and the flow parameters are found from the known velocity distribution on the axis of symmetry.The inverse problem of Laval nozzle theory was considered in 1908 by Meyer [1], who expanded the velocity potential into a series in powers of the Cartesian coordinates and constructed the subsonic and supersonic solutions in the vicinity of the center of the nozzle. Taylor [2] used a similar method to construct a flowfield which is subsonic but has local supersonic zones in the vicinity of the minimal section. Frankl [3] and Fal'kovich [4] studied the flow in the vicinity of the nozzle center in the hodograph plane. Their solution, just as the Meyer solution, made it possible to obtain an idea of the structure of the transonic flow in the vicinity of the center of the nozzle.A large number of studies on transonic flow in the vicinity of the center of the nozzle have been made using the method of small perturbations. The approximate equation for the transonic velocity potential in the physical plane, obtained in [3–6], has been studied in detail for the plane and axisymmetric cases. In [7] Ryzhov used this equation to study the question of the formation of shock waves in the vicinity of the center of the nozzle, and conditions were formulated for the plane and axisymmetric cases under which the flow will not contain shock waves. However, none of the solutions listed above for the inverse problem of Laval nozzle theory makes it possible to calculate the flow in the subsonic and transonic parts of the nozzles with large gradients of the gasdynamic parameters along the normal to the axis of symmetry.Among the studies devoted to the numerical calculation of the flow in the subsonic portion of the Laval nozzle we should note the study of Alikhashkin et al., and the work of Favorskii [9], in which the method of integral relations was used to solve the direct problem for the plane and axisymmetric cases.The present paper provides a numerical solution of the inverse problem of Laval nozzle theory. A stable difference scheme is presented which permits analysis with a high degree of accuracy of the subsonic, transonic, and supersonic flow regions. The result of the calculations is a series of nozzles with rectilinear and curvilinear transition surfaces in which the flow is significantly different from the one-dimensional flow. The flowfield in the subsonic and transonic portions of the nozzles is studied. Several asymptotic solutions are obtained and a comparison is made of these solutions with the numerical solution.The author wishes to thank G. D. Vladimirov for compiling the large number of programs and carrying out the calculations on the M-20 computer.     

Performance evaluation of an inverse integral equation method applied to turbomachine cascades

An improved formulation of the inverse integral equation method proposed in Reference 1 is presented which allows, in particular, a well-posed problem to be ensured. The corresponding computation code is tested in an exhaustive manner for axial and radial compressor and turbine cascades. The agreement between the velocity field obtained with the inverse method and that resulting from a direct calculation is examined for subsonic, transonic and supersonic flows. Accuracy and reliability of the solution to the boundary condition problem are excellent for the subsonic and transonic flows. However, for the supersonic flow, the application of the method seems to be limited by the use of elementary solutions of the Laplace operator.     

Direct and inverse problems of flow over a wing of finite span in the linear formulation

A correspondence between the solutions of the direct and the inverse problem for wing theory is established for a wing of finite span in the framework of linear theory on the basis of solution of a wave equation in Volterra form for supersonic flow and solution of the Laplace equation in the form of Green's formula for subsonic flow. For the direct problem in the case of supersonic flow an expression is derived for finding the load on the wing with maximal allowance for the wing geometry. In the inverse problem for supersonic and subsonic flows, expressions are derived for finding the wing geometry from given values of the load on the wing and the variation of the load along the span of the wing. The solution of the inverse problem is presented in the form of integrals that converge for interior points of the wing surface in the sense of the Cauchy principal value, the wing surface being represented as a vortex surface of mutually orthogonal vortex lines. The conditions of finiteness of the velocities on the edges are discussed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 114–125, September–October, 1979.     

Numerical calculations of supersonic underexpanded jets in a cocurrent flow with the use of parabolized Navier-Stokes equations

Results of a numerical study of three-dimensional supersonic jets propagating in a cocurrent flow are described. Averaged parabolized Navier-Stokes equations are solved numerically on the basis of a developed scheme, which allows calculations in supersonic and subsonic flow regions to be performed in a single manner. A jet flow with a cocurrent flow Mach number 0.05 ⩽ M_∞ ⩽ 7.00 is studied, and its effect on the structure of the mixing layer is demonstrated. The calculated results are compared with available experimental and numerical data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 54–63, May–June, 2008.     

The shock–vortex interaction patterns affected by vortex flow regime and vortex models

We have used a third-order essentially non-oscillatory method to obtain numerical shadowgraphs for investigation of shock–vortex interaction patterns. To search different interaction patterns, we have tested two vortex models (the composite vortex model and the Taylor vortex model) and as many as 47 parametric data sets. By shock–vortex interaction, the impinging shock is deformed to a S-shape with leading and lagging parts of the shock. The vortex flow is locally accelerated by the leading shock and locally decelerated by the lagging shock, having a severely elongated vortex core with two vertices. When the leading shock escapes the vortex, implosion effect creates a high pressure in the vertex area where the flow had been most expanded. This compressed region spreads in time with two frontal waves, an induced expansion wave and an induced compression wave. They are subsonic waves when the shock–vortex interaction is weak but become supersonic waves for strong interactions. Under a intermediate interaction, however, an induced shock wave is first developed where flow speed is supersonic but is dissipated where the incoming flow is subsonic. We have identified three different interaction patterns that depend on the vortex flow regime characterized by the shock–vortex interaction.      

Stability and nonlinear wavy regimes in downward film flows on a corrugated surface

The linear and nonlinear stability of downward viscous film flows on a corrugated surface to freesurface perturbations is analyzed theoretically. The study is performed with the use of an integral approach in ranges of parameters where the calculated results and the corresponding solutions of Navier-Stokes equations (downward wavy flow on a smooth wall and waveless flow along a corrugated surface) are in good agreement. It is demonstrated that, for moderate Reynolds numbers, there is a range of corrugation parameters (amplitude and period) where all linear perturbations of the free surface decay. For high Reynolds numbers, the waveless downward flow is unstable. Various nonlinear wavy regimes induced by varying the corrugation amplitude are determined. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 110–120, January–February, 2007.     

Development of slow disturbances in a hypersonic boundary layer

The mechanisms of development of slow time-dependent disturbances in the wall region of a hypersonic boundary layer are established and a diagram of the disturbed flow patterns is plotted; the corresponding nonlinear boundary value problem is formulated for each of these regimes. It is shown that the main factors that form the disturbed flow are the gas enthalpy near the body surface, the local viscous-inviscid interaction level, and the type, either subsonic or supersonic, of the boundary layer as a whole. Numerical and analytical solutions are obtained in the linear approximation. It is established that enhancement of the local viscous-inviscid interaction or an increased role for the main supersonic region of the boundary layer makes the disturbed flow by and large “supersonic”: the upstream propagation of the disturbances becomes weaker, while their downstream growth is amplified. Contrariwise, local viscous-inviscid interaction attenuation or an increased role for the main subsonic region of the boundary layer has the opposite effect. Surface cooling favors an increased effect of the main region of the boundary layer while heating favors an increased wall region effect. It is also found that in the regimes considered disturbances travel from the turbulent flow region downstream of the disturbed region under consideration counter to the oncoming flow, which may be of considerable significance in constructing the nonlinear stability theory.     

Calculation of hydrogen-air combustion behind detached shock wave for supersonic flow past a sphere

Many of the published theoretical studies of quasi-one-dimensional flows with combustion have been devoted to combustion in a nozzle, wake, or streamtube behind a normal shock wave [1–6].Recently, considerable interest has developed in the study of two-dimensional problems, specifically, the effective combustion of fuel in a supersonic air stream.In connection with experimental studies of the motion of bodies in combustible gas mixtures using ballistic facilities [7–9], the requirement has arisen for computer calculations of two-dimensional supersonic gas flow past bodies in the presence of combustion.In preceding studies [10–12] the present author has solved the steady-state problem under very simple assumptions concerning the structure of the combustion zone in a detonation wave.In the present paper we obtain a numerical solution of the problem of supersonic hydrogen-air flow past a sphere with account for the nonequilibrium nature of eight chemical reactions. The computations encompass only the subsonic and transonic flow regions.The author thanks G. G. Chernyi for valuable comments during discussion of the article.     

Local subsonic zones in supersonic jet flows

The results of the experimental investigation of supersonic turbulent jets with local subsonic zones of forward and reverse flow exhausting into the ambient atmosphere or an outer stream, either parallel or transverse to the jet, are presented. Some gasdynamic features of the flows containing these zones, which have not been adequately addressed in the literature, are analyzed. Thus, supersonic flows with back pressure, e.g., highly overexpanded and underexpanded jet flows, and those upstream and downstream of a jet on the leeward side of a cone in a supersonic gas stream, are studied. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 143–150, January–February, 1998.     

DEVELOPMENT OF A THREE-DIMENSIONAL VISCOUS AEROELASTIC SOLVER FOR NONLINEAR PANEL FLUTTER

A new three-dimensional (3-D) viscous aeroelastic solver for nonlinear panel flutter is developed in this paper. A well-validated full Navier–Stokes code is coupled with a finite-difference procedure for the von Karman plate equations. A subiteration strategy is employed to eliminate lagging errors between the fluid and structural solvers. This approach eliminates the need for the development of a specialized, tightly coupled algorithm for the fluid/structure interaction problem. The new computational scheme is applied to the solution of inviscid two-dimensional panel flutter problems for subsonic and supersonic Mach numbers. Supersonic results are shown to be consistent with the work of previous researchers. Multiple solutions at subsonic Mach numbers are discussed. Viscous effects are shown to raise the flutter dynamic pressure for the supersonic case. For the subsonic viscous case, a different type of flutter behavior occurs for the downward deflected solution with oscillations occurring about a mean deflected position of the panel. This flutter phenomenon results from a true fluid/structure interaction between the flexible panel and the viscous flow above the surface. Initial computations have also been performed for inviscid, 3-D panel flutter for both supersonic and subsonic Mach numbers.     

Supersonic flow past a blunt wedge

A study is made of the asymptotic solution of the problem of flow past a blunt wedge by a uniform supersonic stream of perfect gas. By separation of variables it is shown that at large distances the disturbance of the flow is damped exponentially. In the case of subsonic flow behind the shock wave the exponent of the leading correction term in the expansion of the shock front is calculated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 137–140, July–August, 1984.     

The blunt-body problem in a nonuniform hypersonic viscous gas flow

The laws of heat transfer associated with the interaction of underexpanded supersonic gas jets and obstacles or blunt bodies have been investigated, for example, in [1–3]. Similar problems of nonuniform flow occur when bodies move in the wake behind other bodies; however, in this case the laws of heat transfer have so far received little attention [4–8]. It has been established that for a certain Reynolds number and flow nonuniformity parameters a zone of reverse-circulatory flow develops near the front of the blunt body. However, the conditions of transition to separated flow have not been determined. This paper presents a self-similar solution of the equations of the viscous shock layer near the stagnation line in supersonic flow past an axisymmetric blunt body located behind another body. On the basis of this solution a separationless flow criterion is proposed. The effect of the nonuniformity and the Reynolds number on the shock standoff distance, the convective heat flux and the friction drag of the blunt body is investigated. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 120–125, November–December, 1986. In conclusion the authors wish to thank I. G. Eremeitsev for useful suggestions and G. A. Tirskii for discussing their work.     

Investigation of a Gas Flow with Solid Particles in a Supersonic Nozzle

A two-phase flow with high Reynolds numbers in the subsonic, transonic, and supersonic parts of the nozzle is considered within the framework of the Prandtl model, i.e., the flow is divided into an inviscid core and a thin boundary layer. Mutual influence of the gas and solid particles is taken into account. The Euler equations are solved for the gas in the flow core, and the boundary-layer equations are used in the near-wall region. The particle motion in the inviscid region is described by the Lagrangian approach, and trajectories and temperatures of particle packets are tracked. The behavior of particles in the boundary layer is described by the Euler equations for volume-averaged parameters of particles. The computed particle-velocity distributions are compared with experiments in a plane nozzle. It is noted that particles inserted in the subsonic part of the nozzle are focused at the nozzle centerline, which leads to substantial flow deceleration in the supersonic part of the nozzle. The effect of various boundary conditions for the flow of particles in the inviscid region is considered. For an axisymmetric nozzle, the influence of the contour of the subsonic part of the nozzle, the loading ratio, and the particle diameter on the particle-flow parameters in the inviscid region and in the boundary layer is studied. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 6, pp. 65–77, November–December, 2005.     

Computation of base flow characteristics with uniform blowing

Blowing at bluff body base was considered under different conditions and for small amount of blowing this problem was solved using dividing streamline model [1]. The effect of supersonic blowing on the flow characteristics of the external supersonic stream was studied in [2–4]. The procedure and results of the solution to the problem of subsonic blowing of a homogeneous fluid at the base of a body in supersonic flow are discussed in this paper. Analysis of experimental results (see, e.g., [5]) shows that within a certain range of blowing rate the pressure distribution along the viscous region differs very little from the pressure in the free stream ahead of the base section. In this range the flow in the blown subsonic jet and in the mixing zones can be described approximately by slender channel flow. This approximation is used in the computation of nozzle flows with smooth wall inclination [6, 7]. On the other hand, boundary layer equations are used to compute separated stationary flows with developed recirculation regions [8] in order to describe the flow at the throat of the wake. The presence of blowing has significant effect on the flow structure in the base region. An increasing blowing rate reduces the size of the recirculation region [9] and increases base pressure. This leads to a widening of the flow region at the throat, usually described by boundary-layer approximations. At a certain blowing rate the recirculation region completely disappears which makes it possible to use boundary-layer equations to describe the flow in the entire viscous region in the immediate neighborhood of the base section.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 76–81, January–February, 1984.     

One method of profiling short plane nozzles

One of the possible methods is considered for profiling short plane nozzles for aerodynamic tubes. The nozzle has a straight sonic line, which allows the subsonic and supersonic sections to be constructed separately. The problem is solved numerically in the plane of a hodograph. In the subsonic region, Dirichlet's problem is formulated for Chaplygin's equation in a rectangle, one side of which is the sonic line. At the present time, two approaches have been defined in papers on calculations of a Laval nozzle, associated with the solution of the so-called “direct” and “inverse” problems (one has in mind a study of the flow in the interconnected region of sub- and supersonic flow). The direct problem determines the flow field in the case of a previously specified contour of the channel wall, the shape of which from technical considerations is obtained with certain geometry conditions. The direct problem can be applied in the construction of the Laval nozzle, if the contour of the inlet section of the channel (generally speaking, quite arbitrary) is chosen so successfully that neither shock compressions nor breakaway zones result in the flow. Although a strictly mathematical theory of the direct problem of the Laval nozzle is only being developed at present, there are still very effective numerical methods for its solution [1, 2]. In the inverse problem (which, by definition, is a problem of profiling), the contour of the nozzle is found with respect to a specified velocity distribution on the axis of symmetry. It is assumed that this quite arbitrary dependence can be selected from the condition of the absence of breakaway zones and shock compressions in the nozzle. By its formulation, the inverse problem is Cauchy's problem which, as is well-known, is incorrect in the classical sense in the ellipticity region — the subsonic section of the nozzle. At present, there are also efficient methods of solving the inverse nozzle problem [3], by interpreting it as an arbitrarily correct problem. Difficulties can arise in the inverse problem, in the provision of short (and, consequently, steep) nozzles because of the sharp increase of the error in the calculation. Together with the stated problems, a procedure can be evolved which is associated with the solution of the correctly posed problem for Chaplygin's equation in the plane of the hodograph. This approach is convenient in that it succeeds a priori in fulfilling the important condition of monotonicity of the velocity at the wall, ensuring (in the absence of shock compressions) nonseparability of the streamline flow at any Reynold's numbers.