Pore structures and transport properties of sandstone

To quantitatively analyze the macroscopic properties of the flow in porous media by means of the continuum approach, detailed information (velocity and pressure fields) on the microscopic scale is necessary. In this paper, the numerical solution for incompressible, Newtonian flow in a diverging-converging representative unit cell (RUC) is presented. A new solution procedure for the problem is introduced. A review of the accuracy of the computational method is given.Nomenclature A _ff ^* area of entrance and exit of RUC - A _fs ^* interfacial area between the fluid and solid phases - d throat diameter of RUC (m) - D pore diameter of RUC (m) - i, j unit vector for RUC - L ^* wave length of a unit cell - L _p pore length of RUC (m) - L _t throat length of RUC (m) - n unit outwardly directed vector for the fluid phase - p ^* fluid pressure - ^* cross-sectional mean pressure - _en ^* entrance cross-sectional mean pressure - Re _d Reynolds number - x ^*, r^* cylindrical coordinates - u ^*, v^* velocity - u _cl ^* centerline velocity - _d mean velocity at the throat of RUC (m/s) - _D mean velocity at the large segment of RUC (m/s) Greek viscosity coefficient (Ns/m²) - _p excess momentum loss factor defined in (4.1) - fluid density (kg/m³) - ^* stream function - ^* vorticity - dimensionless circulation defined in (2.7) Symbols - the mean value - * dimensionless quantities     

Hydro-mechanical aspects of the sand production problem

This paper examines the hydro-mechanical aspect of the sand production problem and sets the basic frame of the corresponding mathematical modelling. Accordingly, piping and surface erosion effects are studied on the basis of mass balance and particle transport considerations as well as Darcy's law. The results show that surface erosion is accompanied by high changes of porosity and permeability close to the free surface. Quantities which can be measured in experiment, like the amount of produced solids or fluid discharge, can be used in an inverse way to determine the constitutive parameters of the problem.Notation dV Volume element - dV _s Volume of solids pt - dV _v Volume of voids - dV _ff Volume of fluid phase - dV _fs Volume of fluidized-particles - Volume of mixture - dM _s Mass of solids - dM _ff Mass of fluid phase - d _{M fs} Mass of fluidized-particles - Mass of mixture - s Density of solids - f Density of fluid - ff Density of fluid phase - fs Density of fluidized-particles - Density of mixture - _i ^ff Velocity of fluid - _i ^fs Velocity of fluidized-particles - _i ^s Velocity of solids - Velocity of mixture - q ^ff Volume-discharge of fluid - q ^fs Volume-discharge of fluidized-particles - Volume-discharge of mixture - m ^ff Mass-discharge of fluid - m ^fs Mass-discharge of fluidized-particles - Mass-discharge of mixture - _er Rate of mass-eroded - _dep Rate of mass-deposited - Mass generation term - dS _i Surface element - Pore-surface element - D _IJ Tensor of mechanical dispersion - x _i Location - t Time - Porosity - c Transport concentration - c _cr Critical value ofc - p Fluid-pressure - k Permeability coefficient - _k Kinematic viscosity - Spatial frequency of erosion starter points     

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     

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     

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     

Transient flow towards a well in an aquifer including the effect of fluid inertia

Transient propagation of weak pressure perturbations in a homogeneous, isotropic, fluid saturated aquifer has been studied. A damped wave equation for the pressure in the aquifer is derived using the macroscopic, volume averaged, mass conservation and momentum equations. The equation is applied to the case of a well in a closed aquifer and analytical solutions are obtained to two different flow cases. It is shown that the radius of influence propagates with a finite velocity. The results show that the effect of fluid inertia could be of importance where transient flow in porous media is studied.List of symbols b Thickness of the aquifer, m - c ₀ Wave velocity, m/s - k Permeability of the porous medium, m² - n Porosity of the porous medium - p( ,t) Pressure, N/m² - Q Volume flux, m³/s - r Radial coordinate, m - r _w Radius of the well, m - s Transform variable - S Storativity of the aquifer - S _d(r, t) Drawdown, m - t Time, s - T Transmissivity of the aquifer, m²/s - ( ,t) Velocity of the fluid, m/s - Coordinate vector, m - z Vertical coordinate, m - Coefficient of compressibility, m²/N - Coefficient of fluid compressibility, m²/N - Relaxation time, s - (r, t) Hydraulic potential, m - Dynamic viscosity of the fluid, Ns/m² - Dimensionless radius - Density of the fluid, Ns²/m⁴ - (, ) Dimensionless drawdown - Dimensionless time - , x Dummy variables - ₀, ₁ Auxilary functions     

Thermal effects in the capillary rheometer

Summary The effect of viscous heating in a capillary rheometer is analysed for a power-law fluid by means of a perturbation expansion based upon a boundary-layer-core structure. This expansion is found to complement the eigenfunction series solution obtained by earlier investigators. A similar analysis is presented for the work-of-expansion effect. These two thermal effects are superimposed together with a third perturbation effect due to the pressure dependence of viscosity.On the basis of the present theory, earlier work in this area is discussed and, in some cases, apparent inaccuracies or inconsistencies are pointed out. A means is indicated for correcting data on the basis of the present theory.

---

Zusammenfassung Es wird der Effekt der Erwärmung einer Potenzflüssigkeit infolge viskoser Reibung in einem Kapillar-Rheometer mittels einer Störungsrechnung untersucht, die auf der Unterteilung der Strömung in eine Grenzschicht und einen Kern basiert. Diese Störungsentwicklung ergänzt eine früher von anderen Autoren gefundene Reihenentwicklung mit Hilfe von Eigenfunktionen. Eine ähnliche Untersuchung wird für die thermische Ausdehnungsarbeit durchgeführt. Diese beiden thermischen Effekte sind zusammen einem dritten Störeffekt superponiert, der von der Druckabhängigkeit der Viskosität herrührt.Aufgrund der vorgelegten Theorie werden verschiedene auf diesem Gebiet früher durchgeführte Arbeiten diskutiert, und es werden in einigen Fällen offensichtliche Ungenauigkeiten und Folgewidrigkeiten aufgedeckt. Schließlich wird eine Methode zur Korrektur von Meßdaten mit Hilfe der vorliegenden Theorie angegeben.

---

Nomenclature a tube radius - b ; evaluated atT ₀ andp = 0 when used in perturbation expansion - C _p specific heat - f - f ^* - h defined by eq. [15] - k thermal conductivity - L tube length - m defined by eq. [8] - m ₀ m(T₀, 0) - n power-law index - p pressure - Pe C _p W a/k Peclet number - Pr C _pa/k Prandtl number - Q volumetric flow rate - Q ₀ unperturbed value ofQ in specified-p formulation - r radial coordinate - Re W a/ _a Reynolds number - T temperature - T ₀ inlet temperature - u radial velocity component - u ₀ 0 unperturbed radial velocity - w axial velocity component - w ₀ /W(1 – ) unperturbed axial velocity - W Q/(a ²) average axial velocity - W ₀ Q ₀/(a ²) - z axial coordinate - (3n + 1)/n - ^* ; evaluated atT ₀ andp = 0 when used in perturbation expansion - 4^1-n - ^* - (n + 1)/n - ... shear rate - 4W/a apparent shear rate - p total pressure drop - T _a W ²/k characteristic temperature difference - T _b total bulk-temperature rise - ^* T - r/a - shear viscosity - _a m₀ - (1 –)/ ^1/3 - —p/z - ₀ ... unperturbed value of - z-averaged value of - µ ^{n + 1}/n - z/(a Pe) - _L L/(a Pe) - mass density - _w shear stress at wall - streamfunction - ^*T₀ (absolute temperature scale) - ( )₁ leading-order effect due to viscous heating - ( ) ₁ ^* leading-order effect due to work-of-expansion Note: in specified-p formulation,W gets replaced byW ₀ in definition of Pe, Re, and. With 7 figures and 7 tables     

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     

A three-parameter model describing the experimental relation between shear stress and shear rate for laminar flow of aqueous polymer solutions

Summary A three-parameter model is introduced to describe the shear rate — shear stress relation for dilute aqueous solutions of polyacrylamide (Separan AP-30) or polyethylenoxide (Polyox WSR-301) in the concentration range 50 wppm – 10,000 wppm. Solutions of both polymers show for a similar rheological behaviour. This behaviour can be described by an equation having three parameters i.e. zero-shear viscosity ₀, infinite-shear viscosity , and yield stress ₀, each depending on the polymer concentration. A good agreement is found between the values calculated with this three-parameter model and the experimental results obtained with a cone-and-plate rheogoniometer and those determined with a capillary-tube rheometer.

---

Zusammenfassung Der Zusammenhang zwischen Schubspannung und Schergeschwindigkeit von strukturviskosen Flüssigkeiten wird durch ein Modell mit drei Parametern beschrieben. Mit verdünnten wäßrigen Polyacrylamid-(Separan AP-30) sowie Polyäthylenoxidlösungen (Polyox WSR-301) wird das Modell experimentell geprüft. Beide Polymerlösungen zeigen im untersuchten Schergeschwindigkeitsbereich von ein ähnliches rheologisches Verhalten. Dieses Verhalten kann mit drei konzentrationsabhängigen Größen, nämlich einer Null-Viskosität ₀, einer Grenz-Viskosität und einer Fließgrenze ₀ beschrieben werden. Die Ergebnisse von Experimenten mit einem Kegel-Platte-Rheogoniometer sowie einem Kapillarviskosimeter sind in guter Übereinstimmung mit den Werten, die mit dem Drei-Parameter-Modell berechnet worden sind.

---

a Pa^–1 physical quantity defined by:a = {1 – ( / ₀)}/ ₀ - c l concentration (wppm) - D m capillary diameter - L m length of capillary tube - P Pa pressure drop - R m radius of capillary tube - u m s^–1 average velocity - v _r m s^–1 local axial velocity at a distancer from the axis of the tube - shear rate (–dv _r /dr) - local shear rate in capillary flow - s^–1 wall shear rate in capillary flow - Pa s dynamic viscosity - _a Pa s apparent viscosity defined by eq. [2] - ( _a ) Pa s apparent viscosity in capillary tube at a distanceR from the axis - ₀ Pa s zero-shear viscosity defined by eq. [4] - Pa s infinite-shear viscosity defined by eq. [5] - l ratior/R - kg m^– density - Pa shear stress - ₀ Pa yield stress - _r Pa local shear stress in capillary flow - _R Pa wall shear stress in capillary flow _R = (PR/2L) - _v m³ s^–1 volume rate of flow With 8 figures and 1 table     

Parameter determination for thermodynamical models of the pseudoelastic behaviour of shape memory alloy by resistivity measurements and infrared thermography

Two thermodynamical models of pseudoelastic behaviour of shape memory alloys have been formulated. The first corresponds to the ideal reversible case. The second takes into account the hysteresis loop characteristic of this shape memory alloys.Two totally independent techniques are used during a loading-unloading tensile test to determine the whole set of model parameters, namely resistivity and infrared thermography measurements. In the ideal case, there is no difficulty in identifying parameters.Infrared thermography measurements are well adapted for observing the phase transformation thermal effects.Notations 1 austenite 2 martensite - ₍₎ Macroscopic infinitesimal strain tensor of phase - ₍₂₎ ^f Traceless strain tensor associated with the formation of martensite phase - Macroscopic infiniesimal strain tensor - Macroscopic infinitesimal strain tensor deviator - _f Trace - Equivalent strain - ^pe Macroscopic pseudoelastic strain tensor - x Distortion due to parent (austenite =1)product (martensite =2) phase transformation (traceless symmetric second order tensor) - M Total mass of a system - M() Total mass of phase - V Total volume of a system - V() Total volume of phase - z=M(2)/M Weight fraction of martensite - 1-z=M(1)/M Weight fraction of austenite - u ₀ ^* () Specific internal energy of phase (=1,2) - s ₀ ^* () Specific internal entropy of phase - Specific configurational energy - Specific configurational entropy - ₀ ^f (T) Driving force for temperature-induced martensitic transformation at stress free state ( ₀ ^f T) = T ^*–Ts ^*) - Kirchhoff stress tensor - Kirchhoff stress tensor deviator - Equivalent stress - Cauchy stress tensor - Mass density - K Bulk moduli (K ₀=K) - L Elastic moduli tensor (order 4) - E Young modulus - Energetic shear (₀ = ) - Poisson coefficient - M _s ^o (M _F ^o ) Martensite start (finish) temperature at stress free state - A _s ^o (A _F ^o ) Austenite start (finish) temperature at stress free state - C _v Specific heat at constant volume - k Conductivity - Pseudoelastic strain obtained in tensile test after complete phase transformation (AM) (unidimensional test) - ₀ Thermal expansion tensor - r Resistivity - 1MPa 10⁶ N/m ² - ⁽⁾ Specific free energy of phase - _n Specific free energy at non equilibrium (R model) - _n ^eq Specific free energy at equilibrium (R model) - _n ^v Volumic part of ^eq - Specific free energy at non equilibrium (R _L model) - _conf Specific coherency energy (R _L model) - _c Specific free energy at constrained equilibria (R _L model) - _it (T) Coherency term (R _L model)     

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

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     

Eine stationäre Meßmethode zur Bestimmung der Relaxationszeit viskoelastischer Fluide bei räumlicher Beanspruchung

Zusammenfassung Es wird gezeigt, daß bei Kenntnis der Fließkurve viskoelastischer Flüssigkeiten allein aus der Drehmomentkennlinie des stationär betriebenen Kugel-Kugel-Rheometers eine Relaxationszeit der räumlichen Beanspruchung bestimmt werden kann. Ausgehend von derColeman-Nollschen Entwicklungsschreibweise der rheologischen Zustandsfunktion wird das Geschwindigkeitsfeld als Potenzreihe der Kreisfrequenz bis zur 3. Ordnung bestimmt und zur Drehmomentbeziehung integriert.Messungen an einigen Versuchssubstanzen bestätigen die Tauglichkeit der entwickelten Methode.Häufig verwendete Formelzeichen –a N/m² isotroper Druckanteil - m/s Geschwindigkeitsvektor - e ₁₄ Integrationskonstanten - f _i() Geometriefunktionen - m vektorielle Feldfunktion - ms vektorielle Feldfunktion - ms² vektorielle Feldfunktion - k _i() Geometriefunktionen - t ₀ s Relaxationszeit der räumlichen Beanspruchung - m/s Geschwindigkeitsvektor erster Ordnung - m/s Geschwindigkeitsvektor zweiter Ordnung - m/s Geschwindigkeitsvektor dritter Ordnung - D 1/s Deformationsgeschwindigkeitstensor - 1/s², 1/s³ korotationale, zeitliche Ableitung vonD - 1 Einheitstensor - M Nm Antriebsmoment der rotierenden Kugel - M _i Nm Teilmomente - R m Kugelradius - R _G m Hohlkugelradius - S N/m² Spannungstensor - W 1/s Rotationsgeschwindigkeitstensor - ₁ N s/m² Stoffparameter 1. Ordnung - ₂, ₃ N s²/m² Stoffparameter 2. Ordnung - ₄, ₅, ₆ N s³/m² Stoffparameter 3. Ordnung - RadienverhältnisR/R _G - ₀ N s/m² Anfangsviskosität - kg/m³ Dichte der Flüssigkeit - 1/s Kreisfrequenz der rotierenden Kugel Vorgetragen auf dem 6. Internationalen Rheologie-Kongreß in Lyon-Frankreich vom 4.–8. September 1972.Jetzt: BASF-AG, LudwigshafenMit 4 Abbildungen     

The dynamics of solid-solid phase transitions 2. Incoherent interfaces

Incoherent phase transitions are more difficult to treat than their coherent counterparts. The interface, which appears as a single surface in the deformed configuration, is represented in its undeformed state by a separate surface in each phase. This leads to a rich but detailed kinematics, one in which defects such as vacancies and dislocations are generated by the moving interface. In this paper we develop a complete theory of incoherent phase transitions in the presence of deformation and mass transport, with phase interface structured by energy and stress. The final results are a complete set of interface conditions for an evolving incoherent interface.Frequently used symbols A_i,C_i generic subsurface of S_t - B_i undeformed phase-i region - C configurational bulk stress, Eshelby tensor - F deformation gradient - G inverse deformation gradient - H relative deformation gradient - J bulk Jacobian of the deformation - ¯K, K_i total (twice the mean) curvature of and S_i - Lin (U, V) linear transformations from U into V - Lin⁺ linear transformations of ³ with positive determinant - Orth⁺ rotations of ³ - Q^a external bulk mass supply of species a - ¯S bulk Cauchy stress tensor - S bulk Piola-Kirchhoff stress tensor - S_i undeformed phase i interface - U_i relative velocity of S_i - Unim⁺ linear transformations of ³ with unit determinant - ¯V, V_i normal velocity of and S_i - intrinsic edge velocity of S and A _i S - W_i volume flow across the phase-i interface - X material point - b external body force - e internal bulk configurational force - f_i external interfacial force (configurational) - ¯g external interfacial force (deformational) - grad, div spatial gradient and divergence - gradient and divergence on - h relative deformation - h^a, diffusive mass flux of species a and list of mass fluxes - ¯m outward unit normal to a spatial control volume - ¯n, n_i unit normal to and S_i - n subspace of ³ orthogonal to n - ¯q^a external interfacial mass supply of species a - s ......... - ¯v, v_i compatible velocity fields of and S_i - ¯w, w_i compatible edge velocity fields for and A_i - x spatial point - y_i deformation or motion of phase i - y^. material velocity - generic subsurfaces of - , _i deformed body and deformed phase-i region - () energy supplied to by mass transport - symmetry group of the lattice - _i, surface jacobians - lattice - () power expended on - spatial control volume - S deformed phase interface - lattice point density - interfacial power density - , A total surface stress - C configurational surface stress for phase 1 (material) - ¯C_i configurational surface stress (spatial) - F_i tangential deformation gradient - G_i inverse tangential deformation gradient - H incoherency tensor - ¯1(x), 1_i(X) inclusions of ¯n(x) and n _i (X) into ³ - K configurational surface stress for phase 2 (material) - ¯L, l_i curvature tensor of and S_i - ¯P(x), P_i(X) projections of ³ onto ¯n(x) and n_i (X) - ¯S, S deformational surface stress (spatial and material) - ¯a, a normal part of total surface stress - c normal part of configurational surface stress for phase 1 (material) - e_i internal interfacial configurational force - ¯v, v_i unit normal to and A _i - (x),_i(X) projections of ³ onto ¯n(x) and n _i (X) - _i normal internal force (material) - bulk free energy - slip velocity - _i=(–1)ⁱ _i ......... - ^a, chemical potential of species a and list of potentials - ^a, bulk molar density of species a and list of molar densities - _i normal internal force (spatial) - surface tension - , _i effective shear - referential-to-spatial transform of field - interfacial energy - grand canonical potential - l unit tensor in ³ - x, vector and tensor product in ³ - (...)^., _t(...) material and spatial time derivative - , Div material gradient and divergence - gradient and divergence on S_i - (...), (...) normal time derivative following and S_i - (...) limit of a bulk field asx ,x_i - [...],...> jump and average of a bulk field across the interface - (...)_ext extension of a surface tensor to ³ - tangential part of a vector (tensor) on and S_i     

Two-phase flow in heterogeneous porous media III: Laboratory experiments for flow parallel to a stratified system

Two-phase flow in stratified porous media is a problem of central importance in the study of oil recovery processes. In general, these flows are parallel to the stratifications, and it is this type of flow that we have investigated experimentally and theoretically in this study. The experiments were performed with a two-layer model of a stratified porous medium. The individual strata were composed of Aerolith-10, an artificial: sintered porous medium, and Berea sandstone, a natural porous medium reputed to be relatively homogeneous. Waterflooding experiments were performed in which the saturation field was measured by gamma-ray absorption. Data were obtained at 150 points distributed evenly over a flow domain of 0.1 × 0.6 m. The slabs of Aerolith-10 and Berea sandstone were of equal thickness, i.e. 5 centimeters thick. An intensive experimental study was carried out in order to accurately characterize the individual strata; however, this effort was hampered by both local heterogeneities and large-scale heterogeneities.The theoretical analysis of the waterflooding experiments was based on the method of large-scale averaging and the large-scale closure problem. The latter provides a precise method of discussing the crossflow phenomena, and it illustrates exactly how the crossflow influences the theoretical prediction of the large-scale permeability tensor. The theoretical analysis was restricted to the quasi-static theory of Quintard and Whitaker (1988), however, the dynamic effects described in Part I (Quintard and Whitaker 1990a) are discussed in terms of their influence on the crossflow.Roman Letters A interfacial area between the -region and the -region contained within V, m² - a vector that maps onto , m - b vector that maps onto , m - b vector that maps onto , m - B second order tensor that maps onto , m² - C second order tensor that maps onto , m² - E energy of the gamma emitter, keV - f fractional flow of the -phase - g gravitational vector, m/s² - h characteristic length of the large-scale averaging volume, m - H height of the stratified porous medium , m - i unit base vector in the x-direction - K local volume-averaged single-phase permeability, m² - K - {K}, large-scale spatial deviation permeability - { K} large-scale volume-averaged single-phase permeability, m² - K ^* large-scale single-phase permeability, m² - K ^** equivalent large-scale single-phase permeability, m² - K local volume-averaged -phase permeability in the -region, m² - K local volume-averaged -phase permeability in the -region, m² - K - {K } , large-scale spatial deviation for the -phase permeability, m² - K ^* large-scale permeability for the -phase, m² - l thickness of the porous medium, m - l characteristic length for the -region, m - l characteristic length for the -region, m - L length of the experimental porous medium, m - characteristic length for large-scale averaged quantities, m - n outward unit normal vector for the -region - n outward unit normal vector for the -region - n unit normal vector pointing from the -region toward the -region (n = - n ) - N number of photons - p pressure in the -phase, N/m² - p ₀ reference pressure in the -phase, N/m² - local volume-averaged intrinsic phase average pressure in the -phase, N/m² - large-scale volume-averaged pressure of the -phase, N/m² - large-scale intrinsic phase average pressure in the capillary region of the -phase, N/m² - - , large-scale spatial deviation for the -phase pressure, N/m² - p_c , capillary pressure, N/m² - p _c capillary pressure in the -region, N/m² - p capillary pressure in the -region, N/m² - {p _c } ^c large-scale capillary pressure, N/m² - q -phase velocity at the entrance of the porous medium, m/s - q -phase velocity at the entrance of the porous medium, m/s - S_wi irreducible water saturation - S /, local volume-averaged saturation for the -phase - S _i initial saturation for the -phase - S _r residual saturation for the -phase - S ^* { }^*/}^*, large-scale average saturation for the -phase - S saturation for the -phase in the -region - S saturation for the -phase in the -region - t time, s - v -phase velocity vector, m/s - v local volume-averaged phase average velocity for the -phase, m/s - {v } large-scale averaged velocity for the -phase, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - V large-scale averaging volume, m³ - y position vector relative to the centroid of the large-scale averaging volume, m - {y}^c large-scale average of y over the capillary region, m Greek Letters local porosity - local porosity in the -region - local porosity in the -region - local volume fraction for the -phase - local volume fraction for the -phase in the -region - local volume fraction for the -phase in the -region - {}^* { }^*+{ }^*, large-scale spatial average volume fraction - { }^* large-scale spatial average volume fraction for the -phase - mass density of the -phase, kg/m³ - mass density of the -phase, kg/m³ - viscosity of the -phase, N s/m² - viscosity of the -phase, Ns/m² - V /V , volume fraction of the -region ( + =1) - V /V , volume fraction of the -region ( + =1) - attenuation coefficient to gamma-rays, m^-1 - -      

Stochastic response of structures with small geometric imperfections

Summary A probabilistic model of the geometric imperfections of a real structure is proposed, in order to provide a general theory of the stochastic response of structures in presence of small random deviations from the perfect scheme. The main statistical measures of the stochastic response are derived and an application to the study of a particular conservative elastic system is developed.

---

Sommario Si propone una teoria generale della risposta probabilistica di strutture, in presenza di piccole deviazioni aleatorie dei dati iniziali rispetto allo schema geometrico perfetto. Si deducono le principali proprietà statistiche della risposta della struttura a sollecitazioni esterne deterministiche, e si sviluppa una applicazione riguardante il comportamento aleatorio di un particolare sistema elastico conservativo.

---

List of symbols element of the sample space of events - _kn random variables modelling the structural imperfections - P(o) probability density of random variables - random imperfection of the unloaded structure - u additional displacement of the loaded structure - u_o deterministic fundamental solution for the perfect structure - difference between the additional displacement of the loaded structure and the deterministic fundamental solution for the perfect structure - V₁=u₁ buckling mode of the perfect structure - _i intrinsic coordinates of the structure - suitable measure of the magnitude of the random imperfections - scalar geometric variable representing the internal product - random imperfection divided by - single scalar variable denoting the magnitude of the prescribed loads - potential energy of the structure - potential energy of the perfect structure - difference between and - _c lowest critical load - _s real local maximum for the magnitude of the prescribed loads - _c divided by _S - E{} expected value of a random variable - ² variance of a random variable - , random variables defined by Eq. (21)     

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.     

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     

Flow in porous media I: A theoretical derivation of Darcy's law

Stokes flow through a rigid porous medium is analyzed in terms of the method of volume averaging. The traditional averaging procedure leads to an equation of motion and a continuity equation expressed in terms of the volume-averaged pressure and velocity. The equation of motion contains integrals involving spatial deviations of the pressure and velocity, the Brinkman correction, and other lower-order terms. The analysis clearly indicates why the Brinkman correction should not be used to accommodate ano slip condition at an interface between a porous medium and a bounding solid surface.The presence of spatial deviations of the pressure and velocity in the volume-averaged equations of motion gives rise to aclosure problem, and representations for the spatial deviations are derived that lead to Darcy's law. The theoretical development is not restricted to either homogeneous or spatially periodic porous media; however, the problem ofabrupt changes in the structure of a porous medium is not considered.Roman Letters A interfacial area of the - interface contained within the macroscopic system, m² - A _e area of entrances and exits for the -phase contained within the macroscopic system, m² - A interfacial area of the - interface contained within the averaging volume, m² - A ^* interfacial area of the - interface contained within a unit cell, m² - Ae area of entrances and exits for the -phase contained within a unit cell, m² - B second order tensor used to represent the velocity deviation (see Equation (3.30)) - b vector used to represent the pressure deviation (see Equation (3.31)), m^–1 - d distance between two points at which the pressure is measured, m - g gravity vector, m/s² - K Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the -phase (see Figure 2), m - characteristic length scale for the -phase (see Figure 2), m - n unit normal vector pointing from the -phase toward the -phase (n =–n ) - n _e unit normal vector for the entrances and exits of the -phase contained within a unit cell - p pressure in the -phase, N/m² - p intrinsic phase average pressure for the -phase, N/m² - p –p , spatial deviation of the pressure in the -phase, N/m² - r ₀ radius of the averaging volume and radius of a capillary tube, m - v velocity vector for the -phase, m/s - v phase average velocity vector for the -phase, m/s - v intrinsic phase average velocity vector for the -phase, m/s - v –v , spatial deviation of the velocity vector for the -phase, m/s - V averaging volume, m³ - V volume of the -phase contained within the averaging volume, m³ Greek Letters V/V, volume fraction of the -phase - mass density of the -phase, kg/m³ - viscosity of the -phase, Nt/m² - arbitrary function used in the representation of the velocity deviation (see Equations (3.11) and (B1)), m/s - arbitrary function used in the representation of the pressure deviation (see Equations (3.12) and (B2)), s^–1     

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