A technique for evaluation of three-dimensional behavior in turbulent boundary layers using computer augmented hydrogen bubble-wire flow visualization

A system is described which allows the recreation of the three-dimensional motion and deformation of a single hydrogen bubble time-line in time and space. By digitally interfacing dualview video sequences of a bubble time-line with a computer-aided display system, the Lagrangian motion of the bubble-line can be displayed in any viewing perspective desired. The u and v velocity history of the bubble-line can be rapidly established and displayed for any spanwise location on the recreated pattern. The application of the system to the study of turbulent boundary layer structure in the near-wall region is demonstrated.List of Symbols Reynolds number based on momentum thickness u /v - t⁺ nondimensional time - u shear velocity - u local streamwise velocity, x-direction - u ⁺ nondimensional streamwise velocity - v local normal velocity, -direction - x ⁺ nondimensional coordinate in streamwise direction - ⁺ nondimensional coordinate normal to wall - ₊ ^wire wire nondimensional location of hydrogen bubble-wire normal to wall - z ⁺ nondimensional spanwise coordinate - momentum thickness - v kinematic viscosity - _W wall shear stress     

Momentum balance in highly distorted turbulent pipe flows

An in depth study into the development and decay of distorted turbulent pipe flows in incompressible flow has yielded a vast quantity of experimental data covering a wide range of initial conditions. Sufficient detail on the development of both mean flow and turbulence structure in these flows has been obtained to allow an implied radial static pressure distribution to be calculated. The static pressure distributions determined compare well both qualitatively and quantitatively with earlier experimental work. A strong correlation between static pressure coefficient Cp and axial turbulence intensity is demonstrated.List of symbols C _p static pressure coefficient = (pw-p)/1/2 - D pipe diameter - K turbulent kinetic energy - (r, , z) cylindrical polar co-ordinates. / 0 - R, y pipe radius, distance measured from the pipe wall - U, V axial and radial time mean velocity components - mean value of u - u, u/ , / - u, , w fluctuating velocity components - axial, radial turbulence intensity - turbulent shear stress - u friction velocity, (u ² = ₀/p) - ₀ wall shear stress - ^* boundary layer thickness A version of this paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984     

Mathematical model of the entry length in the flow of suspensions of solid particles in a gas

A mathematical model consisting of equations of mass and momentum and for the velocity field has been used for computing the entry length of the flow of non-Newtonian fluids in laminar, transition and turbulent regions. Experimental data measured in a vertical flow of a suspension of solid particles in air have been used for verifying the predictions. n flow index for laminar flow - Re Reynolds number defined for the flow of the carrier medium - q exponent for turbulent flow - ratio of core radius with a flat velocity profile to pipe radius - _c ratio of the axial component of local velocity in the core to mean velocity - w mean flow velocity - ratio of axial distance from the pipe entrance to the pipe radius - ratio of the entrance length to the pipe radius - relative mass fraction of particles - ratio of the distance from the pipe wall to the pipe radius - coefficient of pressure loss due to friction     

LDV measurements of a turbulent air-solid two-phase flow in a 90° bend

Simultaneous measurements of the mean streamwise and radial velocities and the associated Reynolds stresses were made in an air-solid two-phase flow in a square sectioned (10×10 cm) 90° vertical to horizontal bend using laser Doppler velocimetry. The gas phase measurements were performed in the absence of solid particles. The radius ratio of the bend was 1.76. The results are presented for two different Reynolds numbers, 2.2×10⁵ and 3.47×10⁵, corresponding to mass ratios of 1.5×10^–4 and 9.5×10^–5, respectively. Glass spheres 50 and 100 m in diameter were employed to represent the solid phase. The measurements of the gas and solid phase were performed separately. The streamwise velocity profiles for the gas and the solids crossed over near the outer wall with the solids having the higher speed near the wall. The solid velocity profiles were quite flat. Higher negative slip velocities are observed for the 100 m particles than those for the 50 gm particles. At angular displacement =0°, the radial velocity is directed towards the inner wall for both the 50 and 100 m particles. At =30° and 45°, particle wall collisions cause a clear change in the radial velocity of the solids in the region close to the outer wall. The 100 m particle trajectories are very close to being straight lines. Most of the particle wall collisions occur between the =30° and 60° stations. The level of turbulence of the solids was higher than that of the air.List of symbols D hydraulic diameter (100 mm) - De Dean number,De = - mass flow rate - number of particles per second (detected by the probe volume) - r radial coordinate direction - r _i radius of curvature of the inner wall - r ₀ radius of curvature of the outer wall - r ^* normalized radial coordinate, - R mean radius of curvature - Re Reynolds number, - R _r radius ratio, - U ,U _z mean streamwise velocity - U _r ,U _y mean radial velocity - U _b bulk velocity - , _z rms fluctuating streamwise velocity - _r , _y rms fluctuating radial velocity - -r shear stress component - z-y shear stress component - x spanwise coordinate direction - x ^* normalized spanwise coordinate, - y radial coordinate direction in straight ducts - y ^* normalized radial coordinate in straight ducts, - z streamwise coordinate direction in straight ducts - z ^* normalized streamwise coordinate in straight ducts, Greek symbols streamwise coordinate direction - kinematic viscosity of air     

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 flow regime in the turbulent near wake of a hemisphere

The work presented is a wind tunnel study of the near wake region behind a hemisphere immersed in three different turbulent boundary layers. In particular, the effect of different boundary layer profiles on the generation and distribution of near wake vorticity and on the mean recirculation region is examined. Visualization of the flow around a hemisphere has been undertaken, using models in a water channel, in order to obtain qualitative information concerning the wake structure.List of symbols C _p pressure coefficient, - D diameter of hemisphere - n vortex shedding frequency - p pressure on model surface - p ₀ static pressure - Re Reynolds number, - St Strouhal number, - U, V, W local mean velocity components - mean freestream velocity inX direction - U ^* shear velocity, - u, v, w velocity fluctuations inX, Y andZ directions - X Cartesian coordinate in longitudinal direction - Y Cartesian coordinate in lateral direction - Z Cartesian coordinate in direction perpendicular to the wall - it* boundary layer displacement thickness, - diameter of model surface roughness - elevation angleI - O boundary layer momentum thickness, - _w wall shearing stress - dynamic viscosity of fluid - density of fluid - streamfunction - _x longitudinal component of vorticity, - _y lateral component of vorticity, - _z vertical component of vorticity, This paper was presented at the Ninth symposium on turbulence, University of Missouri-Rolla, October 1–3, 1984     

Laminar source flow between parallel porous disks

An analysis is presented for laminar source flow between infinite parallel porous disks. The solution is in the form of a perturbation from the creeping flow solution. Expressions for the velocity, pressure, and shear stress are obtained and compared with the results based on the assumption of creeping flow.Nomenclature a half distance between disks - radial coordinate - r dimensionless radial coordinate, /a - axial coordinate - z dimensionless axial coordinate, /a - radial coordinate of a point in the flow - R dimensionless radial coordinate of a point in the flow, /a - velocity component in radial direction - u =a/, dimensionless velocity component in radial direction - velocity component in axial direction - v = a/}, dimensionless velocity component in axial direction - static pressure - p = (a ²/ ²), dimensionless static pressure - =p(r, z)–p(R, z), dimensionless pressure drop - V magnitude of suction or injection velocity - Q volumetric flow rate of the source - Re source Reynolds number, Q/4a - reduced Reynolds number, Re/r ² - critical Reynolds number - R _w wall Reynolds number, Va/ - viscosity - density - =/, kinematic viscosity - shear stress at upper disk - ₀ = (a ²/ ²), dimensionless shear stress at upper disk - shear stress ratio, ₀/( ₀)_inertialess - u = , dimensionless average radial velocity - u/u, ratio of radial velocity to average radial velocity - dimensionless stream function     

Velocity field in isothermal turbulent bubbly gas-liquid flow through a pipe

Velocity field was measured by laser Doppler velocimetry in isothermal, turbulent bubbly gas-liquid flow through a 26.6 mm inner diameter vertical pipe. The measurements were made about 33 diameters downstream from the pipe entrance, gas injection being just upstream of the entrance. The gas phase radial distribution at the measurement plane exhibited influence of the injection device in that higher gas fraction existed in the central region of the pipe. For comparison, velocity field was also measured in isothermal, turbulent single-phase liquid flow through the same pipe at the same axial plane. Measured were the radial distributions of liquid mean axial and radial velocities, axial and radial turbulent intensities, and axial Reynolds shear stress. The radial distributions of gas bubble mean axial velocity and axial velocity fluctuation intensity were also measured by LDV. A dualsensor fiberoptic probe was used at the same time to measure the radial distributions of gas fraction, bubble mean axial velocity and size slightly downstream of the LDV measurement plane.List of Symbols an average gas bubble diameter - f, f _TP friction factor, friction factor for gas-liquid flow - k _L liquid turbulent kinetic energy - , gas, liquid mass flow rate - R inner radius of pipe - r, {sitR}^* radial coordinate; nondimensional radial coordinate (=r/R) - Re _L liquid Reynolds number - U _G mean axial velocity of gas bubble - U _L mean axial velocity of liquid - U _LO mean axial velocity for flow at the total mass velocity with properties of the liquid phase - u _L ⁺ nondimensional mean axial velocity of liquid in wall coordinate - friction velocity - axial velocity fluctuation intensity of liquid - axial velocity fluctuation intensity of gas bubbles - V_L mean radial velocity of liquid - v _L radial velocity fluctuation intensity of liquid - (uv)_L single-point cross-correlation between axial and radial velocity fluctuations of liquid ( axial Reynolds shear stress) - T _in mean liquid temperature at test section inlet - x flow quality - y normal distance from wall - y ⁺ nondimensional normal distance from wall in wall coordinate (=yu/v_L) - _G gas phase residence time fraction - _L rate of dissipation in the liquid - _L Kolmogorov length scale in the liquid - _L liquid kinematic viscosity - _L characteristic turbulence length scale in the liquid - _G, _L density of gas, liquid - _m gas-liquid mixture density This work was partly supported by National Science Foundation, Thermal Transport and Thermal Processing Program, Chemical and Thermal Systems Division, under Grant No. CTS-9411898.     

Optimal cupolas of uniform strength: Spherical M-shells and axisymmetric T-shells

Summary The first part of this paper is concerned with the optimal design of spherical cupolas obeying the von Mises yield condition. Five different load combinations, which all include selfweight, are investigated. The second part of the paper deals with the optimal quadratic meridional shape of cupolas obeying the Tresca yield condition, considering selfweight plus the weight of a non-carrying uniform cover. It is established that at long spans some non-spherical Tresca cupolas are much more economical than spherical ones.

---

Optimale Kuppeln gleicher Festigkeit: Kugelschalen und axialsymmetrische Schalen

Übersicht Im ersten Teil dieser Arbeit wird der optimale Entwurf sphärischer Kuppeln behandelt, wobei die von Misessche Fließbewegung zugrunde gelegt wird. Fünf verschiedene Lastkombinationen werden untersucht. Der zweite Teil befaßt sich mit der optimalen quadratischen Form des Meridians von Kuppeln, die der Fließbedingung von Tresca folgen.

---

List of Symbols a_k, b_k, c_k, A_k, B_k, C_k coefficients used in series solutions - A, B constants in the nondimensional equation of the meridional curve - normal component of the load per unit area of the middle surface - meridional and circumferential forces per unit width - radial pressure per unit area of the middle surface, - skin weight per unit area of the middle surface, - vertical external load per unit horizontal area, - base radius, - R radius of convergence - s - cupola thickness, - u, w subsidiary functions for quadratic cupolas - vertical component of the load per unit area of middle surface - resultant vertical force on a cupola segment - structural weight of cupola, - combined weight of cupola and skin, - distance from the axis of rotation, - vertical distance from the shell apex, - z auxiliary variable in series solutions - specific weight of structural material of cupola - radius of the middle surface, - uniaxial yield stress - meridional stress, - circumferential stress, - _a, _b, _c, _d, _e subsidiary variables used in evaluating the meridional stress - auxiliary function used in series solutions This paper constitutes the third part of a study of shell optimization which was initiated and planned by the late Prof. W. Prager     

Evaluation of the performance characteristics of a thermal transient anemometer

The Goertler instability of a hypersonic boundary layer and its influence on the wall heat transfer are experimentally analyzed. Measurements, made in a wind tunnel by means of a computerized infrared (IR) imaging system, refer to the flow over two-dimensional concave walls. Wall temperature maps (that are interpreted as surface flow visualizations) and spanwise heat transfer fluctuations are presented. Measured vortices wavelengths are correlated to non-dimensional parameters and compared with numerical predictions from the literature.List of symbols c _p Specific heat coefficient at constant pressure of the free stream - F Input (true) image - F ₀ Fourier number - Restored image - G Recorded (degraded) image - G Goertler number based on the boundary layer thickness, as defined by Eq. (3) - H System transfer function - M Mach number - Pr Prandtl number - p ₀ Stagnation pressure - Exchanged net heat flux - Convective heat flux - Radiative heat flux - r Recovery factor - Re _m Unit Reynolds number - Re _x Local Reynolds number based on the distance from the leading edge - Re Local Reynolds number based on the boundary layer thickness - Curvature radius - St Stanton number, as defined by Eq. (7) - T _aw Adiabatic wall temperature - T _w Wall temperature - T ₀ Stagnation temperature - t Time - V Free stream velocity - x Streamwise spatial coordinate - y Normal-to-wall spatial coordinate - z Spanwise spatial coordinate - Thermal diffusivity coefficient - Disturbance wavenumber - Non dimensional wavenumber - Boundary layer thickness - Goertler number based on the vortices wavelength - Vortices wavelength - Free stream density - Disturbance total amplification, as defined by Eq. (3) - Disturbance (spatial) growth rate - Non-dimensional growth rate - Perturbation amplitude of a generic quantity - Perturbation amount     

Apparent thickening behavior of dilute polystyrene solutions in extensional flows

An experimental investigation was undertaken to study the apparent thickening behavior of dilute polystyrene solutions in extensional flow. Among the parameters investigated were molecular weight, molecular weight distribution, concentration, thermodynamic solvent quality, and solvent viscosity. Apparent relative viscosity was measured as a function of wall shear rate for solutions flowing from a reservoir through a 0.1 mm I.D. tube. As increased, slight shear thinning behavior was observed up until a critical wall shear rate was exceeded, whereupon either a large increase in or small-scale thickening was observed depending on the particular solution under study. As molecular weight or concentration increased, decreased and, the jump in above , increased. increased as thermodynamic solvent quality improved. These results have been interpreted in terms of the polymer chains undergoing a coil-stretch transition at . The observation of a drop-off in at high (above ) was shown to be associated with inertial effects and not with chain fracture due to high extensional rates.     

Effects of couple-stresses on the dispersion of a soluble matter in a pipe flow of blood

Summary An analysis of the effects of couple-stresses on the effective Taylor diffusion coefficient has been carried out with the help of two non-dimensional parameters based on the concentration of suspensions and , a constant associated with the couple-stresses. It is observed that the concentration distribution increases with increasing or The effective Taylor diffusion coefficient falls rapidly with increasing when is negative.

---

Zusammenfassung Der Einfluß der Momentenspannungen auf den effektiven Taylorschen Diffusionskoeffizienten wird untersucht. Dabei treten zwei dimensionslose Parameter and auf: Der erste bezieht sich auf die Suspensionskonzentration, der zweite kennzeichnet die Momentenspannungen. Man findet, daß die Verteilungsgeschwindigkeit mit wachsendem oder zunimmt. Dagegen fällt der Taylorsche Diffusionskoeffizient bei wachsendem stark ab, wenn negativ ist.

---

a Tube radius - C Concentration - C _i Body moment vector - C ₀ Concentration at the axis of the tube - C _m Mean concentration - D Molecular diffusion coefficient - d _ij Symmetric part of velocity gradient - F Function of and characterising effective Taylor diffusion coefficient - f _i Body force vector - H A function of and - K ₂ Integration constant - K ^* Effective Taylor diffusion coefficient - k Radius of gyration of a unit cuboid with its sides normal to the spatial axes - I _n Modified Bessel's function ofnth order - L Length of the tube over which the concentration is spread - M Function ofH and - M _ij Couple stress tensor - P Function of - p Fluid pressure - Q Volume rate of the transport of the solute across a section of the tube - r Radial distance from the axis of the tube - T _ij Stress tensor - t Time coordinate - T _ij ^A Antisymmetric part of the stress tensor - u Relative fluid velocity - Average velocity - v _i Velocity vector - Fluid velocity at any point of the tube - v ₀ ⁿ Velocity of Newtonian fluid at the axis of the tube - _i Vorticity vector - x Axial coordinate - x ₁ Relative axial coordinate - z Non-Dimensional radial coordinate - Density - _ij Symmetric part of the stress tensor - µ Viscosity of the fluid - µ _ij Deviatoric part ofM _ij - , Constants associated with couple-stress With 3 figures     

On the continuum modelling of porous media containing fluid: a molecular viewpoint with particular attention to scale

Mass conservation and linear momentum balance relations for a porous body and any fluid therein, valid at any given length scale in excess of nearest-neighbour molecular separations, are established in terms of local weighted averages of molecular quantities. The mass density field for the porous body at a given scale is used to identify its boundary at this scale, and a porosity field is defined for any pair of distinct length scales. Specific care is paid to the interpretation of the stress tensor associated with each of the body and fluid at macroscopic scales, and of the force per unit volume each exerts on the other. Consequences for the usual microscopic and macroscopic viewpoints are explored.Nomenclature material system; Section 2.1. - porous body (example of a material system); Sections 2.1, 3.1, 4.1 - fluid body (example of a material system); Sections 2.1, 3.1, 4.1 - weighting function; Sections 2.1, 2.3 - ,h weighting function corresponding to spherical averaging regions of radius and boundary mollifying layer of thicknessh; Section 3.2 - Euclidean space; Section 2.1 - V space of all displacements between pairs of points in; Section 2.1 - mass density field corresponding to; (2.3)₁ - ^P , ^f mass density fields for , ; (4.1) - P momentum density field corresponding to; (2.3)₂ - v velocity field corresponding to; (2.4) - S _r (X) interior of sphere of radiusr with centre at pointx; (3.3) - boundary ofany region - region in which ^p > 0 with = ,h; (3.1) - subset of whose points lie at least+h from boundary of ; (3.4) - abbreviated versions of ; Section 3.2, Remark 4 - strict interior of ; (3.7) - analogues of for fluid system ; Section 3.2 - general version of corresponding to any choice of weighting function; (4.6) - interfacial region at scale; (3.8) - ₀ scale of nearest-neighbour separations in ; Section 3.2. Remark 1 - porosity field at scales ( ₁; ₂); (3.9) - pore space at scales ( ₁; ₂); (3.12)     

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     

Das Dispersionsmodell

Zusammenfassung Ein Vergleich im Frequenzbereich zeigt, daß bei der Berechnung der Verweilzeitverteilung mit dem Dispersionsmodell das endlich begrenzte System für Péclet-Zahlen Pe > 10 mit guter Näherung durch ein einseitig unbegrenztes System und für Pe > 50 durch ein beidseitig unbegrenztes System ersetzt werden kann.

---

The dispersion model. A comparison of approximations

A comparison in the frequency domain shows that for the determination of the residence time distribution with the dispersion model the finitely restricted system may be substituted with good approximation for Peclet numbers Pe > 10 by a one-side unrestricted system and for Pe > 50 by a both-side unrestricted system.

---

Bezeichnungen A Rohrquerschnitt - A=A mittlerer Strömungsquerschnitt in der Schüttschicht - Konzentration (Partialdichte) der Bezugskomponente i - Bezugskonzentration nach Gl. (5) - c_i Konzentration (Dichte) der reinen Bezugskomponente i - D axialer Dispersionskoeffizient - E Fehler im Frequenzbereich nach Gl. (36) - G(,) Übertragungsfunktion - G(,i) Frequenzgang - L Länge der Schüttschicht - m Masse - Massenstrom - Péclet-Zahl - s Laplace-Variable - t Zeit - t Impulsbreite - v Strömungsgeschwindigkeit im leeren Rohr - mittlere axiale Strömungsgeschwin digkeit in der Schüttschicht - V=AL Zwischenraumvolumen der Schüttschicht - x Ortskoordinate - (t) Dirac-Stoss - Porosität - dimensionslose Zeit - dimensionslose Konzentration - Laplace-Transformierte der Konzentration - Fourier-Transformierte der Konzentration - dimensionslose Ortskoordinate - =s dimensionslose Laplace-Variable - mittlere Verweilzeit - Kreisfrequenz - = dimensionslose Kreisfrequenz Indices A Ausgang - D Dispersion - E Eingang - i Bezugskomponente - K Konvektion Mitteilung Nr. 44 des Institutes für Mess-und Regel-technik der Eidgenössischen Technischen Hochschule Zürich (Vorsteher: Prof. Dr. P. Profos)     

Dimensional analysis of pore scale and field scale immiscible displacement

A basic re-examination of the traditional dimensional analysis of microscopic and macroscopic multiphase flow equations in porous media is presented. We introduce a macroscopic capillary number which differs from the usual microscopic capillary number Ca in that it depends on length scale, type of porous medium and saturation history. The macroscopic capillary number is defined as the ratio between the macroscopic viscous pressure drop and the macroscopic capillary pressure. can be related to the microscopic capillary number Ca and the LeverettJ-function. Previous dimensional analyses contain a tacit assumption which amounts to setting = 1. This fact has impeded quantitative upscaling in the past. Our definition for , however, allows for the first time a consistent comparison between macroscopic flow experiments on different length scales. Illustrative sample calculations are presented which show that the breakpoint in capillary desaturation curves for different porous media appears to occur at 1. The length scale related difference between the macroscopic capillary number for core floods and reservoir floods provides a possible explanation for the systematic difference between residual oil saturations measured in field floods as compared to laboratory experiment.     

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.     

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     

Unsteady flow in an annulus between two concentric rotating spheres

An analysis is presented for the unsteady laminar flow of an incompressible Newtonian fluid in an annulus between two concentric spheres rotating about a common axis of symmetry. A solution of the Navier-Stokes equations is obtained by employing an iterative technique. The solution is valid for small values of Reynolds numbers and acceleration parameters of the spheres. In applying the results of this analysis to a rotationally accelerating sphere, a virtual moment of intertia is introduced to account for the local inertia of the fluid.Nomenclature R _i radius of the inner sphere - R _o radius of the outer sphere - radial coordinate - r dimensionless radial coordinate, - meridional coordinate - azimuthal coordinate - time - t dimensionless time, - Re _i instantaneous Reynolds number of the inner sphere, _i R _k ² / - Re _o instantaneous Reynolds number of the outer sphere, _o R _o ² / - radial velocity component - V _r dimensionless radial velocity component, - meridional velocity component - V dimensionless meridional velocity component, - azimuthal velocity component - V dimensionless azimuthal velocity component, - viscous torque - T dimensionless viscous torque, - viscous torque at surface of inner sphere - T _i dimensionless viscous torque at surface of inner sphere, - viscous torque at surface of outer sphere - T _o dimensionless viscous torque at surface of outer sphere, - externally applied torque on inner sphere - T _p,i dimensionless applied torque on inner sphere, - moment of inertia of inner sphere - Z _i dimensionless moment of inertia of inner sphere, - virtual moment of inertia of inner sphere - Z _i,v dimensionless virtual moment of inertia of inner sphere, - virtual moment of inertia of outer sphere - _i instantaneous angular velocity of the inner sphere - _o instantaneous angular velocity of the outer sphere - density of fluid - viscosity of fluid - kinematic viscosity of fluid,/ - radius ratio,R _i/R _o - swirl function, - dimensionless swirl function, - stream function - dimensionless stream function, - _i acceleration parameter for the inner sphere, - _o acceleration parameter for the outer sphere, - shear stress - _r dimensionless shear stress,      

Yield stress: A time-dependent property and how to measure it

This paper reviews the different aspects of the yield stress phenomenon and attempts a synthesis of knowledge. Yield stress can be probed using constant shear stress or shear rate. The magnitude of the result depends on the time allowed to determine whether the sample has developed continuous flow or has ceased flowing. It is closely associated with creep, stress growth and thixotropic breakdown and recovery, and the characteristic times of these transient responses play a part in yield stress measurement. In thixotropic fluids, yield stress is a function of structure and hence of time. In simple thixotropy, the yield stress derived from the equilibrium flow curve is the same as that for the fully built-up structure. But in many materials, the static yield stress obtained after prolonged rest is much higher than the dynamic yield stress from the equilibrium flow curve. This is associated with the phenomenon in which the equilibrium flow curve bends upwards as the shear rate is reduced to very low values. The paper also reviews the many methods that can be used to measure yield stress. It is pointed out that the choice of observation time or shear rate to use should be related to the characteristic time of the flow process to which the result is to be applied. Examples discussed are solids suspension capability of fluids, levelling and sagging, pipeline flow and start-up power requirement of mixers. CS constant structure - D diameter of Gun Rheometer tube - EFC equilibrium flow curve - G measured torque - L length of Gun Rheometer tube - P _min minimum pressure to cause flow - t time - form factor for shear stress - - _y - shear rate - a particular value of shear rate - reference shear rate - test shear rate - shear stress - _y yield stress - _yd dynamic yield stress - _ys static yield stress - 0 initial value after speed change - e equilibrium Paper presented at the British Society of Rheology Conference on New Techniques in Experimental Rheology, University of Reading, 9–11 September 1985.