Impacts of local dispersion and first-order decay on solute transport in randomly heterogeneous porous media

Stochastic subsurface transport theories either disregard local dispersion or take it to be constant. We offer an alternative Eulerian-Lagrangian formalism to account for both local dispersion and first-order mass removal (due to radioactive decay or biodegradation). It rests on a decomposition of the velocityv into a field-scale componentv , which is defined on the scale of measurement support, and a zero mean sub-field-scale componentv _s , which fluctuates randomly on scales smaller than. Without loss of generality, we work formally with unconditional statistics ofv _s and conditional statistics ofv . We then require that, within this (or other selected) working framework,v _s andv be mutually uncorrelated. This holds whenever the correlation scale ofv is large in comparison to that ofv _s . The formalism leads to an integro-differential equation for the conditional mean total concentration c which includes two dispersion terms, one field-scale and one sub-field-scale. It also leads to explicit expressions for conditional second moments of concentration cc. We solve the former, and evaluate the latter, for mildly fluctuatingv by means of an analytical-numerical method developed earlier by Zhang and Neuman. We present results in two-dimensional flow fields of unconditional (prior) mean uniformv . These show that the relative effect of local dispersion on first and second moments of concentration dies out locally as the corresponding dispersion tensor tends to zero. The effect also diminishes with time and source size. Our results thus do not support claims in the literature that local dispersion must always be accounted for, no matter how small it is. First-order decay reduces dispersion. This effect increases with time. However, these concentration moments c and cc of total concentrationc, which are associated with the scale below, cannot be used to estimate the field-scale concentrationc directly. To do so, a spatial average over the field measurement scale is needed. Nevertheless, our numerical results show that differences between the ensemble moments ofc and those ofc are negligible, especially for nonpoint sources, because the ensemble moments ofc are already smooth enough.     

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²     

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

The rheology of polymeric liquid crystals

The molecular theory of Doi has been used as a framework to characterize the rheological behavior of polymeric liquid crystals at the low deformation rates for which it was derived, and an appropriate extension for high deformation rates is presented. The essential physics behind the Doi formulation has, however, been retained in its entirety. The resulting four-parameter equation enables prediction of the shearing behavior at low and high deformation rates, of the stress in extensional flows, of the isotropic-anisotropic phase transition and of the molecular orientation. Extensional data over nearly three decades of elongation rate (10^–2–10¹) and shearing data over six decades of shear rate (10^–2–10⁴) have been correlated using this analysis. Experimental data are presented for both homogeneous and inhomogeneous shearing stress fields. For the latter, a 20-fold range of capillary tube diameters has been employed and no effects of system geometry or the inhomogeneity of the flow-field are observed. Such an independence of the rheological properties from these effects does not occur for low molecular weight liquid crystals and this is, perhaps, the first time this has been reported for polymeric lyotropic liquid crystals; the physical basis for this major difference is discussed briefly. A Semi-empirical constant in eq. (18), N/m² - c rod concentration, rods/m³ - c ^* critical rod concentration at which the isotropic phase becomes unstable, rods/m³ - C interaction potential in the Doi theory defined in eq. (3) - d rod diameter, m - D semi-empirical constant in eq. (19), s^–1 - D _r lumped rotational diffusivity defined in eq. (4), s^–1 - rotational diffusivity of rods in a concentrated (liquid crystalline) system, s^–1 - D _ro rotational diffusivity of a dilute solution of rods, s^–1 - f distribution function defining rod orientation - F tensorial term in the Doi theory defined in eq. (7) (or eq. (19)), s^–1 - G tensorial term in the Doi theory defined in eq. (8) - K _B Boltzmann constant, 1.38 × 10^–23 J/K-molecule - L rod length, m - S scalar order parameter - S tensor order parameter defined in eq. (5) - t time, s - T absolute temperature, K - u unit vector describing the orientation of an individual rod - rate of change ofu due to macroscopic flow, s^–1 - v fluid velocity vector, m/s - v velocity gradient tensor defined in eq. (9), s^–1 - V mean field (aligning) potential defined in eq. (2) - x coordinate direction, m - Kronecker delta (= 0 if = 1 if = ) - _r ratio of viscosity of suspension to that of the solvent at the same shear stress - _s solvent viscosity, Pa · s - ^* viscosity at the critical concentrationc ^*, Pa · s - v ₁, v₂ numerical factors in eqs. (3) and (4), respectively - deviatoric stress tensor, N/m² - volume fraction of rods - ⁰ constant in eq. (16) - ^* volume fraction of rods at the critical concentrationc ^* - average over the distribution functionf(u, t) (= d ²u f(u, t)) - gradient operator - d ²u integral over the surface of the sphere (|u| = 1)     

On the cross-effects of coupled macroscopic transport equations in porous media

In this paper, the derivation of macroscopic transport equations for this cases of simultaneous heat and water, chemical and water or electrical and water fluxes in porous media is presented. Based on themicro-macro passage using the method of homogenization of periodic structures, it is shown that the resulting macroscopic equations reveal zero-valued cross-coupling effects for the case of heat and water transport as well as chemical and water transport. In the case of electrical and water transport, a nonsymmetrical coupling was found.Notations b mobility - c concentration of a chemical - D rate of deformation tensor - D molecular diffusion coefficient - D _ij ^eff macroscopic (or effective) diffusion tensor - electric field - E ₀ initial electric field - k _ij molecular tensor - j, j ^*, current densities - K _ij macroscopic permeability tensor - l characteristic length of the ERV or the periodic cell - L characteristic macroscopic length - L _ijkl coupled flows coefficients - n _i unit outward vector normal to - p pressure - q _t ,q _t ⁺ , heat fluxes - q _c ,q _c ⁺ , chemical fluxes - s specific entropy or the entropy density - S entropy per unit volume - t time variable - t _ij local tensor - T absolute temperature - v _i velocity - V ₀ initial electric potential - V electric potential - x macroscopic (or slow) space variable - y microscopic (or fast) space variable - _i local vectorial field - _i local vectorial field - electric charge density on the solid surface - , bulk and shear viscosities of the fluid - _ij local tensor - _ij local tensor - _i local vector - _ij molecular conductivity tensor - _ij ^eff effective conductivity tensor - homogenization parameter - fluid density - ₀ ion-conductivity of fluid - _ij dielectric tensor - _i ¹ , _i ² , _i ³ local vectors - ⁴ local scalar - _S solid volume in the periodic cell - _L volume of pores in the periodic cell - boundary between _S and _L - ^s rate of entropy production per unit volume - total volume of the periodic cell - _l volume of pores in the cell On leave from the Politechnika Gdanska; ul. Majakowskiego 11/12, 80-952, Gdask, Poland.     

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     

Optical co-ordinates in general relativity 2 — the problem of distance

Summary Let denote the congruence of null geodesics associated with a given optical observer inV ₄. We prove that determines a unique collection of vector fieldsM₍₎ ( =1, 2, 3) and ₍₀₎ overV ₄, satisfying a weak version of Killing's conditions.This allows a natural interpretation of these fields as the infinitesimal generators of spatial rotations and temporal translation relative to the given observer. We prove also that the definition of the fieldsM₍₎ and ₍₀₎ is mathematically equivalent to the choice of a distinguished affine parameter f along the curves of, playing the role of a retarded distance from the observer.The relation between f and other possible definitions of distance is discussed.

---

Sommario Sia la congruenza di geodetiche nulle associata ad un osservatore ottico assegnato nello spazio-tempoV ₄. Dimostriamo che determina un'unica collezione di campi vettorialiM₍₎ ( =1, 2, 3) e ₍₀₎ inV ₄ che soddisfano una versione in forma debole delle equazioni di Killing. Ciò suggerisce una naturale interpretazione di questi campi come generatori infinitesimi di rotazioni spaziali e traslazioni temporali relative all'osservatore assegnato. Dimostriamo anche che la definizione dei campiM₍₎, ₍₀₎ è matematicamente equivalente alla scelta di un parametro affine privilegiato f lungo le curve di, che gioca il ruolo di distanza ritardata dall'osservatore. Successivamente si esaminano i legami tra f ed altre possibili definizioni di distanza in grande.

---

---

Work performed in the sphere of activity of: Gruppo Nazionale per la Fisica Matematica del CNR.     

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      

LDA measurements in low Reynolds number turbulent boundary layer

LDA measurements of the mean velocity in a low Reynolds number turbulent boundary layer allow a direct estimate of the friction velocity U from the value of /y at the wall. The trend of the Reynolds number dependence of / is similar to the direct numerical simulations of Spalart (1988).     

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     

Transcendentally small transversality in the rapidly forced pendulum

The rapidly forced pendulum equation with forcing sin((t/), where =<₀^p,p = 5, for ₀, sufficiently small, is considered. We prove that stable and unstable manifolds split and that the splitting distanced(t) in the ( ,t) plane satisfiesd(t) = sin(t/) sech(/2) +O( ₀ exp(–/2)) (2.3a) and the angle of transversal intersection,, in thet = 0 section satisfies 2 tan/2 = 2S _s = (/2) sech(/2) +O(( ₀ /) exp(–/2)) (2.3b) It follows that the Melnikov term correctly predicts the exponentially small splitting and angle of transversality. Our method improves a previous result of Holmes, Marsden, and Scheuerle. Our proof is elementary and self-contained, includes a stable manifold theorem, and emphasizes the phase space geometry.     

The stokes problems for a conducting fluid with a suspension of particles

The Stokes problems of an incompressible, viscous, conducting fluid with embedded small spherical particles over an infinite plate, set into motion in its plane by impulse and by oscillation, in the presence of a transverse magnetic field, are studied. The velocities of the fluid and of the particles and the wall shear stress are obtained. The stress is found to increase due to the particles and the magnetic field, with the effect of the particles diminishing as the field strength is increased.Nomenclature H ₀ strength of the imposed magnetic field - k density ratio of particles to fluid (per unit volume of flow field) - m _e ² H ₀ ² / - t time - y co-ordinate normal to the plate - u fluid velocity - v particle velocity - _e magnetic permeability of the fluid - kinematic viscosity of the fluid - electric conductivity of the fluid - fluid density - particle relaxation time - frequency of oscillation of the plate     

The Gibbs-Thompson relation within the gradient theory of phase transitions

This paper discusses the asymptotic behavior as 0+ of the chemical potentials associated with solutions of variational problems within the Van der Waals-Cahn-Hilliard theory of phase transitions in a fluid with free energy, per unit volume, given by ²¦¦²+ W(), where is the density. The main result is that is asymptotically equal to E/d+o(), with E the interfacial energy, per unit surface area, of the interface between phases, the (constant) sum of principal curvatures of the interface, and d the density jump across the interface. This result is in agreement with a formula conjectured by M. Gurtin and corresponds to the Gibbs-Thompson relation for surface tension, proved by G. Caginalp within the context of the phase field model of free boundaries arising from phase transitions.     

On the prediction of riblet performance with engineering turbulence models

The paper reports the outcome of a numerical study of fully developed flow through a plane channel composed of ribleted surfaces adopting a two-equation turbulence model to describe turbulent mixing. Three families of riblets have been examined: idealized blade-type, V-groove and a novel U-form that, according to computations, achieves a superior performance to that of the commercial V-groove configuration. The maximum drag reduction attained for any particular geometry is broadly in accord with experiment though this optimum occurs for considerably larger riblet heights than measurements indicate. Further explorations bring out a substantial sensitivity in the level of drag reduction to the channel Reynolds number below values of 15 000 as well as to the thickness of the blade riblet. The latter is in accord with the trends of very recent, independent experimental studies.Possible shortcomings in the model of turbulence are discussed particularly with reference to the absence of any turbulence-driven secondary motions when an isotropic turbulent viscosity is adopted. For illustration, results are obtained for the case where a stress transport turbulence model is adopted above the riblet crests, an elaboration that leads to the formation of a plausible secondary motion sweeping high momentum fluid towards the wall close to the riblet and thereby raising momentum transport.Nomenclature c _f Skin friction coefficient - c _f Skin friction coefficient in smooth channel at the same Reynolds number - k Turbulent kinetic energy - K ⁺ k/ _w - h Riblet height - S Riblet width - H Half height of channel - Re Reynolds number = volume flow/unit width/ - Modified turbulent Reynolds number - R _t turbulent Reynolds numberk ²/ - P _k Shear production rate ofk, _t (U _i /x _j + U _j /x _i ) U _i /x _j - dP/dz Streamwise static pressure gradient - U _i Mean velocity vector (tensor notation) - U Friction velocity, _w/ where _w=–H dP/dz - W Mean velocity - W _b Bulk mean velocity through channel - y ⁺ yU /v. Unless otherwise stated, origin is at wall on trough plane of symmetry - Kinematic viscosity - _t Turbulent kinematic viscosity - Turbulence energy dissipation rate - Modified dissipation rate – 2(k ^1/2/x _j )² - Density - _k , Effective turbulent Prandtl numbers for diffusion ofk and      

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.     

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     

The very slow flow of a Powell-Eyring fluid around a sphere

Summary The very slow flow of a Powell-Eyring type non-Newtonian fluid around a sphere is investigated by a variational technique. The result, a correction factor that is applied to the Stokes' equation, is given as a plot and as an equation which is empirically fit to the plot. Also, a comparison of the very slow flows of a simplified viscoelastic Oldroyd fluid and the Powell-Eyring fluid is made which indicates that in a certain restricted region of the very slow flows, both models give essentially the same results. The Oldroyd and Powell-Eyring model parameters are interrelated by forcing both models to fit the same tube flow viscosity data.Nomenclature B dimensionless quantity, v /R - C dimensionless second invariant - c ₁ constant determined by variational method - D dimensionless variational integral - D _2j , D _j+k position-independent variables used in specification of trial functions - E _2j , E _j+k position-independent variables used in specification of trial functions - f friction factor - f _corr friction factor correction - F _drag drag force on sphere - g, g ₀, g ₁ general trial function; first and second terms in the general trial function - G, H terms in the expression for C - j index - J variational integral - k index - K term in the expression for C - p, q integers - r integer, radial coordinate - R radius of sphere - Re Reynolds number - Re ₀ Reynolds number at point of zero shear rate - Re Reynolds number at infinite distance from sphere - Re _NN Reynolds number based on variable part of viscosity - u, v dimensionless position coordinates - V volume considered - v _i ith velocity component - v _r , v , v _z velocity components in the r, , and z-directions - v approach velocity of the fluid - x/ parameter in Powell-Eyring model - x _i i-position coordinate - parameter in Powell-Eyring model - rate of deformation - , _c , _N , ₀ coefficient of viscosity; cross viscosity; parameter in Powell-Eyring model; viscosity in limit of zero shear rate - spherical coordinate - , _ij rate of deformation tensor; ij-component of rate of deformation tensor - ₁, ₂ parameters in Oldroyd model - Newtonian viscosity - ₁, ₂ parameters in Oldroyd model - dimensionless radial coordinate, r/R - second invariant - fluid density - spherical coordinate - stream function     

Correlations between porous medium flow data and turbulent tube flow results for aqueous hydroxypropylguar solutions (Part 2)

Turbulent tube flow and the flow through a porous medium of aqueous hydroxypropylguar (HPG) solutions in concentrations from 100 wppm to 5000 wppm is investigated. Taking the rheological flow curves into account reveals that the effectiveness in turbulent tube flow and the efficiency for the flow through a porous medium both start at the same onset wall shear stress of 1.3 Pa. The similarity of the curves = ( _w ) and = ( _w ), respectively, leads to a simple linear relation / =k, where the constantk or proportionality depends uponc. This offers the possibility to deduce (for turbulent tube flow) from (for flow through a porous medium). In conjunction with rheological data, will reveal whether, and if yes to what extent, drag reduction will take place (even at high concentrations).The relation of our treatment to the model-based Deborah number concept is shown and a scale-up formula for the onset in turbulent tube flow is deduced as well.     

Modellvorstellung zur turbulenten Filmkondensation strömenden Dampfes

Zusammenfassung Der Wärmeübergang bei turbulenter Film kondensation strömenden Dampfes an einer waagerechten ebenen Platte wurde mit Hilfe der Analogie zwischen Impuls-und Wärmeaustausch untersucht. Zur Beschreibung des Impulsaustausches im Film wurde ein Vierbereichmodell vorgestellt. Nach diesem Modell wird die wellige Phasengrenze als starre rauhe Wand angesehen. Die Abhängigkeit einer Schubspannungs-Nusseltzahl von der Film-Reynoldszahl und Prandtlzahl wurde berechnet und dargestellt.

---

A model for turbulent film condensation of flowing vapour

The heat transfer in turbulent film condensation of flowing vapour on a horizontal flat plate was investigated by means of the analogy between momentum and heat transfer. To describe the momentum transfer in the film a four-region model was presented. With this model the wavy interfacial surface is treated as a stiff rough wall. A shear Nusselt number has been calculated and represented as a function of film Reynolds number and Prandtl number.

---

Formelzeichen a Temperaturleitkoeffizient - k Mischungswegkonstante - k _s äquivalente Sandkornrauhigkeit - Nu _x lokale Schubspannungs-Nusseltzahl,Nu _x=x_xv/uw - Pr Prandtlzahl,Pr=v/a - Pr _t turbulente Prandtlzahl,Pr _t =_m/_q - q Wärmestromdichte _q - R Wärmeübergangswiderstand - R_f Wärmeübergangswiderstand des Films - Re _F Reynoldszahl der Filmströmung - T Temperatur - U, V Geschwindigkeitskomponenten des Dampfes in waagerechter und senkrechter Richtung - u, Geschwindigkeitskomponenten des Kondensats in waagerechter und senkrechter Richtung - V Querschwankungsgeschwindigkeit des Kondensats und des Dampfes - u _/gtD Schubspannungsgeschwindigkeit an der Phasengrenze für die Dampfgrenzschicht, uD =(/)^1/2 - u _F Schubspannungsgeschwindigkeit an der Phasengrenze für den Kondensatfilm,u _F =(/)^1/2 - u _w Schubspannungsgeschwindigkeit an der Wand der Kühlplatte,u _w =(w/)^1/2 - y Wandabstand - x Wärmeübergangskoeffizient - gemittelte Kondensatfilmdicke - _s Dicke der zähen Schicht der Filmströmung an der welligen Phasengrenze - ₄ Dicke der zähen Schicht der Filmströmung an der gemittelten glatten Phasengrenze - Wärmeleitzahl - dynamische Viskosität - v kinematische Viskosität - Dichte - Oberflächenspannung - _w Wandschubspannung - Schubspannung an der Phasengrenzfläche - _m turbulente Impulsaustauschgröße - _q turbulente Wärmeaustauschgröße Indizes d Wert des Dampfes - w Wert an der Wand - x lokaler Wert inx - Wert an der Phasengrenze Stoffgrößen ohne Index gelten für das Kondensat