The virtual origin of riblets in an open channel flow

In this paper, a method using the mean velocity profiles for the buffer layer was developed for the estimation of the virtual origin over a riblets surface in an open channel flow. First, the standardized profiles of the mixing length were estimated from the velocity measurement in the inner layer, and the location of the edge of the viscous layer was obtained. Then, the virtual origins were estimated by the best match between the measured velocity profile and the equations of the velocity profile derived from the mixing length profiles. It was made clear that the virtual origin and the thickness of the viscous layer are the function of the roughness Reynolds number. The drag variation coincided well with other results.Nomenclature f _r skin friction coefficient - f _ro skin friction coefficient in smooth channel at the same flow quantity and the same energy slope - g gravity acceleration - H water depth from virtual origin to water surface - H ⁺ u*H/ - H false water depth from top of riblets to water surface - H ⁺ u*H/ - I _e streamwise energy slope - I _b bed slope - k riblet height - k ⁺ u*k/ - l mixing length - l _s standardized mixing length - Q flow quantity - Re Reynolds number volume flow/unit width/v - s riblet spacing - u mean velocity - u* friction velocity = - u* false friction velocity = - y distance from virtual origin - y distance from top of riblet - y ₀ distance from top of riblet to virtual origin - y _v distance from top of riblet to edge of viscous layer - y ⁺ u*y/ - y ⁺ u*y/ - y ₀ ⁺ u*y ₀/ - u ⁺ u*y/ - shifting coefficient for standardization - thickness of viscous layer=y ₀+y - ⁺ u*/ - ⁺ u*/ - eddy viscosity - ridge angle - v kinematic viscosity - density - shear stress     

Laminar source flow between parallel porous disks

An analysis is presented for laminar source flow between infinite parallel porous disks. The solution is in the form of a perturbation from the creeping flow solution. Expressions for the velocity, pressure, and shear stress are obtained and compared with the results based on the assumption of creeping flow.Nomenclature a half distance between disks - radial coordinate - r dimensionless radial coordinate, /a - axial coordinate - z dimensionless axial coordinate, /a - radial coordinate of a point in the flow - R dimensionless radial coordinate of a point in the flow, /a - velocity component in radial direction - u =a/, dimensionless velocity component in radial direction - velocity component in axial direction - v = a/}, dimensionless velocity component in axial direction - static pressure - p = (a ²/ ²), dimensionless static pressure - =p(r, z)–p(R, z), dimensionless pressure drop - V magnitude of suction or injection velocity - Q volumetric flow rate of the source - Re source Reynolds number, Q/4a - reduced Reynolds number, Re/r ² - critical Reynolds number - R _w wall Reynolds number, Va/ - viscosity - density - =/, kinematic viscosity - shear stress at upper disk - ₀ = (a ²/ ²), dimensionless shear stress at upper disk - shear stress ratio, ₀/( ₀)_inertialess - u = , dimensionless average radial velocity - u/u, ratio of radial velocity to average radial velocity - dimensionless stream function     

Turbulent boundary layer over a smooth wall with widely separated transverse square cavities

Laser-Doppler velocimetry (LDV) measurements and flow visualizations are used to study a turbulent boundary layer over a smooth wall with transverse square cavities at two values of the momentum thickness Reynolds number (R =400 and 1300). The cavities are spaced 20 element widths apart in the streamwise direction. Flow visualizations reveal a significant communication between the cavities and the overlying shear layer, with frequent inflows and ejections of fluid to and from cavities. There is evidence to suggest that quasi-streamwise near-wall vortices are responsible for the ejections of fluid out of the cavities. The wall shear stress, which is measured accurately, increases sharply immediately downstream of the cavity. This increase is followed by a sudden decrease and a slower return to the smooth wall value. Integration of the wall shear stress in the streamwise direction indicates that there is an increase in drag of 3.4% at bothR .Nomenclature C _f skin friction coefficient - C _fsw friction coefficient for a continuous smooth wall - k height of the cavity - k ⁺ ku / - R Reynolds number based on momentum thickness (U ₁ /v) - Rx Reynolds number based on streamwise distance (U ₁ x/) - s streamwise distance between two cavities - t time - t ⁺ tu ₂ / - U ₁ freestream velocity - mean velocity inx direction - u,v,w rms turbulent intensities inx,y andz directions - u local friction velocity - u _sw friction velocity for a continuous smooth wall - w width of the cavity - x streamwise co-ordinate measured from the downstream edge of the cavity - y co-ordinate normal to the wall - z spanwise co-ordinate - y ⁺ yu / - boundary layer thickness - ₀ boundary layer thickness near the upstream edge of the cavity - _i thickness of internal layer - kinematic viscosity of water - ⁺ zu / - momentum thickness     

Analysis of suddenly started laminar flow in the entrance region of a circular tube

Suddenly started laminar flow in the entrance region of a circular tube, with constant inlet velocity, is investigated analytically by using integral momentum approach. A closed form solution to the integral momentum equation is obtained by the method of characteristics to determine boundary layer thickness, entrance length, velocity profile, and pressure gradient.Nomenclature M(, , ) a function - N(, , ) a function - p pressure - p* p/1/2U ², dimensionless pressure - Q(, , ) a function - R radius of the tube - r radial distance - Re 2RU/, Reynolds number - t time - U inlet velocity, constant for all time, uniform over the cross section - u velocity in the boundary layer - u* u/U, dimensionless velocity - u ₁ velocity in the inviscid core - x axial distance - y distance perpendicular to the axis of the tube - y* y/R, dimensionless distance perpendicular to the axis - boundary layer thickness - * displacement thickness - /R, dimensionless boundary layer thickness - momentum thickness - absolute viscosity of the fluid - /, kinematic viscosity of the fluid - x/(R Re), dimensionless axial distance - density of the fluid - tU/(R Re), dimensionless time - _w wall shear stress     

A method for the analysis of recirculating turbulent flows in a cavity

The complex fluid-dynamic aspects of a turbulent recirculating flow in a cavity with axial throughflow, and a rotating wall, were investigated by adopting a simple procedure for evaluating the turbulent stresses. The flow field was divided into two regions, a core and a wall region respectively. A wall function was adopted in the zones near to the solid boundaries, while a constant eddy diffusivity was assumed, in the core, following the indications of computed heat transfer coefficients in comparison with existing experimental data. The distributions of the stream function and of the tangential velocity are presented for a range of the rotational Reynolds number of the rotating wall and of the Reynolds number of the throughflow.

---

Turbulente Rezirkulationsströmung in einem Hohlraum

Zusammenfassung Die komplizierten fluiddynamischen Aspekte einer turbulenten Rezirkulationsströmung in einem Hohlraum mit axialem Durchfluß und einer rotierenden Wand werden unter Verwendung einer vereinfachten Methode zur Berechnung der turbulenten Spannungen betrachtet. Das Strömungsfeld wird in einen Kern und einen Wandbereich aufgeteilt. Für die wandnahen Zonen wird eine Wandfunktion angenommen, während im Kern mit konstanter Wirbeldiffusivität gerechnet wird, was durch den Vergleich berechneter mit gemessenen Wärmeübergangskoeffizienten gerechtfertigt erscheint. Verteilungen der Stromfunktion und der tangentialen Geschwindigkeit sind für einen bestimmten Bereich der Reynoldszahlen für die Wandrotation und der für den Durchfluß angegeben.

---

Nomenclature L axial length of enclosure - P dimensionless pressure, p^*2 - p static pressure - R dimensionless radial coordinate, r/r^* - r radial coordinate - r^* reference length, equal to r_O for enclosure - r_i radii of inlet and exit apertures - Re Reynolds number, v^*r^*/ - Re_i pipe Reynolds number, ¯v_zi(2r_i)/ - Re_t turbulent Reynolds number, Re(/) - Re rotational Reynolds number, r ₀ ² / - t dimensionless time,t/(r^*/v^*) - t time - V_r, V, Vz dimensionless velocity components, V_r/v^*, v, v_z/v^* - v_i turbulent fluctuation of the i-component of velocity - v_r, v, v_z velocity components - v^* reference velocity, equal to ¯v_zi for enclosure - X coordinate along a wall, x/r^* - Y coordinate normal to a wall, y/r^* - Z dimensionless axial coordinate, z/r^* - z axial coordinate - eddy diffusivity for momentum - dynamic viscosity - kinematic viscosity - density - shear stress - dimensionless shear stress, /v^*2 - dimensionless stream function, /r*²v*² - stream function - angular velocity - tangential vorticity component - ()_eff effective - ()_l laminar - ()_t turbulent - mean over the time     

Thermal free convection from a solid sphere

Summary Thermal free convection from a sphere has been studied by melting solid benzene spheres in excess liquid benzene (Pr=8,3; 10⁸<Gr<10⁹). Overall heat transfer as well as local heat transfer were investigated. For the effect of cold liquid produced by the melting a correction has been applied. Results are compared with those obtained by other workers who used alternative experimental methods.Nomenclature coefficient of heat transfer - d characteristic length, here diameter of sphere - thermal conductivity - g acceleration of free fall - cubic expansion coefficient - T temperature difference between wall and fluid at infinity - kinematic viscosity - density - c specific heat capacity - a thermal diffusivity (=/c) - D diffusion coefficient - Nu dimensionless Nusselt number (=d/) - Nu* the analogous number for mass transfer (=kd/D) - mean value of Nusselt number - Gr dimensionless Grashof number (=gd ³T/ ²) - Gr* the analogous number for mass transfer (=gd ³x/ ²) - Pr dimensionless Prandtl number (=/a) - Sc dimensionless Schmidt number (=/D)     

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      

Über den Einfluß von Querkrümmung und Dichte bei inhomogenen turbulenten Freistrahlen

Zusammenfassung Zur Berechnung turbulenter Strömungen wird das k--Modell im Ansatz für die turbulente Scheinzähigkeit erweitert, so daß es den Querkrümmungs- und Dichteeinfluß auf den turbulenten Transportaustausch erfaßt. Die dabei zu bestimmenden Konstanten werden derart festgelegt, daß die bestmögliche Übereinstimmung zwischen Berechnung und Messung erzielt wird. Die numerische Integration der Grenzschichtgleichungen erfolgt unter Verwendung einer Transformation mit dem Differenzenverfahren vom Hermiteschen Typ. Das erweiterte Modell wird auf rotationssymmetrische Freistrahlen veränderlicher Dichte angewendet und zeigt Übereinstimmung zwischen Rechnung und Experiment.

---

On the influence of transvers-curvature and density in inhomogeneous turbulent free jets

The prediction of turbulent flows based on the k- model is extended to include the influence of transverse-curvature and density on the turbulent transport mechanisms. The empirical constants involved are adjusted such that the best agreement between predictions and experimental results is obtained. Using a transformation the boundary layer equations are solved numerically by means of a finite difference method of Hermitian type. The extended model is applied to predict the axisymmetric jet with variable density. The results of the calculations are in agreement with measurements.

---

Bezeichnungen Wirbelabsorptionskoeffizient - c_i Massenkonzentration der Komponente i - c_D, c_L, c, c₁, c₂ Konstanten des Turbulenzmodells - d Düsendurchmesser - E bezogene Dissipationsrate - f bezogene Stromfunktion - f Korrekturfunktion für die turbulente Scheinzähigkeit - j turbulenter Diffusionsstrom - k Turbulenzenergie - k_i Schrittweite in -Richtung - K dimensionslose Turbulenzenergie - L turbulentes Längenmaß - M_i Molmasse der Komponente i - p Druck - allgemeine Gaskonstante - r Querkoordinate - r_0,5 Halbwertsbreite der Geschwindigkeit - r_0,5c Halbwertsbreite der Konzentration - T Temperatur - u Geschwindigkeitskomponente in x-Richtung - v Geschwindigkeitskomponente in r-Richtung - x Längskoordinate - y allgemeine Funktion - Y_i diskreter Wert der Funktion y - Relaxationsfaktor für Iteration - turbulente Dissipationsrate - transformierte r-Koordinate - kinematische Zähigkeit - Exponent - transformierte x-Koordinate - Dichte - _k, Konstanten des Turbulenzmodells - Schubspannung - allgemeine Variable - Stromfunktion - Turbulente Transportgröße Indizes 0 Strahlanfang - m auf der Achse - r mit Berücksichtigung der Krümmung - t turbulent - mit Berücksichtigung der Dichte - im Unendlichen - Schwankungswert oder Ableitung einer Funktion - – Mittelwert Herrn Professor Dr.-Ing. R. Günther zum 70. Geburtstag gewidmet     

Boundary layer effect on thermal NO decomposition behind incident shock waves

The thermal decomposition of nitric oxide (diluted in Argon) has been measured behind incident shock waves by means of IR diode laser absorption spectroscopy. In two independent runs the diode laser was tuned to the=0 =1² _3/2 R(18.5)-rotational vibrational transition and the=1 =2² _3/2 R(20.5)-rotational vibrational transition of nitric oxide, respectively. These two transitions originating from the vibrational ground state (=0) and the first excited vibrational state (=1) were selected in order to probe the homogeneity along the absorption path. The measured NO decomposition could satisfactorily be described by a chemical reaction mechanism after taking into account boundary layer corrections according to the theory of Mirels. The study forms a further proof of Mirels' theory including his prediction of the laminar-turbulent transition. It also shows, that the inhomogeneities from the boundary layer do not affect the IR linear absorption markedly.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Effects of blowing or suction on the laminar boundary layer flow over a rotating cone

The effects of suction or blowing at the surface of a rotating cone in a quiet fluid on the skin friction and heat transfer are described. The equations which govern the fluid motion and thermal energy transfer are transformed by the boundary layer approximations and the resulting equations are solved under the condition that the suction or blowing velocity varies as x ⁿ (x: distance measured from the apex of the cone, n: arbitrary constant). The solutions are obtained as a perturbation from the basic laminar flow of an incompressible viscous fluid over the impermeable rotating cone. Detailed numerical calculations are performed for the case of an isothermal rotating cone with uniform blowing or suction, i.e. n=0, the Prandtl number being 0.72. Results are given for the shear stress, heat transfer and velocity and temperature fields. It is shown from the analysis that suction sharply increases the circumferential shear stress and the heat transfer at the surface.Nomenclature c proportional constant - C _fx dimensionless skin friction factor, _x /(V ²) - C _fx0 dimensionless skin friction factor for an impermeable cone - C _fy dimensionless circumferential skin friction factor, _y /(V ²) - C _fy0 dimensionless circumferential skin friction factor for an impermeable cone - c _p specific heat at constant pressure - f _k function of - g _k function of - h heat transfer coefficient, q/(T _w–T ) - k thermal conductivity of fluid - n arbitrary constant - Nu _x local Nusselt number, hx/k - Nu _x0 local Nusselt number for an impermeable cone - Pr Prandtl number - q heat transfer rate - r radius of a circular cross section of the cone, x sin - R _x Reynolds number, Vx/ - T temperature - T _w surface temperature of the cone - T temperature of the surrounding fluid - u fluid velocity in x-direction - v fluid velocity in y-direction - V circumferential velocity at the cone surface, r - w fluid velocity in z-direction - x coordinate along meridional section - y coordinate along a circular cross section - z coordinate perpendicular to both x and y - perturbation parameter, cx ⁿ /(x sin )^1/2 - dimensionless z-coordinate, z( sin /)^1/2 - _k function of - kinematic viscosity - density of fluid - _x skin friction in x-direction - _y circumferential skin friction - stream function - angular speed of the cone     

A note on thermal stresses in hollow cylinders

Conclusions The extensions of Trostel's solutions derived in this paper may be employed tor general over wide conditions with a resulting error less than 3%, the stress values being too small in magnitude by this amount. The error decreases as the variation of physical properties of the media decreases or as 1/2. M. M. Stanii, Lectures in Mathematical Elasticity during summer semester 1958, Purdue University.     

Flow of second-order fluids over an enclosed rotating disc

Summary This paper is devoted to a study of the flow of a second-order fluid (flowing with a small mass rate of symmetrical radial outflow m, taken negative for a net radial inflow) over a finite rotating disc enclosed within a coaxial cylinderical casing. The effects of the second-order terms are observed to depend upon two dimensionless parameters ₁ and ₂. Maximum values ₁ and ₂ of the dimensionless radial distances at which there is no recirculation, for the cases of net radial outflow (m>0) and net radial inflow (m<0) respectively, decrease with an increase in the second-order effects [represented by T(=₁+₂)]. The velocities at ₁ and ₂ as well as at some other fixed radii have been calculated for different T and the associated phenomena of no-recirculation/recirculation discussed. The change in flow phenomena due to a reversal of the direction of net radial flow has also been studied. The moment on the rotating disc increases with T.Nomenclature , , z coordinates in a cylindrical polar system - z ₀ distance between rotor and stator (gap length) - =/z ₀, dimensionless radial distance - =z/z ₀, dimensionless axial distance - _s = _s/z₀, dimensionless disc radius - V =(u, v, w), velocity vector - dimensionless velocity components - uniform angular velocity of the rotor - , p fluid density and pressure - P =p/(₂ z ₀₂ ² , dimensionless pressure - ₁, ₂, ₃ kinematic coefficients of Newtonian viscosity, elastico-viscosity and cross-viscosity respectively - ₁, _{2 2}/z ₀ ² , resp. ₃/z ₀ ² , dimensionless parameters representing the ratio of second-order and inertial effects - m = , mass rate of symmetrical radial outflow - l a number associated with induced circulatory flow - Rm =m/(z ₀₁), Reynolds number of radial outflow - R _l =l/(z ₀₁), Reynolds number of induced circulatory flow - Rz =z ₀ ² /₁, Reynolds number based on the gap - ₁, ₂ maximum radii at which there is no recirculation for the cases Rm>0 and Rm<0 respectively - ₁(T), ₂(T) ₁ and ₂ for different T - U _1(T) ⁽⁺⁾ = dimensionless radial velocity, Rm>0 - V _1(T) ⁽⁺⁾ = , dimensionless transverse velocity, Rm>0 - U _2(T) ^(–) = , dimensionless radial velocity, Rm=–Rn<0, m=–n - V _2(T) ^(–) = , dimensionless transverse velocity, Rm<0 - C _m moment coefficient     

Heat transfer from laminar flow in ducts with elliptic cross-section

Summary The cooling of a hot fluid in laminar Newtonian flow through cooled elliptic tubes has been calculated theoretically. Numerical data have been computed for the two values 1.25 and 4 of the axial ratio of the elliptic cross-section . For =1.25 the influence of non-zero thermal resistance between outmost fluid layer and isothermal surroundings has also been investigated. Special attention has been given to the distribution of heat flux around the perimeter; when increases the flux varies more with the position at the circumference. This positional dependence becomes less pronounced, however, as the (position-independent) thermal resistance of the wall increases.Flattening of the conduit, while maintaining its cross-sectional area constant, improves the cooling. Comparison with rectangular pipes shows that this improvement is not as marked with elliptic as with rectangular pipes.Nomenclature A _k =A _{m, n} coefficients of expansion (6) - a, b half-axes of ellipse, b<a - a _p =a _{r, s} coefficients of representation (V) - D hydraulic diameter, = 4S/P; S = cross-sectional area, P = perimeter - D _e equivalent diameter, according to (13) - n coordinate (outward) normal to the tube wall - T temperature of fluid - T _i temperature of fluid at the inlet - T _s temperature of surroundings - v ₀ mean velocity of fluid - v _z longitudinal velocity of fluid - x, y carthesian coordinates coinciding with axes of ellipse - z coordinate in flow direction - , dimensionless half-axes of ellipse, =a/D and =b/D - _t heat transfer coefficient from fluid at bulk temperature to surroundings; equation (11) - _w heat transfer coefficient at the wall; equation (3) - axial ratio of ellipse, = a/b = / - , , , dimensionless coordinates; =x/D, =y/D, =z/D, =n/D - dimensionless temperature, = (T–T _s)/(T _i–T _s) - ₀ cup-mixing mean value of ; equation (10) - thermal conductivity of fluid - _m,n = _k eigenvalue - c volumetric heat capacity of fluid - _{m, n} = _k = _k eigenfunction; equations (6) and (I) - Nu total Nusselt number, = _t D/ - Nusselt number at large distance from the inlet - Nu _w wall Nusselt number, = _w D/, based on _w - Pé Péclet number, = ₀ Dc/     

A mixing length model for turbulent flow in constant cross section ducts

In this paper, a study is made of the damping influence of the wall on turbulent fluid flow. By considering the oscillation of the whole of the boundary, van Driest's original hypothesis has been extended to obtain the wall damping factor in flow in a duct of constant cross section. The damping factor is used in conjunction with mixing length expressions to obtain the velocity field. Particular examples considered are plane parallel flow and axisymmetric flow in a pipe and in an annulus.

---

Ein Modell für die Mischungslänge von turbulenten Strömungen in Rohren mit konstantem Querschnitt

Zusammenfassung In dieser Arbeit wurde der dämpfende Wandeinfluß in turbulenten Strömungen untersucht. Unter Berücksichtigung der Schwingungen in der gesamten Grenzschicht wurde die ursprüngliche Theorie von van Driest erweitert und ein Dämpfungsfaktor an der Wand in Rohrströmungen mit konstantem Querschnitt ermittelt. Dieser Dämpfungsfaktor diente in Verbindung mit Ausdrücken für die Mischungslänge zur Bestimmung des Geschwindigkeitsfeldes. Ausgewählte Beispiele waren die ebene Parallelströmung sowie die Zylinderströmung in einem Rohr und einem Ringspalt.

---

Nomenclature A, A^*(=Au/v) Parameter defined in text - b, b^*(=bu/v) semi-width of parallel plate channel - c(= 1/A) parameter defined in text - E[, /2] complete elliptic integral of the second kind - d damping factor - F, G, H functions - l, l^*(=/v) mixing length - M_O, _O functions - r, r^*(=ru/v) radius - A real part of function - R, S, T, U functions - u, u^*(=u/u) velocity in flow direction Z - friction velocity - x, y, z co-ordinates (z in flow direction) - y^*(=yu/v) non-dimensional wall distance - fluid density - , _eff kinematic viscosity, effective kinematic viscosity - phase angle, or polar coordinate angle - shear stress - (=r/r_W) radius ratio - angular velocity Suffixes w wall value - far from a wall     

Laminar source flow between parallel porous disks with equal suction and injection

An analysis is presented for laminar source flow between parallel stationary porous disks with suction at one of the disks and equal injection at the other. The solution is in the form of an infinite series expansion about the solution at infinite radius, and is valid for all suction and injection rates. Expressions for the velocity, pressure, and shear stress are presented and the effect of the cross flow is discussed.Nomenclature a distance between disks - A, B, ..., J functions of R _w only - F static pressure - p dimensionless static pressure, p(a ²/ ²) - Q volumetric flow rate of the source - r radial coordinate - r dimensionless radial coordinate, r/a - R radial coordinate of a point in the flow region - R dimensionless radial coordinate of a point in the flow region, R - Re source Reynolds number, Q/2a - R _w wall Reynolds number, Va/ - reduced Reynolds number, Re/r ² - critical Reynolds number - velocity component in radial direction - u dimensionless velocity component in radial direction, a/ - average radial velocity, Q/2a - u dimensionless average radial velocity, Re/r - ratio of radial velocity to average radial velocity, u/u - velocity component in axial direction - v dimensionless velocity component in axial direction, v - V magnitude of suction or injection velocity - z axial coordinate - z dimensionless axial coordinate, z a - viscosity - density - kinematic viscosity, / - shear stress at lower disk - shear stress at upper disk - ₀ dimensionless shear stress at lower disk, - ₁ dimensionless shear stress at upper disk, - dimensionless stream function     

The effect of non-steady wall shear on pressure surges in conduits carrying viscous liquids

Summary A study is made of the attenuation of pressure surges in a two-dimension a channel carrying a viscous liquid when a valve at the downstream end is suddenly closed. The analysis differs from previous work in this area by taking into account the transient nature of the wall shear, which in the past has been assumed as equivalent to that existing in steady flow. For large values of the frictional resistance parameter the transient wall shear analysis results in less attenuation than given by the steady wall shear assumption.Nomenclature c /, ft/sec - e base of natural logarithms - F(x, y) integration function, equation (38) - (x) mean value of F(x, y) - g local acceleration of gravity, ft/sec² - h width of conduit, ft - k (2m–1)^{2 2} L/h ² c, equation (35) - k* 12L/h ² c, frictional resistance parameter, equation (46) - L length of conduit, ft - m positive integer - n positive integer - p pressure, lb/ft² - p ₀ constant pressure at inlet of conduit, lb/ft² - P pressure plus elevation head, p+gz, equation (4) - mean value of P over the conduit width h - P ₀ p ₀+gz ₀, lbs/ft² - R frictional resistance coefficient for steady state wall shear, lb sec/ft⁴ - s positive integer; also, condensation, equation (6) - t time, sec - t ct/L, dimensionless time - u x component of fluid velocity, ft/sec - u _m mean velocity in conduit, equation (12), ft/sec - u ₀(y) velocity profile in Poiseuille flow, equation (19), ft/sec - transformed velocity - U initial mean velocity in conduit - x distance along conduit, measured from valve (fig. 1), ft - x x/L, dimensionless distance - y distance normal to conduit wall (fig. 1), ft - y y/h, equation (25) - z elevation, measured from arbitrary datum, ft - z ₀ elevation of constant pressure source, ft - isothermal bulk compression modulus, lbs/ft² - _n , equation (37) - _n (2n–1)/2, equation (36) - viscosity, slugs/ft sec - / = kinematic viscosity, ft²/sec - density of fluid, slugs/ft³ - ₀ density of undisturbed fluid, slugs/ft³ - ø angle between conduit and vertical (fig. 1) The research upon which this paper is based was supported by a grant from the National Science Foundation.     

Maximum density effects on thermal instability of horizontal laminar boundary layers

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

On the effects of a periodic, spanwise variation of suction upon asymptotic suction flow

Summary The effects of superposing streamwise vorticity, periodic in the lateral direction, upon two-dimensional asymptotic suction flow are analyzed. Such vorticity, generated by prescribing a spanwise variation in the suction velocity, is known to play an important role in unstable and turbulent boundary layers. The flow induced by the variation has been obtained for a freestream velocity which (i) is steady, (ii) oscillates periodically in time, (iii) changes impulsively from rest. For the oscillatory case it is shown that a frequency can exist which maximizes the induced, unsteady wall shear stress for a given spanwise period. For steady flow the heat transfer to, or from a wall at constant temperature has also been computed.Nomenclature (x, y, z) spatial coordinates - (u, v, w) corresponding components of velocity - (, , ) corresponding components of vorticity - t time - stream function for v and w - v _w mean wall suction velocity - nondimensional amplitude of variation in wall suction velocity - characteristic wavenumber for variation in direction of z - T temperature - P pressure - density - coefficient of kinematic viscosity - coefficient of thermal diffusivity - (/v _w)² - frequency of oscillation of freestream velocity - nondimensional amplitude of freestream oscillation - /v _w ² - z z - y –v _w y/ - v _w ² t/4 - /v _w - U ₀ characteristic freestream velocity - u/U ₀ - coefficient of viscosity - _w wall shear stress - Prandtl number (/) - q heat transfer to wall - T _w wall temperature - T (T _w–T)/(T _w–)     

The investigation on the properties and structures of starting vortex flow past a backward-facing step by WBIV technique

An experimental investigation of a starting vortex flow around a backward-facing step was conducted in a water channel. The properties and structures of the flow were investigated by qualitative flow visualization using the hydrogen bubble method and by quantitative velocity and vorticity measurements using White-light Bubble Image Velocimetry (WBIV) — a newly developed PIV method. Some invariant properties and 4-stage structures of starting vortex flow were observed.List of symbols a flow acceleration during starting stage - h height of backward-facing step - d _v dimensionless vortex size - t time - t dimensionless time - U free uniform velocity - u, v streamwise and spanwise velocity components respectively - Re Reynolds number based on a and h - x, y streamwise and spanwise coordinates respectively in flow field - x _c , y _c dimensionless vortex center position - vorticity - ov dimensionless vorticity - _max maximum vorticity - ov _max dimensionless maximum vorticity - circulation - dimensionless circulation - kinematic viscosity This work was supported by the CNSF Grant 1939 100-1-3