Effect of radial and rotational lubricant inertia on the pressure distribution in externally pressurised thrust bearings

Expressions are obtained for the pressure distribution in an externally pressurised thrust bearing for the condition when one bearing surface is rotated. The influence of centripetal acceleration and the combined effect of rotational and radial inertia terms are included in the analysis. Rotation of the bearing causes the lubricant to have a velocity component in an axial direction towards the rotating surface as it spirals radially outwards between the bearing surfaces. This results in an increase in the pumping losses and a decrease in the load capacity of the bearing. A further loss in the performance of the bearing is found when the radial inertia term, in addition to the rotational inertia term is included in the analysis.Nomenclature r, z, cylindrical co-ordinates - V _r, V , V _z velocity components in the r, and z directions respectively - U, X, W representative velocities - coefficient of viscosity - p static pressure at radius r - p mean static pressure at radius r - Q volume flow per unit time - 2h lubricant film thickness - density of the lubricant - r ₂ outside radius of bearing = D/2 - angular velocity of bearing - R dimensionless radius = r/h - P dimensionless pressure = h ³ p/Q - Re channel Reynolds number = Q/h     

Mixed convection on vertical plate for fluids of any prandtl number

A mixed convection parameter=(Ra) ^1/4/(Re)^1/2, with=Pr/(1+Pr) and=Pr/(1 +Pr)^1/2, is proposed to replace the conventional Richardson number, Gr/Re², for combined forced and free convection flow on an isothermal vertical plate. This parameter can readily be reduced to the controlling parameters for the relative importance of the forced and the free convection,Ra ^1/4/(Re ^1/2 Pr ^1/3) forPr 1, and (RaPr)^1/2/(RePr ^1/2 forPr 1. Furthermore, new coordinates and dependent variables are properly defined in terms of, so that the transformed nonsimilar boundary-layer equations give numerical solutions that are uniformly valid over the entire range of mixed convection intensity from forced convection limit to free convection limit for fluids of any Prandtl number from 0.001 to 10,000. The effects of mixed convection intensity and the Prandtl number on the velocity profiles, the temperature profiles, the wall friction, and the heat transfer rate are illustrated for both cases of buoyancy assisting and opposing flow conditions.

---

Mischkonvektion an einer vertikalen Platte für Fluide beliebiger Prandtl-Zahl

Zusammenfassung Für die kombinierte Zwangs- und freie Konvektion an einer isothermen senkrechten Platte wird ein Mischkonvektions-Parameter=( Ra) ^1/4 (Re)^1/2, mit=Pr/(1 +Pr) und=Pr/(1 +Pr)^1/2 vorgeschlagen, den die gebräuchliche Richardson-Zahl, Gr/Re², ersetzen soll. Dieser Parameter kann ohne weiteres auf die maßgebenden Kennzahlen für den relativen Einfluß der erzwungenen und der freien Konvektion reduziert werden,Ra ^1/4/(Re ^1/2 Pr ^1/3) fürPr 1 und (RaPr)^1/4/(RePr)^1/2 fürPr 1. Weiterhin werden neue Koordinaten und abhängige Variablen als Funktion von definiert, so daß für die transformierten Grenzschichtgleichungen numerische Lösungen erstellt werden können, die über den gesamten Bereich der Mischkonvektion, von der freien Konvektion bis zur Zwangskonvektion, für Fluide jeglicher Prandtl-Zahl von 0.001 bis 10.000 gleichmäßig gültig sind. Der Einfluß der Intensität der Mischkonvektion und der Prandtl-Zahl auf die Geschwindigkeitsprofile, die Temperaturprofile, die Wandreibung und den Wärmeübergangskoeffizienten werden für die beiden Fälle der Strömung in und entgegengesetzt zur Schwerkraftrichtung dargestellt.

---

Nomenclature C _f local friction coefficient - C _p specific heat capacity - f reduced stream function - g gravitational acceleration - Gr local Grashoff number,g T _w –T )x³/v² - Nu local Nusselt number - Pr Prandtl number,v/ - Ra local Rayleigh number,g T _w –T x ³/( v) - Re local Reynolds number,u x/v - Ri Richardson number,Gr/Re ² - T fluid temperature - T _w wall temperature - T free stream temperature - u velocity component in thex direction - u free stream velocity - v velocity component in they direction - x vertical coordinate measuring from the leading edge - y horizontal coordinate Greek symbols thermal diffusivity - thermal expansion coefficient - mixed convection parameter (Ra)1/4/Re)^1/2 - pseudo-similarity variable,(y/x) - ₀ conventional similarity variable,(y/x)Re ^1/2 - dimensionless temperature, (T–T T _W –T - unified mixed-flow parameter, [(Re) ^1/2 + (Ra)^1/4] - dynamic viscosity - kinematic viscosity - stretched streamwise coordinate or mixed convection parameter, [1 + (Re)^1/2/(Ra) ^1/4]^–1=/(1 +) - density - Pr/(1 + Pr) _w wall shear stress - stream function - Pr/(l+Pr)^1/3 This research was supported by a grand from the National Science Council of ROC     

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

A nonequilibrium theory of a slurry

A nonequilibrium theory of a slurry is developed and its practical use is illustrated by a simple stability analysis. Here a slurry is defined as a deformable continuum consisting of a liquid phase containing in suspension a large number of small solid particles which have formed by solidification from the liquid. The liquid is assumed to consist of two components and the solid to contain only one of the two. Consequently, the process of change of phase requires redistribution of material on the scale of the solid particles. This process is assumed to take a finite amount of time, requiring a nonequilibrium macroscopic theory. This theory contains four thermodynamic variables, three to represent the equilibrium state of the binary system and a fourth measuring the departure from thermodynamic equilibrium. The process of microscale diffusion of material is parameterized in the macroscale theory, leading to a Landau-type relaxation term in the equation of evolution of the fourth variable. The theory is simplified to yield a Boussinesq-like set of governing equations. Their practical use is illustrated by analyzing the stability of a simple steady solution of the equations and the effects of a non-zero relaxation time are discussed. A novel instability mechanism involving sedimentation of particles, previously found to occur in the equilibrium case, is found to persist in nonequilibrium, but disappears in the limit of no change of phase.Key to symbols a, b, c thermodynamic coefficients; see (3.36)–(3.38) - sedimentation coefficient; see (5.18) - C _p specific heat; see (3.24) - C _p ^de specific heat of the slurry; see (3.28) and (3.30) - c radius of solid particle (in §4) - D, D diffusive coefficients; see (3.40) and (3.41) - material diffusivity in liquid phase - D ^* modified diffusion coefficient; see (5.15) - d thermodynamic coefficient; see (3.39) - E specific internal energy - f, g, h thermodynamic coefficients; see (3.36)–(3.38) - g acceleration of gravity - reduced gravity; see (5.10) - i total diffusive flux vector of constituent 1 - i diffusive flux vector of constituent 1 in the liquid phase - j diffusive flux vector of solid phase - k thermal conductivity - k entropy flux vector - k _T, k_T thermodiffusion coefficients; see (3.40) and (3.41) - L latent heat of solidification per unit mass; see (3.7) and (3.24) - m wave number - m ^s rate of creation of mass of solid per unit volume through solidification - m ₁ ^s rate of creation of mass of solid constituent 1 per unit volume through solidification - mass rate of freezing per unit area per unit time - N number of solid particles per unit volume - p pressure - p _H hydrostatic component of pressure - p _m mechanical pressure - p ₁ dynamic component of pressure - q heat flux vector - Q _D rate of regeneration of heat through diffusive fluxes - Q _M rate of regeneration of heat through phase-change processes - Q _v rate of regeneration of heat through viscosity - Q vector defined by (3.16) - r heat externally supplied per unit mass (in §3); spherical radial coordinate (in §4) - S specific entropy of slurry - change of specific entropy with mass fraction of constituent 1; also change of chemical potential of liquid phase with temperature barring change of phase - change of chemical potential of liquid phase with temperature in phase equilibrium; see (3.28) and (3.30) - T temperature - t time - t ₀ relaxation time; see (5.30) - u barycentric velocity - u _H horizontal perturbation velocity - V sedimentation speed - w _a upward speed of simple state; see (6.5) and (6.12) - z upward vertical coordinate - upward unit vector - thermal expansion coefficient barring change of phase; see (3.23) - > ^* thermal expansion coefficient in phase equilibrium; see (3.27) and (3.30) - modified thermal expansion coefficient; see (5.1) and (5.4) - isothermal compressibility of slurry barring change of phase; see (3.23) - ^* isothermal compressibility of slurry in phase equilibrium; see (3.27) and (3.30) - dimensionless measure of departure from liquidus equilibrium; see (5.2) - _a deviation from phase equilibrium in simple state; see (6.6) and (6.13) - vertical wave number - volume expansion per unit mass upon melting; see (3.6) - change of chemical potential of liquid phase with pressure; see (3.25) - change of chemical potential of liquid phase with pressure for slurry; see (3.29) and (3.30) - compositional gradient in the static state; see (6.15) - vector defined by (3.35) - constant of integration; see (6.7) and (6.8) - coefficient defined by (6.23) - nonequilibrium expansion coefficient; see (5.1) and (5.4) - thermal diffusivity; =k/C _p - modified thermal diffusivity; see (5.33) - relaxation rate to phase equilibrium; see (2.2) - ₁ relaxation rate to solid-composition equilibrium; see (2.3) - sedimentation coefficient; see (4.29) - horizontal wave number vector - sedimentation coefficient; see (4.30) - ^L , ^s chemical potential of constituent 1 relative to constituent 2 in liquid and solid phase per unit mass; see (2.6) - change of chemical potential of liquid with liquid composition; see (3.8) - coefficient defined by (3.10) - kinematic shear viscosity - total mass fraction of constituent 1 (i.e., solute) - ^L, ^s mass fraction of constituent 1 in liquid and solid phases - density of slurry - _s density of solid phase - - - , growth rate of disturbance - stress tensor - deviatoric stress tensor - dimensionless temperature; see (5,3) - _a constant of integration; see (6.7) - mass fraction of solid phase in slurry - _b vertical gradient of mass fraction of solid; see (6.1) - dimensionless measure of _b; see (6.22) - _c temporal gradient of mass fraction of solid; see (6.1) - specific Gibbs free energy; see (3.13) - _L,_s specific Gibbs free energy of liquid and solid phases; see (2.12) - measure of departure from liquidus equilibrium; see (2.14) - measure of departure from solidus equilibrium; see (2.5) - spherical polar coordinate (in §4); see (4.20); wave angle (in §6); see (6.38)     

Eine analytische Näherungslösung für die eingefrorene laminare Grenzschichtströmung eines binären Gemisches

Zusammenfassung Für die eingefrorene laminare Grenzschichtströmung eines teilweise dissoziierten binären Gemisches entlang einer stark gekühlten ebenen Platte wird eine analytische Näherungslösung angegeben. Danach läßt sich die Wandkonzentration als universelle Funktion der Damköhler-Zahl der Oberflächenreaktion angeben. Für das analytisch darstellbare Konzentrationsprofil stellt die Damköhler-Zahl den Formparameter dar. Die Wärmestromdichte an der Wand bestehend aus einem Wärmeleitungs- und einem Diffusionsanteil wird angegeben und diskutiert. Das Verhältnis beider Anteile läßt sich bei gegebenen Randbedingungen als Funktion der Damköhler-Zahl ausdrücken.

---

An analytical approximation for the frozen laminar boundary layer flow of a binary mixture

An analytical approximation is derived for the frozen laminar boundary layer flow of a partially dissociated binary mixture along a strongly cooled flat plate. The concentration at the wall is shown to be a universal function of the Damkohler-number for the wall reaction. The Damkohlernumber also serves as a parameter of shape for the concentration profile which is presented in analytical form. The heat transfer at the wall depending on a conduction and a diffusion flux is derived and discussed. The ratio of these fluxes is expressed as a function of the Damkohler-number if the boundary conditions are known.

---

Formelzeichen A Atom - A₂ Molekül - C Konstante in Gl. (20) - c₁=1/(2C) Konstante in Gl. (35) - c_p spezifische Wärme bei konstantem Druck - D binärer Diffusionskoeffizient - Ec=u ² /(2h_f) Eckert-Zahl - h spezifische Enthalpie - h_t=h+u²/2 totale spezifische Enthalpie - h _A ⁰ spezifische Dissoziationsenthalpie - K_w Reaktionsgeschwindigkeitskonstante der heterogenen Wandreaktion - 1= /( ) Champman-Rubesin-Parameter - Le=Pr/Sc Lewis-Zahl - M Molmasse - p statischer Druck - Pr= c_pf/ Prandtl-Zahl - q_w Wärmestromdichte an der Wand - q_cw, q_dw Wärmeleitungsbzw. Diffusionsanteil der Wärmestromdichte an der Wand - universelle Gaskonstante - R=/(2M_a) individuelle Gaskonstante der molekularen Komponente - Re_x= u x/ Reynolds-Zahl - Sc=/( D) Schmidt-Zahl - T absolute Temperatur - T_d=h _A ⁰ /R charakteristische Dissoziationstemperatur - u, v x- und y-Komponenten der Geschwindigkeit - U=u/u normierte x-Komponente der Geschwindigkeit - x, y Koordinaten parallel und senkrecht zur Platte Griechische Symbole - =A/ Dissoziationsgrad - Grenzschichtdicke - ₂ Impulsverlustdicke - Damköhler-Zahl der Oberflächenreaktion - =T/T normierte Temperatur - =y/ normierter Wandabstand - Wärmeleitfähigkeit - dynamische Viskosität - , ^* Ähnlichkeitskoordinaten - Dichte - Schubspannung Indizes A auf ein Atom bezogen - M auf ein Molekül bezogen - f auf den eingefrorenen Zustand bezogen - w auf die Wand bezogen - auf den Außenrand der Grenzschicht bezogen     

Representation of the rheological properties of polymer melts in terms of complex fluidity

In dynamic rheological experiments melt behavior is usually expressed in terms of complex viscosity ^* () or complex modulusG ^* (). In contrast, we attempted to use the complex fluidity ^* () = 1/µ ^* () to represent this behavior. The main interest is to simplify the complex-plane diagram and to simplify the determination of fundamental parameters such as the Newtonian viscosity or the parameter of relaxation-time distribution when a Cole-Cole type distribution can be applied. ^* () complex shear viscosity - () real part of the complex viscosity - () imaginary part of the complex viscosity - G ^* () complex shear modulus - G() storage modulus in shear - G() loss modulus in shear - J ^* () complex shear compliance - J() storage compliance in shear - J() loss compliance in shear - shear strain - rate of strain - angular frequency (rad/s) - shear stress - loss angle - ^* () complex shear fluidity - () real part of the complex fluidity - () imaginary part of the complex fluidity - ₀ zero-viscosity - ₀ average relaxation time - h parameter of relaxation-time distribution     

Streamwise pseudo-vortical structures and associated vorticity in the near-wall region of a wall-bounded turbulent shear flow

Streamwise pseudo-vortical motions near the wall in a fully-developed two-dimensional turbulent channel flow are clearly visualized in the plane perpendicular to the flow direction by a sophisticated hydrogen-bubble technique. This technique utilizes partially insulated fine wires, which generate hydrogen-bubble clusters at several distances from the wall. These flow visualizations also supply quantitative data on two instantaneous velocity components, and w, as well as the streamwise vorticity, _x . The vorticity field thus obtained shows quasi-periodicity in the spanwise direction and also a double-layer structure near the wall, both of which are qualitatively in good agreement with a pseudo-vortical motion model of the viscous wall-region.List of symbols C _i ,c _i ,d _i constants in Eqs. (2), (3) and (4) - H channel width (m) - Re _H Reynolds number (= U _c H/) - Re Reynolds number (= U _c /) - T period (s) - t time (s) - U mean streamwise velocity (m/s) - U _c center-line velocity (m/s) - u friction velocity (m/s) - u, , w velocity fluctuations (m/s) - x, y, z coordinates (m) - ^* displacement thickness (m) - momentum thickness (m) - mean low-speed streak spacing (m) - kinematic viscosity (m²/s) - phase difference - _x streamwise vorticity fluctuation (1/s) - ( )⁺ normalized by u and - (^–) root mean square value - (^–) statistical average This paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984     

Normal forms for random diffeomorphisms

Given a dynamical system (,, ,) and a random diffeomorphism (): ^{d d} with fixed point at x=0. The normal form problem is to construct a smooth near-identity nonlinear random coordinate transformation h() to make the random diffeomorphism ()=h()^–1() h() as simple as possible, preferably linear. The linearization D(, 0)=:A() generates a matrix cocycle for which the multiplicative ergodic theorem holds, providing us with stochastic analogues of eigenvalues (Lyapunov exponents) and eigenspaces. Now the development runs pretty much parallel to the deterministic one, the difference being that the appearance of turns all problems into infinite-dimensional ones. In particular, the range of the homological operator is in general not closed, making the conceptof-normal form necessary. The stochastic versions of resonance and averaging are developed. The case of simple Lyapunov spectrum is treated in detail.     

The dynamic and steady shear properties of synovial fluid and of the components making up synovial fluid

Steady-shear and dynamic properties of a pooled sample of cattle synovial fluid have been measured using techniques developed for low viscosity fluids. The rheological properties of synovial fluid were found to exhibit typical viscoelastic behaviour and can be described by the Carreau type A rheological model. Typical model parameters for the fluid are given; these may be useful for the analysis of the complex flow problems of joint lubrication.The two major constituents, hyaluronic acid and proteins, have been successfully separated from the pooled sample of synovial fluid. The rheological properties of the hyaluronic acid and the recombined hyaluronic acid-protein solutions of both equal and half the concentration of the constituents found in the original synovial fluid have been measured. These properties, when compared to those of the original synovial fluid, show an undeniable contribution of proteins to the flow behaviour of synovial fluid in joints. The effect of protein was found to be more prominent in hyaluronic acid of half the normal concentration found in synovial fluid, thus providing a possible explanation for the differences in flow behaviour observed between synovial fluid from certain diseased joints compared to normal joint fluid.Nomenclature A Ratio of angular amplitude of torsion head to oscillation input signal - G Storage modulus - G Loss modulus - I Moment of inertia of upper platen — torsion head assembly - K Restoring constant of torsion bar - N ₁ First normal-stress difference - R Platen radius - S ⁽ⁱ⁾ Geometric factor in the dynamic property analysis - t ₁ Characteristic time parameter of the Carreau model - X, Y Carreau model parameters - Z () Reimann Zeta function of - Carreau model parameter - Shear rate - Apparent steady-shear viscosity - ^* Complex dynamic viscosity - Dynamic viscosity - Imaginary part of the complex dynamic viscosity - ₀ Zero-shear viscosity - ₀ Cone angle - Carreau model characteristic time - Density of fluid - Shear stress - Phase difference between torsion head and oscillation input signals - ₀ Zero-shear rate first normal-stress coefficient - Oscillatory frequency     

Three-component,time-resolved velocity statistics in the wall region of a turbulent pipe flow

Three-component, coincident, time-resolved velocity measurements were obtained in the near wall region, y ⁺ < 100, of a fully developed turbulent pipe flow. The measurements were conducted in the ARL/PSU glycerin tunnel at a Reynolds number (Re _h), based on pipe radius and centerline velocity, of 6436 and an Re of approximately 730. The reported data include velocity statistics up to fourth order, Reynolds stresses and three component, coincident turbulent velocity spectral estimates. The current data are generally in quite good agreement with the fully developed channel flow direct numerical simulation (DNS) results of Antonia et al. (1992) at Re 700 - 700. The accuracy of the current experimental data and the very good agreement with the DNS results provides evidence for the accuracy of the DNS solutions and thus Antonia's conclusions of very near wall, y ⁺ < 20, Re dependence on turbulent velocity statistics. The very good agreement between the low Re rectangular channel flow DNS results and the low Re flat plate turbulent boundary layer statistics of Karlsson and Johansson (1988) suggests that for y ⁺ < 30 statistics of similar flows of differing geometry may be compared on the basis of equal Re . The current data are available on disk or by anonymous ftp by the first author.     

Some remarks on differential equations of quadratic type

In this paper we study differential equations of the formx(t) + x(t)=f(x(t)), x(0)=x ₀ C HereC is a closed, bounded convex subset of a Banach spaceX,f(C) C, and it is often assumed thatf(x) is a quadratic map. We study the differential equation by using the general theory of nonexpansive maps and nonexpansive, non-linear semigroups, and we obtain sharp results in a number of cases of interest. We give a formula for the Lipschitz constant off: C C, and we derive a precise explicit formula for the Lipschitz constant whenf is quadratic,C is the unit simplex inR ⁿ, and thel ₁ norm is used. We give a new proof of a theorem about nonexpansive semigroups; and we show that if the Lipschitz constant off: CC is less than or equal to one, then lim_tf(x(t))–x(t)=0 and, if {x(t):t 0} is precompact, then lim_tx(t) exists. Iff¦C=L¦C, whereL is a bounded linear operator, we apply the nonlinear theory to prove that (under mild further conditions on C) lim_t f(x(t))–x(t)=0 and that lim_t x(t) exists if {x(t):t 0} is precompact. However, forn 3 we give examples of quadratic mapsf of the unit simplex ofR ⁿ into itself such that lim_t x(t) fails to exist for mostx ₀ C andx(t) may be periodic. Our theorems answer several questions recently raised by J. Herod in connection with so-called model Boltzmann equations.     

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     

Optical co-ordinates in general relativity 2 — the problem of distance

Summary Let denote the congruence of null geodesics associated with a given optical observer inV ₄. We prove that determines a unique collection of vector fieldsM₍₎ ( =1, 2, 3) and ₍₀₎ overV ₄, satisfying a weak version of Killing's conditions.This allows a natural interpretation of these fields as the infinitesimal generators of spatial rotations and temporal translation relative to the given observer. We prove also that the definition of the fieldsM₍₎ and ₍₀₎ is mathematically equivalent to the choice of a distinguished affine parameter f along the curves of, playing the role of a retarded distance from the observer.The relation between f and other possible definitions of distance is discussed.

---

Sommario Sia la congruenza di geodetiche nulle associata ad un osservatore ottico assegnato nello spazio-tempoV ₄. Dimostriamo che determina un'unica collezione di campi vettorialiM₍₎ ( =1, 2, 3) e ₍₀₎ inV ₄ che soddisfano una versione in forma debole delle equazioni di Killing. Ciò suggerisce una naturale interpretazione di questi campi come generatori infinitesimi di rotazioni spaziali e traslazioni temporali relative all'osservatore assegnato. Dimostriamo anche che la definizione dei campiM₍₎, ₍₀₎ è matematicamente equivalente alla scelta di un parametro affine privilegiato f lungo le curve di, che gioca il ruolo di distanza ritardata dall'osservatore. Successivamente si esaminano i legami tra f ed altre possibili definizioni di distanza in grande.

---

---

Work performed in the sphere of activity of: Gruppo Nazionale per la Fisica Matematica del CNR.     

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.     

New integral estimates for deformations in terms of their nonlinear strains

If u is a bi-Lipschitzian deformation of a bounded Lipschitz domain in ⁿ (n2), we show that the L ^P norm (p1, pn) of a certain nonlinear strain function e(u) associated with u dominates the distance in L ^q (q= np/(n–p) if p if p>n) from u to a suitably chosen rigid motion of ⁿ . This work extends that of F. John, who proved corresponding estimates for p}>1 under the hypothesis that u has uniformly small strain. We also obtain a bound for the oscillation of Du in L ². These estimates are apparently the first to apply with no a priori pointwise hypotheses upon the strain of u. In ³ the integral e(u) ² d³ is dominated by typical hyperelastic energy functionals proposed in the literature for modeling the behavior of rubber; thus the case n=3, p=2 gives the first general bound for the deformations of such materials in terms of the associated nonlinear elastic work.     

Investigations of mean and turbulent flow characteristics of a two dimensional plane diffuser

This paper reports the investigation of mean and turbulent flow characteristics of a two-dimensional plane diffuser. Both experimental and theoretical details are considered. The experimental investigation consists of the measurement of mean velocity profiles, wall static pressure and turbulence stresses. Theoretical study involves the prediction of downstream velocity profiles and the distribution of turbulence kinetic energy using a well tested finite difference procedure. Two models, viz., Prandtl's mixing length hypothesis and k- model of turbulence, have been used and compared. The nondimensional static pressure distribution, the longitudinal pressure gradient, the pressure recovery coefficient, percentage recovery of static pressure, the variation of U _max/U _bar along the length of the diffuser and the blockage factor have been valuated from the predicted results and compared with the experimental data. Further, the predicted and the measured value of kinetic energy of turbulence have also been compared. It is seen that for the prediction of mean flow characteristics and to evaluate the performance of the diffuser, a simple turbulence model like Prandtl's mixing length hypothesis is quite adequate.List of symbols C ₁ , C ₂ ,C turbulence model constants - F _x body force - k kinetic energy of turbulence - l _m mixing length - L length of the diffuser - u, v, w rms value of the fluctuating velocity - u, v, w turbulent component of the velocity - mean velocity in the x direction - _A average velocity at inlet - U _bar average velocity in any cross section - U _max maximum velocity in any cross section - V mean velocity in the y direction - W local width of the diffuser at any cross section - x, y coordinates - dissipation rate of turbulence - _m eddy diffusivity - Von Karman constant - mixing length constant - _l laminar viscosity - _eff effective viscosity - v kinematic viscosity - density - _k effective Schmidt number for k - effective Schmidt number for - stream function - non dimensional stream function     

Relations between steady flow and oscillatory shear measurements

Summary Results are given of a comparison between dynamic oscillatory and steady shear flow measurements with some polymer melts. Comparison of the steady shear flow viscosity,, with the absolute value of the dynamic viscosity, ¦¦, at equal values of the shear rate,q, and the circular frequency,, has shown the relation thatCox andHerz had found empirically to be substantially correct.Further, the coefficients of the normal stress differences obtained by streaming birefringence techniques have been compared with 2G () · ^{– 2} in the same range of shear rates as covered by the viscosity measurements (G is the real part of the dynamic shear modulus). Two polystyrenes with narrow molecular weight distribution showed the same shift factor along the orq axis for the normal stress coefficients with respect to 2G () · ^{– 2} and the steady shear flow viscosities with respect to the real part of the dynamic viscosity,. For two polyethylenes the results are not so conclusive owing to the smallness of the shift factor found. An empirical equation is proposed predicting the main normal stress difference from dynamic measurements only.

---

Zusammenfassung Die Ergebnisse von Messungen unter erzwungenen Schwingungen und stationärer Scherströmung an einigen Polymerschmelzen werden miteinander verglichen. Der Vergleich der stationären Viskosität mit der absoluten dynamischen Viskosität ¦¦ bei gleichen Werten des Strömungsgradientenq und der Kreisfrequenz zeigt die Gültigkeit der empirischen Beziehung vonCox undHerz.Weiter wurden die Koeffizienten der Normalspannungsdifferenzen, welche durch Messung der Strömungsdoppelbrechung erhalten wurden, mit 2G() · ^–2 verglichen, und zwar wiederum bei gleichen Werten vonq und, wobeiG die Speicherkomponente des dynamischen Schubmoduls ist. Zwei Polystyrole mit enger Molekulargewichtsverteilung zeigen die gleiche Verschiebung entlang der-oderq-Achse für die Normalspannungskoeffizienten in bezug auf2G()· ^–2 und für die stationären Scherviskositäten in bezug auf den Realteil der dynamischen Viskosität. Für zwei Polyäthylene sind die Ergebnisse weniger signifikant, da die entsprechenden Verschiebungen zu klein waren. Eine empirische Beziehung zwischen den Hauptnormalspannungsdifferenzen und den dynamischen Meßwerten wird vorgeschlagen.

---

---

Paper presented at the British Society of Rheology Conference, held at Shrivenham, from 9th–12th September, 1968.     

Der Spannungs- und Verschiebungszustand drehsymmetrischer Membrane mit beliebig gekrümmter Meridiankurve

Zusammenfassung Die Einführung von Zylinderkoordinaten (x, r, ) in die Gleichgewichtsbedingungen der Schnittkräfte bzw. in die Beziehungen zwischen Verzerrung und Verschiebungen am differentialen Schalenabschnitt ermöglicht die Berechnung des Spannungs- und Verschiebungszustandes von drehsymmetrischen Membranen mit beliebig gekrümmter Meridiankurve auf die Integration einer einfachen, linearen partiellen Differentialgleichung zweiter Ordnung für eine charakteristische FunktionF bzw. zurückzuführen. Eine geschlossene Lösung und damit eine Darstellung der Schnittkräfte und Verschiebungen durch explizite Formeln ist bei harmonischer Belastung cosn für zwei Funktionsgruppen=x ² und=x ^–3 möglich. Im Sonderfall der drehsymmetrischen und der antimetrischen Belastung mitn=0 undn=1 gelten die Gleichungen der Schnitt- und Verschiebungsgrößen für eine beliebige Meridianfunktion=(). Die Betrachtungen der Randbedingungen offener Schalen bei harmonischer Belastung geben über die infinitesimalen Deformationen einer drehsymmetrischen Membran mit überall negativer Krümmung Aufschluß.     

A method for the analysis of recirculating turbulent flows in a cavity

The complex fluid-dynamic aspects of a turbulent recirculating flow in a cavity with axial throughflow, and a rotating wall, were investigated by adopting a simple procedure for evaluating the turbulent stresses. The flow field was divided into two regions, a core and a wall region respectively. A wall function was adopted in the zones near to the solid boundaries, while a constant eddy diffusivity was assumed, in the core, following the indications of computed heat transfer coefficients in comparison with existing experimental data. The distributions of the stream function and of the tangential velocity are presented for a range of the rotational Reynolds number of the rotating wall and of the Reynolds number of the throughflow.

---

Turbulente Rezirkulationsströmung in einem Hohlraum

Zusammenfassung Die komplizierten fluiddynamischen Aspekte einer turbulenten Rezirkulationsströmung in einem Hohlraum mit axialem Durchfluß und einer rotierenden Wand werden unter Verwendung einer vereinfachten Methode zur Berechnung der turbulenten Spannungen betrachtet. Das Strömungsfeld wird in einen Kern und einen Wandbereich aufgeteilt. Für die wandnahen Zonen wird eine Wandfunktion angenommen, während im Kern mit konstanter Wirbeldiffusivität gerechnet wird, was durch den Vergleich berechneter mit gemessenen Wärmeübergangskoeffizienten gerechtfertigt erscheint. Verteilungen der Stromfunktion und der tangentialen Geschwindigkeit sind für einen bestimmten Bereich der Reynoldszahlen für die Wandrotation und der für den Durchfluß angegeben.

---

Nomenclature L axial length of enclosure - P dimensionless pressure, p^*2 - p static pressure - R dimensionless radial coordinate, r/r^* - r radial coordinate - r^* reference length, equal to r_O for enclosure - r_i radii of inlet and exit apertures - Re Reynolds number, v^*r^*/ - Re_i pipe Reynolds number, ¯v_zi(2r_i)/ - Re_t turbulent Reynolds number, Re(/) - Re rotational Reynolds number, r ₀ ² / - t dimensionless time,t/(r^*/v^*) - t time - V_r, V, Vz dimensionless velocity components, V_r/v^*, v, v_z/v^* - v_i turbulent fluctuation of the i-component of velocity - v_r, v, v_z velocity components - v^* reference velocity, equal to ¯v_zi for enclosure - X coordinate along a wall, x/r^* - Y coordinate normal to a wall, y/r^* - Z dimensionless axial coordinate, z/r^* - z axial coordinate - eddy diffusivity for momentum - dynamic viscosity - kinematic viscosity - density - shear stress - dimensionless shear stress, /v^*2 - dimensionless stream function, /r*²v*² - stream function - angular velocity - tangential vorticity component - ()_eff effective - ()_l laminar - ()_t turbulent - mean over the time     

Viscometric properties of viscous theta solutions of monodisperse polystyrene

Theta solvents for polystyrene are prepared from high-viscosity blends of styrene and low-molecular-weight polystyrene, and then used to make dilute solutions with monodisperse polystyrene solutes of high-M = 2.3, 6.0, 9.0, 18.0 · 10⁵. A Weissenberg rheogoniometer is used to measure the non-Newtonian viscosity as a function of shear stress, for low values, and also the complex viscosity components and as functions of frequency. A capillary viscometer is used for high- measurements of(). Viscometric properties, at room temperature, are analyzed as functions of high-molecular-weight solute concentrationc with parameters of constant or to obtain [()], [ ()], and [ ()]. Such a collection of data has apparently not previously been available for polymers in theta solvents (in which Gaussian chain statistics prevail). Also unique is the achievement of high stress ( = 2 10⁴ Pa) at low shear rate, by virtue of high solvent viscosity which is not characteristic of other known theta solvents.