Some self-similar solutions of the Leibenzon equation

Self-similar one-dimensional solutions of the Leibenzon equation c²_t= _zz ^k (z 0, k 2) are considered. Approximate solutions are constructed for the two cases in which the initial value = ₁ = const > 0 and on the boundary either a constant value = ₂ < ₁ is maintained or the flow (directed outwards) is given. In the first problem the dependence of the boundary flow on the governing parameters is determined. A characteristic property of the types of motion in question is the existence near the boundary of a region, expanding with time, in which the flow is almost independent of the coordinate.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 145–150, September–October, 1991.     

The measurement of flow-induced surface displacement on a compliant surface by optical holographic interferometry

Flow-induced surface displacements that form on a singlelayer passive isotropic viscoelastic compliant surface as a result of the interaction with a turbulent boundary layer are measured by non-intrusive optical holographic interferometry in connection with an interactive fringe-processing system. The purpose for developing this method is to obtain the instantaneous topographic features of a whole field of the displacements of the compliant surface. Information about dimensions of the foot prints of the turbulence on the compliant surface are obtained in the form of line contours and isometric phase maps. These experimental data are essential in order to determine statistical measures of the random topography of the compliant surface. Furthermore, by coupling with the simultaneous measurements of the turbulence field, the physics of the alteration of the turbulent boundary layer by the undulating surface can then be better understood.List of symbols A light amplitude - A ₀ light amplitude of object beam - A _r light amplitude of reference beam - _c light speed - d _e diameter of laser beam - d _p diameter of pinhole - d _s spatial frequency - E exposure energy - f focal length - I _in light intensity at photographic plate - I ₀ light intensity of object beam - I _r light intensity of reference beam - k wave number, 2 / - L _in light amplitude records on the emulsion surface - L ₀ light ray of object beam - L _r light ray of reference beam - R Reynolds number based on momentum thickness, V / - t time - t _e exposure time - T ₀ uniform background light transmittance - T _re reconstructed light through hologram - T _t amplitude transmittance of hologram - V freestream flow speed - radian frequency, k c - x distance from the leading edge of the flat plate - z position - slope - optical path difference - phase - ₀ phase of object beam - _r phase of reference beam - wavelength - momentum thickness - ₀ angle between the object beam with respect to the normal of the photographic plate - _r angle between the reference beam with respect to the normal of the photographic plate - kinematic viscosity of water     

Hydrodynamic lubrication of ball and socket joints

Summary As a step towards formation of design rules for ball and socket joints, a solution of the lubrication equations for an ideal joint (in which the socket completely surrounds the ball) under steady loading and rotation is presented. The results are compared with those for an infinitely wide plain bearing. The friction for the ball is found to be higher, and the load capacity lower, than that of the plain bearing under similar conditions.Notation h lubricant film thickness - p pressure - absolute viscosity - average of v, w, over film thickness h - v, w components of fluid velocity in , directions - U , U components of relative velocity of surfaces in , directions - R mean radius of spheres - H difference between ball and socket radii - c eccentricity ratio - angular speed of ball - B =pH ²/6cR ² - P load component along Oy - N load component along Oz - M driving moment to spin around Ox - moment coefficient - _p , P _p suffix p refers to unit width of plain bearing     

Rheology of carbon black suspensions. III. Sol-gel transition system

Linear viscoelastic properties of carbon black (CB) suspensions with various CB volume fractions () in a rosin-modified phenol resin type varnish (Varnish-1) were investigated at various temperatures (T). The CB/Varnish-1 suspensions exhibited a sol-gel transition on an increase in , and the _gel value at the gelation point decreased with increasing T. This T dependence of _gel, being opposite to the dependence seen for usual gelling systems, can be related to a phenol resin type polymeric component included in the Varnish-1. At low T, this polymeric component appeared to be rather well solvated in the Varnish-1 thereby allowing the gelation due to bare attraction between the CB particles at large . In contrast, at high T, the polymeric component appeared to have been less solvated, as evidence from a moderate failure of the time-temperature superposition of pure Varnish-1 and a decrease of its elasticity (in a shifted frequency scale) with increasing T. This less solvated polymeric component would have been adsorbed on the CB particles, thereby allowing the agglomeration of the particles at small _gel at high T.     

Viscous properties of calcium carbonate filled polymer melts

Summary The viscous properties of calcium carbonate filled polyethylene and polystyrene melts were examined. The relative vircosity _r defined in the previous paper gave an asymtptotic value( _r)_l in the range of the shear stress below 10⁵ dyne/cm².( _r)_l of the calcium carbonate filled system was higher than that of the glass beads or glass balloons filled system at the same volume fraction of the filler. Maron-Pierce equation with ₀ = 0.44 was able to approximate the( _r)_l — relationship. However, it was deduced here that the high value of( _r)_l of calcium carbonyl filled system was due to the apparent increase of and this increase was attributed to the fixed polymer layer formed on the powder particle. By assuming the particle as a sphere with a diameter of 2 µm, the thickness of the fixed polymer layer was estimated as about 0.17 µm. The yield stress estimated from the Casson's plots increased exponentially with.

---

Zusammenfassung Es wurden die viskosen Eigenschaften von Polyäthylen-und Polystyrol-Schmelzen untersucht, die mit Kalziumkarbonat-Teilchen gefüllt waren. Für die relative Viskosität _r, wie sie in einer vorangegangenen Veröffentlichung definiert worden war, ergab sich bei Schubspannungen unterhalb 10⁵ dyn/cm² ein asymptotischer Wert( _r)_l. Dieser war bei den mit Kalziumkarbonat gefüllten Schmelzen höher als bei Schmelzen, die bis zur gleichen Volumenkonzentration mit Glaskugeln oder Glasballons gefüllt waren. Die ( _r) _l —-Abhängigkeit ließ sich durch eine Gleichung nachMaron und Pierce mit ₀ = 0,44 beschreiben. Es wurde jedoch geschlossen, daß der hohe( _r)_l-Wert der mit Kalziumkarbonat gefüllten Schmelzen auf eine scheinbare Zunahme von zurückzuführen ist, verursacht durch eine feste Polymerschicht auf der Teilchenoberfläche. Unter Annahme kugelförmiger Teilchen mit einem Durchmesser von 2 µm ließ sich die zugeordnete Schichtdicke zu 0,17 µm abschätzen. Die mittels der Casson-Beziehung geschätzte Fließspannung ergab eine exponentielle-Abhängigkeit.

---

---

With 7 figures and 1 table     

A new interpretation of viscosity and yield stress in dense slurries: Coal and other irregular particles

Illinois coal was ground and wet-sieved to prepare three powder stocks whose particle-size distributions were characterized. Three suspending fluids were used (glycerin, bromonaphthalene, Aroclor), with viscosities _s that differed by a factor of 100 and with very different chemistries, but whose densities matched that of the coal. Suspensions were prepared under vacuum, with coal volume fractions that ranged up to 0.46. Viscosities were measured in a cone-and-plate over a shear rate range 10^–3–10² s^–1. Reduced viscosity _r = /_s is correlated in the high-shear limit ( ) with/ _M, where _M is the maximum packing fraction for the high-shear microstructure, to reveal the roles of size distribution and suspending fluid character. A new model that invokes the stress-dependence of _M is found to correlate _r well under non-Newtonian conditions with simultaneous prediction of yield stress at sufficiently high; a critical result is that stress and not governs the microstructure and rheology. Numerous experimental anomalies provide insight into suspension behavior.     

Time-resolved, 3D, laser-induced fluorescence measurements of fine-structure passive scalar mixing in a tubular reactor

A three-dimensional, time-resolved, laser-induced fluorescence (3D-LIF) technique was developed to measure the turbulent (liquid-liquid) mixing of a conserved passive scalar in the wake of an injector inserted perpendicularly into a tubular reactor with Re=4,000. In this technique, a horizontal laser sheet was traversed in its normal direction through the measurement section. Three-dimensional scalar fields were reconstructed from the 2D images captured at consecutive, closely spaced levels by means of a high-speed CCD camera. The ultimate goal of the measurements was to assess the downstream development of the 3D scalar fields (in terms of the full scalar gradient vector field and its associated scalar energy dissipation rate) in an industrial flow with significant advection velocity. As a result of this advection velocity, the measured 3D scalar field is artificially skewed during a scan period. A method to correct for this skewing was developed, tested and applied. Analysis of the results show consistent physical behaviour.List of symbols  A  Deformation tensor - D_t, D_f  Reactor and injector diameter - L_x, L_y, L_z  Dimensions of the 3D-LIF measurement volume - N_x, N_y, N_z  Number of data samples per measurement volume - Re_m  Reynolds number based on mean velocity - Sc  Schmidt number - f  Focal length - f_c,lens, f_c,array  Cut-off frequency for camera lens and sensor array - f, f  Marginal probability density function for and - f  Joint probability density function of and -  Temporal separation of the 2D data planes -  Temporal resolution of the measurement volume -  Spatial resolution of the measurement volume - ,  Deformation angle and deformation, where =tan -  Fluid energy dissipation rate - ,  Strain limited vorticity and scalar diffusion layers -  Scalar concentration - , _B Kolmogorov and Batchelor length scale - ,  Spherical angles of the scalar gradient vector, -  Kinematic viscosity - _e^–2 Half-thickness (1/e²) of the laser sheet - , _a Kolmogorov and Kolmogorov advection time scales -  Scalar energy dissipation rate -  Scalar diffusivity - 2D, 3D Two- and three-dimensional - DNS Direct numerical simulation - LIF Laser-induced fluorescence - SED Scalar energy dissipation rate - TR Tubular reactor

---

Thermal entrance region heat transfer for MHD laminar flow in parallel-plate channels with unequal wall temperatures

The problem of thermal entry heat transfer for Hartmann flow in parallel-plate channels with uniform but unequal wall temperatures considering viscous dissipation, Joule heating and axial conduction effects is approached by the eigenfunction expansion method. The series expansion coefficients for the nonorthogonal eigenfunctions are obtained by using a method for nonorthogonal series described by Kantorovich and Krylov [21]. Numerical results are obtained for the case with entrance condition parameter _o=1 and open circuit condition K=1. The parametric values of Ha=0, 2, 6, 10 and Br=0, –1 are considered for Hartmann and Brinkman numbers, respectively.

---

Zusammenfassung Das Problem der Wärmeübertragung im thermischen Einlauf einer Hartmannströmung im ebenen Spalt mit einheitlichen, aber ungleichen Wandtemperaturen wurde unter Berücksichtigung viskoser Dissipation, Joulescher Heizung und axialer Wärmeleitung mit Hilfe einer Entwicklung nach Eigenfunktionen behandelt. Die Koeffizienten der Entwicklung für nichtorthogonale Eigenfunctionen wurde nach einer Methode für nichtorthogonale Reihen nach Kantorovicz und Krylow [21] berechnet. Numerische Ergebnisse werden für den Eintrittsparameter _o=1 und die Bedingung für den offenen Stromkreis K=1 erhalten. Die Parameterwerte Ha=0, 2, 6, 10 und Br=0, –1 werden für die jeweiligen Werte der Hartmann- und der Brinckman-Zahl betrachtet.

---

Nomenclature a one-half of channel height - ¯B,B₀ magnetic field Induction vector and magnitude of applied magnetic field - Br Brinkman number, _f U_m ²/(k_c) - C_n,D_n coefficients in the series expansion of _e, see eq. 16 - c_p specific heat at constant pressure - ,E₀ electric field intensity vector and component - E_n,O_n even and odd eigenfunctions - Ha Hartmann number, (/_f)^1/2 B_o a - h₁,h₂ local heat transfer coefficients at lower and upper plates - ¯J,J_y electric current density vector and component - K external loading parameter, E_o/(B_o U_m) - k thermal conductivity - Nu₁, Nu₂ local Nusselt numbers, h₁,a/k and h₂a/k, respectively - P fluid pressure - Pe Peclet number, PrRe - Pr Prandtl number, C_{p f}/k - q₁,q₂ rates of heat transfer per unit area,^–k(T/Z)Z=–a'^–k(T/Z) Z=a respectively - Re Reynolds number, U_ma/u_f - T,T₀,T₁,T₂ fluid temperature, uniform entrance temperature, uniform but different lower and upper plate temperatures, respectively - T_b,T_m bulk temperature and (T₁+T₂)/2 - U,U_m,u axial, mean and dimensionless velocities, respectively - ¯V velocity vector - X,Z axial and transverse coordinates - x,z dimensionless coordinates - _n,_n even and odd eigenvalues - ,₀,_b dimensionless fluid, entrance and bulk temperatures, respectively - _c,_e,_f characteristic temperature difference (T₂-T_m), and dimensionless fluid temperatures, defined by eq. (10) - _e,_f magnetic permeability and viscosity of fluid - fluid density - electric conductivity - viscous dissipation function - (1-)/2     

Dynamic-mechanical properties of a bitumen-silica composite

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

---

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

---

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

---

---

With 9 figures and 2 tables     

Topological entropy of a dynamical system on the space of one-dimensional maps

We study properties of the topological entropy of the map F: f , C(I), generated by a fixed continuous map f C(I) of an interval of the straight line. In particular, we show that the topological entropy h(F) > 0 if and only if h(f) > 0.Translated from Neliniini Kolyvannya, Vol. 7, No. 2, pp. 180–187, April–June, 2004.     

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.     

Parameter identification in a soil with constant diffusivity

We consider infiltration into a soil that is assumed to have hydraulic conductivity of the form K = K = K_se^h and water content of the form = K – _r. Here h denotes capillary pressure head while K_s, , and _r represent soil specific parameters. These assumptions linearize the flow equation and permit a closed form solution that displays the roles of all the parameters appearing in the hydraulic function K and . We assume K_s and _r to be known. A measurement of diffusivity fixes the product of and resulting in a parameter identification problem for one parameter. We show that this parameter identification problem, in some cases, has a unique solution. We also show that, in some cases, this parameter identification problem can have multiple solutions, or no solution. In addition it is shown that solutions to the parameter identification problem can be very sensitive to small changes in the problem data.     

A rotating hot-wire technique for spatial sampling of disturbed and manipulated duct flows

The documentation and control of flow disturbances downstream of various open inlet contractions was the primary focus with which to evaluate a spatial sampling technique. An X-wire probe was rotated about the center of a cylindrical test section at a radius equal to one-half that of the test section. This provided quasi-instantaneous multi-point measurements of the streamwise and azimuthal components of the velocity to investigate the temporal and spatial characteristics of the flowfield downstream of various contractions. The extent to which a particular contraction is effective in controlling ingested flow disturbances was investigated by artificially introducing disturbances upstream of the contractions. Spatial as well as temporal mappings of various quantities are presented for the streamwise and azimuthal components of the velocity. It was found that the control of upstream disturbances is highly dependent on the inlet contraction; for example, reduction of blade passing frequency noise in the ground testing of jet engines should be achieved with the proper choice of inlet configurations.List of symbols K _uv correlation coefficient= - P percentage of time that an azimuthal fluctuating velocity derivative dv/d is found - U streamwise velocity component U=U (, t) - V azimuthal or tangential velocity component due to flow and probe rotation V=V (, t) - mean value of streamwise velocity component - U _m resultant velocity from and - mean value of azimuthal velocity component induced by rotation - u fluctuating streamwise component of velocity u=u(, t) - v fluctuating azimuthal component of velocity v = v (, t) - u phase-averaged fluctuating streamwise component of velocity u=u(0) - v phase-averaged fluctuating azimuthal component of velocity v=v() - û average of phase-averaged fluctuating streamwise component of velocity (u()) over cases I-1, II-1 and III-1 û = û() - average of phase-averaged fluctuating azimuthal component of velocity (v()) over cases I-1, II-1 and III-1 - u fluctuating streamwise component of velocity corrected for non-uniformity of probe rotation and/or phase-related vibration u = u(0, t) - v fluctuating azimuthal component of velocity corrected for non-uniformity or probe rotation and/or phase-related vibration v=v (, t) - u ² rms value of corrected fluctuating streamwise component of velocity - rms value of corrected fluctuating azimuthal component of velocity - phase or azimuthal position of X-probe     

Similar solutions and the asymptotics of filtration equations

In this paper we consider the asymptotic behavior of solutions of the quasilinear equation of filtration as t. We prove that similar solutions of the equation u _t = (u )_xx asymptotically represent solutions of the Cauchy problem for the full equation u _t = [(u)]_xx if (u) is close to u for small u.     

Modelling of heat transfer by conduction in a transition region between a porous medium and an external fluid

We study the modelling of purely conductive heat transfer between a porous medium and an external fluid within the framework of the volume averaging method. When the temperature field for such a system is classically determined by coupling the macroscopic heat conduction equation in the porous medium domain to the heat conduction equation in the external fluid domain, it is shown that the phase average temperature cannot be predicted without a generally negligible error due to the fact that the boundary conditions at the interface between the two media are specified at the macroscopic level.Afterwards, it is presented an alternative modelling by means of a single equation involving an effective thermal conductivity which is a function of point inside the interfacial region.The theoretical results are illustrated by means of some numerical simulations for a model porous medium. In particular, temperature fields at the microscopic level are presented.Roman Letters _sf interfacial area of thes-f interface contained within the macroscopic system m² - A _sf interfacial area of thes-f interface contained within the averaging volume m² - C _p mass fraction weighted heat capacity, kcal/kg/K - g vector that maps to _s , m - h vector that maps to _f , m - K _eff effective thermal conductivity tensor, kcal/m s K - l _s,l _f microscopic characteristic length m - L macroscopic characteristic length, m - n _fs outwardly directed unit normal vector for thef-phase at thef-s interface - n outwardly directed unit normal vector at the dividing surface. - R ₀ REV characteristic length, m - T _i macroscopic temperature at the interface, K - error on the external fluid temperature due to the macroscopic boundary condition, K - T ^* macroscopic temperature field obtained by solving the macroscopic Equation (3), K - V averaging volume, m³ - V _s,V _f volume of the considered phase within the averaging volume, m³. - _mp volume of the porous medium domain, m³ - _ex volume of the external fluid domain, m³ - _s , _f volume of the considered phase within the volume of the macroscopic system, m³ - dividing surface, m² - x, z spatial coordinates Greek Letters _s, _f volume fraction - ratio of the effective thermal conductivity to the external fluid thermal conductivity - ^* macroscopic thermal conductivity (single equation model) kcal/m s K - _s, _f microscopic thermal conductivities, kcal/m s K - spatial average density, kg/m³ - microscopic temperature, K - ^* microscopic temperature corresponding toT ^*, K - spatial deviation temperature K - error in the temperature due to the macroscopic boundary conditions, K - ^* _i macroscopic temperature at the interface given by the single equation model, K - spatial average - ^s , ^f intrinsic phase average.     

The Conjugacy of Stochastic and Random Differential Equations and the Existence of Global Attractors

We consider stochastic differential equations in d-dimensional Euclidean space driven by an m-dimensional Wiener process, determined by the drift vector field f₀ and the diffusion vector fields f₁,...,f_m, and investigate the existence of global random attractors for the associated flows . For this purpose is decomposed into a stationary diffeomorphism given by the stochastic differential equation on the space of smooth flows on R^d driven by m independent stationary Ornstein Uhlenbeck processes z¹,...,z^m and the vector fields f₁,...,f_m, and a flow generated by the nonautonomous ordinary differential equation given by the vector field (_t/x)^–1[f₀(_t)+ _i=1 ¹ f_i(_t)z _t ⁱ ]. In this setting, attractors of are canonically related with attractors of . For , the problem of existence of attractors is then considered as a perturbation problem. Conditions on the vector fields are derived under which a Lyapunov function for the deterministic differential equation determined by the vector field f₀ is still a Lyapunov function for , yielding an attractor this way. The criterion is finally tested in various prominent examples.     

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     

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

Squeezing flows of Newtonian liquid films an analysis including fluid inertia

A theoretical study is made of the flow behavior of thin Newtonian liquid films being squeezed between two flat plates. Solutions to the problem are obtained by using a numerical method, which is found to be stable for all Reynolds numbers, aspect ratios, and grid sizes tested. Particular emphasis is placed on including in the analysis the inertial terms in the Navier-Stokes equations.Comparison of results from the numerical calculation with those from Ishizawa's perturbation solution is made. For the conditions considered here, it is found that the perturbation series is divergent, and that in general one must use a numerical technique to solve this problem.Nomenclature a half of the distance, or gap, between the two plates - a ₀ the value of a at time t=0 - adot da/dt - ä d² a/dt ² - d³ a/dt ³ - a ⁱ components of a contravariant acceleration vector - f unknown function of z ₀ and t defined in (6) - f ⁱ function defined in (9) f ¹=r ₀ g(z ₀, t) f ²= ₀ f ³=f(z ₀, t) - F force applied to the plates - g unknown function of z ₀ and t defined in (6) - g g/z ₀ - h grid dimension in the z ₀ direction (see Fig. 5) - Christoffel symbol - i, j, k, l indices - k grid dimension in the t direction (see Fig. 5) - r radial coordinate direction defined in Fig. 1 - r ₀ radial convected coordinate - R radius of the circular plates - t time - v _r fluid velocity in the r direction - v _z fluid velocity in the z direction - v fluid velocity in the direction - x ⁱ cylindrical coordinate x ¹=r x²= x³=z - z vertical coordinate direction defined in Fig. 1 - z ₀ vertical convected coordinate - tangential coordinate direction - ₀ tangential convected coordinate - viscosity - kinematic viscosity, / - ⁱ convected coordinate ¹=r₀ ²=₀ ³=z₀ - density     

Existence of L 1-Connections Between Equilibria of a Semilinear Parabolic Equation

We give a necessary and sufficient condition for the existence of L ¹-connections between equilibria of a semilinear parabolic equation. By an L ¹-connection from an equilibrium ^– to an equilibrium ⁺ we mean a function u(, t) which is a classical solution on the interval (–, T) for some T and blows up at t = T but continues to exist in the space L ¹ for t [T, ) and satisfies u(, t) ^± (in a suitable sense) as t ±. The main tool in our analysis is the zero number.