Conical wing of maximum lift-drag ratio in a supersonic gas flow

The calculation of supersonic flow past three-dimensional bodies and wings presents an extremely complicated problem, whose solution is made still more difficult in the case of a search for optimum aerodynamic shapes. These difficulties made it necessary to simplify the variational problems and to use the simplest dependences, such as, for example, the Newton formula [1–3]. But even in such a formulation it is only possible to obtain an analytic solution if there are stringent constraints on the thickness of the body, and this reduces the three-dimensional problem for the shape of a wing to a two-dimensional problem for the shape of a longitudinal profile. The use of more complicated flow models requires the restriction of the class of considered configurations. In particular, paper [4] shows that at hypersonic flight velocities a wing whose windward surface is concave can have the maximum lift-drag ratio. The problem of a V-shaped wing of maximum lift-drag ratio is also of interest in the supersonic velocity range, where the results of the linear theory of [5] or the approximate dependences of the type of [6] can be used.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 128–133, May–June, 1986.We note in conclusion that this analysis is valid for those flow regimes for which there are no internal shock waves in the shock layer near the windward side of the wing.     

Calculation of the flow of an ideal fluid past a wing of finite thickness

The problem of irrotational flow past a wing of finite thickness and finite span can be reduced by Green's formula to the solution of a system of Fredholm equations of the second kind on the surface of the wing [1]. The wake vortex sheet is represented by a free vortex surface. Besides panel methods (see, for example, [2]) there are also methods of approximate solution of this problem based on a preliminary discretization of the solution along the span of the wing in which the two-dimensional integral equations are reduced to a system of one-dimensional integral equations [1], for which numerical methods of solution have already been developed [3–6]. At the same time, a discretization is also realized for the wake vortex sheet along the span of the wing. In the present paper, this idea of numerical solution of the problem of irrotational flow past a wing of finite span is realized on the basis of an approximation of the unknown functions which is piecewise linear along the span. The wake vortex sheet is represented by vortex filaments [7] in the nonlinear problem. In the linear problem, the sheet is represented both by vortex filaments and by a vortex surface. Examples are given of an aerodynamic calculation for sweptback wings of finite thickness with a constriction, and the results of the calculation are also compared with experimental results.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 124–131, October–December, 1981.     

Experimental study of vortex flow about low aspect ratio wings and circular cones at M = 2

Characteristic of the flow about wings of low aspect ratio with subsonic leading edges and bodies of revolution at angles of attack is the formation of spiral vortices as a result of rolling-up of the transverse flow, which separates near the wing leading edge and on the lateral generator of the body. The vortices, concentrated in a pair of free vortex cores, interact with the boundary layer, causing a complex flow pattern on the surface of the model in question.There are several methods which make it possible to study the flow about the model. Pickups may be used to measure the pressure field or the velocity field near the model. This technique has found wide application and was used for studying the flow pattern about wings and bodies of revolution at both subsonic and supersonic speeds (see, for example, [1–4]). However, this method is very tedious and, in addition, the probes always introduce disturbances into the flow, particularly for supersonic speeds.A visual picture of the vortex flow may be obtained in a towing basin by adding to the water metal powder in the suspended state, or by introducing filaments of colored liquid [1, 5].The vapor screen [6] and smoke [3] methods are also used for flow visualization.The boundary layer flow on the model may be studied with the aid of oil or evaporting coatings. These methods have been used in [1, 7] to study flow about wings and in [8] to study flow about circular cones.According to the studies presented in [9] of an electric discharge with the application of high voltage to electrodes located in an air stream, a stable glow occurs as a result of the prebreakdown discharge.The properties of the prebreakdown discharge have been used by the authors of the present paper to study visually the vortex flows (high voltage electric discharge method). This technique was used to obtain the trajectories of the vortex trails for low aspect ratio wings and circular cones mounted at various angles of attack in a stream with Mach number M=2 and Reynolds number R=0.9·10⁶.In conclusion the author wishes to thank B. V. Kalachev, R. V. Bertyn, and E. D. Korolev for assistance in carrying out the experiments.     

Determination of the vorticity on a wing of small aspect ratio in a hypersonic stream

The theory of a thin shock layer [1–3] is used to obtain a formula for calculating the component of the vorticity in the direction of the flow on a wing of small aspect ratio in a hypersonic gas stream. It is shown that for definite shapes of the wing and flow regimes zones may occur with large local values of the vorticity, which, as is well known, have a significant influence on the structure of the flow field.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 175–178, September–October, 1980.     

Asymptotic solution to the problem of viscous flow in the neighborhood of the axis of a vortex sheet

An asymptotic solution is constructed to the problem of the flow of a viscous incompressible fluid in the neighborhood of the axis of a vortex sheet generated by flow separation from sharp edges of a delta wing of small aspect ratio at large values of the Reynolds number and small angles of attack.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 57–65, January–February, 1984.     

Thin weakly curved wings with maximum aerodynamic quality

Lifting wings that only slightly disturb the supersonic gas flow are considered. The plan shape and thickness distribution of the wing and the free-stream parameters are given. The flow problem is solved within the framework of the Prandtl model. The outer potential flow is determined in accordance with the linear theory. The turbulent boundary layer is found by the method of plane sections with allowance for the three-dimensional inviscid flow pattern. A numerical model of the flow is constructed in the class of piecewise-constant functions on characteristic calculation grids [1]. The variational problem of finding the weakly curved middle surface of the wing giving maximum aerodynamic quality is reduced, by analogy with [2], to a problem of nonlinear programming and is solved by the gradient projection method.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 165–168, July–August, 1991.     

On the theory of a wing with small aspect ratio in a hypersonic flow

A study is made of flow over three-dimensional wings of small aspect ratio with shape close to that of a flat delta-shaped wing. The obtained results make it possible to estimate the influence of the plan shape of the leading edge and the curvature of the wing on the pattern of the flow over its windward surface and on the corresponding gas-dynamic functions.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 112–117, July–August, 1980.     

The flow structure near the windward side of V-shaped wings with an attached shock wave on the leading edges

The results of a numerical calculation of a symmetric flow of supersonic gas with the Mach number M=3 past the windward side of V-shaped wings with an opening angle =40° and apex angles =30, 45, and 90° are given. The possibility of the ascent of one or two Ferri points from the break point of the transverse contour of the wing is discovered and explained. It is shown that conical flow near wings of finite length need not exist in flow regimes corresponding to angles of attack at which a Ferri point ascends, while at angles of attack smaller and larger than a certain interval, conical flow will exist. The investigation is conducted by means of a numerical method of stabilization with an artificial viscosity. The longitudinal coordinate, relative to which the steady system of equations is hyperbolic, played the part of the time variable, usual for methods of stabilization. The numerical method constructed using the scheme of [1] is described in [2] and was successfully applied to the calculation of different regimes of supersonic flow past conical wings with supersonic leading edges [2–6]. In. the present investigation the calculation algorithm of [2] is modified and makes it possible to realize motion with respect to the parameter a, this being particularly important for the stabilization of the solution in the calculation of flow regimes for which regions with a total velocity Mach number close to unity arise in the flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 122–131, January–February, 1986.     

Laminar boundary layer on swept-back wings of infinite span at an angle of attack

A study is made of the flow of a compressible gas in a laminar boundary layer on swept-back wings of infinite span in a supersonic gas flow at different angles of attack. The surface is assumed to be either impermeable or that gas is blown or sucked through it. For this flow and an axisymmetric flow an analytic solution to the problem is obtained in the first approximation of an integral method of successive approximation. For large values of the blowing or suction parameters, asymptotic solutions are found for the boundary layer equations. Some results of numerical solution of the problem obtained by the finite-difference method are given for wings of various shapes in a wide range of angles characterizing the amount by which the wings are swept back and also the blowing or suction parameters. A numerical solution is obtained for the equations of the three-dimensional mixing layer formed in the case of strong blowing of gas from the surface of the body. The analytic and numerical solutions are compared and the regions of applicability of the analytic expressions are estimated. On the basis of the solutions obtained in the present paper and studies of other authors a formula is proposed for the calculation of the heat fluxes to a perfectly catalytic surface of swept-back wings in a supersonic flow of dissociated and ionized air at different angles of attack. Flow over swept-back wings at zero angle of attack has been considered earlier (see, for example, [1–4]) in the theory of a laminar boundary layer. In [5], a study was made of flow over swept-back wings at nonzero angle of attack at small and moderate Reynolds numbers in the framework of the theory of a hypersonic viscous shock layer.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 27–39, May–June, 1980.We thank G. A. Tirskii for a helpful discussion of the results.     

Hypersonic flow of gas over a slender wing of variable shape

In a formulation analogous to [1–3], a study is made of the flow of a uniform homogeneous hypersonic ideal gas over the windward side of a slender wing whose surface profile depends on the time. The problem is solved by the thin shock layer method [4]. The bow shock is assumed to be attached to the leading edge of the wing at at least one point. The corrections of the first approximation to the main Newtonian flow are found. For wings of finite aspect ratio, when the bow shock is attached along the whole of the leading edge of the wing, computational formulas are obtained for determining the parameters of the gas in the shock layer.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 94–101, July–August, 1979.     

Numerical study of steady supersonic separated flow past spatial lifting systems

At the present, spatial lifting systems are usually calculated numerically using linear approximation. However, the practical application of such systems at moderate and large angles of incidence requires new approaches that allow for various nonlinear effects such as large disturbances, flow separation, and jumps in entropy across shock waves. The existing investigations [3, 4] generally cover only simple systems (bodies of revolution, wings, and so on). Here, a numerical method is proposed for investigating supersonic flows past complicated spatial systems. The method extends and continues the well-known methods widely used to solve analogous problems in subsonic aerodynamics [5, 6]. Some examples of the computation of the aerodynamic parameters for flows past wings and spatial lifting systems are also given.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 142–148, May–June, 1993.     

Computation of the two-dimensional viscous incompressible fluid flow with separation at the nlet edge around a cascade

A complex flow consisting of an outer inviscid stream, a dead-water separation domain, and a boundary layer, which interact strongly, is formed in viscous fluid flows with separation at the streamlined profile with high Re numbers. Different jet and vortex models of separation flow are known for an inviscid fluid; numerical, asymptotic, and integral methods [1–3] are used for a viscous fluid. The plane, stationary, turbulent flow through a turbine cascade by a constant-density fluid without and with separation from the inlet edge of the profile and subsequent attachment of the stream to the profile (a short, slender separation domain) is considered in this paper.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 34–44, May–June, 1978.     

Mathematical simulation of unsteady flow past a wing with jets

The various approximate approaches to the investigation of the unsteady aerodynamic characteristics of an airfoil with jet flap [1–3] are applicable only for an airfoil, low jet intensity, and low oscillation frequencies. In the present paper, the method of discrete vortices [4] is generalized to the case of unsteady flow past a wing with jets and arbitrary shape in plan. The problem is solved in the linear formulation; the conditions used are standard: no flow through the wing and jet, finite velocities at the trailing edges where there is no jet, and also a dynamical condition on the jet. The wing and jet are assumed to be thin and the medium inviscid and incompressible.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 139–144, May–June, 1982.     

Calculation of the flow around wings of arbitrary planform over a wide range of angles of attack

A numerical method is described for the calculation of the distributed and total aerodynamic characteristics of a thin wing of any planform. We use only the generally accepted hypotheses-smoothness of flow around the wing and the Chaplygin-Zhukovskii condition of finite velocity at the trailing edges. The medium is considered ideal and incompressible.The development of a nonlinear theory for the wing of small aspect ratio in a compressible medium is one of the most important and difficult problems of wing theory. It has long attracted the attention of the aerodynamiscists. Chaplygin touched on this question in his 1913 report On the vortex theory of the finite span wing, presented to the Moscow Mathematical Society. Several interesting ideas and schemes were proposed by Golubev (see, for example, [2]). The first adequately correct and effective attempt to determine theoretically the nonlinear variation of wing normal force with angle of attack was that of Bollay [3]. In this work he studied rectangular wings of very small aspect ratio. The circulation variation law along the span was taken to be constant, and along the chord it was taken the same as for a flat plate of infinite span. It was also assumed that the centerlines of the free vortices trailing from the wing tips are straight lines and form the same angle with the plane of the wing. The magnitude of this angle was calculated from the average value of the relative velocity. The boundary condition at the wing was satisfied at a single point.In several later studies [4–8] attempts have been made to extend this approach somewhat. In [7] the circulation variation law along the wing chord is calculated, and the boundary conditions are satisfied more exactly. However, attempts to convert to the study of wings of more complex planform, when the circulation can no longer be considered constant along the span, are hydrodynamically incorrect [5, 6, 8]. In these studies schemes are used in which with smooth flow around the wing the free vortices stand off from the wing surface. The angles which the vortex centerlines form with the wing surface are assumed or are calculated on the basis of very arbitrary hypotheses.In the present paper the vortex layer which simulates the wing surface, just as in the linear theory [9, 10], is replaced by a system of discrete vortices. The free vortices away from the wing then are discrete curvilinear vortex filaments. Each of them is replaced by a series of rectilinear vortex segments. The number of bound and free discrete vortices may be increased without limit. The position of the free vortex segments is determined in the computation process, which is carried out sequentially for a series of angles of attack , beginning with 0 when the linear theory scheme holds. We note that the question of accounting for the effect of the leading-edge free vortex sheet is not considered here, although the method described may also be used to obtain results for this problem.The proposed method turned out to be very general, flexible and convenient for the digital computer. It permits studying the practical convergence of the solution, and also permits obtaining not only the total and distributed characteristics of the wing of arbitrary planform, but also studying such delicate questions as the rollup of the vortex sheet behind the wing.The author wishes to thank O. N. Sokolov and T. M. Muzychenko for the example calculations.     

Some aspects of calculating subsonic flow over wings of complex planform

The use of wings of complex planform is characteristic of the present stage of development of aviation; with discontinuities along the leading and trailing edges; with curved edges; with variable geometry (by pivoting the wing panels). This article considers some aspects of the calculation of the over-all and distributed aerodynamic characteristics of such wings for low and high subsonic velocities. The methods, based on the lifting surface scheme and the use of discrete vortical singularities, enable quite efficient and reliable digital computation of the flow about these wings at moderate angles of attack. For steady motion of the wing a further development of the method of [1] is obtained, for harmonic oscillations an extension of [2] is obtained, and for aperiodic motions of the wing and gust inputs a modification of the method of [3] is found.The author wishes to thank T. M. Muzychenko and N. G. Lavrenko for carrying out the calculations of the examples.     

Bifurcation of developed flow in rectangular channels rotating about the transverse axis

The linear problem of the stability of Poiseuille flow between two infinite plates rotating about an axis parallel to the plates and normal to the flow direction was studied in [1, 2]. It was established that the flow is least stable with respect to disturbances in the form of standing waves known as Taylor eddies. The experimental data of [1, 3] and the results [4] of a numerical integration of the Navier-Stokes equations for channels with cross sections highly elongated in the direction of the axis of rotation are in good agreement with the conclusions of the linear theory. In the case of channels with a side ratio of the order of unity, which is of greater practical interest, the primary flow becomes essentially three-dimensional and evolves with variation of the governing criteria: the Reynolds and Rossby numbers. Obviously, this seriously complicates the use of the methods of the linear theory. The effect of the ratio of the sides of the cross section on the stability of the primary flow regime was studied experimentally in [3]. The present article describes the results of an investigation of the problem based on a numerical method of integrating the nonlinear Navier-Stokes equations. Moreover, an asymptotic estimate of the stability limit of the primary regime, based on a local condition of inviscid instability of rotating flows, is presented.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 27–32, September–October, 1985.     

Investigation of the influence of the transverse sweepback of delta wings on their vortex structures and aerodynamic characteristics in the case of flow separation at low subsonic velocities

The results are given of theoretical and experimental investigations into the influence of transverse sweepback of delta wings of small aspect ratio on their lifting properties and the stability of the vortex structure in the case of flow separation. Theoretical investigations have been made using the method of discrete vortices. The physical experiment included weight tests and also investigations on the breakup of vortex filaments by means of aerometric and schlieren methods. The results of the investigations are given and the calculated and experimental data are compared.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 173–175, September–October, 1982.     

Wings with optimal characteristics in hypersonic flow

We present some exact results and results obtained by the local variation method [1–3] from numerical solutions of variational problems on the symmetric wing of minimum drag and the wing of maximum L/D in hypersonic flow; a modification of the method of local variations is proposed for the numerical solution of variational problems with isoperimetric constraints. The Newton method [4] was used to calculate the pressure distribution over the wing surface.The author wishes to thank M. N. Kogan and O. G. Fridlender for many helpful discussions.     

高速小展弦比机翼颤振的微分求积法分析

引入微分求积法,分析高速小展弦比机翼的气动弹性问题。将小展弦比机翼等效为悬臂板,基于一阶活塞气动力理论建立机翼颤振偏微分方程,采用微分求积法将偏微分方程转化为常微分方程,根据频率重合理论对颤振问题进行求解。分析了机翼的固有频率及颤振速度,并与有限元软件计算结果进行比较,误差在2％以内,很好的验证了微分求积法求解小展弦比机翼颤振问题的有效性。分析了机翼面积、展弦比及厚度对颤振速度的影响,结果表明,小展弦比机翼的颤振速度受结构尺寸的影响较大,颤振速度随面积和展弦比的增大而减小,随机翼厚度的增大而增大。     

Vortex structure of separated flows on model wings at low freestream velocities

Flow past model wings is experimentally investigated in a subsonic wind tunnel at large angles of attack at which the laminar boundary layer separates near the leading edge of the wing (flow stall). The object of the study was the flow structure within the separation zone. The carbon-oil visualization of surface streamlines used in the experiments showed that in the separation zone there exist one or more pairs of large-scale vortices rotating in the wing plane. Certain general properties of the vortex structures in the separation zone are found to exist, whereas the flow patterns may differ depending on the model aspect ratio, the yaw angle, and other factors.