The effect of background particulates on the delayed transition of a heated 9:1 ellipsoid

An experimental study was done to quantify the effects of a variety of background particulates on the delayed laminar-turbulent transition of a thermally stabilized boundary layer in water. A Laser-Doppler Velocimeter system was used to measure the location of boundary layer transition on a 50 mm diameter, 9:1 fineness ratio ellipsoid. The ellipsoid had a 0.15 m RMS surface finish. Boundary layer transition locations were determined for length Reynolds numbers ranging from 3.0 × 10⁶ to 7.5 × 106. The ellipsoid was tested in three different heating conditions in water seeded with particles of four distinct size ranges. For each level of boundary layer heating, measurements of transition were made for clean water and subsequently, water seeded with 12.5 m, 38.9 m, 85.5 m and 123.2 m particles, alternately. The three surface heating conditions tested were no heating, T = 10°C and T = 15°C where T is the difference between the inlet model heating water temperature, T _i, and free stream water temperature, T . The effects of particle concentration were studied for 85.5 m and 123.2 m particulates.The results of the study can be summarized as follows. The 12.5 m and 38.9 m particles has no measurable effect on transition for any of the test conditions. However, transition was significantly affected by the 85.5 m and 123.2 m particles. Above a length Reynolds number of 4 × 10⁶ the boundary layer transition location moved forward on the body due to the effect of the 85.5 m particles for all heating conditions. The largest percentage changes in transition location from clean water, were observed for 85.5 m particles seeded water.Transition measurements made with varied concentrations of background particulates indicated that the effect of the 85.5 m particles on the transition of the model reached a plateau between 2.65 particulates/ml concentration and 4.2 particles/ml. Measurements made with 123.3 m particles at concentrations up to 0.3 part/ml indicated no similar plateau.     

A fourfold generalization of Peaucellier's inversion cell

New inversors are proposed which are generalizations of the well-known inversor of Peaucellier. It appears that a kite in the Peaucellier cell is replaceable by an arbitrary 4-bar linkage (abhk) whereas the direction and length of the straight line, produced by the inversor, can be manipulated through the particular choice of the relative polar coordinates of a vertex A of the triangular input link. Formulas are derived for practical inversors with a revolving input link. The ones selected are basically governed by the choice of two transmission angles, ₁ and ₃, by the length L of the acquired line, as well as by its direction represented by the angle /2- comprised between the line L and the frame.

---

Sommario Vengono proposti nuovi inversori quali generalizzazioni del ben noto inversore di Peaucellier. Si mostra che un quadrilatero isoscele nella cella di Peaucellier è sostituibile da un'arbitraria connessione a quattro barre, mentre la direzione e la lunghezza della retta, prodotta dall'inversore, possono essere manipolate con la particolare scelta delle relative coordinate polari di un vertice A del collegamento triangolare di input. Vengono derivate formule per inversori funzionali con un collegamento di input rotante. Quelli solezionati sono principalmente governati dalla scelta di due angoli di trasmissione, ₁ e ₃, dalla lunghezza L della linea ottenuta, così come dalla sua direzione rappresentata dall' angolo /2- compreso tra la linea L ed il riferimento.

     

Turbulent dispersion of droplets for air flow in a pipe

An optical technique was used to study the dispersion of 50 m, 90 m and 150 m droplets downstream of a source located in the center of a vertical pipe through which turbulent air is flowing. A turbulent dispersion coefficient, _P, and a mean-square of the fluctuations in the turbulent velocity, v _p ² , are determined from the change of the measured mean-square displacement of the droplets over the pipe cross section with time. The interesting aspect of the experiments is that they explored conditions where the inertia of the particles is believed to be a much more important effect than that of the crossing of trajectories associated with the inequality of the average velocities of the particles and the fluid.     

On the viscosity of suspensions of solid spheres

A cell theory is used to derive the dependence of the zero-shear-rate viscosity on volume concentration for a suspension of uniform, solid, neutrally buoyant spheres. This result reduces to Einstein's solution at infinite dilution and to Frankel and Acrivos's expression in the limit as the concentration approaches its maximum value. Good agreement is found between the solution and the available data for the entire concentration range, provided that the maximum concentration is determined from the viscosity data themselves.Nomenclature a radius of sphere - d the distance separating the sphere surfaces measured parallel to the line connecting the sphere centers - E energy dissipation rate in one-half the liquid volume separating the spheres - E _cell total energy dissipation rate in the cell - E _homogeneous energy dissipation rate in the cell of a hypothetical one phase fluid - E _interaction energy dissipation rate in the cell due to sphere interactions - E _sphere energy dissipation rate in the cell due to the sphere at the cell center - F force on one sphere - h minimum separation distance between two spheres - J (1/2)d = one-half the distance separating the sphere surfaces measured parallel to the line connecting the sphere centers - p pressure - W velocity of one sphere in squeezing flow between two spheres relative to the midpoint of the line connecting the sphere centers - _i unit vectors in thei-th direction - elongation rate - viscosity of the suspending fluid - _r */ = relative viscosity - * viscosity of the suspension - the total stress tensor - the part of the total stress tensor that vanishes at equilibrium - volume fraction of spheres     

The Reissner–Sagoci Problem for a Non-homogeneous Half-space with a Surface Constraint

This paper deals with the problem of twisting a non-homogeneous, isotropic, half-space by rotating a circular part of its boundary surface (0 r < a, z = 0) through a given angle. A ring (a < r < b, z = 0) outside the circle is stress-free and the remaining part (r > b, z = 0) is rigidly clamped. The shear modulus is assumed to vary with the cylindrical coordinates, r, z by the relation (z) = ₁(c + z), c 0 where ₁, c and are real constants. Expressions for some quantities of physical importance, such as torque applied at the surface of the disk and stress intensity factors, are obtained. The effects of non-homogeneity on torque and stress intensity factor are illustrated graphically.     

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.     

States of classical statistical mechanical systems of infinitely many particles. I

We study a general mathematical model of a classical system of infinitely many point particles. The space X of infinite particle configurations is equipped with a natural topology as well as a measurable structure related to it. It is also connected with a family {X _A } of local spaces of finite configurations indexed by bounded open sets A in the one-particle space E. A theorem analogous to Kolmogoroff's fundamental theorem for stochastic processes is proved, according to which a consistent family { _A } of local probability measures _A defined on the X _A gives rise to a unique probability measure on X. We also study the problem of integral representation for positive linear forms defined over some linear space of real functions on X. We prove that a positive linear form F(f), defined for functions f in the class C+P, admits a uniquely determined integral representation F(f)= f () d, where is a probability measure over X.     

Asymmetric snap-buckling of a column restrained by a stiff wire

Summary A relatively simple example of asymmetric snap-through buckling in a continuous structure is the nonlinear problem of a cantilevered column restrained at its tip by a stiff wire, which is inclined at an acute angle to the column centerline, and loaded at its tip by a force perpendicular to the centerline. A parameter called , which is the nondimensional ratio of the flexural rigidity of the column to the combined extensional stiffness of the wire and the column centerline, determines the essential features of the buckling. If is zero, or is small compared to unity, the bending of the column is small enough to justify the use of linear bending theory for the column. Hence, even though the constraint is nonlinear, the solution to this problem is obtained in closed form. The critical point for the structure is found to be an asymmetric branching point for =0, while for positive, the critical point is a snap-through type. The effect of is similar to that induced by initial imperfections in more complex structures. For very small , the critical load is markedly decreased from the value for =0. Moreover, the graph of the load vs. tip deflection has the appearance of having an acute discontinuity in slope at the critical point for very small, although it is actually found that the graph has a horizontal tangent there.

---

Sommario Un esempio relativamente semplice di collasso asimmetrico per carico di punta in una struttura continua è il problema non lineare di una colonna incastrata ad una estremità e vincolata all'altra estremità da un filo rigido che forma un angolo acuto con l'asse della colonna e caricata a quella estremità con una forza perpendicolare allo stesso asse. Un parametro chiamato , che è il rapporto adimensionale fra la rigidità flessionale della colonna e la rigidezza longitudinale del filo e dell'asse della colonna, determina le caratteristiche essenziali del cedimento. Se è zero o è piccolo rispetto all'unità, l'inflessione della colonna è sufficientemente piccola per giustificare l'uso della teoria lineare di inflessione per la colonna. Di conseguenza anche se il vincolo non è lineare la soluzione del problema è ottenuta in forma compatta. Il punto critico della struttura si trova nel punto di biforcazione asimmetrica per =0, mentre per positivo, il punto critico rappresenta un punto di collasso. L'effetto di è simile a quello prodotto da imperfezioni iniziali in strutture più complesse. Per molto piccolo il carico critico è notevolmente ridotto rispetto al valore per =0. Inoltre il grafico del carico in funzione della curvatura all'estremità sembra avere una netta discontinuità nella pendenza dal punto di biforcazione per molto piccolo benchè, in realtà si trovi che il grafico ha lì una tangente orizzontale.

---

This research was supported by the Advanced Research Projects Agency (Project DEFENDER) and was monitored by the U.S. Army Research Office, Durham, under Contract DA-31-124-ARO-D-257.     

Singular solutions of some quasilinear elliptic equations

We study isolated singularities of the quasilinear equation in an open set of ^N , where 1 < p N, p -1 q < N(p — 1)/ (N -p). We prove that, for any positive solution, if a singularity at the origin is not removable then either or u(x)/(x) any positive constant as x 0 where is the fundamental solution of the p-harmonic equation: . Global positive solutions are also classified.     

On natural convection from a vertical plate with a prescribed surface heat flux in porous media

This paper presents a theoretical and numerical investigation of the natural convection boundary-layer along a vertical surface, which is embedded in a porous medium, when the surface heat flux varies as (1 +x ²)), where is a constant andx is the distance along the surface. It is shown that for > -1/2 the solution develops from a similarity solution which is valid for small values ofx to one which is valid for large values ofx. However, when -1/2 no similarity solutions exist for large values ofx and it is found that there are two cases to consider, namely < -1/2 and = -1/2. The wall temperature and the velocity at large distances along the plate are determined for a range of values of .Notation g Gravitational acceleration - k Thermal conductivity of the saturated porous medium - K Permeability of the porous medium - l Typical streamwise length - q _w Uniform heat flux on the wall - Ra Rayleigh number, =gK(q _w /k)l/(v) - T Temperature - Too Temperature far from the plate - u, v Components of seepage velocity in the x and y directions - x, y Cartesian coordinates - Thermal diffusivity of the fluid saturated porous medium - The coefficient of thermal expansion - An undetermined constant - Porosity of the porous medium - Similarity variable, =y(1+x ) ^/3/x ^1/3 - A preassigned constant - Kinematic viscosity - Nondimensional temperature, =(T – T )Ra^1/3 k/qw - Similarity variable, = =y(log_e x)^1/3/x ^2/3 - Similarity variable, =y/x ^2/3 - Stream function     

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     

Zur Berechnung der Filmverdunstung von Kohlenwasserstoffen in einem Heißluftstrom

The results of an analytical approximation method to predict the film vaporization are compared with the predictions of a finite difference method of Hermitian type. The analytically estimated rate of vaporization of different hydrocarbons, which is the most important value for practical applications, deviates only a few percents from the numerically estimated value.

---

Zur Berechnung der Filmverdunstung von Kohlenwasserstoffen in einem Heißluftstrom

Zusammenfassung Es wird ein Näherungsverfahren zur Berechnung der Filmverdunstung dargestellt, bei dem eine vollständige Lösung der miteinander gekoppelten Grenzschichtgleichungen entfallt. Die nach dieser analytischen Methode ermittelte Verdunstung verschiedener Kohlenwasserstoffe wird mit Werten verglichen, die nach einem Differenzenverfahren vom Hermiteschen Typ berechnet wurden. Es zeigt sich, daß die analytisch berechnete Verdunstungsrate, die für praktische Anwendungen wichtigste Größe, nur wenige Prozent von dem numerisch ermittelten Wert abweicht.

---

Formelzeichen c _ew Konzentrationsdifferenz c_1e -c _1w - c _i Massenkonzentration der Komponentei - c_p, c_pi spezifische Wärmekapazität bei konstantem Druck des Gemisches — der Komponentei - D ₁₂ binärer Diffusionskoeffizient - f dimensionslose Stromfunktion - f dimensionslose Geschwindigkeit - g () allgemeine Funktion - m ₁ Massenstromdichte der Komponente 1 - m ^* dimensionslose Massenstromdichte, G1. (4.8) - M, M_i Molgewicht, — der Komponentei - P, P _i Druck, Partialdruck der Komponentei - Pr Prandtlzahl,C _p/ - q Wärmestromdichte - r ₁ Verdampfungswärme - R allgemeine Gaskonstante - Sc Schmidtzahl/D ₁₂ - T absolute Temperatur - u Geschwindigkeitskomponente inx-Richtung - v Geschwindigkeitskomponente iny-Richtung - x Längskoordinate - y Querkoordinate - z dimensionslose Konzentration - dimensionslose Funktion/ _{e e} - transformierte Koordinatey - dimensionslose Temperatur (T-T _w)/(T_e-T_w) - Wärmeleitfähigkeit des Gemisches - Zähigkeit des Gemisches - transformierte Koordinate - Dichte des Gemisches - Stromfunktion Indizes e am Außenrand der Grenzschicht - i Stoffi - w an der Filmoberfläche - 1, 2 Komponente 1, 2 - () Ableitung ()/ _n      

Analysis of suddenly started laminar flow in the entrance region of a circular tube

Suddenly started laminar flow in the entrance region of a circular tube, with constant inlet velocity, is investigated analytically by using integral momentum approach. A closed form solution to the integral momentum equation is obtained by the method of characteristics to determine boundary layer thickness, entrance length, velocity profile, and pressure gradient.Nomenclature M(, , ) a function - N(, , ) a function - p pressure - p* p/1/2U ², dimensionless pressure - Q(, , ) a function - R radius of the tube - r radial distance - Re 2RU/, Reynolds number - t time - U inlet velocity, constant for all time, uniform over the cross section - u velocity in the boundary layer - u* u/U, dimensionless velocity - u ₁ velocity in the inviscid core - x axial distance - y distance perpendicular to the axis of the tube - y* y/R, dimensionless distance perpendicular to the axis - boundary layer thickness - * displacement thickness - /R, dimensionless boundary layer thickness - momentum thickness - absolute viscosity of the fluid - /, kinematic viscosity of the fluid - x/(R Re), dimensionless axial distance - density of the fluid - tU/(R Re), dimensionless time - _w wall shear stress     

Die Temperaturabhängigkeit der Oberflächenspannung von reinen Kältemitteln vom Tripelpunkt bis zum kritischen Punkt

Zusammenfassung Die Oberflächenspannung von sechs reinen Substanzen — SF₆, CCl₃F, CCl₂F₂, CClF₃, CBrF₃ und CHClF₂ — wurde mit Hilfe einer modifizierten Kapillarmethode gemessen. Die zur Berechnung der Oberflächenspannung erforderlichen Sättigungsdichten und wurden teils aus vorhandenen Zustandsgleichungen, teils aus ebenfalls gemessenen Brechungsindizes bestimmt. Die Temperaturabhängigkeit der Oberflächenspannung läßt sich durch einen erweiterten Ansatz nach van der Waals =_O (T_c-T)(1+...) darstellen, wobei bei einfachen Stoffen ein eingliedriger, bei polaren und assoziierenden Stoffen ein zweigliedriger Ansatz notwendig und ausreichend ist. Für den kritischen Exponenten der Oberflächenspannung wurde ein von der molekularen Substanz weitgehend unabhängiger Wert von =1.284±0.005 gefunden.

---

Temperature dependence of surface tension of pure refrigerants from triple point up to the critical point

The surface tension of six fluids (SF₆, CCl₃F, CCl₂F₂, CClF₃, CBrF₃, CHClF₂) have been measured by means of a modified capillary rise method. The liquid vapor densities, which are needed to calculate the surface tension, have partly been determined by means of refractive indices simultaneously measured in the same apparatus. The temperature dependence of the surface tension is described by an extended van der Waals power law =_O(T_c-T)(1+...). For simple fluids one term and for polar and associating fluids two terms are necessary and sufficient. The critical exponent is found to be 1.284 ± 0.005 and nearly independent of the molecular structure.

---

Formelzeichen a² Laplace-Koeffizient - a Parameter - B_O, B_on Koeffizient der Koexistenzkurve - g Erdbeschleunigung - H Höhe, kapillare Steighöhe - LL Lorentz-Lorenz-Funktion oder Refraktionskonstante - M molare Masse - M Zahl der Meßwerte - N Zahl der unbekannten Parameter - n Brechungsindex - p Druck - R,r Radius - s Entropie - SD Standardabweichung - T, t Temperatur - u innere Energie Griechische Formelzeichen Exponent des Laplace-Koeffizienten - Exponent der Koexistenzkurve - 2. Exponent der Oberflächenspannung - Wellenlänge des Lichts - Exponent der Oberflächenspannung - _D Dipolmoment - , Dichte der Flüssigkeit bzw. des Dampfes - Oberflächenspannung - reduzierte Temperatur (1-T/T_c) - ² gewichtete Varianz Indizes c kritischer Zustand - D Differenz - m Mittelwert - T Isotherme - t Zustand am Tripelpunkt - S Zustand am Schmelzpunkt - bezogen auf Oberfläche     

Elastohydrodynamic lubrication in a plane slider bearing

Summary The present work deals with the case of a two-dimensional slider bearing with a rigid pad and an elastic bearing. Fluid viscosity is assumed to be only a pressure function. We determined the bearing deformation, the pressure distribution and the load capacity at different values of the inclination angle of the slider, with a numerical integration of the system consisting of the elasticity and Reynolds equations. The results show that, with an iso-viscous fluid, bearing elasticity causes a load capacity decrease. Instead bearing elasticity together with the variation of fluid viscosity due to pressure causes a load capacity greater than that of the iso-viscous case (=0).

---

Sommario Il presente lavoro studia il problema della coppia prismatica lubrificata con pattino rigido di allungamento infinito e cuscinetto deformabile; si suppone che la viscosità del fluido sia funzione della sola pressione. Il sistema di equazioni, costituito dall'equazione di Reynolds e dall'equazione dell'elasticità, è stato risolto numericamente, determinando la deformazione del cuscinetto, andamento della pressione e la capacità di carico per diversi valori dell'inclinazione del pattino. I risultati dimostrano che, con fluido isoviscoso, la deformabilità del cuscinetto determina una riduzione della capacità di carico. Se si considera, invece, effetto combinato dell'elasticità del cuscinetto e della variazione della viscosità del fluido, la capacità di carico risulta maggiore di quella che si ottiene con fluido isoviscoso (=0).

---

Nomenclature /L - /L - x/L - x/L - - ¯C CZ/h ₁ - E elasticity modulus - h film thickness - H elastic deformation of the bearing - h ₁ minimum film thickness - h ₂ inlet thickness - inclination of the pad - h Z/h ₁ - HZ/h ₁ - L pad length - viscosity - ₀ viscosity with no over-pressure - p over pressure - p P _ec-P _rc where:ec=elastic caserc=rigid case - P h ₁ ² /6₀VL - h ₂/h ₁=1+L/h ₁ - FV bearing velocity - W load capacity per unit width - Wh ₂ ¹ /6₀ VL ² - Z E h ₃ ¹ /12 ₀ VL ² A first version of this paper was presented at the 7th National AIMETA congress, held at Trieste, October 2–5, 1984. This work was supported by C.N.R.     

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     

Multiple scattering of anisotropic radiation through a layer of clouds

A method by which the transport of anisotropic radiative multiple scattering can be predicted is developed in this paper. A one-dimensional integral intensity model and a three-dimensional diffusion intensity model are both constructed. The former provides a closed-form solution, while the latter involves successive approximation and Gauss's quadrature. On the basis of these methods, the reflection and transmission of solar radiation in a homogeneous cloud layer are computed. The results differ from those for isotropic and Rayleigh scattering assumptions and illustrate the effects on transmission and reflectivity of optical thickness, wavelength, incidence angle, and albedo of single scattering.Nomenclature D ⁺ transmitted diffusion radiation intensity [W/cm² sr · m] - D ^– reflected diffusion radiation intensity [W/cm² sr · m] - I pencil of radiation or specific intensity [W/cm² sr · m] - I ₀ solar irradiance [W/cm²] - K extinction cross-section or total cross-section, + - s (u, ), unit scattered radiation vector - s ₀ (u ₀, ₀), unit incident radiation vector - t optical thickness - u cosine of the viewing angle, , which is measured from the vertical - u ₀ cosine of the angle of incident, ₀, which is measured from the vertical - absorption cross-section - scattering function - absorption coefficient - scattering angle, s · s ₀ - scattering cross-section - scattered azimuthal angle - ₀ incident azimuthal angle - a sphere - a solid angle - ₀ albedo of single scattering,      

Non-persistence of Homoclinic Connections for Perturbed Integrable Reversible Systems

The dynamics of an analytic reversible vector field (X,) is studied in with one real parameter close to 0; X=0 is a fixed point. The differential D_x (0,0) generates an oscillatory dynamics with a frequency of order 1—due to two simple, opposite eigenvalues lying on the imaginary axis—and it also generates a slow dynamics which changes from a hyperbolic type—eigenvalues are —to an elliptic type—eigenvalues are —as passes trough 0. The existence of reversible homoclinic connections to periodic orbits is known for such vector fields. In this paper we study a particular subclass of such vector fields, obtained by small reversible perturbations of the normal form. We give an explicit condition on the perturbation, generically satisfied, which prevents the existence of a homoclinic connections to 0 for the perturbed system. The normal form system of any order admits a reversible homoclinic connection to 0, which then does not survive under perturbation of higher order. It will be seen that normal form essentially decouples the hyperbolic and elliptic part of the linearization to any chosen algebraic order. However, this decoupling does not persist arbitrary reversible perturbation, which finally causes the appearance of small amplitude oscillations.     

Remarks on the bifurcation of solutions of a nonlinear eigenvalue problem

Conclusions We have investigated solutions of equation (3) when ² is an eigenvalue of the linearized operator (13) and when it is not. In Section 4 we have shown that for 0 and ² = _i ² we have exactly two nontrivial solutions which bifurcate to the right of _i ² ; these solutions are shown to exist in an interval ( _i ² , _i ² + ₀). The method of Section 3 may then be used to extend these two solutions to the right of _i ² + ₀ providing that ²= _i ² + ₀ is not an eigenvalue of the linear operator (13) evaluated at = ± ₁. Either a solution can be uniquely extended, or there exists a value of ²where the bifurcation method must be applied again³.While the method used here gives the exact number of solutions bifurcating from _i ² , other problems remain open; for example, it is still not proven that the two bifurcating branches have i zeros, as is the case for Hammerstein operators with oscillation kernels [4]. The conjecture of Odeh and Tadjbakhsh that there are exactly 2(i+1) nontrivial solutions in the interval _i ² < _i ^+1/2 remains un-answered, although it would be proven if one could show that there is no secondary bifurcation as in the cases of Kolodner [7] and Coffman [8].     

Spin-up in a saturated porous medium

It is shown that, on the Brinkman model, spin-up is confined to boundary layers whose thickness is of order k ^1/2, and the spin-up is established in a time of order k/, where k, , and denote permeability, density, porosity and dynamic viscosity, respectively.