Heat transfer between gas-liquid spray stream flowing perpendicularly to the row of the cylinders

Experimental investigation and analysis of heat transfer process between a gas-liquid spray flow and the row of smooth cylinders placed in the surface perpendicular to the flow has been performed. Among others, there was taken into account in the analysis the phenomenon of droplets bouncing and omitting the cylinder as well as the phenomenon of the evaporation process from the liquid film surface.In the experiments test cylinders were used, which were placed between two other cylinders standing in the row.From the experiments and the analysis the conclusion can be drawn that the heat transfer coefficients values for a row of the cylinders are higher than for a single cylinder placed in the gasliquid spray flow.

---

Wärmeübergang an eine senkrecht anf eine Zylinderreihe auftreffende Gas-Flüssigkeits-Sprüh-Strömung

Zusammenfassung Es wurden Messungen und theoretische Analysen des Wärmeübergangs zwischen einer Gas-FlüssigkeitsSprüh-Strömung und den glatten Oberflächen einer Zylinderreihe durchgeführt, die senkrecht zum Sprühstrahl angeordnet waren. Dabei wurde in der Analyse unter anderem das Phänomen betrachtet, daß die Tropfen die Zylinderwand treffen und verfehlen können und daß sich ein Verdampfungsprozeß aus dem flüssigen Film an der Zylinderoberfläche einstellt.Gemessen wurde an einem zwischen zwei Randzylindern befindlichen Zylinder.Die Experimente und die Analyse gestatten die Schlußfolgerung, daß der Wärmeübergangskoeffizient für eine Zylinderreihe höher ist als für einen einzelnen Zylinder in der Sprühströmung.

---

Nomenclature a distance between axes of cylinders, m - c _l specific heat capacity of liquid, J/kg K - c _g specific heat capacity of gas, J/kg K - D cylinder diameter, m - g _l mass velocity of liquid, kg/m²s - ¯k average volume ratio of liquid entering film to amount of liquid directed at the cylinder in gas-liquid spray flow, dimensionless - k() local volume ratio of liquid entering film to amount of liquid directed at the cylinder in gas-liquid spray flow, dimensionless - L specific latent heat of vaporisation, J/kg - m mass fraction of water in gas-liquid spray flow, dimensionless - M constant in Eq. (9) - p pressure, Pa - p _g statical pressure of gas, Pa - p _w pressure of gas on the cylinder surface, Pa - p external pressure on the liquid film surface, Pa - r cylindrical coordinate, m - R radius of cylinder, m - T temperature, K, °C - T _l, t_l liquid temperature in the gas-liquid spray, K, °C - T _w,t_w temperature of cylinder surface, K, °C - T temperature of gas-liquid film interface, K - U liquid film velocity, m/s - w gas velocity on cylinder surface, m/s - w _g gas velocity in free stream, m/s - W _l liquid vapour mass ratio in free stream, dimensionless - W liquid vapour mass ratio at the edge of a liquid film, dimensionless - x coordinate, m - y coordinate, m - z complex variable, dimensionless - average heat transfer coefficient, W/m²K - local heat transfer coefficient, W/m² K - average heat transfer coefficient between cylinder surface and gas, W/m² K - _g, local heat transfer coefficient between cylinder surface and gas, W/m² K - mass transfer coefficient, kg/m²s - liquid film thickness, m - _lg dynamic diffusion coefficient of liquid vapour in gas, kg/m s - pressure distribution function on a cylinder surface - function defined by Eq. (3) - _l liquid dynamic viscosity, kg/m s - _g gas dynamic viscosity, kg/m s - cylindrical coordinate, rad, deg - _l thermal conductivity of liquid, W/m K - _g thermal conductivity of gas, W/m K - mass transfer driving force, dimensionless - _l density of liquid, kg/m³ - _g density of gas, kg/m³ - _w shear stress on the cylinder surface, N/m² - _w shear stress exerted by gas at the liquid film surface, N/m² - air relative humidity, dimensionless - T -T _w - _w =T _wT_l Dimensionless parameters I= enhancement factor of heat transfer - m ^*=M _l/M_g molar mass of liquid to the molar mass of gas ratio - Nu _g= D/ _g gas Nusselt number - Pr _g=c _{g g}/_g gas Prandtl number - Pr _l=c_lll liquid Prandtl number - ¯r=(r–R)/ dimensionless coordinate - Re _g=w_gD _g/_g gas Reynolds number - Re _g,max=w_g,max D _g/_g gas Reynolds number calculated for the maximal gas velocity between the cylinders - Sc=m ^* _g/_l–g Schmidt number =/R dimensionless film thickness     

The dynamic and steady shear properties of synovial fluid and of the components making up synovial fluid

Steady-shear and dynamic properties of a pooled sample of cattle synovial fluid have been measured using techniques developed for low viscosity fluids. The rheological properties of synovial fluid were found to exhibit typical viscoelastic behaviour and can be described by the Carreau type A rheological model. Typical model parameters for the fluid are given; these may be useful for the analysis of the complex flow problems of joint lubrication.The two major constituents, hyaluronic acid and proteins, have been successfully separated from the pooled sample of synovial fluid. The rheological properties of the hyaluronic acid and the recombined hyaluronic acid-protein solutions of both equal and half the concentration of the constituents found in the original synovial fluid have been measured. These properties, when compared to those of the original synovial fluid, show an undeniable contribution of proteins to the flow behaviour of synovial fluid in joints. The effect of protein was found to be more prominent in hyaluronic acid of half the normal concentration found in synovial fluid, thus providing a possible explanation for the differences in flow behaviour observed between synovial fluid from certain diseased joints compared to normal joint fluid.Nomenclature A Ratio of angular amplitude of torsion head to oscillation input signal - G Storage modulus - G Loss modulus - I Moment of inertia of upper platen — torsion head assembly - K Restoring constant of torsion bar - N ₁ First normal-stress difference - R Platen radius - S ⁽ⁱ⁾ Geometric factor in the dynamic property analysis - t ₁ Characteristic time parameter of the Carreau model - X, Y Carreau model parameters - Z () Reimann Zeta function of - Carreau model parameter - Shear rate - Apparent steady-shear viscosity - ^* Complex dynamic viscosity - Dynamic viscosity - Imaginary part of the complex dynamic viscosity - ₀ Zero-shear viscosity - ₀ Cone angle - Carreau model characteristic time - Density of fluid - Shear stress - Phase difference between torsion head and oscillation input signals - ₀ Zero-shear rate first normal-stress coefficient - Oscillatory frequency     

Darcy-Forchheimer natural,forced and mixed convection heat transfer in non-Newtonian power-law fluid-saturated porous media

The governing equation for Darcy-Forchheimer flow of non-Newtonian inelastic power-law fluid through porous media has been derived from first principles. Using this equation, the problem of Darcy-Forchheimer natural, forced, and mixed convection within the porous media saturated with a power-law fluid has been solved using the approximate integral method. It is observed that a similarity solution exists specifically for only the case of an isothermal vertical flat plate embedded in the porous media. The results based on the approximate method, when compared with existing exact solutions show an agreement of within a maximum error bound of 2.5%.Nomenclature A cross-sectional area - b _i coefficient in the chosen temperature profile - B ₁ coefficient in the profile for the dimensionless boundary layer thickness - C coefficient in the modified Forchheimer term for power-law fluids - C ₁ coefficient in the Oseen approximation which depends essentially on pore geometry - C _i coefficient depending essentially on pore geometry - C _D drag coefficient - C _t coefficient in the expression forK ^* - d particle diameter (for irregular shaped particles, it is characteristic length for average-size particle) - f _p resistance or drag on a single particle - F _R total resistance to flow offered byN particles in the porous media - g acceleration due to gravity - g _x component of the acceleration due to gravity in thex-direction - Grashof number based on permeability for power-law fluids - K intrinsic permeability of the porous media - K ^* modified permeability of the porous media for flow of power-law fluids - l _c characteristic length - m exponent in the gravity field - n power-law index of the inelastic non-Newtonian fluid - N total number of particles - Nux,_D,F local Nusselt number for Darcy forced convection flow - Nux,_D-F,F local Nusselt number for Darcy-Forchheimer forced convection flow - Nux,_D,M local Nusselt number for Darcy mixed convection flow - Nux,_D-F,M local Nusselt number for Darcy-Forchheimer mixed convection flow - Nux,_D,N local Nusselt number for Darcy natural convection flow - Nux,_D-F,N local Nusselt number for Darcy-Forchheimer natural convection flow - pressure - p exponent in the wall temperature variation - _Pe c characteristic Péclet number - _Pe x local Péclet number for forced convection flow - _Pe x modified local Péclet number for mixed convection flow - _Ra c characteristic Rayleigh number - _Ra x local Rayleigh number for Darcy natural convection flow - _Ra x local Rayleigh number for Darcy-Forchheimer natural convection flow - Re convectional Reynolds number for power-law fluids - Reynolds number based on permeability for power-law fluids - T temperature - T _e ambient constant temperature - T w,_ref constant reference wall surface temperature - T _w(X) variable wall surface temperature - T _w temperature difference equal toT w,_ref–T _e - T ₁ term in the Darcy-Forchheimer natural convection regime for Newtonian fluids - T ₂ term in the Darcy-Forchheimer natural convection regime for non-Newtonian fluids (first approximation) - T _N term in the Darcy/Forchheimer natural convection regime for non-Newtonian fluids (second approximation) - u Darcian or superficial velocity - u ₁ dimensionless velocity profile - u _e external forced convection flow velocity - u _s seepage velocity (local average velocity of flow around the particle) - u _w wall slip velocity - U c _M characteristic velocity for mixed convection - U c _N characteristic velocity for natural convection - x, y boundary-layer coordinates - x ₁,y ₁ dimensionless boundary layer coordinates - X coefficient which is a function of flow behaviour indexn for power-law fluids - effective thermal diffusivity of the porous medium - shape factor which takes a value of/4 for spheres - shape factor which takes a value of/6 for spheres - ₀ expansion coefficient of the fluid - _T boundary-layer thickness - T ₁ dimensionless boundary layer thickness - porosity of the medium - similarity variable - dimensionless temperature difference - coefficient which is a function of the geometry of the porous media (it takes a value of 3 for a single sphere in an infinite fluid) - ₀ viscosity of Newtonian fluid - ^* fluid consistency of the inelastic non-Newtonian power-law fluid - constant equal toX(2 ^2–n )/ - density of the fluid - dimensionless wall temperature difference     

Entrance flow of a yield-power law fluid

In the present work we have obtained the numerical solution of the momentum equation for a Yield-Pseudoplastic power-law fluid flowing in the entrance region of a tube. The accuracy of the numerical results is checked by comparing the asymptotic values of friction coefficients and velocity profiles with the corresponding results from the analytical solutions for the fully-developed region. The results of the entrance flow solution for the power-law exponent equal to unity (Bingham fluid) are also in agreement with the numerical solution for a Bingham fluid. Detailed results are presented for wide ranges of yield numbers and power law exponents.

Nomenclature

Nomenclature a constant - D diameter - F dimensionless pressure gradient in (4.3) - f _x friction factor in (5.1) - f _app total friction factor in (5.2) - K entrance pressure drop coefficient - n power law exponent - p pressure - r radial co-ordinate - R radius of a tube - Re Reynolds number (5.3) - s rate of shear, u/r - u axial velocity - average velocity - v velocity in radius direction - x axial co-ordinate - y normal co-ordinate - Y yield number in (4.4) - z dimensionless axial distance =(x/D)/Re - z ₁ 1/z Greek Symbols plug flow radius in (4.6) - _eff effective viscosity - density - shear stress - _y yield stress - dimensionless stream function     

The swirling radial jet

The similarity solution of the radial turbulent jet with weak swirl is discussed and a new solution of the radial turbulent jet with swirl is proposed without restrictions assumed in the weak swirl solution.Nomenclature e swirl parameter - k experimental constant - l non-negative constant - M, M , N, P integral invariants - q velocity component in -direction - q _max maximum velocity component in -direction - u radial velocity component - u _max maximum radial velocity component - v axial velocity component - w peripheral velocity component - w _max maximum peripheral velocity component - x radial coordinate - y transverse coordinate - angle introduced in (28) - characteristic width of a jet - (x, y) similarity variable (scaled x and y coordinate) - molecular kinematic viscosity - _T eddy kinematic viscosity - tangential coordinate - fluid density - turbulent shear stress in -direction - _xy , _y components of turbulent shear stress tensor - (x, y) stream function     

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     

Injection moulding of thermoplastics: Freezing of variable-viscosity fluids.

The injection moulding of thermoplastics involves, during mould filling, flows of hot polymer melts into mould networks, the walls of which are so cold that frozen layers form on them. An analytical study of such flows is presented here for the case when the Graetz and Nahme numbers are large and the Pearson number is small. Thus the flows are developing and temperature differences due to heat generation by viscous dissipation are sufficiently large to cause significant variations in viscosity (but the difference between the entry temperature of the polymer to a specific part of the mould network and the melting temperature of the polymer is not). Br Brinkman number - Gz Graetz number - h half-height of channel or disc - h ^* half-height of polymer melt region in channel or disc - L length of channel or pipe - m viscosity shear-rate exponent - Na Nahme number - p pressure - P pressure drop - Pe Péclet number - Pn Pearson number - Q volumetric flowrate - r radial coordinate in pipe or disc - R radius of pipe - Re Reynolds number - R _i inner radius of disc - R _o outer radius of disc - R ^* radius of polymer melt region in pipe - T temperature - T _ad adiabatic temperature rise - T _e entry polymer melt temperature - T _m melting temperature of polymer - T _max maximum temperature - T ₀ reference temperature - T _w wall temperature - flow-average temperature rise - u _r radial velocity in pipe or disc - u _x axial velocity in channel - u _y transverse velocity in channel or disc - u _z axial velocity in pipe - w width of channel - x axial coordinate in channel or modified radial coordinate in disc - y transverse coordinate in channel or disc - z axial coordinate in pipe - thermal conductivity of molten polymer - thermal conductivity of frozen polymer - scaled dimensionless axial coordinate in channel or pipe or radial coordinate in disc - ₀ undetermined integration constant - heat capacity of molten polymer - viscosity temperature exponent - dimensionless transverse coordinate in channel or disc - ^* dimensionless half-height of polymer melt region in channel or disc - H ^* scaled dimensionless half-height of polymer melt region in channel or disc or radius of polymer melt region in pipe - dimensionless temperature - ^* dimensionless wall temperature - scaled dimensionless temperature - numerical constant - µ viscosity of molten polymer - µ ₀ consistency of molten polymer - dimensionless pressure gradient - scaled dimensionless pressure gradient - density of molten polymer - dimensionless radial coordinate in pipe or disc - _i dimensionless inner radius of disc - ^* dimensionless radius of polymer melt region in pipe - dimensionless streamfunction - scaled dimensionless streamfunction - dummy variable - streamfunction - similarity variable - similarity variable     

Velocity and concentration measurements in a model diesel engine

Laser Doppler anemometry and Rayleigh scattering have been used to quantify the velocity and concentration fields after the start of injection in a model diesel engine motored at 200 rpm in the absence of compression. Fuel injection was simulated by a transient jet of vapour Freon-12 initiated at 40 degrees before top-dead-centre through a nozzle incorporated into the centre of a permanently open intake valve. Swirl was induced by means of 60 degree vanes located in the inlet, port. The piston configurations comprised a flat and a re-entrant piston-bowl.The results indicate that for the two nozzle geometries investigated the mass flux decays faster than momentum with nearly constant decay rates along the centreline. The nozzle with the larger exit diameter and wider jet angle gave rise to slower decay of both mass and momentum with associated lower velocity and concentration fluctuations.List of symbols D ₀ nozzle diameter - r radial coordinate - mean axial velocity - mean axial velocity at the centreline - ₀ mean axial velocity at the nozzle exit - rms of axial velocity fluctuations - mean concentration (mole fraction) - mean concentration at the nozzle exit - rms of concentration fluctuations - x axial coordinate A version of this paper was presented at the ASME Winter Annual Meeting of 1984 and printed in AMD, Vol. 66     

Squeezing flows of Newtonian liquid films an analysis including fluid inertia

A theoretical study is made of the flow behavior of thin Newtonian liquid films being squeezed between two flat plates. Solutions to the problem are obtained by using a numerical method, which is found to be stable for all Reynolds numbers, aspect ratios, and grid sizes tested. Particular emphasis is placed on including in the analysis the inertial terms in the Navier-Stokes equations.Comparison of results from the numerical calculation with those from Ishizawa's perturbation solution is made. For the conditions considered here, it is found that the perturbation series is divergent, and that in general one must use a numerical technique to solve this problem.Nomenclature a half of the distance, or gap, between the two plates - a ₀ the value of a at time t=0 - adot da/dt - ä d² a/dt ² - d³ a/dt ³ - a ⁱ components of a contravariant acceleration vector - f unknown function of z ₀ and t defined in (6) - f ⁱ function defined in (9) f ¹=r ₀ g(z ₀, t) f ²= ₀ f ³=f(z ₀, t) - F force applied to the plates - g unknown function of z ₀ and t defined in (6) - g g/z ₀ - h grid dimension in the z ₀ direction (see Fig. 5) - Christoffel symbol - i, j, k, l indices - k grid dimension in the t direction (see Fig. 5) - r radial coordinate direction defined in Fig. 1 - r ₀ radial convected coordinate - R radius of the circular plates - t time - v _r fluid velocity in the r direction - v _z fluid velocity in the z direction - v fluid velocity in the direction - x ⁱ cylindrical coordinate x ¹=r x²= x³=z - z vertical coordinate direction defined in Fig. 1 - z ₀ vertical convected coordinate - tangential coordinate direction - ₀ tangential convected coordinate - viscosity - kinematic viscosity, / - ⁱ convected coordinate ¹=r₀ ²=₀ ³=z₀ - density     

Simultaneous measurements of velocity and temperature fluctuations in thermal boundary layer in a drag-reducing surfactant solution flow

The mechanism of turbulent heat transfer in the thermal boundary layer developing in the channel flow of a drag-reducing surfactant solution was studied experimentally. A two-component laser Doppler velocimetry and a fine-wire thermocouple probe were used to measure the velocity and temperature fluctuations simultaneously. Two layers of thermal field were found: a high heat resistance layer with a high temperature gradient, and a layer with a small or even zero temperature gradient. The peak value of was larger for the flow with the drag-reducing additives than for the Newtonian flow, and the peak location was away from the wall. The profile of was depressed in a similar manner to the depression of the profile of in the flow of the surfactant solution, i.e., decorrelation between v and compared with decorrelation between u and v. The depression of the Reynolds shear stress resulted in drag reduction; similarly, it was conjectured that the heat transfer reduction is due to the decrease in the turbulent heat flux in the wall-normal direction for a flow with drag-reducing surfactant additives.List of symbols ensemble averaged value - (·)⁺ normalized by the inner wall variables - (·) root-mean-square value - C concentration of cetyltrimethyl ammonium chloride (CTAC) solution - c _p heat capacity - D hydraulic diameter - f friction factor - H channel height - h heat transfer coefficient - j _H Colburn factor - l length - Nu Nusselt number, h - Pr Prandtl number, c _p/ - q _w wall heated flux - Re Reynolds number, U _b/ - T temperature - T _b bulk temperature - T _i inlet temperature - T _w wall temperature - T friction temperature, q _w /c _p u - U local time-mean streamwise velocity - U ₁ velocity signals from BSA1 - U ₂ velocity signals from BSA2 - U _b bulk velocity - u streamwise velocity fluctuation - u1 velocity in abscissa direction in transformed coordinates - u friction velocity, - v wall-normal velocity fluctuation - v1 velocity in ordinate direction in transformed coordinates - var(·) variance - x streamwise direction - y wall-normal direction - z spanwise direction - _j junction diameter of fine-wire TC - _w wire diameter of fine-wire TC - angle of principal axis of joint probability function p(u,v) - _f heat conduction of fluid - _w heat conduction of wire of fine-wire TC - kinematic viscosity - local time-mean temperature difference, T _w –T - temperature fluctuation - standard deviation - density - _w wall shear stress     

Eddy viscosity in decaying swirl flow in a pipe

Prediction of heat transfer coefficient for swirling flows can be made provided the values of the eddy viscosity are available. In the present work the axial and tangential velocity fields are surveyed in a pipe for the determination of eddy viscosity. The data thus obtained were utilised to determine the influence of the axial Reynolds number and swirl number on the eddy viscosity. An empirical relationship is suggested to determine the eddy viscosity as a function of Reynolds number and swirl number.Nomenclature A _T angular momentum, equation (10) - a coefficient, equation (1) - b coefficient, equation (1) - D pipe diameter - f friction factor - F(y) initial condition function, equation (8) - J ₀ Bessel's function of the first kind of order zero - J ₁ Bessel's function of the first kind of order one - R pipe radius - Re Reynolds number, u _av D/ - r radial coordinate - S _n swirl number, equation (6) - (S _n )_in swirl number at the inlet of the test pipe - u axial velocity - u _av mean axial velocity in pipe - W non-dimensional local tangential velocity, w/u _av - w tangential velocity - X non-dimensional axial coordinate, x/D - x axial coordinate - y non-dimensional radial coordinate, r/R - z non-dimensional parameter, 4(1+/)/Re(x/D) - kinematic eddy viscosity - _n eigenvalues, equation (7) - kinematic viscosity - density     

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     

Temperature variations across the lubricant film in hydrodynamic lubrication

Summary Temperature variations across the lubricant film in hydrodynamic lubrication have been taken into account. The consequent variations of viscosity cannot be neglected for high Prandtl number lubricants.Nomenclature A constant defined by eq. (13a) - b h ₀/h _L - c heat capacity - h film thickness - h ₀ h at X=0 - h _L h at X=L - H ₀ local heat transfer coefficient defined by (23a) - K thermal conductivity - L length of bearing - m* defined by (9) - Nu Nusselt number defined by (23b) - P pressure - Pe Peclet number (Re) (Pr) - q ₀ slider surface heat flux - q ₀ ^* dimensionless heat flux defined by (20) - Q ₀ slider surface total heat transfer - Q ₀ ^* dimensionless total heat transfer defined by (21) - Re Reynolds number Vh ₀/ - Pr Prandtl number c/ - T temperature - T ₀ slider surface temperature - T _B film bulk temperature - u longitudinal velocity - v transverse velocity - V slider velocity - x longitudinal coordinate - y transverse coordinate - x/h ₀ - u/V - h/h ₀ - y/h - density - viscosity - v/V - - - _S dimensionless stationary surface temperature - dimensionless average stationary surface temperature - _B dimensionless film bulk temperature     

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     

Mean and turbulent velocity measurements of supersonic mixing layers

The behavior of supersonic mixing layers under three conditions has been examined by schlieren photography and laser Doppler velocimetry. In the schlieren photographs, some large-scale, repetitive patterns were observed within the mixing layer; however, these structures do not appear to dominate the mixing layer character under the present flow conditions. It was found that higher levels of secondary freestream turbulence did not increase the peak turbulence intensity observed within the mixing layer, but slightly increased the growth rate. Higher levels of freestream turbulence also reduced the axial distance required for development of the mean velocity. At higher convective Mach numbers, the mixing layer growth rate was found to be smaller than that of an incompressible mixing layer at the same velocity and freestream density ratio. The increase in convective Mach number also caused a decrease in the turbulence intensity ( _u/U).List of symbols a speed of sound - b total mixing layer thickness between U ₁ – 0.1 U and U ₂ + 0.1 U - f normalized third moment of u-velocity, f u³/(U)³ - g normalized triple product of u² , g u²/(U)³ - h normalized triple product of u ², h u ²/(U)³ - l _u axial distance for similarity in the mean velocity - l _u axial distance for similarity in the turbulence intensity - M Mach number - M _c convective Mach number (for ₁ = ₂), M _c (U ₁ – U ₂)/(a ₁ + a ₂) - P static pressure - r freestream velocity ratio, r U ₂/U ₁ - Re unit Reynolds number, Re U/ - s freestream density ratio, s ₂/₁ - T _t total temperature - u instantaneous streamwise velocity - u deviation of u-velocity, uu – U - U local mean streamwise velocity - U ₁ primary freestream velocity - U ₂ secondary freestream velocity - average of freestream velocities, (U ₁ + U ₂)/2 - U freestream velocity difference, U U ₁ – U ₂ - instantaneous transverse velocity - v deviation of -velocity, – V - V local mean transverse velocity - x streamwise coordinate - y transverse coordinate - y ₀ transverse location of the mixing layer centerline - ensemble average - ratio of specific heats - boundary layer thickness (y-location at 99.5% of free-stream velocity) - similarity coordinate, (y – y ₀)/b - compressible boundary layer momentum thickness - viscosity - density - standard deviation - dimensionless velocity, (U – U ₂)/U - 1 primary stream - 2 secondary stream A version of this paper was presented at the 11th Symposium on Turbulence, October 17–19, 1988, University of Missouri-Rolla     

On laminar flow through a uniformly porous pipe

Numerous investigations ([1] and [4–9]) have been made of laminar flow in a uniformly porous circular pipe with constant suction or injection applied at the wall. The object of this paper is to give a complete analysis of the numerical and theoretical solutions of this problem. It is shown that two solutions exist for all values of injection as well as the dual solutions for suction which had been noted by previous investigators. Analytical solutions are derived for large suction and injection; for large suction a viscous layer occurs at the wall while for large injection one solution has a viscous layer at the centre of the channel and the other has no viscous layer anywhere. Approximate analytic solutions are also given for small values of suction and injection.

Nomenclature

General r distance measured radially - z distance measured along axis of pipe - u velocity component in direction of z increasing - v velocity component in direction of r increasing - p pressure - density - coefficient of kinematic viscosity - a radius of pipe - V velocity of suction at the wall - r ²/a ² - R wall or suction Reynolds number, Va/ - f() similarity function defined in (6) - u ₀() eigensolution - U(0) a velocity at z=0 - K an arbitrary constant - B _K Bernoulli numbers Particular Section 5 perturbation parameter, –2/R - ² a constant, –K - x / - g(x) f()/ Section 6 perturbation parameter, –R/2 - ² a constant, –K - g() f() - g _c ()=g() near centre of pipe - * point where g()=0 Section 7 2/R - ² K - t (1–)/ - w(t, ) [1–f(t)]/ - ₀, ₁ constants - g() f()– ₀ - ₀/ - ₀ a constant - * point where f()=0     

The rheology of pigment dispersions

Summary A kinetic model is developed to relate the measured shear stress in a dispersion with the rate of deformation, and with the level of structure caused by the competing effects of flocculation and deflocculation.The model parameters are determined from experimental data obtained from equilibrium and transient oscillatory shear, using dispersions of a pigment in three different oil-based media. It is found that the model can successfully describe the flow behaviour of the dispersions under all three types of deformation, and account for different concentrations and temperatures.

---

Zusammenfassung Es wird ein kinetisches Modell entwickelt, das die in einer Dispersion gemessene Schubspannung mit der Deformationsgeschwindigkeit in Beziehung setzt unter Berücksichtigung der im Wettbewerb stehenden Flokkulations-und Deflokkulationseffekte.Die Modellparameter werden mit Hilfe experimenteller Daten bestimmt, die zu stationären und oszillatorischen Scherströmungen unter Einbeziehung wechselnder Beanspruchungsarten erhalten wurden. Dabei wurden Pigment-Dispersionen in drei verschiedenen Medien auf Öl-Basis verwendet. Man findet, daß das Modell das Fließverhalten der Dispersionen unter allen betrachteten Deformationstypen sowie bei den verschiedenen angewandten Konzentrationen und Temperaturen erfolgreich zu beschreiben vermag.

---

a, b experimental constants - c dispersion concentration by weight - d ₃ mean volume to surface diameter of a floc. - f (·,·) function defined by eq. [20] - g(·) a function of volume fraction - k Boltzmann's constant - n _s number of floccules containings flocs, per unit volume - ratio of the number of floccules containings flocs, per unit volume, to the total number of flocs per unit volume - t present time - B, D derived constants - E ₀,E ₁ activation energies for viscous flow of a dispersion at low and high shear rates - E _m activation energy for viscous flow of the medium - E _f ,E _d activation energies for flocculation and deflocculation of a dispersion - G shear rate - K ₁,K ₂ model constants - M parameter related to the total number of flow units - N number of flocs per unit volume - R, R ₀ the ratio of the rate function for deflocculation to the rate function for flocculation, and its value in the absence of shear - T absolute temperature - ₀, ₀ constants - (·,·),(·) rate functions for flocculation and deflocculation - , , v model parameters - , ₀, ₁ viscosity of a dispersion, and its value at low and high shear rates - _r the viscosity of a dispersion of floccules each containingr flocs - dynamic viscosity of a dispersion - µ, µ ₀ viscosity of a medium at temperatureT, and in the limit of high temperature - _p , _f , _F volume fraction of pigment, of flocs, of floccules - measured shear stress - non-dimensional time - a characteristic time for flocculation With 8 figures and 7 tables     

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     

Analysis of a method and a system to determine heat transfer and properties of fluids in the absence of temperature measurements and a modified Fourier Equation

The technique to determine by capacitance measurements heat transfer, thermal transport and dielectric properties of fluids introduced recently is now analyzed for a simple system of spherical geometry. The temperature distribution under programmed heat input to a fluid annulus between solid walls is computed by finite difference method for the determination of the capacitance time function of the arrangement. A system of heavy wall structure and heated long enough will produce a capacitance-time curve which is a function of thermal conductivity only. Thermal diffusivity is of influence in thin wall systems. The capacitance change of a heavy wall arrangement is related to the thermal conductivity of the test fluid by a modified Fourier equation. This equation describes the heat flow through the fluid layer but includes the thermal expansion of the solid walls. The change of geometry with T is therefore accounted for. For other multicomposite structures the Fourier equation must be further modified by including the thermal expansion of all materials of the structure and possibly also their compressibilities.

---

Zusammenfassung Die kürzlich eingeführte Methode der Bestimmung von Wärmeübergang, thermischen Transport und dielektrischen Größen mittels Kapazitäts-Zeit-Messung wird analysiert für ein einfaches kugeliges System. Die Temperaturverteilung in der Flüssigkeit im Kugelspalt zwischen zwei festen Körpern wird für konstante Wärmezufuhr von außen mittels der Differenzmethode bestimmt und daraus die Kapazitäts-Zeit-Funktion ermittelt. Es wird gezeigt, daß die Kapazitäts-Zeit-Kurve nur eine Funktion der Wärmeleitzahl ist für den Fall dickwandiger Anordnungen. Für dünnwandige Systeme wird sie auch abhängig von der Temperaturleitzahl. Es wird eine modifizierte Fourier-Gleichung eingeführt, die den Wärmetransport durch die Flüssigkeit beschreibt, dabei aber die Änderung der Geometrie der Schicht berücksichtigt, die sich wegen der thermischen Ausdehnung der festen Wände bei der Einstellung der Temperaturdifferenz ergibt. Für andere mehrschichtige Körper muß die Fourier-Gleichung weiterhin modifiziert werden durch Berücksichtigung der thermischen Ausdehnungskoeffizienten aller beteiligten Materialien und möglicherweise auch ihrer Kompressibilitäten.

---

Nomenclature A average cross-sectional area of fluid layer - A coefficient matrix - B matrix defined by Eq. (20) - B₀ geometric constant of fluid layer (A/L) at reference temperature - C capacitance of arrangement - C_i, C_r capacitance of layer of fluid i and reference fluid at temperature T - capacitances at reference temperature - C_H, c_l specific heats of outer and inner wall - FA...FE constants defined in Eqs. (13 ... 17) - L thickness of fluid layer - M_H, M_L mass of outer and inner wall - P power input to the system - R constant defined by Eq. (24) - T temperature - T_ref reference temperature - T (O, t), T (L, t) temperatures of outer and inner wall at time t - T _i ⁿ , T _i+0 ^n+m temperatures at location i and time n (m=number of t's; 0=number of x's) - T temperature difference across fluid layer - T apparent temperature difference - t_h, T_l temperature increases of outer and inner wall - T_max temperature change of system from one to another thermal equilibrium condition a thermal diffusivity - k, k_i, k_r thermal qonductivity of fluids and of fluid i and reference fluid - q heat flow through fluid layer - ^rh,^rl inner radius of outer wall and outer radius of inner wall - ^rOH,^rOL radii at reference temperature - t time - t time interval - x coordinate - ¯x vector of unknown T_i ⁿ⁺¹ - x length interval Greek symbols linear thermal expansion coefficient - _H, _L linear thermal expansion coefficient of materials of outer and inner wall - dielectric constant - _i, _ref dielectric constant of fluid i and reference fluid - ₀ permittivity of free space - multiplyer of conduction Eq. (7) in finite difference form - time needed to establish quasi-steady state conditions in the system heated by a constant power input In honor of Prof. Dr. E. Schmidt to his 80th Birthday.     

On the effect of magnetic field on viscous lifting and drainage of conducting fluid

This paper concerns with obtaining the solution of the problem of viscous lifting and drainage of a thin liquid film clinging to a vertical plane surface moving with a velocityf(t) in the presence of a transverse magnetic field. Specializing to the case when the surface moves with a constant acceleration, it has been found that the film thickness, for large magnetic fields, increases with the increase in magnetic field.Nomenclature a acceleration of the plate - A non-dimensional acceleration, =a/g - B magnetic induction vector - B ₀ applied magnetic field - f(t) any function oft - Laplace transform off(t) - g gravitational acceleration - h film thickness - H non-dimensional film thickness, =h(g/ ²)^1/3 - J current density vector - k (/)^1/2 B ₀ - M k( ²/g)^1/3 - n summation index - q mass flow rate - Q non-dimensional mass flow rate, =q/ - t time - T non-dimensional time, =t(g ²/)^1/3 - Laplace transform ofv(x, t) - V fluid velocity vector, =[0,v(x, t), 0] - (x, y, z) space coordinates - Y non-dimensionaly-coordinate, =y(g/ ²)^1/3 Greek symbols _n (n+1/2) - conductivity - density - kinematic viscosity