Simulation method in the kinetic theory of gases

We develop a method that is based on use of the simulation equation vf/v=vl ^–1(f ₀–f), where v is the modulus of the molecular velocity andl is the mean free-path length. A number of general properties of the model are clarified and the transition to the limit of a continuous medium is discussed in detail.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 110–115, November–December, 1973.     

On the Global Stability of One Nonlinear Difference Equation

We give exact sufficient conditions for the global stability of the zero solution of the difference equation x _{n + 1} = qx _n + f _n (x _n , ..., x _{n – k} ), n , where the nonlinear functions f _n satisfy the conditions of negative feedback and sublinear growth.__________Translated from Neliniini Kolyvannya, Vol. 7, No. 4, pp. 487–494, October–December, 2004.     

Flow structure in motion of a spherical drop in a fluid medium at intermediate Reynolds numbers

The problem of flow of a viscous fluid around a spherical drop has been examined for the limiting case of small and large Reynolds numbers in several investigations (see [1–3], for instance; there is a detailed review of various approximate solutions in [4]). For the intermediate range of Reynolds numbers (approximately 1Re100), where numerical integration of the complete Navier-Stokes equations is necessary, there are solutions of special cases of the problem —flow of air around a solid sphere [5–7], a gas bubble [8, 9], and water drops [10]. The present paper deals with flow around a spherical drop at intermediate Reynolds numbers up to Re=200 for arbitrary values of the ratio of dynamic viscosities =₁/₂ inside and outside the drop. It is shown that a return flow can arise behind the drop in flow without separation. In such conditions the circulatory flow inside the drop breaks up. An approximate formula for the drag coefficient of the drop is given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 8–15, January–February, 1976.We thank L. A. Galin, G. I. Petrov, L. A. Chudov, and participants in the seminars led by them for useful discussions.     

Theory of a laval nozzle for a two-phase mixture containing particles of small lag

A two-velocity and two-temperature model is considered for a continuous medium in relation to the flow of a mixture of gas and particles in the subsonic, transsonic, and supersonic parts of a Laval nozzle. It is assumed that the particles are small, and hence that the coefficients f and ^q, which define the interaction with the gas, are large (these coefficients are inversely proportional to the square of the particle radius for a Stokes mode of flow). This means that the velocity or thermal lag of the particles relative to the gas is small. The solution is sought as expansions with respect to the small parameters ₁=1/f and ₂=1/^q.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 89–100, July–August, 1973.     

Flow of second-order fluids over an enclosed rotating disc

Summary This paper is devoted to a study of the flow of a second-order fluid (flowing with a small mass rate of symmetrical radial outflow m, taken negative for a net radial inflow) over a finite rotating disc enclosed within a coaxial cylinderical casing. The effects of the second-order terms are observed to depend upon two dimensionless parameters ₁ and ₂. Maximum values ₁ and ₂ of the dimensionless radial distances at which there is no recirculation, for the cases of net radial outflow (m>0) and net radial inflow (m<0) respectively, decrease with an increase in the second-order effects [represented by T(=₁+₂)]. The velocities at ₁ and ₂ as well as at some other fixed radii have been calculated for different T and the associated phenomena of no-recirculation/recirculation discussed. The change in flow phenomena due to a reversal of the direction of net radial flow has also been studied. The moment on the rotating disc increases with T.Nomenclature , , z coordinates in a cylindrical polar system - z ₀ distance between rotor and stator (gap length) - =/z ₀, dimensionless radial distance - =z/z ₀, dimensionless axial distance - _s = _s/z₀, dimensionless disc radius - V =(u, v, w), velocity vector - dimensionless velocity components - uniform angular velocity of the rotor - , p fluid density and pressure - P =p/(₂ z ₀₂ ² , dimensionless pressure - ₁, ₂, ₃ kinematic coefficients of Newtonian viscosity, elastico-viscosity and cross-viscosity respectively - ₁, _{2 2}/z ₀ ² , resp. ₃/z ₀ ² , dimensionless parameters representing the ratio of second-order and inertial effects - m = , mass rate of symmetrical radial outflow - l a number associated with induced circulatory flow - Rm =m/(z ₀₁), Reynolds number of radial outflow - R _l =l/(z ₀₁), Reynolds number of induced circulatory flow - Rz =z ₀ ² /₁, Reynolds number based on the gap - ₁, ₂ maximum radii at which there is no recirculation for the cases Rm>0 and Rm<0 respectively - ₁(T), ₂(T) ₁ and ₂ for different T - U _1(T) ⁽⁺⁾ = dimensionless radial velocity, Rm>0 - V _1(T) ⁽⁺⁾ = , dimensionless transverse velocity, Rm>0 - U _2(T) ^(–) = , dimensionless radial velocity, Rm=–Rn<0, m=–n - V _2(T) ^(–) = , dimensionless transverse velocity, Rm<0 - C _m moment coefficient     

Asymptotic solution to the problem of hypersonic flow over bodies in the neighborhood of the separation point of a thin shock layer

A study is made of the problem of hypersonic flow of an inviscid perfect gas over a convex body with continuously varying curvature. The solution is sought in the framework of the asymptotic theory of a strongly compressed gas [1–4] in the limit M when the specific heat ratio tends to 1. Under these assumptions, the disturbed flow is situated in a thin shock layer between the body and the shock wave. At the point where the pressure found by the Newton-Buseman formula vanishes there is separation of the flow and formation of a free layer next to the shock wave [1–4]. The singularity of the asymptotic expansions with respect to the parameter ₁ = ( –1)/( + 1) associated with separation of the strongly compressed layer has been investigated previously by various methods [3–9]. Local solutions to the problem valid in the neighborhood of the singularity have been obtained for some simple bodies [3–7]. Other solutions [7, 9] eliminate the singularity but do not give the transition solution entirely. In the present paper, an asymptotic solution describing the transition from the attached to the free layer is constructed for a fairly large class of flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 99–105, January–February, 1982.     

Supersonic air flow past a sphere with account for vibrational relaxation

A numerical solution is obtained for the problem of air flow past a sphere under conditions when nonequilibrium excitation of the vibrational degrees of freedom of the molecular components takes place in the shock layer. The problem is solved using the method of [1]. In calculating the relaxation rates account was taken of two processes: 1) transition of the molecular translational energy into vibrational energy during collision; 2) exchange of vibrational energy between the air components. Expressions for the relaxation rates were computed in [2]. The solution indicates that in the state far from equilibrium a relaxation layer is formed near the sphere surface. A comparison is made of the calculated values of the shock standoff with the experimental data of [3].Notation uV_max, vV_max velocity components normal and tangential to the sphere surface - V_max maximal velocity - P V _max ² pressure - density - TT temperature - e_viRT vibrational energy of the i-th component per mole (i=–O₂, N₂) - =rb–1 shock wave shape - a _f the frozen speed of sound - HRT/m gas total enthalpy     

Flow structure of wake behind a rotationally oscillating circular cylinder

Flow around a circular cylinder oscillating rotationally with a relatively high forcing frequency has been investigated experimentally. The dominant parameters affecting this experiment are the Reynolds number (Re), oscillation amplitude (θ_A), and frequency ratio F_R=f_f/f_n, where f_f is the forcing frequency and f_n is the natural frequency of vortex shedding. Experiments were carried out under conditions of Re=4.14×10³, 0°θ_A60° and 0.0F_R2.0. Rotational oscillation of the cylinder significantly modified the flow structure in the near-wake. Depending on the frequency ratio F_R, the cylinder wake showed five different flow regimes, each with a distinct wake structure. The vortex formation length and the vortex shedding frequency were greatly changed before and after the lock-on regime where vortices shed at the same frequency as the forcing frequency. The lock-on phenomenon always occurred at F_R=1.0 and the frequency range of the lock-on regime expanded with increasing oscillation amplitude θ_A. In addition, the drag coefficient was reduced when the frequency ratio F_R was less than 1.0 (F_R<1.0) while fixing the oscillation amplitude at θ_A=30°. When the oscillation amplitude θ_A was used as a control parameter at a fixed frequency ratio F_R=1.0 (lock-on regime), the drag reduction effect was observed at all oscillation amplitudes except for the case of θ_A=30°. This type of active flow control method can be used effectively in aerodynamic applications while optimizing the forcing parameters.     

On the Construction of a Solution of a Degenerate Linear System of Differential Equations in the Neighborhood of an Irregular Singular Point

We investigate the asymptotics of the general solution of the linear system of differential equations with irregular singular point

Optimum lifting wings with specified longitudinal trim characteristics

The case of an infinitely slender wing that slightly disturbs a supersonic ideal gas flow is considered. The plan form and the free-stream Mach number M are given. The optimum surface of the wing y=g(x, z) is determined as a result of finding a bounded function of the local angles of attack _M=g(x, z)/x that minimizes the drag coefficient c_x for given values of the lift coefficient c_y and the pitching moment coefficient m_z. The problem is solved in the class of piecewise-constant functions for wings of complex geometry [1].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 185–189, July–August, 1987.     

Interactions in general continua with microstructure

This paper studies the behavior of the one dimensional Broadwell model of a discrete three velocity gas on a bounded domain 0 x 1 with specularly reflective boundary condition at x = 0, 1. For smooth initial data we show that the initial boundary value problem possesses a unique smooth solution which tends as t to a free state consisting of traveling waves f _1e (x – ct), f _2e (x + ct), f _3e (x) where each f _ie is 2-periodic. The convergence is in the weak* topology of an appropriate Orlicz-Banach state space. No smallness assumptions are made on the data.In memory of Ronald J. DiPerna     

Three-dimensional hypersonic rarefied gas flow round bodies of revolution

A proposed method of studying three-dimensional rarefied gas flow around a body of revolution is based on the numerical solution of model kinetic equations. By way of example, the problem is considered of hypersonic flow round an ellipsoid of revolution whose velocity vector forms an angle of ₀ with the axis of symmetry of the body and is located in the plane of symmetry. A study is made of the effect of the angle of attack, surface temperature and Knudsen number on the aerodynamic characteristics of the body.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti 1 Gaza, No. 1, pp. 184–186, January–February, 1986.     

Two-phase geothermal flows with conduction and the connection with Buckley-Leverett theory

Two-phase flows of boiling water and steam in geothermal reservoirs satisfy a pair of conservation equations for mass and energy which can be combined to yield a hyperbolic wave equation for liquid saturation changes. Recent work has established that in the absence of conduction, the geothermal saturation equation is, under certain conditions, asymptotically identical with the Buckley-Leverett equation of oil recovery theory. Here we summarise this work and show that it may be extended to include conduction. In addition we show that the geothermal saturation wave speed is under all conditions formally identical with the Buckley-Leverett wave speed when the latter is written as the saturation derivative of a volumetric flow.Roman Letters C(P, S,q) geothermal saturation wave speed [ms^–1] (14) - c _t (P, S) two-phase compressibility [Pa^–1] (10) - D(P, S) diffusivity [m s^–2] (8) - E(P, S) energy density accumulation [J m^–3] (3) - g gravitational acceleration (positive downwards) [ms^–2] - h _w (P),h _w (P) specific enthalpies [J kg^–1] - J _M (P, S,P) mass flow [kg m^–2 s^–1] (5) - J _E (P, S,P) energy flow [J m^–2s^–1] (5) - k absolute permeability (constant) [m²] - k _w (S),k _s (S) relative permeabilities of liquid and vapour phases - K formation thermal conductivity (constant) [Wm^–1 K^–1] - L lower sheetC<0 in flow plane - m, c gradient and intercept - M(P, S) mass density accumulation [kg m^–3] (3) - O flow plane origin - P(x,t) pressure (primary dependent variable) [Pa] - q volume flow [ms^–1] (6) - S(x, t) liquid saturation (primary dependent variable) - S ^*(x,t) normalised saturation (Appendix) - t time (primary independent variable) [s] - T temperature (degrees Kelvin) [K] - T _sat(P) saturation line temperature [K] - TdT _sat/dP saturation line temperature derivative [K Pa^–1] (4) - T _c ,T _D convective and diffusive time constants [s] - u _w (P),u _s (P),u _r (P) specific internal energies [J kg^–1] - U upper sheetC > 0 in flow plane - U(x,t) shock velocity [m s^–1] - x spatial position (primary independent variable) [m] - X representative length - x, y flow plane coordinates - z depth variable (+z vertically downwards) [m] Greek Letters _P , _S remainder terms [Pa s^–1], [s^–1] - double-valued saturation region in the flow plane - h =h _s –h _w latent heat [J kg^–1] - = _w – _s density difference [kg m^–3] - line envelope - =D _K /D ₀ diffusivity ratio - porosity (constant) - _w (P), _s (P), _t (P, S) dynamic viscosities [Pa s] - v _w (P),v _s (P) kinematic viscosities [m²s^–1] - v ₀ =kh/KT kinematic viscosity constant [m² s^–1] - ₀ =v ₀ dynamic viscosity constant [m² s^–1] - _w (P), _s (P) density [kg m^–3] Suffixes r rock matrix - s steam (vapour) - w water (liquid) - t total - av average - 0 without conduction - K with conduction     

Flow of a liquid metal coolant past a heat generating cylinder

The coupled problem whereby a solid heat generating cylinder is being cooled in steady state by a coolant in potential flow is investigated. An analytical technique for determining the temperature distributions in the solid and the fluid is presented. Numerical studies for six Péclet numbers (0.9<Pe<11.3) and three thermal conductivity ratios (0.31<K<3.1) were carried out.The surface hot-spot temperature and center temperature are presented graphically as functions of the Péclet number with the thermal conductivity ratio as a parameter. The average Nusselt number is found to be proportional to the Péclet number to approximately the one-half power. For the special case of constant surface temperature (uncoupled problem), the variation of local Nusselt number with angle measured from the forward stagnation point is in excellent agreement with the result presented by Grosh and Cess [6].Nomenclature ce _m(, –q) Mathieu function, periodic - D _n Fourier coefficient for solid temperature distribution - E _n Fourier coefficient for fluid temperature distribution - E() a term defined by equation (12), degree - F() a term defined by equation (13), degree - Fek _m(z, –q) modified Mathieu function, non-periodic - Fek _m(z, –q) derivative of Fek _m(z, –q) - h local heat transfer coefficient, energy/time area degree - average heat transfer coefficient, energy/time area degree - h _m mean heat-transfer coefficient, energy/time area degree - k _f thermal conductivity of fluid, energy/time length degree - K thermal conductivity ratio, k _f/k _s - k _s thermal conductivity of solid, energy/time length degree - Nu local Nusselt number, 2Rh/k _f - average Nusselt number defined by equation (55) - (Nu)_m mean Nusselt number defined by equation (57) - Pe Péclet number, 2RU/ _f - Q rate of heat generation per unit volume, energy/time volume - q parameter of Mathieu function, (Pe/4)² - q normal heat flux, energy/time area - R cylinder radius, length - Re Reynolds number, 2R/ - r radial position variable, length - T temperature, degree - T ₀ constant surface temperature, degree - T temperature of fluid at infinity, degree - T _e temperature at center of cylinder, degree - T _f temperature of fluid, degree - T _s temperature of solid, degree - T _w surface temperature, degree - surface hot-spot temperature, degree - reduced temperature, (T–)/E(1) - U approach velocity of flowing fluid, length/time - v velocity component in the direction, length/time - v _r velocity component in the r direction, length/time - z logarithm of Greek symbols _f thermal diffusivity of the fluid, (length)²/time - reduced radius, r/R - angular position variable measured from the trailing stagnation point, radians - kinematic viscosity, (length)²/time - angular position variable measured from the forward stagnation point, degree     

Some problems of nonisothermal steady flow of a viscous fluid

The paper presents solutions to the problems of plane Couette flow, axial flow in an annulus between two infinite cylinders, and flow between two rotating cylinders. Taking into account energy dissipation and the temperature dependence of viscosity, as given by Reynolds's relation =₀ exp (–T) (₀, =const). Two types of boundary conditions are considered: a) the two surfaces are held at constant (but in general not equal) temperatures; b) one surface is held at a constant temperature, the other surface is insulated.Nonisothermal steady flow in simple conduits with dissipation of energy and temperature-dependent viscosity has been studied by several authors [1–11]. In most of these papers [1–6] viscosity was assumed to be a hyperbolic function of temperature, viz. =_m 1/1+²(T–T_m.Under this assumption the energy equation is linear in temperature and can he easily integrated. Couette flow with an exponential viscosity-temperature relation. =₀ e ^–T (₀, =const), (0.1) was studied in [7, 8]. Couette flow with a general (T) relation was studied in (9).Forced flow in a plane conduit and in a circular tube with a general (T) relation was studied in [10]. In particular, it has been shown in [10] that in the case of sufficiently strong dependence of viscosity on temperature there can exist a critical value of the pressure gradient, such that a steady flow is possible only for pressure gradients below this critical value.In a previous work [11] the authors studied Polseuille flow in a circular tube with an exponential (T) relation. This thermohydrodynamic problem was reduced to the problem of a thermal explosion in a cylindrical domain, which led to the existence of a critical regime. The critical conditions for the hydrodynamic thermal explosion and the temperature and velocity profiles were calculated.In this paper we treat the problems of Couette flow, pressureless axial flow in an annulus, and flow between two rotating cylinders taking into account dissipation and the variation of viscosity with temperature according to Reynolds's law (0.1). The treatment of the Couette flow problem differs from that given in [8] in that the constants of integration are found by elementary methods, whereas in [8] this step involved considerable difficulties. The solution to the two other problems is then based on the Couette problem.     

The use of sweep-frequency excitation for unsteady pressure measurement

The use of sweep-frequency excitation for rapid measurement of time-dependent pressures on wind-tunnel models is examined. Results obtained from two different wind-tunnels covering the Mach number range from 0.2 to 0.85, and a wide range of flow conditions, are compared with measurements made using the slower, traditional method of discrete-frequency excitation. It is concluded that the sweep-frequency excitation method can reduce testing time in certain flow conditions with no significant loss in accuracy.List of symbols M Mach number - p broadband rms local static pressure - q 12u ² (dynamic pressure) - R(Cp/) real (in-phase) part of oscillatory Cp/ - I(Cp/) imaginary (in-quadrature) part of oscillatory Cp/ - x/c chord station - wing incidence - canard or wing oscillatory amplitude (plotted in radians unless otherwise stated) - spanwise station - _c canard static incidence - _c canard effective incidence ( _c = 1.89 + _c –0.6) - (T) function of time - ² coherence function The coherence function between two signals x(f), y(f) is defined as - where - G _xy (f) = cross spectral density function between x and y - G _xx (f) = auto spectral density function of x - G _yy (f) = auto spectral density function of y - f = frequency     

Convective heat transfer in the thermal entrance region of parallel-flow noncircular duct heat exchanger arrays

Convective heat transfer properties of a hydrodynamically fully developed flow, thermally developing flow in a parallel-flow, and noncircular duct heat exchanger passage subject to an insulated boundary condition are analyzed. In fact, due to the complexity of the geometry, this paper investigates in detail heat transfer in a parallel-flow heat exchanger of equilateral-triangular and semicircular ducts. The developing temperature field in each passage in these geometries is obtained seminumerically from solving the energy equation employing the method of lines (MOL). According to this method, the energy equation is reformulated by a system of a first-order differential equation controlling the temperature along each line.Temperature distribution in the thermal entrance region is obtained utilizing sixteen lines or less, in the cross-stream direction of the duct. The grid pattern chosen provides drastic savings in computing time. The representative curves illustrating the isotherms, the variation of the bulk temperature for each passage, and the total Nusselt number with pertinent parameters in the entire thermal entry region are plotted. It is found that the log mean temperature difference (T _LM), the heat exchanger effectiveness, and the number of transfer units (NTU) are 0.247, 0.490, and 1.985 for semicircular ducts, and 0.346, 0.466, and 1.345 for equilateral-triangular ducts.

---

Konvektiver Wärmeübergang im thermischen Einlaufgebiet von Gleichstromwärmetauschern mit nichtkreisförmigen Strömungskanälen

Zusammenfassung Die Untersuchung bezieht sich auf das konvektive Wärmeübertragungsverhalten eines Gleichstromwärmetauschers mit nichtkreisförmigen Strömungskanälen bei hydraulisch ausgebildetet, thermisch einlaufender Strömung unter Aufprägung einer adiabaten Randbedingung. Zwei Fälle komplizierter Geometrie, nämlich Kanäle mit gleichseitig dreieckigen und halbkreisförmigen Querschnitten, werden bezüglich des Wärmeübergangsverhaltens bei Gleichstromführung eingehend analysiert. Das sich entwickelnde Temperaturfeld in jedem Kanal von der eben spezifizierten Querschnittsform wird halbnumerisch durch Lösung der Energiegleichung unter Einsatz der Linienmethode (MOL) erhalten. Dieser Methode entsprechend erfolgt eine Umformung der Energiegleichung in ein System von Differentialgleichungen erster Ordnung, welches die Temperaturverteilung auf jeder Linie bestimmt.Die Temperaturverteilung im Einlaufgebiet wird unter Vorgabe von 16 oder weniger Linien über dem Kanalquerschnitt erhalten, wobei die gewählte Gitteranordnung drastische Einsparung an Rechenzeit ergibt. Repräsentative Kurven für das Isothermalfeld, den Verlauf der Mischtemperatur für jeden Kanal und die Gesamt-Nusseltzahl als Funktion relevanter Parameter im gesamten Einlaufgebiet sind in Diagrammform dargestellt. Es zeigt sich, daß die mittlere logarithmische Temperaturdifferenz (T _LM), der Wärmetauscherwirkungsgrad und die Anzahl der Übertragungseinheiten (NTU) folgende Werte annehmen: 0,247, 0,490 und 1,985 für halbkreisförmige Kanäle sowie 0,346, 0,466 und 1,345 für gleichseitig dreieckige Kanäle.

---

Nomenclature A cross sectional area [m²] - a characteristic length [m] - C _c specific heat of cold fluid [J kg^–1 K^–1] - C _h specific heat of hot fluid [J kg^–1 K^–1] - C _p specific heat [J kg^–1 K^–1] - C _r specific heat ratio,C _r=C _c/C_h - D _h hydraulic diameter of duct [m] - f friction factor - k thermal conductivity of fluid [Wm^–1 K^–1] - L length of duct [m] - m mass flow rate of fluid [kg s^–1] - N factor defined by Eq. (20) - NTU number of transfer units - Nu _{x, T} local Nusselt number, Eq. (19) - P perimeter [m] - p pressure [KN m^–2] - Pe Peclet number,RePr - Pr Prandtl number,/ - Q _T total heat transfer [W], Eq. (13) - Q ideal heat transfer [W], Eq. (14) - Re Reynolds number,D _h/ - T temperature [K] - T _b bulk temperature [K] - T _e entrance temperature [K] - T _w circumferential duct wall temperature [K] - u, U dimensional and dimensionless velocity of fluid,U=u/u - , dimensional and dimensionless mean velocity of fluid - w generalized dependent variable - X dimensionless axial coordinates,X=D _h ² /a ² x* - x, x* dimensional and dimensionless axial coordinate,x*=x/D _hPe - y, Y dimensional and dimensionless transversal coordinates,Y=y/a - z, Z dimensional and dimensionless transversal coordinates,Z=z/a Greek symbols thermal diffusivity of fluid [m² s^–1] - * right triangular angle, Fig. 2 - independent variable - T _LM log mean temperature difference of heat exchanger - effectiveness of heat exchanger - generalized independent variable - dimensionless temperature - _b dimensionless bulk temperature - dynamic viscosity of fluid [kg m^–1 s^–1] - kinematic viscosity of fluid [m² s^–1] - density of fluid [kg m^–3] - heat transfer efficiency, Eq. (14) - generalized dependent variable     

Spectral and correlation characteristics of the turbulent boundary layer on an extended flexible cylinder

In marine geophysical seismological prospecting extensive use is made of towed receiving systems consisting of extended flexible cylinders containing acoustic sensors over which the water flows in the longitudinal direction. The boundary layer pressure fluctuations on the cylinder surface are picked up by the sensors as hydrodynamic noise. This paper is concerned with the study of the energy and spacetime characteristics of the pressure fluctuations in the turbulent boundary layer on an extended flexible cylinder in a longitudinal flow. The pressure fluctuations on the cylinder surface have been investigated experimentally for Re_X=(2–4)·10⁷, a dimensionless diameter of the pressure fluctuation sensors d _p ⁺ =d_pu^*/=500, where d_p is the sensor diameter, u^* the dynamic viscosity, and the viscosity coefficient, and frequencies 0.02311.259 (=^*/U). The spectral and correlation characteristics of the pressure fluctuations on the surface of the flexible cylinder are found to differ from the corresponding characteristics for a rigid cylinder [1–4].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i aza, No, 5, pp. 49–54, September–October, 1989.     

Streamwise pseudo-vortical structures and associated vorticity in the near-wall region of a wall-bounded turbulent shear flow

Streamwise pseudo-vortical motions near the wall in a fully-developed two-dimensional turbulent channel flow are clearly visualized in the plane perpendicular to the flow direction by a sophisticated hydrogen-bubble technique. This technique utilizes partially insulated fine wires, which generate hydrogen-bubble clusters at several distances from the wall. These flow visualizations also supply quantitative data on two instantaneous velocity components, and w, as well as the streamwise vorticity, _x . The vorticity field thus obtained shows quasi-periodicity in the spanwise direction and also a double-layer structure near the wall, both of which are qualitatively in good agreement with a pseudo-vortical motion model of the viscous wall-region.List of symbols C _i ,c _i ,d _i constants in Eqs. (2), (3) and (4) - H channel width (m) - Re _H Reynolds number (= U _c H/) - Re Reynolds number (= U _c /) - T period (s) - t time (s) - U mean streamwise velocity (m/s) - U _c center-line velocity (m/s) - u friction velocity (m/s) - u, , w velocity fluctuations (m/s) - x, y, z coordinates (m) - ^* displacement thickness (m) - momentum thickness (m) - mean low-speed streak spacing (m) - kinematic viscosity (m²/s) - phase difference - _x streamwise vorticity fluctuation (1/s) - ( )⁺ normalized by u and - (^–) root mean square value - (^–) statistical average This paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984     

Analysis of ellipsoid compressed by two axial concentrated forces at two ends

Integral equation method and photoelastic experiment are used for the stress analysis of an axial compressive ellipsoid. Let the concentrated forces and the centers of compression, with symmetrical unknown intensive functions X₁(c)=X₁(–c) and X₂(c)=X₂(–c) respectively, be distributed symmetrically to =0 plane along the axis z(=–c) in [a,) and [–a,–) of the elastic space, in addition to a pair of equal and opposite axial forces acting on z=a and z=–a. We can reduce the problem of an axial compressive ellipsoid to two coupled Fredholm integral equations of the first kind. Furthermore, numerical calculation is then made. Two photo-elastic models of ellipsoid were analysed by Freezing and Cutting method, and the results, in which ₂ is quite nearly to those obtained by integral equation method, had been used in the analysis of the data of compressive rock specimens.