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²     

Determination of permeability tensors for two-phase flow in homogeneous porous media: Theory

In this paper we continue previous studies of the closure problem for two-phase flow in homogeneous porous media, and we show how the closure problem can be transformed to a pair of Stokes-like boundary-value problems in terms of pressures that have units of length and velocities that have units of length squared. These are essentially geometrical boundary value problems that are used to calculate the four permeability tensors that appear in the volume averaged Stokes' equations. To determine the geometry associated with the closure problem, one needs to solve the physical problem; however, the closure problem can be solved using the same algorithm used to solve the physical problem, thus the entire procedure can be accomplished with a single numerical code.Nomenclature a a vector that maps V onto , m^-1. - A a tensor that maps V onto . - A area of the - interface contained within the macroscopic region, m². - A area of the -phase entrances and exits contained within the macroscopic region, m². - A area of the - interface contained within the averaging volume, m². - A area of the -phase entrances and exits contained within the averaging volume, m². - Bo Bond number (= (=(–)g²/). - Ca capillary number (= v/). - g gravitational acceleration, m/s². - H mean curvature, m^-1. - I unit tensor. - permeability tensor for the -phase, m². - viscous drag tensor that maps V onto V. - ^* dominant permeability tensor that maps onto v , m². - ^* coupling permeability tensor that maps onto v , m². - characteristic length scale for the -phase, m. - l characteristic length scale representing both and , m. - L characteristic length scale for volume averaged quantities, m. - n unit normal vector directed from the -phase toward the -phase. - n unit normal vector representing both n and n . - n unit normal vector representing both n and n . - P pressure in the -phase, N/m². - p superficial average pressure in the -phase, N/m². - p intrinsic average pressure in the -phase, N/m². - p –p , spatial deviation pressure for the -phase, N/m². - r ₀ radius of the averaging volume, m. - r position vector, m. - t time, s. - v fluid velocity in the -phase, m/s. - v superficial average velocity in the -phase, m/s. - v intrinsic average velocity in the -phase, m/s. - v –v , spatial deviation velocity in the -phase, m/s. - V volume of the -phase contained within the averaging volmue, m³. - averaging volume, m³. Greek Symbols V /, volume fraction of the -phase. - viscosity of the -phase, Ns/m². - density of the -phase, kg/m³. - surface tension, N/m. - (v +v ^T ), viscous stress tensor for the -phase, N/m².     

Filament stretching equation in Lagrangian coordinates

The general integral of the very simple equation ^21/n/() was found to describe the cross sectional area of filaments of isothermal power law fluids while in transient stretching where is time and is the initial location of fluid molecules at time = 0 given as the distance from a reference point fixed in space. Any such stretching transient given as a solution of the above equation is physically realizable subject to the restrictions > 0 and/ < 0.     

Mean and turbulent velocity measurements of supersonic mixing layers

The behavior of supersonic mixing layers under three conditions has been examined by schlieren photography and laser Doppler velocimetry. In the schlieren photographs, some large-scale, repetitive patterns were observed within the mixing layer; however, these structures do not appear to dominate the mixing layer character under the present flow conditions. It was found that higher levels of secondary freestream turbulence did not increase the peak turbulence intensity observed within the mixing layer, but slightly increased the growth rate. Higher levels of freestream turbulence also reduced the axial distance required for development of the mean velocity. At higher convective Mach numbers, the mixing layer growth rate was found to be smaller than that of an incompressible mixing layer at the same velocity and freestream density ratio. The increase in convective Mach number also caused a decrease in the turbulence intensity ( _u/U).List of symbols a speed of sound - b total mixing layer thickness between U ₁ – 0.1 U and U ₂ + 0.1 U - f normalized third moment of u-velocity, f u³/(U)³ - g normalized triple product of u² , g u²/(U)³ - h normalized triple product of u ², h u ²/(U)³ - l _u axial distance for similarity in the mean velocity - l _u axial distance for similarity in the turbulence intensity - M Mach number - M _c convective Mach number (for ₁ = ₂), M _c (U ₁ – U ₂)/(a ₁ + a ₂) - P static pressure - r freestream velocity ratio, r U ₂/U ₁ - Re unit Reynolds number, Re U/ - s freestream density ratio, s ₂/₁ - T _t total temperature - u instantaneous streamwise velocity - u deviation of u-velocity, uu – U - U local mean streamwise velocity - U ₁ primary freestream velocity - U ₂ secondary freestream velocity - average of freestream velocities, (U ₁ + U ₂)/2 - U freestream velocity difference, U U ₁ – U ₂ - instantaneous transverse velocity - v deviation of -velocity, – V - V local mean transverse velocity - x streamwise coordinate - y transverse coordinate - y ₀ transverse location of the mixing layer centerline - ensemble average - ratio of specific heats - boundary layer thickness (y-location at 99.5% of free-stream velocity) - similarity coordinate, (y – y ₀)/b - compressible boundary layer momentum thickness - viscosity - density - standard deviation - dimensionless velocity, (U – U ₂)/U - 1 primary stream - 2 secondary stream A version of this paper was presented at the 11th Symposium on Turbulence, October 17–19, 1988, University of Missouri-Rolla     

Flow in porous media II: The governing equations for immiscible,two-phase flow

The Stokes flow of two immiscible fluids through a rigid porous medium is analyzed using the method of volume averaging. The volume-averaged momentum equations, in terms of averaged quantities and spatial deviations, are identical in form to that obtained for single phase flow; however, the solution of the closure problem gives rise to additional terms not found in the traditional treatment of two-phase flow. Qualitative arguments suggest that the nontraditional terms may be important when / is of order one, and order of magnitude analysis indicates that they may be significant in terms of the motion of a fluid at very low volume fractions. The theory contains features that could give rise to hysteresis effects, but in the present form it is restricted to static contact line phenomena.Roman Letters (, = , , and ) 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² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - g gravity vector, m²/s - H mean curvature of the- interface, m^–1 - H area average of the mean curvature, m^–1 - H – H , deviation of the mean curvature, m^–1 - I unit tensor - K Darcy's law permeability tensor, m² - K permeability tensor for the-phase, m² - K viscous drag tensor for the-phase equation of motion - K viscous drag tensor for the-phase equation of motion - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - n unit normal vector pointing from the-phase toward the-phase (n = –n ) - p _c p – P , capillary pressure, N/m² - 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, m - t time, s - 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² - surface tension of the- interface, N/m - viscous stress tensor for the-phase, N/m² - / kinematic viscosity, m²/s     

Wärmeübergang turbulenter kompressibler Grenz-schichtströmungen mit starken negativen Druckgradienten (Wärmeübergang in einer Lavaldüse)

Zusammenfassung Es wurde eine zweidimensionale, aus ebener Wand und Konturwand bestehende Lavaldüse gebaut, um an den Wänden (unterteilt in gut isolierte und einzeln gekühlte Segmente) genaue Wärmeübergangs messungen durchzuführen. Alle Versuchspunkte zeigen bei geringer Streubreite eine eindeutige Abhängigkeit der Stanton-Zahl von einer aus neu entwickelten Temperaturprofilen gewonnenen Differenz aus dimensionslosen Kennzahlen für Dissipation und Druckgradient. Bei gekühlter Grenzschicht und zunehmendem negativen Druckgradienten nimmt in Theorie und Versuch bei geringer Dissipation der Wärmeübergang ab, bei stärkerer Dissipation zu.

---

Compressible turbulent boundary layer heat transfer with strong favourable pressure gradients (heat transfer in a convergent — Divergent nozzle)

A two-dimensional convergent-divergent nozzle with a plane and a contoured wall was built to perform exact measurements of heat transfer to the walls divided in well insulated single-cooled segments. For all experimental data within a small error-band the Stanton-number definitely depends on a difference of dimensionless numbers for dissipation and pressure gradient found from new developed temperature profiles. Theoretical and experimental results show that in a cooled boundary layer with increasing favourable pressure gradient heat transfer decreases with small dissipation and increases with greater dissipation.

---

Formelzeichen Wärmeübergangsparameter - c_f=_w/(·u ² örtlicher Reibungsbeiwert - c_p spezifische Wärmekapazität - d⁺= u ² /(c_p·T_q) Dissipationsparameter - 1_t turbulenter Mischungsweg - m Massenstrom im Windkanal - p Druck - p⁺= –K/c_f ^3/2 Druckgradientenparameter - q Wärmestromdichte - r Rückgewinnfaktor - u Geschwindigkeit in Hauptstromrichtung - Schubspannungsgeschwin digkeit - u⁺=u/u dimensionslose Geschwin digkeit - x Lauflänge an der Wand - y Abstand senkrecht zur Wand - y⁺=· u_y/ dimensionslose y-Koordinate - C Wärmeübergangsbeiwert - Euler-Zahl - H_s Höhe des engsten Quer-Schnitts - K Beschleunigungsparameter - Ma Machzahl - Pr Prandtl-Zahl - Pr_t turbulente Prandtl-Zahl - R Radius - Reynolds-Zahl - Stanton-Zahl - T absolute Temperatur - T_ei Eigentemperatur - Wärmestromdichtetemperatur - T⁺=T/T_w dimensionslose Temperatur - T^*=l/b⁺ · (1 –T⁺) dimensionslose Temperatur - _ei Wärmeübergangskoeffizient bezogen auf Eigentemperatur - dynamische Zähigkeit - =T – 273,15K Celsius-Temperatur - v kinematische Zähigkeit - Dichte - Schubspannung - Differenz - Wärmeübergangsparameter Indices Freistrom - a Austritt - e Eintritt - k kühlseitige Wand lam laminar - O Gesamt-, Ruhew Wand - x auf Lauflänge x bezogen Auszug aus der vom Fachbereich Maschinenwesen der Technischen Universität München zur Erlangung des akademischen Grades eines Doktor-Ingenieur genehmigten Dissertation über Wärmeübergang turbulenter kompressibler Grenzschichtströmungen mit starken negativen Druckgradienten (Wärmeübergang in einer Lavaldüse) des Diplom-Ingenieurs W. Winkler. Berichterstatter Prof. Dr.-Ing. U. Grigull und Prof. Dr.-Ing. E. Truckenbrodt. — Die Dissertation wurde am 17.12.1975 bei der Technischen Universität München und am 1.4.1976 durch deren Fachbereich Maschinenwesen angenommen. Tag der Promotion 13.5.1976früher: Lehrstuhl A für Thermodynamik Technische Universität München     

The dynamic performance of the Weissenberg rheogoniometer III. Large amplitude oscillatory shearing — harmonic analysis

The harmonic content of the nonlinear dynamic behaviour of 1% polyacrylamide in 50% glycerol/water was studied using a standard Model R 18 Weissenberg Rheogoniometer. The Fourier analysis of the Oscillation Input and Torsion Head motions was performed using a Digital Transfer Function Analyser.In the absence of fluid inertia effects and when the amplitude of the (fundamental) Oscillation Input motion _I is much greater than the amplitudes of the Fourier components of the Torsion Head motion _Tn empirical nonlinear dynamic rheological propertiesG _n (, ⁰),G _n (, ⁰) and/or _n (, ⁰), _n (, ⁰) may be evaluated without a-priori-knowledge of a rheological constitutive equation. A detailed derivation of the basic equations involved is presented.Cone and plate data for the third harmonic storage modulus (dynamic rigidity)G ₃ (, ⁰), loss modulusG ₃ (, ⁰) and loss angle ₃ (, ⁰) are presented for the frequency range 3.14 × 10^–2 1.25 × 10² rad/s at two strain amplitudes, _CP ⁰ = 2.27 and 4.03. Composite cone and plate and parallel plates data for both the third and fifth harmonic dynamic viscosities ₃ (, ⁰), _S (, ⁰) and dynamic rigiditiesG ₃ (, ⁰),G ₅ (, ⁰) are presented for strain amplitudes in the ranges 1.10 _CP ⁰ 4.03 and 1.80 _PP ⁰ 36 for a single frequency, = 3.14 × 10^–1 rad/s. Good agreement was obtained between the results from both geometries and the absence of significant fluid inertia effects was confirmed by the superposition of the data for different gap widths.     

Laser velocimetry measurements in a horizontal gas-solid pipe flow

This paper presents laser measurements of particle velocities in a horizontal turbulent two-phase pipe flow. A phase Doppler particle analyzer, (PDPA), was used to obtain particle size, velocity, and rms values of velocity fluctuations. The particulate phase consisted of glass spheres 50 m in diameter with the volume fraction of the suspension in the range _p=10^-4 to _p=10^-3. The results show that the turbulence increases with particle loading.List of symbols a particle diameter - C _va velocity diameter cross-correlation - d pipe diameter - Fr ² Froude number - g gravitational constant - p(a) Probability density of the particle diameter - Re pipe Reynolds number based on the friction velocity - T characteristic time scale of the energy containing eddies - T _L integral scale of the turbulence sampled along the particle path - u, U, u characteristic fluid velocities: fluctuating, mean and friction - v characteristic velocity of the paricle fluctuations - f expected value of any random variable f - f¦g expected value of f given a value of the random variable g - _p particle volume fraction - _p particle response time - absolute fluid viscosity - v kinematic fluid viscosity - _p, _f densities, particle and fluid - _a ² particle diameter variance - _va ² velocity variance due to the particle diameter variance - _vT ² total particle velocity variance - _vt ² particle velocity variance due to the response to the turbulent field     

Free convection heat transfer from an isothermal plate with arbitrary inclination

This paper presents a new formulation for the laminar free convection from an arbitrarily inclined isothermal plate to fluids of any Prandtl number between 0.001 and infinity. A novel inclination parameter is proposed such that all cases of the horizontal, inclined and vertical plates can be described by a single set of transformed equations. Moreover, the self-similar equations for the limiting cases of the horizontal and vertical plates are recovered from the transformed equations by setting=0 and=1, respectively. Heated upward-facing plates with positive and negative inclination angles are investigated. A very accurate correlation equation of the local Nusselt number is developed for arbitrary inclination angle and for 0.001 Pr .

---

Wärmeübertragung bei freier Konvektion an einer isothermen Platte mit beliebiger Neigung

Zusammenfasssung Diese Untersuchung stellt eine neue Formulierung der laminaren freien Konvektion von Flüssigkeiten mit einer Prandtl-Zahl zwischen 0,001 und unendlich an einer beliebig schräggestellten isothermen Platte dar. Ein neuer Neigungsparameter wird eingeführt, so daß alle Fälle der horizontalen, geneigten oder vertikalen Platte von einem einzigen Satz transformierter Gleichungen beschrieben werden können. Die unabhängigen Gleichungen für die beiden Fälle der horizontalen and vertikalen Platte wurden für=0 und=1 aus den transformierten Gleichungen wieder abgeleitet. Es wurden erwärmte aufwärtsgerichtete Platten mit positiven und negativen Neigungswinkeln untersucht. Eine sehr genaue Gleichung wurde für die lokale Nusselt-Zahl bei beliebigen Neigungswinkeln und für 0,001 Pr entwickelt.

---

Nomenclature C _p specific heat - f reduced stream function - g gravitational acceleration - Gr local Grashof number,g(T _w –T _w ) x³/v² - h local heat transfer coefficient - k thermal conductivity - n constant exponent - Nu local Nusselt number,hx/k - p pressure - Pr Prandtl number, v/ - Ra local Rayleigh number,g(T _w –T )J x³/v - T fluid temperature - T _w wall temperature - T temperature of ambient fluid - u velocity component in x-direction - v velocity component in y-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - (Ra¦sin¦)^1/4/( Ra cos()^1/5 - pseudo-similarity variable, (y/) - dimensionless temperature, (T–T )/(T _w–T ) - ( Ra cos)^1/5+(Rasin)^1/4 - v kinematic viscosity - 1/[1 +(Ra cos)^1/5/( Ra¦sin)^1/4] - density of fluid - Pr/(1+Pr) - _w wall shear stress - angle of plate inclination measured from the horizontal - stream function - dimensionless dynamic pressure     

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

Existence theorem for a minimum problem with free discontinuity set

We study the variational problem Where is an open set in ⁿ ,n2gL ^q () L (), 1q<+, O<, <+ andH _n–1 is the (n–1)-dimensional Hausdorff Measure.     

Robustness of Exponential Dichotomies in Infinite-Dimensional Dynamical Systems

In this paper we examine the issue of the robustness, or stability, of an exponential dichotomy, or an exponential trichotomy, in a dynamical system on an Banach space W. These two hyperbolic structures describe long-time dynamical properties of the associated time-varying linearized equation _t +A=B(t) , where the linear operator B(t) is the evaluation of a suitable Fréchet derivative along a given solution in the set K in W. Our main objective is to show, under reasonable conditions, that if B(t)=B(, t) depends continuously on a parameter and there is an exponential dichotomy, or exponential trichotomy, at a value ₀, then there is an exponential dichotomy, or exponential trichotomy, for all near ₀.We present several illustrations indicating the significance of this robustness property.     

Reservoir transport equations by compositional approach

The objective of this paper is to present an overview of the fundamental equations governing transport phenomena in compressible reservoirs. A general mathematical model is presented for important thermo-mechanical processes operative in a reservoir. Such a formulation includes equations governing multiphase fluid (gas-water-hydrocarbon) flow, energy transport, and reservoir skeleton deformation. The model allows phase changes due to gas solubility. Furthermore, Terzaghi's concept of effective stress and stress-strain relations are incorporated into the general model. The functional relations among various model parameters which cause the nonlinearity of the system of equations are explained within the context of reservoir engineering principles. Simplified equations and appropriate boundary conditions have also been presented for various cases. It has been demonstrated that various well-known equations such as Jacob, Terzaghi, Buckley-Leverett, Richards, solute transport, black-oil, and Biot equations are simplifications of the compositional model.Notation List B reservoir thickness - B formation volume factor of phase - Cⁱ mass of component i dissolved per total volume of solution - C ⁱ mass fraction of component i in phase - C heat capacity of phase at constant volume - C_p heat capacity of phase at constant pressure - D ⁱ hydrodynamic dispersion coefficient of component i in phase - D_MTf thermal liquid diffusivity for fluid f - F = F(x, y, z, t) defines the boundary surface - f_p fractional flow of phase - g gravitational acceleration - H_p enthalpy per unit mass of phase - J_p volumetric flux of phase - k_rf relative permeability to fluid f - k₀ absolute permeability of the medium - M_p ⁱ mass of component i in phase - n porosity - N rate of accretion - P_f pressure in fluid f - p_ca capillary pressure between phases and =p-p - Rⁱ rate of mass transfer of component i from phase to phase - Rⁱ _source source rate of component i within phase - S saturation of phase - s gas solubility - T temperature - t time - U displacement vector - u velocity in the x-direction - v velocity in the y-direction - V volume of phase - V_s velocity of soil solids - W_i body force in coordinate direction i - x horizontal coordinate - z vertical coordinate Greek Letters _p volumetric coefficient of compressibility - _T volumetric coefficient of thermal expansion - _ij Kronecker delta - volumetric strain - _m thermal conductivity of the whole matrix - internal energy per unit mass of phase - gf suction head - density of phase - _ij tensor of total stresses - _ij tensor of effective stresses - volumetric content of phase - f viscosity of fluid f     

Über die Stabilität viskometrischer Strömungen

Zusammenfassung Die Stabilität der ebenen Couette- und der ebenen Poiseuille-Strömung nicht-newtonscher Fluide wird für kleine Störungen in der viskometrischen Ebene untersucht. Der Einfluß der Relaxationszeit der Störungen wird vernachlässigt. Es wird gezeigt, daß die ebene Couette-Strömung unabhängig von der ReZahl instabil wird, fallsd(N)/d > 4 _>d gilt. Hier bedeuten die Schergeschwindigkeit,N den ersten Normalspannungskoeffizienten, die Viskosität und _d die differentielle Viskosität ( _d =d/d). Das gleiche Kriterium gilt mit den Daten an der Kanalwand auch für die Poiseuille-Strömung. In diesem Fall oszillieren die Eigenfunktionen in einer sehr dünnen, wandnahen Schicht und klingen im Flüssigkeitsinnern sehr rasch ab.

---

Summary The stability of plane Couette and plane Poiseuille flow of a non-Newtonian fluid is investigated for small perturbations in the viscometric plane. The influence of the relaxation time of the perturbations is neglected. It is shown that plane Couette flow will become unstable independently of Reynolds number ifd(N)/d > 4 _d holds. Here are the rate of shear velocity,N the first normal stress coefficient, the viscosity and _d the differential viscosity ( _d =d/d). The same criterion holds also for plane Poiseuille flow with the data taken at the wall. In this case the eigenfunctions are oscillating in a very thin layer near the wall and decaying very rapidly in the inner region of the flow field.

Mit 11 Abbildungen     

Flow visualization of compressible vortex structures using density gradient techniques

Mathematical results are derived for the schlieren and shadowgraph contrast variation due to the refraction of light rays passing through two-dimensional compressible vortices with viscous cores. Both standard and small-disturbance solutions are obtained. It is shown that schlieren and shadowgraph produce substantially different contrast profiles. Further, the shadowgraph contrast variation is shown to be very sensitive to the vortex velocity profile and is also dependent on the location of the peak peripheral velocity (viscous core radius). The computed results are compared to actual contrast measurements made for rotor tip vortices using the shadowgraph flow visualization technique. The work helps to clarify the relationships between the observed contrast and the structure of vortical structures in density gradient based flow visualization experiments.Nomenclature a Unobstructed height of schlieren light source in cutoff plane, m - c Blade chord, m - f Focal length of schlieren focusing mirror, m - C _T Rotor thrust coefficient, T/( ² R ⁴) - I Image screen illumination, Lm/m ² - l Distance from vortex to shadowgraph screen, m - n _b Number of blades - p Pressure,N/m ² - p Ambient pressure, N/m ² - r, , z Cylindrical coordinate system - r _c Vortex core radius, m - Non-dimensional radial coordinate, (r/r _c ) - R Rotor radius, m - Tangential velocity, m/s - Specific heat ratio of air - Circulation (strength of vortex), m ²/s - Non-dimensional quantity, ² 8²p r _c ² - Refractive index of fluid medium - ₀ Refractive index of fluid medium at reference conditions - Gladstone-Dale constant, m ³/kg - Density, kg/m ³ - Density at ambient conditions, kg/m ³ - Non-dimensional density, (/ ) - Rotor solidity, (n _b c/ R) - Rotor rotational frequency, rad/s     

Modell zur Beschreibung des Fließens wandgleitender Substanzen durch Düsen

Zusammenfassung Zur Analyse des Fließens einer direkt an der Düsenwand gleitendenOstwald-deWaele-Flüssigkeit (Potenzgesetz) wird ein Modell entwickelt, das die rheologischen Vorgänge tribologisch, d. h. analog derCoulombschen Reibung fester Körper beschreibt.Es zeigt sich, daß in der Düse zwei Bereiche zu unterscheiden sind: ein Haftbereich in der Nähe des Düseneinlaufs und ein am Düsenaustritt liegender Gleitbereich. Die Länge des Gleitbereichs, der Verlauf des Drucks und der Schubspannung längs der Düse sowie die Änderung des Geschwindigkeitsprofils im Gleitbereich werden ermittelt.Überschreitet die Wandschubspannung einen kritischen Betrag, so entsteht am Düsenende ein labiler Bereich, in dem der Betrag der Wandschubspannung sprunghaft auf einen kleineren Wert sinken kann. Der von verschiedenen Autoren gefundene Sprung in der Fließkurve bestimmter Polymerschmelzen kann damit grundsätzlich erklärt werden.

---

Summary Starting from theCoulomb Friction Law for solids, a theoretical model is developed for the pressure flow of a viscous power-law fluid with slip at the wall.It is shown that two flow regions exist in the die: a first region at the upstream part of the die, where the fluid sticks to the wall; and a second region at the downstream part of the die, where the fluid slips at the wall. The length of the slip region, the development of pressure and shear stress along the die as well as the change of the velocity distribution are given for the slip region.For shear stresses above a critical value, an instability region is found at the exit of the die. In this region, a sudden decrease of shear stress can occur. This seems to explain the discontinuity in the flow curve reported by several investigators.

---

F Querschnittsfläche der Kapillaren - Volumendurchsatz - K _R Reibkraft - L Düsenlänge - m Stoffwert (Fließexponent) - N Normalkraft - p hydrostatischer Druck - p _L Druck am Düsenende - p ₁ Druck an der Übergangsstelle Haften-Gleiten - p ₀ Druck vor der Düse - p _0H Druck vor der Düse im Falle des Wandhaftens - r Radius - R Düsenradius - v _g Gleitgeschwindigkeit - v _z Strömungsgeschwindigkeit inz-Richtung - z Koordinate in Strömungsrichtung - z ₁ Längskoordinate der Übergangsstelle Haften-Gleiten - Schergeschwindigkeit - Stoffwert - Viskosität - µ Gleitkoeffizient - µ _H Haftkoeffizient - Dichte - dimensionsloser Radiusr/R - _rz Schubspannung in der Flüssigkeit - _rz ^(R) Wandschubspannung in der Flüssigkeit - ₀ Stoffwert - _wg Wandschubspannung im Falle des Gleitens - _wH Haftschubspannung an der Wand Auszugsweise vorgetragen auf der Jahrestagung der Deutschen Rheologen in Berlin vom 28.–30. April 1975.Mit 10 Abbildungen     

On surface waves in generalized thermoelasticity

Knowles' representation theorem for harmonically time-dependent free surface waves on a homogeneous, isotropic elastic half-space is extended to include harmonically time-dependent free processes for thermoelastic surface waves in generalized thermoelasticity of Lord and Shulman and of Green and Lindsay.r , , r , , .This work was done when author was unemployed.     

Mixed convection from a horizontal plate to fluids of any Prandtl number

Laminar mixed convection over a horizontal plate with uniform wall temperature or uniform wall heat flux is analyzed by introducing proper buoyancy parameters and transformation variables for fluids of any Prandtl number between 0.001 and 10,000. Both cases of buoyancy assisting and opposing flow conditions are investigated. For the buoyancy-assisting case, the obtained numerical results are very accurate over the entire range of mixed convection intensity from pure forced convection limit to pure free convection limit. For the buoyancy-opposing case, solutions are obtained from the forced convection limit to the point of breakdown.

---

Mischkonvektion an einer horizontalen Platte für Fluide mit beliebiger Prandtl-Zahl

Zusammenfassung Es wurde laminare Mischkonvektion an einer horizontalen Platte mit einheitlicher Wandtemperatur oder einheitlicher Wandwärmestromdichte bei Einführung zweckmäßiger Auftriebsparameter und Transformationsvariablen für Fluide mit beliebiger Prandtl-Zahl zwischen 0,001 und 10 000 untersucht. Es wurden die Fälle der Strömung entgegen und in Richtung der Auftriebskraft untersucht. Für den Fall der Strömung in Richtung der Auftriebskraft wurden sehr genaue numerische Ergebnisse für den gesamten Bereich der gemischten Konvektion von rein erzwungener Konvektion bis zu rein freier Konvektion erhalten. Für den Fall der Strömung entgegen der Auftriebsrichtung wurden Lösungen für erzwungene Konvektion bis zum Umkehrpunkt erhalten.

---

Nomenclature C _f local friction coefficient - f reduced stream function - g gravitational acceleration - Gr local Grashof number for UWT,g (T _w –T )x ³/ ² - Gr* local Grashof number for UHF,g q _w x ⁴/k ² - m =10 for UWT; and =6 for UHF - n =5 for UWT; and =3 for UHF - Nu local Nusselt number - p pressure - Pr Prandtl number,/ - q _w wall heat flux - Ra local Rayleigh number for UWT,Gr Pr - Ra* local Rayleigh number for UHF,Gr*Pr - Re local Reynolds number,u x/ - T fluid temperature - T _w wall temperature - T free-stream temperature - u velocity component inx-direction - u free-stream velocity - v velocity component iny-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - =0 for UWT; and =1 for UHF - buoyancy parameter, =( Ra)^1/5/( Re)^1/2 for UWT; and =( Ra*)^1/6/( Re)^1/2 for UHF - pseudo-similarity variable, (y/x) - dimensionless temperature, =(T–T )/(T _w –T ) for UWT; and =(T–T )/(q _w x/k) for UHF - =[( Re)^1/2+( Ra)^1/5] for UWT; and =[( Re)^1/2+( Ra*)^1/6] for UHF - dynamic viscosity - kinematic viscosity - /(1+) - dimensionless pressure - density - Pr/(1+Pr) - _w wall shear stress,(u/y) _y=0 - stream function - Pr/(1+Pr)^1/3     

The effect of background particulates on the delayed transition of a heated 9:1 ellipsoid

An experimental study was done to quantify the effects of a variety of background particulates on the delayed laminar-turbulent transition of a thermally stabilized boundary layer in water. A Laser-Doppler Velocimeter system was used to measure the location of boundary layer transition on a 50 mm diameter, 9:1 fineness ratio ellipsoid. The ellipsoid had a 0.15 m RMS surface finish. Boundary layer transition locations were determined for length Reynolds numbers ranging from 3.0 × 10⁶ to 7.5 × 106. The ellipsoid was tested in three different heating conditions in water seeded with particles of four distinct size ranges. For each level of boundary layer heating, measurements of transition were made for clean water and subsequently, water seeded with 12.5 m, 38.9 m, 85.5 m and 123.2 m particles, alternately. The three surface heating conditions tested were no heating, T = 10°C and T = 15°C where T is the difference between the inlet model heating water temperature, T _i, and free stream water temperature, T . The effects of particle concentration were studied for 85.5 m and 123.2 m particulates.The results of the study can be summarized as follows. The 12.5 m and 38.9 m particles has no measurable effect on transition for any of the test conditions. However, transition was significantly affected by the 85.5 m and 123.2 m particles. Above a length Reynolds number of 4 × 10⁶ the boundary layer transition location moved forward on the body due to the effect of the 85.5 m particles for all heating conditions. The largest percentage changes in transition location from clean water, were observed for 85.5 m particles seeded water.Transition measurements made with varied concentrations of background particulates indicated that the effect of the 85.5 m particles on the transition of the model reached a plateau between 2.65 particulates/ml concentration and 4.2 particles/ml. Measurements made with 123.3 m particles at concentrations up to 0.3 part/ml indicated no similar plateau.     

Normal forms for random diffeomorphisms

Given a dynamical system (,, ,) and a random diffeomorphism (): ^{d d} with fixed point at x=0. The normal form problem is to construct a smooth near-identity nonlinear random coordinate transformation h() to make the random diffeomorphism ()=h()^–1() h() as simple as possible, preferably linear. The linearization D(, 0)=:A() generates a matrix cocycle for which the multiplicative ergodic theorem holds, providing us with stochastic analogues of eigenvalues (Lyapunov exponents) and eigenspaces. Now the development runs pretty much parallel to the deterministic one, the difference being that the appearance of turns all problems into infinite-dimensional ones. In particular, the range of the homological operator is in general not closed, making the conceptof-normal form necessary. The stochastic versions of resonance and averaging are developed. The case of simple Lyapunov spectrum is treated in detail.