On the compactness of quasi-conforming element spaces and the convergence of quasi-conforming element method

In this paper, the compactness of quasi-conforming element spaces and the—convergence of quasi-conforming element method are discussed. The well-known Rellich compactness theorem is generalized to the sequences of quasi-conforming element spaces with certain properties, and the generalized Poincare inequality. The generalized Friedrichs inequality and the generalized inequality of Poincare-Friedrichs are proved true for them. The error estimates are also given. It is shown that the quasi-conforming element method is convergent if the quasi-conforming element spaces have the approximability and the strong continuity, and satisfy the rank condition of element and pass the test IPT. As practical examples, 6-parameter, 9-paramenter, 12-paramenter, 15-parameter, 18-parameter and 21-paramenter quasi-conforming elements are shown to be convergent, and their L²²()-errors are O(h), O(h), O(h ² ), O(h ² ), O(h ), and O(h ⁴ ) respectively.     

Analysis of melt spinning transients in Lagrangian coordinates

Transients in melt spinning of isothermal power law and Newtonian fluids were found to be governed by an extremely simple partial differential equation ² ( ^1/n )/() = 0 in Lagrangian coordinates where is the cross-sectional area,n the power law exponent, the time and the the time at which a fluid molecule constituting the spinline left the spinneret. The general integral ^1/n =f() +g () of the above governing equation containing two arbitrary functions represents physically attainable spinline transients. Hitherto unknown analytical transient solutions of the above governing equation were obtained for the response of isothermal constant tension spinlines to a stepwise change in tension, spinneret hole area, extrusion speed or extrusion viscosity and for the starting transient in gravitational spinning. Linearized perturbation solutions and the stability limit of the spinline derived from the above new found nonlinear solutions were in agreement with previous findings and the above nonlinear response of the spinline to a step increase in the spinneret hole area was found to be equivalent to Orowan's tandem cylinder model of dent growth in filament stretching.     

Turbulent boundary layer over a smooth wall with widely separated transverse square cavities

Laser-Doppler velocimetry (LDV) measurements and flow visualizations are used to study a turbulent boundary layer over a smooth wall with transverse square cavities at two values of the momentum thickness Reynolds number (R =400 and 1300). The cavities are spaced 20 element widths apart in the streamwise direction. Flow visualizations reveal a significant communication between the cavities and the overlying shear layer, with frequent inflows and ejections of fluid to and from cavities. There is evidence to suggest that quasi-streamwise near-wall vortices are responsible for the ejections of fluid out of the cavities. The wall shear stress, which is measured accurately, increases sharply immediately downstream of the cavity. This increase is followed by a sudden decrease and a slower return to the smooth wall value. Integration of the wall shear stress in the streamwise direction indicates that there is an increase in drag of 3.4% at bothR .Nomenclature C _f skin friction coefficient - C _fsw friction coefficient for a continuous smooth wall - k height of the cavity - k ⁺ ku / - R Reynolds number based on momentum thickness (U ₁ /v) - Rx Reynolds number based on streamwise distance (U ₁ x/) - s streamwise distance between two cavities - t time - t ⁺ tu ₂ / - U ₁ freestream velocity - mean velocity inx direction - u,v,w rms turbulent intensities inx,y andz directions - u local friction velocity - u _sw friction velocity for a continuous smooth wall - w width of the cavity - x streamwise co-ordinate measured from the downstream edge of the cavity - y co-ordinate normal to the wall - z spanwise co-ordinate - y ⁺ yu / - boundary layer thickness - ₀ boundary layer thickness near the upstream edge of the cavity - _i thickness of internal layer - kinematic viscosity of water - ⁺ zu / - momentum thickness     

An analytical solution for the velocity and shear rate distribution of non-ideal Bingham fluids in concentric cylinder viscometers

An analytical solution is presented for the calculation of the flow field in a concentric cylinder viscometer of non-ideal Bingham-fluids, described by the Worrall-Tuliani rheological model. The obtained shear rate distribution is a function of the a priori unknown rheological parameters. It is shown that by applying an iterative procedure experimental data can be processed in order to obtain the proper shear rate correction and the four rheological parameters of the Worrall-Tuliani model as well as the yield surface radius. A comparison with Krieger's correction method is made. Rheometrical data for dense cohesive sediment suspensions have been reviewed in the light of this new method. For these suspensions velocity profiles over the gap are computed and the shear layer thicknesses were found to be comparable to visual observations. It can be concluded that at low rotation speeds the actually sheared layer is too narrow to fullfill the gap width requirement for granular suspensions and slip appears to be unavoidable, even when the material is sheared within itself. The only way to obtain meaningfull measurements in a concentric cylinder viscometer at low shear rates seems to be by increasing the radii of the viscometer. Some dimensioning criteria are presented.Notation A, B Integration constants - C Dimensionless rotation speed = µ/_y - c = 2µ - d = ₀ ²–2c_y - f() = (–₀)²+2c(–_y) - r Radius - r _b Bob radius - r _c Cup radius - r _y Yield radius - r ₀ Stationary surface radius - r Rotating Stationary radius - Y ₀ Shear rate parameter = /µ Greek letters Shear rate - = (r _y /r _b )²– 1 - µ Bingham viscosity - µ₀ Initial differential viscosity - µ µ₀-µ - Rotation speed - Angular velocity - Shear stress - _b Bob shear stress - _B Bingham stress - _y (True) yield stress - ₀ Stress parameter = _B +µ Y ₀ - _B - _y      

Laser velocimetry measurements in a horizontal gas-solid pipe flow

This paper presents laser measurements of particle velocities in a horizontal turbulent two-phase pipe flow. A phase Doppler particle analyzer, (PDPA), was used to obtain particle size, velocity, and rms values of velocity fluctuations. The particulate phase consisted of glass spheres 50 m in diameter with the volume fraction of the suspension in the range _p=10^-4 to _p=10^-3. The results show that the turbulence increases with particle loading.List of symbols a particle diameter - C _va velocity diameter cross-correlation - d pipe diameter - Fr ² Froude number - g gravitational constant - p(a) Probability density of the particle diameter - Re pipe Reynolds number based on the friction velocity - T characteristic time scale of the energy containing eddies - T _L integral scale of the turbulence sampled along the particle path - u, U, u characteristic fluid velocities: fluctuating, mean and friction - v characteristic velocity of the paricle fluctuations - f expected value of any random variable f - f¦g expected value of f given a value of the random variable g - _p particle volume fraction - _p particle response time - absolute fluid viscosity - v kinematic fluid viscosity - _p, _f densities, particle and fluid - _a ² particle diameter variance - _va ² velocity variance due to the particle diameter variance - _vT ² total particle velocity variance - _vt ² particle velocity variance due to the response to the turbulent field     

Modellvorstellung zur turbulenten Filmkondensation strömenden Dampfes

Zusammenfassung Der Wärmeübergang bei turbulenter Film kondensation strömenden Dampfes an einer waagerechten ebenen Platte wurde mit Hilfe der Analogie zwischen Impuls-und Wärmeaustausch untersucht. Zur Beschreibung des Impulsaustausches im Film wurde ein Vierbereichmodell vorgestellt. Nach diesem Modell wird die wellige Phasengrenze als starre rauhe Wand angesehen. Die Abhängigkeit einer Schubspannungs-Nusseltzahl von der Film-Reynoldszahl und Prandtlzahl wurde berechnet und dargestellt.

---

A model for turbulent film condensation of flowing vapour

The heat transfer in turbulent film condensation of flowing vapour on a horizontal flat plate was investigated by means of the analogy between momentum and heat transfer. To describe the momentum transfer in the film a four-region model was presented. With this model the wavy interfacial surface is treated as a stiff rough wall. A shear Nusselt number has been calculated and represented as a function of film Reynolds number and Prandtl number.

---

Formelzeichen a Temperaturleitkoeffizient - k Mischungswegkonstante - k _s äquivalente Sandkornrauhigkeit - Nu _x lokale Schubspannungs-Nusseltzahl,Nu _x=x_xv/uw - Pr Prandtlzahl,Pr=v/a - Pr _t turbulente Prandtlzahl,Pr _t =_m/_q - q Wärmestromdichte _q - R Wärmeübergangswiderstand - R_f Wärmeübergangswiderstand des Films - Re _F Reynoldszahl der Filmströmung - T Temperatur - U, V Geschwindigkeitskomponenten des Dampfes in waagerechter und senkrechter Richtung - u, Geschwindigkeitskomponenten des Kondensats in waagerechter und senkrechter Richtung - V Querschwankungsgeschwindigkeit des Kondensats und des Dampfes - u _/gtD Schubspannungsgeschwindigkeit an der Phasengrenze für die Dampfgrenzschicht, uD =(/)^1/2 - u _F Schubspannungsgeschwindigkeit an der Phasengrenze für den Kondensatfilm,u _F =(/)^1/2 - u _w Schubspannungsgeschwindigkeit an der Wand der Kühlplatte,u _w =(w/)^1/2 - y Wandabstand - x Wärmeübergangskoeffizient - gemittelte Kondensatfilmdicke - _s Dicke der zähen Schicht der Filmströmung an der welligen Phasengrenze - ₄ Dicke der zähen Schicht der Filmströmung an der gemittelten glatten Phasengrenze - Wärmeleitzahl - dynamische Viskosität - v kinematische Viskosität - Dichte - Oberflächenspannung - _w Wandschubspannung - Schubspannung an der Phasengrenzfläche - _m turbulente Impulsaustauschgröße - _q turbulente Wärmeaustauschgröße Indizes d Wert des Dampfes - w Wert an der Wand - x lokaler Wert inx - Wert an der Phasengrenze Stoffgrößen ohne Index gelten für das Kondensat     

On properties of solutions of functional differential equations continuously differentiable on (0, +∞)

We establish the saddle-point property of the system of functional differential equations (t) = Ax(t) + Bx((t)) + C ((t)) + f (x(t), x((t))), (0) = 0.Translated from Neliniini Kolyvannya, Vol. 7, No. 3, pp. 302–310, July–September, 2004.     

An experimental study on effects of solid concentration and ζ-potential on pseudoplastic flow of suspensions

The pseudoplastic flow of suspensions, alumina or styrene-acrylamide copolymer particles in water or an aqueous solution of glycerin has been studied by the step-shear-rate method. The relation between the shear rate,D, and the shear stress,, in the step-shear-rate measurements, where the state of dispersion was considered to be constant, was expressed as = AD ^1/2 +CD. The effective solid volume fraction,ø _F, andA were dependent on the shear rate and expressed byø _F =aD ^b andA = D . Combining the above relations, the steady flow curve was expressed by = D ^{1/2 +} + ₀ (1 – a D ^b/0.74)^–1.85 D, where ₀ is the viscosity of the medium.With an increase in solid volume fraction and a decreases in the absolute value of the-potential, the flow behavior of the suspensions changed from Newtonian ( = = b = 0), slightly pseudoplastic ( = b = 0), pseudoplastic ( = 0) to a Bingham-like behavior.The change in viscosity of the medium had an effect on the change in the effective volume fraction.     

Rheological aspects of the hardening of water-borne alkali silicate zinc-rich paints

Summary The kinetics of the reaction between zinc dust and alkali silicate in water-borne silicate zinc-rich paints was investigated on the basis of rheological measurements. The paints studied were formulated with the same zinc dust/silicate vehicle ratio. Different zinc dusts were employed. The vehicle was a litiumsodium silicate solution. Formulation included a rheological agent to prevent rapid settling of suspended zinc dust.Rheological tests were carried out at various time intervals after paint mixing. Each sample was subjected, stepwise, to increasing and decreasing shear rate sequences, each shear rate being applied until stress steady values were attained. Shear rates were within the range 2.35 to 1700 s^–1; subsequently, hysteresis cycles were traced. Equilibrium data were fitted with the Bingham equation.Yield values and ultimate viscosities obtained at the various test time intervals were compared. Both parameters were found to increase with increasing time after mixing. Their increase in time brings out the fact that two successive processes take place; accordingly, structural hypotheses were suggested taking into account the modification set up in the system by the zincsilicate reaction.

---

Zusammenfassung Die Kinetik der Reaktion zwischen Zinkstaub und Alkalisilikat in wäßrigen Silikat-Zink-Farben wurde mit Hilfe rheologischer Messungen erforscht. Die untersuchten Farben wurden stets mit dem gleichen Verhältnis von Zinkstaub zu Silikat-Träger formuliert. Dabei wurden verschiedene Zinkstäube eingesetzt. Der Träger bestand aus einer Lithium-Natrium-Silikat-Lösung. In der Formulierung war außerdem ein rheologischer Wirkstoff, der das schnelle Absetzen des suspendierten Zinkstaubes verhindern sollte, enthalten.Die rheologischen Messungen wurden zu verschiedenen Zeiten nach dem Farbmischen durchgeführt. Jede Probe wurde einer auf- und absteigenden Stufenfolge von Schergeschwindigkeiten unterworfen, wobei die Scherung solange konstant gehalten wurde, bis sich ein stationärer Wert eingestellt hatte. Die Schergeschwindigkeiten lagen dabei in einem Bereich zwischen 2,35 und 1700 s^–1. Anschließend wurden Hysteresisschleifen aufgenommen.Die Gleichgewichtswerte ließen sich mit der Bingham-Gleichung beschreiben. Fließgrenzen und Endviskositäten, die nach verschiedenen Zeitdauern gemessen worden waren, wurden miteinander verglichen. Beide Parameter nahmen mit der Zeit zu. Daraus ist erkennbar, daß zwei aufeinanderfolgende Prozesse stattfinden. Dementsprechend werden Strukturmodelle vorgeschlagen, welche die Veränderungen in Rechnung stellen, die als Folge der Zink-Silikat-Reaktionen in dem System stattfinden.

---

List of symbols shear rate - ÿ shear acceleration - _r relative viscosity - _R reduced viscosity - ₀ continuous phase viscosity - , _,0, _,t ultimate shear rate, at the beginning of the reaction, at timet - shear stress - _R reduced shear stress - ₀, _0,0, _0,t yield stress, at the beginning of the reaction, at timet - , _max disperse phase volume fraction, maximum volume fraction - a, a , a shift factor, viscosity shift factor, shear stress shift factor - a _1,, a_2, see eq. [4] - k ₁,k₂ reaction rate constants - K _E see eq. [5] - t time - t _C characteristic time Abbreviations MSD mean square deviation (=( _exp — _calc)²/(n — K);n = number of experimental points,K = number of parameters) Paper presented at the Joint Conference of the British, Italian and Netherlands Societies of Rheology, Amsterdam, April 18–20, 1979.With 6 figures and 3 tables     

Image processing of hydrogen bubble flow visualization for determination of turbulence statistics and bursting characteristics

A technique is described which employs automated image processing of hydrogen-bubble flow visualization pictures to establish local, instantaneous velocity profile information. Hydrogen bubble flow visualization sequences are recorded using a high-speed video system and then digitized, stored, and evaluated by a VAX 11/780 computer. Employing special smoothing and gradient detection algorithms, individual bubble-lines are computer identified, which allows local velocity profiles to be constructed using time-of-flight techniques. It is demonstrated how this techniques may be used to 1) determine local velocity behavior as a function of position and time, 2) evaluate time-averaged turbulence properties, and 3) correlate probe-type turbulent burst detection techniques with the corresponding visualization data.List of symbols Re Reynolds number based on momentum thickness, u / - t ⁺ nondimensional time tu ² / - T VITA variance averaging time period - u shear velocity = - u local instantaneous streamwise velocity,x-direction - u local fluctuating streamwise velocity,x-direction - u ⁺ nondimensional streamwise velocity, /u - local normal velocity,y-direction - w local spanwise velocity,z-direction - x ⁺ nondimensional coordinate in streamwise direction xu /v - y ⁺ nondimensional coordinate normal to wall, yu /v Greek momentum thickness, - kinematic viscosity - _w wall shear stress This paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984     

A three-parameter model describing the experimental relation between shear stress and shear rate for laminar flow of aqueous polymer solutions

Summary A three-parameter model is introduced to describe the shear rate — shear stress relation for dilute aqueous solutions of polyacrylamide (Separan AP-30) or polyethylenoxide (Polyox WSR-301) in the concentration range 50 wppm – 10,000 wppm. Solutions of both polymers show for a similar rheological behaviour. This behaviour can be described by an equation having three parameters i.e. zero-shear viscosity ₀, infinite-shear viscosity , and yield stress ₀, each depending on the polymer concentration. A good agreement is found between the values calculated with this three-parameter model and the experimental results obtained with a cone-and-plate rheogoniometer and those determined with a capillary-tube rheometer.

---

Zusammenfassung Der Zusammenhang zwischen Schubspannung und Schergeschwindigkeit von strukturviskosen Flüssigkeiten wird durch ein Modell mit drei Parametern beschrieben. Mit verdünnten wäßrigen Polyacrylamid-(Separan AP-30) sowie Polyäthylenoxidlösungen (Polyox WSR-301) wird das Modell experimentell geprüft. Beide Polymerlösungen zeigen im untersuchten Schergeschwindigkeitsbereich von ein ähnliches rheologisches Verhalten. Dieses Verhalten kann mit drei konzentrationsabhängigen Größen, nämlich einer Null-Viskosität ₀, einer Grenz-Viskosität und einer Fließgrenze ₀ beschrieben werden. Die Ergebnisse von Experimenten mit einem Kegel-Platte-Rheogoniometer sowie einem Kapillarviskosimeter sind in guter Übereinstimmung mit den Werten, die mit dem Drei-Parameter-Modell berechnet worden sind.

---

a Pa^–1 physical quantity defined by:a = {1 – ( / ₀)}/ ₀ - c l concentration (wppm) - D m capillary diameter - L m length of capillary tube - P Pa pressure drop - R m radius of capillary tube - u m s^–1 average velocity - v _r m s^–1 local axial velocity at a distancer from the axis of the tube - shear rate (–dv _r /dr) - local shear rate in capillary flow - s^–1 wall shear rate in capillary flow - Pa s dynamic viscosity - _a Pa s apparent viscosity defined by eq. [2] - ( _a ) Pa s apparent viscosity in capillary tube at a distanceR from the axis - ₀ Pa s zero-shear viscosity defined by eq. [4] - Pa s infinite-shear viscosity defined by eq. [5] - l ratior/R - kg m^– density - Pa shear stress - ₀ Pa yield stress - _r Pa local shear stress in capillary flow - _R Pa wall shear stress in capillary flow _R = (PR/2L) - _v m³ s^–1 volume rate of flow With 8 figures and 1 table     

Über Differentialgleichungen für Wärme- und Stoffaustausch von hygroskopischen Stoffen

Zusammenfassung Die beiden Differentialgleichungssysteme vonKrischer undLykow werden miteinander verglichen. Dabei ergibt sich, daß die in der deutschen und russischen Literatur angewandten mathematischen Modelle der Trocknung von kapillarporösen Körpern praktisch übereinstimmen. Es werden die Transformationsgleichungen der dimensionslosen Kenngrößen angegeben, die die Beziehungen zwischen den beiden Systemen herstellen.

---

The differential equations ofKrischer andLuikow for unsteady internal heat and mass transfer in the porous medium are compared. It is shown, that the mathematical models for drying in the German and Russian literature are equivalent. The transform relations of the non-dimensional parameters between the two models are given.

---

Formelzeichen nach Krischer z laufende Koordinate in Strömungsrichtung in m - R kennzeichnende Abmessung des Körpers in m - t Zeit in h - Raumgewicht bei mittlerer Feuchtigkeit in kg/m³ - _w Teilgewicht des Wassers in 1 m³ Trockengut in kg/m³ - _wa Anfangsfeuchtigkeit in kg/m³ - _D Dampfdichte in kg/m³ - _L Luftvolumen je m³ Trocknungsgut in m³/m³ - Temperatur in °C - _u Umgebungstemperatur in °C - _a Anfangstemperatur in °C - r Verdampfungswärme in kcal/kg - q _E Wärmeentwicklung in kcal/m³ h - c spezifische Wärmekapazität des Trockengutes in kcal/kg grd - Wärmeleitfähigkeit in kcal/m h grd - Feuchtigkeitsleitzahl des Trockengutes in m²/h - wirksame Diffusionszahl von Wasserdampf in Luft in m²/h - Diffusionswiderstandszahl des Trockengutes — - Konstante — - Konstante in kg/m³ grd Formelzeichen nach Lykow X=r/R dimensionslose Koordinate des Körpers;r laufende Koordinate in m;R kennzeichnende Abmessung in m; - Fo=a/R ² Fourier-Zahl;a Temperaturleitzahl in m²/h; Zeit in h - T(X, Fo)=t(r, )– 0/t dimensionslose Temperatur des Körpers im Punkt mit KoordinateX für den ZeitpunktFo;t(r, ) Temperatur in °C; 0 mittlere Anfangstemperatur in °C; t ein vorher angenommener Temperaturunterschied in grd - (X, Fo)= ₀–u(r, )/u dimensionsloses Potential des Stoffübergangs im Punkt mit KoordinateX für den ZeitpunktFo;u(r, ) Feuchtigkeitsgehalt des Trockengutes in kg/kg; ₀ mittlerer Anfangsfeuchtigkeitsgehalt in kg/kg; u ein vorher angenommener Unterschied des Feuchtigkeitsgehalts in kg/kg - Ko= u/c t Kosowitsch-Zahl; Verdampfungswärme in kcal/kg;c spezifische Wärmekapazität des Trockengutes in kcal/kg - Ko*=Ko modifizierte Kosowitsch Zhal; Kenngröße der Zustandsänderung - Pn= t/u Posnowsche Zahl;=a _T ^m /a _m Thermogradientkoeffizient in 1/grd;a _T ^m thermische Stoffübergangszahl (charakterisiert den Stoffstrom unter der Einwirkung von Temperaturgradienten) in m²/h grd;a _m Stoffübergangszahl (charakterisiert den Stoffstrom unter der Einwirkung von Feuchtigkeitsgradienten) in m²/h - Lu=a _m/a Lykowsche Zahl - Ki _q=q _q ()·R/ _q t dimensionsloser Wärmestrom (Kirpitschew-Zahl);q _q() Wärmestrom durch die Körperoberfläche in kcal/m²; _q Wärmeleitfähigkeit in kcal/m² h grd - Ki _m=q _m ()·R/a _{m 0} u dimensionsloser Stoffstrom;q _m() Stoffstrom durch die Körperoberfläche in kg/m² h; ₀ Wichte des Trockengutes in kg/m³     

Application of the Chapman-Enskog method to the case of a binary two-temperature gas mixture

An interesting property of the flows of a binary mixture of neutral gases for which the molecular mass ratio =m/M1 is that within the limits of the applicability of continuum mechanics the components of the mixture may have different temperatures. The process of establishing the Maxwellian equilibrium state in such a mixture divides into several stages, which are characterized by relaxation times _i which differ in order of magnitude. First the state of the light component reaches equilibrium, then the heavy component, after which equilibrium between the components is established [1]. In the simplest case the relaxation times differ from one another by a factor of ^*.Here the mixture component temperature difference relaxation time _T /, where is the relaxation time for the light component. If 1, 1, so that _T ~1, then for the characteristic hydrodynamic time scale t~1 the relative temperature difference will be of order unity. In the absence of strong external force fields the component velocity difference is negligibly small, since its relaxation time _vt1.In the case of a fully ionized plasma the Chapman-Enskog method is quite easily extended to the case of the two-temperature mixture [3], since the Landau collision integral is used, which decomposes directly with respect to . In the Boltzmann cross collision integral, the quantity appears in the formulas relating the velocities before and after collision, which hinders the decomposition of this integral with respect to , which is necessary for calculating the relaxation terms in the equations for temperatures differing from zero in the Euler approximation [4] (the transport coefficients are calculated considerably more simply, since for their determination it is sufficient to account for only the first (Lorentzian [5]) terms of the decomposition of the cross collision integrals with respect to ). This led to the use in [4] for obtaining the equations of the considered continuum mixture of a specially constructed model kinetic equation (of the Bhatnagar-Krook type) which has an undetermined degree of accuracy.In the following we use the Boltzmann equations to obtain the equations of motion of a two-temperature binary gas mixture in an approximation analogous to that of Navier-Stokes (for convenience we shall term this approximation the Navier-Stokes approximation) to determine the transport coefficients and the relaxation terms of the equations for the temperatures. The equations in the Burnett approximation, and so on, may be obtained similarly, although this derivation is not useful in practice.     

Analysis of suddenly started laminar flow in the entrance region of a circular tube

Suddenly started laminar flow in the entrance region of a circular tube, with constant inlet velocity, is investigated analytically by using integral momentum approach. A closed form solution to the integral momentum equation is obtained by the method of characteristics to determine boundary layer thickness, entrance length, velocity profile, and pressure gradient.Nomenclature M(, , ) a function - N(, , ) a function - p pressure - p* p/1/2U ², dimensionless pressure - Q(, , ) a function - R radius of the tube - r radial distance - Re 2RU/, Reynolds number - t time - U inlet velocity, constant for all time, uniform over the cross section - u velocity in the boundary layer - u* u/U, dimensionless velocity - u ₁ velocity in the inviscid core - x axial distance - y distance perpendicular to the axis of the tube - y* y/R, dimensionless distance perpendicular to the axis - boundary layer thickness - * displacement thickness - /R, dimensionless boundary layer thickness - momentum thickness - absolute viscosity of the fluid - /, kinematic viscosity of the fluid - x/(R Re), dimensionless axial distance - density of the fluid - tU/(R Re), dimensionless time - _w wall shear stress     

A direct formulation of correlations for wall-proximity corrections of hot-wire measurements

Correlations for corrections to hot-wire data for the effects of wall proximity within the viscous sublayer are usually presented in the form u/u = F (y u /). The application of such correlations requires a prior knowledge of the wall shear stress; alternatively, the correlation must be used in an iterative fashion. It is shown in the present note that any such correlation may be recast with no loss of generality in the explicit form u/u _m = f (y u _m/), which is more convenient for use.List of symbols u difference between measured and actual velocities, u _m – u - u _m measured velocity - u shear velocity, - u ⁺ on-dimensional velocity, u/u - y distance from wall - y ⁺ non-dimensional distance from wall, y u / - fluid density - fluid kinematic viscosity - _s wall shear stress     

The Preston tube in adiabatic compressible flow

A new procedure for the reduction of Preston tube data is introduced, based on the van Driest transformation. It appears to give results agreeing with the better calibration experiments, although a significant assumption in its derivation is violated.List of Symbols M _s Mach number sensed by Preston tube - M Friction Mach number (=u/wall sound speed) - R Gas constant - T _w Wall temperature - d Diameter of Preston tube - h Height of effective centre of Preston tube - p Preston tube pressure difference reading - p _i Equivalent incompressible Preston tube reading - p _w Wall pressure - r Recovery factor (=0.896) - u Friction velocity (=[_w/wall density]^1/2) - Empirical constant allowing for departure from Crocco temperature-velocity correlation (=0.975) - Specific heat ratio - Fluid kinematic viscosity - _w Wall shear-stress     

Flow stability of a Grad-model liquid layer on an inclined plane

A study is presented of the flow of stability of a Grad-model liquid layer [1, 2] flowing over an inclined plane under the influence of the gravity force.It is assumed that at every point of the considered material continuum, along with the conventional velocity vector v, there is defined an angular velocity vector , the internal moment stresses are negligibly small, and in the general case the force stress tensor _kj is asymmetric. The model is characterized by the usual Newtonian viscosity , the Newtonian rolling viscosity _r, and the relaxation time = J/4 _r, where J is a scalar constant of the medium with dimensions of moment of inertia per unit mass, is the density. It is assumed that the medium is incompressible, the coefficients , _r, J are constant [2].The exact solution of the equations of motion, corresponding to flow of a layer with a plane surface, coincides with the solution of the Navier-Stokes equations in the case of flow of a layer of Newtonian fluid. The equations for three-dimensional periodic disturbances differ considerably from the corresponding equations for the problem of the flow stability of a layer of a Newtonian medium. It is shown that the Squire theorem is valid for parallel flows of a Grad liquid.The flow stability of the layer with respect to long-wave disturbances is studied using the method of sequential approximations suggested in [3, 4].     

Viscoelastic properties of semidilute poly(methyl methacrylate) solutions

Viscoelastic properties were examined for semidilute solutions of poly(methyl methacrylate) (PMMA) and polystyrene (PS) in chlorinated biphenyl. The number of entanglement per molecule, N, was evaluated from the plateau modulus, G _N . Two time constants, _s and ₁, respectively, characterizing the glass-to-rubber transition and terminal flow regions, were evaluated from the complex modulus and the relaxation modulus. A time constant _k supposedly characterizing the shrink of an extended chain, was evaluated from the relaxation modulus at finite strains. The ratios ₁/ _s and _k / _s were determined solely by N for each polymer species. The ratio ₁/ _s was proportional to N ^4.5 and N ^3.5 for PMMA and PS solutions, respectively. The ratio _k / _s was approximately proportional to N ^2.0 in accord with the prediction of the tube model theory, for either of the polymers. However, the values for PMMA were about four times as large as those for PS. The result is contrary to the expectation from the tube model theory that the viscoelasticity of a polymeric system, with given molecular weight and concentration, is determined if two material constants _s and G _N are known.     

Wandeffekte bei dispers-plastischen Materialien

Zusammenfassung Die Wandgleitgeschwindigkeit von dispers-plastischen Gemischen aus Kaolinpulver und Paraffinöl wird nach der Drei-Spalte-Methode für die Couette-Strömung mit einem Searle-Rheometer ermittelt. Sie steigt zunächst mit zunehmender Schubspannung an, erreicht ein Maximum, fällt mit weiter steigender Schubspannung wieder ab und wird schließlich sogar negativ. Eine negative Wandgleitgeschwindigkeit ist natürlich physikalisch unmöglich. Dispersplastische Gemische aus Kaolinpulver und Paraffinöl zeigen also ein komplizierteres Wandverhalten als reines Wandgleiten.Zur Deutung dieses komplizierten Wandeffektes wird eine Modellvorstellung entwickelt. Wichtig ist hierbei, daß eine zunehmende Wandgleitgeschwindigkeit auftritt, bevor eine starke Scherströmung im Innern des Strömungsfeldes einsetzt. Mit beginnender Scherströmung führen die plättchenförmigen dispersen Teilchen auf Grund von Zusammenstößen seitliche Schwankungsbewegungen um die makroskopisch wahrnehmbaren Bahnkurven aus.Diese Teilchenbewegungen führen zur Zerstörung der zunächst beim Wandgleiten sich ausbildenden Mikrostrukturen an der Wand. Daher kann die Wandgleitgeschwindigkeit trotz steigender Wandschubspannung abnehmen. Die Behinderung der seitlichen Partikelbewegungen an der Wand — die dispersen Teilchen können sich auf der Wand abstützen — führt bei weiter steigender Schergeschwindigkeit im Innern des Strömungsfeldes makroskopisch zu einer Versteifung des Materials in Wandnähe. Damit können negative Werte der sog. Wandgleitgeschwindigkeit — man spricht besser von einer integralen Wandfunktion — sowie bestimmte experimentelle Befunde bei der Druckabhängigkeit und bei der Temperaturabhängigkeit der rheologischen Eigenschaften und des Wandeffektes erklärt werden.Die experimentellen Untersuchungen beschränken sich im wesentlichen auf den Wandeffekt an schwach gekrümmten Wänden in Couette-Spalten, an denen ein Krümmungseinfluß auf den Wandeffekt mit großer Wahrscheinlichkeit vernachlässigbar ist. Die Auswirkung eines Krümmungseinflusses auf die rheometrischen Meßergebnisse wird jedoch diskutiert. Die aus rheometrischen Messungen bestimmbare integrale Wandfunktion liefert im Fall des komplizierten Wandeffektes noch keine vollständige Information über das Wandverhalten.

---

The wall slip velocity of disperse plastic mixtures of kaolin powder and paraffin oil is determined by the so-called three-gap method for Couette flow with a Searle rheometer. At the start it grows with increasing shear stress, reaches a maximum, then decreases with further increases in shear stress and finally becomes negative. From a physical point of view, negative wall-slip-velocities are impossible. Thus it is concluded that disperse plastic mixtures of kaolin powder and paraffin oil show a more complicated wall effect than pure wall slip.In order to explain this complicated wall effect a model of the microstructure near the wall is developed: It is essential that increasing wall slip velocity occurs before the start of shear flow in the interior of the flow field. With shear flow the slab-like disperse particles perform lateral fluctuations around their macroscopically perceptible flow paths. These are caused by collisions between the particles. These lateral particle movements destroy the microstructure at the wall which was built up by pure wall slip. Therefore the wall slip velocity may decrease inspite of increasing wall shear stress. One may then assume a suppression of lateral particle movement at the wall with further increases in the shear in the interior of the flow field which will cause some kind of stiffening of the material near the wall. This assumption can explain the negative values of the so-called slip velocity (which is better termed an integral wall function) as well as some effects in connection with the pressure and temperature dependence of the flow function and integral wall function.The experimental investigations are confined to slowly curved walls in Couette gaps, where an influence of wall curvature on the wall effect may be neglected, but the influence of wall curvature on the wall effect is discussed theoretically. The integral wall function which can be determined from rheometric measurements does not yield complete information on the complicated wall effect.

---

f() Schubspannungsfunktion - Schubspannungsfunktion in Wandnähe - h axiale Erstreckung eines Couette-Spaltes - M _d übertragenes Drehmoment in der Couette-Strömung - R kleinster Krümmungsradius einer Wand an einer Stelle - R _w Radius einer zylindrischen Wand - R _a, R_i Radien von Außen- und Innenzylinder eines Couette-Spaltes - R ₁, R₂, R₃ Radien eines Drei-Spalte-Couette-Systems - R _w1, R_w2 Radien von zwei Rohren - Volumenstrom in einer Rohrströmung - Volumenströme durch zwei verschiedene Rohre bei gleicher Wandschubspannung - v _w (_w) Wandgleitgeschwindigkeit - Winkel zwischen Wandschubspannung und der Richtung, in der die Wand am schwächsten gekrümmt ist - =(R_a/R_i)² quadratisches Radienverhältnis - (_w) Dicke der vom komplizierteren Wandeffekt beeinflußten Wandschicht - Dicke eines Gleitfilms bei Wandgleiten - _w Schubspannungsänderung in der Wandschicht (_w) - f(_w, ) Wandfunktion - Wandabstand - ø _w (_w) integrale Wandfunktion bei vernachlässigbarer Wandkrümmung und vernachlässigbarer Schubspannungsänderung in der Wandschicht (_w) - ø _Couette ( _w, ₂) integrale Wandfunktion der Couette-Strömung - ø _Rohr ( _w, R_w) integrale Wandfunktion der Rohrströmung - ø _Couette ^* ( _w, R₂) experimentell ermittelte Wandfunktion der Couette-Strömung - ø _Rohr ^* ( _w, R_w, R_w2) experimentell ermittelte Wandfunktion der Rohrströmung - ₁, ₂ größter bzw. kleinster Krümmungsradius einer Wand - _w Wandschubspannung - _a, _i Wandschubspannung am Außen- bzw. Innenzylinder eines Couette-Spaltes - ₂ Wandschubspannung in einem Drei-Spalte-Couette-System am mittleren RadiusR ₂ - Schubspannung - Winkelgeschwindigkeitsdifferenz zwischen Außen- und Innenzylinder eines Couette-Spaltes - _I (M_d), _II (M_d), _III(M_d) Winkelgeschwindigkeitsdifferenzen an einem Drei-Spalte-Couette-System als Funktionen des übertragenen Momentes     

Stability of a radially bounded rotating fluid heated from below

The stability of a bounded rotating cylinder of fluid heated from below is treated mathematically under the assumptions of stationary onset and axisymmetry. Critical Rayleigh numbers are computed by Galerkin's method as a function of the Taylor number and cylinder aspect ratio for Taylor numbers,10⁶. The constraining effect of the side walls is shown to decrease with either increasing or increasing radius/height ratios. For>10⁶, most cylinders, excluding extremely tall ones, will appear infinite in horizontal extent as far as stability characteristics are concerned. The form of the motion at onset is discussed in relation to previous work.