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)     

Heat conduction in a cylinder crossed by an electric current with skin effect

The steady periodic temperature distribution in an infinitely long solid cylinder crossed by an alternating current is evaluated. First, the time dependent and non-uniform power generated per unit volume by Joule effect within the cylinder is determined. Then, the dimensionless temperature distribution is obtained by analytical methods in steady periodic regime. Dimensionless tables which yield the amplitude and the phase of temperature oscillations both on the axis and on the surface of copper or nichrome cylindrical electric resistors are presented.

---

Wärmeleitung in einem stromdurchflossenen Zylinder unter Berücksichtigung des Skin-Effektes

Zusammenfassung Es wird die periodische Temperaturverteilung für den eingeschwungenen Zustand in einem unendlich langen, von Wechselstrom durchflossenen Vollzylinder ermittelt. Zuerst erfolgt die Bestimmung der zeitabhängigen, nichgleichförmigen Energiefreisetzung pro Volumeneinheit des Zylinders infolge Joulescher Wärmeentwicklung und anschließend die Ermittlung der quasistationären Temperaturverteilung auf analytischem Wege. Amplitude und Phasenverzögerung der Temperaturschwingungen werden für die Achse und die Oberfläche eines Kupfer- oder Nickelchromzylinders tabellarisch in dimensionsloser Form mitgeteilt.

---

Nomenclature A integration constant introduced in Eq. (2) - ber, bei Thomson functions of order zero - Bi Biot numberhr ₀/ - c speed of light in empty space - c ₁,c ₂ integration constants introduced in Eq. (46) - c _p specific heat at constant pressure - E electric field - E _z component ofE alongz - E time independent part ofE, defined in Eq. (1) - f function ofs and defined in Eq. (11) - g function ofs and defined in Eq. (37) - h convection heat transfer coefficient - H magnetic field - i imaginary uniti=(–1)^1/2 - I electric current - I _eff effective electric currentI _eff=I/2^1/2 - Im imaginary part of a complex number - J _n Bessel function of first kind and ordern - J electric current density - q _g power generated per unit volume - time average of the power generated per unit volume - time averaged power per unit length - r radial coordinate - R electric resistance per unit length - r ₀ radius of the cylinder - Re real part of a complex number - s dimensionless radial coordinates=r/r ₀ - s, s integration variables - t time - T temperature - time averaged temperature - T _f fluid temperature outside the boundary layer - time average of the surface temperature of the cylinder - u, functions ofs, and defined in Eqs. (47) and (48) - W Wronskian - x position vector - x real variable - Y _n Bessel function of second kind and ordern - z unit vector parallel to the axis of the cylinder - z axial coordinate - · modulus of a complex number - equal by definition Greek symbols amplitude of the dimensionless temperature oscillations - electric permittivity - dimensionless temperature defined in Eq. (16) - ₀, ₁, ₂ functions ofs defined in Eq. (22) - thermal conductivity - dimensionless parameter=(2)^1/2 - magnetic permeability - ₀ magnetic permeability of free space - function of defined in Eq. (59) - dimensionless parameter=c _p/() - mass density - electric conductivity - dimensionless time=t - phase of the dimensionless temperature oscillations - function ofs:= ₁+i ₂ - angular frequency - dimensionless parameter=()^1/2 r ₀     

On natural convection from a vertical plate with a prescribed surface heat flux in porous media

This paper presents a theoretical and numerical investigation of the natural convection boundary-layer along a vertical surface, which is embedded in a porous medium, when the surface heat flux varies as (1 +x ²)), where is a constant andx is the distance along the surface. It is shown that for > -1/2 the solution develops from a similarity solution which is valid for small values ofx to one which is valid for large values ofx. However, when -1/2 no similarity solutions exist for large values ofx and it is found that there are two cases to consider, namely < -1/2 and = -1/2. The wall temperature and the velocity at large distances along the plate are determined for a range of values of .Notation g Gravitational acceleration - k Thermal conductivity of the saturated porous medium - K Permeability of the porous medium - l Typical streamwise length - q _w Uniform heat flux on the wall - Ra Rayleigh number, =gK(q _w /k)l/(v) - T Temperature - Too Temperature far from the plate - u, v Components of seepage velocity in the x and y directions - x, y Cartesian coordinates - Thermal diffusivity of the fluid saturated porous medium - The coefficient of thermal expansion - An undetermined constant - Porosity of the porous medium - Similarity variable, =y(1+x ) ^/3/x ^1/3 - A preassigned constant - Kinematic viscosity - Nondimensional temperature, =(T – T )Ra^1/3 k/qw - Similarity variable, = =y(log_e x)^1/3/x ^2/3 - Similarity variable, =y/x ^2/3 - Stream function     

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 slot jet impinging on a concave semicylindrical surface

A systematic procedure has been laid out for assessment of fluid flow and heat transfer parameters for a slot jet impinging on a concave semicylindrical surface. Based on Walz's modifications of the Karman-Pohlhausen integral method, expressions have been derived for evaluation of the momentum thickness, boundary layer thickness and the displacement thickness at the stagnation point. The work then has been extended for the estimation of thermal boundary layer thickness and local heat transfer coefficients. A correlation has been 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.

---

Berechnung des Wärmeübergangs im Staupunkt eines Strahles, der aus einer rechteckigen öffnung auf eine konkave halbzylindrische Fläche auftrifft

Zusammenfassung Es wurde eine systematische Prozedur für die Abschätzung von Strömungs- und Wärmeübergangsparametern für einen Strahl, der auf eine konkave halbzylindrische Fläche auftrifft, aufgestellt. Basierend auf Walz's Modifikationen der Karman-Pohlhausen Integral-Methode, wurden Ausdrücke für die Berechnung der Impulsdicke, der Grenzschichtdicke und die Versetzungsdicke am Staupunkt abgeleitet. Die Arbeit wurde dann auf die Abschätzung der thermischen Grenzschichtdicke und der lokalen Wärmeübertragungskoeffizienten ausgedehnt. Es wird eine Beziehung für die Nusselt-Zahl am Staupunkt als eine Funktion der Reynolds-Zahl für verschiedene dimensionslose Abstände von der Austrittsfläche des Schlitzes bis zur Aufprallfläche aufgestellt.

---

Nomenclature c _p specific heat at constant pressure - 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 _.5 Nusselt number based on impinging jet quantities =h _0.50/k - Nu _.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 - p(x*) static pressure at nondimensional distancex* from the stagnation point - Re _J jet Reynolds number =U _J W/ - Re _0.5 Reynolds number based on impinging jet quantities =u _{m0 0.50}/ - T temperature - T* nondimensional temperature =(T–T _W)/(T _J–T _W) - T _a room temperature - T _J jet temperature - T _W wall temperature - u velocity component inx andx directions - 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 - W jet nozzle width - x* nondimensional coordinate starting at the stagnation point =x/2 _0.50 - x, y rectangular cartesian coordinates - y coordinate normal to the wall and 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 - _.50 jet half width at the plane of impingement if the plate were not there - d_.5 free jet (half width) thickness whereu=u _m/2 - _m inner boundary layer thickness atu =u _m - _T thermal boundary layer thickness - nondimensional coordinate normal to wall =y/ _m - _T nondimensional coordinate normal to wall =y/ _T - Pohlhausen's form parameter - dynamic viscosity - kinematic viscosity = / - fluid density - momentum thickness - ₀ momentum thickness at the stagnation point     

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      

Transport in ordered and disordered porous media V: Geometrical results for two-dimensional systems

In this paper we continue the geometrical studies of computer generated two-phase systems that were presented in Part IV. In order to reduce the computational time associated with the previous three-dimensional studies, the calculations presented in this work are restricted to two dimensions. This allows us to explore more thoroughly the influence of the size of the averaging volume and to learn something about the use of anon-representative region in the determination of averaged quantities.

Nomenclature

Roman Letters A interfacial area of the– interface associated with the local closure problem, m² - a _i i=1, 2, gaussian probability distribution used to locate the position of particles - l unit tensor - characteristic length for the-phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1,2,3 lattice vectors, m - m convolution product weighting function - m _V special convolution product weighting function associated with a unit cell - n _i i=1, 2 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume,V, m³ - x position of the centroid of an averaging area, m - x ₀ reference position of the centroid of an averaging area, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters V /V, volume average porosity - _{a i} standard deviation ofa _i - _r standard deviation ofr - intrinsic phase average of      

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.     

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     

Velocity measurements in a plane turbulent air jet at moderate Reynolds numbers

An experimental investigation of the moderate Reynolds number plane air jets was undertaken and the effect of the jet Reynolds number on the turbulent flow structure was determined. The Reynolds number, which was defined by the jet exit conditions, was varied between 1000 and 7000. Other initial conditions, such as the initial turbulence intensity, were kept constant throughout the experiments. Both hot-wire and laser Doppler anemometry were used for the velocity measurements. In the moderate Reynolds number regime, the turbulent flow structure is in transition. The average size and the number of the large scale of turbulence (per unit length of jet) was unaffected by the Reynolds number. A broadening of the turbulent spectra with increasing Reynolds number was observed. This indicated that there is a decrease in the strength of the large eddies resulting from a reduction of the relative energy available to them. This diminished the jet mixing with the ambient as the Reynolds number increased. Higher Reynolds numbers led to lower jet dilution and spread rates. On the other hand, at higher Reynolds numbers the dependence of jet mixing on Reynolds number became less significant as the turbulent flow structure developed into a self-preserving state.List of symbols b _u velocity half-width of the jet - C _u, C _u,0 constants defining the velocity decay rate - D nozzle width - E _u one dimensional power spectrum of velocity fluctuations - f frequency - K _u, K _u,0 constants defining the jet spread rate - k wavenumber (2f/U) - L longitudinal integral scale - R ₁₁ correlation function - r separation distance - Re jet Reynolds number (U ₀ D/v) - St Strouhal number (fD/U ₀) - t time - U axial component of the mean velocity - U _m mean velocity on the jet axis - U ₀ mean velocity at the jet exit - u the rms of u - u fluctuating component of the axial velocity - V lateral component of the mean velocity - fluctuating component of the lateral velocity - x axial distance from the nozzle exit - y lateral distance from the jet axis - z spanwise distance from the jet axis - v kinematic viscosity - time lag A version of this paper was presented as paper no. 86-0038 at the AIAA 24th Aerospace Sciences Meeting, Reno NV, USA, January 1986     

Experimental investigation of impinging jet arrays

We report on measurements of the velocity field and turbulence fluctuations in a hexagonal array of circular jets, impinging normally on a plane wall, using particle image velocimetry (PIV). Results for mean velocity and turbulent stresses are presented in various horizontal and vertical planes. From the measurements, we have identified some major features of impinging jet arrays and we discuss their mutual interaction, collision on the plate, and consequent backwash, which generate recirculating motion between the jets. The length of the jet core, the production of turbulence kinetic energy, and the model of the exhaust mechanisms for spent fluid are also discussed. The measurements indicated that the interaction between the self-induced cross flow and the wall jets resulted in the formation of a system of horseshoe-type vortices that circumscribe the outer jets of the array. The instantaneous snapshots of the velocity field reveal some interesting features of the flow dynamics, indicating a breakdown of some of the jets before reaching the plate, which may have consequences on the distribution of the instantaneous heat transfer.List of symbols D_m Nozzle diameter in multiple jet array nozzle plate (m) - D_s Pipe diameter in single jet rig (m) - H Distance between nozzle and impingement plate (m) - k Turbulent kinetic energy (m²/s²) - L Pipe length (m) - P_k Production of turbulent kinetic energy (m²/s³) - P_uu , P_vv Normal components of P_k (m²/s³) - P_uv Shear component of P_k (m²/s³) - s Pitch (m) - U_bulk Surface-averaged exit velocity (single jet) (m/s) - U_CL Center line jet exit velocity (jet array), m/s - u, v Mean velocity components in x and y directions (m/s) - u, v, w Instantaneous velocity in x, y, and z directions (m/s) - u, v, w Velocity fluctuation in x, y, and z directions (m/s) - u², v², w² Reynolds normal stress components (m²/s²) - uv Reynolds shear stress component (m²/s²) - x, z Coordinates parallel to impingement plate (m) - y Coordinate perpendicular to impingement plate (m)     

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      

Effect of measuring volume length on two-component laser velocimeter measurements in a turbulent boundary layer

The effects of finite measuring volume length on laser velocimetry measurements of turbulent boundary layers were studied. Four different effective measuring volume lengths, ranging in spanwise extent from 7 to 44 viscous units, were used in a low Reynolds number (Re=1440) turbulent boundary layer with high data density. Reynolds shear stress profiles in the near-wall region show that u v strongly depends on the measuring volume length; at a given y-position, u v decreases with increasing measuring volume length. This dependence was attributed to simultaneous validations on the U and V channels of Doppler bursts coming from different particles within the measuring volume. Moments of the streamwise velocity showed a slight dependence on measuring volume length, indicating that spatial averaging effects well known for hot-films and hot-wires can occur in laser velocimetry measurements when the data density is high.List of symbols time-averaged quantity - u wall friction velocity, ( _w /)^1/2 - v kinematic viscosity - d _p pinhole diameter - l _eff spanwise extent of LDV measuring volume viewed by photomultiplier - l ⁺ non-dimensional length of measuring volume, l _eff u /v - y ⁺ non-dimensional coordinate in spanwise direction, y u /v - z ⁺ non-dimensional coordinate in spanwise direction, z u /v - U ⁺ non-dimensional mean velocity, /u - u instantaneous streamwise velocity fluctuation, U &#x2329;U - v instantaneous normal velocity fluctuation, V–V - u RMS streamwise velocity fluctuation, u ^21/2 - v RMS normal velocity fluctuation, v ^21/2 - Re Reynolds number based on momentum thickness, U 0/v - R _uv cross-correlation coefficient, u v/u v - R₁₂(0, 0, z) two point correlation between u and v with z-separation, <u(0, 0, 0) v (0, 0, z)>/<u(0, 0, 0) v (0, 0, 0)> - N rate at which bursts are validated by counter processor - T Taylor time microscale, u (dv/dt²)^–1/2     

Cycle-to-cycle variation effects on turbulent shear stress measurements in pulsatile flows

Accurate evaluation of turbulent velocity statistics in pulsatile flows is important in estimating potential damage to blood constituents from prosthetic heart valves. Variations in the mean flow from one cycle to the next can result in artificially high estimates. Here we demonstrate a procedure using a digital, low-pass filter to remove the cycle-to-cycle variation from turbulence statistics. The results show that cycle-to-cycle variations can significantly affect estimates of turbulent Reynolds stress and should be either eliminated or demonstrated to be small when reporting pulsatile flow results.List of symbols D inside diameter of aortic valve - R radius of model aorta - t time window - t time - T period of cycle - T duration of outflow pulse from ventricle - U instantaneous axial velocity - U _L low-pass axial velocity - U _p mean periodic axial velocity - U ensemble averaged axial velocity - uv ensemble-average turbulent velocity product - u root-mean-square of turbulent axial velocity - U _max maximum, ensemble-averaged axial velocity - V instantaneous radial velocity - y vertical distance from aorta centerline - z axial distance downstream of prosthetic heart valve This paper was presented at the Tenth Symposium on Turbulence, University of Missouri-Rolla, Sept. 22–24, 1986     

Numerical computation of high Rayleigh number natural convection and prediction of hot radiator induced room air motion

The two-dimensional stationary turbulent buoyant flow and heat transfer in a cavity at high Rayleigh numbers was computed numerically. The k– turbulence model was used. The time-averaged equations for momentum, energy and continuity, which are coupled to the turbulence equations, were solved using a finite difference formulation. In order to validate the computer code, a comparison exercise was carried out. The test results are in good agreement with the internationally accepted benchmark solution. Grid-refinement shows the necessity of a very fine grid at high Rayleigh numbers with especially small grid-distances in the near-wall region. The computed boundary layer velocity profiles are in excellent agreement with available experimental data. The local heat transfer in the turbulent part of the boundary layers is predicted 20% too high. Computations were carried out for the natural convective flow in a room induced by a hot radiator and a cold window. Various radiator configurations and types of thermal boundary conditions were applied including thermal radiation interaction between surfaces.Nomenclature a thermal diffusivity (m²/s) - C constant in _t expression - D cavity dimensions (m) - g acceleration of gravity (m/s²) - G _k production/destruction of k by buoyancy (kg/ms³) - h enthalpy (J/kg) - IX index of grid point - k turbulent kinetic energy (m²/s²) - m dimensionless stratification parameter - Nu overall Nusselt number - Nu _y local Nusselt number - NX total number of grid points - p pressure (N/m²) - P _k production of k by shear stress (kg/ms³) - Q heat flux through wall (W/m) - Ra overall Rayleigh number - Ra _y local Rayleigh number - Re _t turbulent Reynolds number - S source term in -equation (kg/ms⁴) - S source term for - T _c, T _h temperatures of cold and hot walls (K) - T _s (y) stratification temperature on vertical mid-line (K) - T ₀ mean cavity temperature (K) - u, v horizontal and vertical velocity components (m/s) - u ₀ Brunt-Vaisälä velocity scale (m/s) - x, y horizontal and vertical coordinates (m) - non-linearity parameter for grid - coefficient of thermal expansion (l/K) - jet angle (°) - diffusivity for - S dissipation rate for turbulent kinetic energy (m²/s³) - variable to be solved - thermal conductivity (W/mK) - , _t kinematic and eddy viscosities (m²/s) - stream function (kg/ms) - density (kg/m³) - _k, , _t constants in k– model     

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     

Constituency measurements in the mixing region of a cross flow jet using a laser velocimeter

The velocities in the mixing region of a cross flow jet injected into a freestream were studied in detail with a laser velocimeter. Three jet to freestream momentum ratios were used (3.1, 8.1, 16.2). By purposely seeding the jet and freestream separately (as well as both simultaneously), marking the fluid was feasible. Thus, determining the velocities that emanated from the different streams was possible. By methodically analyzing the three sets of dependent data, the size and location of the mixing region was determined. The mixing regions for the three momentum ratios were found to be of different sizes and at different locations. By proper scaling, however, the regions for the three momentum ratios were found to collapse to one scaled region. Because of the intermittent behavior of the mixing, conventional turbulence models for such mixing may not be applicable; however, detailed velocities and turbulence quantities are included for benchmarking predictions.List of symbols B slot width - H channel height - MR momentum ratio, jet to free stream = _j V _j ²/ U ² - Re _H Reynolds number, U H/v - U free stream velocity - u axial velocity - u rms of axial velocity fluctuation - v transverse velocity - v rms of transverse velocity fluctuation - V _j slot exit transverse velocity - x axial direction (Fig. 3) - x _c x-center of mixing region - scaled value of x, = x/B - y transverse direction (Fig. 3) - y _c y-center of mixing region - scaled value of y, = y/ MRB - _x mixing region width in x-direction - _y mixing region width in y-direction - scaled mixing region width in x-direction, = _x /B - scaled mixing region width in y-direction, = _y / MRB - free stream density - _j slot exit density - v kinematic viscosity of freestream This research was sponsored in part by the Fulbright Commission (Bonn, Germany), the Institut für Thermische Strömungsmaschinen, Universität Karlsruhe (Karlsruhe, Germany), and the Rotating Machinery and Controls Industrial Research Program, University of Virginia (Charlottesville, VA, USA)     

Approximate solutions for drag coefficient of bubbles moving in shear-thinning elastic fluids

The drag coefficient for bubbles with mobile or immobile interface rising in shear-thinning elastic fluids described by an Ellis or a Carreau model is discussed. Approximate solutions based on linearization of the equations of motion are presented for the highly elastic region of flow. These solutions are in reasonably good agreement with the theoretical predictions based on variational principles and with published experimental data. C _D Drag coefficient - E ^* Differential operator [E ^{* 2} = ²/² + (sin/ ²)/(1/sin /)] - El Ellis number - F _D Drag force - K Consistency index in the power-law model for non-Newtonian fluid - n Flow behaviour index in the Carreau and power-law models - P Dimensionless pressure [=(p – p ₀)/₀ (U /R)] - p Pressure - R Bubble radius - Re ₀ Reynolds number [= 2R U /₀] - Re Reynolds number defined for the power-law fluid [= (2R) ⁿ U ^2–n /K] - r Spherical coordinate - t Time - U Terminal velocity of a bubble - u Velocity - Wi Weissenberg number - Ellis model parameter - Rate of deformation - Apparent viscosity - ₀ Zero shear rate viscosity - Infinite shear rate viscosity - Spherical coordinate - Parameter in the Carreau model - ^* Dimensionless time [=/(U /R)] - Dimensionless length [=r/R] - Second invariant of rate of deformation tensors - ^* Dimensionless second invariant of rate of deformation tensors [=/(U /R)²] - Second invariant of stress tensors - ^* Dimensionless second invariant of second invariant of stress tensor [= / ₀ ² (U /R)²] - Fluid density - Shear stress - ^* Dimensionless shear stress [=/ ₀ (U /R)] - _1/2 Ellis model parameter - ₁ ^2/* Dimensionless Ellis model parameter [= _1/2/ ₀(U /R)] - Stream function - ^* Dimensionless stream function [=/U R ²]     

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     

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