Deformation to breaking of deep water gravity waves

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

The dynamic performance of the Weissenberg rheogoniometer III. Large amplitude oscillatory shearing — harmonic analysis

The harmonic content of the nonlinear dynamic behaviour of 1% polyacrylamide in 50% glycerol/water was studied using a standard Model R 18 Weissenberg Rheogoniometer. The Fourier analysis of the Oscillation Input and Torsion Head motions was performed using a Digital Transfer Function Analyser.In the absence of fluid inertia effects and when the amplitude of the (fundamental) Oscillation Input motion _I is much greater than the amplitudes of the Fourier components of the Torsion Head motion _Tn empirical nonlinear dynamic rheological propertiesG _n (, ⁰),G _n (, ⁰) and/or _n (, ⁰), _n (, ⁰) may be evaluated without a-priori-knowledge of a rheological constitutive equation. A detailed derivation of the basic equations involved is presented.Cone and plate data for the third harmonic storage modulus (dynamic rigidity)G ₃ (, ⁰), loss modulusG ₃ (, ⁰) and loss angle ₃ (, ⁰) are presented for the frequency range 3.14 × 10^–2 1.25 × 10² rad/s at two strain amplitudes, _CP ⁰ = 2.27 and 4.03. Composite cone and plate and parallel plates data for both the third and fifth harmonic dynamic viscosities ₃ (, ⁰), _S (, ⁰) and dynamic rigiditiesG ₃ (, ⁰),G ₅ (, ⁰) are presented for strain amplitudes in the ranges 1.10 _CP ⁰ 4.03 and 1.80 _PP ⁰ 36 for a single frequency, = 3.14 × 10^–1 rad/s. Good agreement was obtained between the results from both geometries and the absence of significant fluid inertia effects was confirmed by the superposition of the data for different gap widths.     

Loss of real characteristics for models of three-phase flow in a porous medium

In this paper we examine the generalized Buckley-Leverett equations governing threephase immiscible, incompressible flow in a porous medium, in the absence of gravitational and diffusive/dispersive effects. We consider the effect of the relative permeability models on the characteristic speeds in the flow. Using a simple idea from projective geometry, we show that under reasonable assumptions on the relative permeabilities there must be at least one point in the saturation triangle at which the characteristic speeds are equal. In general, there is a small region in the saturation triangle where the characteristic speeds are complex. This is demonstrated with the numerical results at the end of the paper.Symbols and Notation a, b, c, d entries of Jacobian matrix - A, B, C, D coefficients in Taylor expansion of _t, _v, _a - det J determinant of matrix J - dev J deviator of matrix J - J Jacobian matrix - L linear term in Taylor expansion for J near (s _v, s_a) = (0, 1) - m slope of r ₊ - p pressure - r± eigenvectors of Jacobian matrix - R real line - S intersection of saturation triangle with circle of radius centered at (1, 0) - S intersection of saturation triangle with circle of radius centered at (0, 1) - s _l, s_v, s_a saturations of phases (liquid, vapor, aqua) - tr J trace of matrix J - v _l , v _v , v _a phase flow rates (Darcy velocities) - v _T total flow rate - X, Y, Z entries of dev J - smooth closed curve inside saturation triangle - saturation triangle - _l, _v, _a phase density times gravitational acceleration times resevoir dip angle - K total permeability - _l, _v, _a three-phase relative permeabilities - _lv>, _la liquid phase relative permeabilities from two-phase data - _l, _v, _a mobilities of phases - _T total mobility - _l Corey mobility - _l, _v, _a phase viscosities - _± eigenvalues of Jacobian matrix - porosity Supported in part by National Science Foundation grant No. DMS-8701348, by Air Force Office of Scientific Research grant No. AFOSR-87-0283, and by Army Research Office grant No. DAAL03-88-K-0080.This work was performed under the auspices of the U.S. Department of Energy by the Lawrence Livermore National Laboratory under contract No. W-7405-Eng-48.     

Diffusion in anisotropic porous media

An experimental system was constructed in order to measure the two distinct components of the effective diffusivity tensor in transversely isotropic, unconsolidated porous media. Measurements were made for porous media consisting of glass spheres, mica particles, and disks made from mylar sheets. Both the particle geometry and the void fraction of the porous media were determined experimentally, and theoretical calculations for the two components of the effective diffusivity tensor were carried out. The comparison between theory and experiment clearly indicates that the void fraction and particle geometry are insufficient to characterize the process of diffusion in anisotropic porous media. Roman Letters A interfacial area between - and -phases for the macroscopic system, m² - A _e area of entrances and exits of the -phase for the macroscopic system, m² - A interfacial area contained within the averaging volume, m² - a characteristic length of a particle, m - b average thickness of a particle, m - c _A concentration of species A, moles/m³ - c _o reference concentration of species A, moles/m³ - c _A intrinsic phase average concentration of species A, moles/m³ - c _a c _A–c _A, spatial deviation concentration of species A, moles/m³ - C c _A/c ₀, dimensionless concentration of species A - binary molecular diffusion coefficient, m²/s - D _eff effective diffusivity tensor, m²/s - D _xx component of the effective diffusivity tensor associated with diffusion parallel to the bedding plane, m²/s - D _yy component of the effective diffusivity tensor associated with diffusion perpendicular to the bedding plane, m²/s - D _eff effective diffusivity for isotropic systems, m₂/s - f vector field that maps c _A on to c _a , m - h depth of the mixing chamber, m     

Marangoni-convection in a rotating liquid container

In a partially filled and constantly spinning container in zerogravity condition there arises under the action of an axial temperature gradient a thermo-capillary convection. This so-called Marangoni convection has been treated analytically for a directly imposed temperature gradient upon the free liquid surface and also for a constant but different temperature at the upper and lower disc wall. The streamfunction and circulation have been obtained, from which the velocity distribution could be determined.

---

Marangoni-Konvektion in einer in einem Behälter rotierenden Flüssigkeit

Zusammenfassung Durch das Vorhandensein eines axialen Temperaturgradienten ergibt sich in einem mit konstanter Geschwindigkeit rotierenden teilweise mit Flüssigkeit gefüllten Behälter eine thermalkapillare Korrelation. Diese sogenannte Marangoni-Konvektion wird analytisch behandelt für eine lineare axiale und eine beliebige axiale Temperaturverteilung auf der Flüssigkeitsoberfläche. Stromfunktion und Zirkulation werden analytisch bestimmt. Daraus ergeben sich die Geschwindigkeitsverteilungen in radialer, zirkumferentialer und axialer Richtung.

---

Nomenclature a radius of cylindrical container - b radius to free liquid surface - h height of container - I _m, K_m Modified Besselfunktions of first and second kind and orderm - k _j roots of bi-cubic equation (24 b) - k=b/a diameter ratio of location of free liquid surface and container wall - r, , z polar cylindrical coordinates - T(r, z) temperature distribution of liquid - u, v, w radial-, circumferential-, and axial velocity of the liquid, resp. - thermal expansion coefficient - dynamic viscosity of liquid - =/ kinematic viscosity - density of liquid - surface tension of liquid - _r , _rz shear stresses - (r, z) circulation - (r, z) stream function - ₀ speed of spin of container about axis of symmetry     

Stress analysis of the cylinder subjected to arbitrarily distributed axisymmetrical loads

Summary Stress analysis has been carried out for a finite cylinder subjected to arbitrarily distributed axisymmetrical surface loads. Direct stress _x in the axial direction is assumed to be of the form _x = ₀+r ₁ +r ₂ where ₀ to ₂ are functions of x. Using the equations of equilibrium and compatibility the other direct stresses and the shearing stress are expressed by ₁ and ₂. Fundamental equations governing ₁ and ₂ are introduced using the variational principle of complementary energy. From the results of the present analysis it is evident that the boundary conditions can be satisfied completely even for the case where the external forces are specified in complicated form, and that more accurate solutions can easily be obtained by introducing additional terms in _x.

---

Spannungsanalyse für den Zylinder unter axialsymmetrischer Last in beliebiger Verteilung

Übersicht Für einen endlichen Zylinder unter axialsymmetrischer Oberflächenlast in beliebiger Verteilung werden die Spannungen ermittelt. Die Normalspannung in Axialrichtung wird in der Form _x = ₀+r ₁ +r ₂ angesetzt mit ₀, ₁, ₂ als Funktionen von x. Mit Hilfe der Gleichgewichtsund Verträglichkeitsbedingungen werden die anderen Normalspannungen und die Schubspannung durch ₁ und ₂ ausgedrückt. Über das Variationsprinzip für die Komplementärenergie werden die grundlegenden Gleichungen für ₁ und ₂ eingeführt. Die Ergebnisse zeigen, daß die Randbedingungen selbst für komplizierte Belastungsarten vollständig erfüllbar sind und mit zusätzlichen Termen in _x mühelos noch genauere Lösungen bestimmt werden können.

     

A noninvasive method for the measurement of flow-induced surface displacement of a compliant surface

A noninvasive optical method is described which allows the measurement of the vertical component of the instantaneous displacement of a surface at one or more points. The method has been used to study the motion of a passive compliant layer responding to the random forcing of a fully developed turbulent boundary layer. However, in principle, the measurement technique described here can be used equally well with any surface capable of scattering light and to which optical access can be gained. The technique relies on the use of electro-optic position-sensitive detectors; this type of transducer produces changes in current which are linearly proportional to the displacement of a spot of light imaged onto the active area of the detector. The system can resolve displacements as small as 2 m for a point 1.8 mm in diameter; the final output signal of the system is found to be linear for displacements up to 200 m, and the overall frequency response is from DC to greater than 1 kHz. As an example of the use of the system, results detailing measurements obtained at both one and two points simultaneously are presented.List of symbols C _t elastic transverse wave speed = (G/)^1/2 - d ⁺ spot diameter normalized by viscous length scale - G frequency average of G() - G() shear storage modulus - G() shear loss modulus - l. viscous length scale = v/u _* - N total number of sampled data values - r separation vector for 2-point measurements = (, ) - rms root-mean-square value - R momentum thickness Reynolds number = U _t8/v - t time - u (y) mean streamwise component of velocity in boundary layer - u _* friction velocity = (t _w/)^1/2 - U free-stream velocity - x, y, z longitudinal, normal and spanwise directions - y _o undisturbed surface position - vertical component of compliant surface displacement - ₉₉ boundary layer thickness for which u(y) = 0.99 U _t8 - _l viscous sublayer thickness 5 l _* - frequency average of G()/ - boundary layer momentum thicknes = - fluid dynamic viscosity - v fluid kinematic viscosity = / - , longitudinal, spanwise components of separation vector r - fluid density - time delay - _w wall shear stress     

Hausdorff Dimension of Invariant Sets for Random Dynamical Systems

Suppose X() is a compact random set, invariant with respect to a continuously differentiable random dynamical system (RDS) on a separable Hilbert space. It is shown that the Hausdorff dimension dim _H (X()) is an invariant random variable, and it is bounded by d, provided the RDS contracts d-dimensional volumes exponentially fast. Both exponential decrease of d-volumes as well as the approximation of the RDS by its linearization are assumed to hold uniformly in . The results are applied to reaction diffusion equations with additive noise and to two-dimensional Navier–Stokes equations with bounded real noise.     

Green's functions,bi-linear forms,and completeness of the eigenfunctions for the elastostatic strip and wedge

The bi-harmonic Green's functionG(r,r) for the infinite strip region -1y1, -<x<, with the boundary conditionsG=G/y ony=±1, is obtained in integral form. It is shown thatG has an elegant bi-linear series representation in terms of the (Papkovich-Fadle) eigenfunctions for the strip. This representation is then used to show that any function bi-harmonic in arectangle, and satisfying the same boundary conditions asG, has a unique representation in the rectangle as an infinite sum of these eigenfunctions. For the case of the semi-infinite strip, we investigate conditions on sufficient to ensure that is exponentially small asx. In particular it is proved that this is so, solely under the condition that be bounded asx.A corresponding pattern of results is established for the wedge of general angle. The Green's function is obtained in integral form and expressed as a bilinear series of the (Williams) eigenfunctions. These eigenfunctions are proved to be complete for all functions bi-harmonic in anannular sector (and satisfying the same boundary conditions as the Green's function). As an application it is proved that if an elastostatic field exists in a corner region with free-free boundaries, and with either (i) the total strain energy bounded, or (ii) the displacement field bounded, then this field has a unique representation as a sum of those Williams eigenfunctions whichindividually posess the properties (i), (ii).The methods used here extend to all other linear homogeneous boundary conditions for these geometries.On leave of absence at the University of British Columbia, Vancouver, B.C. Canada, during 1977–79. This work was supported in part by N.R.C. grants Nos. A9259 and A9117.     

Rheology of molten polymers

Summary TheCross equation describes the flow of pseudoplastic liquids in terms of an upper and a lower Newtonian viscosity corresponding to infinite and zero shear, and ₀, and of a third material constant related to the mechanism of rupture of linkages between particles in the intermediate, non-Newtonian flow regime, Calculation of of bulk polymers is important, since it cannot be determined experimentally. The equation was applied to the melt flow data of two low density polyethylenes at three temperatures.Using data in the non-Newtonian region covering 3 decades of shear rate to extrapolate to the zero-shear viscosity resulted in errors amounting to about onethird of the measured ₀ values. The extrapolated upper Newtonian viscosity was found to be independent of temperature within the precision of the data, indicating that it has a small activation energy.The ₀ values were from 100 to 1,400 times larger than the values at the corresponding temperatures.The values of were large compared to the values found for colloidal dispersions and polymer solutions, but decreased with increasing temperature. This shows that shear is the main factor in reducing chain entanglements, but that the contribution of Brownian motion becomes greater at higher temperatures.

---

Zusammenfassung Die Gleichung vonCross beschreibt das Fließverhalten von pseudoplastischen Flüssigkeiten durch drei Konstante: Die obereNewtonsche Viskosität (bei sehr hohen Schergeschwindigkeiten), die untereNewtonsche Viskosität ₀ (bei Scherspannung Null), und eine Materialkonstante, die vom Brechen der Bindungen zwischen Partikeln im nicht-Newtonschen Fließbereich abhängt. Die Berechnung von ist wichtig für unverdünnte Polymere, wo man sie nicht messen kann.Die Gleichung wurde auf das Fließverhalten der Schmelzen von zwei handelsüblichen Hochdruckpolyäthylenen bei drei Temperaturen angewandt. Die Werte von ₀, durch Extrapolation von gemessenen scheinbaren Viskositäten im Schergeschwindigkeitsbereich von 10 bis 4000 sec^–1 errechnet, wichen bis 30% von den gemessenen ₀-Werten ab. Die Aktivierungsenergie der war so klein, daß die-Werte bei den drei Temperaturen innerhalb der Genauigkeit der Extrapolation anscheinend gleich waren. Die ₀-Werte waren 100 bis 1400 mal größer als die-Werte.Im Verhältnis zu kolloidalen Dispersionen und verdünnten Polymerlösungen war das der Schmelzen groß, nahm aber mit steigender Temperatur ab. Deshalb wird die Verhakung der Molekülketten hauptsächlich durch Scherbeanspruchung vermindert, aber der Beitrag derBrownschen Bewegung nimmt mit steigender Temperatur zu.

     

Temperature-dependent heat sources or sinks in a stagnation point flow

Summary An analytical study has been made to determine the heat transfer characteristics of a stagnation point flow in which there are temperature-dependent heat sources or sinks. Results have been obtained for both strong and weak sources or sinks for a Prandtl number of 0.7. An analytical method, applicable to all Prandtl numbers, was utilized which circumvented the need for extensive numerical solutions and which, at the same time, provided a closed-form representation for the heat transfer. A few numerical solutions were carried out to verify the method.Nomenclature a _i constants depending on Prandtl number - c _p specific heat at constant pressure - f dimensionless velocity variable - g function defined by equation (13) - g _n functions of (n=1, 2, 3,...) - k thermal conductivity - Pr Prandtl number, c _p /k - q heat transfer rate per unit area at surface - Q heat flux parameter, q/k(u ₁/)^1/2 - S rate of heat generation or removal per unit volume (divided by c _p ) - T static temperature; T _w , wall temperature; T , free-stream temperature - u ₁ proportionality constant for free-stream velocity - U free-stream velocity - v normal velocity component - x coordinate measuring distance along surface from stagnation point - y coordinate measuring distance normal to surface - heat generation parameter, equation (3) - dimensionless normal coordinate, - dimensionless temperature - _n functions of (n=1, 2, 3,...) - absolute viscosity - kinematic viscosity - density     

Conditions of simulation of stagnation point heat transfer from a high-enthalpy flow

The problem of local simulation of stagnation point heat transfer to a blunt body is solved within the framework of boundary layer theory on the assumption that the simulation subsonic high-enthalpy flow is in equilibrium outside the boundary layer on the model, while the parameters of the natural flow are in equilibrium at the outer edge of the boundary layer on the body. The parameters of the simulating subsonic flow are expressed in terms of the total enthalpyH ₀, the stagnation point pressurep _w and the velocityV ₁ for the natural free-stream flow in the form of universal functions of the dimensionless modeling coefficients=R _m ^* /R _b ^* ( .<1),=V ₁/2H ₀ ( .<1) whereR _m ^* and R _b ^* are the effective radii of the model and the body at their stagnation points. Approximate conditions for modeling the heat transfer from a high-enthalpy (including hypersonic) flow to the stagnation point on a blunt body by means of hyposonic (M1) flows, corresponding to the case ²1, are obtained. The possibilities of complete local simulation of hypersonic nonequilibrium heat transfer to the stagnation point on a blunt body in the hyposonic dissociated air jets of a VGU-2 100-kilowatt induction plasma generator [4, 5] are analyzed.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.1, pp. 172–180, January–February, 1993.     

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

Buoyancy Sustained by Viscous Dissipation

The quasi-parallel regime of a Darcy–Boussinesq boundary-layer flow over a permeable vertical flat plate adjacent to a fluid saturated porous medium is considered. Quasi-parallel means here a plane flow with a constant transversal velocity v=–v ₀ directed perpendicularly towards the vertical surface, where a lateral suction with the same velocity –v ₀ is applied. The plate is held at a constant temperature T _w which coincides with the ambient temperature T of the fluid. The heat released by viscous dissipation induces a density gradient in the fluid. Thus, although T _w=T , a thermal convection occurs. The steady regime of this self-sustaining buoyant flow has been examined in detail. Wall jet-like profiles with a continuous but finite spectrum of the momentum flow have been found. These self-sustaining buoyant jets show a universal behavior, that is, there exist certain length, velocity and temperature scales such that the flow characteristics become independent of the (constant) material properties of the fluid and the porous medium as well.     

Das erweiterte Korrespondenzgesetz für die dynamische Idealviskosität und die Idealwärmeleitfähigkeit reiner Stoffe

Zusammenfassung Zur Berechnung der dynamischen Idealviskosität Ideal (T) und der Idealwärmeleitfähigkeit ideal (T) benötigt man die kritische TemperaturT _kr, das kritische spezifische Volum _kr, die MolmasseM, den kritischen Parameter _kr und die molare isochore WärmekapazitätC _v(T). Sowohl das theoretisch, als auch das empirisch abgeleitete erweiterte Korrespondenzgesetz ergeben eine für praktische Zwecke ausreichende Genauigkeit für die Meßwertwiedergabe, die bei den assoziierenden Stoffen und den Quantenstoffen jedoch geringer ist als bei den Normalstoffen.

---

The extended correspondence law for the ideal dynamic viscosity and the ideal thermal conductivity of pure substances

For the calculation of the ideal dynamic viscosity Ideal (T) and the ideal thermal conductivity ideal (T) the critical temperatureT _kr, the critical specific volumev _kr, the molecular massM, the critical parameter _kr, and the molar isochoric heat capacityC _v(T) is needed. Not only the theoretically determined but also the empirically determined extended correspondence law gives for practical use a good representation of the measured data, which for the associating substances and the quantum substances is not so good as for the normal substances.

     

Elongational flow behavior of viscoelastic liquids: Modelling bubble dynamics with viscoelastic constitutive relations

Summary Earlier parts of this series have described a technique based on the collapse of single bubbles in the fluids for studying the elongational rheology of viscoelastic solutions and melts of moderate viscosities ( ₀ > 10²p) at relatively high strain rates . The present paper describes the modelling of bubble collapse with both rate and integral type constitutive relations using a body coordinate system. Predictions of the stress at the bubble wall as a function of time during collapse from a BKZ model and a modified corotational Maxwell model compared favorably with experimental data for two polymer solutions, 1% polyacrylamide in water/glycerine and 2% hydroxypropyl cellulose in water.

---

Zusammenfassung In vorangehenden Veröffentlichungen dieser Reihe wurde eine Methode beschrieben, mit Hilfe derer man aus dem Zerfall von einzelnen Blasen in einer Flüssigkeit auf die Dehn-Rheologie viskoelastischer Lösungen und Schmelzen mittlerer Viskosität ( ₀ > 10² P) bei relativ hohen Dehngeschwindigkeiten schließen kann. Die vorliegende Untersuchung beschreibt Modelle des Blasenzerfalls mit Hilfe von Stoffgleichungen sowohl vom rate- als auch vom Integral-Typ, wobei ein körperfestes Koordinatensystem benutzt wird. Die Voraussagen der Spannung an der Blasenwand als Funktion der Zeit während des Zerfalls bei Verwendung eines BKZ- und eines modifizierten korotatorischen Maxwell-Modells zeigen eine recht gute Übereinstimmung mit experimentellen Werten, die an zwei Polymerlösungen, nämlich einer 1%igen Polyacrylamid-Lösung in einer Wasser-Glycerin-Mischung und einer 2%igen wäßrigen Hydropropylcellulose, erhalten worden sind.

---

Nomenclature a material constant - b material constant - g metric tensor, space coordinates - m material constant - n material constant - p pressure - P _G pressure within bubble - P _R pressure outside bubble at the wall - P pressure far away from the bubble - R bubble radius - dR/dt - R ₀ initial bubble radius - t time - u velocity - U potential function - Y R/R ₀ Greek symbols covariant body metric tensor - surface tension - rate of deformation matrix, II -second invariant of - strain rate - ₀ zero shear rate viscosity - _e elongational viscosity - _ef effective viscosity - ¹, ², ³ coordinates in body system - ^{1 1}/R ₀ ³ - body stress tensor - density - space stress tensor - relaxation time - _ef effective relaxation time - bubble pressure function, defined in eq. [19] - vorticity tensor With 11 figures and 1 table     

Combined natural and forced convection in a laminar wall jet along a vertical plate with uniform surface heat flux

The steady state flow and heat transfer characteristics of the combined natural and forced convection in a two dimensional, laminar, incompressible wall jet over a vertical wall are obtained for constant wall heat flux boundary condition. The velocity and temperature distribution are assumed to be power series, where the zeroth term corresponds to that for a plane wall jet in the absence of buoyancy effects. Numerical results for the momentum and thermal series functions are presented for a Prandtl number of 0.73. Wall values of the momentum and thermal series functions are presented for Prandtl numbers ranging from 0.01 to 1000.Nomenclature Gr* modified Grashof number - k thermal conductivity - Nu Nusselt number - Pr Prandtl number - q _w heat flux at the wall - Re Reynolds number - T temperature - u velocity component in x-direction - v velocity component in y-direction - x co-ordinate along the plane wall - y co-ordinate normal to the wall - () gamma function - non-dimensional co-ordinate defined in (6) - non-dimensional temperature - dynamic viscosity - kinematic viscosity - non-dimensional co-ordinate defined in (6) - density - w values at the wall - values at large distances away from the wall     

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     

Flow of a simple non-newtonian fluid past a sphere

Summary Creeping flow past a sphere is solved for a limiting case of fluid behaviour: an abrupt change in viscosity.List of Symbols d _ij Component of rate-of-deformation tensor - F _d Drag force exerted on sphere by fluid - G ^(d) Coefficients in expression for _ij in terms of d _ij - G _YOJK ^(d) Coefficients in power series representing G ^(d) - R Radius of sphere - r Spherical coordinate - V Velocity of fluid very far from sphere - v _i Component of the velocity vector - x Dimensionless radial distance, r/R - x _i Rectangular Cartesian coordinate - Dimensionless quantity defined by (26) - ^(d) Potential defined by (7) - Value of x denoting border between Regions 1 and 2 as a function of - ₁, ₂ Lower and upper limiting viscosities defined by (10) - Spherical coordinate - * Value of for which =1 - Value of denoting border between regions 1 and 2 as a function of x - Newtonian viscosity - _ij Component of the stress tensor - Spherical coordinate - ₁, ₂ Stream functions defined by (12) and (14) - Second and third invariants of the stress tensor and of the rate-of-deformation tensor, defined by (3)     

Method of mass-average quantities for the boundary layer in an outer flow with transverse nonuniformity

An engineering method is proposed for calculating the friction and heat transfer through a boundary layer in which a nonuniform distribution of the velocity, total enthalpy, and static enthalpy is specified across the streamlines at the initial section x₀. Such problems arise in the vortical interaction of the boundary layer with the high-entropy layer on slender blunt bodies, with sudden change of the boundary conditions for an already developed boundary layer (temperature jump, surface discontinuity), and in wake flow past a body, etc.Notation x, y longitudinal and transverse coordinates - u,, H, h gas velocity, stream function, total and static enthalpy - p,,, pressure, density, viscosity, Prandtl number - , q friction and thermal flux at the body surface - r(x), (x) body surface shape and boundary layer thickness - V, M freestream velocity and Mach number - u⁽⁰⁾(x₀,), H⁽⁰⁾(x₀,), h⁽⁰⁾(x₀,) parameter distributions at initial section - u⁽⁰⁾(x,), h⁽⁰⁾(x,), h⁽⁰⁾(x,) profiles of quantities in outer flow in absence of friction and heat transfer at the surface of the body The indices v=0, 1 relate to plane and axisymmetric flows - , w, b, relate to quantities at the outer edge of the inner boundary layer, at the body surface in viscid and nonviscous flows, and in the freestream, respectively. The author wishes to thank O. I. Gubanov, V. A. Kaprov, I. N. Murzinov, and A. N, Rumynskii for discussions and assistance in this study.