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     

Shear viscosity behavior of highly concentrated suspensions at low and high shear-rates

A method is proposed for determining the shear viscosity behavior of highly concentrated suspensions at low and high shear-rates through the use of a formulation that is a function of three parameters signifying the effects of particle size distribution. These parameters are the intrinsic viscosity [], a parametern that reflects the level of particle association at the initiation of motion and the maximum packing concentration _m. The formulation reduces to the modified Eilers equation withn = 2 for high shear rates. An analytical method was used for the calculation of maximum packing concentration which was subsequently correlated with the experimental values to account for the surface induced interaction of particles with the fluid. The calculated values of viscosities at low and high shear-rates were found to be in good agreement with various experimental data reported in literature. A brief discussion is also offered on the reliability of the methods of measuring the maximum packing concentration. _r = /₀ relative viscosity of the suspension - volumetric concentration of solids - k _n coefficient which characterizes a specific effect of particle interactions - _m maximum packing concentration - _r,0 relative viscosity at low shear-rates - [] intrinsic viscosity - n, n parameter that reflects the level of particle interactions at low and high shear-rates, respectively - _r, relative viscosity at high shear-rates - (_m)_s, (_m)_i, (_m)_l packing factors for small, intermediate and large diameter classes - v _s, v_i, v_l volume fractions of small, intermediate and large diameter classes, respectively - _si, _sl coefficient to be used in relating a smaller to an intermediate and larger particle group, respectively - _is, _il coefficient to be used in relating an intermediate to a smaller and larger particle group, respectively - _ls, _li coefficient to be used in relating a larger to a smaller and intermediate particle group, respectively - _m0 maximum packing concentration for binary mixtures - _m,e measured maximum packing concentration - _m,c calculated maximum packing concentration     

Laser velocimetry measurements in a horizontal gas-solid pipe flow

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

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

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

A light transmittance probe for interfacial area measurements of dispersed two phase flows with small particles

An optical probe measuring interfacial area () by light attenuation has been designed with a special emphasis on flows with sub-millimetric particles. It permits measurements in liquid-liquid or gas-liquid dispersions without need of introducing empirical correcting factors for the standard exponential decay law of light intensity while keeping an extended application range. This probe was successfully tested with an air-glass particle flow, the parameters of which were carefully determined basically by hold-up methods. The volume fraction of the dispersed phase was varied between 0.05% and 5%, and the particle size between 10 m and 300 m.List of symbols D diameter of spherical particle - D _S Sauter diameter - E ₀ irradiance on a surface perpendicular to light propagation 226E;=(1/l) averaged density function along y axis - f density function of a dispersion - f ₁, f ₂ focal length of the lenses L ₁, L ₂ - g granulometry function of a powder (probability density) - h granulometry function of a powder (unnormalized) - I ₀, I light beam intensity respectively before and after passing through the dispersion - j volumetric powder flow - K ₁, K ₂, K ₃ dimensionless constants - l optical path length of the beam in the dispersion - L experimental pipe width along x axis - m mass of a sample - n optical index of the continuous phase - p _a, p ₀ respectively slope of _a and ₀ straight line - r distance between particles - S _d scattering cross-section - V volume of dispersion - averaged particles velocity - x, y, z spatial coordinates - interfacial area - _a absolute interfacial area (by unit volume of dispersion) - ₀ interfacial area measured by light attenuation method - _d angle (around the initial direction of light propagation) within which a particle diffracts - _dr detector aperture angle - light wavelength - _d scattering cross section by unit volume of dispersion - light beam diameter - ₁, ₂ L₁, and L₂ lenses diameters - local volumetric fraction of dispersed phase - averaged fraction of dispersed phase along x axis - ₂ averaged fraction of dispersed phase along x and y axis - volumetric mass of particles     

Hydrodynamic models as applied to the investigation of magnetic plasma stability

In the present paper magnetohydrodynamic models are employed to investigate the stability of an inhomogeneous magnetic plasma with respect to perturbations in which the electric field may be regarded as a potential field (rot E 0). A hydrodynamic model, actually an extension of the well-known Chew-Goldberg er-Low model [1], is used to investigate motions transverse to a strong magnetic field in a collisionless plasma. The total viscous stress tensor is given; this includes, together with magnetic viscosity, the so-called inertial viscosity.Ordinary two-fluid hydrodynamics is used in the case of strong collisions=. It is shown that the collisional viscosity leads to flute-type instability in the case when, collisions being neglected, the flute mode is stabilized by a finite Larmor radius. A treatment is also given of the case when epithermal high-frequency oscillations (not leading immediately to anomalous diffusion) cause instability in the low-frequency (drift) oscillations in a manner similar to the collisional electron viscosity, leading to anomalous diffusion.Notation f particle distribution function - E electric field component - H₀ magnetic field - density - V particle velocity - e charge - m, M electron and ion mass - _i, _e ion and electron cyclotron frequencies - viscous stress tensor - P pressure - r_i Larmor radius - P pressure tensor - t time - frequency - T temperature - collision frequency - collision time - j current density - _i, _e ion and electron drift frequencies - k_x, k_y, k_z wave-vector components - n₀ particle density - g acceleration due to gravity. The authors are grateful to A. A. Galeev for valuable discussion.     

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     

Reservoir transport equations by compositional approach

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

Transport in ordered and disordered porous media V: Geometrical results for two-dimensional systems

In this paper we continue the geometrical studies of computer generated two-phase systems that were presented in Part IV. In order to reduce the computational time associated with the previous three-dimensional studies, the calculations presented in this work are restricted to two dimensions. This allows us to explore more thoroughly the influence of the size of the averaging volume and to learn something about the use of anon-representative region in the determination of averaged quantities.

Nomenclature

Roman Letters A interfacial area of the– interface associated with the local closure problem, m² - a _i i=1, 2, gaussian probability distribution used to locate the position of particles - l unit tensor - characteristic length for the-phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1,2,3 lattice vectors, m - m convolution product weighting function - m _V special convolution product weighting function associated with a unit cell - n _i i=1, 2 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume,V, m³ - x position of the centroid of an averaging area, m - x ₀ reference position of the centroid of an averaging area, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters V /V, volume average porosity - _{a i} standard deviation ofa _i - _r standard deviation ofr - intrinsic phase average of      

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)     

Effect of measuring volume length on two-component laser velocimeter measurements in a turbulent boundary layer

The effects of finite measuring volume length on laser velocimetry measurements of turbulent boundary layers were studied. Four different effective measuring volume lengths, ranging in spanwise extent from 7 to 44 viscous units, were used in a low Reynolds number (Re=1440) turbulent boundary layer with high data density. Reynolds shear stress profiles in the near-wall region show that u v strongly depends on the measuring volume length; at a given y-position, u v decreases with increasing measuring volume length. This dependence was attributed to simultaneous validations on the U and V channels of Doppler bursts coming from different particles within the measuring volume. Moments of the streamwise velocity showed a slight dependence on measuring volume length, indicating that spatial averaging effects well known for hot-films and hot-wires can occur in laser velocimetry measurements when the data density is high.List of symbols time-averaged quantity - u wall friction velocity, ( _w /)^1/2 - v kinematic viscosity - d _p pinhole diameter - l _eff spanwise extent of LDV measuring volume viewed by photomultiplier - l ⁺ non-dimensional length of measuring volume, l _eff u /v - y ⁺ non-dimensional coordinate in spanwise direction, y u /v - z ⁺ non-dimensional coordinate in spanwise direction, z u /v - U ⁺ non-dimensional mean velocity, /u - u instantaneous streamwise velocity fluctuation, U &#x2329;U - v instantaneous normal velocity fluctuation, V–V - u RMS streamwise velocity fluctuation, u ^21/2 - v RMS normal velocity fluctuation, v ^21/2 - Re Reynolds number based on momentum thickness, U 0/v - R _uv cross-correlation coefficient, u v/u v - R₁₂(0, 0, z) two point correlation between u and v with z-separation, <u(0, 0, 0) v (0, 0, z)>/<u(0, 0, 0) v (0, 0, 0)> - N rate at which bursts are validated by counter processor - T Taylor time microscale, u (dv/dt²)^–1/2     

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     

Transport in ordered and disordered porous media IV: Computer generated porous media for three-dimensional systems

In the method of volume averaging, the difference between ordered and disordered porous media appears at two distinct points in the analysis, i.e. in the process of spatial smoothing and in the closure problem. In theclosure problem, the use of spatially periodic boundary conditions isconsistent with ordered porous media and the fields under consideration when the length-scale constraint,r ₀L is satisfied. For disordered porous media, spatially periodic boundary conditions are an approximation in need of further study.In theprocess of spatial smoothing, average quantities must be removed from area and volume integrals in order to extractlocal transport equations fromnonlocal equations. This leads to a series of geometrical integrals that need to be evaluated. In Part II we indicated that these integrals were constants for ordered porous media provided that the weighting function used in the averaging process contained thecellular average. We also indicated that these integrals were constrained by certain order of magnitude estimates for disordered porous media. In this paper we verify these characteristics of the geometrical integrals, and we examine their values for pseudo-periodic and uniformly random systems through the use of computer generated porous media.

Nomenclature

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 _i i=1, 2, 3 gaussian probability distribution used to locate the position of particles - I unit tensor - L general characteristic length for volume averaged quantities, m - L characteristic length for , m - L characteristic length for , m - characteristic length for the -phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1, 2, 3 lattice vectors, m - m convolution product weighting function - m _v special convolution product weighting function associated with the traditional volume average - n _i i=1, 2, 3 integers used to locate the position of particles - 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 - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - r position vector, m - r _m support of the weighting functionm, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume,, m³ - x positional vector locating the centroid of an averaging volume, m - x ₀ reference position vector associated with the centroid of an averaging volume, m - y position vector locating points 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 - /L, small parameter in the method of spatial homogenization - standard deviation ofa _i - _r standard deviation ofr - _r intrinsic phase average of      

Nonlocal dispersion in media with continuously evolving scales of heterogeneity

General nonlocal diffusive and dispersive transport theories are derived from molecular hydrodynamics and associated theories of statistical mechanical correlation functions, using the memory function formalism and the projection operator method. Expansion approximations of a spatially and temporally nonlocal convective-dispersive equation are introduced to derive linearized inverse solutions for transport coefficients. The development is focused on deriving relations between the frequency-and wave-vector-dependent dispersion tensor and measurable quantities. The resulting theory is applicable to porous media of fractal character.Nomenclature C _v (t) particle velocity correlation function - C _v ,(t) particle fluctuation velocity correlation function - C _j (x,t) current correlation function - D(x,t) dispersion tensor - D(x,t) fluctuation dispersion tensor - f ₀(x,p) equilibrium phase probability distribution function - f(x, p;t) nonequilibrium phase probability distribution function - G(x,t) conditional probability per unit volume of finding a particle at (x,t) given it was located elsewhere initially - (k,t) Fourier transform ofG(x,t) - G(x,t) fluctuation conditional probability per unit volume of finding a particle at (x,t) given it was located elsewhere initially - k wave vector - K(t) memory function - L Liouville operator - m mass - p(t) particle momentum coordinate - P = (0)( , (0)) projection operator - Q =I-P projection operator - s real Laplace space variable - S(k, ) time-Fourier transform of(k,t) - t time - v(t) particle velocity vector - v(t) particle fluctuation velocity vector - V phase space velocity - time-Fourier variable - ^(itn)(k) frequency moment of(k,t) - x(t) particle displacement coordinate - x(t) particle displacement fluctuation coordinate - friction coefficient - (t) normalized correlation function General Functions () Dirac delta function - () Gamma function Averages ₀ Equilibrium phase-space average - Nonequilibrium phase-space average - (,) L ² inner product with respect tof ₀     

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²     

Transport in ordered and disordered porous media I: The cellular average and the use of weighting functions

In this work we consider transport in ordered and disordered porous media using singlephase flow in rigid porous mediaas an example. We defineorder anddisorder in terms of geometrical integrals that arise naturally in the method of volume averaging, and we show that dependent variables for ordered media must generally be defined in terms of thecellular average. The cellular average can be constructed by means of a weighting function, thus transport processes in both ordered and disordered media can be treated with a single theory based on weighted averages. Part I provides some basic ideas associated with ordered and disordered media, weighted averages, and the theory of distributions. In Part II a generalized averaging procedure is presented and in Part III the closure problem is developed and the theory is compared with experiment. Parts IV and V provide some geometrical results for computer generated porous media.Roman Letters A interfacial area of the- interface contained within the macroscopic region, m² - A_e area of entrances and exits for the-phase contained within the macroscopic system, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor, m² - L general characteristic length for volume averaged quantities, m - characteristic length (pore scale) for the-phase - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - N unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p0 reference pressure in the-phase, N/m² - p traditional intrinsic volume averaged pressure, N/m² - r0 radius of a spherical averaging volume, m - r position vector, m - r position vector locating points in the-phase, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - v velocity vector in the-phase, m/s - v traditional superficial volume averaged velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V/V, volume average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

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²     

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     

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     

Maximum density effects on thermal instability of horizontal laminar boundary layers

Linear stability theory is used to investigate the onset of longitudinal vortices in laminar boundary layers along horizontal semi-infinite flat plates heated or cooled isothermally from below by considering the density inversion effect for water using a cubic temperature-density relationship. The analysis employs non-parallel flow model incorporating the variation of the basic flow and temperature fields with the streamwise coordinate as well as the transverse velocity component in the disturbance equations. Numerical results for the critical Grashof number Gr _L ^* =Gr _X ^* /Re _X< ^{Emphasis>/3/2} are presented for thermal conditions corresponding to –0.5 ₁–2.0 and –0.8 ₂1.2.Nomenclature a wavenumber, 2/ - D operator, d/d - F (f–f)/2 - f dimensionless stream function - g gravitational acceleration - G eigenvalue, Gr _L/Re_L - Gr _L Grashof number based on L - Gr _X Grashof number based on X - L characteristic length, (X/U)^1/2 - M number of divisions in y direction - P pressure - Pr Prandtl number, / - p dimensionless pressure, P/( ² /Re _L) - Re _L, Re_X Reynolds numbers, (U L/)=Re _X< ^1/2 and (U), respectively - T temperature - U, V, W velocity components in X, Y, Z directions - u, v, w dimensionless perturbation velocities, (U, V, W)/U - X, Y, Z rectangular coordinates - x, y, z dimensionless coordinates, (X, Y, Z)/L - thermal diffusivity - coefficient of thermal expansion - ₁, ₂ temperature coefficients for density-temperature relationship - similarity variable, Y/L=y - dimensionless temperature disturbance, /T - dimensionless wavelength of vortex rolls, 2/a - ₁, ₂ thermal parameters defined by equation (12) - kinematic viscosity - density - dimensionless basic temperature, (T _b T )/T - –1 - T temperature difference, (T _w–T ) - * critical value or dimensionless disturbance amplitude - prime, disturbance quantity or differentiation with respect to - b basic flow quantity - max value at a density maximum - w value at wall - free stream condition