Transient Free Convection from a Horizontal Surface in a Porous Medium Subjected to a Sudden Change in Surface Heat Flux

This paper presents an analytical and numerical study of transient free convection from a horizontal surface that is embedded in a fluid-saturated porous medium. It is assumed that for time steady state velocity and temperature fields are obtained in the boundary-layer which occurs due to a uniform flux dissipation rate q ₁ on the surface. Then, at the heat flux on the surface is suddenly changed to q ₂ and maintained at this value for . Firstly, solutions which are valid for small and large are obtained. The full boundary-layer equations are then integrated step-by-step for the transient regime from the initial unsteady state ( ) until such times at which this forward marching approach is no longer well posed. Beyond this time no valid solutions could be obtained which matched the final solution from the forward integration to the steady state profiles at large times .     

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     

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)     

Inhomogeneous ‘longitudinal’ plane waves in a deformed elastic material

By definition, a homogeneous isotropic compressible Hadamard material has the property that an infinitesimal longitudinal homogeneous plane wave may propagate in every direction when the material is maintained in a state of arbitrary finite static homogeneous deformation. Here, as regards the wave, homogeneous means that the direction of propagation of the wave is parallel to the direction of eventual attenuation; and longitudinal means that the wave is linearly polarized in a direction parallel to the direction of propagation. In other words, the displacement is of the form u = ncos k(n · x – ct), where n is a real vector. It is seen that the Hadamard material is the most general one for which a longitudinal inhomogeneous plane wave may also propagate in any direction of a predeformed body. Here, inhomogeneous means that the wave is attenuated, in a direction distinct from the direction of propagation; and longitudinal means that the wave is elliptically polarized in the plane containing these two directions, and that the ellipse of polarization is similar and similarly situated to the ellipse for which the real and imaginary parts of the complex wave vector are conjugate semi-diameters. In other words, the displacement is of the form u = {S exp i(S · x – ct)}, where S is a complex vector (or bivector). Then a Generalized Hadamard material is introduced. It is the most general homogeneous isotropic compressible material which allows the propagation of infinitesimal longitudinal inhomogeneous plane circularly polarized waves for all choices of the isotropic directional bivector. Finally, the most general forms of response functions are found for homogeneously deformed isotropic elastic materials in which longitudinal inhomogeneous plane waves may propagate with a circular polarization in each of the two planes of central circular section of the ⁿ -ellipsoid, where is the left Cauchy-Green strain tensor corresponding to the primary pure homogeneous deformation.     

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     

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     

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.     

On the closure problem for Darcy's law

In a previous derivation of Darcy's law, the closure problem was presented in terms of an integro-differential equation for a second-order tensor. In this paper, we show that the closure problem can be transformed to a set of Stokes-like equations and we compare solutions of these equations with experimental data. The computational advantages of the transformed closure problem are considerable.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the averaging volume, m² - B second-order tensor used to respresent the velocity deviation - b vector used to represent the pressure deviation, m^–1 - C second-order tensor related to the permeability tensor, m^–2 - D second-order tensor used to represent the velocity deviation, m² - d vector used to represent the pressure deviation, m - g gravity vector, m/s² - I unit tensor - K C ^–1,–D, Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - l _i i=1, 2, 3, lattice vectors, m - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of the-phase - p pressure in the-phase, N/m ² - p intrinsic phase average pressure, N/m² - p –p , spatial deviation of the pressure in the-phase, N/m² - r position vector locating points in the-phase, m - r ₀ radius of the averaging volume, m - t time, s - v velocity vector in the-phase, m/s - v intrinsic phase average velocity in the-phase, m/s - v phase average or Darcy velocity in the \-phase, m/s - v –v , spatial deviation of the velocity in the-phase m/s - V averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ Greek Letters V /V volume fraction of the-phase - mass density of the-phase, kg/m³ - viscosity of the-phase, Nt/m²     

Evaluation of the dynamic-mechanical properties of solids from a cooperative model for stress relaxation

For many solid materials the stress relaxation process obeys the universal relationF = – (d/d lnt)_max = (0.1 ± 0.01) ( ₀ — _i ), regardless of the structure of the material. Here denotes the stress,t the time, ₀ the initial stress of the experiment and _i the internal stress. A cooperative model accounting for the similarity in relaxation behaviour between different materials was developed earlier. Since this model has a spectral character, the concepts of linear viscoelasticity are used here to evaluate the corresponding prediction of the dynamic mechanical properties, i.e. the frequency dependence of the storageE () and lossE () moduli. Useful numerical approximations ofE () andE () are also evaluated. It is noted that the universal relation in stress relaxation had a counterpart in the frequency dependence ofE (). The theoretical prediction of the loss factor for high-density polyethylene is compared with experimental results. The agreement is good.     

Hypersonic flow over a narrow delta plate at large angles of attack

The problem of hypersonic flow over a flat delta plate with a high sweepback anglex at angles of attack close to /2 is solved using a numerical algorithm based on transition to the conical solution. The existence of conical flow at /2 with the velocity vector directed towards the apex of the plate is established. Values ofC _p/sin² and the thickness of the shock layer in the plane of symmetry of the plate are given as functions of the hypersonic similarity parameterk=tan tanx. A comparison of the calculated and experimental data shows that they are in good agreement.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.5, pp. 183–185, September–October, 1992.     

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.     

Uniqueness and the maximum principle for quasilinear elliptic equations with quadratic growth conditions

In this paper, we show that the maximum principle holds for quasilinear elliptic equations with quadratic growth under general structure conditions.Two typical particular cases of our results are the following. On one hand, we prove that the equation (1) {ie77-01} where {ie77-02} and {ie77-03} satisfies the maximum principle for solutions in H ¹()L(), i.e., that two solutions u ₁, u ₂H¹() L() of (1) such that u ₁u₂ on , satisfy u ₁u₂ in . This implies in particular the uniqueness of the solution of (1) in H ₀ ¹ ()L().On the other hand, we prove that the equation (2) {ie77-04} where fH^–1() and g(u)>0, g(0)=0, satisfies the maximum principle for solutions uH¹() such that g(u)¦Du|{²L¹(). Again this implies the uniqueness of the solution of (2) in the class uH ₀ ¹ () with g(u)¦Du|{²L¹().In both cases, the method of proof consists in making a certain change of function u=(v) in equation (1) or (2), and in proving that the transformed equation, which is of the form (3) {ie77-05}satisfies a certain structure condition, which using ((v₁ -v ₂)⁺)ⁿ for some n>0 as a test function, allows us to prove the maximum principle.     

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

A generalized Norton-Hoff model and the Prandtl-Reuss law of plasticity

The aim of this article is to study the quasistatic evolution of a three-dimensional elastic-perfectly plastic solid which satisfies the Prandtl-Reuss law. The evolution of the field of stresses -which solves a time dependent variational inequality — and that of the field of displacements u, have been described in previous works [15], [26], [35], [36], [37] but it was not shown there that and u satisfy indeed the Prandtl-Reuss constitutive law. In this article we find and u in a class of functions which are sufficiently regular for the Prandtl-Reuss law to make sense and we prove that and u satisfy the constitutive law. This result is attained by considering the elastic-perfectly plastic model as the limit of a family of elastic-visco-plastic models like those of Norton and Hoff. The Norton-Hoff type models which we introduce depend on a viscosity parameter > 0; we study the perturbed models (i.e. > 0 fixed) and then we pass to the limit 0.Dedicated to James Serrin on the occasion of his 60^th Birthday     

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     

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–)     

Heat conduction in a plane wall subject to periodic temperature changes

Analytical solutions for the heat conduction in a plane wall with periodic temperature variations at the wall surface are presented. Series and asymptotic developments of these solutions are deduced. The results are important for the calculation of the heat transfer in rotary kilns or other rotaring units.

---

Die Wärmeleitung in einer ebenen Wand mit periodischen Temperaturänderungen

Zusammenfassung Es werden analytische Lösungen für die Wärmeleitung in einer ebenen Wand mit periodischen Temperaturänderungen an ihrer Oberfläche mitgeteilt. Reihen- und asymptotische Entwicklungen dieser Lösungen werden abgeleitet. Die Ergebnisse sind wichtig für die Berechnung des Wärmetransportes in Drehrohröfen oder ähnlichen Maschinen.

---

Nomenclature a ² =/ C m²/s thermal diffusivity, Eq. (1) - C J/kg K specific heat - F K initial temperature of the wall, Eq. (4) - F m² surface of the wall - G Green's function, Eq. (10) - G₁ Green's function, Eq. (12) - h m thickness of the wall - H Heaviside function, Eq. (5) - k constant, Eq. (25) - k _x constant, Eq. (25) - Q J total energy, Eq. (17) - Q ^u J total energy from temperatureU, Eq. (18) - Q ^v J total energy from temperatureV, Eq. (19) - s s^–1 Laplace variable - t s time - t ₁ s heating time, Eq. (5) - t ₂ s period, Eq. (5) - T K temperature of the wall - T _i K surface temperature of the wall - T ₁ K surface temperature of the wall during the heating time - T ₂ K surface temperature of the wall during the cooling time - U K temperature of the wall defined in problem 1 - V, K temperature of the wall defined in problem 2 - x m coordinate - ₀ W/m²K overall heat transfer coefficient, Eq. (31) - ₁₀ W/m² K overall heat transfer coefficient, Eq. (32) - ₂₀ W/m² K overall heat transfer coefficient, Eq. (33) - Dirac Delta function - s^–1/2 parameter, Eq. (6) - W/mK thermal conductivity - kg/m³ specific mass - dimensionless time, Eq. (34) - Riemann Zeta function surfix Laplace transformed variable     

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     

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