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     

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

Effects of MHD free convection and mass transfer on the flow past an oscillating infinite coaxial vertical circular cylinder

The effects of MHD free convection and mass transfer are taken into account on the flow past oscillating infinite coaxial vertical circular cylinder. The analytical expressions for velocity, temperature and concentration of the fluid are obtained by using perturbation technique.

---

Einwirkungen von freier MHD-Konvektion und Stoffübertragung auf eine Strömung nach einem schwingenden unendlichen koaxialen vertikalen Zylinder

Zusammenfassung Die Einwirkungen der freien MHD-Konvektion und Stoffübertragung auf eine Strömung nach einem schwingenden, unendlichen, koaxialen, vertikalen Zylinder wurden untersucht. Die analytischen Ausdrücke der Geschwindigkeit, Temperatur und Fluidkonzentration sind durch die Perturbationstechnik erhalten worden.

---

Nomenclature C _p Specific heat at constant temperature - C the species concentration near the circular cylinder - C _w the species concentration of the circular cylinder - C the species concentration of the fluid at infinite - * dimensionless species concentration - D chemical molecular diffusivity - g acceleration due to gravity - Gr Grashof number - Gm modified Grashof number - K thermal conductivity - Pr Prandtl number - r _a ,r _b radius of inner and outer cylinder - a, b dimensionless inner and outer radius - r,r coordinate and dimensionless coordinate normal to the circular cylinder - Sc Schmidt number - t time - t dimensionless time - T temperature of the fluid near the circular cylinder - T _w temperature of the circular cylinder - T temperature of the fluid at infinite - u velocity of the fluid - u dimensionless velocity of the fluid - U ₀ reference velocity - z,z coordinate and dimensionless coordinate along the circular cylinder - coefficient of volume expansion - * coefficient of thermal expansion with concentration - dimensionless temperature - H ₀ magnetic field intensity - coefficient of viscosity - _e permeability (magnetic) - kinematic viscosity - electric conductivity - density - M Hartmann number - dimensionless skin-friction - frequency - dimensionless frequency     

On the closure problem for Darcy's law

In a previous derivation of Darcy's law, the closure problem was presented in terms of an integro-differential equation for a second-order tensor. In this paper, we show that the closure problem can be transformed to a set of Stokes-like equations and we compare solutions of these equations with experimental data. The computational advantages of the transformed closure problem are considerable.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the averaging volume, m² - B second-order tensor used to respresent the velocity deviation - b vector used to represent the pressure deviation, m^–1 - C second-order tensor related to the permeability tensor, m^–2 - D second-order tensor used to represent the velocity deviation, m² - d vector used to represent the pressure deviation, m - g gravity vector, m/s² - I unit tensor - K C ^–1,–D, Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - l _i i=1, 2, 3, lattice vectors, m - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of the-phase - p pressure in the-phase, N/m ² - p intrinsic phase average pressure, N/m² - p –p , spatial deviation of the pressure in the-phase, N/m² - r position vector locating points in the-phase, m - r ₀ radius of the averaging volume, m - t time, s - v velocity vector in the-phase, m/s - v intrinsic phase average velocity in the-phase, m/s - v phase average or Darcy velocity in the \-phase, m/s - v –v , spatial deviation of the velocity in the-phase m/s - V averaging volume, m³ - V volume of the-phase contained in 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²     

Flow about a circular cylinder with a single large-scale surface perturbation

An experimental study of the flow around a cylinder with a single straight perturbation was conducted in a wind tunnel. With this bluff body, positioned in a uniform crossflow, the vortex shedding frequency and other flow characteristics could be manipulated.The Strouhal number has been shown to be a function of the perturbation angular position, _p , as well as the perturbation size and Reynolds number. As much as a 50% change in Strouhal number could be achieved, simply by changing _p by 1°. The perturbation size compared to the local boundary layer thickness, , was varied from approximately 1 to about 20 . The Reynolds number was varied from 10,000 to 40,000. For perturbation sizes approximately 5 to 20 and Reynolds numbers of 20,000 to 40,000, a consistent Strouhal number variation with _p was observed.A detailed investigation of the characteristic Strouhal number variation has shown that varying _p had a significant influence on the boundary layer separation and transition to turbulence. These significant changes occurring in the boundary layer have been shown to cause variations in the spacing between the shear layers, base pressure, drag, lift, and the longitudinal spacing between the vortices in the vortex street.List of Symbols a longitudinal spacing of vortices in the vortex street - C _d drag coefficient - C _dc drag coefficient corrected for blockage effect - C _l lift coefficient - C _p pressure coefficient, p _i –p /q - C _pb base pressure coefficient - C _pbc base pressure coefficient corrected for blockage effect - d perturbation diameter - d ^* spacing between the shear layers; defined as conditionally averaged spacing between points in the shear layers corresponding to 0.99u _max/U - D cylinder diameter; diameter of the circumscribing circle for a cable - f _v vortex shedding frequency - H wind tunnel test section cross-sectional width - L spanwise length of the cylinder - p _i tap pressure - p free stream static pressure - q free stream dynamic pressure - Re Reynolds number based on cylinder diameter - rms root-mean-square - S Strouhal number, f _v D/U - S _max maximum value of S - S _min minimum value of S - t time - u _c vortex convection velocity - u _max maximum velocity in the shear layer - U free stream velocity - U _c free stream velocity, corrected for blockage effect - x streamwise dimension referenced from the back of the cylinder - z lateral wake dimension, i.e., perpendicular to the free stream velocity vector and cylinder axis, referenced from the cylinder axis - x spacing between two hot wire probes aligned streamwise - phase difference between two hot wire probes aligned streamwise - boundary layer thickness - angle from stagnation point in degrees - _p perturbation angular position - _{b p} where S drops back to about the S of a cylinder - _c critical angle, angular position where S drops steeply with 1° change in - _{m p} where S was minimum - _{r p} after S recovers from drop in value - _{t p} where S starts to increase from about the S of a cylinder     

Heat conduction in a cylinder crossed by an electric current with skin effect

The steady periodic temperature distribution in an infinitely long solid cylinder crossed by an alternating current is evaluated. First, the time dependent and non-uniform power generated per unit volume by Joule effect within the cylinder is determined. Then, the dimensionless temperature distribution is obtained by analytical methods in steady periodic regime. Dimensionless tables which yield the amplitude and the phase of temperature oscillations both on the axis and on the surface of copper or nichrome cylindrical electric resistors are presented.

---

Wärmeleitung in einem stromdurchflossenen Zylinder unter Berücksichtigung des Skin-Effektes

Zusammenfassung Es wird die periodische Temperaturverteilung für den eingeschwungenen Zustand in einem unendlich langen, von Wechselstrom durchflossenen Vollzylinder ermittelt. Zuerst erfolgt die Bestimmung der zeitabhängigen, nichgleichförmigen Energiefreisetzung pro Volumeneinheit des Zylinders infolge Joulescher Wärmeentwicklung und anschließend die Ermittlung der quasistationären Temperaturverteilung auf analytischem Wege. Amplitude und Phasenverzögerung der Temperaturschwingungen werden für die Achse und die Oberfläche eines Kupfer- oder Nickelchromzylinders tabellarisch in dimensionsloser Form mitgeteilt.

---

Nomenclature A integration constant introduced in Eq. (2) - ber, bei Thomson functions of order zero - Bi Biot numberhr ₀/ - c speed of light in empty space - c ₁,c ₂ integration constants introduced in Eq. (46) - c _p specific heat at constant pressure - E electric field - E _z component ofE alongz - E time independent part ofE, defined in Eq. (1) - f function ofs and defined in Eq. (11) - g function ofs and defined in Eq. (37) - h convection heat transfer coefficient - H magnetic field - i imaginary uniti=(–1)^1/2 - I electric current - I _eff effective electric currentI _eff=I/2^1/2 - Im imaginary part of a complex number - J _n Bessel function of first kind and ordern - J electric current density - q _g power generated per unit volume - time average of the power generated per unit volume - time averaged power per unit length - r radial coordinate - R electric resistance per unit length - r ₀ radius of the cylinder - Re real part of a complex number - s dimensionless radial coordinates=r/r ₀ - s, s integration variables - t time - T temperature - time averaged temperature - T _f fluid temperature outside the boundary layer - time average of the surface temperature of the cylinder - u, functions ofs, and defined in Eqs. (47) and (48) - W Wronskian - x position vector - x real variable - Y _n Bessel function of second kind and ordern - z unit vector parallel to the axis of the cylinder - z axial coordinate - · modulus of a complex number - equal by definition Greek symbols amplitude of the dimensionless temperature oscillations - electric permittivity - dimensionless temperature defined in Eq. (16) - ₀, ₁, ₂ functions ofs defined in Eq. (22) - thermal conductivity - dimensionless parameter=(2)^1/2 - magnetic permeability - ₀ magnetic permeability of free space - function of defined in Eq. (59) - dimensionless parameter=c _p/() - mass density - electric conductivity - dimensionless time=t - phase of the dimensionless temperature oscillations - function ofs:= ₁+i ₂ - angular frequency - dimensionless parameter=()^1/2 r ₀     

The finite element solutions of laminar flow and combined convection of air in a staggered or an in-line tube-bank

A finite element method is used to solve the full Navier-Stokes and energy equations for the problems of laminar combined convection from three isothermal heat horizontal cylinders in staggered tube-bank and four isothermal heat horizontal cylinders in in-line tube-bank. The variations of surface shear stress, pressure and Nusselt number are obtained over the entire cylinder surface including the zone beyond the separation point. The predicted values of total, pressure and friction drag coefficients, average Nusselt number and the plots of velocity flow fields and isotherms are also presented.

---

Die Finite-Elemente-Lösung von laminarer Strömung und kombinierter Konvektion von Luft in einer versetzten oder fluchtenden Rohranordnung

Zusammenfassung Eine Methode der finiten Elemente wird zur Lösung der vollständigen Navier-Stokes- und der Energiegleichung für die Probleme der laminaren kombinierten Konvektion an drei isothermen geheizten horizontalen Zylindern in versetzter Rohranordnung sowie für vier isotherme geheizte horizontale Zylinder in fluchtender Anordnung verwendet.Die Veränderung der Wandschubspannung, des Druckes und der Nusselt-Zahl werden für die gesamte Zylinderoberfläche, einschließlich des Bereiches nach dem Ablösepunkt, bestimmt. Die Werte des gesamten Widerstandsbeiwertes aufgrund von Druck und Reibung, die durchschnittliche Nusselt-Zahl und die Diagramme des Geschwindigkeitsfeldes und der Isothermen werden ebenfalls aufgezeigt.

---

Nomenclature C specifie heat - C _D total drag coefficient - C _f friction drag coefficient - C _p pressure drag coefficient - D diameter of cylinder,L=2R ₀ - G, g gravitational acceleration - Gr Grashof number, g(T_w–T )D ³/v ² - h local heat transfer coefficient - K thermal conductivity - L spacing between the centers of cylinder - M _l shape function - N _i shape function - Nu, local and average Nusselt numbers - P dimensionless pressure, p^*/u ² - p ^*,p pressure, free stream pressure - Pe Peclet number,RePr - Pr Prandtl number, c/K - Ra Rayleigh number,Gr Pr - Re Reynolds number,Du /v - R ₀ radius of cylinder - T temperature - T _w temperature on cylinder surface with fixed value - T free stream temperature - v dimensionless x-direction component of velocity,v ^*/u - u ^* x-direction component of velocity - u free stream velocity - v dimensionless Y-direction component of velocity,v ^*/u - v ^* Y-direction component of velocity - X x-direction axis - x dimensionless x-direction coordinate,x ^*/D - x^* x-direction coordinate - Y Y-direction axis - y dimensionless Y-direction coordinate,y ^*/D - y ^* Y-direction coordinate Greek symbols coefficient of volumetric thermal expansion - plane angle - dynamic viscosity - kinematic viscosity, / - density of fluid - _w dimensionless surface shear stress, _* ^w /u ² - skw/* surface shear stress - dimensionless temperature,      

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³     

Nonstationary vibration of a rotating shaft with nonlinear spring characteristics during acceleration through a critical speed (A critical speed of a summed-and-differential harmonic oscillation)

Nonstationary vibration of a flexible rotating shaft with nonlinear spring characteristics during acceleration through a critical speed of a summed-and-differential harmonic oscillation was investigated. In numerical simulations, we investigated the influence of the angular acceleration , the initial angular position of the unbalance _n and the initial rotating speed on the maximum amplitude. We also performed experiments with various angular accelerations. The following results were obtained: (1) the maximum amplitude depends not only on but also on _n and : (2) when the initial angular position _n changes. the maximum amplitude varies between two values. The upper and lower bounds of the maximum amplitude do not change monotonously for the angular acceleration: (3) In order to always pass the critical speed with finite amplitude during acceleration. the value of must exceed a certain critical value.Nomenclature O-xyz rectangular coordinate system - , ₁, ₁ inclination angle of rotor and its projections to thexy- andyz-planes - I _r polar moment of inertia of rotor - I diametral moment of inertia of rotor - i _r ratio ofI _r toI - dynamic unbalance of rotor - directional angle of fromx-axis - c damping coefficient - spring constant of shaft - N _nt ,N _nt nonlinear terms in restoring forees in ₁ and ₁ directions - ₄ representative angle - a small quantity - V. V _u .V _N potential energy and its components corresponding to linear and nonlinear terms in the restoring forees - directional angle - _n coefficients of asymmetrical nonlinear terms - _n coefficients of symmetrical nonlinear terms - coefficients of asymmetrical nonlinear terms experessed in polar coordinates - coefficients of symmetrical nonlinear terms expressed in polar coordinates - rotating speed of shaft - t time - _n initial angular position of att=0 - p natural frequency - p ₁.p _t natural frequencies of forward and backward precessions - , ₁, ₁ total phases of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - , ₁, ₁ phases of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - P, R _t ,R _b amplitudes of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - difference between phases ( = _f – _u) - acceleration of rotor - initial rotating speed - t _t ,r _b amplitudes of nonstationary oscillation during acceleration - (r _t )_max, (r _b )_max maximum amplitudes of nonstationary oscillation during acceleration - (r ₁ ¹ )_max, (r _b ¹ )_max maximum value of angular acceleration of non-passable case - ₀ critical value over which the rotor can always pass the critical speed - p ₁,p ₂,p ₃,p ₄ natural frequencies of experimental apparatus     

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

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

On forced vibration of anisotropic cylinders

Summary The dynamic response of a circular cylinder with thick walls of transverse curvilinear isotropy subjected to a uniformly distributed pressure varying periodically with time is analyzed by means of the Laplace transformation, and the exact solution is obtained in closed form. The previously obtained solutions for forced vibrations with isotropy, and free vibrations with transverse curvilinear isotropy are included as special cases of the general results reported here.Nomenclature t time - r, , z cylindrical coordinates - _ii components of normal strain - _ii components of normal stress - u radial displacement - c _ij elastic constant - mass density - c ² c ₁₁/ - ² c ₂₂/c ₁₁ - a, b inner, outer radius of the cylinder - , A, B constants - forced angular frequency - function defined by (9) - p, real, complex variables - constant defined by (14) - real number - , Lamé elastic constants - J (x) Bessel function of first kind - Y (x) Bessel function of second kind - I (x) modified Bessel function of first kind - K (x) modified Bessel function of second kind     

Approximate solutions for drag coefficient of bubbles moving in shear-thinning elastic fluids

The drag coefficient for bubbles with mobile or immobile interface rising in shear-thinning elastic fluids described by an Ellis or a Carreau model is discussed. Approximate solutions based on linearization of the equations of motion are presented for the highly elastic region of flow. These solutions are in reasonably good agreement with the theoretical predictions based on variational principles and with published experimental data. C _D Drag coefficient - E ^* Differential operator [E ^{* 2} = ²/² + (sin/ ²)/(1/sin /)] - El Ellis number - F _D Drag force - K Consistency index in the power-law model for non-Newtonian fluid - n Flow behaviour index in the Carreau and power-law models - P Dimensionless pressure [=(p – p ₀)/₀ (U /R)] - p Pressure - R Bubble radius - Re ₀ Reynolds number [= 2R U /₀] - Re Reynolds number defined for the power-law fluid [= (2R) ⁿ U ^2–n /K] - r Spherical coordinate - t Time - U Terminal velocity of a bubble - u Velocity - Wi Weissenberg number - Ellis model parameter - Rate of deformation - Apparent viscosity - ₀ Zero shear rate viscosity - Infinite shear rate viscosity - Spherical coordinate - Parameter in the Carreau model - ^* Dimensionless time [=/(U /R)] - Dimensionless length [=r/R] - Second invariant of rate of deformation tensors - ^* Dimensionless second invariant of rate of deformation tensors [=/(U /R)²] - Second invariant of stress tensors - ^* Dimensionless second invariant of second invariant of stress tensor [= / ₀ ² (U /R)²] - Fluid density - Shear stress - ^* Dimensionless shear stress [=/ ₀ (U /R)] - _1/2 Ellis model parameter - ₁ ^2/* Dimensionless Ellis model parameter [= _1/2/ ₀(U /R)] - Stream function - ^* Dimensionless stream function [=/U R ²]     

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

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

Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, 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 a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

A three-parameter model describing the behaviour of a viscoelastic liquid in a tangential annular flow

A three-parameter model describing the shear rate-shear stress relation of viscoelastic liquids and in which each parameter has a physical significance, is applied to a tangential annular flow in order to calculate the velocity profile and the shear rate distribution. Experiments were carried out with a 5000 wppm aqueous solution of polyacrylamide and different types of rheometers. In a shear-rate range of seven decades (5 10^–3 s^–1 < < 1.2 10⁵ s^–1) a good agreement is obtained between apparent viscosities calculated with our model and those measured with three different types of rheometers, i.e. Couette rheometers, a cone-and-plate rheogoniometer and a capillary tube rheometer. a physical quantity defined by:a = {1 – ( / ₀)}/ ₀ (Pa^–1) - C constant of integration (1) - r distancer from the center (m) - r ₁,r ₂ radius of the inner and outer cylinder (m) - v _r local tangential velocity at a distancer from the center (v _r = _r r) (m s^–1) - v ₂ local tangential velocity at a distancer ₂ from the center (m s^–1) - shear rate (s^–1) - local shear rate (s^–1) - ₁ wall shear rate at the inner cylinder (s^–1) - dynamic viscosity (Pa s) - _a apparent viscosity (_a = / ) (Pa s) - _a1 apparent viscosity at the inner cylinder (Pa s) - ₀ zero-shear viscosity (Pa s) - infinite-shear viscosity (Pa s) - shear stress (Pa) - _r local shear stress at a distancer from the center (Pa) - ₀ yield stress (Pa) - ₁, ₂ wall shear-stress at the inner and outer cylinder (Pa) - _r local angular velocity (s^–1) - ₂ angular velocity of the outer cylinder (s^–1)     

The value of the critical exponent for reaction-diffusion equations in cones

Let D R ^N be a cone with vertex at the origin i.e., D = (0, )x where S ^N–1 and x D if and only if x = (r, ) with r=¦x¦, . We consider the initial boundary value problem: u _t = u+u ^p in D×(0, T), u=0 on Dx(0, T) with u(x, 0)=u ₀(x) 0. Let ₁ denote the smallest Dirichlet eigenvalue for the Laplace-Beltrami operator on and let ₊ denote the positive root of (+N–2) = ₁. Let p ^* = 1 + 2/(N + ₊). If 1 < p < p ^*, no positive global solution exists. If p>p ^*, positive global solutions do exist. Extensions are given to the same problem for u _t=+¦x¦ u ^p .This research was supported in part by the Air Force Office of Scientific Research under Grant # AFOSR 88-0031 and in part by NSF Grant DMS-8 822 788. The United States Government is authorized to reproduce and distribute reprints for governmental purposes not withstanding any copyright notation therein.     

Spectral and correlation characteristics of the turbulent boundary layer on an extended flexible cylinder

In marine geophysical seismological prospecting extensive use is made of towed receiving systems consisting of extended flexible cylinders containing acoustic sensors over which the water flows in the longitudinal direction. The boundary layer pressure fluctuations on the cylinder surface are picked up by the sensors as hydrodynamic noise. This paper is concerned with the study of the energy and spacetime characteristics of the pressure fluctuations in the turbulent boundary layer on an extended flexible cylinder in a longitudinal flow. The pressure fluctuations on the cylinder surface have been investigated experimentally for Re_X=(2–4)·10⁷, a dimensionless diameter of the pressure fluctuation sensors d _p ⁺ =d_pu^*/=500, where d_p is the sensor diameter, u^* the dynamic viscosity, and the viscosity coefficient, and frequencies 0.02311.259 (=^*/U). The spectral and correlation characteristics of the pressure fluctuations on the surface of the flexible cylinder are found to differ from the corresponding characteristics for a rigid cylinder [1–4].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i aza, No, 5, pp. 49–54, September–October, 1989.     

Transport in ordered and disordered porous media III: Closure and comparison between theory and experiment

In this paper we examine the closure problem associated with the volume averaged form of the Stokes equations presented in Part II. For both ordered and disordered porous media, we make use of a spatially periodic model of a porous medium. Under these circumstances the closure problem, in terms of theclosure variables, is independent of the weighting functions used in the spatial smoothing process. Comparison between theory and experiment suggests that the geometrical characteristics of the unit cell dominate the calculated value of the Darcy's law permeability tensor, whereas the periodic conditions required for thelocal form of the closure problem play only a minor role.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² - A interfacial area of the- interface associated with the local closure problem, m² - A _p surface area of a particle, m² - b vector used to represent the pressure deviation, m^–1 - B ⁰ B+I, a second order tensor that maps v _m ontov - B second-order tensor used to represent the velocity deviation - d _p 6V _p/A_p, effective particle diameter, m - d a vector related to the pressure, m - D a second-order tensor related to the velocity, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor calculated on the basis of a spatially periodic model, m² - 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 characteristic length for the volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - characteristic length (pore scale) for the-phase - _i i=1, 2, 3 lattice vectors, m - weighting function - m(-y) , convolution product weighting function - 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 - 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 _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function - r position vector, m - r position vector locating points in the-phase, m. - V averaging volume, m³ - B 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 _m superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - v traditional superficial volume averaged velocity, m/s - v – v _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²     

Take-off threshold velocity of solid particles lying under a turbulent boundary layer

The take-off of free solid particles by wind erosion has been investigated. A preliminary bibliographic study (Foucaut and Stanislas 1995) has enabled an analysis of the Bagnold criterion (1941) to be made and its accordance with the results obtained by White (1982) to be shown, thus leading to a semi-empirical model. This initial study led to a more judicious choice of the primary parameters, allowing a more physical representation of the threshold velocity. In the present study, the criterion validation was carried out in a specific boundary layer wind tunnel, by means of a direct measurement of the threshold velocity. The basic idea was to increase the wind tunnel velocity slowly and linearly and to perform an optical detection of the first take-offs.List of symbols D _p particle diameter - D _Pref reference diameter - _p D _p /D _pref - g gravity - H shape factor = (/) - h ⁺ roughness parameter ( =u D _p /2) - R _f =u D _p / - R _x u _e x/v - û u _e /v - u ⁺ flow velocity ( =u/u ) - u _e external velocity - u _elim external threshold velocity - û _e u _e /u _elim - u friction velocity - u ^* threshold friction velocity - u _tref reference velocity - ^* u ^*/u ^ref - y ⁺ yu / - _p ( _p –)/ - boundary layer thickness - displacement thickness - dynamic viscosity - v kinematic viscosity - momentum thickness - fluid density - _p particle density - standard deviation of particle diameter - shear stress Research carried out at Ecole des Mines de Douai