Flow around λ wings with a bend in the leading edge

In recent years a considerable number of studies have been published on flow around wings at high supersonic velocities. The researches have been conducted in two directions: there are studies of hypersonic flow around wings of traditional shape and a search is carried out for new types of lay-out which possess optimal aerodynamic characteristics. The second direction relates to the numerous studies of flow around wings with shaped transverse cross sections [1–7]. The calculation of the aerodynamic quality of a shaped delta wing composed of plane surfaces on the basis of the relationships on an oblique shock [1, 2], from the results of experiments on the pressure distribution and from weight tests [3, 4], showed that the shaped wing has a higher quality than the plane delta wing.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 171–175, January–February, 1985.     

Contribution to the theory of the high compression ratio gas ejector with cylindrical mixing chamber

Ways of improving the operation of a gas ejector with a high compression ratio are investigated. The conditions for obtaining the maximal compression ratio at the critical operating regime of the gas ejector are studied theoretically and experimentally with account for mixing of the supersonic injecting and subsonic ejected streams ahead of the choking section. The principles for the rational utilization of the effect of stream mixing in the ejector ahead of the choking section are indicated; the use of these principles permits a several-fold increase of the compression ratio of the supersonic ejector. A theory is given for the critical regime of the gas ejector with uniformly perforated nozzle, and the hydraulic parameters of the required wall perforationss are determined. It is shown that perforation as a hydraulic factor can improve significantly the parameters of the sonic ejector in the critical regime.The foundations of modem gas ejector theory were developed by Khristianovich [1, 2]. In these studies he established the relationship between the parameters of the flow at the end of the mixing chamber (section 3, p ₀ is the total pressure, is the reduced velocity) and the parameters of the ejecting (section 1, p ₀ ,) and the ejected (p₀₁,₁) flows with account for compressibility for the ejector with a cylindrical mixing chamber (Fig. 1a). The ejector theory [1, 2] (see also [3, 4]) is given in the hydraulic approximation: the flow at the end of the mixing chamber is assumed uniform, flow friction on the mixing chamber walls is neglected. The use of the gasdynamic functions [5–9] made it possible to obtain computational equations for the ejector in a convenient form and to extend them to the case of mixing of gases with different thermophysical properties. We note that for subsonic velocities of the ejecting and ejected flows the system of ejector equations [1, 2] is supplemented by the condition of equality of the static pressures p=P₁ at the stream contact section 1.The results of extensive experimental studies of subsonic ejectors are in good agreement with the results of this theory.For sonic or supersonic velocity of the ejecting gas (=1) the condition p=p₁ is not satisfied in the general case. Fundamental for the development of ejector theory was the establishment by Millionshchikov and Ryabinkov in 1948 of the existence of a critical operating regime of the supersonic ejector [7, 10]. They showed that the limiting operating regimes of the gas ejector for high pressure differentials ==p ₀ /p₀₁ are determined by the conditions for the choking of the ejected jet by the expanding supersonic ejecting flow. With the occurrence of the critical regime the velocity of the ejected jet at the choking section (section 2, Fig. 1a) reaches the speed of sound (=1); this limits the further increase of the pressure ratio and the ejector compression ratio =p ₀ /p ₀ for a given ejection coefficient k (k is the ratio of the ejected and ejecting gas flow rates). The relationships between these flow parameters at sections 2 and 1 supplement the system of ejector equations and permit determining its critical characteristics.Millionshchikov and Ryabinkov showed that for moderate values of the pressure ratio good agreement of the theoretical and experimental ejector characteristics are given by the assumption of constant static pressure p₂=const at section 2 (Fig. 1a).The limit of the applicability of the theory based on the condition p₂= = const, was studied experimentally by Lyzhin [10].The theory of the critical regime of the gas ejector was developed in 1953 in studies of Nikol'skii, Shustov, Vasil'ev, Taganov, and Mezhirov [10, 11]. Nikol'skii showed that the condition of constant static pressure at the choking section is not in agreement with the momentum equation.For a more rigorous theoretical determination of the critical ejector regime he proposed joining between sections 1 and 2 (Fig. 1a) the calculation of the ejecting jet using the method of characteristics and the hydraulic calculation of the ejected jet; example calculations were made by Nikol'skii and Shustov. Taganov and Mezhirov suggested a method for calculating the ejector critical regime using a linear distribution of the pressure in the supersonic ejecting jet (at the choking section 2).A simple and successful method for calculating the ejector critical regime was given by Vasil'ev, who used the hydraulic representation of the ejecting and ejected flows in the choking section; both flows are assumed uniform at section 2, the static pressures in these flows in the general case are different and are determined by the momentum equation. A similar theory for the ejector critical regime was developed independently in [12, 13], and the theory with account for the supersonic ejecting flow (ahead of the choking section) was developed using the method of characteristics in [14].It should be noted that the results of the calculations of the critical characteristics of the ejectors using all three of these methods were practically indentical and in good agreement with experiment for large and moderate values of the ejection coefficients. We emphasize that in the theories of the ejector critical regime the flow mixing between sections 1 and 2 is neglected.The critical regime theory imposes significant limitations on the possible characteristics of the gas ejector, first of all, on the achievable compression ratio =p ₀ /p ₀ . Thus, from the data of [10], even for a pressure ratio =1000 the maximal theoretical value of the compression ratio for the supersonic ejector does not exceed 40 (see in Fig. 2 the limiting ejector characteristics based on the critical regime theory); for the sonic air ejector (=1) the theoretical value of 3.5 (see Fig. 9b on p. 26). Therefore it is important to analyze the methods for influencing the critical regime parameters in order to determine ways to improve the operation of the gas ejector with a high compression ratio.     

Method of calculating the longitudinal,lateral, and crossed aerodynamic derivatives of an aircraft at subsonic velocities

It is impossible to take account of the numerous, sometimes contradictory, requirements in the design of modern aircraft without carrying out broad parametric investigations. Therefore, there is a need for effective and sufficiently exact methods of computer calculation, making it possible to carry out a large volume of investigations in a short period of time, even in the design stage of aircraft. The present article sets forth a method for calculating the longitudinal, lateral, and crossed aerodynamic derivatives and the induced resistance, taking account of the thickness of the carrying elements and suspensions of an aircraft moving at a subsonic velocity and performing harmonic vibrations with a small frequency. The method allows of a sufficiently exact calculation of the aerodynamic characteristics of an aircraft of any given complex configuration. A limitation is the development of critical phenomena (breakaway of the flow and a zone of supersonic flows).Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 76–88, March–April, 1978.The author thanks N. N. Tyunin for kindly furnishing the experimental data on m _x ^y + m _x ^* .     

Indentation of the semi-infinite elastic solid by a concave rigid punch

A solution is obtained for the relationship between load, displacement and inner contact radius for an axisymmetric, spherically concave, rigid punch, indenting an elastic half-space. Analytic approximations are developed for the limiting cases in which the ratio of the inner and outer radii of the annular contact region is respectively small and close to unity. These approximations overlap well at intermediate values. The same method is applied to the conically concave punch and to a punch with a central hole. , , . , . . .     

Supersonic flow past axisymmetric flared cones

Systematic data on the determination of the aerodynamic characteristics of axisymmetric bodies with a break in the generating line (Fig. 1a, b) in supersonic flow at zero angle of attack are presented in [1, 2, and others]. A characteristic feature of the flow past such bodies is the appearance of an extensive separation zone dec in the region of the break in the generator when the break angle exceeds some minimum value _min, which for a turbulent boundary layer depends basically on the Mach number M at the body surface ahead of the separation zone. In this case, compression waves which change into the oblique compression shocks dc and cc, emanate both from the beginning of the separation zone (point c) and from the end of it (point d). These shocks, intersecting at the point c, form the triple shock configuration acd and acc for which we introduce the notationac[c, d]. The maximum value (_max) of the generator break angle is limited by the possibility of the existence of an attached compression shock, dc. According to these data a change in the generator break angle for the range _minmax of the angle does not disrupt the nature of the flow in the separation zone, but only alters the size of this zone.We shall examine the flow past cones with values of the generator break angles (_max) for which the attached shock dc cannot exist.     

Effect of angle of attack on supersonic flow past axisymmetric blunt bodies with surface injection

The effect of angles of attack in the interval 0 40° on the flow pattern and the aerodynamic characteristics of a body of power-law shape (equation of the generator in the cylindrical coordinate system r=zⁿ, n=0.125) is investigated for supersonic flow without injection and with intense subsonic localized injection from the surface. As a result of numerical calculations it is established that the use of Newton's theory for determining the coordinates of the gas stagnation point behind the shock in flow past an impermeable body of the shape in question leads to serious errors, and an expression for determining the location of this point is given. It is shown that for three-dimensional flow the flow pattern and the surface pressure distribution are sharply different from the case=0. It is established that on the parameter interval in question intense injection considerably reduces the aerodynamic drag without loss of static stability, which is important in connection with the solution of the problem of gas-stable aircraft control.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 95–101, September–October, 1987.     

General expressions for unsteady forces in a lattice of profiles

A large number of studies have been devoted to the unsteady flow of a viscid incompressible fluid past a lattice of thin profiles and the determination of the resulting aerodynamic forces and moments. For example, in the particular case of the motion of a lattice with stagger with zero phase shift of the oscillations between neighboring profiles, Haskind [1] determined the unsteady lift force and moment. Popescu [2] suggested expressions for the force and moment in the case when =0 and =0, using the method of conformal mapping. Samoilovich [3] obtained equations for the unsteady lift force and moment by the method of the acceleration potential for phase shift =0 and = of the oscillations between neighboring profiles. Musatov [4] used an electronic digital computer to calculate the overall unsteady aerodynamic characteristics of a grid by the vortex method, taking into account the amplitude of the oscillations and the initial circulation for =m (m1). Gorelov [5] determined the coefficients of the over-all unsteady aerodynamic force and moment of a profile in a lattice with the stagger and any value of =m. He used a method based on the unsteady flow past an isolated profile with subsequent account for the interference of the profiles in the lattice.In the following we find general expressions for the unsteady lift force and moment acting on a lattice moving in an incompressible fluid with the constant velocity U. These formulas generalize the known formulas for the isolated profile [6]. The profiles of a staggered grid (Section 1) are considered to be thin and slightly curved, and perform oscillations with a phase shift of the oscillations between neighboring profiles. The method of separation of singularities is used to obtain the solution in closed form. The coefficients of the expansion of the complex velocity in a series in the derivatives of a function are calculated. An integral equation relative to the unknown tangential velocity component in the wake is derived (Section 2), and its analytic solution is given (Section 3). For =0 the solution coincides with the solution obtained earlier in [7]. Expressions are obtained for the forces and moments (Section 4) in the form of four terms. The first two terms determine the force and moment for motion with constant circulation, and the last two determine these characteristics for motion with variable circulation. The suction force acting at the leading edges of the profiles is found in a general form. Particular cases of closely and widely spaced lattices are considered. Computational results are presented.     

On the Ellipticity of the Middle Surface of a Shell and its Application to the Asymptotic Analysis of ‘Membrane’ Shells

We consider a surface S = (), where ² is a bounded, connected, open set with a smooth boundary and : ³ is a smooth map; let () denote the components of the two-dimensional linearized strain tensor of S and let 0 with length 0 > 0. We assume the the norm ,|| ()||0, in the space V0() = { H¹() × H¹() × L²(); = 0 on 0 } is equivalent to the usual product norm on this space. We then establish that this assumption implies that the surface S is uniformly elliptic and that we necessarily have 0 = .     

Study of supersonic flow past axisymmetric bodies of different form

A method is proposed in [1] for calculating supersonic flow past smooth bodies. The present article presents a computational scheme and calculational formulas for determining the gasdynamic functions at the nodes of the lines=const. A comparison is made of certain of the results obtained with the results of other studies [2, 3]. Results are also presented of the calculation of the flow of a perfect gas past ellipsoids of revolution (=2.3) and inverted cones with spherical and ellipsoidal blunting.The author wishes to thank G. F. Telenin for his guidance in the present study.     

Thin-Walled Beams: the Case of the Rectangular Cross-Section

In this paper we present an asymptotic analysis of the three-dimensional problem for a thin linearly elastic cantilever =×(0,l) with rectangular cross-section of sides and ², as goes to zero. Under suitable assumptions on the given loads, we show that the three-dimensional problem converges in a variational sense to the classical one-dimensional model for extension, flexure and torsion of thin-walled beams. Mathematics Subject Classifications (2000) 474K20, 74B10, 49J45.     

Hypersonic flow past blunt-edged delta wings

A method is proposed for calculating hypersonic ideal-gas flow past blunt-edged delta wings with aspect ratios = 100–200. Systematic wing flow calculations are carried out on the intervals 6 M 20, 0 20, 60 80; the results are analyzed in terms of hypersonic similarity parameters.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 175–179, September–October, 1990.     

Gasdynamics of a plasma in a multiple-mirror magnetic trap with “point” mirrors

Equations are derived for the gasdynamics of a dense plasma confined by a multiple-mirror magnetic field. The limiting cases of large and small mean free paths have been analyzed earlier: ₀ and k, where is the length of an individual mirror machine, ₀ is the size of the mirror, and k is the mirror ratio. The present work is devoted to a study of the intermediate range of mean free paths ₀ k. It is shown that in this region of the parameters the process of expansion of the plasma has a diffusional nature, and the coefficients of transfer of the plasma along the magnetic field are calculated.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 14–19, November–December, 1974.The authors thank D. D. Ryutov for the statement of the problem and interest in the work.     

Heteroclinic orbits in singular systems: A unifying approach

We consider singularly perturbed systems , such that=f(, _o, 0). _o ^m , has a heteroclinic orbitu(t). We construct a bifurcation functionG(, ) such that the singular system has a heteroclinic orbit if and only ifG(, )=0 has a solution=(). We also apply this result to recover some theorems that have been proved using different approaches.     

Concentration dependence of the critical viscoelastic properties of gelatin at the gel point

Gelatin gel properties have been studied through the evolution of the storage [G()] and the loss [G()] moduli during gelation or melting near the gel point at several concentrations. The linear viscoelastic properties at the percolation threshold follow a power-law G()G() and correspond to the behavior described by a rheological constitutive equation known as the Gel Equation. The critical point is characterized by the relation: tan = G/G = cst = tan ( · /2) and it may be precisely located using the variations of tan versus the gelation or melting parameter (time or temperature) at several frequencies. The effect of concentration and of time-temperature gel history on its variations has been studied. On gelation, critical temperatures at each concentration were extrapolated to infinite gel times. On melting, critical temperatures were determined by heating step by step after a controlled period of aging. Phase diagrams [T = f(C)] were obtained for gelation and melting and the corresponding enthalpies were calculated using the Ferry-Eldridge relation. A detailed study of the variations of A with concentration and with gel history was carried out. The values of which were generally in the 0.60–0.72 range but could be as low as 0.20–0.30 in some experimental conditions, were compared with published and theoretical values.     

Methods of Small and Large δ in the Nonlinear Dynamics – A Comparative Analysis

New asymptotic approaches for dynamical systems containing a power nonlinear term x ⁿ are proposed and analyzed. Two natural limiting cases are studied: n 1 + , 1 and n . In the firstcase, the 'small method' (SM)is used and its applicability for dynamical problems with the nonlinearterm sin as well as the usefulness of the SMfor the problem with small denominators are outlined. For n , a new asymptotic approach is proposed(conditionally we call it the 'large method' –LM). Error estimations lead to the followingconclusion: the LM may be used, even for smalln, whereas the SM has a narrow application area. Both of the discussed approaches overlap all values ofthe parameter n.     

Heat transfer and equilibrium temperature of a sphere in a supersonic rarefied gas flow

Some results are presented of experimental studies of the equilibrium temperature and heat transfer of a sphere in a supersonic rarefied air flow.The notations D sphere diameter - u, , T,,l, freestream parameters (u is velocity, density, T the thermodynamic temperature,l the molecular mean free path, the viscosity coefficient, the thermal conductivity) - T₀ temperature of the adiabatically stagnated stream - T_e mean equilibrium temperature of the sphere - T_w surface temperature of the cold sphere (T_we) - mean heat transfer coefficient - _e air thermal conductivity at the temperature T_e - P Prandtl number - M Mach number     

The extension of Beltrami's ellipsoid to anisotropic bodies

Summary The spectral decomposition of the compliance, stiffness, and failure tensors for transversely isotropic materials was studied and their characteristic values were calculated using the components of these fourth-rank tensors in a Cartesian frame defining the principal material directions. The spectrally decomposed compliance and stiffness or failure tensors for a transversely isotropic body (fiber-reinforced composite), and the eigenvalues derived from them define in a simple and efficient way the respective elastic eigenstates of the loading of the material. It has been shown that, for the general orthotropic or transversely isotropic body, these eigenstates consist of two double components, ₁ and ₂ which are shears (₂ being a simple shear and ₁, a superposition of simple and pure shears), and that they are associated with distortional components of energy. The remaining two eigenstates, with stress components ₃, and ₄, are the orthogonal supplements to the shear subspace of ₁ and ₂ and consist of an equilateral stress in the plane of isotropy, on which is superimposed a prescribed tension or compression along the symmetry axis of the material. The relationship between these superimposed loading modes is governed by another eigenquantity, the eigenangle .The spectral type of decomposition of the elastic stiffness or compliance tensors in elementary fourth-rank tensors thus serves as a means for the energy-orthogonal decomposition of the energy function. The advantage of this type of decomposition is that the elementary idempotent tensors to which the fourth-rank tensors are decomposed have the interesting property of defining energy-orthogonal stress states. That is, the stress-idempotent tensors are mutually orthogonal and at the same time collinear with their respective strain tensors, and therefore correspond to energy-orthogonal stress states, which are therefore independent of each other. Since the failure tensor is the limiting case for the respective _x, which are eigenstates of the compliance tensor S, this tensor also possesses the same remarkable property.An interesting geometric interpretation arises for the energy-orthogonal stress states if we consider the projections of _x in the principal₃D stress space. Then, the characteristic state ₂ vanishes, whereas stress states ₁, ₃ and ₄ are represented by three mutually orthogonal vectors, oriented as follows: The ₃ and ₄ lie on the principal diagonal plane (₃₁₂) with subtending angles equaling (–/2) and (-), respectively. On the positive principal ₃-axis, is the eigenangle of the orthotropic material, whereas the ₁-vector is normal to the (₃₁₂)-plane and lies on the deviatoric -plane. Vector ₂ is equal to zero.It was additionally conclusively proved that the four eigenvalues of the compliance, stiffness, and failure tensors for a transversely isotropic body, together with value of the eigenangle , constitute the five necessary and simplest parameters with which invariantly to describe either the elastic or the failure behavior of the body. The expressions for the _x-vector thus established represent an ellipsoid centered at the origin of the Cartesian frame, whose principal axes are the directions of the ₁-, ₃- and ₄-vectors. This ellipsoid is a generalization of the Beltrami ellipsoid for isotropic materials.Furthermore, in combination with extensive experimental evidence, this theory indicates that the eigenangle alone monoparametrically characterizes the degree of anisotropy for each transversely isotropic material. Thus, while the angle for isotropic materials is always equal to _i = 125.26° and constitutes a minimum, the angle || progressively increases within the interval 90–180° as the anisotropy of the material is increased. The anisotropy of the various materials, exemplified by their ratiosE _L/2G_L of the longitudinal elastic modulus to the double of the longitudinal shear modulus, increases rapidly tending asymptotically to very high values as the angle approaches its limits of 90 or 180°.     

Analysis of turbulent boundary layer interaction with an external supersonic flow behind a step

An integral method of analyzing turbulent flow behind plane and axisymmetric steps is proposed, which will permit calculation of the pressure distribution, the displacement thickness, the momentum-loss thickness, and the friction in the zone of boundary layer interaction with an external ideal flow. The characteristics of an incompressible turbulent equilibrium boundary layer are used to analyze the flow behind the step, and the parameters of the compressible boundary layer flow are connected with the parameters of the incompressible boundary layer flow by using the Cowles-Crocco transformation.A large number of theoretical and experimental papers devoted to this topic can be mentioned. Let us consider just two [1, 2], which are similar to the method proposed herein, wherein the parameter distribution of the flow of a plane nearby turbulent wake is analyzed. The flow behind the body in these papers is separated into a zone of isobaric flow and a zone of boundary layer interaction with an external ideal flow. The jet boundary layer in the interaction zone is analyzed by the method of integral relations.The flow behind plane and axisymmetric steps is analyzed on the basis of a scheme of boundary layer interaction with an external ideal supersonic stream. The results of the analysis by the method proposed are compared with known experimental data.Notation x, y longitudinal and transverse coordinates - X, Y transformed longitudinal and transverse coordinates - , ^*, ^** boundary layer thickness, displacement thickness, momentum-loss thickness of a boundary layer - , ^*, ^** layer thickness, displacement thickness, momentum-loss thickness of an incompressible boundary layer - u, velocity and density of a compressible boundary layer - U, velocity and density of the incompressible boundary layer - , stream function of the compressible and incompressible boundary layers - , dynamic coefficient of viscosity of the compressible and incompressible boundary layers - r₁ radius of the base part of an axisymmetric body - r radius - R transformed radius - M Mach number - friction stress - p pressure - a speed of sound - s enthalpy - v Prandtl-Mayer angle - P Prandtl number - P_t turbulent Prandtl number - r₂ radius of the base sting - b step depth - =0 for plane flow - =1 for axisymmetric flow Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 33–40, May–June, 1971.In conclusion, the authors are grateful to M. Ya. Yudelovich and E. N. Bondarev for useful comments and discussions.     

Flow through a helically coiled annulus

In this paper the flow is studied of an incompressible viscous fluid through a helically coiled annulus, the torsion of its centre line taken into account. It has been shown that the torsion affects the secondary flow and contributes to the azimuthal component of velocity around the centre line. The symmetry of the secondary flow streamlines in the absence of torsion, is destroyed in its presence. Some stream lines penetrate from the upper half to the lower half, and if is further increased, a complete circulation around the centre line is obtained at low values of for all Reynolds numbers for which the analysis of this paper is valid, being the ratio of the torsion of the centre line to its curvature.Nomenclature A =constant - a outer radius of the annulus - b unit binormal vector to C - C helical centre line of the pipe - D rL - g 1000 - K Dean number=Re² - L 1+r sin - M (L ²+ ² r ²)^1/2 - n unit normal vector to C - P, P pressure and nondimensional pressure - p ₀, p pressures of O(1) and O() - Re Reynolds number=aW ₀/ - (r, , s), (r, , s) coordinates and nondimensional coordinates - nonorthogonal unit vectors along the coordinate directions - r ₀ radius of the projection of C - t unit tangent vector to C - V _r, V , V _s velocity components along the nonorthogonal directions - V_r, V, V _s nondimensional velocity components along - W ₀ average velocity in a straight annulus Greek symbols , curvature and nondimensional curvature of C - U, V, W lowest order terms for small in the velocity components along the orthogonal directions t - _r, , _s first approximations to V _r , V, V _s for small - =/=/ - kinematic viscosity - density of the fluid - , torsion and nondimensional torsion of C - , stream function and nondimensional stream function - nondimensional streamfunction for U, V - a inner radius of the annulus After this paper was accepted for publication, a paper entitled On the low-Reynolds number flow in a helical pipe, by C.Y. Wang, has appeared in J. Fluid. Mech., Vol 108, 1981, pp. 185–194. The results in Wangs paper are particular cases of this paper for =0, and are also contained in [9].     

Solutions for the anisotropic elastic wedge at critical wedge angles

The classical solution for an isotropic elastic wedge loaded by uniform tractions on the sides of the wedge becomes infinite everywhere in the wedge when the wedge angle 2 equals , 2 or 2_* where tan 2_* = 2_*. When the wedge is loaded by a concentrated couple at the wedge apex the solution also becomes infinite at 2 = 2_*. A similar situation occurs when the wedge is anisotropic except that 2_* is governed by a different equation and depends on material properties. Solutions which do not become infinite everywhere in the wedge are available for isotropic elastic wedges. In this paper we present solutions for the anisotropic elastic wedge at critical wedge angles. The main feature of the solutions obtained here is that they are in a real form even though Stroh's complex formalism is employed.