Nonlinear Taylor stability of viscoelastic fluids

Using Stuart's energy method, the torque on the inner cylinder, for a second order fluid, in the supercritical regime is calculated. It is found that when the second normal stress difference is negative, the flow is more stable than for a Newtonian fluid and the torque is reduced. If the second normal stress difference is positive, then the flow is more stable and there is no torque reduction. Experimental data related to the present work are discussed.Nomenclature a amplitude of the fundamentals - A _ij ⁽¹⁾ , A _ij ⁽²⁾ first and second Rivlin-Ericksen tensors - d r ₂–r ₁ - D d/dx - E - F - g _ij metric tensor - G torque on the inner cylinder in the supercritical regime - h height of the cylinders - k ₀ /d ² - k ₁ /d ² - I ₁ - I ₂ - I ₃ - I ₄ - r ₁, r ₂ radii of inner and outer cylinders respectively - r ₀ 1/2(r ₁+r ₂) - R Reynolds number ₁ r ₁ d/ ₀ - R _c critical Reynolds number - T Taylor number r _{1 1} ² d ^{3 2}/ ₀ ² *) - T _c critical Taylor number - u ₁, v ₁, w ₁ Fundamentals of the disturbance - u _i , v _i , w _i , (i>1) harmonics - mean velocity (not laminar velocity) - u –u ₁/ar _{1 1} - v v ₁/Rar _{1 1} - x (r–r ₀)/d - , material constants - 0 viscosity - wave number d - density - ₁ angular velocity of inner cylinder - tilde denotes complex conjugate     

Reduced and Generalized Stokes Resolvent Equations in Asymptotically Flat Layers, Part I: Unique Solvability

We study the generalized Stokes equations in asymptotically flat layers, which can be considered as compact perturbations of an infinite (flat) layer Besides standard non-slip boundary conditions, we consider a mixture of slip and non-slip boundary conditions on the upper and lower boundary, respectively. In the first part we prove the unique solvability in L^q-Sobolev spaces, 1 < q < , by extending the known results in the case of an infinite layer ₀ via a perturbation argument to asymptotically flat layers which are sufficiently close to ₀. Combining this result with standard cut-off techniques and the parametrix constructed in the second part, we prove the unique solvability for an arbitrary asymptotically flat layer. Moreover, we show equivalence of unique solvability of the generalized and the reduced Stokes resolvent equations, which is essential for the second part of this contribution.     

Flow through a helically coiled annulus

In this paper the flow is studied of an incompressible viscous fluid through a helically coiled annulus, the torsion of its centre line taken into account. It has been shown that the torsion affects the secondary flow and contributes to the azimuthal component of velocity around the centre line. The symmetry of the secondary flow streamlines in the absence of torsion, is destroyed in its presence. Some stream lines penetrate from the upper half to the lower half, and if is further increased, a complete circulation around the centre line is obtained at low values of for all Reynolds numbers for which the analysis of this paper is valid, being the ratio of the torsion of the centre line to its curvature.Nomenclature A =constant - a outer radius of the annulus - b unit binormal vector to C - C helical centre line of the pipe - D rL - g 1000 - K Dean number=Re² - L 1+r sin - M (L ²+ ² r ²)^1/2 - n unit normal vector to C - P, P pressure and nondimensional pressure - p ₀, p pressures of O(1) and O() - Re Reynolds number=aW ₀/ - (r, , s), (r, , s) coordinates and nondimensional coordinates - nonorthogonal unit vectors along the coordinate directions - r ₀ radius of the projection of C - t unit tangent vector to C - V _r, V , V _s velocity components along the nonorthogonal directions - V_r, V, V _s nondimensional velocity components along - W ₀ average velocity in a straight annulus Greek symbols , curvature and nondimensional curvature of C - U, V, W lowest order terms for small in the velocity components along the orthogonal directions t - _r, , _s first approximations to V _r , V, V _s for small - =/=/ - kinematic viscosity - density of the fluid - , torsion and nondimensional torsion of C - , stream function and nondimensional stream function - nondimensional streamfunction for U, V - a inner radius of the annulus After this paper was accepted for publication, a paper entitled On the low-Reynolds number flow in a helical pipe, by C.Y. Wang, has appeared in J. Fluid. Mech., Vol 108, 1981, pp. 185–194. The results in Wangs paper are particular cases of this paper for =0, and are also contained in [9].     

Flow of second-order fluids over an enclosed rotating disc

Summary This paper is devoted to a study of the flow of a second-order fluid (flowing with a small mass rate of symmetrical radial outflow m, taken negative for a net radial inflow) over a finite rotating disc enclosed within a coaxial cylinderical casing. The effects of the second-order terms are observed to depend upon two dimensionless parameters ₁ and ₂. Maximum values ₁ and ₂ of the dimensionless radial distances at which there is no recirculation, for the cases of net radial outflow (m>0) and net radial inflow (m<0) respectively, decrease with an increase in the second-order effects [represented by T(=₁+₂)]. The velocities at ₁ and ₂ as well as at some other fixed radii have been calculated for different T and the associated phenomena of no-recirculation/recirculation discussed. The change in flow phenomena due to a reversal of the direction of net radial flow has also been studied. The moment on the rotating disc increases with T.Nomenclature , , z coordinates in a cylindrical polar system - z ₀ distance between rotor and stator (gap length) - =/z ₀, dimensionless radial distance - =z/z ₀, dimensionless axial distance - _s = _s/z₀, dimensionless disc radius - V =(u, v, w), velocity vector - dimensionless velocity components - uniform angular velocity of the rotor - , p fluid density and pressure - P =p/(₂ z ₀₂ ² , dimensionless pressure - ₁, ₂, ₃ kinematic coefficients of Newtonian viscosity, elastico-viscosity and cross-viscosity respectively - ₁, _{2 2}/z ₀ ² , resp. ₃/z ₀ ² , dimensionless parameters representing the ratio of second-order and inertial effects - m = , mass rate of symmetrical radial outflow - l a number associated with induced circulatory flow - Rm =m/(z ₀₁), Reynolds number of radial outflow - R _l =l/(z ₀₁), Reynolds number of induced circulatory flow - Rz =z ₀ ² /₁, Reynolds number based on the gap - ₁, ₂ maximum radii at which there is no recirculation for the cases Rm>0 and Rm<0 respectively - ₁(T), ₂(T) ₁ and ₂ for different T - U _1(T) ⁽⁺⁾ = dimensionless radial velocity, Rm>0 - V _1(T) ⁽⁺⁾ = , dimensionless transverse velocity, Rm>0 - U _2(T) ^(–) = , dimensionless radial velocity, Rm=–Rn<0, m=–n - V _2(T) ^(–) = , dimensionless transverse velocity, Rm<0 - C _m moment coefficient     

About the Regularized Navier–Stokes Equations

The first goal of this paper is to study the large time behavior of solutions to the Cauchy problem for the 3-dimensional incompressible Navier–Stokes system. The Marcinkiewicz space L^3, is used to prove some asymptotic stability results for solutions with infinite energy. Next, this approach is applied to the analysis of two classical regularized Navier–Stokes systems. The first one was introduced by J. Leray and consists in mollifying the nonlinearity. The second one was proposed by J.-L. Lions, who added the artificial hyper-viscosity (–)^{/ 2}, > 2 to the model. It is shown in the present paper that, in the whole space, solutions to those modified models converge as t toward solutions of the original Navier–Stokes system.     

Momentum balance in highly distorted turbulent pipe flows

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

Reduced and Generalized Stokes Resolvent Equations in Asymptotically Flat Layers, Part II: H ∞-Calculus

We study the generalized Stokes equations in asymptotically flat layers, which can be considered as compact perturbations of an infinite (flat) layer Besides standard non-slip boundary conditions, we consider a mixture of slip and non-slip boundary conditions on the upper and lower boundary, respectively. In this second part, we use pseudodifferential operator techniques to construct a parametrix to the reduced Stokes equations, which solves the system in L^q-Sobolev spaces, 1 < q < , modulo terms which get arbitrary small for large resolvent parameters . This parametrix can be analyzed to prove the existence of a bounded H-calculus of the (reduced) Stokes operator.     

The virtual origin of riblets in an open channel flow

In this paper, a method using the mean velocity profiles for the buffer layer was developed for the estimation of the virtual origin over a riblets surface in an open channel flow. First, the standardized profiles of the mixing length were estimated from the velocity measurement in the inner layer, and the location of the edge of the viscous layer was obtained. Then, the virtual origins were estimated by the best match between the measured velocity profile and the equations of the velocity profile derived from the mixing length profiles. It was made clear that the virtual origin and the thickness of the viscous layer are the function of the roughness Reynolds number. The drag variation coincided well with other results.Nomenclature f _r skin friction coefficient - f _ro skin friction coefficient in smooth channel at the same flow quantity and the same energy slope - g gravity acceleration - H water depth from virtual origin to water surface - H ⁺ u*H/ - H false water depth from top of riblets to water surface - H ⁺ u*H/ - I _e streamwise energy slope - I _b bed slope - k riblet height - k ⁺ u*k/ - l mixing length - l _s standardized mixing length - Q flow quantity - Re Reynolds number volume flow/unit width/v - s riblet spacing - u mean velocity - u* friction velocity = - u* false friction velocity = - y distance from virtual origin - y distance from top of riblet - y ₀ distance from top of riblet to virtual origin - y _v distance from top of riblet to edge of viscous layer - y ⁺ u*y/ - y ⁺ u*y/ - y ₀ ⁺ u*y ₀/ - u ⁺ u*y/ - shifting coefficient for standardization - thickness of viscous layer=y ₀+y - ⁺ u*/ - ⁺ u*/ - eddy viscosity - ridge angle - v kinematic viscosity - density - shear stress     

Steady States of Anisotropic Generalized Newtonian Fluids

We consider the stationary flow of a generalized Newtonian fluid which is modelled by an anisotropic dissipative potential f. More precisely, we are looking for a solution of the following system of nonlinear partial differential equations

The Extinction Behavior of the Solutions for a Class of Reaction-Diffusion Equations

1　ＩｎｔｒｏｄｕｃｔｉｏｎａｎｄＬｅｍｍａｓＴｈｅｒｅａｒｅｍａｎｙｒｅｓｕｌｔｓａｂｏｕｔｅｘｉｓｔｅｎｃｅ (ｇｌｏｂａｌｏｒｌｏｃａｌ)ａｎｄａｓｙｍｐｔｏｔｉｃｂｅｈａｖｉｏｒｏｆｓｏｌｕｔｉｏｎｓｆｏｒｒｅａｃｔｉｏｎ＿ｄｉｆｆｕｓｉｏｎｅｑｕａｔｉｏｎｓ[1- 9].Ｂｙｔｈｅａｉｄｓｏｆｒｅｓｕｌｔｓ[2 ,3]ｏｆｅｑｕａｔｉｏｎ ｕ/ ｔ=Δｕ-λ｜ｕ｜γ- 1ｕｗｉｔｈｉｎｉｔｉａｌ＿ｂｏｕｎｄａｒｙｖａｌｕｅｓ,ｐａｐｅｒ [4 ]ｓｔｕｄｉｅｄｔｈｅｐｒｏｂｌｅｍｏｆ ｕ/ ｔ=Δｕ-λ｜ｅβｔｕ｜γ- …     

Inertial Manifolds on Squeezed Domains

Let be an arbitrary smooth bounded domain in and > 0 be arbitrary. Squeeze by the factor in the y-direction to obtain the squeezed domain = {(x,y)(x,y)}. In this paper we study the family of reaction-diffusion equations

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

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

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

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³     

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.     

Laminar source flow between parallel porous disks

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

Asymptotic analysis of linearly elastic shells: ‘Generalized membrane shells’

We consider a family of linearly elastic shells indexed by their half-thickness , all having the same middle surface % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaadofacqGH9aqpcqaHvpGAcaGGOaGafqyYdCNbaebacaGGPaaa% aa!4317!\[S = \varphi (\bar \omega )\], with % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiabew9aQjaacQdacuaHjpWDgaqeaiabgkOimlaadkfadaahaaWc% beqaaiaaikdaaaGccqGHsgIRcaWGsbWaaWbaaSqabeaacaaIZaaaaa% aa!4812!\[\varphi :\bar \omega \subset R^2 \to R^3 \], and clamped along a portion of their lateral face whose trace on S is % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiabew9aQjaacIcacqaHZoWzdaWgaaWcbaGaaGimaaqabaGccaGG% Paaaaa!41EB!\[\varphi (\gamma _0 )\], where % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaacIcacqaHZoWzdaWgaaWcbaGaaGimaaqabaGccaGGPaaaaa!401F!\[(\gamma _0 )\] is a fixed portion of with length % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaacIcacqaHZoWzdaWgaaWcbaGaaGimaaqabaGccaGGPaGaeyOp% a4JaaGimaaaa!41E1!\[(\gamma _0 ) > 0\]. Let % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaacIcacqaHZoWzdaWgaaWcbaGaeqySdeMaeqOSdigabeaakiaa% cIcacqaH3oaAcaGGPaGaaiykaaaa!45AA!\[(\gamma _{\alpha \beta } (\eta ))\] be the linearized strain tensor of S. We make an essential geometric and kinematic assumption, according to which the semi-norm % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaacYhacqGHflY1caGG8bWaa0baaSqaaiabeM8a3bqaaiaad2ea% aaaaaa!4345!\[| \cdot |_\omega ^M \] defined by % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaacYhacqaH3oaAcaGG8bWaa0baaSqaaiabeM8a3bqaaiaad2ea% aaGccqGH9aqpdaGadeqaamaaqababaGaaiiFaiaacYhaaSqaaiabeg% 7aHfrbbjxAHXgaiuaacaWFSaGaeqOSdigabeqdcqGHris5aOGaeq4S% dCMaeqySdeMaeqOSdiMaaiikaiabeE7aOjaacMcacaGG8bGaaiiFam% aaDaaaleaacaWGmbWaaWbaaWqabeaacaaIYaaaaSGaaiikaiabeM8a% 3jaacMcaaeaacaaIYaaaaaGccaGL7bGaayzFaaWaaWbaaSqabeaaca% aIXaGaai4laiaaikdaaaaaaa!61F1!\[|\eta |_\omega ^M = \left\{ {\sum\nolimits_{\alpha ,\beta } {||} \gamma \alpha \beta (\eta )||_{L^2 (\omega )}^2 } \right\}^{1/2} \] is a norm over the space % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaadAfacaGGOaGaeqyYdCNaaiykaiabg2da9maacmqabaGaeq4T% dGMaeyicI4SaamisamaaCaaaleqabaGaaGymaaaakiaacIcacqaHjp% WDcaGGPaGaai4oaiabeE7aOjabg2da9iaab+gacaqGUbGaeq4SdC2a% aSbaaSqaaiaabcdaaeqaaaGccaGL7bGaayzFaaaaaa!5361!\[V(\omega ) = \left\{ {\eta \in H^1 (\omega );\eta = {\text{on}}\gamma _{\text{0}} } \right\}\], excluding however the already analyzed membrane shells, where % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiabeo7aNnaaBaaaleaacaqGWaaabeaakiabg2da9iabgkGi2kab% eM8a3baa!42F8!\[\gamma _{\text{0}} = \partial \omega \] and S is elliptic. This new assumption is satisfied for instance if % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiabeo7aNnaaBaaaleaacaqGWaaabeaakiabgcMi5kabgkGi2kab% eM8a3baa!43B9!\[\gamma _{\text{0}} \ne \partial \omega \] and S is elliptic, or if S is a portion of a hyperboloid of revolution.We then show that, as 0, the averages % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaamaalaaabaGaaGymaaqaaiaaikdadaahaaWcbeqaaiabew7aLbaa% aaGcdaWdXaqaaiaadwhadaqhaaWcbaGaamyAaaqaaiabew7aLbaaki% aabsgacaWG4bWaa0baaSqaaiaaiodaaeaacqaH1oqzaaaabaGaeyOe% I0IaeqyTdugabaGaeqyTduganiabgUIiYdaaaa!4E28!\[\frac{1}{{2^\varepsilon }}\int_{ - \varepsilon }^\varepsilon {u_i^\varepsilon {\text{d}}x_3^\varepsilon } \] across the thickness of the shell of the covariant components u _i of the displacement of the points of the shell strongly converge in the completion V _M ^#() of V() with respect to the norm % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaacYhacqGHflY1caGG8bWaa0baaSqaaiabeM8a3bqaaiaad2ea% aaaaaa!4345!\[| \cdot |_\omega ^M \], toward the solution of a generalized membrane shell problem. This convergence result also justifies the recent formal asymptotic approach of D. Caillerie and E. Sanchez-Palencia.The limit problem found in this fashion is sensitive, according to the terminology recently introduced by J.L. Lions and E. Sanchez-Palencia, in the sense that it possesses two unusual features: it is posed in a space that is not necessarily contained in a space of distributions, and its solution is highly sensitive to arbitrarily small smooth perturbations of the data.Under the same assumption, we also show that the average % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaamaalaaabaGaaGymaaqaaiaaikdadaahaaWcbeqaaiabew7aLbaa% aaGcdaWdXaqaaiaadwhadaqhaaWcbaGaamyAaaqaaiabew7aLbaaki% aabsgacaWG4bWaa0baaSqaaiaaiodaaeaacqaH1oqzaaaabaGaeyOe% I0IaeqyTdugabaGaeqyTduganiabgUIiYdaaaa!4E28!\[\frac{1}{{2^\varepsilon }}\int_{ - \varepsilon }^\varepsilon {u_i^\varepsilon {\text{d}}x_3^\varepsilon } \] where % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaadwhadaahaaWcbeqaaiabew7aLbaakiabg2da9iaacIcacaWG% 1bWaa0baaSqaaiaadMgaaeaacqaH1oqzaaGccaGGPaaaaa!452C!\[u^\varepsilon = (u_i^\varepsilon )\], and the solution % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiabe67a4naaCaaaleqabaGaeqyTdugaaOGaeyicI4SaamOvamaa% BaaaleaacaWGlbaabeaakiaacIcacqaHjpWDcaGGPaaaaa!465B!\[\xi ^\varepsilon \in V_K (\omega )\] of Koiter's equations have the same principal part as 0 in the same space V _M () as above. For such generalized membrane shells, the two-dimensional shell model of W.T. Koiter is thus likewise justified.We also treat the case where % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaacYhacqGHflY1caGG8bWaa0baaSqaaiabeM8a3bqaaiaad2ea% aaaaaa!4345!\[| \cdot |_\omega ^M \] is no longer a norm over V(), but is a norm over the space % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaadAfadaWgaaWcbaGaam4saaqabaGccaGGOaGaeqyYdCNaaiyk% aiabg2da9maacmqabaGaeq4TdGMaeyypa0JaaiikaiabeE7aOnaaBa% aaleaacaWGPbaabeaakiaacMcacqGHiiIZcaWGibWaaWbaaSqabeaa% caaIXaaaaOGaaiikaiabeM8a3jaacMcacqGHxdaTcaWGibWaaWbaaS% qabeaacaaIYaaaaOGaaiikaiabeM8a3jaacMcacaGG7aGaeq4TdG2a% aSbaaSqaaiaadMgaaeqaaOGaeyypa0JaeyOaIy7aaSbaaSqaaiaadA% haaeqaaOGaeq4TdG2aaSbaaSqaaiaaiodaaeqaaOGaeyypa0JaaGim% aiGac+gacaGGUbGaeq4SdC2aaSbaaSqaaiaaicdaaeqaaaGccaGL7b% GaayzFaaaaaa!68B8!\[V_K (\omega ) = \left\{ {\eta = (\eta _i ) \in H^1 (\omega ) \times H^2 (\omega );\eta _i = \partial _v \eta _3 = 0\operatorname{on} \gamma _0 } \right\}\], thus also excluding the already analyzed flexural shells. Then a convergence theorem can still be established, but only in the completion of the quotient space V()/V ₀() with repect to % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaacYhacqGHflY1caGG8bWaa0baaSqaaiabeM8a3bqaaiaad2ea% aaaaaa!4345!\[| \cdot |_\omega ^M \], where % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaadAfadaWgaaWcbaGaaGimaaqabaGccaGGOaGaeqyYdCNaaiyk% aiabg2da9maacmqabaGaeq4TdGMaeyicI4SaamOvaiaacIcacqaHjp% WDcaGGPaGaai4oaiabeo7aNjabeg7aHjabek7aIjaacIcacqaH3oaA% caGGPaGaeyypa0JaaeimaiaabMgacaqGUbGaeqyYdChacaGL7bGaay% zFaaaaaa!5997!\[V_0 (\omega ) = \left\{ {\eta \in V(\omega );\gamma \alpha \beta (\eta ) = {\text{0in}}\omega } \right\}\].These convergence results, together with those that we already obtained for membrane and flexural shells, jointly with B. Miara in the second case, thus constitute an asymptotic analysis of linearly elastic shells in all possible cases.     

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.     

Symplectic approach to the theory of quantized fields. II

We prove that the set D of vector fields on the configuration space B of a field whose 1-parameter groups locally associated are groups of fibre-preserving transformations of B that leave invariant that field in the sense of variational theory, is a Lie algebra with respect to ordinary addition, multiplication by real numbers and Lie brackets. We see that this Lie algebra structure can be carried over to the corresponding set of Noether invariants, which then becomes a Lie algebra in a natural way.Further, we define the n-form of Poincaré-Cartan of a field, and we use it to generalize the Lie algebras D and in a reasonable way. The algebras D and are subalgebras of the new Lie algebras D and introduced. A main result in this connection is the following: the differential d of the n-form of Poincaré-Cartan is –(d+f), where (, d+f) are the field equations on the vertical bundle B.The symplectic manifold of solutions associated with a field is introduced in a formal way and the former Lie algebras D, , D, are interpreted on this manifold. In imitation of the case of analytical dynamics, the main results in this direction are: a) Every vector field of the Lie algebra D defines, in a canonical way, a vector field on the manifold of solutions such that its polar 1-form with respect to the symplectic metric ₂ is the differential of its corresponding Noether invariant, and b) the Lie bracket [, ] of two Noether invariants , is the Noether invariant given by ₂(D, D), where D, D are the vector fields on the manifold of solutions defined, in the sense a), by two infinitesimal generators of , , respectively. This will allow us to regard the Lie algebra as the analogous object in field theory to the Poisson algebra of analytic dynamics.We apply the general formalism to the relativistic theory of non-linear scalar fields, and we compare our results with the formalism developed by I. Segal for this case.     

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