Reynolds number dependence of the drag coefficient for laminar flow through electrically-heated photoetched screens

Experiments are performed to measure the drag coefficient of electrically-heated screens. Square-pattern 80 mesh and 100 mesh screens of 50.8 m-wide wires photoetched from 50.8 m thick Inconel sheets are examined. Ambient air is passed through these screens at upstream velocities yielding wire-width Reynolds numbers from 2 to 35, and electrical current is passed through the screens to generate heat fluxes from o to 0.17 MW/m², based on the total screen area. The dependence of the drag coefficient on Reynolds number and heat flux is determined for these two screens by measuring pressure drops across the screens for a variety of conditions in these ranges. In all cases, heating is found to increase the drag coefficient above the unheated value. A correlation relating the heated drag coefficient to the unheated drag coefficient is developed based on the idea that the main effect of heating at these levels is to modify the Reynolds number through modifying the viscosity. This correlation is seen to reproduce the experimental results closely.List of Symbols A total screen cross sectional area - C fitting coefficient, near unity - c _D heated drag coefficient - c _{D, 0} unheated drag coefficient - C _p air specific heat at constant pressure - D photoetched wire width, sheet thickness - h _s stagnation point heat-transfer coefficient - k air thermal conductivity - M distance between adjacent wires - O open area fraction - p air pressure - p air pressure drop across screen - Pr Prandtl number for air, c _p/k - Q total electrical power to screen - R radius of curvature at stagnation point - Re _D wire width Reynolds number, UD/ - T air temperature - U air speed upstream of screen - air specific heat ratio - air density - air viscosity - exponent in temperature power law for viscosity - () quantity () evaluated at heated screen temperature The authors thank John Lewin and Bob Meyer for their assistance in the design and fabrication of the heated screen test facility and Tom Grasser for his help in performing the experiments. This work was performed at Sandia National laboratories, supported by the U.S. Department of Energy under contract number DE-AC04-94AL85000.     

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].     

Finite difference analysis of plane Poiseuille and Couette flow developments

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

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

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

Nonlinear analysis of the forced response of a beam with three mode interaction

An analysis is presented for the primary resonance of a clamped-hinged beam, which occurs when the frequency of excitation is near one of the natural frequencies,_n . Three mode interaction (₂ 3₁ and _{3 1} + 2₂) is considered and its influence on the response is studied. The case of two mode interaction (₂ 3₁) is also considered to compare it with the case of three mode interaction. The straight beam experiencing mid-plane stretching is governed by a nonlinear partial differential equation. By using Galerkin's method the governing equation is reduced to a system of nonautonomous ordinary differential equations. The method of multiple scales is applied to solve the system. Steady-state responses and their stability are examined. Results of numerical investigations show that there exists no significant difference between both modal interactions' influences on the responses.     

Chaos and integrability in a nonlinear wave equation

We consider the parametrized family of equations _tt ,u- _xx u-au+u ₂ ² u=O,x(0,L), with Dirichlet boundary conditions. This equation has finite-dimensional invariant manifolds of solutions. Studying the reduced equation to a four-dimensional manifold, we prove the existence of transversal homoclinic orbits to periodic solutions and of invariant sets with chaotic dynamics, provided that =2, 3, 4,.... For =1 we prove the existence of infinitely many first integrals pairwise in involution.     

The value of the critical exponent for reaction-diffusion equations in cones

Let D R ^N be a cone with vertex at the origin i.e., D = (0, )x where S ^N–1 and x D if and only if x = (r, ) with r=¦x¦, . We consider the initial boundary value problem: u _t = u+u ^p in D×(0, T), u=0 on Dx(0, T) with u(x, 0)=u ₀(x) 0. Let ₁ denote the smallest Dirichlet eigenvalue for the Laplace-Beltrami operator on and let ₊ denote the positive root of (+N–2) = ₁. Let p ^* = 1 + 2/(N + ₊). If 1 < p < p ^*, no positive global solution exists. If p>p ^*, positive global solutions do exist. Extensions are given to the same problem for u _t=+¦x¦ u ^p .This research was supported in part by the Air Force Office of Scientific Research under Grant # AFOSR 88-0031 and in part by NSF Grant DMS-8 822 788. The United States Government is authorized to reproduce and distribute reprints for governmental purposes not withstanding any copyright notation therein.     

Aspect ratio effect on flow-induced forces on circular cylinders in a cross-flow

Finite-span circular cylinders with two different aspect ratios, placed in a cross-flow, are investigated experimentally at a cylinder Reynolds number of 46,000. Simultaneous measurements of the flow-induced unsteady forces on the cylinders and the stream velocity in the wake are carried out. These results together with mean drag measurements along the span and available literature data are used to evaluate the flow mechanisms responsible for the induced unsteady forces and the effect of aspect ratio on these forces. The coherence of vortex shedding along the span of the cylinder is partially destroyed by the separated flow emanating from the top and by the recirculating flow behind the cylinder. As a result, the fluctuating lift decreases drastically. Based on the data collected, it is conjectured that the fluctuating recirculating flow behind the cylinder is the flow mechanism responsible for the unsteady drag and causes it to increase beyond the fluctuating lift. The fluctuating recirculating flow is a direct consequence of the unsteady separated flow. The unsteady forces vary along the span, with lift increasing and drag decreasing towards the cylinder base. When the cylinder span is large compared to the wall boundary layer thickness, a submerged two-dimensional region exists near the base. As the span decreases, the submerged two-dimensional region becomes smaller and eventually vanishes. Altogether, these results show that fluctuating drag is the dominant unsteady force in finite-span cylinders placed in a cross-flow. Its characteristic frequency is larger than that of the vortex shedding frequency.List of symbols a span of active element on cylinder, = 2.5 cm - C _D local rms drag coefficient, 2D/ U ² da - C _L local rms lift coefficient, 2l/ U ² da - C _D local mean drag coefficient, 2D/ U ² da - C _D spanwise-averaged C _D for finite-span cylinder - (C _D ) _2D spanwise-averaged mean drag coefficient for two-dimensional cylinder - C _p pressured coefficient - -(C _p ) _b pressure coefficient at = - d diameter of cylinder, = 10.2 cm - D fluctuating component of instantaneous drag - D local rms of fluctuating drag - D local mean drag - E _D power spectrum of fluctuating drag, defined as - E _L power spectra of fluctuating lift, defined as - f _D dominant frequency of drag spectrum - f _L dominant frequency of lift spectrum - f _u dominant frequency of velocity spectrum - h span of cylinder - H height of test section, = 30.5 cm - L fluctuating component of instantaneous lift - L local rms of fluctuating lift - R _Du () cross-correlation function of streamwise velocity and local drag, - R _Lu () cross-correlation function of stream wise velocity and local lift, - Re Reynolds number, U d/y - S _L Strouhal number based on f _L ,f _L d/U - S _D Strouhal number based on f _D ,f _D d/U - S _u Strouhal number based on f _u , f _u d/U - t time - u fluctuating component of instantaneous streamwise velocity - U mean streamwise velocity - mean stream velocity upstream of cylinder - x streamwise distance measured from axis of cylinder - y transverse distance measured from axis of test section - z spanwise distance measured from cylinder base - angular position on cylinder circumference measured from forward stagnation - kinematic viscosity of air - density of air - time lag in cross-correlation function - _D normalized spectrum of fluctuating drag - _L normalized spectrum of fluctuating lift     

New concepts in relative permeabilities at high capillary numbers for surfactant flooding

Flooding oil reservoirs with surfactant solutions can increase the amount of oil that can be recovered. Macroscopic modelling of the process requires relative permeabilities to be functions of saturation and capillary number. With only limited experimental data, relative permeabilities have usually been assumed to be linear functions of saturation at high capillary numbers. The experimental data is reviewed, some of which suggest that this assumption is not necessarily correct. The basis for the assumption is therefore reviewed and it is concluded that the linear model corresponds to microscopically segregated flow in the porous medium. Based on new but equally plausible complementary assumptions about the flow pattern, a mixed flow model is derived. These models are then shown to be limiting cases of a droplet model which represents the mixing scale within the porous medium and gives a physical basis for interpolating between the models. The models are based on physical concepts of flow in a porous medium and so the approach described here represents a significant improvement in the understanding of high capillary number flow. This is shown by the fact that fewer parameters are needed to describe experimental data.Notation A total cross-sectional area assigned to capillary bundle - A ⁽ⁱ⁾ physical cross-sectional area of tube i - c ⁽ⁱ⁾ ordered configurational label for droplets in tube i - c configuration label for tube i (order not considered) - D defined by Equation (26) - E(...) expectation value with respect to the trinomial distribution - S _r ⁽⁾ fractional flow of phase - k absolute permeability - k _r relative permeability of phase - k _r ⁰ endpoint relative permeability of phase - L capillary tube length in bundle model - m ⁽ⁱ⁾ number of droplets of phase a occupying tube i - n exponent for phase a in Equation (2) - N number of droplets in bundle model - N _c capillary number - p pressure - p(c') probability of configuration c - Q ⁽ⁱ⁾ total volume flow rate in tube i - S saturation of phase - S flowing saturation of phase - S _r residual saturation of phase - S _r ⁽⁾ saturations when fractional flow of phase is 1 in the case of varying residual saturations for three-phase flow ( ) - t _c residence time for droplet configuration c - v ⁽ⁱ⁾ total fluid velocity in bundle tube i - , phase label - p pressure differential across capillary bundle - ⁽ⁱ⁾ tube conductivity defined by Equation (7) - viscosity of phase - interfacial tension - gradient operator - ... average over tube droplet configurations     

A power-law fluid flow past a porous sphere

Summary A model has been developed for the flow of a non-Newtonian fluid past a porous sphere. The drag force exerted on a porous sphere moving in a power-law fluid is obtained by an approximate solution of equations of motion in the creeping flow regime. It is predicted that the effect of the pseudoplastic anomaly on the drag force is more pronounced at large porosity parameters.

---

Zusammenfassung Es wird ein Modell für die Strömung einer nichtnewtonschen Flüssigkeit längs einer porösen Kugel entwickelt. Die auf die in einer Ostwald-DeWaele-Flüssigkeit bewegte Kugel ausgeübte Reibungskraft wird durch eine Näherungslösung der Bewegungsgleichungen für schleichende Strömung gewonnen. Man findet, daß der Einfluß der Abweichung vom newtonschen Verhalten um so ausgeprägter wird, je größer die Porosität ist.

---

A, B, C, D a, b, c, d coefficients in eqs. [10] and [18] - F _D drag force - K consistency index in power-law model - k ₁ ,k ₂ coefficients defined by eq. [18] - m porosity parameter - n flow index in power-law model - P pressure - P ^* dimensionless pressure defined by eq. [4] - P pressure difference - R radius of porous sphere - r radial distance from the center of the sphere - U velocity of uniform stream - u _i velocity component - u _i ^* dimensionless velocity component defined by eq. [4] - Y drag force correction factor defined by eq. [27] - _ij rate of deformation tensor - _ij ^* dimensionless rate of deformation tensor defined by eq. [4] - , spherical coordinates - dimensionless radial distance defined by eq. [4] - second invariant of rate of deformation tensor - ^* dimensionless second invariant of rate of deformation tensor defined by eq. [4] - _ij stress tensor - _ij ^* dimensionless stress tensor defined by eq. [4] - stream function - ^* dimensionless stream function defined by eq. [4] - i inside the surface of the sphere - o outside the surface of the sphere With 1 figure and 1 table     

Shear viscosity behavior of highly concentrated suspensions at low and high shear-rates

A method is proposed for determining the shear viscosity behavior of highly concentrated suspensions at low and high shear-rates through the use of a formulation that is a function of three parameters signifying the effects of particle size distribution. These parameters are the intrinsic viscosity [], a parametern that reflects the level of particle association at the initiation of motion and the maximum packing concentration _m. The formulation reduces to the modified Eilers equation withn = 2 for high shear rates. An analytical method was used for the calculation of maximum packing concentration which was subsequently correlated with the experimental values to account for the surface induced interaction of particles with the fluid. The calculated values of viscosities at low and high shear-rates were found to be in good agreement with various experimental data reported in literature. A brief discussion is also offered on the reliability of the methods of measuring the maximum packing concentration. _r = /₀ relative viscosity of the suspension - volumetric concentration of solids - k _n coefficient which characterizes a specific effect of particle interactions - _m maximum packing concentration - _r,0 relative viscosity at low shear-rates - [] intrinsic viscosity - n, n parameter that reflects the level of particle interactions at low and high shear-rates, respectively - _r, relative viscosity at high shear-rates - (_m)_s, (_m)_i, (_m)_l packing factors for small, intermediate and large diameter classes - v _s, v_i, v_l volume fractions of small, intermediate and large diameter classes, respectively - _si, _sl coefficient to be used in relating a smaller to an intermediate and larger particle group, respectively - _is, _il coefficient to be used in relating an intermediate to a smaller and larger particle group, respectively - _ls, _li coefficient to be used in relating a larger to a smaller and intermediate particle group, respectively - _m0 maximum packing concentration for binary mixtures - _m,e measured maximum packing concentration - _m,c calculated maximum packing concentration     

On the boundary conditions at the macroscopic level

We study the problem of the boundary conditions specified at the boundary of a porous domain in order to solve the macroscopic transfer equations obtained by means of the volume-averaging method. The analysis is limited to the case of conductive transport but the method can be extended to other cases. A numerical study enables us to illustrate the theoretical results in the case of a model porous medium. Roman Letters _sf interfacial area of the s-f interface contained within the macroscopic system m² - A _sf interfacial area of the s-f interface contained within the averaging volume m² - C _p mass fraction weighted heat capacity, kcal/kg/K - d _s , d _f microscopic characteristic length m - g vector that maps to _s, m - h vector that maps to _f , m - K _eff effective thermal conductivity tensor, kcal/m s K - l REV characteristic length, m - L macroscopic characteristic length, m - n _fs outwardly directed unit normal vector for the f-phase at the f-s interface - n _e outwardly directed unit normal vector at the dividing surface - T ^* macroscopic temperature field obtained by solving the macroscopic equation (3), K - V averaging volume, m³ - V _s , V _f volume of the considered phase within the averaging volume, m³ - volume of the macroscopic system, m³ - _s , _f volume of the considered phase within the volume of the macroscopic system, m³ - dividing surface, m² Greek Letters _s , _f volume fraction - ratio of thermal conductivities - _s , _f thermal conductivities, kcal/m s K - spatial average density, kg/m³ - microscopic temperature, K - ^* microscopic temperature corresponding to T ^* , K - spatial deviation temperature K - error on the temperature due to the macroscopic boundary conditions, K - spatial average - ^s , ^f intrinsic phase average     

The theory of a vibrating-rod viscometer

The paper presents a complete theory for a viscometer based upon the principle of a circular-section rod, immersed in a fluid, performing transverse oscillations perpendicular to its axis. The theory is established as a result of a detailed analysis of the fluid flow around the rod and is subject to a number of criteria which subsequently constrain the design of an instrument. Using water as an example it is shown that a practical instrument can be designed so as to enable viscosity measurement with an accuracy of ±0.1%, although it is noted that many earlier instruments failed to satisfy one or more of the newly-established constraints.Nomenclature A, D constants in equation (46) - A _m , B _m , C _m , D _m constants in equations (50) and (51) - A _j , B _j constants in equation (14) - a _j ⁺ , a _j ^– wavenumbers given by equation (15) - C _f drag coefficient defined in equation (53) - c speed of sound - D _b drag force of fluid b - D ₀ coefficient of internal damping - E extensional modulus - f(z) initial deformation of rod - f(), F _m () functions of defined in equation (41) - F force in the rod - force per unit length near t=0 - F dimensionless force per unit length near t=0 - g _m amplitude of transient force - G modulus of rigidity - h, h* functions defined by equations (71) and (72) - H functions defined by equation (69) and (70) - I second moment of area - I _0,1, J _0,1, K _0,1 modified Bessel functions - k, k functions defined in equations (2) - L half-length of oscillator - Ma Mach number - m _b added mass per unit length of fluid b - m _s mass per unit length of solid - n _j eigenvalue defined in equations (15) and (16) - R radius of rod - R _c radius of container - r radial coordinate - T tension - T _visc temperature rise due to heat generation by viscous dissipation - t time - v _r , v radial and angular velocity components - y lateral displacement - y ₀ initial lateral displacement - y ₁, y ₂ successive maximum lateral displacement - z axial coordinate - dimensionless tension - dimensionless mass of fluid - dimensionless drag of fluid - amplification factor - logarithmic decrement in a fluid - _a , _b logarithmic decrement in fluids a and b - ₀ logarithmic decrement in vacuo - _j logarithmic decrement in mode j in a fluid - spatial resolution of amplitude - _v voltage resolution - r, , , _s, , increments in R, , , _s , , - dimensionless amplitude of oscillation - dimensionless axial coordinate - angular coordinate - _f thermal conductivity of fluid - viscosity of fluid - viscosity of fluid calculated on assumption that * - _a , _b viscosity of fluids a and b - _m constants in equation (10) - dimensionless displacement - _j j the component of - density of fluid - _a , _b density of fluids a and b - _s density of tube or rod material - dimensionless radial coordinate - * dimensionless radius of container - dimensionless times - spatial component of defined in equation (11) - _j , _tm jth, mth component of - dimensionless streamfunction - ⁰, ¹ components of in series expansion in powers of - streamfunction - dimensionless frequency (based on ) - angular frequency - ₀ angular frequency in absence of fluid and internal damping - _j angular frequency in mode j in a fluid - _a , _b frequencies in fluids a and b     

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

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

---

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

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

---

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

Heat and mass transfer in laminar flow between parallel porous plates

Summary The problem of heat transfer in a two-dimensional porous channel has been discussed by Terrill [6] for small suction at the walls. In [6] the heat transfer problem of a discontinuous change in wall temperature was solved. In the present paper the solution of Terrill for small suction at the walls is revised and the whole problem is extended to the cases of large suction and large injection at the walls. It is found that, for all values of the Reynolds number R, the limiting Nusselt number Nu increases with increasing R.Nomenclature stream function - 2h channel width - x, y distances measured parallel and perpendicular to the channel walls respectively - U velocity of fluid at x=0 - V constant velocity of fluid at the wall - =y/h nondimensional distance perpendicular to the channel walls - f() function defined in equation (1) - coefficient of kinematic viscosity - R=Vh/ suction Reynolds number - density of fluid - C _p specific heat at constant pressure - K thermal conductivity - T temperature - x=x ₀ position where temperature of walls changes - T ₀, T ₁ temperature of walls for x<x ₀, x>x ₀ respectively - = (T – T ₁)/T ₀ – T ₁) nondimensional temperature - =x/h nondimensional distance along channel - R ^* = Uh/v channel Reynolds number - Pr = C _p/K Prandtl number - _n eigenvalues - B _n() eigenfunctions - B _n ⁽⁰⁾ , () eigenfunctions for R=0 - B ₀ ⁽ⁱ⁾ , B ₀ ⁽ⁱⁱ⁾ ... change in eigenfunctions when R0 and small - K _n constants given by equation (13) - h heat transfer coefficient - Nu Nusselt number - _m mean temperature - C _n constants given by equation (18) - perturbation parameter - B _0i () perturbation approximations to B ₀() - Q = B ₀/ ₀ derivative of eigenfunction with respect to eigenvalue - z nondimensional distance perpendicular to the channel walls - F(z) function defined by (54)     

Eine analytische Näherungslösung für die eingefrorene laminare Grenzschichtströmung eines binären Gemisches

Zusammenfassung Für die eingefrorene laminare Grenzschichtströmung eines teilweise dissoziierten binären Gemisches entlang einer stark gekühlten ebenen Platte wird eine analytische Näherungslösung angegeben. Danach läßt sich die Wandkonzentration als universelle Funktion der Damköhler-Zahl der Oberflächenreaktion angeben. Für das analytisch darstellbare Konzentrationsprofil stellt die Damköhler-Zahl den Formparameter dar. Die Wärmestromdichte an der Wand bestehend aus einem Wärmeleitungs- und einem Diffusionsanteil wird angegeben und diskutiert. Das Verhältnis beider Anteile läßt sich bei gegebenen Randbedingungen als Funktion der Damköhler-Zahl ausdrücken.

---

An analytical approximation for the frozen laminar boundary layer flow of a binary mixture

An analytical approximation is derived for the frozen laminar boundary layer flow of a partially dissociated binary mixture along a strongly cooled flat plate. The concentration at the wall is shown to be a universal function of the Damkohler-number for the wall reaction. The Damkohlernumber also serves as a parameter of shape for the concentration profile which is presented in analytical form. The heat transfer at the wall depending on a conduction and a diffusion flux is derived and discussed. The ratio of these fluxes is expressed as a function of the Damkohler-number if the boundary conditions are known.

---

Formelzeichen A Atom - A₂ Molekül - C Konstante in Gl. (20) - c₁=1/(2C) Konstante in Gl. (35) - c_p spezifische Wärme bei konstantem Druck - D binärer Diffusionskoeffizient - Ec=u ² /(2h_f) Eckert-Zahl - h spezifische Enthalpie - h_t=h+u²/2 totale spezifische Enthalpie - h _A ⁰ spezifische Dissoziationsenthalpie - K_w Reaktionsgeschwindigkeitskonstante der heterogenen Wandreaktion - 1= /( ) Champman-Rubesin-Parameter - Le=Pr/Sc Lewis-Zahl - M Molmasse - p statischer Druck - Pr= c_pf/ Prandtl-Zahl - q_w Wärmestromdichte an der Wand - q_cw, q_dw Wärmeleitungsbzw. Diffusionsanteil der Wärmestromdichte an der Wand - universelle Gaskonstante - R=/(2M_a) individuelle Gaskonstante der molekularen Komponente - Re_x= u x/ Reynolds-Zahl - Sc=/( D) Schmidt-Zahl - T absolute Temperatur - T_d=h _A ⁰ /R charakteristische Dissoziationstemperatur - u, v x- und y-Komponenten der Geschwindigkeit - U=u/u normierte x-Komponente der Geschwindigkeit - x, y Koordinaten parallel und senkrecht zur Platte Griechische Symbole - =A/ Dissoziationsgrad - Grenzschichtdicke - ₂ Impulsverlustdicke - Damköhler-Zahl der Oberflächenreaktion - =T/T normierte Temperatur - =y/ normierter Wandabstand - Wärmeleitfähigkeit - dynamische Viskosität - , ^* Ähnlichkeitskoordinaten - Dichte - Schubspannung Indizes A auf ein Atom bezogen - M auf ein Molekül bezogen - f auf den eingefrorenen Zustand bezogen - w auf die Wand bezogen - auf den Außenrand der Grenzschicht bezogen     

Flow visualization and analysis of aerodynamically detuned supersonic rotors

The effect of aerodynamic detuning on the supersonic steady and unsteady blade passage flow field is experimentally investigated on a free surface water table by means of color Schlieren and shadowgraph flow visualization techniques. Two aerodynamic detuning mechanisms are considered: (1) alternate circumferential spacing of adjacent airfoils; (2) the replacement of alternate airfoils with splitters. The steady flow visualization demonstrates the significant effect of aerodynamic detuning on the passage flow field and, in particular, the shock wave-airfoil surface intersection locations. The unsteady flow visualization studies show the importance of the interblade phase angle. A mathematical model is also described and utilized to demonstrate the enhanced aeroelastic stability associated with the altered cascade passage shock wave structure due to these aerodynamic detuning mechanisms.List of symbols a dimensionless perturbation sonic velocity - C airfoil chord - I _x mass moment of inertia - k reduced frequency, k = C/2 - K spring constant - M _R dimensionless unsteady aerodynamic moment - M Mach number - P _split split splitter circumferential spacing - P _start splitter leading edge location - R reference full chord airfoil - R _s reference splitter - Sp _chord splitter chord length - u dimensionless perturbation chordwise velocity - v dimensionless perturbation normal velocity - amplitude of oscillation - interblade phase angle - level of aerodynamic detuning - undamped natural torsional frequency - ₀ reference frequency - flutter frequency     

Translatory accelerating motion of a circular disk in a viscous fluid

The exact solution of the equation of motion of a circular disk accelerated along its axis of symmetry due to an arbitrarily applied force in an otherwise still, incompressible, viscous fluid of infinite extent is obtained. The fluid resistance considered in this paper is the Stokes-flow drag which consists of the added mass effect, steady state drag, and the effect of the history of the motion. The solutions for the velocity and displacement of the circular disk are presented in explicit forms for the cases of constant and impulsive forcing functions. The importance of the effect of the history of the motion is discussed.Nomenclature a radius of the circular disk - b one half of the thickness of the circular disk - C dimensionless form of C ₁ - C ₁ magnitude of the constant force - D fluid drag force - f(t) externally applied force - F() dimensionaless form of applied force - F ₀ initial value of F - g gravitational acceleration - H() Heaviside step function - k magnitude of impulsive force - K dimensionless form of k - M a dimensionless parameter equals to (1+37#x03C0;_s/4_f) - S displacement of disk - t time - t ₁ time of application of impulsive force - u velocity of the disk - V dimensionless velocity - V ₀ initial velocity of V - V _t terminal velocity - parameter in (13) - parameter in (13) - (t) Dirac delta function - ratio of b/a - () function given in (5) - dynamical viscosity of the fluid - kinematic viscosity of the fluid - _f fluid density - _s mass density of the circular disk - dimensionless time - _i dimensionless form of t _i - dummy variable - dummy variable     

The theory of a vibrating-rod densimeter

The paper presents a theory of a device for the accurate determination of the density of fluids over a wide range of thermodynamic states. The instrument is based upon the measurement of the characteristics of the resonance of a circular section tube, or rod, performing steady, transverse oscillations in the fluid. The theory developed accounts for the fluid motion external to the rod as well as the mechanical motion of the rod and is valid over a defined range of conditions. A complete set of working equations and corrections is obtained for the instrument which, together with the limits of the validity of the theory, prescribe the parameters of a practical design capable of high accuracy.Nomenclature A, B, C, D constants in equation (60) - A _j , B _j constants in equation (18) - a _j ⁺ , a _j ^– wavenumbers given by equation (19) - C _f drag coefficient defined in equation (64) - C _f ^/0 , C _f ^/1 components of C _f in series expansion in powers of - c speed of sound - D _b drag force of fluid b - D ₀ coefficient of internal damping - E extensional modulus - force per unit length - F _j ⁺ , F _j ^– constants in equation (24) - f, g functions of defined in equations (56) - G modulus of rigidity - I second moment of area - K constant in equation (90) - k, k constants defined in equations (9) - L half-length of oscillator - Ma Mach number - m _a mass per unit length of fluid a - m _b added mass per unit length of fluid b - m _s mass per unit length of solid - n _j eigenvalue defined in equation (17) - P power (energy per cycle) - P _a , P _b power in fluids a and b - p pressure - R radius of rod or outer radius of tube - R _c radius of container - R _i inner radius of tube - r radial coordinate - T tension - T _visc temperature rise due to heat generation by viscous dissipation - t time - v _r , v radial and angular velocity components - y lateral displacement - z axial coordinate - dimensionless tension - _a dimensionless mass of fluid a - _b dimensionless added mass of fluid b - _b dimensionless drag of fluid b - dimensionless parameter associated with - ₀ dimensionless coefficient of internal damping - dimensionless half-width of resonance curve - dimensionless frequency difference defined in equation (87) - spatial resolution of amplitude - R, , , _s , increments in R, , , _s , - dimensionless amplitude of oscillation - dimensionless axial coordinate - ratio of to - _a , _b ratios of to for fluids a and b - angular coordinate - parameter arising from distortion of initially plane cross-sections - _f thermal conductivity of fluid - dimensionless parameter associated with - viscosity of fluid - _a , _b viscosity of fluids a and b - dimensionless displacement - _j jth component of - density of fluid - _a , _b density of fluids a and b - _s density of tube or rod material - density of fluid calculated on assumption that * - dimensionless radial coordinate - * dimensionless radius of container - dimensionless times - _{rr rr}, _r radial normal and shear stress components - spatial component of defined in equation (13) - _j jth component of - dimensionless streamfunction - ⁰, ¹ components of in series expansion in powers of - phase angle - _r phase difference - _ra , _rb phase difference for fluids a and b - streamfunction - _j jth component defined in equation (22) - dimensionless frequency (based on ) - _a , _b dimensionless frequency in fluids a and b - _s dimensionless frequency (based on _s ) - angular frequency - ₀ resonant frequency in absence of fluid and internal damping - _r resonant frequency in absence of internal fluid - _ra , _rb resonant frequencies in fluids a and b - dimensionless frequency - dimensionless frequency when _a vanishes - dimensionless frequencies when _a vanishes in fluids a and b - dimensionless resonant frequency when _a , _b, _b and ₀ vanish - dimensionless resonant frequency when _a , _b and _b vanish - dimensionless resonant frequency when _b and _b vanish - dimensionless frequencies at which amplitude is half that at resonance     

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²