Effect of Porous Coatings on Stability of Hypersonic Boundary Layers

The influence of ultrasoundabsorbing coatings on stability of hypersonic boundary layers is considered. Two types of coatings were used in experiments: feltmetal with a random porous microstructure and a sheet perforated by blind cylindrical microchannels. The experiments were performed in a wind tunnel at a Mach number M = 5.95 on sharp cones with a 7° apex halfangle. Evolution of natural disturbances and artificially induced wave packets in the boundary layer was studied with the help of hotwire anemometry. Spatial characteristics of artificial disturbances were obtained. It is demonstrated that such coatings exert a stabilizing effect on secondmode disturbances.     

On the measurement of lateral velocity derivatives in turbulent flows

Direct numerical simulation data for the lateral velocity derivative u/y at the centreline of a fully developed turbulent channel flow provide reasonable support for Wyngaard's analysis of the error involved in measuring this quantity using parallel hot wires. Numerical data in the wall region of the channel flow also provide a useful indication of how to select the separation between the wires. Justification for this choice is obtained by comparing several measured statistics of u/y with the corresponding numerical data.     

Chaos and integrability in a nonlinear wave equation

We consider the parametrized family of equations _tt ,u- _xx u-au+u ₂ ² u=O,x(0,L), with Dirichlet boundary conditions. This equation has finite-dimensional invariant manifolds of solutions. Studying the reduced equation to a four-dimensional manifold, we prove the existence of transversal homoclinic orbits to periodic solutions and of invariant sets with chaotic dynamics, provided that =2, 3, 4,.... For =1 we prove the existence of infinitely many first integrals pairwise in involution.     

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.     

Approximate method for calculating turbulent boundary layer with positive pressure gradient

A single-parameter integral method is proposed for calculating the turbulent boundary layer with positive pressure gradient which makes it possible to calculate the friction, thermal flux, and layer thickness both ahead of the separation point and in some region behind the separation point.Notation u velocity - density - ^* displacement thickness - ^** momentum thickness - energy thickness - M Mach number - r radius - dynamic viscosity - c_p specific heat at constant pressure - Reynolds number based on initial boundary layer thickness - P Prandtl number - p₁ static pressure at point of initial interaction - p₂ static pressure at pressureplateau - p₀ stagnation pressure - T₀ stagnation temperature - I enthalpy - T_e recovery temperature - T_w ⁰ temperature factor - H form parameter - r₁ recovery coefficient Indices 0 denotes initial section of boundary layer - 1 parameters taken at edge of boundary layer - w parameters taken at the wall temperature - * parameters referred to flow on a flat plate with =0     

On the evolutions of the discontinuities of orders ⩾ m in solutions for quasi-linear systems of PDEs of order m,or for the linearizations of these systems

Summary One considers a system L[u]=0 of PDEs, quasi-linear (according to [1]) and of order m, which possesses a bicharacteristic line , as it happens in the hyperbolic case. For v=0, , –m (>0) let u(^v) be a discontinuity wave of order m+v that solves the system above and whose discontinuity hypersurface includes . The corresponding transport equations along are considered. Furthermore some interesting cases are pointed out, in which these equations turn out to be mutually equivalent in a suitable sense. Some theorems are stated to compare the transport equations for the discontinuities of the above kinds, that are connected with the systems d^hL[u]/dt^h=0 (h=0, , –m) and/or the linearization of the system L[u]=0 around any regular solution of it.

---

Sommario Si considera un sistema L[u]=0 di equazioni alle derivate parziali, quasi lineare (secondo [1]) e di ordine m, il quale sia dotato di qualche bicaratteristica , come accade nel caso iperbolico. Per v=0, , –m(>0) sia u(^v) un'onda di discontinuità di ordine m+v risolvente il detto sistema e avente ipersuperficie di discontinuità contenente Si considerano le relative equazioni di trasporto lungo e si determinano casi interessanti in cui queste equazioni sono mutuamente equivalenti in senso opportuno. Si stabiliscono teoremi di confronto per il trasporto delle discontinuità del tipo suddetto, relative ai sistemi d^hL[u]/dt^h=0 (h=0, , –m) e/o alla linearizazione del sistema L[u]=0 attorno a qualche sua soluzione regolare.

     

Calculation of the pressure on the side surface of bodies consisting of a spherical segment joined to an inverted cone in a supersonic stream

The results of investigations of inviscid flow over inverted cones with nose consisting of a spherical segment were published for the first time in Soviet literature in [1–4]. In the present paper, a numerical solution to this problem is obtained using the improved algorithms of [5, 6], which have proved themselves well in problems of exterior flow over surfaces with positive angles of inclination to the oncoming flow. It is shown that the Mach number 2 M , equilibrium and nonequilibrium physicochemical transformations in air (H = 60 km, V = 7.4 km/sec, R₀ = 1 m), and the angle of attack 0 40° influence the investigated pressure distributions. A comparison of the results of the calculations with drainage experiments for M = 6, = 0-25° confirms the extended region of applicability of the developed numerical methods. Also proposed is a simple correlation of the dependence on the Mach number in the range 1.5 M of the shape of the shock wave near a sphere in a stream of ideal gas with adiabatic exponent = 1.4.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 178–183, January–February, 1981.     

Free convection from a horizontal surface with mass diffusion

The flow due to free convection and mass diffusion of the steady laminar flow of an incompressible fluid on upper side of a hot horizontal surface has been studied. Similar solutions have been obtained for the boundary layer equations under certain restricting assumptions. The effects of mass diffusion on the free convection flow field have been analyzed. Numerical solutions have been obtained for some particular cases of gases with Prandtl number and Schmidt number close to unity.

---

Zusammenfassung Es wird die freie Konvektionsströmung und die Diffusion in der inkompressiblen laminaren und stationären Grenzschicht auf der Oberseite einer geheizten Platte behandelt. Die Grenzschichtgleichungen für dieses Problem wurden für diejenigen Fälle gelöst, die zu ähnlichen Lösungen führen. Der Einfluß der Diffusion auf die Strömung wurde untersucht. Numerische Beispiele für einige Fälle mit Prandtl-Zahl und Schmidt-Zahl nahe Eins sind berechnet und angegeben worden.

---

Nomenclature b Width of plate perpendicular to the plane of paper - C Concentration of the diffusing species in the main gas - C ₁, C₂ Equation constants, defined in text - C _p Specific heat at constant pressure - D Coefficient of mass diffusion - f _w Flow parameter, dimensionless - h Heat transfer coefficient - l Length of plate from which heat transfer occurs - m ₁ Molecular weight of the ambient gas - m ₂ Molecular weight of the diffusing gas - p Fluid pressure defined in text - u, v Fluid velocity components inx andy-directions respectively. - T Absolute temperature at any point in the free convection field - T ₀ Absolute temperature of the ambient atmosphere - Gr Grashof number,(l ³g T/v²) - Gr_m Diffusion Grashof number, - Gr_TM Gr + Gr_M - Nu Average nusselt number, - Pr Prandtl number,/ - Sc Schmidt number,/D h - Thermal diffusivity, - Coefficient of cubical expansion - Dynamic viscosity - Fluid density - Kinematic viscosity - ()_w Value of quantity at the horizontal surface - ()(x) Local value of the quantity - ()_t Refers to diffusing fluid     

Two-dimensional transonic vortex flows of an ideal gas

Equations are obtained for two-dimensional transonic adiabatic (nonisoenergetic and nonisoentropic) vortex flows of an ideal gas, using the natural coordinates (=const is the family of streamlines, and =const is the family of lines orthogonal to them). It is not required that the transonic gas flow be close to a uniform sonic flow (the derivation is given without estimates). Solutions are found for equations describing vortex flows inside a Laval nozzle and near the sonic boundary of a free stream.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 105–109, September–October, 1973.     

An application of non-singularity BEM to plate bending, buckling and natural vibration problems

In this paper the fundamental solution of the singular governing equation of plate static bending is taken as the Green's function, which can satisfy the governing equation precisely in the plate region. Based on the principle of superposition, let the function values on the plate boundary, induced by a set of the Green's function sources (including the known sources in the plate region and the unknown sources in the fictitious region), satisfy the prescribed conditions on specially chosen boundary matching points, and the corresponding semi-analytical and semi-numerical solution can be obtained, which is free from the restraint of boundary forms and boundary conditions. The more matching points there are on the boundary, the better the accuracy of results is. Finally, in static bending problems a set of linear algebraic equations has to be computed; in buckling problems the minimum value of buckling eigenvalue equation has to be found; in natural vibration problems the eigenvalues of the frequency equation have to be calculated. Numerical examples are given and the results are compared with those by the analytical method and other methods. It can be seen that they are very close to each other.     

Nonlinear analysis of the forced response of a beam with three mode interaction

An analysis is presented for the primary resonance of a clamped-hinged beam, which occurs when the frequency of excitation is near one of the natural frequencies,_n . Three mode interaction (₂ 3₁ and _{3 1} + 2₂) is considered and its influence on the response is studied. The case of two mode interaction (₂ 3₁) is also considered to compare it with the case of three mode interaction. The straight beam experiencing mid-plane stretching is governed by a nonlinear partial differential equation. By using Galerkin's method the governing equation is reduced to a system of nonautonomous ordinary differential equations. The method of multiple scales is applied to solve the system. Steady-state responses and their stability are examined. Results of numerical investigations show that there exists no significant difference between both modal interactions' influences on the responses.     

Hopf bifurcation in a three-dimensional system

In this paper,Liapunor-Schmidl reduction and singularity theory are employed to discuss Hopf and degenerate Hopf bifureations in global parametric region in a three-dimensional system x=-βx+y, y=-x-βy(1-kz), z=β[α(1-z)-ky²], The conditions on existence and stability are given.     

Nonlinear rheological behavior of a concentrated spherical silica suspension

Linear and nonlinear viscoelastic properties were examined for a 50 wt% suspension of spherical silica particles (with radius of 40 nm) in a viscous medium, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The effective volume fraction of the particles evaluated from zero-shear viscosities of the suspension and medium was 0.53. At a quiescent state the particles had a liquid-like, isotropic spatial distribution in the medium. Dynamic moduli G^* obtained for small oscillatory strain (in the linear viscoelastic regime) exhibited a relaxation process that reflected the equilibrium Brownian motion of those particles. In the stress relaxation experiments, the linear relaxation modulus G(t) was obtained for small step strain (0.2) while the nonlinear relaxation modulus G(t, ) characterizing strong stress damping behavior was obtained for large (>0.2). G(t, ) obeyed the time-strain separability at long time scales, and the damping function h() (–G(t, )/G(t)) was determined. Steady flow measurements revealed shear-thinning of the steady state viscosity () for small shear rates (< ^–1; = linear viscoelastic relaxation time) and shear-thickening for larger (>^–1). Corresponding changes were observed also for the viscosity growth and decay functions on start up and cessation of flow, + (t, ) and _– (t, ). In the shear-thinning regime, the and dependence of ₊(t,) and _–(t,) as well as the dependence of () were well described by a BKZ-type constitutive equation using the G(t) and h() data. On the other hand, this equation completely failed in describing the behavior in the shear-thickening regime. These applicabilities of the BKZ equation were utilized to discuss the shearthinning and shear-thickening mechanisms in relation to shear effects on the structure (spatial distribution) and motion of the suspended particles.Dedicated to the memory of Prof. Dale S. Parson     

Shock Waves in a Fluid with Bubbles Containing Evaporating Drops of a Liquefied Gas

The propagation and reflection of one-dimensional plane unsteady waves and pulses in a mixture of a fluid with two-phase bubbles containing evaporating drops is investigated. A significant effect of unsteady evaporation of the drops in the zone ahead of the shock wave on the wave propagation is demonstrated. The evaporation of the drops results in a pressure increase ahead of the wave and the shock wave as it were climbs to increasing pressure level. In contrast to bubbly fluids with single-phase bubbles, in a fluid with two-phase bubbles, at a fixed phase volume fraction, a decrease in bubble size results in an increase rather than a decrease of the oscillation amplitude. The wave reflection from a solid wall is essentially nonlinear and the maximum pressure attained at the wall is several times greater than the incident-wave intensity.     

Lp-estimates for the nonlinear spatially homogeneous Boltzmann equation

This paper studies L^p-estimates for solutions of the nonlinear, spatially homogeneous Boltzmann equation. The molecular forces considered include inverse k^th-power forces with k > 5 and angular cut-off.The main conclusions are the following. Let f be the unique solution of the Boltzmann equation with f(v,t)(1 + ¦v²¦)^(s ₁ ^{+ /p)/2} L¹, when the initial value f ₀ satisfies f ₀(v) 0, f ₀(v) (1 + ¦v¦²)^(s ₁ ^{+ /p)/2} L¹, for some s₁ 2 + /p, and f ₀(v) (1 + ¦v¦²)^s/2 L^p. If s 2/p and 1 < p < , then f(v, t)(1 + ¦v¦²)^{(s s} ₁)^/2 L^p, t > 0. If s >2 and 3/(1+ ) < p < , thenf(v,t) (1 + ¦v¦²)^(s(s ₁ ^{+ 3/p))/2} L^p, t > 0. If s >2 + 2C₀/C₁ and 3/(l + ) < p < , then f(v,t)(1 + ¦v¦²)^s/2 L^p, t > 0. Here 1/p + 1/p = 1, x y = min (x, y), and C₀, C₁, 0 < 1, are positive constants related to the molecular forces under consideration; = (k – 5)/ (k – 1) for k^th-power forces.Some weaker conclusions follow when 1 < p 3/ (1 + ).In the proofs some previously known L-estimates are extended. The results for L^p, 1 < p < , are based on these L-estimates coupled with nonlinear interpolation.     

Motion of large gas bubbles ascending in a liquid

An equation is derived for the ascent velocity of large gas bubbles in a liquid. This velocity is assumed to be governed by the propagation of a wavelike perturbation caused by the bubble in the liquid.Notation w bubble (or drop) velocity - specific gravity - dynamic viscosity - kinematic viscosity - r bubble (or drop) radius - surface tension - coefficient of friction - g gravitational acceleration - D bubble (or drop) diameter - p pressure - c propagation velocity of the wavelike perturbation - wavelength     

Finite difference analysis of plane Poiseuille and Couette flow developments

Summary The problem of flow development from an initially flat velocity profile in the plane Poiseuille and Couette flow geometry is investigated for a viscous fluid. The basic governing momentum and continuity equations are expressed in finite difference form and solved numerically on a high speed digital computer for a mesh network superimposed on the flow field. Results are obtained for the variations of velocity, pressure and resistance coefficient throughout the development region. A characteristic development length is defined and evaluated for both types of flow.Nomenclature h width of channel - L ratio of development length to channel width - p fluid pressure - p ₀ pressure at channel mouth - P dimensionless pressure, p/ ² - P ₀ dimensionless pressure at channel mouth - P pressure defect, P ₀–P - (P)₀ pressure defect neglecting inertia - Re Reynolds number, uh/ - u fluid velocity in x-direction - mean u velocity across channel - u ₀ wall velocity - U dimensionles u velocity u/ - U _c dimensionless centreline velocity - U ₀ dimensionless wall velocity - v fluid velocity in y-direction - V dimensionless v velocity, hv/ - x coordinate along channel - X dimensionless x-coordinate, x/h ² - y coordinate across channel - Y dimensionless y-coordinate, y/h - resistance coefficient, - ₀ resistance coefficient neglecting inertia - fluid density - fluid viscosity     

Supersonic nonequilibrium flow past segmental bodies of a mixture simulating the atmosphere of venus

Calculations of the flow of the mixture 0.94 CO₂+0.05 N₂+0.01 Ar past the forward portion of segmentai bodies are presented. The temperature, pressure, and concentration distributions are given as a function of the pressure ahead of the shock wave and the body velocity. Analysis of the concentration distribution makes it possible to formulate a simplified model for the chemical reaction kinetics in the shock layer that reflects the primary flow characteristics. The density distributions are used to verify the validity of the binary similarity law throughout the shock layer region calculated.The flow of a CO₂+N₂+Ar gas mixture of varying composition past a spherical nose was examined in [1]. The basic flow properties in the shock layer were studied, particularly flow dependence on the free-stream CO₂ and N₂ concentration.New revised data on the properties of the Venusian atmosphere have appeared in the literature [2, 3] One is the dominant CO₂ concentration. This finding permits more rigorous formulation of the problem of blunt body motion in the Venus atmosphere, and attention can be concentrated on revising the CO₂ thermodynamic and kinetic properties that must be used in the calculation.The problem of supersonic nonequilibrium flow past a blunt body is solved within the framework of the problem formulation of [4].Notation V body velocity - shock wave standoff - universal gas constant - ratio of frozen specific heats - hRt/m enthalpy per unit mass undisturbed stream P pressure - density - T temperature - m molecular weight - c_p specific heat at constant pressure - (X) concentration of component X (number of particles in unit mass) - R body radius of curvature at the stagnation point - _j rate of j-th chemical reaction shock layer P V ² pressure - density - TT temperature - mm molecular weight Translated from Izv. AN SSSR. Mekhanika Zhidkosti i Gaza, Vol. 5, No. 2, pp. 67–72, March–April, 1970.The author thanks V. P. Stulov for guidance in this study.     

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.