Transition,turbulence and oscillating flow in a pipe a visual study

Transitional and turbulent oscillatory flow in a rigid pipe with long entry sections was investigated using flow visualization to establish the existence of coherent structures. Flow tracer and high speed motion pictures were used. The simple harmonic motion of a scotch yoke and flywheel linked to a piston and cylinder provided the flow driving force. The camera was convected with the flow by attaching it through a gearing system to the scotch yoke.List of symbols A cross sectional area of flow - C, K constants - D pipe diameter - N _Re,ave Reynolds number based on average velocity (DU _ave /v) - N _Re,p Reynolds number based on maximum oscillatory velocity (DU _max /v) - Reynolds number based on maximum oscillatory velocity and Stokes (boundary) layer thickness (U _max /v) - R pipe radius - U instantaneous velocity in the flow direction - short-term average instantaneous velocity - U ^* friction velocity (U _ave (f/2)^1/2) - U _amp amplitude parameter (U _max /U _ave ) - U _ave average velocity - U _s steady velocity - U _t instantaneous oscillatory velocity - U _max maximum oscillatory velocity ( X _max /T) - u _r , u _z deviations from _r, and _z - y radial coordinate from wall (R–r) - y ⁺ dimensionless radial coordinate from wall (y U*/v) - frequency parameter [R (/v) ^1/2] - Stokes (boundary) layer thickness [C (2 v/)^1/2] - normalized time into cycle - fluid viscosity - v fluid kinematic viscosity (/) - density - angular frequency (2/T) - - overbar, average - sub-c critical value     

Finite difference analysis of plane Poiseuille and Couette flow developments

Summary The problem of flow development from an initially flat velocity profile in the plane Poiseuille and Couette flow geometry is investigated for a viscous fluid. The basic governing momentum and continuity equations are expressed in finite difference form and solved numerically on a high speed digital computer for a mesh network superimposed on the flow field. Results are obtained for the variations of velocity, pressure and resistance coefficient throughout the development region. A characteristic development length is defined and evaluated for both types of flow.Nomenclature h width of channel - L ratio of development length to channel width - p fluid pressure - p ₀ pressure at channel mouth - P dimensionless pressure, p/ ² - P ₀ dimensionless pressure at channel mouth - P pressure defect, P ₀–P - (P)₀ pressure defect neglecting inertia - Re Reynolds number, uh/ - u fluid velocity in x-direction - mean u velocity across channel - u ₀ wall velocity - U dimensionles u velocity u/ - U _c dimensionless centreline velocity - U ₀ dimensionless wall velocity - v fluid velocity in y-direction - V dimensionless v velocity, hv/ - x coordinate along channel - X dimensionless x-coordinate, x/h ² - y coordinate across channel - Y dimensionless y-coordinate, y/h - resistance coefficient, - ₀ resistance coefficient neglecting inertia - fluid density - fluid viscosity     

An experimental study on the two-dimensional opposed wall jet in a turbulent boundary layer: Change in the flow pattern with velocity ratio

Results of the measurement of flow properties in a two-dimensional turbulent wall jet which is injected into the turbulent boundary layer in the direction opposite to that of the boundary layer flow are presented by varying the ratio of the jet issuing velocity to the mainstream velocity . This flow includes the flow separation and the recirculating flow, and so it requires the magnitude and direction of instantaneous velocity be measured. In the present work, a tandem hot-wire probe is manufactured and its characteristics are examined experimentally. With the use of this probe the change in the penetration length of the jet with the velocity ratio is investigated. It is clarified that two regimes of flow patterns exist: in the low velocity ratio the penetration length of the jet changes little with , and in the high velocity ratio it goes far from the nozzle with increasing . Streamlines, turbulence intensity contours and static pressure contours are presented in the two typical velocity ratios corresponding to each of two flow patterns, and they are compared.List of symbols b ₀ nozzle width (Fig. 1) - , e mean and fluctuating output voltages of hot-wire anemometer - p, p mean static pressure, p = p – p _o - p ₀ static pressure in undisturbed mainstream - p _w , p _w mean wall pressure, p _w = p _w – p _o - U ₀ mainstream velocity - U _j jet velocity at the nozzle exit - , u mean and fluctuating velocity components in x-direction - u _e effective cooling velocity - x distance along the wall from nozzle exit - x _pmax, x _pmin positions where the wall pressure has maximum and minimum values respectively - x _s penetration length of jet - y distance from the wall - forward flow fraction - ₁, ₂ yaw and pitch angles of flow direction (Fig. 4) - velocity ratio, = U _j /U _o - density of the fluid - nondimensional stream function The authors wish to express their appreciation to Prof. Toshio Tanaka of Gifu University for his advice in the course of the experiment. They also would like to thank the Research Foundation for the Electrotechnology of Chubu which partly supported this work.     

Laminar flow in the entrance region of a porous-wall channel

Summary The effect of fluid injection at the walls of a two-dimensional channel on the development of flow in the entrance region of the channel has been investigated. The integral forms of the boundary layer equations for flow in the channel were set up for an injection velocity uniformly distributed along the channel walls.With an assumed polynomial of the n-th degree for the one-parameter velocity profile a solution of the above boundary layer equations was obtained by an iteration method. A closed form solution was also obtained for the case when a similar velocity profile was assumed. The agreement between the entrance region velocity profiles of the present analysis for an impermeable-walled channel and of Schlichting¹) and Bodoia and Osterle²) is found to be very good.The results of the analysis show that fluid injection at the channel walls increases the rate of the growth of the boundary layer thickness, and hence reduces considerably the entrance length required for a fully developed flow.Nomenclature h half channel thickness - L entrance length with wall-injection - L ₀ entrance length without wall-injection - p static pressure - p=p/U ₀ ² dimensionless pressure - Re=U ₀ h/ Reynolds number at inlet cross-section - u velocity in the x direction at any point in the channel - =u/U ₀ dimensionless velocity in the x direction at any point in the channel - U _av average velocity at a channel cross-section - U _c center line velocity - U ₀ inlet cross-section velocity - _c =U _c /U ₀ dimensionless center line velocity - v velocity in the y direction at any point in the channel - v ₀ constant injection velocity of fluid at the wall - v=v/v ₀ dimensionless velocity in the y direction at any point in the channel - x distance along the channel wall measured from the inlet cross-section - x=x/hRe dimensionless distance in the x direction - y distance perpendicular to the channel wall - y=y/h dimensionless distance in the y direction - thickness of the boundary layer - =/h dimensionless boundary layer thickness - =/ dimensionless distance within the boundary layer region - =v ₀ h/ injection parameter or injection Reynolds number - kinematic viscosity - 1+ie - mass density of the fluid - parameter defined in (14)     

Measurements of thermally stratified pipe flow using image-processing techniques

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

The effect of surface mass transfer on buoyancy induced flow in a variable porosity medium adjacent to a vertical heated plate

The effect of surface mass transfer on buoyancy induced flow in a variable porosity medium adjacent to a heated vertical plate is studied for high Rayleigh numbers. Similarity solutions are obtained within the frame work of boundary layer theory for a power law variation in surface temperature,T _Wx and surface injectionv _Wx(–1/2). The analysis incorporates the expression connecting porosity and permeability and also the expression connecting porosity and effective thermal diffusivity. The influence of thermal dispersion on the flow and heat transfer characteristics are also analysed in detail. The results of the present analysis document the fact that variable porosity enhances heat transfer rate and the magnitude of velocity near the wall. The governing equations are solved using an implicit finite difference scheme for both the Darcy flow model and Forchheimer flow model, the latter analysis being confined to an isothermal surface and an impermeable vertical plate. The influence of the intertial terms in the Forchheimer model is to decrease the heat transfer and flow rates and the influence of thermal dispersion is to increase the heat transfer rate.

---

Der Effekt des Oberflächenstoffaustausches bei auftriebsinduzierter Strömung in einem variablen porösen Medium, das an eine vertikale, beheizte Platte angrenzt

Zusammenfassung Es wird der Effekt des Oberflächenstoffaustausches in auftriebsinduzierter Strömung in einem variablen porösen Medium, das an eine vertikale, beheizte Platte angrenzt, für große Reynoldszahlen untersucht. Ähnliche Lösungen werden im Rahmen der Grenzschicht-Theorie, durch Variation des Potenzansatzes der Oberflächentemperatur,T _Wx , und der Oberflächengeschwindigkeit,v _Wx(–1/2), erreicht. Die Analyse vereinigt sowohl den Ausdruck, der Porösität und Permeabilität verbindet, als auch den Ausdruck, der Porösität und Wärmeleitfähigkeit miteinander verbindet. Der Einfluß der Temperaturverteilung auf Strömung und Wärmeübergangskennzahlen wird ebenfalls im Detail analysiert. Als Ergebnis der vorliegenden Untersuchung ergibt sich die Tatsache, daß variable Porösität Wärmeübertragungsrate und Betrag der Geschwindigkeit in Wandnähe steigert. Die bestimmenden Gleichungen, sowohl für das Darcysche Strömungsmodell als auch für das Forchheimersche Strömungsmodell, werden mit Hilfe eines implizierten Differenzenschemas gelöst. Die Berechnung wird für die beiden Fälle, isotherme Oberfläche und undurchlässige vertikale Platte, angewandt. Der Einfluß der Terme für die Trägheitskräfte im Forchheimerschen Modell senkt Wärmeübergangs- und Durchgangsrate, wogegen die Wärmeübergangsrate durch den Einfluß der Temperaturverteilung erhöht wird.

---

Nomenclature a constant defined by Eq. (12) - A constant defined by Eq. (12) - B constant defined by Eq. (3) - b _s/_f ratio of thermal conductivities - C constant defined by Eq. (1) - C _P specific heat of the convective fluid - d particle diameter - f dimensionless function defined by Eq. (14) - f _w lateral mass flux parameter - g acceleration due to gravity - k ₀ mean permeability of the mediumk ₀= ₀ ³ d ²/150 (1– ₀)² k ₀=1.75d/(1– ₀) 150 (Inertia parameter) - L length of the source or sink - m mass transfer - n constant defined in Eq. (12) - k (y) permeability of the porous medium - k (y) interial coefficient in the Ergun expression - Gr modified Grashof numberGr=(g k ₀ k ₀ (T _w–))/ ² - R _a Rayleigh number (g k ₀ x T _w–)/ - R _ad modified Rayleigh number (g k ₀ d|T _w–|)/ - N _u Nusselt number - s x/d - Q overall heat transfer rate - T temperature - T _w surface temperature - T ambient fluid temperature - u velocity in vertical direction - v velocity in horizontal direction - x vertical coordinate - y horizontal coordinate Greek symbols ₀ mean thermal diffusivity _f/ C_p - coefficient of thermal expansion - constant defined in Eq. (4) - ratio of particle to bed diameter - _e effective thermal conductivity - f thermal conductivity of fluid - _s thermal conductivity of solid - dimensionless similarity variable in Eq. (13) - value of at the edge of the boundary layer - constant defined in Eq. (1) - _e effective molecular thermal diffusivity - (y) porosity of the medium - ₀ mean porosity of the medium - viscosity of the fluid - ₀ density of the convective fluid - stream function - w condition at the wall - condition at infinity     

LDA measurements in low Reynolds number turbulent boundary layer

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

Prediction of stagnation point heat transfer for a single round jet impinging on a concave hemispherical surface

This paper deals with a systematic procedure for assessment of fluid flow and heat transfer parameters for a single round jet impinging on a concave hemispherical surface. Based on Scholkemeier's modifications of the Karman-Pohlhausen integral method, expressions are derived for evaluation of the momentum thickness, boundary layer thickness and the displacement thickness at the stagnation point. This is followed by the estimation of thermal boundary layer thickness and local heat transfer coefficients. A correlation is presented for the Nusselt number at the stagnation point as a function of the Reynolds number for different non-dimensional distances from the exit plane of the jet to the impingement surface.

---

Bestimmung des Staupunktes bei der Wärmeübertragung für einen einzelnen Strahl, der auf eine konkave halbkugelige Oberfläche trifft

Zusammenfassung Diese Arbeit beschäftigt sich mit dem systematischen Verfahren der Bewertung von Fluidströmungen und Wärmeübertragungsparametern für einen einzelnen runden Strahl, der auf eine konkave halbkugelförmige Oberfläche trifft. Das Verfahren beruht auf Scholkemeiers Modifikation des Karman-Pohlhausen Integrationsverfahrens. Ausdrücke sind für die Berechnung der Impuls-Dicke, der Grenzschichtdicke und der Verschiebungsdicke am Staupunkt hergeleitet worden. Dies ist aus der Berechnung der thermischen Grenzschichtdicke und des lokalen Wärmeübertragungskoeffizienten abgeleitet worden. Es wird eine Gleichung für die Nusselt-Zahl am Staupunkt als Funktion der Reynolds-Zahl für verschiedene dimensionslose Abstände vom Strahlaustrittspunkt bis zum Auftreffpunkt auf die Oberfläche vorgestellt.

---

Nomenclature c _p specific heat at constant pressure - d diameter of single round nozzle - h ₀ heat transfer coefficient at the stagnation point - H distance from the exit plane of the jet to the impingement surface - k thermal conductivity - Nu _0.5 Nusselt number based on impinging jet quantities=h _0.50/k - Nu _{0.5, 0} stagnation point Nusselt number=h _{0 0,50}/k - p pressure - p _a ambient pressure - p ₀ maximum pressure or stagnation pressure - p(x) static pressure at a distancex from the stagnation point - R radius of curvature of the hemisphere - Re _J jet Reynolds number=U _Jd/ - Re _0.5 Reynolds number based on impinging jet quantities=u _{m0 0.50}/ - T temperature - T _a room temperature - T _J jet temperature - T _W wall temperature - u velocity component inx andx directions (Fig. 1) - u _m jet centerline (or maximum) free jet velocity: external (or maximum) boundary layer velocity aty= _m - u _m0 arrival velocity defined as the maximum velocity the free jet would have at the plane of impingement if the plane were not there - U _J jet exit velocity - x* non-dimensional coordinate starting at the stagnation point=x/2 _0.50 - x, y rectangular Cartesian coordinates - y coordinate normal to the wall starting at the wall - ratio of thermal to velocity boundary layer thickness= _T/_m - ₀ ratio of thermal to velocity boundary layer thickness at the stagnation point - * inner layer displacement thickness - _0.50 jet half width at the plane of impingement if the plate were not there - _m inner boundary layer thickness atu=u _m - Pohlhausen's form parameter - dynamic viscosity - kinematic viscosity=/ - fluid density - momentum thickness - ₀ momentum thickness at the stagnation point     

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     

On the effects of uniform high suction on the rotationally symmetric flow of a conducting liquid near a stationary disc in the presence of a transverse magnetic field

In the present paper an attempt has been made to find out effects of uniform high suction in the presence of a transverse magnetic field, on the motion near a stationary plate when the fluid at a large distance above it rotates with a constant angular velocity. Series solutions for velocity components, displacement thickness and momentum thickness are obtained in the descending powers of the suction parameter a. The solutions obtained are valid for small values of the non-dimensional magnetic parameter m (= 4 _e ² H ₀ ² /) and large values of a (a2).Nomenclature a suction parameter - E electric field - E _r , E , E _z radial, azimuthal and axial components of electric field - F, G, H reduced radial, azimuthal and axial velocity components - H magnetic field - H _r , H , H _z radial, azimuthal and axial components of magnetic field - H ₀ uniform magnetic field - H* displacement thickness and momentum thickness ratio, */ - h induced magnetic field - h _r , h , h _z radial, azimuthal and axial components of induced magnetic field - J current density - m nondimensional magnetic parameter - p pressure - P reduced pressure - R Reynolds number - U ₀ representative velocity - V velocity - V _r , V , V _z radial, azimuthal and axial velocity components - w ₀ uniform suction through the disc. - density - electrical conductivity - kinematic viscosity - _e magnetic permeability - a parameter, (/)^1/2 z - a parameter, a - * displacement thickness - momentum thickness - angular velocity     

Dielectric properties of the composite system polyurethane — sodium chloride

Summary The dielectric properties of the composite system polyurethane-sodium chloride have been measured at frequencies between 10^–4 Hz and 3 · 10⁵ Hz in the temperature range from –150 °C up to +90 dgC. Three dielectric loss mechanisms have been found; they are indicated by ₁, ₂ and. The activation energy of the ₁-transition is 35 kcal/mole, that of the-transition 6.7 kcal/mole. The ₂-loss peak was only observed at frequencies of 10³ Hz and higher, forming one broad peak with the ₁-loss peak at lower frequencies. In the composite materials, the- and ₂-loss peaks measured at fixed frequencies were found at the same temperature. The ₂-loss peak, however, was shifted to a lower temperature, due to the sodium chloride filler. Comparison of experimental data of and tan with theoretical predictions concerning the dielectric properties of composite systems showed only partial agreement. The difference mainly consisted in. the temperature shift in the tan-peak of the ₁-transition.

---

Zusammenfassung Die dielektrischen Eigenschaften des Verbundssystems Kochsalz-Polyurethankautschuk wurden im Frequenzgebiet zwischen 10^–4 Hz und 3.10⁵ Hz und im Temperaturbereich von –150 °C bis +90 °C gemessen. Es wurden drei dielektrische Verlustmaxima gefunden, die mit ₁, ₂ und angedeutet werden. Die Aktivierungsenergie des ₁-Überganges beträgt 35 kcal/Mol, die des-Überganges 6.7 kcal/Mol. Das ₂-Maximum konnte nur bei Frequenzen höher als 10³Hz vom ₁-Maximum gesondert erfaßt werden. Die Lage der ₂- und-Maxima war vom Füllgrad unabhängig. Das ₁-Maximum verschiebt sich mit steigendem Füllgrad zu niedrigeren Temperaturen. Die gemessenen Werte von und tan stimmen nur teilweise mit den Aussagen einer Theorie der dielektrischen Eigenschaften von Mischkörpern überein.

     

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     

Asymptotic Nusselt numbers for Newtonian flow through a parallel-plate channel with recycle

General expressions for evaluating the asymptotic Nusselt number for a Newtonian flow through a parallel-plate channel with recycle at the ends have been derived. Numerical results with the ratio of thicknesses as a parameter for various recycle ratios are obtained. A regression analysis shows that the results can be expressed by log Nur^0.83=0.3589 (log)² -0.2925 (log) + 0.3348 forR 3, 0.1 0.9; logNu=0.5982(log)² +0.3755 × 10^–2 (log) +0.8342 forR 10^–2, 0.1 0.9.

---

Asymptotische Nusselt-Zahlen für die Newtonsche Strömung durch einen Kanal aus parallelen Platten mit Rückführung

Zusammenfassung In dieser Untersuchung wurden allgemeine Ausdrücke hergeleitet um die asymptotische Nusselt-Zahl für eine Newtonsche Strömung durch einen Kanal aus parallelen Platten mit Rückführung an den Enden berechnen zu können. Es wurden numerische Ergebnisse mit den Dicken-Verhältnissen, als Parameter für verschiedene Rückführungs-verhältnisse, erhalten. Eine Regressionsanalyse zeigt, daß die Ergebnisse wie folgt ausgedrückt werden können: log Nur^0,83=0,3589 (log)² -0,2925 (log) + 0,3348 fürR 3, 0,1 0,9; logNu=0,5982(log)² +0,3755 × 10^–2 (log) + 0,8342 fürR 10^–2, 0,1 0,9.

---

Nomenclature A₁ shooting value,d(0)/d - A₂ shooting value,d(1)/d - B channel width - Gz Graetz number, U_bW²/L - h _m logarithmic average convective heat transfer coefficient - h _x average local convective heat transfer coefficient - k thermal conductivity - L channel length - Nu average local Nusselt number, 2 h_xW/k - Nu _m logarithmic average Nusselt number, 2h_mW/k - R recycle ratio, reverse volume flow rate divided by input volume flow rate - T temperature of fluid - T _m bulk temperature, Eq. (8) - T ₀ temperature of feed stream - T _s wall temperature - U velocity distribution - U _b reference velocity,V/BW - V input volume flow rate - v dimensionless velocity distribution, U/U_b - W channel thickness - x longitudinal coordinate - y transversal coordinate - Z₁-z₆ functions defined in Eq. (A1) - thermal diffusivity - least squares error, Eq. (A7) - weight, Eqs. (A8), (A9) - dimensionless coordinate,y/W - dimensionless coordinate,x/GzL - function, Eq. (7)     

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     

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     

Some flow characteristics of lobed forced mixers at higher velocity ratios

The flow characteristics of two types of lobed forced mixers, the unscalloped and the scalloped mixers, have been examined at velocity ratios higher than unity, in relation to the variation of mass flux uniformity, the decay of the streamwise vorticity, the variation of turbulent kinetic energy and the growth of the shear layer with distance from the trailing edge. Three trailing edge configurations have also been considered for each type of mixer, namely a square wave, a semi-circular wave and a triangular wave. The analysis showed that the strength of the streamwise vorticity shed at the trailing edge and the subsequent decaying rate with downstream distance are found to be very important in studying the mixing effectiveness of the lobed mixers.List of Symbols C _I normalized streamwise circulation, _s/U _r h tan - _s streamwise circulation - k turbulent kinetic energy = 1/2(u²+v²+w²) - Re Reynolds number, U _r /=2.27×10⁴ - h (=) Lobe height, 33 mm - U ₁, U ₂ mean velocity of the slow and fast streams - U _r reference mean velocity, (U ₁ + U ₂)/2 = 10 m/s - U, u streamwise mean and the corresponding rms velocities - V, v horizontal mean and the corresponding rms velocities - W, w vertical mean and the corresponding rms velocities - x,y, z streamwise, horizontal and vertical directions - A _wake cross-sectional area bounded by the wake region. The wake region boundary is defined at the region bounded by one half of a lobe along the y/ direction and at the locations along the z/ direction where U ₂/U _r2 and U ₁/U _r1<0.95. - nominal lobe wavelength, 33 mm - half of the included divergent angle of the penetration region, 22° - U uniformity factor - momentum thickness Financial supports from the Applied Research Grant is gratefully acknowledged. The contribution of Mr. J. K. L. Teh, Dr. J. H. Yeo and Mr. T. H. Yip to the work presented here are sincerely appreciated.     

Theoretical investigation of the effect of an internal heat source on turbulent liquid metal heat transfer in a channel with uniform heat flux

Summary The effect of an internal heat source on the heat transfer characteristics for turbulent liquid metal flow between parallel plates is studied analytically. The analysis is carried out for the conditions of uniform internal heat generation, uniform wall heat flux, and fully established temperature and velocity profiles. Consideration is given both to the uniform or slug flow approximation and the power law approximation for the turbulent velocity profile. Allowance is made for turbulent eddying within the liquid metal through the use of an idealized eddy diffusivity function. It is found that the Nusselt number is unaffected by the heat source strength when the velocity profile is assumed to be uniform over the channel cross section. In the case of a 1/7-power velocity expression, the Nusselt numbers are lower than those in the absence of internal heat generation, and decrease with diminishing eddy conduction. Nusselt numbers, in the absence of an internal heat source, are compared with existing calculations, and indications are that the present results are adequate for preliminary design purposes.Nomenclature A hydrodynamic parameter - a half height of channel - a ₁ a constant, 1+0.01 Pr Re ^0.9 - a ₂ a constant, 0.01 Pr Re ^0.9 - C _p specific heat at constant pressure - D _h hydraulic diameter of channel, 4a - h heat transfer coefficient, q _w/(t _w–t _b) - I ₁ integral defined by (17) - I ₂ integral defined by (18) - k diffusivity parameter, (1+0.01 Pr Re ^0.9)^1/2 - m exponent in power velocity expression - Nu Nusselt number, hD _h/ - Nu ₀ Nusselt number in absence of internal heat generation - Pr Prandtl number, / - Q heat generation rate per volume - q _w wall heat flux - Re Reynolds number for channel, 2/ - s ratio of heat generation rate to wall heat flux, Qa/q _w - T dimensionless temperature, (t _w–t)/(t _w–t _b) - t fluid temperature, t _w wall temperature, t _b fluid bulk temperature - u fluid velocity in x direction, , fluid mean velocity - x longitudinal coordinate measured from channel entrance - x ⁺ dimensionless longitudinal coordinate, 2(x/a)/Pr Re - y transverse coordinate measured from channel centerline - z transverse coordinate measured from channel wall, a–y - molecular diffusivity of heat, /C _p - dummy variable of integration - dummy variable of integration - _H eddy diffusivity of heat - _M eddy diffusivity of momentum - dummy variable of integration - fluid thermal conductivity - _T dimensionless diffusivity, Pr ( _H/) - fluid kinematic viscosity - dummy variable of integration - fluid density - dummy variable of integration - ratio of eddy diffusivity for heat transfer to that for momentum transfer, _H/ _M - average value of - dimensionless velocity distribution, u/     

On the measurement of normal stress differences in steady shear flow

Summary Compared to the similar pressure-distribution cone-and-plate apparatus of Adams and Lodge (4), the new apparatus' improvements include: temperature control of the cone (as well as the plate); increased stiffening of the frame; four (not three) pressuremeasuring holes in the cone/plate region; inclusion of a pressure-measuring hole on the axis under the cone truncation; exclusive use of a vertical free liquid boundary at the cone rim (without a sea of liquid). Temperature control of the rotating cone and of the fixed plate leads to acceptable temperature uniformity in the test liquid for working temperatures within 10°C or 20°C of ambient; the corresponding interval is about 1°C if the cone temperature control is abandoned. Pressure gradients measured using a Newtonian liquid agree with those calculated using Walters' eq. (3). For a viscoelastic liquid, after subtracting inertial contributions, pressure distributions measured at a given shear rate in the cone/plate region do not change when the gap angle is changed from 2° to 3°, showing that the effects of secondary flow are negligible. Values ofN ₃ =N ₁ + 2N ₂ obtained from the gradients of these distributions are believed to be in error by not more than ±1 Pa, in favorable cases. The present most useful ranges are: 10 to 5000 Pa forN ₃; 0.1 to 200 sec^–1 for shear rate; up to 5 Pa s for viscosity; and 5 to 40°C for temperature. As an application, it is shown that adding 0.1% of a high molecular weight polyisobutylene to a 2% polyisobutylene solution doublesN ₃ and has no detectable effect on the viscosity measured at low shear rates with a Ferranti-Shirley viscometer.

---

Zusammenfassung Im Vergleich zu dem ähnlichen Kegel-Platte-Gerät von Adams und Lodge (4) zur Messung der Druckverteilung wurden an dem neuen Gerät die folgenden Verbesserungen vorgenommen: Temperaturregelung an Kegel und Platte, Versteifung des Rahmens, vier (anstatt drei) Druckmeßlöcher im Kegel-Platte-Bereich, ein zusätzliches Druckmeßloch auf der Achse unter der Kegelstumpf-Deckfläche, ausschließliche Verwendung einer vertikalen freien Grenzfläche der Flüssigkeit am Kegelrand (ohne umgebenden Flüssigkeitssee). Die Temperaturregelung des rotierenden Kegels und der festen Platte führt zu einer ausreichenden Temperaturgleichförmigkeit in der Testflüssigkeit für Betriebstemperaturen, die höchstens um 10–20°C von der Umgebungstemperatur abweichen. Dieses Intervall beträgt dagegen nur etwa 1°C, wenn auf die Temperaturregelung am Kegel verzichtet wird. Für newtonsche Flüssigkeiten entsprechen die gemessenen Druckgradienten den mittels der Gleichung von Walters (3) berechneten. Für viskoelastische Flüssigkeiten zeigen sich bei der Änderung des Spaltwinkels von 2° auf 3° nach Abzug der Trägheitsbeiträge keine Änderungen der bei einer bestimmten Schergeschwindigkeit gemessenen Druckverteilung. Dies zeigt, daß Sekundärströmungseffekte vernachlässigbar sind. Es darf angenommen werden, daß die Werte vonN ₃ =N ₁ + 2N ₂, die man aus den Gradienten dieser Verteilungen erhält, unter günstigen Umständen mit einem Fehler von nicht mehr als ±1 Pa behaftet sind. Gegenwärtig liegen die günstigsten Bereiche bei 10 bis 5000 Pa fürN ₃, 0,1 bis 200 s^–1 für die Schergeschwindigkeit, unterhalb von 5 Pa s für die Viskosität und 5 bis 40°C für die Temperatur. Als Anwendung wird gezeigt, daß ein Zusatz von 0,1% hochmolekularen Polyisobutylens zu einer 2%igen Polyisobutylenlösung den Wert vonN ₃ verdoppelt, aber keinen erkennbaren Einfluß auf die (bei geringen Schergeschwindigkeiten mit einem Ferranti-Shirley-Viskosimeter gemessen) Viskosität hat.

---

udsf unidirectional shear flow - TCP truncated-cone and plate - N ₁,N ₂ 1st and 2nd normal stress differences in udsf - N ₃ N ₁ + 2N ₂ - : = A is defined by the equationA := B - P ^* hole pressurePw – Pm; Pw, Pm = pressures measured by flush transducer and by hole-mounted transducer - t time - , strain rate, shear rate - (P,t) covariant body metric tensor at particleP and timet - ⁱ , _i covariant and contravariant udsf body base vectors (i = 1, 2, 3) - ^–1 inverse of - R, plate radius, cone/plate gap angle - r ₀,h ₀ radius and height of cone truncation - r,, spherical polar coordinates; cone axis = 0; plate surface = /2 - physical components of stress; for a tensile component - cone angular velocity - p on the plate = /2 - ,T, density, absolute temperature, viscosity - P 0.15 ²(r ² –R ²) (inertial contribution) [2.7] - P _ve contribution [2.8] from flow perturbations of viscoelastic origin - r _i i = 1,2,3,4; values ofr at centers of holes in cone/plate region - P _i () pressure change recorded by transducerTi when cone angular velocity goes from zero to - 1/2 {P _i ()+ P _i (–)} (average for 2 senses of rotation) - rim pressure, from least-squares line through four points - Re Reynolds' number:R ²/ - (P,t)/t With 11 figures and 2 tables     

Investigations of mean and turbulent flow characteristics of a two dimensional plane diffuser

This paper reports the investigation of mean and turbulent flow characteristics of a two-dimensional plane diffuser. Both experimental and theoretical details are considered. The experimental investigation consists of the measurement of mean velocity profiles, wall static pressure and turbulence stresses. Theoretical study involves the prediction of downstream velocity profiles and the distribution of turbulence kinetic energy using a well tested finite difference procedure. Two models, viz., Prandtl's mixing length hypothesis and k- model of turbulence, have been used and compared. The nondimensional static pressure distribution, the longitudinal pressure gradient, the pressure recovery coefficient, percentage recovery of static pressure, the variation of U _max/U _bar along the length of the diffuser and the blockage factor have been valuated from the predicted results and compared with the experimental data. Further, the predicted and the measured value of kinetic energy of turbulence have also been compared. It is seen that for the prediction of mean flow characteristics and to evaluate the performance of the diffuser, a simple turbulence model like Prandtl's mixing length hypothesis is quite adequate.List of symbols C ₁ , C ₂ ,C turbulence model constants - F _x body force - k kinetic energy of turbulence - l _m mixing length - L length of the diffuser - u, v, w rms value of the fluctuating velocity - u, v, w turbulent component of the velocity - mean velocity in the x direction - _A average velocity at inlet - U _bar average velocity in any cross section - U _max maximum velocity in any cross section - V mean velocity in the y direction - W local width of the diffuser at any cross section - x, y coordinates - dissipation rate of turbulence - _m eddy diffusivity - Von Karman constant - mixing length constant - _l laminar viscosity - _eff effective viscosity - v kinematic viscosity - density - _k effective Schmidt number for k - effective Schmidt number for - stream function - non dimensional stream function     

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