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

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

Transport in ordered and disordered porous media I: The cellular average and the use of weighting functions

In this work we consider transport in ordered and disordered porous media using singlephase flow in rigid porous mediaas an example. We defineorder anddisorder in terms of geometrical integrals that arise naturally in the method of volume averaging, and we show that dependent variables for ordered media must generally be defined in terms of thecellular average. The cellular average can be constructed by means of a weighting function, thus transport processes in both ordered and disordered media can be treated with a single theory based on weighted averages. Part I provides some basic ideas associated with ordered and disordered media, weighted averages, and the theory of distributions. In Part II a generalized averaging procedure is presented and in Part III the closure problem is developed and the theory is compared with experiment. Parts IV and V provide some geometrical results for computer generated porous media.Roman Letters A interfacial area of the- interface contained within the macroscopic region, m² - A_e area of entrances and exits for the-phase contained within the macroscopic system, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor, m² - L general characteristic length for volume averaged quantities, m - characteristic length (pore scale) for the-phase - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - N unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p0 reference pressure in the-phase, N/m² - p traditional intrinsic volume averaged pressure, N/m² - r0 radius of a spherical averaging volume, m - r position vector, m - r position vector locating points in the-phase, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - v velocity vector in the-phase, m/s - v traditional superficial volume averaged velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V/V, volume average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

Flow visualization of compressible vortex structures using density gradient techniques

Mathematical results are derived for the schlieren and shadowgraph contrast variation due to the refraction of light rays passing through two-dimensional compressible vortices with viscous cores. Both standard and small-disturbance solutions are obtained. It is shown that schlieren and shadowgraph produce substantially different contrast profiles. Further, the shadowgraph contrast variation is shown to be very sensitive to the vortex velocity profile and is also dependent on the location of the peak peripheral velocity (viscous core radius). The computed results are compared to actual contrast measurements made for rotor tip vortices using the shadowgraph flow visualization technique. The work helps to clarify the relationships between the observed contrast and the structure of vortical structures in density gradient based flow visualization experiments.Nomenclature a Unobstructed height of schlieren light source in cutoff plane, m - c Blade chord, m - f Focal length of schlieren focusing mirror, m - C _T Rotor thrust coefficient, T/( ² R ⁴) - I Image screen illumination, Lm/m ² - l Distance from vortex to shadowgraph screen, m - n _b Number of blades - p Pressure,N/m ² - p Ambient pressure, N/m ² - r, , z Cylindrical coordinate system - r _c Vortex core radius, m - Non-dimensional radial coordinate, (r/r _c ) - R Rotor radius, m - Tangential velocity, m/s - Specific heat ratio of air - Circulation (strength of vortex), m ²/s - Non-dimensional quantity, ² 8²p r _c ² - Refractive index of fluid medium - ₀ Refractive index of fluid medium at reference conditions - Gladstone-Dale constant, m ³/kg - Density, kg/m ³ - Density at ambient conditions, kg/m ³ - Non-dimensional density, (/ ) - Rotor solidity, (n _b c/ R) - Rotor rotational frequency, rad/s     

Convective heat transfer within fibrous insulation slabs

A numerical study of convective heat flow within a fibrous insulating slab is presented. The material is treated as an anisotropic porous medium and the variation of properties with temperature is taken into account. Good agreement is obtained with available experimental data for the same geometry.

---

Zusammenfassung Für den konvektiven Wärmestrom in einem faserförmigen Isolierstoff wird eine numerische Berechnung angegeben. Der Stoff wird als anisotropes poröses Medium mit temperaturabhängigen Stoffwerten angesehen. Die Übereinstimmung mit verfügbaren Versuchswerten ist gut.

---

Nomenclature C_p specific heat of the gas at the mean temperature - Da Darcy number=k_y/H² - Gr^* modified Grashof number=gTHk_y/²= (Grashof number) × (Darcy number) - H thickness of the specimen - P gas pressure - Pr^* modified Prandtl number= C_p/_x - Ra^* modified Rayleigh number=Gr^* Pr^* - R_p ratio of permeabilities=k_y/k_x - R_k ratio of conductivities= _y/_x - T absolute temperature of the gas - t₁ absolute temperature of the hot face - T₂ absolute temperature of the cold face - T_m mean temperature of the gas=(T₁+T₂)/2 - k_x specific permeability of the porous medium along the x-direction - k_y specific permeability of the porous medium along the y-direction - p T/T_m - q exponent - r exponent - u gas velocity along the x-direction - v gas velocity along the y-direction - X_* distance along the x-direction - y_* distance along the y-direction - T temperature difference=t₁–T₂ - thermal coefficient of expansion of the gas - _m thermal coefficient of expansion of the gas at the mean temperature - _* T–T_m - dimensionless temperature= ^*/T - _a apparent thermal conductivity of the porous medium along the x-direction - _al local apparent thermal conductivity of the porous medium along the x-direction - _x thermal conductivity of the porous medium along the x-direction in the absence of convection - _y thermal conductivity of the porous medium along the y-direction in the absence of convection - dynamic viscosity of the gas - _m dynamic viscosity of the gas at the mean temperature - kinematic viscosity of the gas - _m kinematic viscosity of the gas at the mean temperature - density of the gas - _m density of the gas at the mean temperature - ^* stream function at any point - dimensionless stream function= ^*/( _m/_m)     

Free convection heat transfer from an isothermal plate with arbitrary inclination

This paper presents a new formulation for the laminar free convection from an arbitrarily inclined isothermal plate to fluids of any Prandtl number between 0.001 and infinity. A novel inclination parameter is proposed such that all cases of the horizontal, inclined and vertical plates can be described by a single set of transformed equations. Moreover, the self-similar equations for the limiting cases of the horizontal and vertical plates are recovered from the transformed equations by setting=0 and=1, respectively. Heated upward-facing plates with positive and negative inclination angles are investigated. A very accurate correlation equation of the local Nusselt number is developed for arbitrary inclination angle and for 0.001 Pr .

---

Wärmeübertragung bei freier Konvektion an einer isothermen Platte mit beliebiger Neigung

Zusammenfasssung Diese Untersuchung stellt eine neue Formulierung der laminaren freien Konvektion von Flüssigkeiten mit einer Prandtl-Zahl zwischen 0,001 und unendlich an einer beliebig schräggestellten isothermen Platte dar. Ein neuer Neigungsparameter wird eingeführt, so daß alle Fälle der horizontalen, geneigten oder vertikalen Platte von einem einzigen Satz transformierter Gleichungen beschrieben werden können. Die unabhängigen Gleichungen für die beiden Fälle der horizontalen and vertikalen Platte wurden für=0 und=1 aus den transformierten Gleichungen wieder abgeleitet. Es wurden erwärmte aufwärtsgerichtete Platten mit positiven und negativen Neigungswinkeln untersucht. Eine sehr genaue Gleichung wurde für die lokale Nusselt-Zahl bei beliebigen Neigungswinkeln und für 0,001 Pr entwickelt.

---

Nomenclature C _p specific heat - f reduced stream function - g gravitational acceleration - Gr local Grashof number,g(T _w –T _w ) x³/v² - h local heat transfer coefficient - k thermal conductivity - n constant exponent - Nu local Nusselt number,hx/k - p pressure - Pr Prandtl number, v/ - Ra local Rayleigh number,g(T _w –T )J x³/v - T fluid temperature - T _w wall temperature - T temperature of ambient fluid - u velocity component in x-direction - v velocity component in y-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - (Ra¦sin¦)^1/4/( Ra cos()^1/5 - pseudo-similarity variable, (y/) - dimensionless temperature, (T–T )/(T _w–T ) - ( Ra cos)^1/5+(Rasin)^1/4 - v kinematic viscosity - 1/[1 +(Ra cos)^1/5/( Ra¦sin)^1/4] - density of fluid - Pr/(1+Pr) - _w wall shear stress - angle of plate inclination measured from the horizontal - stream function - dimensionless dynamic pressure     

Perturbation analysis for periodic heat transfer in radiating fins

A perturbation analysis is presented for periodic heat transfer in radiating fins of uniform thickness. The base temperature is assumed to oscillate around a mean value. The perturbation expansion is carried out in terms of dimensionless amplitude of the base temperature oscillation. The zero-order problem which is nonlinear, and corresponds to the steady state fin behaviour, is solved by quasilinearization. A method of complex combination is used to reduce both the first and the second order problems to two, coupled linear boundary value problems which are subsequently solved by a noniterative numerical scheme. The second-order term is composed of an oscillatory component with twice the frequency of base temperature oscillation and a time-independent term which causes a net change in the steady state values of temperature and heat transfer rate. Within the range of parameters used, the net effect is to decrease the mean temperature and increase the mean heat transfer rate. This is in constrast to the linear case of convecting fins where the mean values are unaffected by base temperature oscillations. Detailed numerical results are presented illustrating the effects of fin parameter N and dimensionless frequency B on temperature distribution, heat transfer rate, and time-average fin efficiency. The time-average fin efficiency is found to reduce significantly at low N and high B.

---

Störungsanalyse für periodische Wärmeübertragung an Strahlungsrippen

Zusammenfassung Eine Störungsanalyse wird für periodische Wärmeübertragung in Strahlungsrippen gleicher Dicke vorgelegt. Die Fußtemperatur wird als um einen Mittelwert schwingend angenommen. Die Störungsentwicklung wird in Termen einer dimensionslosen Amplitude e dieser Schwingung angesetzt. Das Problem nullter Ordnung, das nichtlinear ist und dem stationären Verhalten der Rippe entspricht, wird durch Quasilinearisierung gelöst. Eine Methode der komplexen Kombination wird angewandt, um die Probleme erster und zweiter Ordnung auf zwei gekoppelte Grenzwertprobleme zu reduzieren, die nacheinander nach einem nichtiterativen Schema gelöst werden. Der Term zweiter Ordnung besteht aus einer Schwingungskomponente mit der doppelten Frequenz der Schwingung der Fußtemperatur und einem zeitunabhängigen Term, der eine Nettoänderung der stationären Werte der Temperatur und der Wärmeübertragung verursacht. Im verwendeten Bereich der Parameter tritt eine Abnahme der mittleren Temperatur und eine Zunahme der mittleren Wärmeübertragung auf. Das steht im Gegensatz zum linearen Fall der Konvektionsrippe, bei dem die Mittelwerte durch Schwingungen der Fußtemperatur nicht beeinflußt werden. Detaillierte numerische Ergebnisse zeigen die Einflüsse des Rippenparameters N und der dimensionslosen Frequenz B auf Temperatur Verteilung, Wärmeübertragung und zeitliches Mittel des Rippengütegrades. Dieses zeitliche Mittel nimmt merklich ab bei kleinem N und hohem B.

---

Nomenclature b fin thickness - B dimensionless frequency, L²/ - E emissivity - f₀, f₁ functions of X - g₀, g₁, g₂ functions of X - h₀, h₁, h₂ functions of X - k thermal conductivity - L fin Length - N fin parameter, 2EL²T_bm/bk - q heat transfer rate - Q dimensionless heat transfer rate, qL/kbT_bm - t time - T temperature - T_b fin base temperature - T_S effective sink temperature - T_bm mean fin base temperature - x axial distance - X dimensionless axial distance, x/L - dimensionless amplitude of base temperature (s. Eq.2) - thermal diffusivity - instantaneous fin efficiency - time-average fin efficiency - _ss steady state fin efficiency - dimensionless temperature, T/T_bm - ₀ zero-order approximation - ₁ first-order approximation - ₂ second-order approximation - _2s steady component of ₂ - , ₁, ₂ constants - complex function of X - ₁ real part of - ₂ imaginary part of - complex function of X - ₁ real part of Y - ₂ imaginary part of - dimensionless time, t/L² - frequency of base temperature oscillation     

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

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

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

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

Transport in ordered and disordered porous media II: Generalized volume averaging

In this paper we develop the averaged form of the Stokes equations in terms of weighting functions. The analysis clearly indicates at what point one must choose a media-specific weighting function in order to achieve spatially smoothed transport equations. The form of the weighting function that produces the cellular average is derived, and some important geometrical theorems are presented.Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - A _p surface area of a particle, m² - d _p 6V _p/A_p, effective particle diameter, m - g gravity vector, m/s² - I unit tensor - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p general characteristic length for volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - l characteristic length (pore scale) for the-phase - l _i i=1, 2, 3 lattice vectors, m - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - m _D special convolution product weighting function for disordered media - m _M master convolution product weighting function for ordered and disordered media - n unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p_m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p p p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function, m - r position vector, m - r position vector locating points in the-phase, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - V velocity vector in the-phase, m/s - vm superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - V volume of the-phase contained in the averaging volume, m³ - V _p volume of a particle, m³ - v traditional superficial volume averaged velocity, m/s - v v p _m spatial deviation velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - _m m * . weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m² - V /V, volume fraction of the-phase     

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

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

Mixed convection from a horizontal plate to fluids of any Prandtl number

Laminar mixed convection over a horizontal plate with uniform wall temperature or uniform wall heat flux is analyzed by introducing proper buoyancy parameters and transformation variables for fluids of any Prandtl number between 0.001 and 10,000. Both cases of buoyancy assisting and opposing flow conditions are investigated. For the buoyancy-assisting case, the obtained numerical results are very accurate over the entire range of mixed convection intensity from pure forced convection limit to pure free convection limit. For the buoyancy-opposing case, solutions are obtained from the forced convection limit to the point of breakdown.

---

Mischkonvektion an einer horizontalen Platte für Fluide mit beliebiger Prandtl-Zahl

Zusammenfassung Es wurde laminare Mischkonvektion an einer horizontalen Platte mit einheitlicher Wandtemperatur oder einheitlicher Wandwärmestromdichte bei Einführung zweckmäßiger Auftriebsparameter und Transformationsvariablen für Fluide mit beliebiger Prandtl-Zahl zwischen 0,001 und 10 000 untersucht. Es wurden die Fälle der Strömung entgegen und in Richtung der Auftriebskraft untersucht. Für den Fall der Strömung in Richtung der Auftriebskraft wurden sehr genaue numerische Ergebnisse für den gesamten Bereich der gemischten Konvektion von rein erzwungener Konvektion bis zu rein freier Konvektion erhalten. Für den Fall der Strömung entgegen der Auftriebsrichtung wurden Lösungen für erzwungene Konvektion bis zum Umkehrpunkt erhalten.

---

Nomenclature C _f local friction coefficient - f reduced stream function - g gravitational acceleration - Gr local Grashof number for UWT,g (T _w –T )x ³/ ² - Gr* local Grashof number for UHF,g q _w x ⁴/k ² - m =10 for UWT; and =6 for UHF - n =5 for UWT; and =3 for UHF - Nu local Nusselt number - p pressure - Pr Prandtl number,/ - q _w wall heat flux - Ra local Rayleigh number for UWT,Gr Pr - Ra* local Rayleigh number for UHF,Gr*Pr - Re local Reynolds number,u x/ - T fluid temperature - T _w wall temperature - T free-stream temperature - u velocity component inx-direction - u free-stream velocity - v velocity component iny-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - =0 for UWT; and =1 for UHF - buoyancy parameter, =( Ra)^1/5/( Re)^1/2 for UWT; and =( Ra*)^1/6/( Re)^1/2 for UHF - pseudo-similarity variable, (y/x) - dimensionless temperature, =(T–T )/(T _w –T ) for UWT; and =(T–T )/(q _w x/k) for UHF - =[( Re)^1/2+( Ra)^1/5] for UWT; and =[( Re)^1/2+( Ra*)^1/6] for UHF - dynamic viscosity - kinematic viscosity - /(1+) - dimensionless pressure - density - Pr/(1+Pr) - _w wall shear stress,(u/y) _y=0 - stream function - Pr/(1+Pr)^1/3     

On laminar flow through a uniformly porous pipe

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

Nomenclature

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

Instationäre Messung der Wärmeleitfähigkeit mit optischer Registrierung

Zusammenfassung Die Wärmeleitfähigkeit von Wasser wird im Temperaturbereich von 20 bis 90 °C und bei 1 bar mit einem neuen instationären Absolutverfahren bestimmt. Zur Aufzeichnung des instationären Feldes des Brechungsindex und der Temperatur werden zwei interferometrische Anordnungen benutzt: Die Methode vonMach-Zehnder und das Biprismaverfahren. Die Ergebnisse stehen in guter Übereinstimmung zu den Messungen anderer Autoren, die stationäre Methoden benutzten. Die Unsicherheit der 374 Einzelmessungen wird auf höchstens±1% geschätzt. Damit ist nachgewiesen, daß ein instationäres Meßverfahren mit optischer Registrierung mit den klassischen stationären Verfahren hinsichtlich der Meßunsicherheit konkurrieren kann. Das instationäre Verfahren kommt ohne kalorische Messungen aus und besteht bei optischer Registrierung im wesentlichen aus Längenmessungen.

---

Unsteady-state measurements of the thermal conductivity with optical recording

The thermal conductivity of water in the temperature region from 20 to 90°C and at 1 bar was measured by means of a new unsteady-state absolute method. To record the unsteady-state field of the index of refraction and of the temperature, two interferometric arrangements were used: TheMach-Zehnder and the biprisma methods. The results are in good agreement with measurements of other authors, who had used steady-state methods. The maximum degree of uncertainty of the 374 measurements is estimated to be±1%. Thus it is shown that unsteady-state methods with optical recording can well be compared with classical steady-state methods regarding uncertainties. The method does not require caloric measurements and uses primarly determinations of lengths.

---

Bezeichnungen A=2 a in den Auswerteverfahren gebrauchte Abkürzung, m - B=q _x=0/ in den Auswerteverfahren gebrauchte Abkürzung, grd/m - F Fläche, m² - Fo=a/x ² Fourierzahl - I Strom, A - R Widerstand, - U Spannung, V - T Kelvintemperatur, °K - a=/c _p Temperaturleitfähigkeit, m²/s - Wärmeeindringzahl, Ws^1/2/m²grd - c _p isobare spezifische Wärmekapazität, kJ/kggrd - l Modellänge, m - n Brechungsindex - q Wärmestromdichte, W/m² - t Celsiustemperatur, °C - =1/ spezifisches Volumen, m³/kg - x Wandabstand, Ortskoordinate, m - z Ordnungszahl der Interferenzstreifen - =/_x=0 dimensionslose Übertemperatur - n=n -n _x=0 Brechungsindexdifferenz - Verhältnis der Wärmeeindringzahlen - Übertemperatur, grd - Wärmeleitfähigkeit, W/mgrd - Lichtwellenlänge, m - =1/ Dichte, kg/msu3 - Zeit, s Indices Zustand des Bades, Umgebung - x=0 Wand - z Stelle der z-ten Ordnung - i Laufparameter - err errechneter Wert - mess Meßwert - Bez Bezug Auszug aus der von der Fakultät für Maschinenwesen, und Elektrotechnik der Technischen Hochschule München genehmigten Dissertation von J.Bach.     

Mixed convection flow in narrow vertical ducts

The mixed convection flow in a vertical duct is analysed under the assumption that , the ratio of the duct width to the length over which the wall is heated, is small. It is assumed that a fully developed Poiseuille flow has already been set up in the duct before heat from the wall causes this to be changed by the action of the buoyancy forces, as measured by a buoyancy parameter . An analytical solution is derived for the case when the Reynolds numberRe, based on the duct width, is of 0 (1). This is extended to the case whenRe is 0 (^–1) by numerical integrations of the governing equations for a range of values of representing both aiding and opposing flows. The limiting cases, || 1 andR=Re of 0 (1), andR and both large, with of 0 (R ^1/3) are considered further. Finally, the free convection limit, large with R of 0 (1), is discussed.

---

Mischkonvektion in engen senkrechten Rohren

Zusammenfassung Mischkonvektion in einem senkrechten Rohr wird unter der Voraussetzung untersucht, daß das Verhältnis der Rohrbreite zur Länge, über welche die Wand beheizt wird, klein ist. Es wird angenommen, daß sich bereits eine voll entwickelte Poiseuille-Strömung in dem Rohr eingestellt hat, bevor Antriebskräfte, gemessen mit dem Auftriebsparameter , aufgrund der Wandbeheizung die Strömung verändern. Es wird eine analytische Lösung für den Fall erhalten, daß die mit der Rohrbreite als charakteristische Länge gebildete Reynolds-ZahlRe konstant ist. Dies wird mittels einer numerischen Integration der wichtigsten Gleichungen auf den FallRe =f (^–1) sowohl für Gleich- als auch für Gegenstrom ausgedehnt. Weiterhin werden die beiden Grenzfälle betrachtet, wenn || 1 undR=Re konstant ist, sowieR und beide groß mit proportionalR ^1/3. Schließlich wird der Grenzfall der freien Konvektion, großes mit konstantem R, diskutiert.

---

Nomenclature g acceleration due to gravity - Gr Grashof number - G modified Grashof number - h duct width - l length of the heated section of the duct wall - p pressure - Pr Prandtl number - Q flow rate through the duct - Q ₀ heat transfer on the wally=0 - Q ₁ heat transfer on the wally=1 - Re Reynolds number - R modified Reynolds number - T temperature of the fluid - T ₀ ambient temperature - T applied temperature difference - u, velocity component in thex-direction - v, velocity component in they-direction - x, co-ordinate measuring distance along the duct - y, co-ordinate measuring distance across the duct - buoyancy parameter - ₀ modified buoyancy parameter, ₀=R ^–1/3 - coefficient of thermal expansion - ratio of duct width to heated length, =h/l - (non-dimensional) temperature - _w applied temperature on the wally=0 - kinematic viscosity - density of the fluid - ₀ shear stress on the wally=0 - ₁ shear stress on the wally=1 - stream function     

Flow through a helically coiled annulus

In this paper the flow is studied of an incompressible viscous fluid through a helically coiled annulus, the torsion of its centre line taken into account. It has been shown that the torsion affects the secondary flow and contributes to the azimuthal component of velocity around the centre line. The symmetry of the secondary flow streamlines in the absence of torsion, is destroyed in its presence. Some stream lines penetrate from the upper half to the lower half, and if is further increased, a complete circulation around the centre line is obtained at low values of for all Reynolds numbers for which the analysis of this paper is valid, being the ratio of the torsion of the centre line to its curvature.Nomenclature A =constant - a outer radius of the annulus - b unit binormal vector to C - C helical centre line of the pipe - D rL - g 1000 - K Dean number=Re² - L 1+r sin - M (L ²+ ² r ²)^1/2 - n unit normal vector to C - P, P pressure and nondimensional pressure - p ₀, p pressures of O(1) and O() - Re Reynolds number=aW ₀/ - (r, , s), (r, , s) coordinates and nondimensional coordinates - nonorthogonal unit vectors along the coordinate directions - r ₀ radius of the projection of C - t unit tangent vector to C - V _r, V , V _s velocity components along the nonorthogonal directions - V_r, V, V _s nondimensional velocity components along - W ₀ average velocity in a straight annulus Greek symbols , curvature and nondimensional curvature of C - U, V, W lowest order terms for small in the velocity components along the orthogonal directions t - _r, , _s first approximations to V _r , V, V _s for small - =/=/ - kinematic viscosity - density of the fluid - , torsion and nondimensional torsion of C - , stream function and nondimensional stream function - nondimensional streamfunction for U, V - a inner radius of the annulus After this paper was accepted for publication, a paper entitled On the low-Reynolds number flow in a helical pipe, by C.Y. Wang, has appeared in J. Fluid. Mech., Vol 108, 1981, pp. 185–194. The results in Wangs paper are particular cases of this paper for =0, and are also contained in [9].     

On the universality of the velocity profiles of a turbulent flow in an axially rotating pipe

If a fluid enters an axially rotating pipe, it receives a tangential component of velocity from the moving wall, and the flow pattern change according to the rotational speed. A flow relaminarization is set up by an increase in the rotational speed of the pipe. It will be shown that the tangential- and the axial velocity distribution adopt a quite universal shape in the case of fully developed flow for a fixed value of a new defined rotation parameter. By taking into account the universal character of the velocity profiles, a formula is derived for describing the velocity distribution in an axially rotating pipe. The resulting velocity profiles are compared with measurements of Reich [10] and generally good agreement is found.Nomenclature b constant, equation (34) - D pipe diameter - l mixing length - l ₀ mixing length in a non-rotating pipe - N rotation rate,N=Re /Re _D - p pressure - R pipe radius - Re _D flow-rate Reynolds number, - Re rotational Reynolds number, Re =v _w D/ - Re* Reynolds number based on the friction velocity, Re*=v*R/ - (Re*)₀ Reynolds number based on the friction velocity in a non-rotating pipe - Ri Richardson number, equation (10) - r coordinate in radial direction - dimensionless coordinate in radial direction, - v _r ,v ,v _z time mean velocity components - v _r ,v ,v _z velocity fluctations - v _w tangential velocity of the pipe wall - v* friction velocity, - axial mean velocity - v _ZM maximum axial velocity - dimensionless radial distance from pipe wall, - y ⁺ dimensionless radial distance from pipe wall - y ₁ ⁺ constant - Z rotation parameter,Z =v _w/v _* =N Re _D /2Re* - _m eddy viscosity - ( _m )₀ eddy viscosity in a non-rotating pipe - coefficient of friction loss - von Karman constant - ₁ constant, equation (31) - density - dynamic viscosity - kinematic viscosity     

Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A ^* interfacial area of the- interface contained within a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

Deformation to breaking of deep water gravity waves

The results of laboratory observations of the deformation of deep water gravity waves leading to wave breaking are reported. The specially developed visualization technique which was used is described. A preliminary analysis of the results has led to similar conclusions than recently developed theories. As a main fact, the observed wave breaking appears as the result of, first, a modulational instability which causes the local wave steepness to approach a maximum and, second, a rapidly growing instability leading directly to the breaking.List of symbols L total wave length - H total wave height - crest elevation above still water level - trough depression below still water level - wave steepness =H/L - crest steepness =/L - trough steepness =/L - F ₁ forward horizontal length from zero-upcross point (A) to wave crest - F ₂ backward horizontal length from wave crest to zero-downcross point (B) - crest front steepness =/F ₁ - crest rear steepness =/F ₂ - vertical asymmetry factor=F ₂/F ₁ (describing the wave asymmetry with respect to a vertical axis through the wave crest) - µ horizontal asymmetry factor=/H (describing the wave asymmetry with respect to a horizontal axis: SWL) - T ₀ wavemaker period - L ₀ theoretical wave length of a small amplitude sinusoïdal wave generated at T _inf0 ^sup–1 frequency - ₀ average wave height     

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     

Free convection from a vertical plate with concentrated and distributed thermal sources

An analysis is developed for the laminar free convection from a vertical plate with uniformly distributed wall heat flux and a concentrated line thermal source embedded at the leading edge. We introduce a parameter=(1 +Q _L/Q_w)^–1=(1 + Ra_L/Ra_w)^–1 to describe the relative strength of the two thermal sources; and propose a unified buoyancy parameter=( Ra_L+ Ra_w)^1/5 with=1/(1 +Pr ^–1) to properly scale the dependent and independent variables. The variables are so defined that the resulting nonsimilar boundary-layer equations can describe exactly the buoyancy-induced flow from the dual sources with any relative strength to fluids of any Prandtl number from very small values to infinity. These nonsimilar equations are readily reducible to the self-similar equations of an adiabatic wall plume for=0, and to those of free convection from uniform flux plate for=1. Rigorous finite-difference solutions for fluids of Pr from 0.001 to are obtained over the entire range of from 0 to 1. The effects of both relative source strength and Prandtl number on the velocity profiles, temperature profiles, and the variations of wall temperature, are clearly illustrated.

---

Freie Konvektion an einer vertikalen Platte mit einer konzentrierten und einer gleichmäßig verteilten Wärmequelle

Zusammenfassung Für die freie Konvektion an einer vertikalen Platte mit einer gleichmäßig verteilten Wandwärmestromdichte und einer in der Vorderkante eingebetteten linienförmigen Wärmequelle wird eine Berechnungsmethode entwickelt. Zur Beschreibung der relativen Stärke der beiden Wärmequellen führen wir einen Parameter=(1 + Q_L/Q_w)^–1=(1 + Ra_L/Ra_w)^–1 ein und schlagen einen vereinheitlichten Auftriebsparameter=( Ra _L+ Ra _w)^1/5 mit=1/(1 +Pr ^–1 für die Skalierung der abhängigen und unabhängigen Variablen vor. Die Variablen werden so definiert, daß mit den sich ergebenden unabhängigen Grenzschichtgleichungen die von den beiden Wärmequellen beliebiger Stärke verursachte Auftriebsströmung von Fluiden beliebiger Prandtl-Zahl genau beschrieben werden kann. Diese unabhängigen Gleichungen können ohne weiteres auf die selbstähnlichen Gleichungen für den Fall einer lokalen Wärmezufuhr an einer sonst adiabatischen Wand für=0 und jenen der freien konvektion an einer Platte mit einheitlichem Wärmestrom für=1 zurückgeführt werden. Für Fluide mit der Prandtl-Zahl zwischen 0,001 und Unendlich werden nach der strengen finite Differenzen-Methode Lösungen im Bereich von zwischen 0 und 1 erhalten. Der jeweilige Einfluß der relativen Quellenstärke und der Prandtl-Zahl auf die Geschwindigkeits- und Temperaturprofile sowie die Veränderung der Wandtemperatur werden deutlich dargestellt.

---

Nomenclature C _f friction coefficient - C _p specific heat - f reduced stream function - g gravitational acceleration - k thermal conductivity - L width of the plate - Nu local Nusselt number - Pr Prandtl number - q _w wall heat flux - Q _L heat generated by the line source - Q _w heat released by the uniform-flux wall from 0 tox, q _w Lx - Ra _L local Rayleigh number, g T _L ^* x ³/( ) - Ra _w local Rayleigh number,g T _w ^* w ³/( ) - T fluid temperature - T temperature of ambient fluid - T _L ^* characteristic temperature of the line source,Q _{L/(C p} L) - T _w ^* characteristic temperature of the uniform flux wall, =q _w x/k=Q _w /(C _p L) - u velocity component in then-direction - U₀ dimensionless velocity,u/(/x) Ra _L ^2/5 - U ₁ dimensionless velocity,u/(/x) Ra _w ^2/5 - velocity component in they-direction - x coordinate parallel to the plate - y coordinate normal to the plate - thermal diffusivity - thermal expansion coefficient - pseudo-similarity variable,(y/x) - dimensionless temperature, (T–T )/(T _L ^* +T _w ^* ) - ₀ dimensionless temperature, (Ra_l)^1/5 (T–T )/T _L ^* - ₁ dimensionless temperature, (Ra_w)Ra_w)^1/5 (T–T )/T _w ^* - (Ra _L+Ra_w)^1/5 - kinematic viscosity - (1 +Ra _L/Ra_w)^–1=(1 +T _L ^* /T _w ^* )^–1=(1 + Q_L/Q_w)^–1 - density - Pr/(1 +Pr) - _w wall shear stress - stream function