Magneto-hydrodynamically lubricated externally pressurised bearing with variable film thickness

Summary This paper studies theoretically the use of a conducting lubricant in an externally pressurised bearing with variable film-thickness in the presence of an axial magnetic field. The flow and other characteristics are determined and it is shown that the pressure and load capacity can be increased by increasing the strength of the applied magnetic field at a given flow rate. But at a given feeding pressure the load capacity and pressure do not depend upon the magnetic field. The load capacity of this bearing is greater than that of a bearing having a constant film-thickness. It is also pointed out that the frictional drag on the rotor can be minimised by supplying electrical energy to the system.Nomenclature angle which the rotor surface makes with the stator (see fig. 1) - angular velocity of the rotor - _t terminal voltage between the electrodes - _t.o.c. open circuit voltage - viscosity of the lubricant - conductivity of the lubricant - B ₀ strength of the applied magnetic field - E _r radial component of the electric field - h variable film-thickness - h ₀ minimum film-thickness - I total current - I _s.c. short circuit current - L depth of the recess - M Hartmann number - p pressure - p _e exit pressure - p _i inlet pressure - Q rate of volume flow - r radial coordinate - R ₀ radius of the recess - R outer radius of the stator - R _i internal resistance - T frictional drag - u radial velocity - v tangential velocity - W load capacity - normalised load capacity - z axial coordinate     

Aspect ratio effect on flow-induced forces on circular cylinders in a cross-flow

Finite-span circular cylinders with two different aspect ratios, placed in a cross-flow, are investigated experimentally at a cylinder Reynolds number of 46,000. Simultaneous measurements of the flow-induced unsteady forces on the cylinders and the stream velocity in the wake are carried out. These results together with mean drag measurements along the span and available literature data are used to evaluate the flow mechanisms responsible for the induced unsteady forces and the effect of aspect ratio on these forces. The coherence of vortex shedding along the span of the cylinder is partially destroyed by the separated flow emanating from the top and by the recirculating flow behind the cylinder. As a result, the fluctuating lift decreases drastically. Based on the data collected, it is conjectured that the fluctuating recirculating flow behind the cylinder is the flow mechanism responsible for the unsteady drag and causes it to increase beyond the fluctuating lift. The fluctuating recirculating flow is a direct consequence of the unsteady separated flow. The unsteady forces vary along the span, with lift increasing and drag decreasing towards the cylinder base. When the cylinder span is large compared to the wall boundary layer thickness, a submerged two-dimensional region exists near the base. As the span decreases, the submerged two-dimensional region becomes smaller and eventually vanishes. Altogether, these results show that fluctuating drag is the dominant unsteady force in finite-span cylinders placed in a cross-flow. Its characteristic frequency is larger than that of the vortex shedding frequency.List of symbols a span of active element on cylinder, = 2.5 cm - C _D local rms drag coefficient, 2D/ U ² da - C _L local rms lift coefficient, 2l/ U ² da - C _D local mean drag coefficient, 2D/ U ² da - C _D spanwise-averaged C _D for finite-span cylinder - (C _D ) _2D spanwise-averaged mean drag coefficient for two-dimensional cylinder - C _p pressured coefficient - -(C _p ) _b pressure coefficient at = - d diameter of cylinder, = 10.2 cm - D fluctuating component of instantaneous drag - D local rms of fluctuating drag - D local mean drag - E _D power spectrum of fluctuating drag, defined as - E _L power spectra of fluctuating lift, defined as - f _D dominant frequency of drag spectrum - f _L dominant frequency of lift spectrum - f _u dominant frequency of velocity spectrum - h span of cylinder - H height of test section, = 30.5 cm - L fluctuating component of instantaneous lift - L local rms of fluctuating lift - R _Du () cross-correlation function of streamwise velocity and local drag, - R _Lu () cross-correlation function of stream wise velocity and local lift, - Re Reynolds number, U d/y - S _L Strouhal number based on f _L ,f _L d/U - S _D Strouhal number based on f _D ,f _D d/U - S _u Strouhal number based on f _u , f _u d/U - t time - u fluctuating component of instantaneous streamwise velocity - U mean streamwise velocity - mean stream velocity upstream of cylinder - x streamwise distance measured from axis of cylinder - y transverse distance measured from axis of test section - z spanwise distance measured from cylinder base - angular position on cylinder circumference measured from forward stagnation - kinematic viscosity of air - density of air - time lag in cross-correlation function - _D normalized spectrum of fluctuating drag - _L normalized spectrum of fluctuating lift     

Reynolds number dependence of the drag coefficient for laminar flow through fine-scale photoetched screens

The laminar steady flow downstream of fine-mesh screens is studied. Instead of woven-wire screens, high-uniformity screens are fabricated by photoetching holes into 50.8 m thick Inconel sheets. The resulting screens have minimum wire widths of 50.8 m and inter-wire separations of 254 m and 318 m for the two screens examined. A flow facility has been constructed for experiments with these screens. Air is passed through the screens at upstream velocities yielding wire width Reynolds numbers from 2 to 35. To determine the drag coefficient, pressure drops across the screens are measured using pressure transducers and manometers. Threedimensional flow simulations are also performed. The computational drag coefficients consistently overpredict the experimental values. However, the computational results exhibit sensitivity to the assumed wire cross section, indicating that detailed knowledge of the wire cross section is essential for unambiguous interpretation of experiments using photoetched screens. Standard semi-empirical drag correlations for woven-wire screens do not predict the present experimental results with consistent accuracy.List of symbols A ₁, A ₂ screen aspect ratios - c _d screen drag coefficient - d woven-wire diameter - D photoetched minimum wire width (spanwise) - f woven-wire screen drag function - M distance between adjacent wires - N spectral-element order - o woven-wire open area fraction - O photoetched open area fraction - _p pressure drop across screen - Re _d woven-wire diameter Reynolds number - Re _D photoetched wire width Reynolds number - U fluid velocity upstream of screen - W photoetched sheet thickness (streamwise) - x, y, z spatial coordinates - fluid density - fluid viscosity     

Rheological behaviour of two-phase flow of solid particles in a gas

The present investigation was concerned with the rheological behaviour of dilute suspensions of solid particles in a gas in a vertical cocurrent flow moving upwards. Starting from the experimentally determined dependence of the pressure drop on the concentration of solid particles and the Reynolds number of the carrier medium in the steady flow region, the rheological parameters were estimated using pseudo-shear diagrams. Air was the carrier medium and the dispersed phase was one of six fractions of polypropylene powder and five fractions of glass ballotini. The results show that the investigated two-phase systems have pseudoplastic character which becomes more pronounced with increases in concentration, equivalent diameter and density of solid particles in the flowing suspension. C _d coefficient of particle resistance - d _e equivalent diameter of particles - D column diameter - Fr Froude number - g gravitational acceleration - K rheological parameter - L length - n rheological parameter - p _t pressure drop due to friction - p _m total pressure drop - p _ag pressure drop due to acceleration of the gas phase - p _as pressure drop due to acceleration of the solid phase - p _g hydrostatic pressure of the gas phase - p _s specific effective weight of the dispersed phase - r radius - Re Reynolds number - Re _p Reynolds number of a particle - Re _G generalized Reynolds number - Re _G1 generalized Reynolds number relating to the end of the laminar flow region - Re _G2 generalized Reynolds number relating to the beginning of the turbulent flow region - w _z axial component of velocity - u _t steady free-fall velocity of a single particle - w average velocity - w _g average velocity of the gas phase - w _s average velocity of the dispersed phase of solid particles - relative mass fraction of solid particles - x _s volume fraction of solid particles - _g coefficient of pressure drop due to friction - µ dynamic viscosity - _g density of the gas phase - _m density of the suspension - _s density of solid particles - _ds density of the dispersed phase - _w shear stress at the wall     

Thermal and inertia effects in externally pressurized conical thrust oil bearings

This paper presents theoretical and experimental investigations of thermal and inertia effects on the performance of externally pressurized conical thrust bearings. The analysis, as well as the experimental results, revealed that the increase in oil temperature due to pad rotation has a detrimental effect on the load carrying capacity, while it increases the flow rate. Increasing the speed of rotation, will increase or decrease the bearing load carrying capacity depending on the recess dimensions.Nomenclature c lubricant specific heat - F frictional torque - h film thickness - L load carrying capacity - P pressure - P pressure ratio (P/P ₁) - P ₁ inlet pressure - Q volume flow rate - r radius measured on cone surface - r radius ratio (R/R ₃) - R ₁ supply hole radius - R ₂ recess radius - R ₃ outside radius of bearing - S inertia parameter (0.15 ² R ₃ ² /P ₁) - T temperature - u, v, w velocity components (see Fig. 2) - z coordinate normal to cone surface - lubricant density - lubricant viscosity - 2 cone apex angle - rotational speed - recess depth     

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     

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.     

Viscoelastic swirling flow with free surface in cylindrical chambers

The flow of a viscoelastic liquid driven by the steadily rotating bottom cover of a cylindrical cup is investigated. The flow field and the shape of the free surface are determined at the lowest significant orders of the regular domain perturbation in terms of the angular velocity of the bottom cap. The meridional field superposed on a primary azimuthal field shows a structure of multiple cells. The velocity field and the shape of the free surface are strongly effected by the cylinder aspect ratio and the elasticity of the liquid. The use of this flow configuration as a free surface rheometer to determine the first two Rivlin-Ericksen constants is shown to be promising.Nomenclature R, ,Z Coordinates in the physical domain D - , , Coordinates in the rest stateD ₀ - r, ,z Dimensionless coordinates in the rest stateD ₀ - Angular velocity - Zero shear viscosity - Surface tension coefficient - Density - Dimensionless surface tension parameter - ₁, ₂ The first two Rivlin-Ericksen constants - Stream function - Dimensionless second order meridional stream function - ^* Dimensionless second normal stress function - 2 Dimensionless sum of the first and second normal stress functions - N ₁,N ₂ The first and second normal stress functions - n Unit normal vector - D Stretching tensor - A _n nth order Rivlin-Ericksen tensor - S Extra-stress - u Velocity field - U Dimensionless second order meridional velocity field - V Dimensionless first order azimuthal velocity field - p Pressure - Modified pressure field - P Dimensionless second order pressure field - J Mean curvature - a Cylinder radius - d Liquid depth at rest - D Dimensionless liquid depth at rest - h Free surface height - H Dimensionless free surface height at the second order     

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²     

Buoyancy- and pressure-driven motion in a vertical porous layer: Effects of quadratic drag

The seepage velocity arising from pressure and buoyancy driving forces in a slender vertical layer of fluid-saturated porous media is considered. Quadratic drag (Forcheimer effects) and Brinkman viscous forces are included in the analysis. Parameters are identified which characterize the influence of matrix permeability, quadratic drag and buoyancy. An explicit solution is obtained for pressure-driven flow which illustrates the influence of quadratic drag and the strong boundary layer behavior expected for low permeability media. The experimental data of Givler and Altobelli [2] for water seepage through a high porosity foam is found to yield good agreement with the present analysis. For the case of buoyancy-driven flow, a uniformly valid approximate solution is found for low permeability media. Comparison with the pressure-driven case shows strong similarities in the near-wall region.Nomenclature B function of - d layer thickness - D discriminant defined by Equation (9) - modified Darcy number - F Forcheimer constant - g gravitational acceleration - k porous matrix permeability - m parameter defined by Equation (11) - p pressure - p modified pressure - pressure gradient - R buoyancy parameter - T ₀ nominal layer temperature - u seepage velocity - dimensionless seepage velocity - _c composite approximation - _i boundary layer velocity - _o outer or core flow approximation - _m midplane velocity - U matching velocity - V cross-sectional average velocity - w variable defined by Equation (12) - x, z Cartesian coordinates - , dimensionless Cartesian coordinates - inertia parameter - T layer temperature difference - larger root of cubic given by Equation (8) - fluid dynamic viscosity - _e effective viscosity of fluid saturated medium - variable defined by Equation (18) - 0 fluid density - smaller root of cubic given by Equation (8) - variable defined by Equation (18) - stretched inner coordinate - porosity - function of      

Effect of angular inertia of lubricant in MHD hydrostatic thrust bearings

Summary Circumferential motion of a conducting lubricant in a hydrostatic thrust bearing is caused either by the angular motion of a rotating disk or by the interaction of a radial electric field and an axial magnetic field. Under the assumption that the fluid inertia due to radial motion is negligibly small in comparison with that due to angular motion, it is found analytically that the rotor causes an increase in flow rate and a decrease in load capacity, while both are increased by the application of an electric field in the presence of an axial magnetic field. The critical angular speed of the rotor at which the bearing can no longer support any load is obtained, and the possibility of flow separation in the lubricant is discussed.Nomenclature a recess radius - b outside disk radius - B ₀ magnetic induction of uniform axial magnetic field - E ₀ radial electric field at r=a - E _r radial electric field - h half of lubricant film thickness - M Hartmann number = (B ₀ ² h ²/)^1/2 - P pressure - P ₀ pressure at r=a - P _e pressure at r=b - Q volume flow rate of lubricant - Q ₀ flow rate of a nonrotating bearing without magnetic field - r radial coordinate - r _s position of flow separation on stationary disk - u, v fluid velocity components in radial and circumferential directions, respectively - W load carrying capacity of bearing - W ₀ load capacity of a nonrotating bearing without magnetic field - z axial coordinate - coefficient of viscosity - _e magnetic permeability - fluid density - electrical conductivity - electric potential - angular speed of rotating disk - _c critical rotor speed at which W=0     

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

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

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      

The chaotic response of a fluttering panel: The influence of maneuvering

The influence of maneuvering on the chaotic response of a fluttering buckled plate on an aircraft has been studied. The governing equations, derived using Lagrangian mechanics, include geometric non-linearities associated with the occurrence of tensile stresses, as well as coupling between the angular velocity of the maneuver and the elastic degrees of freedom. Numerical simulation for periodic and chaotic responses are conducted in order to analyze the influence of the pull-up maneuver on the dynamic behavior of the panel. Long-time histories phase-plane plots, and power spectra of the responses are presented. As the maneuver (load factor) increases, the system exhibits complicated dynamic behavior including a direct and inverse cascade of subharmonic bifurcations, intermittency, and chaos. Beside these classical routes of transition from a periodic state to chaos, our calculations suggest amplitude modulation as a possible new mode of transition to chaos. Consequently this research contributes to the understanding of the mechanisms through which the transition between periodic and strange attractors occurs in, dissipative mechanical systems. In the case of a prescribed time dependent maneuver, a remarkable transition between the different types of limit cycles is presented.Nomenclature a plate length - a _r u _r /h - D plate bending stiffness - E modulus of elasticity - g acceleration due to gravity - h plate thickness - j₁,j₂,j₃ base vectors of the body frame of reference - K spring constant - M Mach number - n 1 + ₀/g - N ₁ applied in-plane force - p–p aerodynamic pressure - P pa ⁴/Dh - q ₀/2 - Q _r generalized Lagrangian forces - R rotation matrix - R ₄ N, a ²/D - t time - kinetic energy - u plate deflection - u displacement of the structure - u _r modal amplitude - v₀ velocity - x coordinates in the inertial frame of reference - z coordinates in the body frame of reference - Ka/(Ka+Eh) - - elastic energy - 2qa ³/D - a/_mh - Poisson's ratio - material coordinates - air density - _m plate density - - _r prescribed functions - _r sin(r z/a) - angular velocity - a/v₀ - skew-symmetric matrix form of the angular velocity     

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.     

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     

Reynolds number dependence of the drag coefficient for laminar flow through electrically-heated photoetched screens

Experiments are performed to measure the drag coefficient of electrically-heated screens. Square-pattern 80 mesh and 100 mesh screens of 50.8 m-wide wires photoetched from 50.8 m thick Inconel sheets are examined. Ambient air is passed through these screens at upstream velocities yielding wire-width Reynolds numbers from 2 to 35, and electrical current is passed through the screens to generate heat fluxes from o to 0.17 MW/m², based on the total screen area. The dependence of the drag coefficient on Reynolds number and heat flux is determined for these two screens by measuring pressure drops across the screens for a variety of conditions in these ranges. In all cases, heating is found to increase the drag coefficient above the unheated value. A correlation relating the heated drag coefficient to the unheated drag coefficient is developed based on the idea that the main effect of heating at these levels is to modify the Reynolds number through modifying the viscosity. This correlation is seen to reproduce the experimental results closely.List of Symbols A total screen cross sectional area - C fitting coefficient, near unity - c _D heated drag coefficient - c _{D, 0} unheated drag coefficient - C _p air specific heat at constant pressure - D photoetched wire width, sheet thickness - h _s stagnation point heat-transfer coefficient - k air thermal conductivity - M distance between adjacent wires - O open area fraction - p air pressure - p air pressure drop across screen - Pr Prandtl number for air, c _p/k - Q total electrical power to screen - R radius of curvature at stagnation point - Re _D wire width Reynolds number, UD/ - T air temperature - U air speed upstream of screen - air specific heat ratio - air density - air viscosity - exponent in temperature power law for viscosity - () quantity () evaluated at heated screen temperature The authors thank John Lewin and Bob Meyer for their assistance in the design and fabrication of the heated screen test facility and Tom Grasser for his help in performing the experiments. This work was performed at Sandia National laboratories, supported by the U.S. Department of Energy under contract number DE-AC04-94AL85000.     

Neutrally buoyant bubbles used as flow tracers in air

Research has been performed to determine the accuracy of neutrally buoyant and near-neutrally-buoyant bubbles used as flow tracers in an incompressible potential flowfield. Experimental and computational results are presented to evaluate the quantitative accuracy of neutrally buoyant bubbles using a commercially available helium bubble generation system. A two-dimensional experiment was conducted to determine actual bubble trajectories in the stagnation region of a NACA 0012 airfoil at 0° angle of attack. A computational scheme evaluating the equation of motion for a single bubble was also used to determine the factors which affect a bubble's trajectory. The theoretical and computational analysis have shown that neutrally buoyant bubbles will trace complex flow patterns faithfully in the flowfield of interest. Experimental analysis revealed that the use of bubbles generated by the commercially available system to trace flow patterns should be limited to qualitative measurements unless care is taken to ensure neutral buoyancy.Nomenclature a _c centripetal acceleration - c model chord - c _D bubble drag coefficient - D bubble diameter - g acceleration due to gravity - g _v acceleration due to gravity vector - h trajectory deviation normalization parameter - K nondimensional inertia parameter, - m _f mass of fluid - m _p mass of bubble - p static pressure - r radial distance, bubble radius - R gas constant - Re free-stream Reynolds number, - Re _p bubble slip Reynolds number, - S cross-sectional area of sphere - T temperature - t time - u streamwise velocity component - U free-stream velocity - v _f fluid velocity vector - V _p bubble velocity vector - x _p bubble position vector - y _b bubble trajectory y/c - y _s streamline y/c - model angle of attack - bubble solution surface tension - potential vortex strength - _bfs bubble solution density - fluid density - bubble density - bubble wall thickness - fluid viscosity     

An experimental investigation of separating flow on a convex surface

The mean and turbulent characteristics of an incompressible turbulent boundary layer developing on a convex surface under the influence of an adverse pressure gradient are presented in this paper.The turbulence quantities measured include all the components of Reynolds stresses, auto-correlation functions and power spectra of the three components of turbulence. The results indicate the comparative influence of the convex curvature and adverse pressure gradient which are simultaneously acting on the flow. The investigation provides extensive experimental information which is much needed for a better understanding of turbulent shear flows.Nomenclature a, b constants in equation for velocity defect profile (Fig. 6) - c _f skin-friction coefficient (= _w/F _1/2 U ₁ ² ) - E(k ₁) one-dimensional wave number spectra - f frequency in Hz - G Clauser's equilibrium parameter = (H–1)/H(c _f /2) - H shape parameter (= ₁/ ₂) - k ₁ wave number (=2f/U) - L _u, L _v, L _w length scales of u, v and w fluctuations - p _s static pressure on the measurement surface - p _w reference tunnel wall static pressure - q ² total turbulent kinetic energy - R radius of curvature of the convex surface - R() auto-correlation function - T _u, T _v, T _w time scales of u, v and w fluctuations - U local mean velocity - U ₁ local free stream velocity - U _* friction velocity - u, v, w velocity fluctuations in x, y and z directions respectively - X streamwise coordinate measured along the surface from A (Fig. 1b) - x streamwise coordinate measured along the surface reckoned from station 9 - y coordinate normal to the surface - z spanwise coordinate - ₁/ _w · dp/dx - - boundary layer thickness - ₁ displacement thickness - ₂ momentum thickness - ₃ energy thickness - kinematic viscosity - density - time delay - _w wall shear stress     

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)