Solitary waves in a layer of viscous liquid

Solitary waves in a thin layer of viscous liquid which is running down a vertical surface under the action of gravity are investigated. The existence of such waves was demonstrated in the experiments of [1, 2]. The difficulties that must be faced in a theoretical computation were also noted in these studies. Below a solution of the problem of stationary waves is obtained by the method of expansion in the small parameter in two regions with subsequent matching and also by a numerical integration method. It is shown that in each case a solution of solitary wave type exists along with the single-parameter family of periodic solutions (parameter—the wave number ). On decreasing the wave number, the periodic waves go over into a succession of solitary waves.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 63–66, January–February, 1977.The authors thank L. N. Maurin for helpful discussions and A. M. Tereshchenko for assisting in the computations.     

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.     

Critical layer formation in a stratified medium with mean shear flows

The problem of the excitation of internal waves with a given wave number k and frequency in a stratified medium with shear flows is considered. The internal wave field of the form v(z)exp(–it+ikx) established as t in a medium without dissipation has a singular point at the level z=z₀ (critical level), at which the flow velocity U(z) coincides with the phase velocity /k. Dissipative effects (viscosity and heat conduction) smooth out this singularity. An exact solution of the model equation describing as t and zz₀ the field excited by oscillating sources activated at t=0 is constructed with allowance for dissipation. This makes it possible to describe the limiting steady-state field, determine the critical layer as the neighborhood of the critical level in which dissipation effects are important, and to estimate its width and the rate of convergence to the limiting steady-state regime. The asymptotic behavior of the fields is examined for Ri1, where Ri is the Richardson number. It is shown that when the well-known Miles stability condition Ri>1/4 is satisfied there are no natural oscillations with a critical level.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 82–93, May–June, 1990.     

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.     

Asymptotic behavior of solutions of the Orr-Sommerfeld equations in regions adjacent to two branches of the neutral curve

The asymptotic behavior of the upper and lower branches of the neutral stability curve of a boundary layer found by Lin [1] was determined more accurately by various authors [2–4], who, on the basis of the linearized Navien-Stokes equations, analyzed the higher approximations in the Reynolds number R. In the limit R , neutral perturbations have wavelengths that exceed in order of magnitude the boundary layer thickness. The long-wavelength asymptotic behavior of the Orr-Sommerfeld equation is, in particular, of interest because the characteristic solutions of the linearized equations of free interaction (triple-deck theory) [5–7] are a limiting form of Tollmierr-Schlichting waves in an incompressible fluid with critical layers next to the wall [8–9]. At the same time, the dispersion relation, which is identical to the secular equation of the Orr-Sommerfeld problem, contains an entire spectrum of solutions not considered in the earlier studies [2–4]. The first oscillation mode in the spectrum may be either stable or unstable. In the present paper, solutions are constructed for each of the subregions (including the critical layer) into which the perturbed velocity field in the linear stability problem is divided at large Reynolds numbers. Dispersion relations describing the neighborhood of the upper and lower branches of the neutral curve for the boundary layer are derived. These relations, which contain neutral solutions as a special case, go over asymptotically into each other in the unstable region between the two branches.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 3–11, July–August, 1984.     

Magnetic hydrodynamics of two-dimensional flows in a plasma with an anisotropic pressure

The hydrodynamic equations of Chew, Goldberger, and Low [1] are used to analyze certain types of two-dimensional flows of a plasma with an anisotropic pressure (the pressure along the magnetic field p differs from the pressure across it p). In Sec. 1 the relationships derived in [2] for the transition of plasma state across surfaces of strong discontinuity are invoked to investigate the variation of the hydrodynamic parameters in weak shock waves in the linear approximation. The flow around bodies which only slightly perturb the main flow is investigated in Sec. 2 in the linear approximation. Similar problems for the case of an isotropic pressure are studied in detail in [3–5], for example.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 3–10, November–December, 1970.     

Differential equation for a stratified fluid and its particular solutions

The plane steady motion of a stratified ideal incompressible fluid in a gravity field is examined. Considering that the parameter characterizing the fluid particles — their density ₀ — is constant along the streamline, it is convenient to take the stream function as one of the independent variables and, in view of the presence of the gravity force, the Cartesian coordinate as the other. In this study, on the basis of Lavrent'eva's equation [1, 2, 3], the differential equation is derived for a single unknown function — the vertical displacement of the streamline y(y₀, x), where y₀ is its initial position and x is the horizontal coordinate. The particular solutions corresponding to a wave guide, cnoidal and solitary waves and, moreover, waves of the type corresponding to a smooth ascent to a new level are presented. A similar coordinate system was used in [4], but there the problem was reduced to a system of partial differential equations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 83–87, September–October, 1986.The authors are grateful to A. A. Barmin for discussing their results.     

Stability of convective motion caused by inhomogeneously distributed internal heat sources

The stability of steady convective plane-parallel flow in a vertical layer of viscous incompressible liquid of thickness h is investigated. The motion is caused by heat sources distributed in the liquid with volume density Q = Q₀exp (^–x) (the x axis is taken perpendicular to the boundary layer). The region of instability is determined for various values of the Prandtl number and the parameter N = h characterizing the inhomogeneity of the internal sources. It is shown that with increase in N there is qualitative rearrangement of the stability limit for perturbations of hydrodynamic type and incremental thermal waves.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 140–144, May–June, 1977.     

Stationary travelling waves on a film flowing down an inclined plane

At small flow rates, the study of long-wavelength perturbations reduces to the solution of an approximate nonlinear equation that describes the change in the film thickness [1–3]. Steady waves can be obtained analytically only for values of the wave numbers close to the wave number _n that is neutral in accordance with the linear theory [1, 2]. Periodic solutions were constructed numerically for the finite interval of wave numbers 0.5_{n n} in [4]. In the present paper, these solutions are found in almost the complete range of wave numbers 0 _n that are unstable in the linear theory. In particular, soliton solutions of this equation are obtained. The results were partly published in [5].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 142–146, July–August, 1980.     

Doubly periodic and quasiperiodic wave regimes in a liquid film flowing down an inclined plane,their stability and bifurcations

A nonlinear evolution equation frequently encountered in modeling the behavior of disturbances in various nonconservative media, for example, in problems of the hydrodynamics of liquid film flow, is considered. Wave solutions of this equation, regular in space and both periodic and quasiperiodic in time, branching off from steady and steady-state traveling waves are found numerically. The stability and bifurcations are analyzed for some of the solutions obtained. As a result, a bifurcation chain is found for solutions stable with respect to disturbances of the same spatial period. It is shown that the bifurcations are related to the loss of certain symmetries of the initial solution. It is demonstrated that as the bifurcation parameter increases it is possible to distinguish in the structure of the solutions intervals of quiet behavior and intervals of intense outbursts.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 98–107, July–August, 1992.     

Loss of stability and bifurcation of auto-oscillatory regimes close to Poiseuille flow

Bifurcation of Poiseuille flow in a flat channel is used as an example to analyze the problem of determining variables that permit study of bifurcation of a main steady flow of a viscous incompressible liquid for parameters close to the values of the coordinates of a point on the curve of neutral stability at which the first Lyapunov exponent d₀ vanishes and there is a changeover from subcritical to supercritical bifurcation. For Poiseuille flow, such a point (R₂,₂, where R₂ is the Reynolds number, and ₂ is the wave number, occurs on the lower branch of the neutral curve. In this paper, it is shown by the Lyapunov-Schmidt method that for < ₂ the stable time-periodic solution that bifurcates into the subcritical region loses stability in the case of slight supercriticality, and a fold singularity is formed in the amplitude surface. The nature of this additional bifurcation is determined by the sign of the second Lyapunov exponent d₁. For its calculation, the value of ₂ is fixed, and the bifurcation that occurs when the Reynolds number is changed is considered. A solution is sought in the form of a convergent series in powers of = ((R – R₀)^1/4, = ±1. The condition of solvability, which serves to determine the coefficient of ⁴, makes it possible to determine the value of d₂. This procedure is entirely general and makes it possible to study bifurcation in the neighborhood of a point of degeneracy on the neutral curve in other hydrodynamic problems too.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 41–48, March–April, 1991.     

Gas flow in cylindrical capillaries under intermediate conditions

At present, there are sufficient solutions of the problem of free-molecular gas flow through a short cylindrical channel, for example, [1–3]. In intermediate flow conditions, for Knudsen number Kn 1, solutions have been obtained for the limiting cases: an infinitely long channel [4] and a channel of zero length (an aperture) [5]. However, no solution is known for short channels for Kn 1. The present work reports a calculation by the Monte Carlo method of the macroscopic characteristics of the gas flow through a short cylindrical channel (for various length—radius ratios), taking into account intermolecular collisions.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 187–190, January–February, 1977.     

Stability of flows of the Hamel type in channels with permeable walls

The stability of nonparallel flows of a viscous incompressible fluid in an expanding channel with permeable walls is studied. The fluid is supplied to the channel through the walls with a constant velocity v₀ and through the entrance cross section, where a Hamel velocity profile is assigned. The resulting flow in the channel depends on the ratio of flow rates of the mixing streams. This flow was studied through the solution of the Navier—Stokes equations by the finite-difference method. It is shown that for strong enough injection of fluid through the permeable walls and at a distance from the initial cross section of the channel the flow approaches the vortical flow of an ideal fluid studied in [1]. The steady-state solutions obtained were studied for stability in a linear approximation using a modified Orr—Sommerfeld equation in which the nonparallel nature of the flow and of the channel walls were taken into account. Such an approach to the study of the stability of nonparallel flows was used in [2] for self-similar Berman flow in a channel and in [3] for non-self-similar flows obtained through a numerical solution of the Navier—Stokes equations. The critical parameters *, R*, and Cr* at the point of loss of stability are presented as functions of the Reynolds number R₀, characterizing the injection of fluid through the walls, and the parameter , characterizing the type of Hamel flow. A comparison is made with the results of [4] on the stability of Hamel flows with R₀ = 0.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 125–129, November–December, 1977.The author thanks G.I. Petrov for a discussion of the results of the work at a seminar at the Institute of Mechanics of Moscow State University.     

Taylor-Couette stability analysis for a Doi-Edwards fluid

The stability of Taylor-Couette flow of entangled polymeric solutions to small axisymmetric stationary disturbances is analyzed using the Doi-Edwards constitutive equation in the small gap limit. A previous analysis of Karlsson, Sokolov, and Tanner for the general K-BKZ equation, of which the Doi-Edwards equation is a special case, reduces the problem to one of numerically evaluating seven viscoelastic functions of the shear rate in the gap. Of these seven, only three — two of which are related to the second normal stress difference, and one of them to shear thinning — significantly affect the flow stability. The negative second normal stress difference of the Doi-Edwards fluid stabilizes the flow at low values of the Weissenberg number ₁ , while shear thinning produces strong destabilization at moderate Weissenberg number. Here ₁ is the longest relaxation time. Non-monotonic effects of viscoelasticity on Taylor-Couette stability analogous to those predicted here have been observed in experiments of Giesekus. The extreme shear thinning of the Doi-Edwards fluid is also predicted to produce a large growth in the height of the Taylor cells, a phenomenon that has been seen experimentally by Beavers and Joseph.     

Boundary layer transition on the surface of a delta wing in supersonic flow

Laminar-turbulent transition on the surface of a delta wing has been experimentally investigated in a supersonic wind tunnel at Mach numbers M_t8=3–5. It is shown that when M,=3, Re_L=6.5·10⁶, and =–5.5° much of the upper surface of the wing in the neighborhood of the line of symmetry is occupied by a wedge-shaped region of turbulent flow. In this region the heat fluxes reach the same values as at the heat transfer maxima induced here by separated flows and may exceed the turbulent heat flux level on the windward surface of the wing. Changing the shape of the under surface of the wing from plane to pyramidal leads to acceleration of the boundary layer transition on the under surface.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 87–92, May–June, 1989.     

Kinetic cooling of a moving gas

A study is made of the effects that arise when a moving gas absorbs electromagnetic radiation whose frequency is in resonance with the frequency of the center of a spectral line of a vibrational-rotational transition of molecules of the mixture. It is shown that the variation of the gas-dynamic parameters depends on the relationships between the rates of the stimulated transitions, intramolecular V — V exchange, and V — T relaxation, and the maximal effects are attained in the neighborhood of the sonic point.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 127–138, May–June, 1982.     

Shock heating of plasma in a rapidly increasing magnetic field

The experimental excitation of intense collisionless shock waves (M 5) with subsequent plasma compression by the magnetic field of a shock coil is described. A magnetic plug > 20 kOe is produced in 100 × 10^–9 sec by a current generator, a long line with 250-kV water insulation and a characteristic impedance of l At an initial deuterium-plasma density of 2 × 10¹⁴ cm^–3, shock waves with a front width of 20c/₀₃and a velocity of 5 × 10₇ cm/sec are recorded. The ion energy after the accumulation, determined from the neutron yield, turns out to be 2 ke V. Axial shock waves excited by the plasma flow beneath the shock coil are observed.Translated from Zhurnal Prikladnoi Mekhaniki i Teknicheskoi Fiziki, Vol. 11, No. 2, pp. 28–38, March–April, 1970.The authors thank G. I. Budker and R. Z. Sagdeev for formulating the problem, R. I. Soloukhin for interest in the study, and S. P. Shalamov for construction of the apparatus.     

Messung des ersten Normalspannungskoeffizienten von Kunststoffschmelzen im Bogenspalt

Zusammenfassung Bei einer stationären Schichtenströmung in einem Bogenspalt (azimutale Druckströmung im Ringspalt) bildet sich zwischen Innen- und Außenwand eine Druckdifferenz aus, deren Größe ein Maß für den 1. Normalspannungskoeffizienten der elastischen Flüssigkeit im Spalt ist. Die Strömung läßt sich zur Messung des 1. Normalspannungskoeffizienten verwenden. Der Schergeschwindigkeitsbereich der Messung liegt, wie bei der Kapillarrheometrie zur Bestimmung der Viskosität, zwischen 1 und 1000 s^–1. Die Auswertung der Messungen ist wegen des inhomogenen Scherfeldes relativ kompliziert. In der Arbeit wird ein besonders wirkungsvolles numerisches Auswerteverfahren hergeleitet und auf bestehende Messungen angewendet. Eine Besonderheit des Auswerteverfahrens ist die Freiheit der Wahl des Approximationsansatzes für die Viskositätskurve, während analytische Verfahren meist an einen bestimmten Ansatz gebunden sind. Außerdem braucht, im Gegensatz zu anderen derartigen Verfahren, die Position des schubspannungsfreien Stromfadensr ₀ nicht bestimmt zu werden.

---

Summary The stress in steady viscometric flow of molten polymers is determined by the viscosity and by the two normal stress coefficients ₁ and ₂. The paper describes a method of measuring ₁ by means of steady circumferential shear flow in an annulus. The cylinders are stationary and the fluid flows due to a circumferential pressure gradient. The radial normal stresses at the outer and at the inner wall are different from each other. The pressure-differencep is a measure for the 1. normal stress coefficient of the viscoelastic fluid. Due to the inhomogeneous shear field, the evaluation of ₁ fromp measurements is quite complicated. A powerful numerical method of evaluation has been developed and applied to existing data. The method is not restricted to a special empirical formula for the flow curve (as an analytical method would be) and does not require the knowledge of the positionr ₀ of the stress-free stream line.

---

a Pa s² Stoffparameter des Ansatzes des 1. Normalspannungskoeffizienten, s. Gl. [8] - AR _i — Koeffizient des Druckgefälles in-Richtung (Programm PFEIL) - AU _i — Koeffizient für Integration nach Simpson-Regel (Programm PFEIL) - b s² Stoffparameter des Ansatzes des 1. Normalspannungskoeffizienten - B _i — Koeffizient auf der rechten Seite des linearen Gleichungssystems (Programm PFEIL) - c — Exponent des Ansatzes des 1. Normalspannungskoeffizienten - CL _i CM _i CR _i — Koeffizienten der dimensionslosen Geschwindigkeit in dem linearen Gleichungssystem (Programm PFEIL) - F ₁,F ₂,F ₃ — Ableitungen der Summe der Fehlerquadrate nacha, b undc - G _k — Gewichtsfaktor - h m Spaltweite,r _a –r _i - H — dimensionslose Spaltweite, (r _a –r _i )/r _a - l m Länge des Bogenspaltes, 0,75(r _a +r _i ) - m — Exponent des Potenzansatzes, s. Gl. [13] - n — Dämpfungskonstante - N ₁ Pa 1. Normalspannungsdifferenz, _rr – - N ₂ Pa 2. Normalspannungsdifferenz - p Pa Druck - p Pa Druckgradient in-Richtung - P — dimensionsloser Druckgradient in-Richtung, s. Gl. [14] - p, p _k Pa Normalspannungsdifferenz zwischen Innen- und Außenwand im Bogenspalt, (– p + _rr ) _a – (–p + _rr ) _i - Q — Summe der Fehlerquadrate - r, R= r/r _a m, — Radiusvektor (Koordinate in Gradientenrichtung) - r ₀,R ₀=r ₀/r _a m, — Radius des neutralen Fadens - R — dimensionslose radiale Schrittweite - T, °C Temperatur bzw. Bezugstemperatur - v ms^–1 Geschwindigkeitskomponente in-Richtung - V ,V _,i — dimensionslose Geschwindigkeitskomponente in-Richtung - V _a ,V _k — dimensionslose Geschwindigkeit an der Außen- bzw. Innenwand - v _r ,v _z ms^–1 Geschwindigkeitskomponenten inr-undz-Richtung - ms ^–1 mittlere Geschwindigkeit in-Richtung - z m Koordinate in der indifferenten Richtung - K^–1 Temperaturkoeffizient der Viskosität - s^–1 Schergeschwindigkeit - s^–1 kritische Schergeschwindigkeit der Viskositätskurve, s. Gl. [13] - s^–1 Bezugsschergeschwindigkeit, - — dimensionslose Schergeschwindigkeit - — dimensionslose kritische Schergeschwindigkeit, - Pa s Viskosität - ₀ Pa s Nullviskosität - Pa s Bezugsviskosität, - — Radienverhältnis,r _i /r _a - ₁ Pa s ² 1. Normalspannungskoeffizient - Pa s² mittlerer 1. Normalspannungskoeffizient - ₂ Pa s² 2. Normalspannungskoeffizient - — Koordinate in Strömungsrichtung - Pa Spannung - a an der Außenwand - i, an der Innenwand - i laufender Index inr-Richtung - k Nummer des Meßpunktes - n Anzahl der Meßpunkte - n _i nord für Programm PFEIL - s _i süd für Programm PFEIL Mit 9 Abbildungen und 2 Tabellen     

Unsteady flows of an erosion plasma in irradiation of solid targets by annular beams

In recent times high-pressure physics has made ever wider application of constructions which use convergent shock waves [1–8]. The study of gas dynamic flows with convergent shock waves imposes the need for more careful calculation of the motions of a gas in regions whose dimensions are much less than the characteristic dimensions of the flow. In the present study the numerical method is used to study the gas dynamic phenomena accompanying the irradiation of solid obstacles by annular beams of monochromatic radiation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 179–182, November–December, 1988.In conclusion we note that at very short durations t t_k the solution to the problem is similar to the flow during separation of a gaseous toroid [19].     

Parametric excitation of waves at an interface

Experiments on the parametric excitation of waves at a fluid interface show a strong disagreement with theoretical results [1–3], since the latter do not take into account the influence of the second medium. This proves to be especially important at low frequencies. Thus, for a water-air interface with an excitation frequency = 60 sec^–1 the contribution amounts to 10%,and with = 30 sec^–1, even 20%. In this paper the stability of the interface of two viscous, incompressible fluids of finite depth in a variable gravity field is considered. The problem is put in the linear form by making an expansion with respect to the small viscosity and is solved by taking the Laplace transform with respect to time. A second-order integrodifferential equation with periodic coefficients is obtained for the deviation of the interface from the equilibrium position; its solution is sought by the method of averaging [4]. It is shown that the presence of the second fluid significantly raises the threshold of instability.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 167–170, March–April, 1977.