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

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

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

Singular solutions of some quasilinear elliptic equations

We study isolated singularities of the quasilinear equation in an open set of ^N , where 1 < p N, p -1 q < N(p — 1)/ (N -p). We prove that, for any positive solution, if a singularity at the origin is not removable then either or u(x)/(x) any positive constant as x 0 where is the fundamental solution of the p-harmonic equation: . Global positive solutions are also classified.     

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     

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

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

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     

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.     

Two-phase flow in heterogeneous porous media II: Numerical experiments for flow perpendicular to a stratified system

In this work, we make use of numerical experiments to explore our original theoretical analysis of two-phase flow in heterogeneous porous media (Quintard and Whitaker, 1988). The calculations were carried out with a two-region model of a stratified system, and the parameters were chosen be consistent with practical problems associated with groundwater flows and petroleum reservoir recovery processes. The comparison between theory (the large-scaled averaged equations) and experiment (numerical solution of the local volume averaged equations) has allowed us to identify conditions for which the quasi-static theory is acceptable and conditions for which a dynamic theory must be used. Byquasi-static we mean the following: (1) The local capillary pressure,everywhere in the averaging volume, can be set equal to the large-scale capillary pressure evaluated at the centroid of the averaging volume and (2) the large-scale capillary pressure is given by the difference between the large-scale pressures in the two immiscible phases, and is therefore independent of gravitational effects, flow effects and transient effects. Bydynamic, we simply mean a significant departure from the quasi-static condition, thus dynamic effects can be associated with gravitational effects, flow effects and transient effects. To be more precise about the quasi-static condition we need to refer to the relation between the local capillary pressure and the large-scale capillary pressure derived in Part I (Quintard and Whitaker, 1990). Herep _c ¦_y represents the local capillary pressure evaluated at a positiony relative to the centroid of the large-scale averaging volume, and {p _c }¦_x represents the large-scale capillary pressure evaluated at the centroid.In addition to{p _c } ^c being evaluated at the centroid, all averaged terms on the right-hand side of Equation (1) are evaluated at the centroid. We can now write the equations describing the quasi-static condition as , , This means that the fluids within an averaging volume are distributed according to the capillary pressure-saturation relationwith the capillary pressure held constant. It also means that the large-scale capillary pressure is devoid of any dynamic effects. Both of these conditions represent approximations (see Section 6 in Part I) and one of our main objectives in this paper is to learn something about the efficacy of these approximations. As a secondary objective we want to explore the influence of dynamic effects in terms of our original theory. In that development only the first four terms on the right hand side of Equation (1) appeared in the representation for the local capillary pressure. However, those terms will provide an indication of the influence of dynamic effects on the large-scale capillary pressure and the large-scale permeability tensor, and that information provides valuable guidance for future studies based on the theory presented in Part I.Roman Letters A scalar that maps {}*/t onto - A scalar that maps {}*/t onto - A interfacial area between the -region and the -region contained within, m² - A interfacial area between the -region and the -region contained within, m² - A interfacial area between the -region and the -region contained within, m² - a vector that maps ({}*/t) onto , m - a vector that maps ({}*/t) onto , m - b vector that maps ({p}– g) onto , m - b vector that maps ({p}– g) onto , m - B second order tensor that maps ({p}– g) onto , m² - B second order tensor that maps ({p}– g) onto , m² - c vector that maps ({}*/t) onto , m - c vector that maps ({}*/t) onto , m - C second order tensor that maps ({}*/t) onto , m² - C second order tensor that maps ({}*/t) onto . m² - D third order tensor that maps ( ) onto , m - D third order tensor that maps ( ) onto , m - D second order tensor that maps ( ) onto , m² - D second order tensor that maps ( ) onto , m² - E third order tensor that maps () onto , m - E third order tensor that maps () onto , m - E second order tensor that maps () onto - E second order tensor that maps () onto - p _c =(), capillary pressure relationship in the-region - p _c =(), capillary pressure relationship in the-region - g gravitational vector, m/s² - largest of either or - - - i unit base vector in thex-direction - I unit tensor - K local volume-averaged-phase permeability, m² - K local volume-averaged-phase permeability in the-region, m² - K local volume-averaged-phase permeability in the-region, m² - {K } large-scale intrinsic phase average permeability for the-phase, m² - K –{K }, large-scale spatial deviation for the-phase permeability, m² - K –{K }, large-scale spatial deviation for the-phase permeability in the-region, m² - K –{K }, large-scale spatial deviation for the-phase permeability in the-region, m² - K ^* large-scale permeability for the-phase, m² - L characteristic length associated with local volume-averaged quantities, m - characteristic length associated with large-scale averaged quantities, m - I _i i = 1, 2, 3, lattice vectors for a unit cell, m - l characteristic length associated with the-region, m - ; characteristic length associated with the-region, m - l _H characteristic length associated with a local heterogeneity, m - - n unit normal vector pointing from the-region toward the-region (n =–n ) - n unit normal vector pointing from the-region toward the-region (n =–n ) - p pressure in the-phase, N/m² - p local volume-averaged intrinsic phase average pressure in the-phase, N/m² - {p } large-scale intrinsic phase average pressure in the capillary region of the-phase, N/m² - p local volume-averaged intrinsic phase average pressure for the-phase in the-region, N/m² - p local volume-averaged intrinsic phase average pressure for the-phase in the-region, N/m² - p –{p }, large scale spatial deviation for the-phase pressure, N/m² - p –{p }, large scale spatial deviation for the-phase pressure in the-region, N/m² - p –{p }, large scale spatial deviation for the-phase pressure in the-region, N/m² - P _c p –{p }, capillary pressure, N/m² - {p_c}^c large-scale capillary pressure, N/m² - r ₀ radius of the local averaging volume, m - R ₀ radius of the large-scale averaging volume, m - r position vector, m - , m - S /, local volume-averaged saturation for the-phase - S ^* {}*{}*, large-scale average saturation for the-phaset time, s - t time, s - u , m - U , m² - v -phase velocity vector, 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 } large-scale intrinsic phase average velocity for the-phase in the capillary region of the-phase, m/s - {v } large-scale phase average velocity for the-phase in the capillary region of the-phase, 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 local averaging volume, m³ - V volume of the-phase in, m³ - V large-scale averaging volume, m³ - V capillary region for the-phase within, m³ - V capillary region for the-phase within, m³ - V _c intersection of m³ - V volume of the-region within, m³ - V volume of the-region within, m³ - V () capillary region for the-phase within the-region, m³ - V () capillary region for the-phase within the-region, m³ - V () , region in which the-phase is trapped at the irreducible saturation, m³ - y position vector relative to the centroid of the large-scale averaging volume, m Greek Letters local volume-averaged porosity - local volume-averaged volume fraction for the-phase - local volume-averaged volume fraction for the-phase in the-region - local volume-averaged volume fraction for the-phase in the-region - local volume-averaged volume fraction for the-phase in the-region (This is directly related to the irreducible saturation.) - {} large-scale intrinsic phase average volume fraction for the-phase - {} large-scale phase average volume fraction for the-phase - {}* large-scale spatial average volume fraction for the-phase - –{}, large-scale spatial deviation for the-phase volume fraction - –{}, large-scale spatial deviation for the-phase volume fraction in the-region - –{}, large-scale spatial deviation for the-phase volume fraction in the-region - a generic local volume-averaged quantity associated with 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, N s/m² - interfacial tension of the - phase system, N/m - , N/m - , volume fraction of the-phase capillary (active) region - , volume fraction of the-phase capillary (active) region - , volume fraction of the-region ( + =1) - , volume fraction of the-region ( + =1) - {p }– g, N/m³ - {p }– g, N/m³     

Mixed convection flow in narrow vertical ducts

The mixed convection flow in a vertical duct is analysed under the assumption that , the ratio of the duct width to the length over which the wall is heated, is small. It is assumed that a fully developed Poiseuille flow has already been set up in the duct before heat from the wall causes this to be changed by the action of the buoyancy forces, as measured by a buoyancy parameter . An analytical solution is derived for the case when the Reynolds numberRe, based on the duct width, is of 0 (1). This is extended to the case whenRe is 0 (^–1) by numerical integrations of the governing equations for a range of values of representing both aiding and opposing flows. The limiting cases, || 1 andR=Re of 0 (1), andR and both large, with of 0 (R ^1/3) are considered further. Finally, the free convection limit, large with R of 0 (1), is discussed.

---

Mischkonvektion in engen senkrechten Rohren

Zusammenfassung Mischkonvektion in einem senkrechten Rohr wird unter der Voraussetzung untersucht, daß das Verhältnis der Rohrbreite zur Länge, über welche die Wand beheizt wird, klein ist. Es wird angenommen, daß sich bereits eine voll entwickelte Poiseuille-Strömung in dem Rohr eingestellt hat, bevor Antriebskräfte, gemessen mit dem Auftriebsparameter , aufgrund der Wandbeheizung die Strömung verändern. Es wird eine analytische Lösung für den Fall erhalten, daß die mit der Rohrbreite als charakteristische Länge gebildete Reynolds-ZahlRe konstant ist. Dies wird mittels einer numerischen Integration der wichtigsten Gleichungen auf den FallRe =f (^–1) sowohl für Gleich- als auch für Gegenstrom ausgedehnt. Weiterhin werden die beiden Grenzfälle betrachtet, wenn || 1 undR=Re konstant ist, sowieR und beide groß mit proportionalR ^1/3. Schließlich wird der Grenzfall der freien Konvektion, großes mit konstantem R, diskutiert.

---

Nomenclature g acceleration due to gravity - Gr Grashof number - G modified Grashof number - h duct width - l length of the heated section of the duct wall - p pressure - Pr Prandtl number - Q flow rate through the duct - Q ₀ heat transfer on the wally=0 - Q ₁ heat transfer on the wally=1 - Re Reynolds number - R modified Reynolds number - T temperature of the fluid - T ₀ ambient temperature - T applied temperature difference - u, velocity component in thex-direction - v, velocity component in they-direction - x, co-ordinate measuring distance along the duct - y, co-ordinate measuring distance across the duct - buoyancy parameter - ₀ modified buoyancy parameter, ₀=R ^–1/3 - coefficient of thermal expansion - ratio of duct width to heated length, =h/l - (non-dimensional) temperature - _w applied temperature on the wally=0 - kinematic viscosity - density of the fluid - ₀ shear stress on the wally=0 - ₁ shear stress on the wally=1 - stream function     

On a Boundary-Value Problem of Periodic Type for First-Order Linear Functional Differential Equations

We establish unimprovable sufficient conditions for the unique solvability of the boundary-value problem

Squeezing flows of Newtonian liquid films an analysis including fluid inertia

A theoretical study is made of the flow behavior of thin Newtonian liquid films being squeezed between two flat plates. Solutions to the problem are obtained by using a numerical method, which is found to be stable for all Reynolds numbers, aspect ratios, and grid sizes tested. Particular emphasis is placed on including in the analysis the inertial terms in the Navier-Stokes equations.Comparison of results from the numerical calculation with those from Ishizawa's perturbation solution is made. For the conditions considered here, it is found that the perturbation series is divergent, and that in general one must use a numerical technique to solve this problem.Nomenclature a half of the distance, or gap, between the two plates - a ₀ the value of a at time t=0 - adot da/dt - ä d² a/dt ² - d³ a/dt ³ - a ⁱ components of a contravariant acceleration vector - f unknown function of z ₀ and t defined in (6) - f ⁱ function defined in (9) f ¹=r ₀ g(z ₀, t) f ²= ₀ f ³=f(z ₀, t) - F force applied to the plates - g unknown function of z ₀ and t defined in (6) - g g/z ₀ - h grid dimension in the z ₀ direction (see Fig. 5) - Christoffel symbol - i, j, k, l indices - k grid dimension in the t direction (see Fig. 5) - r radial coordinate direction defined in Fig. 1 - r ₀ radial convected coordinate - R radius of the circular plates - t time - v _r fluid velocity in the r direction - v _z fluid velocity in the z direction - v fluid velocity in the direction - x ⁱ cylindrical coordinate x ¹=r x²= x³=z - z vertical coordinate direction defined in Fig. 1 - z ₀ vertical convected coordinate - tangential coordinate direction - ₀ tangential convected coordinate - viscosity - kinematic viscosity, / - ⁱ convected coordinate ¹=r₀ ²=₀ ³=z₀ - density     

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.     

On laminar flow through a uniformly porous pipe

Numerous investigations ([1] and [4–9]) have been made of laminar flow in a uniformly porous circular pipe with constant suction or injection applied at the wall. The object of this paper is to give a complete analysis of the numerical and theoretical solutions of this problem. It is shown that two solutions exist for all values of injection as well as the dual solutions for suction which had been noted by previous investigators. Analytical solutions are derived for large suction and injection; for large suction a viscous layer occurs at the wall while for large injection one solution has a viscous layer at the centre of the channel and the other has no viscous layer anywhere. Approximate analytic solutions are also given for small values of suction and injection.

Nomenclature

General r distance measured radially - z distance measured along axis of pipe - u velocity component in direction of z increasing - v velocity component in direction of r increasing - p pressure - density - coefficient of kinematic viscosity - a radius of pipe - V velocity of suction at the wall - r ²/a ² - R wall or suction Reynolds number, Va/ - f() similarity function defined in (6) - u ₀() eigensolution - U(0) a velocity at z=0 - K an arbitrary constant - B _K Bernoulli numbers Particular Section 5 perturbation parameter, –2/R - ² a constant, –K - x / - g(x) f()/ Section 6 perturbation parameter, –R/2 - ² a constant, –K - g() f() - g _c ()=g() near centre of pipe - * point where g()=0 Section 7 2/R - ² K - t (1–)/ - w(t, ) [1–f(t)]/ - ₀, ₁ constants - g() f()– ₀ - ₀/ - ₀ a constant - * point where f()=0     

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

Asymptotic behavior of solutions of nonlinear elliptic equations

We study and obtain formulas for the asymptotic behavior as ¦x¦ of C ² solutions of the semilinear equation u=f(x, u), x (*) where is the complement of some ball in ⁿ and f is continuous and nonlinear in u. If, for large x, f is nearly radially symmetric in x, we give conditions under which each positive solution of (*) is asymptotic, as ¦x¦, to some radially symmetric function. Our results can also be useful when f is only bounded above or below by a function which is radially symmetric in x or when the solution oscillates in sign. Examples when f has power-like growth or exponential growth in the variables x and u usefully illustrate our results.     

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     

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²     

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

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

Heteroclinic orbits in singular systems: A unifying approach

We consider singularly perturbed systems , such that=f(, _o, 0). _o ^m , has a heteroclinic orbitu(t). We construct a bifurcation functionG(, ) such that the singular system has a heteroclinic orbit if and only ifG(, )=0 has a solution=(). We also apply this result to recover some theorems that have been proved using different approaches.     

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