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³     

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

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

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

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

Oscillatory free convection heat transfer along a semi-infinite vertical plate

Summary The fluctuating free convection flow along a semi-infinite vertical plate is considered when the plate temperature is of the form T _p –T =(T ₀ –T ) where 0 < 1, denotes the frequency of oscillation and the mean temperature T ₀–T is proportional to ⁿ (0 n < 1). Flow and temperature fields have been obtained by means of two asymptotic expansions. For small values of the frequency parameter , a regular expansion is obtained while for large the method of matched asymptotic expansion is used. It is found that the skin friction and the rate of heat transfer obtained from two expansions overlap satisfactorily for a certain value of . For n=1 the flow governing equations to a semisimilar form, and have been solved by finite difference method. The results obtained from the series and the finite difference methods are in good agreement.

---

Oszillierender Wärmeübergang an einer halbunendlichen senkrechten Platte bei freier Konvektion

Übersicht Betrachtet wird die fluktuierende freie Konvektionsströmung längs einer halbunendlichen senkrechten Platte, deren Temperatur dem Gesetz T _p –T =(T ₀–T ) [1+ sin {ie1-03}] folgt, wobei 0 < 1 gelte, {ie1-04} die Frequenz ist und der Temperatur-Mittelwert T ₀–T proportional zu ⁿ (0 n < 1) ist. Mit Hilfe zweier asymptotischer Entwicklungen werden die Strömungs- und Temperaturfelder gewonnen. Für kleine Werte des Frequenzparameters wird eine gewöhnliche Entwicklung benutzt, für große die Methode angepaßter asymptotischer Entwicklungen. Es stellt sich heraus, daß die Oberflächenreibung und die Wärmeübergangsrate aus zwei Entwicklungen für ein bestimmtes zufriedenstellend aufeinander fallen. Für n=1 werden die Grundgleichungen zueinander ähnlich und werden nach der Finite-Differenzen-Methode gelöst. Die Ergebnisse nach den Reihenentwicklungen und der Finite-Differenzen-Methode stimmen gut überein.

     

On the boundary conditions at the macroscopic level

We study the problem of the boundary conditions specified at the boundary of a porous domain in order to solve the macroscopic transfer equations obtained by means of the volume-averaging method. The analysis is limited to the case of conductive transport but the method can be extended to other cases. A numerical study enables us to illustrate the theoretical results in the case of a model porous medium. Roman Letters _sf interfacial area of the s-f interface contained within the macroscopic system m² - A _sf interfacial area of the s-f interface contained within the averaging volume m² - C _p mass fraction weighted heat capacity, kcal/kg/K - d _s , d _f microscopic characteristic length m - g vector that maps to _s, m - h vector that maps to _f , m - K _eff effective thermal conductivity tensor, kcal/m s K - l REV characteristic length, m - L macroscopic characteristic length, m - n _fs outwardly directed unit normal vector for the f-phase at the f-s interface - n _e outwardly directed unit normal vector at the dividing surface - T ^* macroscopic temperature field obtained by solving the macroscopic equation (3), K - V averaging volume, m³ - V _s , V _f volume of the considered phase within the averaging volume, m³ - volume of the macroscopic system, m³ - _s , _f volume of the considered phase within the volume of the macroscopic system, m³ - dividing surface, m² Greek Letters _s , _f volume fraction - ratio of thermal conductivities - _s , _f thermal conductivities, kcal/m s K - spatial average density, kg/m³ - microscopic temperature, K - ^* microscopic temperature corresponding to T ^* , K - spatial deviation temperature K - error on the temperature due to the macroscopic boundary conditions, K - spatial average - ^s , ^f intrinsic phase average     

Chaos and integrability in a nonlinear wave equation

We consider the parametrized family of equations _tt ,u- _xx u-au+u ₂ ² u=O,x(0,L), with Dirichlet boundary conditions. This equation has finite-dimensional invariant manifolds of solutions. Studying the reduced equation to a four-dimensional manifold, we prove the existence of transversal homoclinic orbits to periodic solutions and of invariant sets with chaotic dynamics, provided that =2, 3, 4,.... For =1 we prove the existence of infinitely many first integrals pairwise in involution.     

New calculations on the Faraday effect in wave guides

Summary The propagation of electromagnetic waves is investigated theoretically for a round wave guide, containing a gyroelectrie-gyromagnetic medium with gyroaxis parallel to the guide in the form of a cylindrical shell of thickness, adjacent to the wall of the guide. An equation is set up, permitting to compute the change in the propagation constant due to the presence of the shell, including terms proportional to ². Assuming only the presence of gyromagnetism, the change ₁ of first order in for TE-waves is determined and is found to be the same fpr right- and left-circular polarization. The second order difference ₂ ⁺ — ₂ ^- for the two senses of polarization, however, appears to have a non-vanishing value which, just like ₁ can be expressed in terms of the radius of the guide, the frequency, the dielectric constant and the elements of the gyromagnetic permeability tensor which characterize the medium of the shell.     

Thin-Walled Beams: the Case of the Rectangular Cross-Section

In this paper we present an asymptotic analysis of the three-dimensional problem for a thin linearly elastic cantilever =×(0,l) with rectangular cross-section of sides and ², as goes to zero. Under suitable assumptions on the given loads, we show that the three-dimensional problem converges in a variational sense to the classical one-dimensional model for extension, flexure and torsion of thin-walled beams. Mathematics Subject Classifications (2000) 474K20, 74B10, 49J45.     

Über Differentialgleichungen für Wärme- und Stoffaustausch von hygroskopischen Stoffen

Zusammenfassung Die beiden Differentialgleichungssysteme vonKrischer undLykow werden miteinander verglichen. Dabei ergibt sich, daß die in der deutschen und russischen Literatur angewandten mathematischen Modelle der Trocknung von kapillarporösen Körpern praktisch übereinstimmen. Es werden die Transformationsgleichungen der dimensionslosen Kenngrößen angegeben, die die Beziehungen zwischen den beiden Systemen herstellen.

---

The differential equations ofKrischer andLuikow for unsteady internal heat and mass transfer in the porous medium are compared. It is shown, that the mathematical models for drying in the German and Russian literature are equivalent. The transform relations of the non-dimensional parameters between the two models are given.

---

Formelzeichen nach Krischer z laufende Koordinate in Strömungsrichtung in m - R kennzeichnende Abmessung des Körpers in m - t Zeit in h - Raumgewicht bei mittlerer Feuchtigkeit in kg/m³ - _w Teilgewicht des Wassers in 1 m³ Trockengut in kg/m³ - _wa Anfangsfeuchtigkeit in kg/m³ - _D Dampfdichte in kg/m³ - _L Luftvolumen je m³ Trocknungsgut in m³/m³ - Temperatur in °C - _u Umgebungstemperatur in °C - _a Anfangstemperatur in °C - r Verdampfungswärme in kcal/kg - q _E Wärmeentwicklung in kcal/m³ h - c spezifische Wärmekapazität des Trockengutes in kcal/kg grd - Wärmeleitfähigkeit in kcal/m h grd - Feuchtigkeitsleitzahl des Trockengutes in m²/h - wirksame Diffusionszahl von Wasserdampf in Luft in m²/h - Diffusionswiderstandszahl des Trockengutes — - Konstante — - Konstante in kg/m³ grd Formelzeichen nach Lykow X=r/R dimensionslose Koordinate des Körpers;r laufende Koordinate in m;R kennzeichnende Abmessung in m; - Fo=a/R ² Fourier-Zahl;a Temperaturleitzahl in m²/h; Zeit in h - T(X, Fo)=t(r, )– 0/t dimensionslose Temperatur des Körpers im Punkt mit KoordinateX für den ZeitpunktFo;t(r, ) Temperatur in °C; 0 mittlere Anfangstemperatur in °C; t ein vorher angenommener Temperaturunterschied in grd - (X, Fo)= ₀–u(r, )/u dimensionsloses Potential des Stoffübergangs im Punkt mit KoordinateX für den ZeitpunktFo;u(r, ) Feuchtigkeitsgehalt des Trockengutes in kg/kg; ₀ mittlerer Anfangsfeuchtigkeitsgehalt in kg/kg; u ein vorher angenommener Unterschied des Feuchtigkeitsgehalts in kg/kg - Ko= u/c t Kosowitsch-Zahl; Verdampfungswärme in kcal/kg;c spezifische Wärmekapazität des Trockengutes in kcal/kg - Ko*=Ko modifizierte Kosowitsch Zhal; Kenngröße der Zustandsänderung - Pn= t/u Posnowsche Zahl;=a _T ^m /a _m Thermogradientkoeffizient in 1/grd;a _T ^m thermische Stoffübergangszahl (charakterisiert den Stoffstrom unter der Einwirkung von Temperaturgradienten) in m²/h grd;a _m Stoffübergangszahl (charakterisiert den Stoffstrom unter der Einwirkung von Feuchtigkeitsgradienten) in m²/h - Lu=a _m/a Lykowsche Zahl - Ki _q=q _q ()·R/ _q t dimensionsloser Wärmestrom (Kirpitschew-Zahl);q _q() Wärmestrom durch die Körperoberfläche in kcal/m²; _q Wärmeleitfähigkeit in kcal/m² h grd - Ki _m=q _m ()·R/a _{m 0} u dimensionsloser Stoffstrom;q _m() Stoffstrom durch die Körperoberfläche in kg/m² h; ₀ Wichte des Trockengutes in kg/m³     

Existence and Uniqueness of Time-Periodic Physically Reasonable Navier-Stokes Flow Past a Body

Let be a three-dimensional exterior domain of class C^2,, 0<<1. Assume that a Navier-Stokes liquid is moving in under the action of a body force F that is time-periodic of period T, and that the velocity of the liquid is zero at spatial infinity. In this paper we show that, if F satisfies suitable conditions, and its norm, in appropriate function spaces, is sufficiently small, there is at least one time-periodic strong solution. Furthermore, the velocity field v of such a solution decays to zero for large |x| as |x|^–1 and its spatial gradient decays as |x|^–2, both uniformly in time. In addition, the pressure p decays like |x|^–2 and its gradient like |x|^–3, for almost all t[0,T]. In the special case where F is time-independent, these solutions are also time-independent and coincide with Finns physically reasonable solutions [4]. Moreover, we show that our time-periodic solutions are unique in a very large class, namely, the class of time-periodic weak solutions satisfying the energy inequality and with corresponding pressure fields verifying mild summability conditions in ×[0,T].     

Remarks on the bifurcation of solutions of a nonlinear eigenvalue problem

Conclusions We have investigated solutions of equation (3) when ² is an eigenvalue of the linearized operator (13) and when it is not. In Section 4 we have shown that for 0 and ² = _i ² we have exactly two nontrivial solutions which bifurcate to the right of _i ² ; these solutions are shown to exist in an interval ( _i ² , _i ² + ₀). The method of Section 3 may then be used to extend these two solutions to the right of _i ² + ₀ providing that ²= _i ² + ₀ is not an eigenvalue of the linear operator (13) evaluated at = ± ₁. Either a solution can be uniquely extended, or there exists a value of ²where the bifurcation method must be applied again³.While the method used here gives the exact number of solutions bifurcating from _i ² , other problems remain open; for example, it is still not proven that the two bifurcating branches have i zeros, as is the case for Hammerstein operators with oscillation kernels [4]. The conjecture of Odeh and Tadjbakhsh that there are exactly 2(i+1) nontrivial solutions in the interval _i ² < _i ^+1/2 remains un-answered, although it would be proven if one could show that there is no secondary bifurcation as in the cases of Kolodner [7] and Coffman [8].     

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.     

Production of temperature fluctuations in grid turbulence: Wiskind's experiment revisited

Measurements have been made in nearly-isotropic grid turbulence on which is superimposed a linearly-varying transverse temperature distribution. The mean-square temperature fluctuations, , increase indefinitely with streamwise distance, in accordance with theoretical predictions, and consistent with an excess of production over dissipation some 50% greater than values recorded in previous experiments. This high level of production has the effect of reducing the ratio,r, of the time scales of the fluctuating velocity and temperature fields. The results have been used to estimate the coefficient,C, in Monin's return-to-isotropy model for the slow part of the pressure terms in the temperature-flux equations. An empirical expression by Shih and Lumley is consistent with the results of earlier experiments in whichr 1.5, C 3.0, but not with the present data where r 0.5, C 1.6. Monin's model is improved when it incorporates both time scales.List of symbols C coefficient in Monin model, Eq. (5) - M grid mesh length - m exponent in power law for temperature variance, x ^m - n turbulence-energy decay exponent,q ² x ^-n - p production rate of - p pressure - q ² - R microscale Reynolds number - r time-scale ratiot/t - T mean temperature - U mean velocity - mean-square velocity fluctuations (turbulent energy components) - turbulent temperature flux - x, y, z spatial coordinates - temperature gradient dT/dy - thermal diffusivity - dissipation rate ofq ²/2 - dissipation rate of - Taylor microscale (₂=5q₂/) - temperature microscale - _v temperature-flux correlation coefficient, /v - dimensionless distance from the grid,x/M     

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

Zur analytischen Beschreibung des Füllstoffeinflusses auf das Fließverhalten von Kunststoff-Schmelzen

From literature some representative equations have been compiled describing the influence of filler on the viscosity of polymer melts. By application of these on the experimental results obtained from GF-SAN it was found that the relative viscosity _R , i.e. the ratio of the viscosities of the filled and unfilled melt, shows a pronounced dependence on the shear rate but not on the shear stress. Defining _R with constant and not with constant (as it is usually done), an analytical approach is possible independent of Further the influence of pressure, temperature and filler content on the zero-shear viscosity of filled polymer melts may be expressed by a modified Arrhenius equation.

     

Existence theorem for a minimum problem with free discontinuity set

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

Surface wave profile measurement by image analysis

A way of measuring the geometrical characteristics of progressive steep water surface waves is to use a visualization technique connected with image analysis processing. In the laboratory, visualization of wave profiles can be realized with quite simple procedures: a previous paper (Bonmarin and Ramamonjiarisoa 1985) has described a technique allowing such a visualization in a large water tank 40 m long, 3.2 m wide and about 1 m deep. This paper reported also on a manual process for analysing the wave pictures obtained. In the present paper, we describe an automated image analysis method which is complementary to the manual process mentioned above. It uses a video technique and allows analysis of a large number of pictures leading to statistical measurements.List of symbols L total wave length - H total wave height - crest elevation above still water level - trough depression below still water level - wave steepness = H/L - crest steepness = /L - F ₁ forward horizontal length from zero-upcross point (A) to wave crest - F ₂ backward horizontal length from wave crest to zero-downcross point (B) - crest front steepness = /F ₁ - crest rear steepness = /F ₂ - vertical asymmetry factor = F ₂/F ₁ (describing the wave's asymmetry with respect to a vertical axis through the wave crest) - horizontal asymmetry factor = /H (describing the wave's asymmetry with respect to a horizontal axis: SWL) - L ₃ vertical asymmetry factor = L ₂/L ₁ (describing the asymmetry between the crest and the trough) - E _p potential energy of the wave - e ₊ ratio between the potential energy located in the crest and the total potential energy of the wave     

An experimental study of the lateral migration of a droplet in a creeping flow

The distribution of droplets in a plane Hagen-Poiseuille flow of dilute suspensions has been measured by a special LDA technique. This method assumes a well defined relation between the velocity of the droplets and their lateral position in the channel. The measurements have shown that the droplet distribution is non-uniform and depends on the viscosity ratio between the droplets and the carrier liquid. The results have been compared with a theory by Chan and Leal describing the lateral migration of suspended droplets.List of symbols a particle radius, m - d half width of the channel, m - Re flow Reynolds number, = 2 _m · d · /µ - flow velocity, m/s - _m flow velocity at the channel axis, m/s - We Weber number, = 2 _m ^Emphasis>/2 · d · / - x distance from center line (x = 0) of the channel, m - non-dimensional distance from the channel center line, _x ^d - y distance along the channel (y = 0 at channel inlet), m - non-dimensional distance along the channel, = y/2d - non-dimensional, normalized distance along the channel, = · _m · µ/ - interfacial tension, N/m - viscosity ratio of dispersed (droplet) phase to viscosity of continuous phase - µ viscosity of continuous phase, Pa · s - density of continuous phase, kg/m³ - phase density difference, kg/m³ Experiments were performed at Max-Planck-Institut, Göttingen     

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

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