Instability of Buckley-Leverett flow in a heterogeneous medium

We study the simultaneous one-dimensional flow of water and oil in a heterogeneous medium modelled by the Buckley-Leverett equation. It is shown both by analytical solutions and by numerical experiments that this hyperbolic model is unstable in the following sense: Perturbations in physical parameters in a tiny region of the reservoir may lead to a totally different picture of the flow. This means that simulation results obtained by solving the hyperbolic Buckley-Leverett equation may be unreliable.Symbols and Notation f fractional flow function varying withs andx - value off outsideI - value off insideI - local approximation off around¯x - f ^–,f ⁺ values of - f _j ⁿ value off atS _j ⁿ andx _j - g acceleration due to gravity [ms^–2] - I interval containing a low permeable rock - k dimensionless absolute permeability - k ^* absolute permeability [m²] - k _c ^* characteristic absolute permeability [m²] - k _ro relative oil permeability - k _rw relative water permeability - L ^* characteristic length [m] - L ₁ the space of absolutely integrable functions - L the space of bounded functions - P _c dimensionless capillary pressure function - P _c ^* capillary pressure function [Pa] - P _c ^* characteristic pressure [Pa] - S similarity solution - S _j ⁿ numerical approximation tos(x_j, t_n) - S ₁, S₂,S ₃ constant values ofs - s water saturation - value ofs at - s ^L left state ofs (wrt. ) - s ^R right state ofs (wrt. ) - s s for a fixed value of in Section 3 - T value oft - t dimensionless time coordinate - t ^* time coordinate [s] - t _c ^* characteristic time [s] - t _n temporal grid point,t _n=n t - v ^* total filtration (Darcy) velocity [ms^–1] - W, , v dimensionless numbers defined by Equations (4), (5) and (6) - x dimensionless spatial coordinate [m] - x ^* spatial coordinate [m] - x _j spatial grid piont,x _j=j x - discontinuity curve in (x, t) space - right limiting value of¯x - left limiting value of¯x - angle between flow direction and horizontal direction - t temporal grid spacing - x spatial grid spacing - length ofI - parameter measuring the capillary effects - argument ofS - _o dimensionless dynamic oil viscosity - _w dimensionless dynamic water viscosity - _c ^* characteristic viscosity [kg m^–1s^–1] - _o ^* dynamic oil viscosity [kg m^–1s^–1] - _w ^* dynamic water viscosity [k gm^–1s^–1] - _o dimensionless density of oil - _w dimensionless density of water - _c ^* characteristic density [kgm^–3] - _o ^* density of oil [kgm^–3] - _w ^* density of water [kgm^–3] - porosity - dimensionless diffusion function varying withs andx - ^* dimensionless function varying with s andx ^* [kg^–1m³s] - _j ⁿ value of atS _j ⁿ andx _j This research has been supported by VISTA, a research cooperation between the Norwegian Academy of Science and Letters and Den norske stats oljeselskap a.s. (Statoil).     

An analytical solution for the velocity and shear rate distribution of non-ideal Bingham fluids in concentric cylinder viscometers

An analytical solution is presented for the calculation of the flow field in a concentric cylinder viscometer of non-ideal Bingham-fluids, described by the Worrall-Tuliani rheological model. The obtained shear rate distribution is a function of the a priori unknown rheological parameters. It is shown that by applying an iterative procedure experimental data can be processed in order to obtain the proper shear rate correction and the four rheological parameters of the Worrall-Tuliani model as well as the yield surface radius. A comparison with Krieger's correction method is made. Rheometrical data for dense cohesive sediment suspensions have been reviewed in the light of this new method. For these suspensions velocity profiles over the gap are computed and the shear layer thicknesses were found to be comparable to visual observations. It can be concluded that at low rotation speeds the actually sheared layer is too narrow to fullfill the gap width requirement for granular suspensions and slip appears to be unavoidable, even when the material is sheared within itself. The only way to obtain meaningfull measurements in a concentric cylinder viscometer at low shear rates seems to be by increasing the radii of the viscometer. Some dimensioning criteria are presented.Notation A, B Integration constants - C Dimensionless rotation speed = µ/_y - c = 2µ - d = ₀ ²–2c_y - f() = (–₀)²+2c(–_y) - r Radius - r _b Bob radius - r _c Cup radius - r _y Yield radius - r ₀ Stationary surface radius - r Rotating Stationary radius - Y ₀ Shear rate parameter = /µ Greek letters Shear rate - = (r _y /r _b )²– 1 - µ Bingham viscosity - µ₀ Initial differential viscosity - µ µ₀-µ - Rotation speed - Angular velocity - Shear stress - _b Bob shear stress - _B Bingham stress - _y (True) yield stress - ₀ Stress parameter = _B +µ Y ₀ - _B - _y      

Elongational flow behavior of viscoelastic liquids: Modelling bubble dynamics with viscoelastic constitutive relations

Summary Earlier parts of this series have described a technique based on the collapse of single bubbles in the fluids for studying the elongational rheology of viscoelastic solutions and melts of moderate viscosities ( ₀ > 10²p) at relatively high strain rates . The present paper describes the modelling of bubble collapse with both rate and integral type constitutive relations using a body coordinate system. Predictions of the stress at the bubble wall as a function of time during collapse from a BKZ model and a modified corotational Maxwell model compared favorably with experimental data for two polymer solutions, 1% polyacrylamide in water/glycerine and 2% hydroxypropyl cellulose in water.

---

Zusammenfassung In vorangehenden Veröffentlichungen dieser Reihe wurde eine Methode beschrieben, mit Hilfe derer man aus dem Zerfall von einzelnen Blasen in einer Flüssigkeit auf die Dehn-Rheologie viskoelastischer Lösungen und Schmelzen mittlerer Viskosität ( ₀ > 10² P) bei relativ hohen Dehngeschwindigkeiten schließen kann. Die vorliegende Untersuchung beschreibt Modelle des Blasenzerfalls mit Hilfe von Stoffgleichungen sowohl vom rate- als auch vom Integral-Typ, wobei ein körperfestes Koordinatensystem benutzt wird. Die Voraussagen der Spannung an der Blasenwand als Funktion der Zeit während des Zerfalls bei Verwendung eines BKZ- und eines modifizierten korotatorischen Maxwell-Modells zeigen eine recht gute Übereinstimmung mit experimentellen Werten, die an zwei Polymerlösungen, nämlich einer 1%igen Polyacrylamid-Lösung in einer Wasser-Glycerin-Mischung und einer 2%igen wäßrigen Hydropropylcellulose, erhalten worden sind.

---

Nomenclature a material constant - b material constant - g metric tensor, space coordinates - m material constant - n material constant - p pressure - P _G pressure within bubble - P _R pressure outside bubble at the wall - P pressure far away from the bubble - R bubble radius - dR/dt - R ₀ initial bubble radius - t time - u velocity - U potential function - Y R/R ₀ Greek symbols covariant body metric tensor - surface tension - rate of deformation matrix, II -second invariant of - strain rate - ₀ zero shear rate viscosity - _e elongational viscosity - _ef effective viscosity - ¹, ², ³ coordinates in body system - ^{1 1}/R ₀ ³ - body stress tensor - density - space stress tensor - relaxation time - _ef effective relaxation time - bubble pressure function, defined in eq. [19] - vorticity tensor With 11 figures and 1 table     

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 -      

Diabatic internal source flow

Summary A study is made of diabatic internal source flow. It is shown that for such a flow self similar solutions to the full equations of motion are possible provided that the transport properties ( and k) are assumed constant. Solutions are presented for several cases representing internal heat generation in a confined supersonic flow.Nomenclature a, b, c, d, e constant exponents - f dimensionless velocity - F dimensionless velocity - g dimensionless temperature - G dimensionless temperature - H ₀ stagnation enthalpy - h dimensionless pressure - k thermal conductivity - M Mach number - M _c centerline Mach number - m dimensionless density - n dimensionless internal heat strength - Pr Prandtl number - p pressure - Q dimensionless internal heat source parameter - strength of internal heat source - R gas constant - Re Reynolds number - r radial coordinate - T temperature - u radial velocity - v tangential velocity - ratio of specific heats - tangential coordinate - coefficient of viscosity - density     

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

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

Convective heat transfer within fibrous insulation slabs

A numerical study of convective heat flow within a fibrous insulating slab is presented. The material is treated as an anisotropic porous medium and the variation of properties with temperature is taken into account. Good agreement is obtained with available experimental data for the same geometry.

---

Zusammenfassung Für den konvektiven Wärmestrom in einem faserförmigen Isolierstoff wird eine numerische Berechnung angegeben. Der Stoff wird als anisotropes poröses Medium mit temperaturabhängigen Stoffwerten angesehen. Die Übereinstimmung mit verfügbaren Versuchswerten ist gut.

---

Nomenclature C_p specific heat of the gas at the mean temperature - Da Darcy number=k_y/H² - Gr^* modified Grashof number=gTHk_y/²= (Grashof number) × (Darcy number) - H thickness of the specimen - P gas pressure - Pr^* modified Prandtl number= C_p/_x - Ra^* modified Rayleigh number=Gr^* Pr^* - R_p ratio of permeabilities=k_y/k_x - R_k ratio of conductivities= _y/_x - T absolute temperature of the gas - t₁ absolute temperature of the hot face - T₂ absolute temperature of the cold face - T_m mean temperature of the gas=(T₁+T₂)/2 - k_x specific permeability of the porous medium along the x-direction - k_y specific permeability of the porous medium along the y-direction - p T/T_m - q exponent - r exponent - u gas velocity along the x-direction - v gas velocity along the y-direction - X_* distance along the x-direction - y_* distance along the y-direction - T temperature difference=t₁–T₂ - thermal coefficient of expansion of the gas - _m thermal coefficient of expansion of the gas at the mean temperature - _* T–T_m - dimensionless temperature= ^*/T - _a apparent thermal conductivity of the porous medium along the x-direction - _al local apparent thermal conductivity of the porous medium along the x-direction - _x thermal conductivity of the porous medium along the x-direction in the absence of convection - _y thermal conductivity of the porous medium along the y-direction in the absence of convection - dynamic viscosity of the gas - _m dynamic viscosity of the gas at the mean temperature - kinematic viscosity of the gas - _m kinematic viscosity of the gas at the mean temperature - density of the gas - _m density of the gas at the mean temperature - ^* stream function at any point - dimensionless stream function= ^*/( _m/_m)     

Sensitivity of three-segment electrodiffusion probes to eddy shedding

The influence of eddy shedding on the instantaneous readings of a three-segment cylindrical electrodiffusion velocity probe was investigated in an immersed jet with a very low turbulence intensity, = 1.2%. The velocity fluctuations measured by the three-segment probe were smaller than 2.6%, and the maximum error in the flow angle estimation was 2. Vortices with the Strouhal frequency were detected by a simple electrodiffusion probe placed downstream of the three-segment probe, but no peaks with this frequency were found on the frequency spectra of the three-segment probe. From the probe response to a stepwise change of the polarization voltage the characteristic times of the transient process were estimated. List of symbols a parameter in Eq. (1) [A s^b m^-b] - A amplitude gain - b parameter in Eq. (1) - c parameter in Eq. (3) [A s^–1/2] - d probe diameter [m] - f frequency [s^–1] - f _s recording frequency [s^–1] - G power spectrum - I _k relative current through k-th segment, Eq. (2) - i total current [A] - i _k current through k-th segment [A] - N number of data samples - Re Reynolds number, - Sr Strouhal number, - t time [s] - t ₀ characteristic transient time [s] - v jet velocity [m s^-1] - v time mean value of velocity [m s^-1] - v _x, y velocity components measured by probe [m s^-1] - var variance, var - dynamic viscosity [Pa s] - density [kg m^-3] - relative deviation, [%] - flow angle, see Fig. 1 - dimensionless frequency For the financial support of this work we express our thanks to the DFG, Bonn. The assistance of Dr. Ondra Wein and Dr. Pavel Mitschka is greatly appreciated.     

Experimentelle Prüfung der Analogie zwischen konvektiver Wärme- und Stoffübertragung bei nichtabgelöster Strömung

Zusammenfassung Stoffübertragung für das System Naphthalin/Luft und Wärmeübertragung an Luft werden an der Platte mit laminarer und turbulenter Grenzschicht, in einem Rechteckkanal im Bereich des thermischen Analufs bei hydraulisch ausgebildeter turbulenter Strömung und am Kreisrohr bzw. Ringspalt bei vollausgebildeter Strömung gemessen. Die bekannten Gesetze bei Wärmeübertragung für Platte, Kreisrohr und Ringspalt in der Schreibweise für Stoffübertragung werden bestätigt. Die Gleichung vonElser für den thermischen Analufvorgang wird den Versuchsergebnissen angepaßt. Der Exponent der Prandtl- bzw. Schmidt-Zahl nimmt im Bereich 0,7<(Pr;Sc)<2,5 je nach Strömungsform Werte zwischen 0,33 und 0,67 an.

---

Mass transfer for the system naphthalin/air and heat transfer with air were measured for the following geometries: a plate with laminar and turbulent boundary-layer, a rectangular channel with fully developed turbulent velocity distribution in the thermal entrance region, a pipe of circular cross-section, an annular both in fully developed turbulent flow. For the plate, pipe and annular, the results of the two measuring methods agree very well and confirm the well known laws of heat transfer.Elser equation for the thermal entrance region is adapted to the results. The exponent of the Prandtl and Schmidt numbers varies in the range of 0,7<(Pr;Sc)<2,5 between the values 0,33 and 0,67 depending on the state of flow.

---

Bezeichnungen A empirische Zahlenkonstante - B empirische Zahlenkonstante - A Austauschgröße für Impuls [kg/ms] - A _q Austauschgröße für Wärme [kg/ms] - A _S Austauschgröße für Stoff [kg/ms] - B Plattenbreite [m] - C _f Widerstandsbeiwert - D Diffusionskoeffizient [m²/s] - F freier Strömungsquerschnitt [m²] - K empirische Zahlenkonstante - K _h Korrekturfaktor für Stefanstrom - L Plattenlänge [m] - M relative Molekülmasse [g/mol] - P Gesamtdruck [N/m²] - R Gaskonstante [Nm/kg grd] - S Oberfläche [m²] - T absolute Temperatur [°K] - U benetzter Umfang [m] - W Strömungswiderstand [N] - a Temperaturleitfähigkeit [m²/s] Exponent - b Exponent - c _p spezifische Wärme [J/kg grd] - d Rohrdurchmesser [m] - hydraulischer Durchmesser [m] - l charakteristische Bezugslänge [m] - Massenstromdichte [kg/m²s] - m Exponent - n Exponent - p Partial- bzw. Dampfdruck [N/m²] - q Wärmestromdichte [W/m²] - t Zeit [s] - u Exponent - mittlere Strömungsgeschwindigkeit [m/s] - y laufende Koordinate [m] - y* mittlere Lauflänge der Grenzschicht [m] - mittlere Wärmeübergangszahl [W/m²grd] - örtliche Wärmeübergangszahl [W/m²grd] - mittlere Stoffübergangszahl [m/s] - örtliche Stoffübergangszahl [m/s] - Temperatur [°C] - Wärmeleitfähigkeit [W/m grd] - kinematische Zähigkeit [m²/s] - Dichte [kg/m³] - Druckverlustbeiwert - P Druckverlust [N/m²] - G Gewichtsverlust [kg] Dimensionslose Kenngrößen Pr=/a Prandtl-Zahl - Sc=/D Schmidt-Zahl - Le=a/D Lewis-Zahl - Pr _t=A /A _q turbulente Prandtl-Zahl - Sc _t=A /A _S turbulente Schmidt-Zahl - Re=ie225-8 ·l/ Reynolds-Zahl - Nu= ·l/ Nusselt-Zahl - Sh _l= · l/D Sherwood-Zahl - St=Nu/Re Pr Stanton-Zahl - St=Sh/ReSc Stanton-Zahl für Stoffübertragung - j _W=St Pr ^1–n Wärmeübertragungskoeffizient - j _S=St Sc ^1–n Stoffübertragungskoeffizient - Tu % Turbulenzgrad Indizes A Stoff eines Zweistoffsystems - L Luft, Plattenlänge - M Mischung, bezogen auf den Massenstrom - N Naphthalin - S Stoffübertragung - W Wand; Wärmeübertragung - a, i außen, innen - 0 Bezugszustand     

Non-isothermal flow of polymer melts in a curved tube

The results of a numerical study (using finite differences) of heat transfer in polymer melt flow is presented. The rheological behaviour of the melt is described by a temperature-dependent power-law model. The curved tube wall is assumed to be at constant temperature. Convective and viscous dissipation terms are included in the energy equation. Velocity, temperature and viscosity profiles, Nusselt numbers, bulk temperatures, etc. are presented for a variety of flow conditions. Br — Brinkman number - c specific heat, J/kg K - De — Dean number - E dimensionless apparent viscosity, eq. (14d) - G dimensionless shear rate, eq. (19) - k parameter of the power-law model, °C^–1, eq. (7) - mass flow rate, kg/s - m ₀ parameter of the power-law model, Pa · s ⁿ , eq. (7) - n parameter of the power-law model, eq. (7) - Nu 2r _p/ — Nusselt number, eqs. (28,31) - p pressure, Pa - Pe — Péclet number - P —(p/)/r _c — pressure gradient, Pa/m - dissipated energy, W, eq. (29) - total energy, W, eq. (30) - r radial coordinate, m - r _c radius of tube-curvature, m, fig. 1 - r _p radius of tube, m, fig. 1 - r _t variable, m, eq. (6) - R dimensionless radial coordinate, eq. (14a) - R _c dimensionlessr _c, eq. (14a) - R _t dimensionlessr _t, eq. (14a) - t temperature, °C - bulk temperature, °C, eq. (27) - t ₀ inlet temperature of the melt, °C - t _w tube wall temperature, °C - T dimensionless temperature, eq. (14c) - T _w dimensionless tube wall temperature - T dimensionless bulk temperature - u ₁ variable, s^–1, eq. (4) - u ₂ variable, s^–1, eq. (5) - U ₁ dimensionlessu ₁, eq. (18) - U ₂ dimensionlessu ₂, eq. (18) - v velocity in-direction, m/s - average velocity of the melt, m/s - V dimensionlessv, eq. (14b) - dimensionless , eq. (15) - z r _c — centre length of the tube, m - Z dimensionlessz, eq. (14e) - heat transfer coefficient, W/m² K - shear rate, s^–1, eq. (8) - — shear rate, s^–1 - apparent viscosity, Pa · s, eq. (7) - ₀ — apparent viscosity, Pa · s - angular coordinate, rad, fig. 1 - thermal conductivity, W/m K - melt density, kg/m³ - axial coordinate, rad, fig. 1 - rate of strain tensor, s^–1, eq. (8) - (—p) pressure drop, Pa     

Non-Linear response of a long,discontinuous fiber/melt system in elongational flows

The authors investigated the transient elongational behavior of a highly-aligned 600% volume fraction long, discontinuous fiber filled poly-ether-ketone-ketone melt with a computer-controlled extensional rheometer at 370°C. Prior experiments at controlled strain rate and stress produced _E ⁺ (t, ) and ^– (t, _E) similar to a shear dominated flow of a non-linear viscoelastic fluid. Stress relaxation following steady extension showed nonlinear effects in the change in stress decay rate with increasing strain rate. Continuous relaxation spectra showed a shift in the spectral peak to smaller values of with increasing strain rate. The Giesekus nonlinear constitutive relation modeled the elongation and stress relaxation with shearing rate at the fiber surface set by a strain rate magnification factor. Suitable for elongation, the model produced insufficient shift in the stress relaxation spectrum to account for the large change in stress decay rate exhibited in the experiments.English alphabet a _r aspect ratio of the fibers or l/d - A ₀ initial uniform cross-section area of the specimen - d fiber diameter - f fiber volume fraction - H() relaxation spectrum found by the method of Ferry and William l length of the fiber - L(t) time function specimen length - L ₀ initial specimen length - r radial coordinate across the shear cell - R _i fiber radius and inner cell dimension - R _o outer cell radius - t time in s - t _max duration of the extension - T _g glass transition temperature of the polymer - v velocity of the moving end of the test specimen - x axial position where is calculated Greek alphabet nonlinearity parameter in the Giesekus relation - axial mass distribution along the specimen major axis - shear strain rate - strain tensor - ₍₁₎ first convected derivative of the strain tensor - ₍₂₎ second convected derivative of the strain tensor - average strain at the end of extension as determined from - extension strain rate - average extension strain rate determined from - ^– transient strain rate under controlled stress, creep, test - _E elongational viscosity - _Eapp apparent elongational viscosity determined from - _E ⁺ transient elongational viscosity - ₀ zero shear rate viscosity - relaxation parameter - ₁ relaxation parameter in either Jeffrey's or Giesekus fluid - ₂ retardation parameter in either Jeffrey's or Giesekus fluid - _max relaxation value at which 99.9% of the H spectrum had occurred - _p relaxation value at which H reaches a maximum - volumetric composite density - _E elongational stress - _E ⁺ transient elongational stress - _E ^– controlled elongational stress, creep stress - _E ^y peak elongational stress in controlled experiment - shear stress at surface of the fiber in a shear cell - _yx simple shear component of the strain rate tensor - stress tensor - ₁ first convected derivative of the stress tensor     

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     

Flow of a simple non-newtonian fluid past a sphere

Summary Creeping flow past a sphere is solved for a limiting case of fluid behaviour: an abrupt change in viscosity.List of Symbols d _ij Component of rate-of-deformation tensor - F _d Drag force exerted on sphere by fluid - G ^(d) Coefficients in expression for _ij in terms of d _ij - G _YOJK ^(d) Coefficients in power series representing G ^(d) - R Radius of sphere - r Spherical coordinate - V Velocity of fluid very far from sphere - v _i Component of the velocity vector - x Dimensionless radial distance, r/R - x _i Rectangular Cartesian coordinate - Dimensionless quantity defined by (26) - ^(d) Potential defined by (7) - Value of x denoting border between Regions 1 and 2 as a function of - ₁, ₂ Lower and upper limiting viscosities defined by (10) - Spherical coordinate - * Value of for which =1 - Value of denoting border between regions 1 and 2 as a function of x - Newtonian viscosity - _ij Component of the stress tensor - Spherical coordinate - ₁, ₂ Stream functions defined by (12) and (14) - Second and third invariants of the stress tensor and of the rate-of-deformation tensor, defined by (3)     

Material functions generated by the complex viscosity

Based on the complex viscosity model various steady-state and transient material functions have been completed. The model is investigated in terms of a corotational frame reference. Also, BKZ-type integral constitutive equations have been studied. Some relations between material functions have been derived. C ^–1 Finger tensor - F[], (F ^–1[]) Fourier (inverse) transform - rate of deformation tensor in corotating frame - h(I, II) Wagner's damping function - J (x) Bessel function - m parameter inh (I, II) - m(s) memory function - m _k, n_k integers (powers in complex viscosity model) - P principal value of the integral - parameter in the complex viscosity model - rate of deformation tensor - shear rates - [], [] incomplete gamma function - (a) gamma function - steady-shear viscosity - ^* complex viscosity - , real and imaginary parts of ^* - ₀ zero shear viscosity - ⁺, ₁ ⁺ stress growth functions - ^–, ₁ ^- stress relaxation functions - (s) relaxation modulus - ₁(s) primary normal-stress coefficient - ø(a, b; z) degenerate hypergeometric function - ₁, ₂ time constants (parameters of ^*) - frequency - extra stress tensor     

Molecular weight distribution from rheological measurements

Summary A simple and reliable relative method to derive the molecular weight distribution of linear polymers is proposed.It is shown that both the zero-shear viscosity, ₀, and the intrinsic viscosity, [], have a logarithmic dependence on the weight average molecular weight, , and the polydispersity, . The coefficients of these relationships can be determined by applying a multiple regression analysis to a series of samples for which andQ are known.By making use of the two established relationships, the determination of andQ for a given polymer sample reduces to the experimental measurement of its ₀ and [].An analysis has been performed to estimate to what extent experimental errors on ₀ and [] affect the calculated molecular weight distribution.It has been found that only the experimental error on [] contributes heavily to the final error on the polydispersity.

---

Zusammenfassung Es wird eine einfache und zuverlässige Relativmethode vorgeschlagen, um die Uneinheitlichkeit linearer Polymere abzuleiten.Es wird gezeigt, daß alle beide, Nullschergradient-viskosität ₀, und Grenzviskositätszahl [], einfach logarithmisch vom Gewichtsmittel des Molekulargewichts , und vom Polymolekularitätsindex , abhängig sind.Die Koeffizienten dieser Beziehungen können mit statistischer Analyse festgesetzt werden, wenn undQ einer Probenreihe bekannt sind.Mit den zwei vorher festgesetzten Beziehungen besteht die Bestimmung von undQ einer gegebenen Polymersprobe nur aus den experimentellen Massen seiner ₀- und []-Werte.Eine Analyse wurde ausgeführt, um die Bedeutung des experimentellen Irrtums über die berechnete Uneinheitlichkeit zu wissen.Es wurde gefunden, daß ein experimenteller Irrtum betreffs [] schwer an endlichem Irrtum der Uneinheitlichkeit teilnimmt.

---

---

With 2 figures and 2 tables     

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     

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

Onset of viscous fingering for miscible liquid-liquid displacements in porous media

The displacement of one fluid by another miscible fluid in porous media is an important phenomenon that occurs in petroleum engineering, in groundwater movement, and in the chemical industry. This paper presents a recently developed stability criterion which applies to the most general miscible displacement. Under special conditions, different expressions for the onset of fingering given in the literature can be obtained from the universally applicable criterion. In particular, it is shown that the commonly used equation to predict the stable velocity ignores the effects of dispersion on viscous fingering.Nomenclature C Solvent concentration - Unperturbed solvent concentration - D _L Longitudinal dispersion coefficient [m²/s] - D _T Transverse dispersion coefficient [m²/s] - g Gravitational acceleration [m/s²] - I _sr Instability number - k Permeability [m²] - K Ratio of transverse to longitudinal dispersion coefficient - L Length of the porous medium [m] - L _x Width of the porous medium [m] - L _y Height of the porous medium [m] - M Mobility ratio - V Superficial velocity [m/s] - V _c Critical velocity [m/s] - V _s Velocity at the onset of instability [m/s] - µ Viscosity [Pa/s] - Unperturbed viscosity [Pa/s] - µ ₀,µ _s Viscosities of oil and solvent, respectively [Pa/s] - Density [kg/m³] - ₀, _s Densities of oil and solvent, respectively [kg/m³] - Porosity - Dimensionless length     

Production of temperature fluctuations in grid turbulence: Wiskind's experiment revisited

Measurements have been made in nearly-isotropic grid turbulence on which is superimposed a linearly-varying transverse temperature distribution. The mean-square temperature fluctuations, , increase indefinitely with streamwise distance, in accordance with theoretical predictions, and consistent with an excess of production over dissipation some 50% greater than values recorded in previous experiments. This high level of production has the effect of reducing the ratio,r, of the time scales of the fluctuating velocity and temperature fields. The results have been used to estimate the coefficient,C, in Monin's return-to-isotropy model for the slow part of the pressure terms in the temperature-flux equations. An empirical expression by Shih and Lumley is consistent with the results of earlier experiments in whichr 1.5, C 3.0, but not with the present data where r 0.5, C 1.6. Monin's model is improved when it incorporates both time scales.List of symbols C coefficient in Monin model, Eq. (5) - M grid mesh length - m exponent in power law for temperature variance, x ^m - n turbulence-energy decay exponent,q ² x ^-n - p production rate of - p pressure - q ² - R microscale Reynolds number - r time-scale ratiot/t - T mean temperature - U mean velocity - mean-square velocity fluctuations (turbulent energy components) - turbulent temperature flux - x, y, z spatial coordinates - temperature gradient dT/dy - thermal diffusivity - dissipation rate ofq ²/2 - dissipation rate of - Taylor microscale (₂=5q₂/) - temperature microscale - _v temperature-flux correlation coefficient, /v - dimensionless distance from the grid,x/M     

The measurement of the material parameters of viscoelastic fluids using a rotating sphere and a rheogoniometer

Summary In this work, measurement of the flow field around a rotating sphere has been used to obtain the material parameters of a second-order Rivlin-Ericksen fluid. Experiments were carried out with a Laser-Doppler anemometer to obtain the velocity distribution and usingGiesekus' analysis, the material parameters for the second-order fluid were obtained.

---

Zusammenfassung In dieser Untersuchung wird die Ausmessung des Strömungsfeldes um eine rotierende Kugel dazu verwendet, um die Stoffparameter einer Rivlin-Ericksen-Flüssigkeit zweiter Ordnung zu erhalten. Die Experimente zur Bestimmung der Geschwindigkeitsverteilung werden mit einem Laser-Doppler-Anemometer durchgeführt, und zur Auswertung der Parameter der Flüssigkeit zweiter Ordnung wird eine Analyse vonGiesekus benutzt.

---

Notations A ₁,A₂ Rivlin-Ericksen tensor - A ² Parameter used in eq. [12] - a Radius of the sphere - B Parameter used in eq. [12] - I Unit tensor - m ₀(₁ – ₂)/a², parameter used by ref. (8) - N ₁,N₂ First and second normal stress difference - p Isotropic pressure - Radial distance from the centre of the rotating body - S ₁,S₂ Stress tensor - v _r,v,v Velocity components in a spherical coordinate system - ₀,₁,₂ Material parameters used in eq. [2] - Shear rate - _a Apparent voscosity - ₀ Zero-shear viscosity - Angle measured from the axis of rotation - Fluid density - Stream function - Shear stress - Angular velocity With 3 figures