Generation of vortices in a partially filled,rapidly rotating cylinder

We describe a system in which vortices are shed from a cylindrical free surface approximately centered in a rotating flow. Shedding is controlled by the parameter =2 g/ ² d, where g, , d denote gravity, rotation rate and the diameter of the free surface. We find vortex shedding for >0.162 and no vortex shedding for < 0.0847. The range depends on the aspect ratio L/d, where L is the column length, in a nonmonotonic fashion. These results are independent of viscosity and surface tension for small values of these parameters.Now at Martin Marietta, Orlando Aerospace, PO Box 5837, Mail Point 150, Orlando, FL 32855, USA     

An analytical approach to coupled heat and moisture transport in soil

The simultaneous diffusion of heat and moisture through soil is described by two coupled partial differential equations in which the diffusion coefficients are highly non-linear functions of the dependent variables. The system has been regarded as analytically intractable for any generality of coupled flow. However, for an asymptotically steady state, the equations show a marked periodic stability. Computer simulation indicates that the behaviour quickly becomes entrained to input boundary periodicity for any initial state, regardless of the detailed functional form of the diffusion coefficients. This property allows an harmonic series solution to be assembled. Factors such as amplitude decay, phase shift and wave form evolution may be evaluated. The solution is adapted to boundary conditions pertaining to arid soils and the results validated against the 1968 field data of Rose and the 1973 experiment by Jackson.Notation gradient operator - divergence operator - _A amplitude of surface moisture content variation - _l volumetric liquid content, m³/m³ - _c value for moisture content, at which vapour diffusivity decays to zero - _M mean of surface moisture content variation - _s saturation value of moisture content - tortuosity factor, m/m - _i eigenvalues of ₀ - hypothetical thermal conductivity, J/m/sec/K - ₀ density of saturated water vapour, kg/m³ - _l density of liquid water, kg/m³ - _v density of water vapour, kg/m³ - surface tension, kg/sec² - matric potential, m - C volumetric heat capacity, J/m³/K - D ^* molecular diffusivity of water vapour in the porous medium, m²/sec - D _atm molecular diffusivity of water vapour in air, m²/sec - D _TV thermally induced vapour diffusivity, m²/sec/K - D _Tl thermally induced liquid diffusivity, m²/sec/K - D _v isothermal vapour diffusivity, m²/sec - D _l isothermal liquid diffusivity, m²/sec - L latent heat of vaporisation, J/kg - P atmospheric pressure at soil surface,Pa - R gas constant of water vapour, J/kg/K - T temperature,K - T _M mean temperature at surface, K - T _A temperature amplitude at surface, K - g acceleration due to gravity, m/sec² - h relative humidity, dimensionless - p partial pressure of water vapour,Pa - q _v water vapour flux, kg/m²/sec - t time, sec - z depth, (measured downwards), m     

Cryogenic transonic wind tunnels and the condensation of nitrogen

A brief tutorial on the need for Reynolds number similarity and the advent of cryogenic transonic wind tunnels is presented. Experimental results of nitrogen condensation in nozzles are collected and related to the flow in the wind tunnels. New theoretical approaches to a solution of the condensation problem in the supersaturated state are proposed.List of symbols a speed of sound - A area - {ovc} wing area/wing span - c _p pressure coefficient, Eq. (12) - G ^* energy of formation of a critical droplet, Eq. (14) - h altitude - J homogeneous nucleation rate, Eq. (13) - k Boltzmann constant - l characteristic length - M Mach number, Eq. (2) - n ^* number of molecules in a critical droplet, Eq. (14) - p static pressure - p ₀ wind tunnel supply pressure - p ⁰ standard pressure - p equilibrium vapor pressure - P wind tunnel fan power - q dynamic pressure - Re Reynolds number, Eq. (1) - t time - T temperature - T ₀ wind tunnel supply temperature - molecular volume - V air speed - ratio of specific heats - dynamic viscosity - v kinematic viscosity - density - surface tension This paper is dedicated to my old friend Eberhard Berger upon his retirement from the Föttinger Institut of the Technical University of Berlin     

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

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

Zur Kapillarviskosimetrie verdünnter Lösungen

Zusammenfassung Es werden allgemeine Dimensionierungsregeln für Kapillarviskosimeter angegeben. Der Einfluß der Oberflächenspannung () der zu untersuchenden Flüssigkeit auf die mittlere hydrostatische Druckhöhe ist für den Fall berechnet, daß das Gefäß zwischen den beiden Meßmarken eine Kugel mit dem Volumen von 50 cm³ ist. Die berechnete mittlere Druckhöhe ist eine lineare Funktionvon /. Daher ist es im Fall desUbbelohde-Viskosimeters möglich, den Kapillarzug durch eine entsprechende Krümmung des hängenden Niveaus im Ausfluß der Kapillare zu kompensieren.

---

Summary General rules to dimension capillary viscosimeters are given. The effect of the surface tension () of the fluid under test on the mean hydrostatic head is calculated for a sphere of 50 cm³ between the measuring marks.The calculated mean hydrostatic head is a linear function of/. Therefore it is possible in the case of theUbbelohde-viscosimeter to compensate the capillary traction by appropriate curvature of the suspended level in the outflow of the capillary.

     

Mixed convection flow in narrow vertical ducts

The mixed convection flow in a vertical duct is analysed under the assumption that , the ratio of the duct width to the length over which the wall is heated, is small. It is assumed that a fully developed Poiseuille flow has already been set up in the duct before heat from the wall causes this to be changed by the action of the buoyancy forces, as measured by a buoyancy parameter . An analytical solution is derived for the case when the Reynolds numberRe, based on the duct width, is of 0 (1). This is extended to the case whenRe is 0 (^–1) by numerical integrations of the governing equations for a range of values of representing both aiding and opposing flows. The limiting cases, || 1 andR=Re of 0 (1), andR and both large, with of 0 (R ^1/3) are considered further. Finally, the free convection limit, large with R of 0 (1), is discussed.

---

Mischkonvektion in engen senkrechten Rohren

Zusammenfassung Mischkonvektion in einem senkrechten Rohr wird unter der Voraussetzung untersucht, daß das Verhältnis der Rohrbreite zur Länge, über welche die Wand beheizt wird, klein ist. Es wird angenommen, daß sich bereits eine voll entwickelte Poiseuille-Strömung in dem Rohr eingestellt hat, bevor Antriebskräfte, gemessen mit dem Auftriebsparameter , aufgrund der Wandbeheizung die Strömung verändern. Es wird eine analytische Lösung für den Fall erhalten, daß die mit der Rohrbreite als charakteristische Länge gebildete Reynolds-ZahlRe konstant ist. Dies wird mittels einer numerischen Integration der wichtigsten Gleichungen auf den FallRe =f (^–1) sowohl für Gleich- als auch für Gegenstrom ausgedehnt. Weiterhin werden die beiden Grenzfälle betrachtet, wenn || 1 undR=Re konstant ist, sowieR und beide groß mit proportionalR ^1/3. Schließlich wird der Grenzfall der freien Konvektion, großes mit konstantem R, diskutiert.

---

Nomenclature g acceleration due to gravity - Gr Grashof number - G modified Grashof number - h duct width - l length of the heated section of the duct wall - p pressure - Pr Prandtl number - Q flow rate through the duct - Q ₀ heat transfer on the wally=0 - Q ₁ heat transfer on the wally=1 - Re Reynolds number - R modified Reynolds number - T temperature of the fluid - T ₀ ambient temperature - T applied temperature difference - u, velocity component in thex-direction - v, velocity component in they-direction - x, co-ordinate measuring distance along the duct - y, co-ordinate measuring distance across the duct - buoyancy parameter - ₀ modified buoyancy parameter, ₀=R ^–1/3 - coefficient of thermal expansion - ratio of duct width to heated length, =h/l - (non-dimensional) temperature - _w applied temperature on the wally=0 - kinematic viscosity - density of the fluid - ₀ shear stress on the wally=0 - ₁ shear stress on the wally=1 - stream function     

Transport in ordered and disordered porous media I: The cellular average and the use of weighting functions

In this work we consider transport in ordered and disordered porous media using singlephase flow in rigid porous mediaas an example. We defineorder anddisorder in terms of geometrical integrals that arise naturally in the method of volume averaging, and we show that dependent variables for ordered media must generally be defined in terms of thecellular average. The cellular average can be constructed by means of a weighting function, thus transport processes in both ordered and disordered media can be treated with a single theory based on weighted averages. Part I provides some basic ideas associated with ordered and disordered media, weighted averages, and the theory of distributions. In Part II a generalized averaging procedure is presented and in Part III the closure problem is developed and the theory is compared with experiment. Parts IV and V provide some geometrical results for computer generated porous media.Roman Letters A interfacial area of the- interface contained within the macroscopic region, m² - A_e area of entrances and exits for the-phase contained within the macroscopic system, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor, m² - L general characteristic length for volume averaged quantities, m - characteristic length (pore scale) for the-phase - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - N unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p0 reference pressure in the-phase, N/m² - p traditional intrinsic volume averaged pressure, N/m² - r0 radius of a spherical averaging volume, m - r position vector, m - r position vector locating points in the-phase, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - v velocity vector in the-phase, m/s - v traditional superficial volume averaged velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V/V, volume average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

Motion of large gas bubbles ascending in a liquid

An equation is derived for the ascent velocity of large gas bubbles in a liquid. This velocity is assumed to be governed by the propagation of a wavelike perturbation caused by the bubble in the liquid.Notation w bubble (or drop) velocity - specific gravity - dynamic viscosity - kinematic viscosity - r bubble (or drop) radius - surface tension - coefficient of friction - g gravitational acceleration - D bubble (or drop) diameter - p pressure - c propagation velocity of the wavelike perturbation - wavelength     

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.     

Evaluating Herschel-Bulkley fluids with the back extrusion (annular pumping) technique

A mathematical model was developed to describe the behavior of Herschel-Bulkley fluids in a back extrusion (annular pumping) device. A technique was also developed to determine the rheological properties (yield stress, flow behavior index, and consistency coefficient) of these fluids. Mathematical terms were expressed in four dimensionless terms, and graphical aids and tables were prepared to facilitate the handling of the expressions.Nomenclature a radius of the plunger, m - dv/dr shear rate, s^–1 - F force applied to the plunger, N - F _b buoyancy force, N - F _cb force corrected for buoyancy, N - F _T recorded force just before the plunger is stopped, N - F _Te recorded force after the plunger is stopped, N - g acceleration due to gravity, m/s² - H(t) momentary height between plunger and container bottom, m - K a/R, dimensionless - L length of annular region, m - L(t) depth of plunger penetration, m - n flow behavior index, dimensionless - p static pressure, Pa - P _L pressure in excess of hydrostatic pressure at the plunger base, Pa - p ₀ pressure at entrance to annulus, Pa - P pressure drop per unit of length, Pa/m - Q total volumetric flow rate through the annulus, m³/s - r radial coordinate, measured from common axis of cylinder forming annulus, m - R radius of outer cylinder of annulus, m - s reciprocal of n, dimensionless - t time, s - T dimensionless shear stress, defined in Eq. (3) - T ₀ dimensionless yield stress, defined in Eq. (4) - T _w dimensionless shear stress at the plunger wall - _p velocity of plunger, m/s - velocity, m/s - mass density of fluid, kg/m³ - Newtonian viscosity, Pa s - P p ₀ –p _L , Pa - consistency coefficient, Pa sⁿ - value of where shear stress is zero - _–, ₊ limits of the plug flow region (Fig. 1) - r/R - shear stress, Pa - _y yield stress, Pa - _w shear stress at the plunger wall, Pa - dimensionless flow rate defined in Eq. (24) - dimensionless velocity defined by Eq. (5) - _–, ₊ dimensionless velocity outside the plug flow region - _max dimensionless maximum velocity in the plug flow region - _p dimensionless velocity at the plunger wall     

Viscoelastic properties of semidilute poly(methyl methacrylate) solutions

Viscoelastic properties were examined for semidilute solutions of poly(methyl methacrylate) (PMMA) and polystyrene (PS) in chlorinated biphenyl. The number of entanglement per molecule, N, was evaluated from the plateau modulus, G _N . Two time constants, _s and ₁, respectively, characterizing the glass-to-rubber transition and terminal flow regions, were evaluated from the complex modulus and the relaxation modulus. A time constant _k supposedly characterizing the shrink of an extended chain, was evaluated from the relaxation modulus at finite strains. The ratios ₁/ _s and _k / _s were determined solely by N for each polymer species. The ratio ₁/ _s was proportional to N ^4.5 and N ^3.5 for PMMA and PS solutions, respectively. The ratio _k / _s was approximately proportional to N ^2.0 in accord with the prediction of the tube model theory, for either of the polymers. However, the values for PMMA were about four times as large as those for PS. The result is contrary to the expectation from the tube model theory that the viscoelasticity of a polymeric system, with given molecular weight and concentration, is determined if two material constants _s and G _N are known.     

A generalized Norton-Hoff model and the Prandtl-Reuss law of plasticity

The aim of this article is to study the quasistatic evolution of a three-dimensional elastic-perfectly plastic solid which satisfies the Prandtl-Reuss law. The evolution of the field of stresses -which solves a time dependent variational inequality — and that of the field of displacements u, have been described in previous works [15], [26], [35], [36], [37] but it was not shown there that and u satisfy indeed the Prandtl-Reuss constitutive law. In this article we find and u in a class of functions which are sufficiently regular for the Prandtl-Reuss law to make sense and we prove that and u satisfy the constitutive law. This result is attained by considering the elastic-perfectly plastic model as the limit of a family of elastic-visco-plastic models like those of Norton and Hoff. The Norton-Hoff type models which we introduce depend on a viscosity parameter > 0; we study the perturbed models (i.e. > 0 fixed) and then we pass to the limit 0.Dedicated to James Serrin on the occasion of his 60^th Birthday     

Linear shear flow past a porous particle

Linear shear flow past a porous spherical particle is studied using a generalized boundary condition proposed by Jones. The torque on a porous sphere rotating in a quiescent fluid is calculated. Streamlines patterns are illustrated for the case of a particle freely suspended in a simple shear flow. These patterns are shown to differ significantly from those associated with an impermeable rigid sphere. Finally, an expression for the effective viscosity of a dilute suspension of porous spherical particles is obtained.Nomenclature A, B dimensionless flow parameter - a radius of the porous sphere - C, E, F constants of integration - d shear strength - d constant rate of deformation of ambient field - e rate of strain tensor - f, g functions of distance - k permeability of the porous medium - n unit normal vector - p pressure - p unit vector - Q coefficient of spherical harmonic - q filter velocity within the porous medium - r polar spherical coordinate - S _p surface of porous particle - S, T, T* coefficients of spherical harmonics - T torque exerted on the particle - u fluid velocity vector - x cartesian coordinates - dimensionless constant - , polar spherical coordinates - dimensionless flow parameter - viscosity of the fluid - stress tensor - rotational velocity of the particle - rotational velocity of the ambient field.     

Instantaneous density measurement in two-dimensional gas flow by high speed differential interferometry

The aim of this paper is to present an experimental set-up using a Wollaston prism differential interferometer producing up to twenty successive short exposure white light interferograms at a high framing rate. It is shown that, through optical component calibration, the interferograms can be analysed to yield the instantaneous density field. This method has been successfully tested in the two-dimensional unsteady flow generated by the interaction of a mixing layer and a cavity.List of symbols h height of the downstream edge of the cavity - H height of backward facing step - M Mach number - t time - t time interval between two successive frames - N frequency - double-prism median plane - birefringence angle - _p pressure fluctuation - C _p pressure coefficient - biprism abscissa corresponding to any colour - ₀ biprism reference abscissa corresponding to background colour - _y deviation of light rays - R radius of curvature of spherical mirror - L virtual distance from the middle of the test section to the spherical mirror - E optical thickness - E _e optical thickness corresponding to background colour - d _E difference of optical thickness - d _x abscissa difference - gas density - ₀ stagnation gas density - _e gas density of background colour     

Flow in vicinity of blunt body stagnation point in diverging hypersonic stream

In the hypersonic thin shock layer approximation for a small ratio k of the densities before and after the normal shock wave the solution of [1] for the vicinity of the stagnation point of a smooth blunt body is extended to the case of nonuniform outer flow. It is shown that the effect of this nonuniformity can be taken into account with the aid of the effective shock wave radius of curvature R_*, whose introduction makes it possible to reduce to universal relations the data for different nonuniform outer flows with practically the same similarity criterion k. The results of the study are compared with numerical calculations of highly underexpanded jet flow past a sphere.Notations x, y a curvilinear coordinate system with axes directed respectively along and normal to the body surface with origin at the forward stagnation point - R radius of curvature of the meridional plane of the body surface - uV, vV., , p V ² respectively the velocity projections on the x, y axes, density, and pressure - and V freestream density and velocity The indices =0 and=1 apply to plane and axisymmetric flows Izv. AN SSSR, Mekhanika Zhidkosti i Gaza, Vol. 5, No. 3, pp. 102–105, 1970.     

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     

Absolute determination of the isochromatic parameter by load-stepping photoelasticity

A new approach to full-field automated photoelasticity is presented in which a circular polariscope is used to enable the isochromatic phase value () to be determined unambiguously and without input of a known isochromatic value obtained using an auxiliary technique. Values of cos are obtained from light-field and dark-field images for three loads of small incremental steps. Using a relatively straight-forward procedure, ramped phase maps for are produced which can be unwrapped using conventional techniques. The resulting distribution of is then found absolutely using information provided by which is the incremental change in the isochromatic phase value between the load steps. The results obtained for disk-in-compression tests presented here in comparison with theoretical solutions demonstrate that the technique is both simple to use and very accurate. A similar approach may be adopted using three wavelengths instead of three load steps.     

Natural convection heat transfer through an enclosed horizontal layer of supercritical carbon dioxide

In natural convection heat transfer through a thin horizontal layer of carbon dioxide, maxima in the equivalent thermal conductivities are obtained in the vicinity of the respective pseudocritical temperatures at pressures of 75.8, 89.6 and 103.4 bar. The maxima are the more pronounced, the closer the critical point is approached.Comparison of experimental results with Nusselt equations shows good agreement except for the immediate vicinity of the pseudocritical temperature.In visual observations a distinct change in flow structure appears in the immediate vicinity of the pseudocritical temperature. A steady state polygon pattern and a boiling-like action could not be observed in this geometry.

---

Zusammenfassung Beim Wärmetransport durch freie Konvektion in einer dünnen waagerechten Schicht von Kohlendioxid ergaben sich Maxima der scheinbaren Wärmeleitfähigkeit in der Nähe der pseudokritischen Temperaturen bei Drükken von 75,8, 89,6 und 103,4 bar. Die Maxima sind um so ausgeprägter, je mehr man sich dem kritischen Punkt nähert.Ein Vergleich der Versuchsergebnisse mit Nusseltbeziehungen ergibt gute Übereinstimmung außer in unmittelbarer Umgebung der pseudokritischen Temperatur. Direkte Beobachtungen der Konvektionsmuster zeigen in unmittelbarer Umgebung der pseudokritischen Temperatur eine deutliche Strukturänderung. Ein stationäres Zellmuster und siedeähnliche Vorgänge konnten in dieser Anordnung nicht beobachtet werden.

---

Nomenclature A area of the heating or cooling plate - C constant in the correlation - g acceleration of gravity - h heat transfer coefficient - k thermal conductivity of fluid in the gap - k _e equivalent thermal conductivity - m, n exponents of dimensionless numbers - q heat flux - T _C,PC absolute temperature; critical C, pseudocritical PC - Gr Grashof numberg ( _h– _c) ³/ ² - Nu Nusselt numberh/k - Pr Prandtl number/ - thermal diffusivity - coefficient of volume expansion - width of gap - _c,h temperature of cooling (c)-, heating (h)-plate - _m arithmetic mean temperature ( _c+ _h)/2 - kinematic viscosity - _c,h fluid density at the temperature of the cooling (c)- or heating (h)-plate - heat flow rate through the gap     

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.     

Dielectric properties of heterogeneous mixtures with a polar constituent

Summary After defining the boundaries for the dielectric constant of a heterogeneous mixture, the behaviour of such a mixture is studied as a function of the frequency, when one of its components is polar. Deviations from a semicircle are to be expected for the function _m =f( _m ) even when the dielectric properties of the polar constituent can be described with a semicircular Cole-Cole-arc. The relaxation time of the mixture is shorter than that of the polar constituent.