Hydrodynamic instability of the contact zone between accelerated gases

An experimental apparatus for investigating Rayleigh-Taylor instability in the transition layer between two gases at accelerations g 10⁵g₀ (g₀ is the acceleration of gravity) is described. The constantly acting acceleration is communicated to the contact zone by the compression wave formed ahead of a flame front. The linear stage of development is investigated together with the effect of the thickness of the contact zone. It is shown that on the interval 0.3 < <- ( is the wavelength of the disturbance at the edge of the contact zone) the rate of growth of the perturbation amplitude 0.5₀, where ₀ is the amplitude growth rate for media separated by an interface with a discontinuous change of density.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 15–21, November–December, 1991.     

Rheology of molten polymers

Summary TheCross equation describes the flow of pseudoplastic liquids in terms of an upper and a lower Newtonian viscosity corresponding to infinite and zero shear, and ₀, and of a third material constant related to the mechanism of rupture of linkages between particles in the intermediate, non-Newtonian flow regime, Calculation of of bulk polymers is important, since it cannot be determined experimentally. The equation was applied to the melt flow data of two low density polyethylenes at three temperatures.Using data in the non-Newtonian region covering 3 decades of shear rate to extrapolate to the zero-shear viscosity resulted in errors amounting to about onethird of the measured ₀ values. The extrapolated upper Newtonian viscosity was found to be independent of temperature within the precision of the data, indicating that it has a small activation energy.The ₀ values were from 100 to 1,400 times larger than the values at the corresponding temperatures.The values of were large compared to the values found for colloidal dispersions and polymer solutions, but decreased with increasing temperature. This shows that shear is the main factor in reducing chain entanglements, but that the contribution of Brownian motion becomes greater at higher temperatures.

---

Zusammenfassung Die Gleichung vonCross beschreibt das Fließverhalten von pseudoplastischen Flüssigkeiten durch drei Konstante: Die obereNewtonsche Viskosität (bei sehr hohen Schergeschwindigkeiten), die untereNewtonsche Viskosität ₀ (bei Scherspannung Null), und eine Materialkonstante, die vom Brechen der Bindungen zwischen Partikeln im nicht-Newtonschen Fließbereich abhängt. Die Berechnung von ist wichtig für unverdünnte Polymere, wo man sie nicht messen kann.Die Gleichung wurde auf das Fließverhalten der Schmelzen von zwei handelsüblichen Hochdruckpolyäthylenen bei drei Temperaturen angewandt. Die Werte von ₀, durch Extrapolation von gemessenen scheinbaren Viskositäten im Schergeschwindigkeitsbereich von 10 bis 4000 sec^–1 errechnet, wichen bis 30% von den gemessenen ₀-Werten ab. Die Aktivierungsenergie der war so klein, daß die-Werte bei den drei Temperaturen innerhalb der Genauigkeit der Extrapolation anscheinend gleich waren. Die ₀-Werte waren 100 bis 1400 mal größer als die-Werte.Im Verhältnis zu kolloidalen Dispersionen und verdünnten Polymerlösungen war das der Schmelzen groß, nahm aber mit steigender Temperatur ab. Deshalb wird die Verhakung der Molekülketten hauptsächlich durch Scherbeanspruchung vermindert, aber der Beitrag derBrownschen Bewegung nimmt mit steigender Temperatur zu.

     

Numerical analysis of the branching of couette flow between concentric cylinders for various wave numbers

The character of stability loss of the circular Couette flow, when the Reynolds number R passes through the critical value R₀, is investigated within a broad range of variation of the wave numbers. The Lyapunov-Schmidt method is used [1, 2]; the boundary-value problems for ordinary differential equations arising in the case of its realization are solved numerically on a computer. It is shown that the branching character substantially depends on the wave number . For all a, excluding a certain interval (₁, ₂), the usual postcritical branching takes place: at a small supercriticality the circular flow loses stability and is softly excited into a secondary stationary flow — stable Taylor vortices. For wave numbers from the interval (₁,₂) a hard excitation of Taylor vortices takes place: at a small subcriticality R=R₀–² the secondary mode is unstable and merges with the Couette flow for 0; however, for a small supercriticality in the neighborhood of a circular flow there exist no stationary modes which are different.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 132–135, May–June, 1976.     

Binary laminar boundary layer on a vertical surface in the presence of combined free and forced convection

Heat and mass transfer at a vertical surface is examined in the case of combined free and forced convection. The boundary layer equations, transformed to ordinary differential equations, contain a parameter that determines the effect of free convection on the forced motion. Criteria are offered for differentiating the free-convection, forced-convection, and combined regimes.Notation x, y coordinates - u, v velocity components - g acceleration of gravity - T temperature - kinematic viscosity - coefficient of thermal expansion - a thermal diffusivity - ₁ partial vapor density - D diffusion coefficient - W₂ mass velocity of air - independent variable - _w shear stress at wall - thermal conductivity - r latent heat of phase transition - , dimensionless temperature and partial vapor density - m^* the complex (m ₁–m _1w )/(1–m(_1w ) - c_p specific heat at constant pressure - G Grashof number - R Reynolds number - P Prandtl number - S Schmidt number     

Investigation of the flow past wings of infinite span with a blunt leading edge in the vorticity interaction regime

An asymptotic analysis of the Navier-Stokes equations is carried out for the case of hypersonic flow past wings of infinite span with a blunt leading edge when 0, Re , and M . Analytic solutions are obtained for an inviscid shock layer and inviscid boundary layer. The results of a numerical solution of the problems of vorticity interaction at the blunt edge and on the lateral surface of the wing are presented. These solutions are compared with the solution of the equations of a thin viscous shock layer and on the basis of this comparison the boundaries of the asymptotic regions are estimated.deceasedTranslated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 120–127, November–December, 1987.     

Flow in porous media I: A theoretical derivation of Darcy's law

Stokes flow through a rigid porous medium is analyzed in terms of the method of volume averaging. The traditional averaging procedure leads to an equation of motion and a continuity equation expressed in terms of the volume-averaged pressure and velocity. The equation of motion contains integrals involving spatial deviations of the pressure and velocity, the Brinkman correction, and other lower-order terms. The analysis clearly indicates why the Brinkman correction should not be used to accommodate ano slip condition at an interface between a porous medium and a bounding solid surface.The presence of spatial deviations of the pressure and velocity in the volume-averaged equations of motion gives rise to aclosure problem, and representations for the spatial deviations are derived that lead to Darcy's law. The theoretical development is not restricted to either homogeneous or spatially periodic porous media; however, the problem ofabrupt changes in the structure of a porous medium is not considered.Roman Letters A interfacial area of the - interface contained within the macroscopic system, m² - A _e area of entrances and exits for the -phase contained within the macroscopic system, m² - A interfacial area of the - interface contained within the averaging volume, m² - A ^* interfacial area of the - interface contained within a unit cell, m² - Ae area of entrances and exits for the -phase contained within a unit cell, m² - B second order tensor used to represent the velocity deviation (see Equation (3.30)) - b vector used to represent the pressure deviation (see Equation (3.31)), m^–1 - d distance between two points at which the pressure is measured, m - g gravity vector, m/s² - K Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the -phase (see Figure 2), m - characteristic length scale for the -phase (see Figure 2), m - n unit normal vector pointing from the -phase toward the -phase (n =–n ) - n _e unit normal vector for the entrances and exits of the -phase contained within a unit cell - p pressure in the -phase, N/m² - p intrinsic phase average pressure for the -phase, N/m² - p –p , spatial deviation of the pressure in the -phase, N/m² - r ₀ radius of the averaging volume and radius of a capillary tube, m - v velocity vector for the -phase, m/s - v phase average velocity vector for the -phase, m/s - v intrinsic phase average velocity vector for the -phase, m/s - v –v , spatial deviation of the velocity vector for the -phase, m/s - V averaging volume, m³ - V volume of the -phase contained within the averaging volume, m³ Greek Letters V/V, volume fraction of the -phase - mass density of the -phase, kg/m³ - viscosity of the -phase, Nt/m² - arbitrary function used in the representation of the velocity deviation (see Equations (3.11) and (B1)), m/s - arbitrary function used in the representation of the pressure deviation (see Equations (3.12) and (B2)), s^–1     

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.     

Critique of the computational Preston tube method

The general validity of the Computational Preston Tube Method (CPM) proposed by Nitsche et al. is disputed. The method involves the determination of both the skin-friction coefficient c _f and the von Karman constant from fitting two near-wall mean velocity measurements to a generalized van Driest family of velocity profiles. It is demonstrated by a detailed examination of measured mean velocity profiles that for at least one flow situation for which success with the CPM has been claimed, the laminar-turbulent transition of a flat-plate boundary layer, the method must either fail completely or lead to non-unique results. It is concluded that the applicability of the CPM is restricted to flows for which the normal law-of-the-wall is valid.     

The mirror relations and nonlinear viscoelasticity of polymer melts

An attempt is made to incorporate into a quasilinear viscoelastic constitutive equation of the Boltzmann superposition type the two mirror relations of Gleissle, as well as his relation between the steady-state first normal-stress difference and the shear viscosity curve. It is shown that the three relations can hold separately within this constitutive model, but not simultaneously, because they require a different nonlinear strain measure, namelyS ₁₂ () = – a ( – 1) (a = 0 for 1,a = 1 for 1) for the mirroring of the viscosities,S ₁₂ () = – a (–k ²/) (a = 0 for k, a = 1 for k) for the mirroring of the first normal-stress coefficients, and for the third relation. Here denotes the shear strain and erf the error function. Experimental data on melts of a low-density polyethylene, a high-density polyethylene and a polypropylene show that the mirror relations are passable approximations, but that the third relation meets reality surprisingly close if the right value ofk is used.     

Influence of vorticity on the heat transfer to blunt bodies in hypersonic flight

When blunt bodies are in hypersonic flight, a high-entropy layer of gas with nonzero vorticity is formed near their surface. The transverse gradients of the entropy, density, and gas velocity in the layer are high, which makes it necessary to take into account its absorption by the boundary layer of finite thickness . This vortex interaction is usually accompanied by an increase in the heat flux q and the frictional stress on the wall compared with their values as calculated in accordance with the classical scheme of a thin boundary layer, when the parameters on the outer edge of the boundary layer are set equal to the inviscid parameters on the body. This effect has been investigated on the side surface of slender (with angle 1 to the undisturbed flow) blunt bodies in a hypersonic stream [1–3]. It is shown in the present paper that the effect can have a stronger and even qualitative influence on the flow over blunt bodies with 1 if the radius of curvature R_s of the detached shock wave on the axis is small compared with the midsection radius R of the body. It is shown that the distributions of the heat fluxes with allowance for the vorticity of the inviscid shock layer are similar in the case of slightly blunt (r₀/R 0) cones with half-angles less than a critical _*.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 50–57, March–April, 1981.     

Instantaneous visualization of surface flows

An experimental system is described for visualizing the surface flow of a wing, using an oil smoke tracer technique. The method leads to the determination of the instantaneous velocity direction at the output of surface injectors. A preliminary investigation is made on a flat plate to optimize the conditions of oil smoke injection. Then, the visualization is performed on the upperside of a sweptback wing in the vicinity of the reattachment of the vortex flow. This visualization technique can be applied to other types of wall flows — separated or not — around various bodies.List of symbols b wing span - c _n normal (to leading edge) chord - c _r streamwise (or root) chord - d diameter of the injectors - distance from the apex along the leading edge - relative distance from the apex along the leading edge ( = /C _d) - sweep angle - e injector geometric parameter (e = d/l) - angle of attack - K injection parameter - l length of the injectors - v kinematic viscosity - P _t, P_s total and static pressure of the flow - P _inj injection pressure - P _r reduced pressure (P _r = (P_inj – P_t)/(P_t – P_s)) - Re flow Reynolds number (Re = V ·c _n/v) - Re _i injector Reynolds number (Re = V ·d/v) - s curvilinear distance along c _d - s relative curvilinear distance along c _d(s = s/c _d) - V infinite upstream flow velocity     

Calculation of heat transfer accompanying hypersonic three-dimensional flow of air over slender blunt-nosed cones

The effective length method [1, 2] has been used to make systematic calculations of the heat transfer for laminar and turbulent boundary layers on slender blunt-nosed cones at small angles of attack ( + 5° in a separationless hypersonic air stream dissociating in equilibrium (half-angles of the cones 0 20°, angles of attack 0 15°, Mach numbers 5 M 25). The parameters of the gas at the outer edge of the boundary layer were taken equal to the inviscid parameters on the surface of the cones. Analysis of the results leads to simple approximate dependences for the heat transfer coefficients.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 173–177, September–October, 1981.     

Waves in a plasma with Hall dispersion

In a magnetohydrodynamic approximation, an investigation is made of the propagation of waves in a plasma, whose characteristic frequency is much less than the collision frequency of the electrons _e ^–1. It is assumed that the magnetic field is sufficiently strong so that the equality _ee1 will be satisfied, where _e is the cyclotron frequency of the rotation of the electrons. With large magnetic Reynolds numbers (R_m1), which are characteristic for many astrophysical problems, this latter condition leads to a need to take account of dispersion effects connected with Hall currents, in the absence of Joule dissipation. The dispersion equation for the propagation of small perturbations is analyzed in the limiting cases of weak dispersion and of a wave propagating along the magnetic field. In the case of weak dispersion, an equation is derived for nonlinear waves. The solutions are found in the form of stationary solitons. The region of such solutions is analyzed. A typical example of a medium with Hall dispersion is an interplanetary plasma, in which the parameter _ee is generally great.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 108–113, May–June, 1974.     

On natural convection from a vertical plate with a prescribed surface heat flux in porous media

This paper presents a theoretical and numerical investigation of the natural convection boundary-layer along a vertical surface, which is embedded in a porous medium, when the surface heat flux varies as (1 +x ²)), where is a constant andx is the distance along the surface. It is shown that for > -1/2 the solution develops from a similarity solution which is valid for small values ofx to one which is valid for large values ofx. However, when -1/2 no similarity solutions exist for large values ofx and it is found that there are two cases to consider, namely < -1/2 and = -1/2. The wall temperature and the velocity at large distances along the plate are determined for a range of values of .Notation g Gravitational acceleration - k Thermal conductivity of the saturated porous medium - K Permeability of the porous medium - l Typical streamwise length - q _w Uniform heat flux on the wall - Ra Rayleigh number, =gK(q _w /k)l/(v) - T Temperature - Too Temperature far from the plate - u, v Components of seepage velocity in the x and y directions - x, y Cartesian coordinates - Thermal diffusivity of the fluid saturated porous medium - The coefficient of thermal expansion - An undetermined constant - Porosity of the porous medium - Similarity variable, =y(1+x ) ^/3/x ^1/3 - A preassigned constant - Kinematic viscosity - Nondimensional temperature, =(T – T )Ra^1/3 k/qw - Similarity variable, = =y(log_e x)^1/3/x ^2/3 - Similarity variable, =y/x ^2/3 - Stream function     

Parameter identification in a soil with constant diffusivity

We consider infiltration into a soil that is assumed to have hydraulic conductivity of the form K = K = K_se^h and water content of the form = K – _r. Here h denotes capillary pressure head while K_s, , and _r represent soil specific parameters. These assumptions linearize the flow equation and permit a closed form solution that displays the roles of all the parameters appearing in the hydraulic function K and . We assume K_s and _r to be known. A measurement of diffusivity fixes the product of and resulting in a parameter identification problem for one parameter. We show that this parameter identification problem, in some cases, has a unique solution. We also show that, in some cases, this parameter identification problem can have multiple solutions, or no solution. In addition it is shown that solutions to the parameter identification problem can be very sensitive to small changes in the problem data.     

On the prediction of riblet performance with engineering turbulence models

The paper reports the outcome of a numerical study of fully developed flow through a plane channel composed of ribleted surfaces adopting a two-equation turbulence model to describe turbulent mixing. Three families of riblets have been examined: idealized blade-type, V-groove and a novel U-form that, according to computations, achieves a superior performance to that of the commercial V-groove configuration. The maximum drag reduction attained for any particular geometry is broadly in accord with experiment though this optimum occurs for considerably larger riblet heights than measurements indicate. Further explorations bring out a substantial sensitivity in the level of drag reduction to the channel Reynolds number below values of 15 000 as well as to the thickness of the blade riblet. The latter is in accord with the trends of very recent, independent experimental studies.Possible shortcomings in the model of turbulence are discussed particularly with reference to the absence of any turbulence-driven secondary motions when an isotropic turbulent viscosity is adopted. For illustration, results are obtained for the case where a stress transport turbulence model is adopted above the riblet crests, an elaboration that leads to the formation of a plausible secondary motion sweeping high momentum fluid towards the wall close to the riblet and thereby raising momentum transport.Nomenclature c _f Skin friction coefficient - c _f Skin friction coefficient in smooth channel at the same Reynolds number - k Turbulent kinetic energy - K ⁺ k/ _w - h Riblet height - S Riblet width - H Half height of channel - Re Reynolds number = volume flow/unit width/ - Modified turbulent Reynolds number - R _t turbulent Reynolds numberk ²/ - P _k Shear production rate ofk, _t (U _i /x _j + U _j /x _i ) U _i /x _j - dP/dz Streamwise static pressure gradient - U _i Mean velocity vector (tensor notation) - U Friction velocity, _w/ where _w=–H dP/dz - W Mean velocity - W _b Bulk mean velocity through channel - y ⁺ yU /v. Unless otherwise stated, origin is at wall on trough plane of symmetry - Kinematic viscosity - _t Turbulent kinematic viscosity - Turbulence energy dissipation rate - Modified dissipation rate – 2(k ^1/2/x _j )² - Density - _k , Effective turbulent Prandtl numbers for diffusion ofk and      

Base-pressure pulsations behind a cylinder and disk in a subsonic stream flow

Base-pressure fluctuations behind a long cylinder (l/d 5–10) and the disk (l/ d 0.0) is investigated experimentally in this paper. The spectral and correlation characteristics of the base-pressure fluctuations behind axisymmetric bodies at a Mach number M 1.0 are generalized on the basis of the data obtained and the results of other authors.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 181–183, January–February, 1977.     

The chaotic response of a fluttering panel: The influence of maneuvering

The influence of maneuvering on the chaotic response of a fluttering buckled plate on an aircraft has been studied. The governing equations, derived using Lagrangian mechanics, include geometric non-linearities associated with the occurrence of tensile stresses, as well as coupling between the angular velocity of the maneuver and the elastic degrees of freedom. Numerical simulation for periodic and chaotic responses are conducted in order to analyze the influence of the pull-up maneuver on the dynamic behavior of the panel. Long-time histories phase-plane plots, and power spectra of the responses are presented. As the maneuver (load factor) increases, the system exhibits complicated dynamic behavior including a direct and inverse cascade of subharmonic bifurcations, intermittency, and chaos. Beside these classical routes of transition from a periodic state to chaos, our calculations suggest amplitude modulation as a possible new mode of transition to chaos. Consequently this research contributes to the understanding of the mechanisms through which the transition between periodic and strange attractors occurs in, dissipative mechanical systems. In the case of a prescribed time dependent maneuver, a remarkable transition between the different types of limit cycles is presented.Nomenclature a plate length - a _r u _r /h - D plate bending stiffness - E modulus of elasticity - g acceleration due to gravity - h plate thickness - j₁,j₂,j₃ base vectors of the body frame of reference - K spring constant - M Mach number - n 1 + ₀/g - N ₁ applied in-plane force - p–p aerodynamic pressure - P pa ⁴/Dh - q ₀/2 - Q _r generalized Lagrangian forces - R rotation matrix - R ₄ N, a ²/D - t time - kinetic energy - u plate deflection - u displacement of the structure - u _r modal amplitude - v₀ velocity - x coordinates in the inertial frame of reference - z coordinates in the body frame of reference - Ka/(Ka+Eh) - - elastic energy - 2qa ³/D - a/_mh - Poisson's ratio - material coordinates - air density - _m plate density - - _r prescribed functions - _r sin(r z/a) - angular velocity - a/v₀ - skew-symmetric matrix form of the angular velocity     

Unsteady supersonic flow around blunt bodies with detached shock wave

The numerical method of calculating the supersonic three-dimensional flow about blunt bodies with detached shock wave presented in [1–3] is applied to the case of unsteady flow. The formulation of the unsteady problem is analogous to that of [4], which assumes smallness of the unsteady disturbances.The paper presents some results of a study of the unsteady flow about blunt bodies over a wide range of variation of the Mach number M=1.50– and dimensionless oscillation frequency l/V=0–1.0. A comparison is made with the results obtained from the Newton theory.     

Hypersonic flow over a narrow delta plate at large angles of attack

The problem of hypersonic flow over a flat delta plate with a high sweepback anglex at angles of attack close to /2 is solved using a numerical algorithm based on transition to the conical solution. The existence of conical flow at /2 with the velocity vector directed towards the apex of the plate is established. Values ofC _p/sin² and the thickness of the shock layer in the plane of symmetry of the plate are given as functions of the hypersonic similarity parameterk=tan tanx. A comparison of the calculated and experimental data shows that they are in good agreement.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.5, pp. 183–185, September–October, 1992.