On physical criteria for the selection of wood for soundboards of musical instruments

In order to develop criteria for the physical evaluation of wood for soundboards of musical instruments, measurements were made of dynamic Young's modulusE, static Young's modulusE, internal frictionQ ^–1 in longitudinal direction, and specific gravity for numerous species of broad-leaved wood. From the results obtained, including those of our previous paper on coniferous wood [1], it was found that the suitability of wood for soundboards could be evaluated by the quantity ofQ ^–1/(E/), and that there were very high correlations betweenQ ^–1/(E/) andE/, and betweenE andE, regardless of wood species. Consequently, it becomes possible to select practically any wood suitable for soundboards by using the value ofE/, which can be measured easily, and it was derived that the relation betweenE/ andQ ^–1 of wood could be expressed by an exponential equation regardless of wood species.     

Evaluation of the dynamic-mechanical properties of solids from a cooperative model for stress relaxation

For many solid materials the stress relaxation process obeys the universal relationF = – (d/d lnt)_max = (0.1 ± 0.01) ( ₀ — _i ), regardless of the structure of the material. Here denotes the stress,t the time, ₀ the initial stress of the experiment and _i the internal stress. A cooperative model accounting for the similarity in relaxation behaviour between different materials was developed earlier. Since this model has a spectral character, the concepts of linear viscoelasticity are used here to evaluate the corresponding prediction of the dynamic mechanical properties, i.e. the frequency dependence of the storageE () and lossE () moduli. Useful numerical approximations ofE () andE () are also evaluated. It is noted that the universal relation in stress relaxation had a counterpart in the frequency dependence ofE (). The theoretical prediction of the loss factor for high-density polyethylene is compared with experimental results. The agreement is good.     

The virtual origin of riblets in an open channel flow

In this paper, a method using the mean velocity profiles for the buffer layer was developed for the estimation of the virtual origin over a riblets surface in an open channel flow. First, the standardized profiles of the mixing length were estimated from the velocity measurement in the inner layer, and the location of the edge of the viscous layer was obtained. Then, the virtual origins were estimated by the best match between the measured velocity profile and the equations of the velocity profile derived from the mixing length profiles. It was made clear that the virtual origin and the thickness of the viscous layer are the function of the roughness Reynolds number. The drag variation coincided well with other results.Nomenclature f _r skin friction coefficient - f _ro skin friction coefficient in smooth channel at the same flow quantity and the same energy slope - g gravity acceleration - H water depth from virtual origin to water surface - H ⁺ u*H/ - H false water depth from top of riblets to water surface - H ⁺ u*H/ - I _e streamwise energy slope - I _b bed slope - k riblet height - k ⁺ u*k/ - l mixing length - l _s standardized mixing length - Q flow quantity - Re Reynolds number volume flow/unit width/v - s riblet spacing - u mean velocity - u* friction velocity = - u* false friction velocity = - y distance from virtual origin - y distance from top of riblet - y ₀ distance from top of riblet to virtual origin - y _v distance from top of riblet to edge of viscous layer - y ⁺ u*y/ - y ⁺ u*y/ - y ₀ ⁺ u*y ₀/ - u ⁺ u*y/ - shifting coefficient for standardization - thickness of viscous layer=y ₀+y - ⁺ u*/ - ⁺ u*/ - eddy viscosity - ridge angle - v kinematic viscosity - density - shear stress     

Mean and turbulent velocity measurements of supersonic mixing layers

The behavior of supersonic mixing layers under three conditions has been examined by schlieren photography and laser Doppler velocimetry. In the schlieren photographs, some large-scale, repetitive patterns were observed within the mixing layer; however, these structures do not appear to dominate the mixing layer character under the present flow conditions. It was found that higher levels of secondary freestream turbulence did not increase the peak turbulence intensity observed within the mixing layer, but slightly increased the growth rate. Higher levels of freestream turbulence also reduced the axial distance required for development of the mean velocity. At higher convective Mach numbers, the mixing layer growth rate was found to be smaller than that of an incompressible mixing layer at the same velocity and freestream density ratio. The increase in convective Mach number also caused a decrease in the turbulence intensity ( _u/U).List of symbols a speed of sound - b total mixing layer thickness between U ₁ – 0.1 U and U ₂ + 0.1 U - f normalized third moment of u-velocity, f u³/(U)³ - g normalized triple product of u² , g u²/(U)³ - h normalized triple product of u ², h u ²/(U)³ - l _u axial distance for similarity in the mean velocity - l _u axial distance for similarity in the turbulence intensity - M Mach number - M _c convective Mach number (for ₁ = ₂), M _c (U ₁ – U ₂)/(a ₁ + a ₂) - P static pressure - r freestream velocity ratio, r U ₂/U ₁ - Re unit Reynolds number, Re U/ - s freestream density ratio, s ₂/₁ - T _t total temperature - u instantaneous streamwise velocity - u deviation of u-velocity, uu – U - U local mean streamwise velocity - U ₁ primary freestream velocity - U ₂ secondary freestream velocity - average of freestream velocities, (U ₁ + U ₂)/2 - U freestream velocity difference, U U ₁ – U ₂ - instantaneous transverse velocity - v deviation of -velocity, – V - V local mean transverse velocity - x streamwise coordinate - y transverse coordinate - y ₀ transverse location of the mixing layer centerline - ensemble average - ratio of specific heats - boundary layer thickness (y-location at 99.5% of free-stream velocity) - similarity coordinate, (y – y ₀)/b - compressible boundary layer momentum thickness - viscosity - density - standard deviation - dimensionless velocity, (U – U ₂)/U - 1 primary stream - 2 secondary stream A version of this paper was presented at the 11th Symposium on Turbulence, October 17–19, 1988, University of Missouri-Rolla     

On the indentation of an inhomogeneous plastic solid

Summary The problem considered here is that of the indentation of a semi infinite, inhomogeneous rigid-plastic solid by a smooth, flat ended punch under conditions of plane strain. It is assumed that the yield stress of the solid k(x, y) has the form k ⁰+k(x, y) where k ⁰ is a constant and is small. A perturbation method of solution developed by Spencer [1] is used, and general results are obtained for arbitrary values of k(x, y). Some particular cases are then considered.     

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].     

Drag reduction by polymer solutions in a riblet-lined pipe

The flow of 3 to 100 wppm aqueous solutions of a polyethyleneoxide polymer,M _w=6.2×;10⁶, was studied in a 10.2 mm i.d. pipe lined with 0.15 mm V-groove riblets, at diametral Reynolds numbers from 300 to 150000. Measurements in the riblet pipe were accompanied by simultaneous measurements in a smooth pipe of the same diameter placed in tandem. The chosen conditions provided turbulent drag reductions from zero to the asymptotic maximum possible. The onset of polymer-induced drag reduction in the riblet pipe occurred at the same wall shear stress, * _w =0.65 N/m², as that in the smooth pipe. After onset, the polymer solutions in the riblet pipe initially exhibited linear segments on Prandtl-Karman coordinates, akin to those seen in the smooth pipe, with specific slope increment . The maximum drag reduction observed in the riblet pipe was independent of polymer concentration and well below the asymptotic maximum drag reduction observed in the smooth pipe. Polymer solution flows in the riblet pipe exhibited three regimes: (i) Hydraulically smooth, in which riblets induced no drag reduction, amid varying, and considerable, polymer-induced drag reduction; this regime extended to non-dimensional riblet heightsh ⁺<5 in solvent andh ⁺<10 in polymer solutions. (ii) Riblet drag reduction, in which riblet-induced flow enhancementR>0; this regime extended from 5<h ⁺<22 in solvent and from 10<h ⁺<30 in the 3 wppm polymer solution, with respective maximaR=0.6 ath ⁺=14 andR=1.6 ath ⁺=21. Riblet drag reduction decreased with increasing polymer concentration and increasing polymer-induced flow enhancement S. (iii) Riblet drag enhancement, whereinR<0; this regime extended for 22<h ⁺<110 in solvent, withR;–2 forh ⁺>70, and was observed in all polymer solutions at highh ⁺, the more so as polymer-induced drag reduction increased, withR<0 for allS>8. The greatest drag enhancement in polymer solutions,R=–7±1 ath ⁺=55 whereS=20, considerably exceeded that in solvent. Three-dimensional representations of riblet- and polymer-induced drag reductions versus turbulent flow parameters revealed a hitherto unknown dome region, 8<h ⁺<31, 0<S<10, 0<R<1.5, containing a broad maximum at (h ⁺,S,R) = (18, 5, 1.5). The existence of a dome was physically interpreted to suggest that riblets and polymers reduce drag by separate mechanisms.     

Bursting and structure of the turbulence in an internal flow manipulated by riblets

The buffer layer of an internal flow manipulated by riblets is investigated. The distributions of the ejection and bursting frequency from the beginning to the middle part of the buffer layer, together with high moments of the fluctuating streamwise velocity,u, and its time derivative are reported. The profiles of the ejection and bursting frequency are determined and compared using three single point detection schemes. The effect of the riblets on the bursting mechanism is found confined in a localized region in the buffer layer. The multiple ejection bursts are more affected than the single ejection bursts. The skewness and flatness factors of theu signal are larger in the manipulated layer than in the standard boundary layer. That, also holds true for the flatness factor of the time derivative, but the Taylor and Liepman scales are not affected. The spectrum of theu signal is altered at the beginning part of the viscous sublayer.Nomenclature u Friction velocity - Viscosity - l _v ;f _v wall scalesv/u ;u ² /v - y Vertical distance to the wall - z Spanwise extent - (⁺) Variable normalized with wall scales - u Velocity;u=Turbulence intensity - h, s Height and width of the riblets - f _e Ejection frequency - f _b Bursting frequency - f _BME Frequency of the Bursts with Multiple Ejection - f _BSE Frequency of Single Ejection Bursts - S andS _du/dt Skewness factor ofu and its time derivative - F _u andF _du/dt Flatness factor ofu and its time derivative Abbreviations SBL Standard (non-manipulated) Boundary Layer - MBL Manipulated Boundary Layer - BME Bursts with Multiple Ejections - BSE Bursts with Single Ejections - VITA Variable Interval Time Averaging technique - u–l u-level technique - mu Modifiedu-level technique     

Dynamic-mechanical properties of a bitumen-silica composite

Summary Dynamic-mechanical properties of bitumen-silica composite materials, measured at room temperature, do not vary with the volume ratio () in a simple manner as do usual bituminous concretes. However,E is a linear function of the interfacial area () between the filler and the binder per unit volume. ThusE = E ₀ +a(), wherea is a constant related to the storage modulus, in the absence of voids and with a void ratio factor. The loss moduli, plotted against () go through a maximum in a similar way as when plotted versus decreasing temperatures.

---

Zusammenfassung Die bei Zimmertemperatur gemessenen dynamischmechanischen Eigenschaften von Bitumen, welches mit Siliziumdioxid gefüllt worden ist, ändern sich nicht in einer so einfachen Weise mit dem Volumenanteil des Füllstoffs, wie es bei gewöhnlichen Asphaltbetonen der Fall ist. Dagegen ergibt sich der SpeichermodulE als eine lineare Funktion der Grenzfläche zwischen Bindemittel und Füllstoff pro Volumeneinheit:E = E ₀ +a(), wobeia eine Konstante bedeutet, die zum einen von dem Speichermodul bei Abwesenheit von Hohlräumen und zum andern von einem durch solche Hohlräume bedingten Faktor abhängt. Der Verlustmodul als Funktion von zeigt ein Maximum ähnlich wie bei der Auftragung gegen die Temperatur.

---

Résumé Les propriétés mécaniques dynamiques, mesurées à la température ambiante, d'un système composite bitume-silice, n'évoluent pas simplement avec la fraction volumique de charges comme dans le cas des enrobés usuels. En effet, le module de conservationE varie linéairement avec l'aire de l'interface liant-charge par unité de volume (). Ainsi la relation suivante a pu être mise en évidence:E = E ₀ +a(), la constantea étant fonction d'un module indépendant du taux de vide et d'un terme relié à ce dernier. Si augmente, les modules de perte passent par une valeur maximale. Ces variations sont semblables à celles que l'on aurait si l'on portait ces modules en fonction de la température.

---

---

With 9 figures and 2 tables     

The effect of non-steady wall shear on pressure surges in conduits carrying viscous liquids

Summary A study is made of the attenuation of pressure surges in a two-dimension a channel carrying a viscous liquid when a valve at the downstream end is suddenly closed. The analysis differs from previous work in this area by taking into account the transient nature of the wall shear, which in the past has been assumed as equivalent to that existing in steady flow. For large values of the frictional resistance parameter the transient wall shear analysis results in less attenuation than given by the steady wall shear assumption.Nomenclature c /, ft/sec - e base of natural logarithms - F(x, y) integration function, equation (38) - (x) mean value of F(x, y) - g local acceleration of gravity, ft/sec² - h width of conduit, ft - k (2m–1)^{2 2} L/h ² c, equation (35) - k* 12L/h ² c, frictional resistance parameter, equation (46) - L length of conduit, ft - m positive integer - n positive integer - p pressure, lb/ft² - p ₀ constant pressure at inlet of conduit, lb/ft² - P pressure plus elevation head, p+gz, equation (4) - mean value of P over the conduit width h - P ₀ p ₀+gz ₀, lbs/ft² - R frictional resistance coefficient for steady state wall shear, lb sec/ft⁴ - s positive integer; also, condensation, equation (6) - t time, sec - t ct/L, dimensionless time - u x component of fluid velocity, ft/sec - u _m mean velocity in conduit, equation (12), ft/sec - u ₀(y) velocity profile in Poiseuille flow, equation (19), ft/sec - transformed velocity - U initial mean velocity in conduit - x distance along conduit, measured from valve (fig. 1), ft - x x/L, dimensionless distance - y distance normal to conduit wall (fig. 1), ft - y y/h, equation (25) - z elevation, measured from arbitrary datum, ft - z ₀ elevation of constant pressure source, ft - isothermal bulk compression modulus, lbs/ft² - _n , equation (37) - _n (2n–1)/2, equation (36) - viscosity, slugs/ft sec - / = kinematic viscosity, ft²/sec - density of fluid, slugs/ft³ - ₀ density of undisturbed fluid, slugs/ft³ - ø angle between conduit and vertical (fig. 1) The research upon which this paper is based was supported by a grant from the National Science Foundation.     

Effects of MHD free convection and mass transfer on the flow past an oscillating infinite coaxial vertical circular cylinder

The effects of MHD free convection and mass transfer are taken into account on the flow past oscillating infinite coaxial vertical circular cylinder. The analytical expressions for velocity, temperature and concentration of the fluid are obtained by using perturbation technique.

---

Einwirkungen von freier MHD-Konvektion und Stoffübertragung auf eine Strömung nach einem schwingenden unendlichen koaxialen vertikalen Zylinder

Zusammenfassung Die Einwirkungen der freien MHD-Konvektion und Stoffübertragung auf eine Strömung nach einem schwingenden, unendlichen, koaxialen, vertikalen Zylinder wurden untersucht. Die analytischen Ausdrücke der Geschwindigkeit, Temperatur und Fluidkonzentration sind durch die Perturbationstechnik erhalten worden.

---

Nomenclature C _p Specific heat at constant temperature - C the species concentration near the circular cylinder - C _w the species concentration of the circular cylinder - C the species concentration of the fluid at infinite - * dimensionless species concentration - D chemical molecular diffusivity - g acceleration due to gravity - Gr Grashof number - Gm modified Grashof number - K thermal conductivity - Pr Prandtl number - r _a ,r _b radius of inner and outer cylinder - a, b dimensionless inner and outer radius - r,r coordinate and dimensionless coordinate normal to the circular cylinder - Sc Schmidt number - t time - t dimensionless time - T temperature of the fluid near the circular cylinder - T _w temperature of the circular cylinder - T temperature of the fluid at infinite - u velocity of the fluid - u dimensionless velocity of the fluid - U ₀ reference velocity - z,z coordinate and dimensionless coordinate along the circular cylinder - coefficient of volume expansion - * coefficient of thermal expansion with concentration - dimensionless temperature - H ₀ magnetic field intensity - coefficient of viscosity - _e permeability (magnetic) - kinematic viscosity - electric conductivity - density - M Hartmann number - dimensionless skin-friction - frequency - dimensionless frequency     

Certain yawed magnetogasdynamic flows

Certain steady yawed magnetogasdynamic flows, in which the magnetic field is everywhere parallel to the velocity field, are related to certain reduced three-dimensional compressible gas flows having zero magnetic field. Under a restriction, the reduced flows are linked, by certain reciprocal relations, to a four parameter class of plane gas flows. In the instance of constant entropy an approximation method is suggested for obtaining magnetogasdynamic flows from the corresponding plane, irrotational gasdynamic flows and examples are given.

Nomenclature

magnetogasdynamic flow variables H magnetic intensity - q fluid velocity - fluid density - p pressure - s entropy - Q _t, H _t component of q, H in the x–y plane - w , h component of q, H perpendicular to the x–y plane reduced gasdynamic flow factor of proportionality - q* fluid velocity - * fluid density - p* pressure - Q _t ^* =u*î+v*, w* components of q* - l arbitrary constant - A _v Alfvén speed - Q _t, , p fluid velocity, density, pressure of the reciprocal gas dynamic flow - L, n, k, arbitrary constants - , velocity potential, stream function - angle made by Q _t, Q _t ^* , and V with the x-axis - adiabatic gas constant - a ²=(–1)/2 constant - M Mach number - W constant value of w* - E approximate constant value of g(p) - * modified potential function - modified velocity coordinate - +i - complex potential of the irrotational flow - B arbitrary constant - V incompressible flow velocity - V modified fluid velocity - X _p, Y _p points on the profile     

Flow of soft solids squeezed between planar and spherical surfaces

The force to squeeze a Herschel–Bulkley material without slip between two approaching surfaces of various curvature is calculated. The Herschel–Bulkley yield stress requires an infinite force to make plane–plane and plane–concave surfaces touch. However, for plane–convex surfaces this force is finite, which suggests experiments to access the mesoscopic thickness region (1–100 m) of non-Newtonian materials using squeeze flow between a plate and a convex lens. Compared to the plane–parallel surfaces that are used most often for squeeze flow, the dependence of the separation h and approach speed V on the squeezing-time is more complicated. However, when the surfaces become close, a simplification occurs and the near-contact approach speed is found to vary as V h⁰ if the Herschel–Bulkley index is n<1/3, and V h^(3n-1)/(2n) if n 1/3. Using both plane–plane and plane–convex surfaces, concordant measurements are made of the Herschel–Bulkley index n and yield stress ₀ for two soft solids. Good agreement is also found between ₀ measured by the vane and by each squeeze-flow method. However, one of the materials shows a limiting separation and a V(h) behaviour not predicted by theory for h<10 m, possibly owing to an interparticle structure of similar lengthscale.     

Symplectic approach to the theory of quantized fields. II

We prove that the set D of vector fields on the configuration space B of a field whose 1-parameter groups locally associated are groups of fibre-preserving transformations of B that leave invariant that field in the sense of variational theory, is a Lie algebra with respect to ordinary addition, multiplication by real numbers and Lie brackets. We see that this Lie algebra structure can be carried over to the corresponding set of Noether invariants, which then becomes a Lie algebra in a natural way.Further, we define the n-form of Poincaré-Cartan of a field, and we use it to generalize the Lie algebras D and in a reasonable way. The algebras D and are subalgebras of the new Lie algebras D and introduced. A main result in this connection is the following: the differential d of the n-form of Poincaré-Cartan is –(d+f), where (, d+f) are the field equations on the vertical bundle B.The symplectic manifold of solutions associated with a field is introduced in a formal way and the former Lie algebras D, , D, are interpreted on this manifold. In imitation of the case of analytical dynamics, the main results in this direction are: a) Every vector field of the Lie algebra D defines, in a canonical way, a vector field on the manifold of solutions such that its polar 1-form with respect to the symplectic metric ₂ is the differential of its corresponding Noether invariant, and b) the Lie bracket [, ] of two Noether invariants , is the Noether invariant given by ₂(D, D), where D, D are the vector fields on the manifold of solutions defined, in the sense a), by two infinitesimal generators of , , respectively. This will allow us to regard the Lie algebra as the analogous object in field theory to the Poisson algebra of analytic dynamics.We apply the general formalism to the relativistic theory of non-linear scalar fields, and we compare our results with the formalism developed by I. Segal for this case.     

Distribution of eddy viscosity and mixing length in a non-Newtonian flow of a solid particle-gas suspension

The present work concerns the study of the radial distribution of eddy viscosity and mixing length and their dependence upon the Reynolds number and the concentration of the solid phase in a non-Newtonian flow of a suspension of solid particles in a gas. The investigated systems have a pseudoplastic character and the deviation from Newtonian behaviour increases with an increase in the concentration of the dispersed phase. Relations are presented for the eddy viscosity and mixing length in the flow of pseudoplastic fluids. From the analysis of results it follows that the mixing length and eddy viscosity increases with an increase in the Reynolds number. In contrast, an increase in the concentration of the solid phase and consequently of the pseudoplasticity causes a decrease in the investigated quantities. The radial distribution of the mixing length and the eddy viscosity is characterized by a maximum, after which the investigated quantities vary only slightly. This enables the area of the core of the turbulent flow to be defined. E Non-dimensional eddy viscosity - K fluid consistency defined by Ostwald-De Waele's formula (power law) - K fluid consistency, eq. (12) - L mixing length - L _t non-dimensional mixing length - N ⁺ position parameter, eq. (3) - n power-law index - n ⁺ Reichardt's position parameter, eq. (4) - n slope of the dependence ln _w = f[ln (8w/D)] - R pipe radius - r radial distance from the pipe axis - Re=D W /µ Reynolds number - U ⁺ non-dimensional local mean axial velocity, eq. (2) - u ^* = w (/8)^0.5 friction velocity - non-dimensional local mean axial velocity - local mean axial velocity - u turbulence velocity component - mean axial velocity at pipe axis - w average velocity over cross-section of pipe - loading ratio of solid to air, i.e. ratio of mass flow rates of solids to air - Y _m ⁺ =Rⁿ u^2–n /K 8^n–1 non-dimensional distance of pipe centre from the wall - y distance from pipe wall - y ⁺ non-dimensional distance from pipe wall, eq. (5) - friction factor - µ _L laminar part of viscosity - µ _t eddy viscosity - density - shear stress - _w shear stress on the wall     

Vibration of rectangular orthotropic plates

Summary An approximate analytical procedure has been given to solve the problem of a vibrating rectangular orthotropic plate, with various combinations of simply supported and clamped boundary conditions. Numerical results have been given for the case of a clamped square plate.Nomenclature 2a, 2b sides of the rectangular plate - h plate thickness - E _x , E _y , E, G elastic constants of te orthotropic material - D _x E _x h ³/12 - D _y E _y h ³/12 - H _xy Eh ³/12+Gh ³/6 D _x , D _y and H _xy are rigidity constants of the orthotropic plate - mass of the plate per unit area - Poisson's ratio - W deflection of the plate - p circular frequency - b/a ratio - X _m , Y _n characteristic functions of the vibrating beam problem - p ² a ² b ²/H _xy the frequency parameter.     

New concepts in relative permeabilities at high capillary numbers for surfactant flooding

Flooding oil reservoirs with surfactant solutions can increase the amount of oil that can be recovered. Macroscopic modelling of the process requires relative permeabilities to be functions of saturation and capillary number. With only limited experimental data, relative permeabilities have usually been assumed to be linear functions of saturation at high capillary numbers. The experimental data is reviewed, some of which suggest that this assumption is not necessarily correct. The basis for the assumption is therefore reviewed and it is concluded that the linear model corresponds to microscopically segregated flow in the porous medium. Based on new but equally plausible complementary assumptions about the flow pattern, a mixed flow model is derived. These models are then shown to be limiting cases of a droplet model which represents the mixing scale within the porous medium and gives a physical basis for interpolating between the models. The models are based on physical concepts of flow in a porous medium and so the approach described here represents a significant improvement in the understanding of high capillary number flow. This is shown by the fact that fewer parameters are needed to describe experimental data.Notation A total cross-sectional area assigned to capillary bundle - A ⁽ⁱ⁾ physical cross-sectional area of tube i - c ⁽ⁱ⁾ ordered configurational label for droplets in tube i - c configuration label for tube i (order not considered) - D defined by Equation (26) - E(...) expectation value with respect to the trinomial distribution - S _r ⁽⁾ fractional flow of phase - k absolute permeability - k _r relative permeability of phase - k _r ⁰ endpoint relative permeability of phase - L capillary tube length in bundle model - m ⁽ⁱ⁾ number of droplets of phase a occupying tube i - n exponent for phase a in Equation (2) - N number of droplets in bundle model - N _c capillary number - p pressure - p(c') probability of configuration c - Q ⁽ⁱ⁾ total volume flow rate in tube i - S saturation of phase - S flowing saturation of phase - S _r residual saturation of phase - S _r ⁽⁾ saturations when fractional flow of phase is 1 in the case of varying residual saturations for three-phase flow ( ) - t _c residence time for droplet configuration c - v ⁽ⁱ⁾ total fluid velocity in bundle tube i - , phase label - p pressure differential across capillary bundle - ⁽ⁱ⁾ tube conductivity defined by Equation (7) - viscosity of phase - interfacial tension - gradient operator - ... average over tube droplet configurations     

The dynamic performance of the Weissenberg rheogoniometer III. Large amplitude oscillatory shearing — harmonic analysis

The harmonic content of the nonlinear dynamic behaviour of 1% polyacrylamide in 50% glycerol/water was studied using a standard Model R 18 Weissenberg Rheogoniometer. The Fourier analysis of the Oscillation Input and Torsion Head motions was performed using a Digital Transfer Function Analyser.In the absence of fluid inertia effects and when the amplitude of the (fundamental) Oscillation Input motion _I is much greater than the amplitudes of the Fourier components of the Torsion Head motion _Tn empirical nonlinear dynamic rheological propertiesG _n (, ⁰),G _n (, ⁰) and/or _n (, ⁰), _n (, ⁰) may be evaluated without a-priori-knowledge of a rheological constitutive equation. A detailed derivation of the basic equations involved is presented.Cone and plate data for the third harmonic storage modulus (dynamic rigidity)G ₃ (, ⁰), loss modulusG ₃ (, ⁰) and loss angle ₃ (, ⁰) are presented for the frequency range 3.14 × 10^–2 1.25 × 10² rad/s at two strain amplitudes, _CP ⁰ = 2.27 and 4.03. Composite cone and plate and parallel plates data for both the third and fifth harmonic dynamic viscosities ₃ (, ⁰), _S (, ⁰) and dynamic rigiditiesG ₃ (, ⁰),G ₅ (, ⁰) are presented for strain amplitudes in the ranges 1.10 _CP ⁰ 4.03 and 1.80 _PP ⁰ 36 for a single frequency, = 3.14 × 10^–1 rad/s. Good agreement was obtained between the results from both geometries and the absence of significant fluid inertia effects was confirmed by the superposition of the data for different gap widths.     

On the localness of the spectral energy transfer in turbulence

By utilizing available experimental data for net energy transfer spectra for homogeneous turbulence, contributions P(, ) to the energy transfer at a wavenumber from various other wavenumbers are calculated. This is done by fitting a truncated power-exponential series in and to the experimental data for the net energy transfer T(), and using known properties of P(, ). Although the contributions P(, ) obtained by using this procedure are not unique, the results obtained by using various assumptions do not differ significantly. It seems clear from the results that for a region where the energy entering a wavenumber band dominates that leaving, much of the energy entering the band comes from wavenumbers which are about an order of magnitude smaller. That is, the energy transfer is rather nonlocal. This result is not significantly dependent on Reynolds number (for turbulence Reynolds numbers based on microscale from 3 to 800). For lower wavenumbers, where more energy leaves than enters a wavenumber band, the energy transfer into the band is more local, but much of the energy then leaves at distant wavenumbers.     

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.