Irrotational outflow of liquid from a stream through a lateral hole of finite depth

An algorithm is constructed for numerical determination of the flow parameters and coefficient of contraction of a jet in the case of irrotational lateral outflow of liquid from a semiinfinite stream through a nozzle of finite depth situated at an arbitrary angle to the mainstream flow. The solution is based on the use of N. E. Zhukovskii's method and the Schwarz-Christoffel formula. The results of calculations for a nozzle situated at an angle = /2 ± , where = /6, are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 162–164, January–February, 1977.     

General expressions for unsteady forces in a lattice of profiles

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

Determination of permeability tensors for two-phase flow in homogeneous porous media: Theory

In this paper we continue previous studies of the closure problem for two-phase flow in homogeneous porous media, and we show how the closure problem can be transformed to a pair of Stokes-like boundary-value problems in terms of pressures that have units of length and velocities that have units of length squared. These are essentially geometrical boundary value problems that are used to calculate the four permeability tensors that appear in the volume averaged Stokes' equations. To determine the geometry associated with the closure problem, one needs to solve the physical problem; however, the closure problem can be solved using the same algorithm used to solve the physical problem, thus the entire procedure can be accomplished with a single numerical code.Nomenclature a a vector that maps V onto , m^-1. - A a tensor that maps V onto . - A area of the - interface contained within the macroscopic region, m². - A area of the -phase entrances and exits contained within the macroscopic region, m². - A area of the - interface contained within the averaging volume, m². - A area of the -phase entrances and exits contained within the averaging volume, m². - Bo Bond number (= (=(–)g²/). - Ca capillary number (= v/). - g gravitational acceleration, m/s². - H mean curvature, m^-1. - I unit tensor. - permeability tensor for the -phase, m². - viscous drag tensor that maps V onto V. - ^* dominant permeability tensor that maps onto v , m². - ^* coupling permeability tensor that maps onto v , m². - characteristic length scale for the -phase, m. - l characteristic length scale representing both and , m. - L characteristic length scale for volume averaged quantities, m. - n unit normal vector directed from the -phase toward the -phase. - n unit normal vector representing both n and n . - n unit normal vector representing both n and n . - P pressure in the -phase, N/m². - p superficial average pressure in the -phase, N/m². - p intrinsic average pressure in the -phase, N/m². - p –p , spatial deviation pressure for the -phase, N/m². - r ₀ radius of the averaging volume, m. - r position vector, m. - t time, s. - v fluid velocity in the -phase, m/s. - v superficial average velocity in the -phase, m/s. - v intrinsic average velocity in the -phase, m/s. - v –v , spatial deviation velocity in the -phase, m/s. - V volume of the -phase contained within the averaging volmue, m³. - averaging volume, m³. Greek Symbols V /, volume fraction of the -phase. - viscosity of the -phase, Ns/m². - density of the -phase, kg/m³. - surface tension, N/m. - (v +v ^T ), viscous stress tensor for the -phase, N/m².     

Effect of the Nature of the Interaction between a Gas and a Surface on the Relations for the Knudsen Layer in the Case of Strong Evaporation

The effect of the temperature accommodation coefficient _T on the relations at the Knudsen layer edge is investigated for strong evaporation using the moment method. An explicit expression for the dimensionless density as a function of the temperature and the Mach number M is obtained for 0 < _T < 1. For _T = 0 the entire solution is obtained in explicit form. It is shown that for = 0 and a condensation coefficient << 1 the temperature outside the Knudsen layer changes sharply as M varies from 0 to a certain value much less than unity after which the temperature ceases to depend on . For the model of specular reflection of the molecules from the surface the density and the temperature outside the Knudsen layer are found in explicit form as functions of the Mach number.     

Heteroclinic orbits in singular systems: A unifying approach

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

Thermodynamically complete equations of state for α-, ε-, and γ-iron

The numerical model of phase transition in iron in stress waves described in [1] contains equations of state with a limited range of applicability. They do not consider thermal excitation of conduction electrons and the presence of and — -triple point on the phase equilibrium curve, the effect of which should appear in shock loading of porous or preheated specimens. The present study will offer thermodynamically complete equations of state for the -, -, -phases of iron, free of these shortcomings.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 112–114, May–June, 1986.     

Flow with separation of an ideal fluid past a wing of infinite aspect ratio with a triangular protrusion

A study is made of the steady flow of an incompressible ideal fluid past a rectangular wing of infinite aspect ratio that has a protrusion on its leading edge. The protrusion is a triangular plate of small aspect ratio. It is assumed that flow separation occurs only from the side edges of the protrusion. An asymptotic solution to the problem is found for = O(), = o(1), where is the half-angle at the tip of the protrusion, and is the angle of attack. The results of numerical calculations are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 65–72, May–June, 1982.We are sincerely grateful to G. P. Svishchev for posing the problem and for interest in the work.     

Stability of stationary traveling waves on the surface of a vertical film of viscous fluid

The stability of stationary traveling waves of the first and second families with respect to infinitesimal perturbations of arbitrary wavelength is subjected to a detailed numerical investigation. The existence of a unique region of stability of the first family is established for wave numbers (₁, ₁) corresponding to the optimal wave regime. There are several regions of stability of the second family ( _k , _k),k=2,3,..., lying close to the local flow rate maxima. In the regions of instability the growth rates of perturbations of the first family are several times greater than for the second family. This difference increases with increase in the Reynolds number. The calculations make it possible to explain a number of experimental observations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 33–41, May–June, 1989.The authors are grateful to V. Ya. Shkadov for his constant interest, and to A. G. Kulikovskii, A. A. Barmin and their seminar participants for useful discussions and suggestions.     

Theoretical derivation of tangential velocity profiles in a flat vortex chamber-influence of turbulence and wall friction

Summary As part of a study on the hydrodynamics of a cyclone separator, a theoretical investigation of the flow pattern in a flat box cyclone (vortex chamber) has been carried out. Expressions have been derived for the tangential velocity profile as influenced by internal friction (eddy viscosity) and wall friction. The most important parameter controlling the tangential velocity profile is = –u ₀ R/(v+ ), where u ₀ is the radial velocity at the outer radius R of the cyclone, the kinematic liquid viscosity and is the kinematic eddy viscosity. For values of greater than about 10 the tangential velocity profile is nearly hyperbolic, for smaller than 1 the tangential velocity even decreases towards the centre. It is shown how and also the wall friction coefficient may be obtained from experimental velocity profiles with the aid of suitable graphs. Because of the close relation between eddy viscosity and eddy diffusion, measurements of velocity profiles in flat box cyclones will also provide information on the eddy motion of particles in a cyclone, a motion reducing its separation efficiency.List of symbols A cross-sectional area of cyclone inlet - h height of cyclone - p static pressure in cyclone - p static pressure difference in cyclone between two points on different radius - r radius in cyclone - r ₁ radius of cyclone outlet - R radius of cyclone circumference - u radial velocity in cyclone - u ₀ radial velocity at circumference of flat box cyclone - v tangential velocity - v ₀ tangential velocity at circumference of flat box cyclone - w axial velocity - z axial co-ordinate in cyclone - friction coefficient in flat box cyclone (for definition see § 5) - ₁ value of friction coefficient for ₁<< ₂ - ₂ value of friction coefficient for ₂<<1 - = - ₁ value of for ₁<< ₂ - ₂ value of for ₂<<1 - thickness of laminar boundary layer - =/h - turbulent kinematic viscosity - ratio of z to h - k ratio of height of cyclone to radius R of cyclone - parameter describing velocity profile in cyclone =–u ₀ R/(+) - kinematic viscosity of fluid - density of fluid - ratio of r to R - ₁ value of at outlet of cyclone - ₂ value of at inner radius of cyclone inlet - _w shear stress at cyclone wall - angular momentum in cyclone/angular momentum in cyclone inlet - ₁ value of at = ₁ - ₂ value of at = ₂     

On surface waves in generalized thermoelasticity

Knowles' representation theorem for harmonically time-dependent free surface waves on a homogeneous, isotropic elastic half-space is extended to include harmonically time-dependent free processes for thermoelastic surface waves in generalized thermoelasticity of Lord and Shulman and of Green and Lindsay.r , , r , , .This work was done when author was unemployed.     

On laminar flow through a uniformly porous pipe

Numerous investigations ([1] and [4–9]) have been made of laminar flow in a uniformly porous circular pipe with constant suction or injection applied at the wall. The object of this paper is to give a complete analysis of the numerical and theoretical solutions of this problem. It is shown that two solutions exist for all values of injection as well as the dual solutions for suction which had been noted by previous investigators. Analytical solutions are derived for large suction and injection; for large suction a viscous layer occurs at the wall while for large injection one solution has a viscous layer at the centre of the channel and the other has no viscous layer anywhere. Approximate analytic solutions are also given for small values of suction and injection.

Nomenclature

General r distance measured radially - z distance measured along axis of pipe - u velocity component in direction of z increasing - v velocity component in direction of r increasing - p pressure - density - coefficient of kinematic viscosity - a radius of pipe - V velocity of suction at the wall - r ²/a ² - R wall or suction Reynolds number, Va/ - f() similarity function defined in (6) - u ₀() eigensolution - U(0) a velocity at z=0 - K an arbitrary constant - B _K Bernoulli numbers Particular Section 5 perturbation parameter, –2/R - ² a constant, –K - x / - g(x) f()/ Section 6 perturbation parameter, –R/2 - ² a constant, –K - g() f() - g _c ()=g() near centre of pipe - * point where g()=0 Section 7 2/R - ² K - t (1–)/ - w(t, ) [1–f(t)]/ - ₀, ₁ constants - g() f()– ₀ - ₀/ - ₀ a constant - * point where f()=0     

Global Behavior of a Nonlinear Quasiperiodic Mathieu Equation

In this paper, we investigate the interaction of subharmonicresonances in the nonlinear quasiperiodic Mathieu equation,x + [ + (cos ₁ t + cos ₂ t)] x + x³ = 0.We assume that 1 and that the coefficient of the nonlinearterm, , is positive but not necessarily small.We utilize Lie transform perturbation theory with elliptic functions –rather than the usual trigonometric functions – to study subharmonic resonances associated with orbits in 2m:1 resonance with a respective driver. In particular, we derive analytic expressions that place conditions on (, , ₁, ₂) at which subharmonic resonance bands in a Poincaré section of action space begin to overlap. These results are used in combination with Chirikov's overlap criterion to obtain an overview of the O() global behavior of equation (1) as a function of and ₂ with ₁, , and fixed.     

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

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

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.     

Flow through a helically coiled annulus

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

Solutions for the anisotropic elastic wedge at critical wedge angles

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

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

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

Solution of regular beam equations in arbitrary emission conditions on a curvilinear surface

The regular beam equations are solved analytically for the case of emission from an arbitrary surface in conditions of total space charge (-mode) and in a given external magnetic field H (§2) for temperature-limited emission (T-mode), in an external magnetic field H (§3); and for emission with nonzero initial velocity (§4). The emitter is taken as the coordinate surface x¹=0 in an orthogonal system x¹ (i = =1,2,3), while the current density J and field on it are given functions j(x², x³), (x², x³. The solution is written as series in (x¹) with coefficients dependent on x², x³, determined from recurrence relations. For emission in the -mode and H 0, =1/3; for temperature-limited emission, =1/2; with nonzero initial velocity, =1. The results are extended to the case of a beam in the presence of a moving background of uniform density (5).     

Up-scaling Flow in Fractured Media: Equivalence Between the Large Scale Averaging Theory and the Continuous Time Random Walk Method

In a recent paper, (McNabb, 1978), we set up a method allowing to compute both the transient and steady-state exchange terms between the matrix and fractured regions of a naturally fractured porous medium using the continuous time random walk method (CTRW). In particular, the exchange coefficient parametrizing the large-scale exchange term was computed on physical grounds using a time integration of the so-called time correlation function corresponding to the particle presence in the fractures. On the other hand, the large scale averaging theory (LSAT) as developed by Quintard and Whitaker (Quintard and Whitaker, 1996) gives another definition for this exchange coefficient . It also provides a so-called closure problem allowing to compute from the solution of a well-defined steady state boundary value problem, to be solved over a representative volume of the high resolution fractured map. The goal of the present paper is to show analytically that both definitions coincide, yielding a unique and well defined value of the coefficient. This provides an unification of two approaches whose respective backgrounds are very different at the first glance.