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     

Stability characteristics of condensate films

The linear stability theory is used to study stability characteristics of laminar condensate film flow down an arbitrarily inclined wall. A critical Reynolds number exists above which disturbances will be amplified. The magnitude of the critical Reynolds number is in all practical situations so small that a laminar gravity-induced condensate film can be expected to be unstable. Several stabilizing effects are acting on the film flow; at an inclined wall these effects are due to surface tension, gravity and condensation mass transfer.

---

Zusammenfassung Mit Hilfe der linearen Stabilitätstheorie werden die Stabilitätseigenschaften laminarer Kondensatfilme an einer geneigten Wand untersucht. Es zeigt sich, daß Kondensatfilme in jedem praktischen Fall ein unstabiles Verhalten aufweisen. Der stabilisierende Einfluß von Oberflächenspannung, Schwerkraft und Stoffübertragung durch Kondensation bewkkt jedoch, daß Störungen in bestimmten Wellenlängenbereichen gedämpft werden.

---

Nomenclature c=c^*/u₀ complex wave velocity, celerity, dimensionless - c^*=c _r ^* + i c _i ^* complex wave velocity, celerity, dimensional - c_p specific heat at constant pressure - g gravitational acceleration - h_fg latent heat - k thermal conductivity of liquid - p^* pressure - p=p^*/u₀ ² dimensionless pressure - Pe=Pr Re^* Peclet number - Pr Prandtl number - Re^*=u₀ / Reynolds number (defined with surface velocity) - S temperature perturbation amplitude - t^* time - t=t^* u₀/ dimensionless time - T temperature - T_s saturation temperature - T_w wall temperature - T=T_s-T_w temperature drop across liquid film - u^*, v^* velocity components - u=u^*/u₀ dimensionless velocity components - v=v^*/u₀ dimensionless velocity components - u₀ surface velocity of undisturbed film flow - v _g ^* vapor velocity - x^*, y^* coordinates - x=x^*/ dimensionless coordinates - y=y^*/ dimensionless coordinates Greek Symbols =^* wave number, dimensionless - ^*=2 /^* wave number dimensional - ^* wave length, dimensional - =^*/ wave length, dimensionless - local thickness of undisturbed condensate film - kinematic viscosity, liquid - density, liquid - _g density vapor - surface tension - = (1 +) film thickness of disturbed film, Fig. 1 - stream function perturbation amplitude - angle of inclination Base flow quantities are denoted by^–, disturbance quantities are denoted by.     

Simultaneous concentration and velocity measurements using combined laser-induced fluorescence and laser Doppler velocimetry: Application to turbulent transport

This paper describes the implementation of an optical technique, allowing to perform concentration and velocity measurements simultaneously and at the same point. This method is based on the coupling of laser-induced fluorescence of rhodamine B, applied to the determination of local concentration, and laser Doppler velocimetry. The method developed provides an accurate measurement of the concentration-velocity cross-correlation. The latter is a parameter linked to the eddy diffusivity tensor of a passive contaminant. This method was tested with a turbulent submerged free jet and it allowed the determination of the mean field of concentration and velocity, the concentration-velocity cross-correlation, and the local eddy diffusivity.List of symbols C molar concentration - c fluctuating part of the concentration - mean value of the concentration - concentration-velocity cross-correlation - D molecular diffusivity - (D _eddy) _ij eddy diffusivity tensor - I _abs absorbed intensity - I _e local incident intensity - K _opt optical constant - N number of samples - r _c half-width radius for the concentration profile - r _v half-width radius for the velocity profile - S _f fluorescence signal - Sc Schmidt number - V _c collection volume - U velocity - U _e flow velocity in the channel - U _i injection velocity Greek symbols kinetic energy dissipation rate - ₁ molar extinction coefficient for the laser radiation (in m²mol^–1l^–1) - ₂ molar extinction coefficient for the fluorescence signal (in m²mol^–1l^–1) - quantum yield - _c Batchelor scale - _k Kolmogorov scale - v kinematic viscosity - normalized values     

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

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

Das erweiterte Korrespondenzgesetz für die dynamische Idealviskosität und die Idealwärmeleitfähigkeit reiner Stoffe

Zusammenfassung Zur Berechnung der dynamischen Idealviskosität Ideal (T) und der Idealwärmeleitfähigkeit ideal (T) benötigt man die kritische TemperaturT _kr, das kritische spezifische Volum _kr, die MolmasseM, den kritischen Parameter _kr und die molare isochore WärmekapazitätC _v(T). Sowohl das theoretisch, als auch das empirisch abgeleitete erweiterte Korrespondenzgesetz ergeben eine für praktische Zwecke ausreichende Genauigkeit für die Meßwertwiedergabe, die bei den assoziierenden Stoffen und den Quantenstoffen jedoch geringer ist als bei den Normalstoffen.

---

The extended correspondence law for the ideal dynamic viscosity and the ideal thermal conductivity of pure substances

For the calculation of the ideal dynamic viscosity Ideal (T) and the ideal thermal conductivity ideal (T) the critical temperatureT _kr, the critical specific volumev _kr, the molecular massM, the critical parameter _kr, and the molar isochoric heat capacityC _v(T) is needed. Not only the theoretically determined but also the empirically determined extended correspondence law gives for practical use a good representation of the measured data, which for the associating substances and the quantum substances is not so good as for the normal substances.

     

Existence theorem for a minimum problem with free discontinuity set

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

A solitary wave in a porous medium

The one-phase Darcy continuity equation, including the quadratic gradient term, is considered. The exact linearization of the equation is found by a functional transformation for an arbitrary spatial dimension in the limit case where the constant fluid compressibility is much more dominant than the constant compressibilities of the reservoir parameters.The equation permits a solution representing a localized wave travelling through a one-dimensional reservoir without changing its form. This is the actual long-time limit of the transient solution for a constant sandface-rate injection of a compressible fluid with a constant compressibility if the fluid is much more compressible than the matrix. A solitary wave solution is not possible for production.A fully developed solitary wave would appear only for very high pressure increases, but the first signs of the emerging solitary wave are detectable at the sandface for moderate pressure increases which can appear under physical reservoir conditions.Latin symbols a Dimensionless wave propagation velocity - A _N Sandface area (N = 0, 1, 2) - c ₁, c ₂ Sums of compressibilities - c _x Generic (generalized) compressibility - c Fluid compressibility - c _h Reservoir height (i.e. bulk volume) compressibility (N = 0, 1) - c _k , c , c Generalized compressibilities - D Spatial reservoir dimensionality (D = 1, 2, 3) - f Fractional change of p _n1 due to nonlinear effects - h Reservoir height (proportional to bulk volume for N = 0, 1) - Horizontal reservoir width (N = 0) - k Reservoir permeability - K _N Constant with dimension of pressure (N = 0, 1, 2) - n Sum index - N Integer variable (N = D – 1) - p Reservoir pressure - p^* Overburden pressure - p _D Dimensionless (scaled) version of p - p ₀ Initial pressure - q Volumetric flow rate referred to sandface - r Radial (or linear) spatial distance from center of well - r _w Well radius - r _e External reservoir radius (or length) from center of well - t Time variable - t _f Injection/production time corresponding to fraction f - T Cole-Hopf-transformed version of dimensionless pressure y - u Rescaled (dimensionless) version of v _D - v Darcy velocity - v _d Dimensionless (scaled) version of v - x Generic symbol in compressibility expression (also used for auxiliary function and for auxiliary variable) - y Rescaled (dimensionless) version of p _D - z Dimensionless (scaled) version of r Greek symbols Coefficient of inertial resistance - Variable in wave solution for y - p _n1 Absolute change in physical sandface pressure due to production or injection - _p Pressure change over (dimensionless) distance behind and far away from front - _r Physical distance at constant time corresponding to - Characteristic (dimensionless) width of solitary wave - Formation porosity - ₁, ₂ Integration constants - Dimensionless (scaled) length of finite reservoir - Fluid viscosity - Fluid density - Dimensionless (scaled) version of t - Wave solution for dimensionless pressure y - Integer variable (±1) distinguishing between production and injection     

Measurements of thermally stratified pipe flow using image-processing techniques

The cross-correlation technique and Laser Induced Fluorescence (LIF) have been adopted to measure the time-dependent and two-dimensional velocity and temperature fields of a stably thermal-stratified pipe flow. One thousand instantaneous and simultaneous velocity and temperature maps were obtained at overall Richardson numberRi = 0 and 2.5, from which two-dimensional vorticity, Reynolds stress and turbulent heat flux vector were evaluated. The quasi-periodic inclined vortices (which connected to the crest) were revealed from successive instantaneous maps and temporal variation of vorticity and temperature. It has been recognized that these vortices are associated with the crest and valley in the roll-up motion.List of symbols A Fraction of the available light collected - C Concentration of fluorescence - D Pipe diameter - I Fluorescence intensity - L Sampling length along the incident beam - I ₀ Intensity of an excitation beam - I _c (T) Calibration curve between temperature and fluorescence intensity - I _ref Reference intensity of fluorescence radiation - Re _b Reynolds number based on bulk velocity,U _b D/v - Ri Overall Richardson number based on velocity difference,gDT/U ² - t Time - t Time interval between the reference and corresponding matrix - T Temperature - T ₁,T ₂ Temperature of lower and upper layer - T ^* Normalized temperature, (T–T ₁)/T - T _c (I) Inverse function of temperature as a function ofI _c - T _ref Reference temperature - T Temperature difference between upper and lower flow,T ₂ –T ₁ - U ₁ Velocity of lower stream - U ₂ Velocity of upper stream - U _b Bulk velocity - U _c Streamwise mean velocity atY/D=0 - U Streamwise velocity difference between upper and lower flow,U ₁–U ₂ - u, v, T Fluctuating component ofU, V, T - U, V Velocity component of X, Y direction - X Streamwise distance from the splitter plate - Y Transverse distance from the centerline of the pipe - Z Spanwise distance from the centerline of the pipe - Quantum yield - Absorptivity - vorticity calculated from a circulation - Kinematic viscosity - circulation     

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     

Finite-difference analysis of MHD Stokes problem for a vertical infinite plate in dissipative fluid

Finite-difference solution of MHD flow past an impulsively started vertical infinite plate in an electrically conducting fluid has been presented on taking into account the viscous dissipative heat. Results for velocity and temperature are shown graphically whereas the numerical values of the skin-friction and the rate of heat transfer are entered in the table. The results are discussed in terms of the parameters M (the Hartmann number), G (the Grashof number, G>0, cooling of the plate by free convection, G<0, heating of the plate by free convection currents), E (the Eckert number) and P (the Prandtl number).Nomenclature B ₀ applied magnetic field - c _p specific heat at constant pressure - g acceleration due to gravity - k thermal conductivity - t time - T temperature of the fluid near the plate - T temperature of the fluid far away from the plate - U ₀ velocity of the plate - u velocity of the fluid - coefficient of volume expansion - kinematic viscosity - scalar electrical conductivity - coefficient of viscosity - density of the fluid     

Boundary layer diagnostics by means of an infrared scanning radiometer

A computerized infrared (IR) scanning radiometer is employed to characterize the boundary layer development over a model wing, having a Göttingen 797 cross-section, by measuring the temperature distribution over its heated surface. The Reynolds analogy is used to relate heat transfer measurements to skin friction. The results show that IR thermography is capable of rapidly detecting location and extent of transition and separation regions of the boundary layer over the whole surface of the tested model wing. Thus, the IR technique appears to be a suitable and effective diagnostic tool for aerodynamic research in wind tunnels.List of symbols c airfoil chord - c _f local skin friction coefficient = 2/( V ²) - c _p specific heat coefficient at constant pressure - h local convective heat transfer coefficient - Nu Nusselt number = h x/ - Nu _c Nusselt number based on airfoil chord = h c/ - Pr Prandtl number c _p / - Q _j wall heat flux due to Joule heating - Q _l heat flux loss - Re Reynolds number V x/ - Re _c Reynolds number based on airfoil chord = V c/ - St Stanton number = h/c _p V - T _w wall temperature - T _aw adiabatic wall temperature - V velocity of the free stream - x chordwise spatial coordinate - angle of attack - thermal conductivity coefficient - dynamic viscosity coefficient - mass density - wall shear stress     

Hydrodynamical changes due to polymer migration in very dilute solutions

The slip hypothesis, based on thermodynamical arguments, has been extended to obtain the flow characteristics of polymer solutions flowing in a nonhomogeneous flow field. An asymptotic analysis, valid for both channel and falling film flows, is presented that predicts the flow enhancement due to polymer migration. Concentration-viscosity coupling is shown to be a critical factor in the hydrodynamic analysis. The analysis, which essentially provides an upper bound on flow enhancement, explicitly accounts for the influence of wall shear stress, initial polymer concentration etc. A comparison with the pertinent experimental data shows reasonable agreement. c concentration - c ₀ concentration in shear-free region - c _i initial concentration - d rate of deformation tensor - g acceleration due to gravity - g ₁ function defined in eq. [13] or [15] - g ₂ function defined in eq. [18] or [20] - H half-channel thickness or film thickness - K gas law constant - L length of the channel or film - q flow rate per unit width - q ^* normalized flow rate - T temperature - v velocity - V mean velocity - y transverse distance - y _c location of solvent layer - _w /µ _s - _w /c ₀ KT - /t convected derivative - dimensionless cenentration,c/c ₀ - _c dimensionless interface concentration - _w dimensionless wall concentration - relaxation time - µ _eff effective viscosity - µ _s solvent viscosity - dimensionless transverse distance,y/H - _c dimensionless interface location - density - stress tensor - _w wall shear stress - c _i KT/ _w - ns no slip NCL-Communication No. 3155     

The motion of a detaching elastic body

Conclusions The qualitative behavior of the displacement (t) and the radius R(t) during the different phases of the motion is illustrated in the diagram of Fig. 6.1.After the first impact at t = 0 the displacement (t) varies according to (5.2). If the first maximum of (t) is higher than ₁ then at time t ₁ the graph of (t) intersects the straight line = _cand detachment first occurs. In the second phase the dependance of on t is expressed by (5.6). The detachment will end at the instant t ₂ when vanishes.The radius R remains equal to R ₀ until (t) reaches the critical value ₁ = _c that is at t = t ₁. After t ₁, R(t) will decrease according to (4.4) up to its final value ₂.A rather unexpected property of the solution is that the greatest elongation is finite for every non-vanishing value of the ratio .To Jerry Ericksen for his 60^th birthday     

Double diffusive convection in the diffusive regime

Numerical solutions to convection of a fluid which is heated, and to which salt is introduced from the bottom, have been obtained. Although different in boundary conditions from the conditions in the study of Huppert and Moore, qualitatively the flow investigated here has many features the same as theirs. The differences are discussed and solutions are given for two fluids whose Prandtl numbers are 1 and 7 and ratio of molecular and thermal diffusivities is 0.1. The fields of temperature and salinity, and the stream patterns are shown in contour plots. The saline and thermal Nusselt numbers are given as functions of the thermal Rayleigh number.Nomenclature A aspect ratio of the enclosure - g gravitational acceleration - H the height of the enclosure - k _S molecular diffusivity of salt - k _T molecular diffusivity of heat - N _S averaged saline Nusselt number - N _T averaged thermal Nusselt number - P transported variable of temperature, salinity, or vorticity - p pressure above its hydrostatic value - Pr Prandtl number - velocity vector - R _S saline Rayleigh number - R _T thermal Rayleigh number - S salinity - S slainity difference between the top and the bottom - S ₀ salinity at the top surface - S _r salinity at the reference state=S₀+S/2 - T temperature - T temperature difference between the top and and the bottom - T ₀ temperature at the top surface - T _r temperature at the reference state=T₀+T/2 - t time - u velocity component in the horizontal direction - transport velocity - v velocity component in the vertical direction - W the width of the enclosure - x horizontal coordinate - x mesh size in the horizontal direction - y vertical coordinate Greek Letters coefficient of thermal expansion - coefficient of volumetric expansion produced by salt - kinematic viscosity - ratio of diffusivities - vorticity - density - _r density at the reference state - density difference between the top and the bottom - streamfunction     

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     

Effects of viscous dissipation on heat transfer in an oscillating flow past a flat plate

Theoretical investigation has been carried out of laminar thermal boundary layer response to harmonic oscillations in velocity associated with a progressive wave imposed on a steady free stream velocity and convected in the free stream direction. Series solutions are derived both to velocity and temperature field and the resulting equations are solved numerically. The functions affecting the temperature field are shown graphically for different values of Prandtl number. It is observed that there is more reduction in the rate of heat transfer for P _r<1 and a rise in the rate of heat transfer for P _r>1 due to the presence of oscillatory free-stream.Nomenclature u, v velocity components in the x and y direction - x, y Cartesian coordinate axes - t time - U, U ₀ instantaneous value of and mean free stream velocity - density of fluid - kinematic viscosity - T, T _w, T temperature of the fluid, wall and free stream fluid - c _p specific heat at constant pressure - thermal diffusivity - amplitude of free stream velocity - frequency - _p non-dimensional temperature (T–T /T _w–T ) - P _r Prandtl number (c _p/K) - E _c Eckert number (U ₀ ² /c _p(T _w–T )) - a parameter ( ) - ₀ boundary layer thickness of the oscillation of a harmonic oscillation of frequency ( ) - ordinary boundary layer thickness ( ) - time-averaged, time-independent external velocity - A, B, C, D, E, K, L, M, N, P functions used in expansion for u and - Nu Nusselt number (hx/k) - T _w–% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8V4rqqrFfpeea0Jc9yq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepGe9fr-xfr-x% frpeWZqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaacIcadaGcaa% qaaiaadAhacaWG4bGaai4laiqadwfagaqeaaWcbeaakiaacMcaaaa!3CA6!\[(\sqrt {vx/\bar U} )\] - k thermal conductivity     

Slit die viscometry at shear rates up to 5 × 106s−1: An analytical correction for small viscous heating errors

If the viscosity can be expressed in the form = (T)f(), the walls are at a constant temperatureT ₀, and the extra stress, velocity and temperature fields are fully developed, then the wall shear rate can be calculated by applying the Weissenberg-Rabinowitsch operator toF _c Q instead of to the flow rateQ, whereF _c is a correction factor which differs from 1 when the temperature field is non-uniform; the isothermal equation relating the wall shear stress and pressure gradient is still valid. For the case in whcih = ₀|| ⁿ /(1 +(T–T ₀)), wheren, ₀, and are independent of shear stress and temperatureT, an exact analytical expression forF _c in terms of the Nahme-Griffith numberNa andn is obtained. Use of this expression gives agreement with data obtained for degassed decalin ( = 2.5 mPa s) from a new miniature slit-die viscometer at shear rates up to 5 × 10⁶s^–1; here, the correction is only 7%,Na is 1.3, andGz, the Graetz number at the die exit, is 119. For a Cannon standard liquidS6 ( = 9 mPa s), agreement extends up to 5 × 10⁵s^–1; at 2×10⁶s^–1 (whereNa = 7.2 andGz = 231), the corrections are 11% (measured) and 36% (calculated).Notation x, y Cartesian coordinates - v _x ,v velocity inx-direction, dimensionless velocity - p _xx ,p _yy normal stress onx- andy-planes - N ₁ first normal stress difference - shear stress ony-planes acting inx-direction - _w value of shear stress at the wall - shear rate, shear rate at the wall - Q, Q flow rate (Eqs. (2.13), (2.15)) - T, T ₀ temperature, temperature at the wall - ø, dimensionless temperature (Eqs. (2.24), (2.25)) - h, w half of die height, width of die - R diameter of a tube - , ₀ viscosity, viscosity atT = T ₀ - viscosity-temperature coefficient - k thermal conductivity - c _p specific heat at constant pressure - n, m dimensionless parameters characterizing shear stress dependence of viscosity - Na Nahme Griffith number (Eq. (2.21)) - Gz Graetz number (Eq. (5.1)) - F _c viscous heating correction factor (Eq. (2.18)) - ( ) a function characterizing temperature dependence of viscosity (Eq. (2.8)) - J _k ( ) Bessel function of the first kind This paper is dedicated to Professor Hanswalter Giesekus on the occasion of his retirement as Editor of Rheologica Acta.     

Stress analysis of the cylinder subjected to arbitrarily distributed axisymmetrical loads

Summary Stress analysis has been carried out for a finite cylinder subjected to arbitrarily distributed axisymmetrical surface loads. Direct stress _x in the axial direction is assumed to be of the form _x = ₀+r ₁ +r ₂ where ₀ to ₂ are functions of x. Using the equations of equilibrium and compatibility the other direct stresses and the shearing stress are expressed by ₁ and ₂. Fundamental equations governing ₁ and ₂ are introduced using the variational principle of complementary energy. From the results of the present analysis it is evident that the boundary conditions can be satisfied completely even for the case where the external forces are specified in complicated form, and that more accurate solutions can easily be obtained by introducing additional terms in _x.

---

Spannungsanalyse für den Zylinder unter axialsymmetrischer Last in beliebiger Verteilung

Übersicht Für einen endlichen Zylinder unter axialsymmetrischer Oberflächenlast in beliebiger Verteilung werden die Spannungen ermittelt. Die Normalspannung in Axialrichtung wird in der Form _x = ₀+r ₁ +r ₂ angesetzt mit ₀, ₁, ₂ als Funktionen von x. Mit Hilfe der Gleichgewichtsund Verträglichkeitsbedingungen werden die anderen Normalspannungen und die Schubspannung durch ₁ und ₂ ausgedrückt. Über das Variationsprinzip für die Komplementärenergie werden die grundlegenden Gleichungen für ₁ und ₂ eingeführt. Die Ergebnisse zeigen, daß die Randbedingungen selbst für komplizierte Belastungsarten vollständig erfüllbar sind und mit zusätzlichen Termen in _x mühelos noch genauere Lösungen bestimmt werden können.

     

Deformation to breaking of deep water gravity waves

The results of laboratory observations of the deformation of deep water gravity waves leading to wave breaking are reported. The specially developed visualization technique which was used is described. A preliminary analysis of the results has led to similar conclusions than recently developed theories. As a main fact, the observed wave breaking appears as the result of, first, a modulational instability which causes the local wave steepness to approach a maximum and, second, a rapidly growing instability leading directly to the breaking.List of symbols L total wave length - H total wave height - crest elevation above still water level - trough depression below still water level - wave steepness =H/L - crest steepness =/L - trough steepness =/L - F ₁ forward horizontal length from zero-upcross point (A) to wave crest - F ₂ backward horizontal length from wave crest to zero-downcross point (B) - crest front steepness =/F ₁ - crest rear steepness =/F ₂ - vertical asymmetry factor=F ₂/F ₁ (describing the wave asymmetry with respect to a vertical axis through the wave crest) - µ horizontal asymmetry factor=/H (describing the wave asymmetry with respect to a horizontal axis: SWL) - T ₀ wavemaker period - L ₀ theoretical wave length of a small amplitude sinusoïdal wave generated at T _inf0 ^sup–1 frequency - ₀ average wave height     

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³