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

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

On a Boundary-Value Problem of Periodic Type for First-Order Linear Functional Differential Equations

We establish unimprovable sufficient conditions for the unique solvability of the boundary-value problem

Pulsatile flow in a tube with wall injection and suction

This paper studies similarity solutions for pulsatile flow in a tube with wall injection and suction. The Navier-Stokes equations are reduced to a system of three ordinary differential equations. Two of the equations represent the effects of suction and injection on the steady flow while the third represents the effects of suction and injection on pulsatile flow. Since the equations for steady flow have been studied previously, the analysis centers on the third equation. This equation is solved numerically and by the method of matched asymptotic expansions. The exact numerical solutions compare well with the asymptotic solutions.The effects of suction and injection on pulsatile flow are the following: a) Small values of suction can cause a resonance-like effect for low frequency pulsatile flow. b) The annular effect still occurs but for large injection or suction the frequency at which this effect becomes dominant depends on the cross-flow Reynolds number. c) The maximum shear stress at the wall is decreased by injection, but may be increased or decreased by suction.Nomenclature a radius of the tube - a ₀ ² i ² - A₀, B₀, C₀, D₀, E₀ constant coefficients appearing in the expression for pressure - b a non-dimensionalized length - b ₀ ² i ^{2 2} - b _k complex coefficients of a power series - B - C ₁, C ₂, D complex constants - d - D _1,2 - f() F(a ^1/2)/aV - f ₀,f ₁,... functions of order one used in asymptotic expansions of f() - F(r) rv _r - g() - G(r) a steady component of velocity in axial direction - h() 4/C₀ a ² H(a ^1/2) - h ₀,h ₁,h ₂,...;l ₀,l ₁,l ₂,... functions of order one used in asymptotic expansions for h() in outer regions - H(r) complex valued function giving unsteady component of velocity - H ₀, H ₁, H ₂, ... K ₀, K ₁, K ₂, ...; L ₀, L ₁, L ₂, ... functions of order one used in asymptotic expansions for h() in inner regions - i - J ₀, J ₁, Y ₀, Y ₁ Bessel functions of first and second kind - k - K Rk/2b ² - O order symbol - p pressure - p ₁(z, t) arbitrary function related to pressure - r radial coordinate - r ₀ (1+16 ^{4 4})^1/4 - R Va/, the crossflow Reynolds number - t time - u() G(r)/V - v _r radial velocity - v _z axial velocity - V constant velocity at which fluid is injected or extracted - z axial coordinate - ² a ²/4 - 4.196 - small parameter; =–2/R (Sect. 4); =–R/2 (Sect. 5); =2/R(Sect. 6) - r ²/a ² - * 0.262 - Arctan (4 ^{2 2}) - , inner variables - kinematic viscosity - b - * zero of g() - density - (r, t) arbitrary function related to axial velocity - frequency     

Natural convection heat transfer through an enclosed horizontal layer of supercritical carbon dioxide

In natural convection heat transfer through a thin horizontal layer of carbon dioxide, maxima in the equivalent thermal conductivities are obtained in the vicinity of the respective pseudocritical temperatures at pressures of 75.8, 89.6 and 103.4 bar. The maxima are the more pronounced, the closer the critical point is approached.Comparison of experimental results with Nusselt equations shows good agreement except for the immediate vicinity of the pseudocritical temperature.In visual observations a distinct change in flow structure appears in the immediate vicinity of the pseudocritical temperature. A steady state polygon pattern and a boiling-like action could not be observed in this geometry.

---

Zusammenfassung Beim Wärmetransport durch freie Konvektion in einer dünnen waagerechten Schicht von Kohlendioxid ergaben sich Maxima der scheinbaren Wärmeleitfähigkeit in der Nähe der pseudokritischen Temperaturen bei Drükken von 75,8, 89,6 und 103,4 bar. Die Maxima sind um so ausgeprägter, je mehr man sich dem kritischen Punkt nähert.Ein Vergleich der Versuchsergebnisse mit Nusseltbeziehungen ergibt gute Übereinstimmung außer in unmittelbarer Umgebung der pseudokritischen Temperatur. Direkte Beobachtungen der Konvektionsmuster zeigen in unmittelbarer Umgebung der pseudokritischen Temperatur eine deutliche Strukturänderung. Ein stationäres Zellmuster und siedeähnliche Vorgänge konnten in dieser Anordnung nicht beobachtet werden.

---

Nomenclature A area of the heating or cooling plate - C constant in the correlation - g acceleration of gravity - h heat transfer coefficient - k thermal conductivity of fluid in the gap - k _e equivalent thermal conductivity - m, n exponents of dimensionless numbers - q heat flux - T _C,PC absolute temperature; critical C, pseudocritical PC - Gr Grashof numberg ( _h– _c) ³/ ² - Nu Nusselt numberh/k - Pr Prandtl number/ - thermal diffusivity - coefficient of volume expansion - width of gap - _c,h temperature of cooling (c)-, heating (h)-plate - _m arithmetic mean temperature ( _c+ _h)/2 - kinematic viscosity - _c,h fluid density at the temperature of the cooling (c)- or heating (h)-plate - heat flow rate through the gap     

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

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

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     

An experimental study of the lateral migration of a droplet in a creeping flow

The distribution of droplets in a plane Hagen-Poiseuille flow of dilute suspensions has been measured by a special LDA technique. This method assumes a well defined relation between the velocity of the droplets and their lateral position in the channel. The measurements have shown that the droplet distribution is non-uniform and depends on the viscosity ratio between the droplets and the carrier liquid. The results have been compared with a theory by Chan and Leal describing the lateral migration of suspended droplets.List of symbols a particle radius, m - d half width of the channel, m - Re flow Reynolds number, = 2 _m · d · /µ - flow velocity, m/s - _m flow velocity at the channel axis, m/s - We Weber number, = 2 _m ^Emphasis>/2 · d · / - x distance from center line (x = 0) of the channel, m - non-dimensional distance from the channel center line, _x ^d - y distance along the channel (y = 0 at channel inlet), m - non-dimensional distance along the channel, = y/2d - non-dimensional, normalized distance along the channel, = · _m · µ/ - interfacial tension, N/m - viscosity ratio of dispersed (droplet) phase to viscosity of continuous phase - µ viscosity of continuous phase, Pa · s - density of continuous phase, kg/m³ - phase density difference, kg/m³ Experiments were performed at Max-Planck-Institut, Göttingen     

The mirror relations and nonlinear viscoelasticity of polymer melts

An attempt is made to incorporate into a quasilinear viscoelastic constitutive equation of the Boltzmann superposition type the two mirror relations of Gleissle, as well as his relation between the steady-state first normal-stress difference and the shear viscosity curve. It is shown that the three relations can hold separately within this constitutive model, but not simultaneously, because they require a different nonlinear strain measure, namelyS ₁₂ () = – a ( – 1) (a = 0 for 1,a = 1 for 1) for the mirroring of the viscosities,S ₁₂ () = – a (–k ²/) (a = 0 for k, a = 1 for k) for the mirroring of the first normal-stress coefficients, and for the third relation. Here denotes the shear strain and erf the error function. Experimental data on melts of a low-density polyethylene, a high-density polyethylene and a polypropylene show that the mirror relations are passable approximations, but that the third relation meets reality surprisingly close if the right value ofk is used.     

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     

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     

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     

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     

Helical flow of molten polymers in a cylindrical annulus

Summary The paper is concerned with an analytical investigation of helical flow of a non-Newtonian fluid through an annulus with a rotating inner cylinder. The shear dependence of viscosity is described by a power law and the temperature dependence by an exponential function.Velocity and temperature profiles, energy input and shear along the stream lines, pressure drop, and torque are presented for the range of input parameters encountered in polymer extrusion.The results of the study can be applied to a mixing element in a screw extruder and for a device to control extrudate temperature and output.Nomenclature a thermal diffusivity [m²/s] - b temperature coefficient [K^–1], see eq. [4] - c heat capacity [J/kg K] - h slot width [m] - I ₁,I ₂,I ₃ invariants of the rate of deformation tensor, see eq. [5] - k thermal conductivity [J/m s K] - l, L = 1/h length of the slot - l _T ,l _K thermal and kinematic entrance length - m power law exponent, see eq. [3] - M torque [m N] - p pressure [N/m²] - P dimensionless pressure gradient, see eq. [24] - P _R,P _RZ dimensionless components of the shear stress tensor, see eq. [25] and eq. [26] - r, R = r/r _wa radial coordinate - r _wa, r_wi outer and inner radius of annulus [m] - t time [s]; dwell time in the annulus - T temperature [K] - v , v_r, V_z velocity components [m/s] - v ₀ angular velocity at inner wall [m/s] - average velocity inz-direction [m/s] - V , V_R, V_Z dimensionless velocity components,v /v₀, v_r/v₀, v_z/v₀ - V _z velocity ratio, helical parameter - Y coordinate inr-direction, see eq. [20] - z, Z = z/h Pe axial coordinate - deformation - rate of deformation tensor [s^–1] - apparent viscosity [N s/m²], see eq. [3] - dimensionless temperature,b (T – T ₀) - azimuth coordinate - ratio of radii,r _wi/r_wa - density [kg/m³] - , kl shear stress tensor [N/m²] - fluidity [m²w/Nw s], see eq. [4] - Gf Griffith number, see eq. [12] - Pe Péclet number, see eq. [13] - Re Reynolds number, - 0 initial state, reference state - equilibrium state - e entrance - wi, wa at surface of inner or outer wall - r, R, z, Z, coordinates - i, j radial and axial position of nodal point in the grid - k, l tensor components Presented at Euromech 37, Napoli 6. 20–23. 1972.With 15 figuresDedicated to Prof. Dr.-Ing. G. Schenkel on his 60th birthday     

An analytical solution for the temperature profiles during injection molding,including dissipation effects

An analytical solution is obtained for the stationary temperature profile in a polymeric melt flowing into a cold cavity, which also takes into account viscous heating effects. The solution is valid for the injection stage of the molding process. Although the analytical solution is only possible after making several (at first sight) rather stringent assumptions, the calculated temperature field turns out to give a fair agreement with a numerical, more realistic approach. Approximate functions were derived for both the dissipation-independent and the dissipation-dependent parts which greatly facilitate the temperature calculations. In particular, a closed-form expression is derived for the position where the maximum temperature occurs and for the thickness of the solidified layer.The expression for the temperature field is a special case of the solution of the diffusion equation with variable coefficients and a source term.Nomenclature a thermal diffusivity [m²/s] - c specific heat [J/kg K] - D channel half-height [m] - L channel length [m] - m 1/ - P pressure [Pa] - T temperature [°C] - T _W wall temperature [°C] - T _i injection temperature [°C] - T _A Br independent part of T - T _B Br dependent part of T - T ^core asymptotic temperature - v _z() axial velocity [m/s] - W channel width [m] - x cross-channel direction [m] - z axial coordinate [m] - (x) gamma function - (a, x) incomplete gamma function - M(a, b, x) Kummer function - small parameter - () temperature function - thermal conductivity [W/mK] - viscosity [Pa · s] - ₀ consistency index - power-law exponent - density [kg/m] - similarity variable Dimensionless variables Br Brinkman number - Gz Graetz number -      

Interaction of an oscillating vortex with a turbulent boundary layer

An oscillating vortex embedded within a turbulent boundary layer was generated experimentally by forcing a periodic lateral translation of a half-delta wing vortex generator. The objective of the experiment was to investigate the possibility that a natural oscillation, or meander, might be responsible for flattened vortex cores observed in previous work, which could also have contaminated previous turbulence measurements. The effect of this forced oscillation was characterized by comparison of measurements of the mean velocities and Reynolds stresses at two streamwise stations, for cases with and without forcing. The Reynolds stresses, especially w, were affected significantly by the forced oscillation, mainly through contributions from the individual production terms, provided the vortex was not too diffuse.List of Symbols a amplitude of forced vortex motion - f frequency of forced vortex generator motion - l vortex generator root chord - L flow length scale - R _Y , R _Z vortex core radial dimensions in vertical and spanwise directions, respectively - R_r vortex circulation Reynolds number R = / - u, v, w instantaneous velocity components in X, Y, Z directions - U, V, W mean velocities; shorthand notation for u, , w - X, Y, Z right-hand Cartesian streamwise, vertical, and spanwise coordinate directions - boundary-layer thickness - overall circulation - air kinematic viscosity - _x streamwise vorticity, _X = W/Y–V/d+t6Z - ( )₀ reference value (measured at X = 10 cm) - ( )c refers to vortex center - ( ) _max maximum value for a particular crossflow plane - ( ) (overbar) time average - ( ) (prime) fluctuating component, e.g., u=U+u     

An analogue study for flame flickering

An analogue experiment is proposed to simulate flame flickering comprising a free ascending column fed on its side with a light gas (helium) emerging from a vertical slot in ambient air. The convective motion of the helium jet is considered to represent the motion of burnt gases of buoyant jet flames. The helium jet is accelerated by buoyancy effects and the flow field is similar to that of burnt gases observed for real buoyant flames. The vertical velocity profile of the steady helium jet is measured at different vertical distances. The unsteady helium jet is also studied by measuring the instability frequency as a function of ambient pressure at different injection flow rates, and by analyzing the tomography images of the helium jet. The instability morphology is the same as that observed on real buoyant flames. We conclude that this type of instability can be approximately characterized by the maximum vertical velocityu _max, and the distance betweenu _max in the helium ascending column andu = o in the ambient air. For this type of instability the local vorticity is proportional to which can be influenced by gravity and ambient pressure. Theoretical prediction of the instability frequency as a function of gravity and ambient pressure has been obtained, and is in good agreement with the experimental results.List of symbols C ₁,C ₂ constants - F instability frequency - F _c critical frequency - F _m the most amplified frequency - F (K, ) function defined in (11) - g gravitational acceleration - g reduced gravity acceleration g(₀-^*)/^* - k real wave number of the disturbance - K reduced wave numberK=2k - K _c reduced wave number of the critical instability mode - K _m nondimensional wavenumber of the most amplified mode - L vertical characteristic length (in x direction) - P ambient pressure - u local vertical buoyant velocity (inx direction) - u _max local maximum vertical velocity - v local velocity component iny direction (horizontal) - V ₀ injection velocity of helium (iny direction) - x vertical distance measured from the leading edge of boundary layer - y horizontal distance measured from the exit plane of the vertical slot - Z(K, ) function defined in equation (11) Greek symbols distance betweenu _max in the helium ascending column andu = o in the ambient air - - wavelength of instability - _c critical wavelength - _m the most amplified wavelength - ^* helium density at slot exit - ₀ ambient air density - ^* helium dynamic viscosity at slot exit - v ^* helium kinematic viscosity at slot exit - complex number presented in disturbancee ^i(kx+t) - _i imaginary part of , representing the amplification rate of disturbance - _r real part of , where ( _r /k) represents the group velocity - reduced complex number of , defined      

Experimental investigation of the creeping motion of a drop in a vertical tube

New experimental data regarding the motion of a drop along the axis of a vertical tube, filled with another highly viscous liquid, are obtained. The experiments are realised with sufficiently large drops for an internal circulation to develop and also for different pairs of fluids; the preponderant role of the gravity on the drop shape and consequently on its terminal velocity is pointed out. Moreover, by means of a visualization technique, details on the flow both inside and outside the drop are given.List of symbols g gravity acceleration - r distance from the drop center - R equivalent radius of the drop, i.e. the radius of the sphere having the same volume as the drop - R _EQ radius of the equatorial section of the drop - R _T tube radius - L _AX half length of the drop - U ₀ terminal velocity of the drop - P _s Poiseuille number= U _{0 e} /4 g R ² - Fr Fronde number = U ₀ ² _e /2 g R - Re Reynolds number = 2 U ₀ R _e / _e - E _o Eötvös number = 4g R ²/ - deformation parameter = _e U ₀/ - apparent density of the suspended liquid= | _i – _e | - _i viscosity of the suspended liquid - _e viscosity of the suspending liquid - drop-to-tube radius ratio = R/R _T - _EQ equatorial drop-to-tube radius ratio = R EQ/R _T - interfacial tension     

Two-phase flow in heterogeneous porous media III: Laboratory experiments for flow parallel to a stratified system

Two-phase flow in stratified porous media is a problem of central importance in the study of oil recovery processes. In general, these flows are parallel to the stratifications, and it is this type of flow that we have investigated experimentally and theoretically in this study. The experiments were performed with a two-layer model of a stratified porous medium. The individual strata were composed of Aerolith-10, an artificial: sintered porous medium, and Berea sandstone, a natural porous medium reputed to be relatively homogeneous. Waterflooding experiments were performed in which the saturation field was measured by gamma-ray absorption. Data were obtained at 150 points distributed evenly over a flow domain of 0.1 × 0.6 m. The slabs of Aerolith-10 and Berea sandstone were of equal thickness, i.e. 5 centimeters thick. An intensive experimental study was carried out in order to accurately characterize the individual strata; however, this effort was hampered by both local heterogeneities and large-scale heterogeneities.The theoretical analysis of the waterflooding experiments was based on the method of large-scale averaging and the large-scale closure problem. The latter provides a precise method of discussing the crossflow phenomena, and it illustrates exactly how the crossflow influences the theoretical prediction of the large-scale permeability tensor. The theoretical analysis was restricted to the quasi-static theory of Quintard and Whitaker (1988), however, the dynamic effects described in Part I (Quintard and Whitaker 1990a) are discussed in terms of their influence on the crossflow.Roman Letters A interfacial area between the -region and the -region contained within V, m² - a vector that maps onto , m - b vector that maps onto , m - b vector that maps onto , m - B second order tensor that maps onto , m² - C second order tensor that maps onto , m² - E energy of the gamma emitter, keV - f fractional flow of the -phase - g gravitational vector, m/s² - h characteristic length of the large-scale averaging volume, m - H height of the stratified porous medium , m - i unit base vector in the x-direction - K local volume-averaged single-phase permeability, m² - K - {K}, large-scale spatial deviation permeability - { K} large-scale volume-averaged single-phase permeability, m² - K ^* large-scale single-phase permeability, m² - K ^** equivalent large-scale single-phase permeability, m² - K local volume-averaged -phase permeability in the -region, m² - K local volume-averaged -phase permeability in the -region, m² - K - {K } , large-scale spatial deviation for the -phase permeability, m² - K ^* large-scale permeability for the -phase, m² - l thickness of the porous medium, m - l characteristic length for the -region, m - l characteristic length for the -region, m - L length of the experimental porous medium, m - characteristic length for large-scale averaged quantities, m - n outward unit normal vector for the -region - n outward unit normal vector for the -region - n unit normal vector pointing from the -region toward the -region (n = - n ) - N number of photons - p pressure in the -phase, N/m² - p ₀ reference pressure in the -phase, N/m² - local volume-averaged intrinsic phase average pressure in the -phase, N/m² - large-scale volume-averaged pressure of the -phase, N/m² - large-scale intrinsic phase average pressure in the capillary region of the -phase, N/m² - - , large-scale spatial deviation for the -phase pressure, N/m² - p_c , capillary pressure, N/m² - p _c capillary pressure in the -region, N/m² - p capillary pressure in the -region, N/m² - {p _c } ^c large-scale capillary pressure, N/m² - q -phase velocity at the entrance of the porous medium, m/s - q -phase velocity at the entrance of the porous medium, m/s - S_wi irreducible water saturation - S /, local volume-averaged saturation for the -phase - S _i initial saturation for the -phase - S _r residual saturation for the -phase - S ^* { }^*/}^*, large-scale average saturation for the -phase - S saturation for the -phase in the -region - S saturation for the -phase in the -region - t time, s - v -phase velocity vector, m/s - v local volume-averaged phase average velocity for the -phase, m/s - {v } large-scale averaged velocity for the -phase, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - V large-scale averaging volume, m³ - y position vector relative to the centroid of the large-scale averaging volume, m - {y}^c large-scale average of y over the capillary region, m Greek Letters local porosity - local porosity in the -region - local porosity in the -region - local volume fraction for the -phase - local volume fraction for the -phase in the -region - local volume fraction for the -phase in the -region - {}^* { }^*+{ }^*, large-scale spatial average volume fraction - { }^* large-scale spatial average volume fraction for the -phase - mass density of the -phase, kg/m³ - mass density of the -phase, kg/m³ - viscosity of the -phase, N s/m² - viscosity of the -phase, Ns/m² - V /V , volume fraction of the -region ( + =1) - V /V , volume fraction of the -region ( + =1) - attenuation coefficient to gamma-rays, m^-1 - -      

The very slow flow of a Powell-Eyring fluid around a sphere

Summary The very slow flow of a Powell-Eyring type non-Newtonian fluid around a sphere is investigated by a variational technique. The result, a correction factor that is applied to the Stokes' equation, is given as a plot and as an equation which is empirically fit to the plot. Also, a comparison of the very slow flows of a simplified viscoelastic Oldroyd fluid and the Powell-Eyring fluid is made which indicates that in a certain restricted region of the very slow flows, both models give essentially the same results. The Oldroyd and Powell-Eyring model parameters are interrelated by forcing both models to fit the same tube flow viscosity data.Nomenclature B dimensionless quantity, v /R - C dimensionless second invariant - c ₁ constant determined by variational method - D dimensionless variational integral - D _2j , D _j+k position-independent variables used in specification of trial functions - E _2j , E _j+k position-independent variables used in specification of trial functions - f friction factor - f _corr friction factor correction - F _drag drag force on sphere - g, g ₀, g ₁ general trial function; first and second terms in the general trial function - G, H terms in the expression for C - j index - J variational integral - k index - K term in the expression for C - p, q integers - r integer, radial coordinate - R radius of sphere - Re Reynolds number - Re ₀ Reynolds number at point of zero shear rate - Re Reynolds number at infinite distance from sphere - Re _NN Reynolds number based on variable part of viscosity - u, v dimensionless position coordinates - V volume considered - v _i ith velocity component - v _r , v , v _z velocity components in the r, , and z-directions - v approach velocity of the fluid - x/ parameter in Powell-Eyring model - x _i i-position coordinate - parameter in Powell-Eyring model - rate of deformation - , _c , _N , ₀ coefficient of viscosity; cross viscosity; parameter in Powell-Eyring model; viscosity in limit of zero shear rate - spherical coordinate - , _ij rate of deformation tensor; ij-component of rate of deformation tensor - ₁, ₂ parameters in Oldroyd model - Newtonian viscosity - ₁, ₂ parameters in Oldroyd model - dimensionless radial coordinate, r/R - second invariant - fluid density - spherical coordinate - stream function     

Peristaltic pumping of a second-order fluid in a planar channel

The peristaltic motion of a non-Newtonian fluid represented by the constitutive equation for a second-order fluid was studied for the case of a planar channel with harmonically undulating extensible walls. A perturbation series for the parameter ( half-width of channel/wave length) obtained explicit terms of 0(²), 0(²Re²) and 0(₁Re²) respectively representing curvature, inertia and the non-Newtonian character of the fluid. Numerical computations were performed and compared to the perturbation analysis in order to determine the range of validity of the terms.Presented at the second conference Recent Developments in Structured Continua, May 23–25, 1990, in Sherbrooke, Québec, Canada