An experimental study of the lateral migration of a droplet in a creeping flow

The distribution of droplets in a plane Hagen-Poiseuille flow of dilute suspensions has been measured by a special LDA technique. This method assumes a well defined relation between the velocity of the droplets and their lateral position in the channel. The measurements have shown that the droplet distribution is non-uniform and depends on the viscosity ratio between the droplets and the carrier liquid. The results have been compared with a theory by Chan and Leal describing the lateral migration of suspended droplets.List of symbols a particle radius, m - d half width of the channel, m - Re flow Reynolds number, = 2 _m · d · /µ - flow velocity, m/s - _m flow velocity at the channel axis, m/s - We Weber number, = 2 _m ^Emphasis>/2 · d · / - x distance from center line (x = 0) of the channel, m - non-dimensional distance from the channel center line, _x ^d - y distance along the channel (y = 0 at channel inlet), m - non-dimensional distance along the channel, = y/2d - non-dimensional, normalized distance along the channel, = · _m · µ/ - interfacial tension, N/m - viscosity ratio of dispersed (droplet) phase to viscosity of continuous phase - µ viscosity of continuous phase, Pa · s - density of continuous phase, kg/m³ - phase density difference, kg/m³ Experiments were performed at Max-Planck-Institut, Göttingen     

Rheological properties of reactive polymer melts.

A new method for describing the rheological properties of reactive polymer melts, which was presented in an earlier paper, is developed in more detail. In particular, a detailed derivation of the equation of a first-order rheometrical flow surface is given and a procedure for determining parameters and functions occurring in this equation is proposed. The experimental verification of the presented approach was carried out using our data for polyamide-6.Notation E Dimensionless reduced viscosity, eq. (34) - E ₀ Newtonian asymptote of the function (36) - E power-law asymptote of the function (36) - E _{= 1} the value ofE at = 1 - k degradation reaction rate constant, s^–1 - k ₁ rate constant of function (t), eq. (26), s^–1 - k ₂ rate constant of function (t), eq. (29), s^–1 - K(t) residence-time-dependent consistency factor, eq. (22) - M _w weight-average molecular weight - M _x x-th moment of the molecular weight distribution - R gas constant - S _x M _x /M _w - t residence time in molten state, s - t _j thej-th value oft, s - T temperature, K - % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xd9vqpe0x% c9q8qqaqFn0dXdir-xcvk9pIe9q8qqaq-xir-f0-yqaqVeLsFr0-vr% 0-vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaaieGaceWFZo% Gbaiaaaaa!3B4E!\[\dot \gamma \] shear rate, s^–1 - _i thei-th value of , s^–1 - _r =1 the value of at = 1, s^–1 - ^* reduced shear rate, eq. (44), s^–1 - dimensionless reduced shear rate, eq. (35) - viscosity, Pa · s - shear-rate and residence-time dependent viscosity, Pa · s - zero-shear-rate degradation curve - degradation curve at - _{t0 (t)} zero-residence-time flow curve - Newtonian asymptote of the RFS - instantaneous flow curve - power-law asymptote of the RFS - _0,0 zero-shear-rate and zero-residence-time viscosity, Pa · s - _E=1 value of viscosity atE=1, Pa · s - ^* reduced viscosity, eq. (43), Pa · s - zero-residence-time rheological time constant, s - density, kg/m³ - (t),(t) residence time functions     

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.     

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

---

Wärmeübertragung in einem axial rotierenden,durchströmten Rohr im Bereich des thermischen Einlaufs

Zusammenfassung Der Einfluß der Rotation auf das Temperaturprofil und die Wärmeübergangszahl einer turbulenten Rohrströmung im Bereich des thermischen Einlaufs wird theoretisch untersucht und mit Meßwerten verglichen. Es wird angenommen, daß das Geschwindigkeitsprofil voll ausgebildet ist. Die Rotation hat aufgrund der radial ansteigenden Zentrifugalkräfte einen ausgeprägten Einfluß auf die Unterdrückung der turbulenten Bewegung. Dadurch verschlechtert sich die Wärmeübertragung mit steigender Rotations-Reynoldszahl und die thermische Einlauflänge nimmt beträchtlich zu.

---

Heat transfer in an axially rotating pipe in the thermal entrance region. Part 1: Effect of rotation on turbulent pipe flow

The effects of rotation on the temperature distribution and the heat transfer to a fluid flowing inside a tube are examined by analysis in the thermal entrance region. The theoretical results are compared with experimental findings. The flow is assumed to have a fully developed velocity profile. Rotation was found to have a very marked influence on the suppression of the turbulent motion because of radially growing centrifugal forces. Therefore, a remarkable decrease in heat transfer with increasing rotational Reynolds number can be observed. The thermal entrance length increases remarkably with growing rotational Reynolds number.

---

Formelzeichen a Temperaturleitzahl - C _n , ,C ₁,C ₃ Konstanten - c _p spezifische Wärme bei konstantem Druck - D Rohrdurchmesser - E Funktion nach Gl. (30) - H _n Eigenfunktionen - l hydrodynamischer Mischungsweg - l _q thermischer Mischungsweg - Massenstrom - N=Re /Re Reynoldszahlenverhältnis - Nu Nusseltzahl - Nu Nusseltzahl für die thermisch voll ausgebildete Strömung - Pr Prandtlzahl - Pr _t turbulente Prandtlzahl - Wärmestromdichte - Re _* Schubspannungsreynoldszahl - R _n Eigenfunktionen - Durchfluß-Reynoldszahl - Re v =D/ Rotations-Reynoldszahl - Ri Richardsonzahl - R Rohrradius - r Koordinate in radialer Richtung - dimensionslose Koordinate in radialer Richtung - T Temperatur - T Temperaturschwankung - T _b bulk temperature - mittlere Axialgeschwindigkeit - v Geschwindigkeit - v Geschwindigkeitsschwankung - turbulenter Wärmestrom - dimensionsloser Wandabstand - =1/6 Konstante - Integrationsvariable - Integrationsvariable - , ₁, ₂, dimensionslose Temperaturen - Wärmeleitzahl - _n Eigenwerte - kinematische Viskosität - Dichte - tangentiale Koordinate - , Hilfsfunktionen Indizes m in der Rohrmitte - r radial - w an der Rohrwand - z axial - 0 am Rohreintritt - 0 ohne Rotation - tangential     

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.     

Simultaneous measurements of velocity and temperature fluctuations in thermal boundary layer in a drag-reducing surfactant solution flow

The mechanism of turbulent heat transfer in the thermal boundary layer developing in the channel flow of a drag-reducing surfactant solution was studied experimentally. A two-component laser Doppler velocimetry and a fine-wire thermocouple probe were used to measure the velocity and temperature fluctuations simultaneously. Two layers of thermal field were found: a high heat resistance layer with a high temperature gradient, and a layer with a small or even zero temperature gradient. The peak value of was larger for the flow with the drag-reducing additives than for the Newtonian flow, and the peak location was away from the wall. The profile of was depressed in a similar manner to the depression of the profile of in the flow of the surfactant solution, i.e., decorrelation between v and compared with decorrelation between u and v. The depression of the Reynolds shear stress resulted in drag reduction; similarly, it was conjectured that the heat transfer reduction is due to the decrease in the turbulent heat flux in the wall-normal direction for a flow with drag-reducing surfactant additives.List of symbols ensemble averaged value - (·)⁺ normalized by the inner wall variables - (·) root-mean-square value - C concentration of cetyltrimethyl ammonium chloride (CTAC) solution - c _p heat capacity - D hydraulic diameter - f friction factor - H channel height - h heat transfer coefficient - j _H Colburn factor - l length - Nu Nusselt number, h - Pr Prandtl number, c _p/ - q _w wall heated flux - Re Reynolds number, U _b/ - T temperature - T _b bulk temperature - T _i inlet temperature - T _w wall temperature - T friction temperature, q _w /c _p u - U local time-mean streamwise velocity - U ₁ velocity signals from BSA1 - U ₂ velocity signals from BSA2 - U _b bulk velocity - u streamwise velocity fluctuation - u1 velocity in abscissa direction in transformed coordinates - u friction velocity, - v wall-normal velocity fluctuation - v1 velocity in ordinate direction in transformed coordinates - var(·) variance - x streamwise direction - y wall-normal direction - z spanwise direction - _j junction diameter of fine-wire TC - _w wire diameter of fine-wire TC - angle of principal axis of joint probability function p(u,v) - _f heat conduction of fluid - _w heat conduction of wire of fine-wire TC - kinematic viscosity - local time-mean temperature difference, T _w –T - temperature fluctuation - standard deviation - density - _w wall shear stress     

Laminar source flow between parallel porous disks with equal suction and injection

An analysis is presented for laminar source flow between parallel stationary porous disks with suction at one of the disks and equal injection at the other. The solution is in the form of an infinite series expansion about the solution at infinite radius, and is valid for all suction and injection rates. Expressions for the velocity, pressure, and shear stress are presented and the effect of the cross flow is discussed.Nomenclature a distance between disks - A, B, ..., J functions of R _w only - F static pressure - p dimensionless static pressure, p(a ²/ ²) - Q volumetric flow rate of the source - r radial coordinate - r dimensionless radial coordinate, r/a - R radial coordinate of a point in the flow region - R dimensionless radial coordinate of a point in the flow region, R - Re source Reynolds number, Q/2a - R _w wall Reynolds number, Va/ - reduced Reynolds number, Re/r ² - critical Reynolds number - velocity component in radial direction - u dimensionless velocity component in radial direction, a/ - average radial velocity, Q/2a - u dimensionless average radial velocity, Re/r - ratio of radial velocity to average radial velocity, u/u - velocity component in axial direction - v dimensionless velocity component in axial direction, v - V magnitude of suction or injection velocity - z axial coordinate - z dimensionless axial coordinate, z a - viscosity - density - kinematic viscosity, / - shear stress at lower disk - shear stress at upper disk - ₀ dimensionless shear stress at lower disk, - ₁ dimensionless shear stress at upper disk, - dimensionless stream function     

Concentration-dependent behaviour of the shear viscosity of coal-fuel oil suspensions

A function correlating the relative viscosity of a suspension of solid particles in liquids to their concentration is derived here theoretically using only general thermodynamic ideas, with out any consideration of microscopic hydrodynamic models. This function ( _r = exp (1/2B ^* C ²)) has a great advantage over the many different functions proposed in literature, for it depends on a single parameter,B ^*, and is therefore concise. To test the validity of this function, a least-squares regression analysis was undertaken of available data on the viscosity and concentration of suspensions of coal particles in fuel oil, which promise to be a useful alternative to fuel oil in the near future. The proposed function was found to accurately describe the concentration-dependent behaviour of the relative viscosity of these suspensions. Furthermore, an attempt was made to obtain information about the factors affecting the value ofB ^*, however the results were only qualitative because of, among other things, the inaccuracy of the viscosity measurements in such highly viscous fluids. shear viscosity of the suspension - ₀ shear viscosity of the Newtonian suspending medium - _r = /₀ relative viscosity - solid volume concentration - c solid weight concentration - _m maximum attainable volume concentration of solids - solid volume concentration at which the relative viscosity of the suspension becomes infinite - c _m maximum attainable solid weight concentration - _s density of the solid phase - _l density of the liquid phase - _m density of the suspension - k _n coefficients of theø-power series expansion of _r - { _j } sets of parameters specifying the thermodynamic state of the solid phase of a suspension - T absolute temperature (K) - f (c, T, _j) formal expression for the relative variation of the viscosity with concentration = [1 / (/c)] _T,j - d median size of the granulometric distribution - _B plastic or Bingham viscosity - K consistency factor - n flow index - g ([c _m –c],T, _j ) function including an asymptotic divergence asc tends toc _m , formally describing the concentration dependent behaviour of the shear viscosity of a suspension - A (T, _j) regression analysis parameters - B (T, _j) regression analysis parameters - B ^* (T, _j ) regression analysis parameters     

Transport in ordered and disordered porous media II: Generalized volume averaging

In this paper we develop the averaged form of the Stokes equations in terms of weighting functions. The analysis clearly indicates at what point one must choose a media-specific weighting function in order to achieve spatially smoothed transport equations. The form of the weighting function that produces the cellular average is derived, and some important geometrical theorems are presented.Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - A _p surface area of a particle, m² - d _p 6V _p/A_p, effective particle diameter, m - g gravity vector, m/s² - I unit tensor - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p general characteristic length for volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - l characteristic length (pore scale) for the-phase - l _i i=1, 2, 3 lattice vectors, m - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - m _D special convolution product weighting function for disordered media - m _M master convolution product weighting function for ordered and disordered media - n unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p_m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p p p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function, m - r position vector, m - r position vector locating points in the-phase, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - V velocity vector in the-phase, m/s - vm superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - V volume of the-phase contained in the averaging volume, m³ - V _p volume of a particle, m³ - v traditional superficial volume averaged velocity, m/s - v v p _m spatial deviation velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - _m m * . weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m² - V /V, volume fraction of the-phase     

Wärmeübertragung in einem axial rotierenden,durchströmten Rohr im Bereich des thermischen Einlaufs

Zusammenfassung Der Einfluß der Rotation auf das Temperaturprofil und die Wärmeübergangszahl einer laminaren Rohrströmung im Bereich des thermischen Einlaufs wird theoretisch untersucht. Es wird angenommen, daß das Geschwindigkeitsprofil voll ausgebildet ist. Die Rotation hat einen destabilisierenden Einfluß auf die Laminarströmung, die umschlägt und turbulent wird. Aufgrund der Anfachung der Turbulenz durch die Rotation verbessert sich die Wärmeübertragung mit steigender Rotations-Reynoldszahl und die thermische Einlauflänge nimmt beträchtlich ab.

---

Heat transfer in an axially rotating pipe in the thermal entrance region. Part 2: Effect of rotation on laminar pipe flow

The effects of tube rotation on the temperature distribution and the heat transfer to a fluid flowing inside a tube are examined by analysis in the thermal entrance region. The flow is assumed to be hydrodynamically fully developed. The rotation has a destabilizing effect on the laminar pipe flow, causing a transition to turbulent flow. Therefore, a remarkable increase in heat transfer with increasing rotational Reynolds number can be observed. The thermal entrance length decreases remarkably with growing rotational Reynolds number.

---

Formelzeichen a Temperaturleitzahl - C _n , ,C ₁,C ₃ Konstanten - c _p spezifische Wärme bei konstantem Druck - D Rohrdurchmesser - E Funktion nach Gl. (16) - H _n Eigenfunktionen - l hydrodynamischer Mischungsweg - l _q thermischer Mischungsweg - N=Re /Re Reynoldszahlenverhältnis - Nu Nusseltzahl - Nu Nusseltzahl für die thermisch voll ausgebildete Strömung - Pr Prandtlzahl - Pr _t turbulente Prandtlzahl - Wärmestromdichte - Re _* Schubspannungsreynoldszahl - R _s Eigenfunktionen - Durchfluß-Reynoldszahl - Re =v D/v Rotations-Reynoldszahl - R Rohrradius - r Koordinate in radialer Richtung - Dimensionslose Koordinate in radialer Richtung - T Temperatur - T _b bulk temperature - mittlere Axialgeschwindigkeit - v Geschwindigkeit - dimensionsloser Wandabstand - Integrationsvariable - Integrationsvariable - , dimensionslose Temperaturen - Wärmeleitzahl - _p Eigenwerte - kinematische Viskosität - Dichte - tangentiale Koordinate - , Hilfsfunktionen Indizes m in der Rohrmitte - r radial - w an der Rohrwand - z axial - 0 am Rohreintritt - 0 ohne Rotation - tangential     

On the universality of the velocity profiles of a turbulent flow in an axially rotating pipe

If a fluid enters an axially rotating pipe, it receives a tangential component of velocity from the moving wall, and the flow pattern change according to the rotational speed. A flow relaminarization is set up by an increase in the rotational speed of the pipe. It will be shown that the tangential- and the axial velocity distribution adopt a quite universal shape in the case of fully developed flow for a fixed value of a new defined rotation parameter. By taking into account the universal character of the velocity profiles, a formula is derived for describing the velocity distribution in an axially rotating pipe. The resulting velocity profiles are compared with measurements of Reich [10] and generally good agreement is found.Nomenclature b constant, equation (34) - D pipe diameter - l mixing length - l ₀ mixing length in a non-rotating pipe - N rotation rate,N=Re /Re _D - p pressure - R pipe radius - Re _D flow-rate Reynolds number, - Re rotational Reynolds number, Re =v _w D/ - Re* Reynolds number based on the friction velocity, Re*=v*R/ - (Re*)₀ Reynolds number based on the friction velocity in a non-rotating pipe - Ri Richardson number, equation (10) - r coordinate in radial direction - dimensionless coordinate in radial direction, - v _r ,v ,v _z time mean velocity components - v _r ,v ,v _z velocity fluctations - v _w tangential velocity of the pipe wall - v* friction velocity, - axial mean velocity - v _ZM maximum axial velocity - dimensionless radial distance from pipe wall, - y ⁺ dimensionless radial distance from pipe wall - y ₁ ⁺ constant - Z rotation parameter,Z =v _w/v _* =N Re _D /2Re* - _m eddy viscosity - ( _m )₀ eddy viscosity in a non-rotating pipe - coefficient of friction loss - von Karman constant - ₁ constant, equation (31) - density - dynamic viscosity - kinematic viscosity     

Slip flow of non-plasticized PVC compounds

Measurements on seven rigid PVC compounds were carried out with a slit rheometer working in combination with an injection moulding machine. Plastication of the compounds occurred in the screw of the plastication unit, which also forced the melt through the die with a controlled forward velocity. The rectangular slit had a length of 90 mm and a widthB of 20 mm. The heightH could be varied between 0.8 and 3.3 mm. Pressures and temperatures were recorded at several positions in and before the die. Measurements were carried out at shear rates from 10 to 2000 s^–1.When the reduced volume output was plotted against the wall shear stress _W , only four compounds showed master curves independent ofH, which is indicative of wall adhesion. In the other cases this plot did not produce such a master curve, but the plot of the mean velocity against _W was independent ofH (slip curve). This indicated that slip flow prevailed with a slip velocityv _G When, in the case of wall slip, the smooth inner surfaces of the die were replaced by surfaces with grooves perpendicular to the direction of flow, slip flow was prevented and the flow curves were shifted to much higher values of _Wc Above a critical value of the wall shear stress ( _Wc ) at which slip flow began, the output became nearly independent of _W . From the measurements made below _Wc a vs. relation for the shear flow could be derived, which was used to calculate the superimposed shear flow . Exact values of the slip velocity were then given by . Wall slip only occurred for compounds with a high shear viscosity, which corresponds to a high molecular weight (K-value).Dedicated to Professor H. Janeschitz-Kriegl on the occasion of his 60th birthday.     

Injection moulding of thermoplastics: Freezing of variable-viscosity fluids.

The injection moulding of thermoplastics involves, during mould filling, flows of hot polymer melts into mould networks, the walls of which are so cold that frozen layers form on them. An analytical study of such flows is presented here for the case when the Graetz and Nahme numbers are large and the Pearson number is small. Thus the flows are developing and temperature differences due to heat generation by viscous dissipation are sufficiently large to cause significant variations in viscosity (but the difference between the entry temperature of the polymer to a specific part of the mould network and the melting temperature of the polymer is not). Br Brinkman number - Gz Graetz number - h half-height of channel or disc - h ^* half-height of polymer melt region in channel or disc - L length of channel or pipe - m viscosity shear-rate exponent - Na Nahme number - p pressure - P pressure drop - Pe Péclet number - Pn Pearson number - Q volumetric flowrate - r radial coordinate in pipe or disc - R radius of pipe - Re Reynolds number - R _i inner radius of disc - R _o outer radius of disc - R ^* radius of polymer melt region in pipe - T temperature - T _ad adiabatic temperature rise - T _e entry polymer melt temperature - T _m melting temperature of polymer - T _max maximum temperature - T ₀ reference temperature - T _w wall temperature - flow-average temperature rise - u _r radial velocity in pipe or disc - u _x axial velocity in channel - u _y transverse velocity in channel or disc - u _z axial velocity in pipe - w width of channel - x axial coordinate in channel or modified radial coordinate in disc - y transverse coordinate in channel or disc - z axial coordinate in pipe - thermal conductivity of molten polymer - thermal conductivity of frozen polymer - scaled dimensionless axial coordinate in channel or pipe or radial coordinate in disc - ₀ undetermined integration constant - heat capacity of molten polymer - viscosity temperature exponent - dimensionless transverse coordinate in channel or disc - ^* dimensionless half-height of polymer melt region in channel or disc - H ^* scaled dimensionless half-height of polymer melt region in channel or disc or radius of polymer melt region in pipe - dimensionless temperature - ^* dimensionless wall temperature - scaled dimensionless temperature - numerical constant - µ viscosity of molten polymer - µ ₀ consistency of molten polymer - dimensionless pressure gradient - scaled dimensionless pressure gradient - density of molten polymer - dimensionless radial coordinate in pipe or disc - _i dimensionless inner radius of disc - ^* dimensionless radius of polymer melt region in pipe - dimensionless streamfunction - scaled dimensionless streamfunction - dummy variable - streamfunction - similarity variable - similarity variable     

Stoffaustausch zwischen zwei fluiden Phasen,von denen sich eine als Strahl in der anderen bewegt

Zusammenfassung Der Übergang eines Stoffes zwischen zwei fluiden Phasen wird betrachtet, von denen sich einer als Strahl in der anderen bewegt. Die Geschwindigkeit der laminar strömenden Phase wird durch eine Gleichung ausgedrückt, die Geschwindigkeitsprofile zwischen der Kolben- und der Rohrströmung kontinuierlich beschreibt. Der Transport des Stoffes im Strahl durch Diffusion in radialer und durch Konvektion in axialer Richtung wird für den isothermen, stationären Fall untersucht. Die das Problem beschreibende Differentialgleichung wird anscheinend erstmals geschlossen gelöst. Die Lösungen beinhalten konfluente hypergeometrische Funktionen. Berechnet werden Eigenwerte, Koeffizienten, örtliche und mittlere Konzentrationsfelder sowie Stoffübergangszahlen.

---

Mass transfer between two fluids, one of the two fluids is moving as jet within the other

The mass transfer between two fluids is calculated, one of the two fluids is moving as a jet within the other. The velocity of the laminar flowing phase is expressed by an equation, which describes continously the velocity profiles from plug flow to tubular flow. For the isothermal, stationary state the transport of substance i by radial diffusion and by axial convection is investigated. It appears to be that the differential equations describing the problem are solved rigorously for the first time. The solutions contain confluent hypergeometrical functions. Results include eigenvalues, coefficients, local and mean concentration fields, mass transfer numbers.

---

Verwendete Zeichen und ihre Bedeutung a - A, A_n Koeffizienten - B, B_n Koeffizienten - c Konzentration, Konstante im Anhang - C_r=0 Mittenkonzentration - c₀ Konzentration in Phase I bis z=0 - c_II Konzentration in Phase II - ¯c mittlere Konzentration, definiert in Gl. (35) - C Koeffizient, definiert in Gl. (A 21) - D Diffusionskoeffizient - Da Damköhlerzahl - E Funktion, gegeben durch Gl. (A 12) - f, f(R) Funktion f von R - f_n, f_n (R) Funktionswerte - g, g(Z) Funktion g von Z - g_n, g_n (Z) Funktionswerte - h(Z) Funktion h von z - H_q Koeffizienten, gegeben durch Gl. (A 10) - j Massenstromdichte - J _k , J_q Besselfunktion der Ordnungk, q - k definiert durch Gl. (A 9) - n laufende Zahl - m laufende Zahl - p laufende Zahl - Pe=Re·Sc Pecletzahl - q laufende Zahl - Q_n Koeffizienten, definiert in Gl. (31) - r radiale Koordinate - r₀ Radius - R r/r₀ - Re=u₀r₀/ Reynoldszahl - S=2r₀z Zylinderfläche - Sc=/D Schmidtzahl - Sh=2r₀ /D Sherwoodzahl - Sherwoodzahl, definiert in Gl. (52) - Sh_u Sherwoodzahl, definiert in Gl. (54) - Sh_z Sherwoodzahl, definiert in Gl. (40) - Sherwoodzahl, definiert in Gl. (45) - t R² - u Geschwindigkeit - u₀ maximale Geschwindigkeit - v - Volumenstrom - w Variable - x Variable - y abhängige Variable - z axiale Koordinate, Lauflänge - Z z/r₀ - Z_Pe dimensionslose Lauflänge, definiert durch Gl. (34) - a_n Koeffizienten, definiert durch Gl. (A 19) - Stoffübergangskoeffizient - Stoffübergangskoeffizient, definiert in Gl. (48) - _u Stoffübergangskoeffizient, definiert in Gl. (49) - _z Stoffübergangskoeffizient, definiert in Gl. (38) - Stoffübergangskoeffizient, definiert in Gl. (44) - definiert in Gl. (A 21) - Gammafunktion - c Konzentrationsdifferenz - m Stoffmenge - Zahl zwischen Null und Eins - laufende Zahl - kinematische Zähigkeit - (v) (t) - konfluente hypergeometrische Funktion - (t) - konfluente hypergeometrische Funktion - , _n Eigenwerte Hochzeichen - * kennzeichnet asymptotische Lösungen     

The unsteady aerodynamics of a first stage stator vane row

The fundamental unsteady aerodynamics on a vane row of an axial flow research compressor stage are experimentally investigated, demonstrating the effects of airfoil camber and steady loading. In particular, the rotor wake generated unsteady surface pressure distributions on the first stage vane row are quantified over a range of operating conditions. These cambered airfoil unsteady data are correlated with predictions from a flat plate cascade inviscid flow model. At the design point, the unsteady pressure difference coefficient data exhibit good correlation with the nonseparated predictions, with the aerodynamic phase lag data exhibiting fair trendwise correlation. The quantitative phase lag differences are associated with the camber of the airfoil. An aft suction surface flow separation region is indicated by the steady state surface static pressure data as the aerodynamic loading is increased. This separation affects the increased incidence angle unsteady pressure data.List of symbols b airfoil semi-chord - C airfoil chord - C _p dynamic pressure coefficient, - _p static pressure coefficient, - i incidence angle - k reduced frequency, - N number of rotor revolutions - p dynamic pressure difference - static pressure difference, - S stator vane circumferential spacing - U _t rotor blade tip speed - u longitudinal perturbation velocity - V absolute velocity - V _axial absolute axial velocity - v transverse perturbation velocity - x _sep location of separation point - inlet angle - inlet air density - blade passing angular frequency     

The effective permeability of a heterogeneous porous medium

The effective (single-phase) permeability of an (infinite) heterogeneous porous medium is studied using a formalism of Green's functions. We give formal expressions for it in the form of a series expansion involving the microscopic random-permeability field many-body correlation functions of higher and higher order.The particular case of a log-normal medium of infinite extent is studied using field-theoretical methods. Using partial series resummation techniques, we derivea formula up to all orders in the local correlations which was first reckoned by many authors by means of a first-order calculation. The formula — which remains an approximation — works whatever the dimensionality of the space, and gives the following simple estimate for the effective permeability in 3 D:K _eff=k ^1/33. The method is general and the approximations can be systematically improved on when more complex situations are studied.Roman Letters D number of dimensions of the space in which the flow takes place - f(r) body force field,N - f(q) Fourier-transformed body-force field, Nm³ - G ₀(r, r) Green's function of the Laplace operator, m^–1 - g(k,r, r) velocity propagator before averaging, m^–1 - G(r, r) velocity propagator after averaging, m^–1 - j(r) a scalar dimensionless field - k(r) local value of the permeability at point r, m² - K _eff effective permeability - K _g geometric average of the local permeability, m² - l typical size of the averaging volume, m - L characteristic length of the porous medium or of the reservoir, m - L(r, r) projection operator, m^–2 - M(r, r) scattering operator, m^–3 - p(r) local value of the pressure, Nm^–2 - p(k,r, r) pressure propagator before averaging, m^–1 - P(r, r) pressure propagator after averaging, m^–1 - r position vector, m - r modulus of vectorr, m - unit vector pointing in the direction ofr - q Fourier wave vector, m^–1 - q modulus of the Fourier wave-vectorq, m^–1 - unit vector pointing in the direction ofq - projector over vector - 1 unit tensor - X(r) a local random variable - ¯X(r) volume averaged local random variable - X (r) ensemble averaged local random variable - V large-scale averaging volume, m³ - Z(j) generating functional of a random field - Z(r,j) modified generating functional of a random field - Z normalization factor Greek Letters ₀ average value of the logarithm of the permeability - (r) fluctuation of the logarithm of permeability at pointr - viscosity of the fluid, Nt/m² - (r–r) two-point correlation function of the fluctuations of the logarithm of the permeability - _k correlation length of the permeability correlation function, m - _u correlation length of the velocity correlation function, m     

An approximate theoretical analysis and experimental verification of turbulent entrance region flow of drag reducing fluids

Summary Entry lengths for pipe flows of moderately drag reducing fluids are determined using momentum integral technique. It is shown theoretically that the entry lengths for drag reducing fluids could be significantly larger than the Newtonian fluids flowing turbulently under otherwise identical conditions. The experimental data from the literature bear out the theoretical calculations.

---

Zusammenfassung Mit Hilfe der Impuls-Methode wird die Einlauflänge in einer Rohrströmung für Flüssigkeiten mit mäßig starker Widerstandsverminderung berechnet. Es wird vorausgesagt, daß die Einlauflänge für derartige Flüssigkeiten erheblich größer sein kann als für newtonsche Flüssigkeiten unter sonst identischen Bedingungen. Aus der Literatur entnommene experimentelle Daten bestätigen diese theoretischen Berechnungen.

---

Nomenclature A ₁ Coefficient in eq. [7] - A Slope of logarithmic velocity profile - a Exponent in eq. [10] - B Intercept function for logarithmic velocity profile - De Deborah number, - f Friction factor - F Function, eq. [30] - G Function given in eq. [11] - Static pressure, dynes/cm² - q Index of power law velocity profile - R Pipe radius, cm - r Radial distance, cm - R Core radius, cm - Re Reynolds number - Axial velocity, cm/s - u _c Core velocity, cm/s - u ⁺ Dimensionless velocity, eq. [5] - u ^* Friction velocity, , cm/s - Radial velocity, cm/s - V Average velocity, cm/s - x Axial distance, cm - x _e Entry length, cm - y Distance from the wall, cm - y ⁺ Dimensionless distance, eq. [5] - y _I ⁺ Dimensionless viscous sublayer thickness - coefficient in eq. [17] - exponent of Reynolds number in eq. [17] - shear rate, s^–1 - boundary layer thickness, cm - _fl fluid relaxation time, s - µ fluid viscosity, gm/cm s - v kinematic viscosity, cm²/s - _l laminar sublayer thickness, dimensionless - fluid density, gm/cm³ - shear stress, dynes/cm² - _w shear stress at the wall, dynes/cm² - ₁, ₂, ₃, ₄ functions in eq. [27] - ~ time averaged quantities - — dimensionless quantity With 3 figures and 1 table     

Transient Free Convection from a Horizontal Surface in a Porous Medium Subjected to a Sudden Change in Surface Heat Flux

This paper presents an analytical and numerical study of transient free convection from a horizontal surface that is embedded in a fluid-saturated porous medium. It is assumed that for time steady state velocity and temperature fields are obtained in the boundary-layer which occurs due to a uniform flux dissipation rate q ₁ on the surface. Then, at the heat flux on the surface is suddenly changed to q ₂ and maintained at this value for . Firstly, solutions which are valid for small and large are obtained. The full boundary-layer equations are then integrated step-by-step for the transient regime from the initial unsteady state ( ) until such times at which this forward marching approach is no longer well posed. Beyond this time no valid solutions could be obtained which matched the final solution from the forward integration to the steady state profiles at large times .     

An analytical solution for the temperature profiles during injection molding,including dissipation effects

An analytical solution is obtained for the stationary temperature profile in a polymeric melt flowing into a cold cavity, which also takes into account viscous heating effects. The solution is valid for the injection stage of the molding process. Although the analytical solution is only possible after making several (at first sight) rather stringent assumptions, the calculated temperature field turns out to give a fair agreement with a numerical, more realistic approach. Approximate functions were derived for both the dissipation-independent and the dissipation-dependent parts which greatly facilitate the temperature calculations. In particular, a closed-form expression is derived for the position where the maximum temperature occurs and for the thickness of the solidified layer.The expression for the temperature field is a special case of the solution of the diffusion equation with variable coefficients and a source term.Nomenclature a thermal diffusivity [m²/s] - c specific heat [J/kg K] - D channel half-height [m] - L channel length [m] - m 1/ - P pressure [Pa] - T temperature [°C] - T _W wall temperature [°C] - T _i injection temperature [°C] - T _A Br independent part of T - T _B Br dependent part of T - T ^core asymptotic temperature - v _z() axial velocity [m/s] - W channel width [m] - x cross-channel direction [m] - z axial coordinate [m] - (x) gamma function - (a, x) incomplete gamma function - M(a, b, x) Kummer function - small parameter - () temperature function - thermal conductivity [W/mK] - viscosity [Pa · s] - ₀ consistency index - power-law exponent - density [kg/m] - similarity variable Dimensionless variables Br Brinkman number - Gz Graetz number -