Theoretical derivation of tangential velocity profiles in a flat vortex chamber-influence of turbulence and wall friction

Summary As part of a study on the hydrodynamics of a cyclone separator, a theoretical investigation of the flow pattern in a flat box cyclone (vortex chamber) has been carried out. Expressions have been derived for the tangential velocity profile as influenced by internal friction (eddy viscosity) and wall friction. The most important parameter controlling the tangential velocity profile is = –u ₀ R/(v+ ), where u ₀ is the radial velocity at the outer radius R of the cyclone, the kinematic liquid viscosity and is the kinematic eddy viscosity. For values of greater than about 10 the tangential velocity profile is nearly hyperbolic, for smaller than 1 the tangential velocity even decreases towards the centre. It is shown how and also the wall friction coefficient may be obtained from experimental velocity profiles with the aid of suitable graphs. Because of the close relation between eddy viscosity and eddy diffusion, measurements of velocity profiles in flat box cyclones will also provide information on the eddy motion of particles in a cyclone, a motion reducing its separation efficiency.List of symbols A cross-sectional area of cyclone inlet - h height of cyclone - p static pressure in cyclone - p static pressure difference in cyclone between two points on different radius - r radius in cyclone - r ₁ radius of cyclone outlet - R radius of cyclone circumference - u radial velocity in cyclone - u ₀ radial velocity at circumference of flat box cyclone - v tangential velocity - v ₀ tangential velocity at circumference of flat box cyclone - w axial velocity - z axial co-ordinate in cyclone - friction coefficient in flat box cyclone (for definition see § 5) - ₁ value of friction coefficient for ₁<< ₂ - ₂ value of friction coefficient for ₂<<1 - = - ₁ value of for ₁<< ₂ - ₂ value of for ₂<<1 - thickness of laminar boundary layer - =/h - turbulent kinematic viscosity - ratio of z to h - k ratio of height of cyclone to radius R of cyclone - parameter describing velocity profile in cyclone =–u ₀ R/(+) - kinematic viscosity of fluid - density of fluid - ratio of r to R - ₁ value of at outlet of cyclone - ₂ value of at inner radius of cyclone inlet - _w shear stress at cyclone wall - angular momentum in cyclone/angular momentum in cyclone inlet - ₁ value of at = ₁ - ₂ value of at = ₂     

Chaos and integrability in a nonlinear wave equation

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

Representation of the rheological properties of polymer melts in terms of complex fluidity

In dynamic rheological experiments melt behavior is usually expressed in terms of complex viscosity ^* () or complex modulusG ^* (). In contrast, we attempted to use the complex fluidity ^* () = 1/µ ^* () to represent this behavior. The main interest is to simplify the complex-plane diagram and to simplify the determination of fundamental parameters such as the Newtonian viscosity or the parameter of relaxation-time distribution when a Cole-Cole type distribution can be applied. ^* () complex shear viscosity - () real part of the complex viscosity - () imaginary part of the complex viscosity - G ^* () complex shear modulus - G() storage modulus in shear - G() loss modulus in shear - J ^* () complex shear compliance - J() storage compliance in shear - J() loss compliance in shear - shear strain - rate of strain - angular frequency (rad/s) - shear stress - loss angle - ^* () complex shear fluidity - () real part of the complex fluidity - () imaginary part of the complex fluidity - ₀ zero-viscosity - ₀ average relaxation time - h parameter of relaxation-time distribution     

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

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

---

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

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

---

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

The contribution of internal degrees of freedom to the non-Newtonian rheology of model polymer fluids

We report non-equilibrium molecular dynamics simulations of rigid and non-rigid dumbbell fluids to determine the contribution of internal degrees of freedom to strain-rate-dependent shear viscosity. The model adopted for non-rigid molecules is a modification of the finitely extensible nonlinear elastic (FENE) dumbbell commonly used in kinetic theories of polymer solutions. We consider model polymer melts — that is, fluids composed of rigid dumbbells and of FENE dumbbells. We report the steady-state stress tensor and the transient stress response to an applied Couerte strain field for several strain rates. We find that the rheological properties of the rigid and FENE dumbbells are qualitatively and quantitatively similar. (The only exception to this is the zero strain rate shear viscosity.) Except at high strain rates, the average conformation of the FENE dumbbells in a Couette strain field is found to be very similar to that of FENE dumbbells in the absence of strain. The theological properties of the two dumbbell fluids are compared to those of a corresponding fluid of spheres which is shown to be the most non-Newtonian of the three fluids considered.Symbol Definition b dimensionless time constant relating vibration to other forms of motion - F force on center of mass of dumbbell - F _i force on bead i of dumbbell - F force between center of masses of dumbbells and - F _ij force between beads i and j - h vector connecting bead to center of mass of dumbbell - H dimensionless spring constant for dumbbells, in units of / ² - I moment of inertia of dumbbell - J general current induced by applied field - k _B Boltzmann's constant - L angular momentum - m mass of bead, (= m/2) - M mass of dumbbell, g - N number of dumbbells in simulation cell - P translational momentum of center of mass of dumbbell - P pressure tensor - P _xy xy component of pressure tensor - Q separation of beads in dumbbell - Q _eq equilibrium extension of FENE dumbbell and fixed extension of rigid dumbbell - Q ₀ maximum extension of dumbbell - r _ij vector connecting beads i and j - r position vector of center of mass dumbbell - R vector connecting centers of mass of two dumbbells - t time - t ^* dimensionless time, in units of m/ - T ^* dimensionless temperature, in units of /k - u potential energy - u velocity vector of flow field - u _x x component of velocity vector - V volume of simulation cell - X general applied field - strain rate, s^–1 - ^* dimensionless shear rate, in units of /m ² - general transport property - Lennard-Jones potential well depth - friction factor for Gaussian thermostat - shear viscosity, g/cms - ^* dimensionless shear viscosity, in units of m/ ² - ^* dimensionless number density, in units of ^–3 - Lennard-Jones separation of minimum energy - relaxation time of a fluid - angular velocity of dumbbell - orientation angle of dumbbell      

Reflection, transmission and excitation of SH-surface waves by a discontinuity in mass-loading on a semi-infinite elastic medium

The scattering of an SH-wave by a discontinuity in mass-loading on a semi-infinite elastic medium is investigated theoretically. The incident wave is either a plane body wave or a plane SH-surface wave. The problem is reduced to a Wiener-Hopf problem for the scattered wave. In this problem the amplitude spectral density of the particle displacement occurs as unknown function. Special attention is given to the numerical values of the surface wave contributions to the scattered field.Nomenclature x ₁, x ₂, x ₃ Cartesian coordinates - , polar coordinates in x ₁, x ₃-plane - volume mass density - surface mass density of mass-loading - , Lamé constants - U scalar wave function, defined by (2.1) - c _S speed of propagation of uniform shear waves in bulk medium (c _S=(/)^1/2) - angular frequency - t time - k _S wave number of uniform shear waves (k _S=/c _S) - reduced specific acoustic impedance of mass-loading (=k _S /) - k _m wave number of SH-surface wave (k _m=k _S(1+ ²)^1/2) - _1,2,3 partial differentiation with respect to x _1,2,3 - ⁱ angle between x ₃-axis and direction of propagation of incident body wave - ⁱ wave number in horizontal direction of incident body wave ( ⁱ=k _S sin( ⁱ)) - ⁱ wave number in vertical direction of incident body wave ( ⁱ=k _S cos( ⁱ)) - C _1,2 complex amplitude of surface wave excited by a body wave - R reflection factor of surface wave, when surface wave is incident - T transmission factor of surface wave, when surface wave is incident - S particle displacement vector The research presented in this paper has been carried out with partial financial support from the Delfts Hogeschoolfonds.     

Free convection from a vertical plate with concentrated and distributed thermal sources

An analysis is developed for the laminar free convection from a vertical plate with uniformly distributed wall heat flux and a concentrated line thermal source embedded at the leading edge. We introduce a parameter=(1 +Q _L/Q_w)^–1=(1 + Ra_L/Ra_w)^–1 to describe the relative strength of the two thermal sources; and propose a unified buoyancy parameter=( Ra_L+ Ra_w)^1/5 with=1/(1 +Pr ^–1) to properly scale the dependent and independent variables. The variables are so defined that the resulting nonsimilar boundary-layer equations can describe exactly the buoyancy-induced flow from the dual sources with any relative strength to fluids of any Prandtl number from very small values to infinity. These nonsimilar equations are readily reducible to the self-similar equations of an adiabatic wall plume for=0, and to those of free convection from uniform flux plate for=1. Rigorous finite-difference solutions for fluids of Pr from 0.001 to are obtained over the entire range of from 0 to 1. The effects of both relative source strength and Prandtl number on the velocity profiles, temperature profiles, and the variations of wall temperature, are clearly illustrated.

---

Freie Konvektion an einer vertikalen Platte mit einer konzentrierten und einer gleichmäßig verteilten Wärmequelle

Zusammenfassung Für die freie Konvektion an einer vertikalen Platte mit einer gleichmäßig verteilten Wandwärmestromdichte und einer in der Vorderkante eingebetteten linienförmigen Wärmequelle wird eine Berechnungsmethode entwickelt. Zur Beschreibung der relativen Stärke der beiden Wärmequellen führen wir einen Parameter=(1 + Q_L/Q_w)^–1=(1 + Ra_L/Ra_w)^–1 ein und schlagen einen vereinheitlichten Auftriebsparameter=( Ra _L+ Ra _w)^1/5 mit=1/(1 +Pr ^–1 für die Skalierung der abhängigen und unabhängigen Variablen vor. Die Variablen werden so definiert, daß mit den sich ergebenden unabhängigen Grenzschichtgleichungen die von den beiden Wärmequellen beliebiger Stärke verursachte Auftriebsströmung von Fluiden beliebiger Prandtl-Zahl genau beschrieben werden kann. Diese unabhängigen Gleichungen können ohne weiteres auf die selbstähnlichen Gleichungen für den Fall einer lokalen Wärmezufuhr an einer sonst adiabatischen Wand für=0 und jenen der freien konvektion an einer Platte mit einheitlichem Wärmestrom für=1 zurückgeführt werden. Für Fluide mit der Prandtl-Zahl zwischen 0,001 und Unendlich werden nach der strengen finite Differenzen-Methode Lösungen im Bereich von zwischen 0 und 1 erhalten. Der jeweilige Einfluß der relativen Quellenstärke und der Prandtl-Zahl auf die Geschwindigkeits- und Temperaturprofile sowie die Veränderung der Wandtemperatur werden deutlich dargestellt.

---

Nomenclature C _f friction coefficient - C _p specific heat - f reduced stream function - g gravitational acceleration - k thermal conductivity - L width of the plate - Nu local Nusselt number - Pr Prandtl number - q _w wall heat flux - Q _L heat generated by the line source - Q _w heat released by the uniform-flux wall from 0 tox, q _w Lx - Ra _L local Rayleigh number, g T _L ^* x ³/( ) - Ra _w local Rayleigh number,g T _w ^* w ³/( ) - T fluid temperature - T temperature of ambient fluid - T _L ^* characteristic temperature of the line source,Q _{L/(C p} L) - T _w ^* characteristic temperature of the uniform flux wall, =q _w x/k=Q _w /(C _p L) - u velocity component in then-direction - U₀ dimensionless velocity,u/(/x) Ra _L ^2/5 - U ₁ dimensionless velocity,u/(/x) Ra _w ^2/5 - velocity component in they-direction - x coordinate parallel to the plate - y coordinate normal to the plate - thermal diffusivity - thermal expansion coefficient - pseudo-similarity variable,(y/x) - dimensionless temperature, (T–T )/(T _L ^* +T _w ^* ) - ₀ dimensionless temperature, (Ra_l)^1/5 (T–T )/T _L ^* - ₁ dimensionless temperature, (Ra_w)Ra_w)^1/5 (T–T )/T _w ^* - (Ra _L+Ra_w)^1/5 - kinematic viscosity - (1 +Ra _L/Ra_w)^–1=(1 +T _L ^* /T _w ^* )^–1=(1 + Q_L/Q_w)^–1 - density - Pr/(1 +Pr) - _w wall shear stress - stream function     

The rheology of polymeric liquid crystals

The molecular theory of Doi has been used as a framework to characterize the rheological behavior of polymeric liquid crystals at the low deformation rates for which it was derived, and an appropriate extension for high deformation rates is presented. The essential physics behind the Doi formulation has, however, been retained in its entirety. The resulting four-parameter equation enables prediction of the shearing behavior at low and high deformation rates, of the stress in extensional flows, of the isotropic-anisotropic phase transition and of the molecular orientation. Extensional data over nearly three decades of elongation rate (10^–2–10¹) and shearing data over six decades of shear rate (10^–2–10⁴) have been correlated using this analysis. Experimental data are presented for both homogeneous and inhomogeneous shearing stress fields. For the latter, a 20-fold range of capillary tube diameters has been employed and no effects of system geometry or the inhomogeneity of the flow-field are observed. Such an independence of the rheological properties from these effects does not occur for low molecular weight liquid crystals and this is, perhaps, the first time this has been reported for polymeric lyotropic liquid crystals; the physical basis for this major difference is discussed briefly. A Semi-empirical constant in eq. (18), N/m² - c rod concentration, rods/m³ - c ^* critical rod concentration at which the isotropic phase becomes unstable, rods/m³ - C interaction potential in the Doi theory defined in eq. (3) - d rod diameter, m - D semi-empirical constant in eq. (19), s^–1 - D _r lumped rotational diffusivity defined in eq. (4), s^–1 - rotational diffusivity of rods in a concentrated (liquid crystalline) system, s^–1 - D _ro rotational diffusivity of a dilute solution of rods, s^–1 - f distribution function defining rod orientation - F tensorial term in the Doi theory defined in eq. (7) (or eq. (19)), s^–1 - G tensorial term in the Doi theory defined in eq. (8) - K _B Boltzmann constant, 1.38 × 10^–23 J/K-molecule - L rod length, m - S scalar order parameter - S tensor order parameter defined in eq. (5) - t time, s - T absolute temperature, K - u unit vector describing the orientation of an individual rod - rate of change ofu due to macroscopic flow, s^–1 - v fluid velocity vector, m/s - v velocity gradient tensor defined in eq. (9), s^–1 - V mean field (aligning) potential defined in eq. (2) - x coordinate direction, m - Kronecker delta (= 0 if = 1 if = ) - _r ratio of viscosity of suspension to that of the solvent at the same shear stress - _s solvent viscosity, Pa · s - ^* viscosity at the critical concentrationc ^*, Pa · s - v ₁, v₂ numerical factors in eqs. (3) and (4), respectively - deviatoric stress tensor, N/m² - volume fraction of rods - ⁰ constant in eq. (16) - ^* volume fraction of rods at the critical concentrationc ^* - average over the distribution functionf(u, t) (= d ²u f(u, t)) - gradient operator - d ²u integral over the surface of the sphere (|u| = 1)     

Transient flow towards a well in an aquifer including the effect of fluid inertia

Transient propagation of weak pressure perturbations in a homogeneous, isotropic, fluid saturated aquifer has been studied. A damped wave equation for the pressure in the aquifer is derived using the macroscopic, volume averaged, mass conservation and momentum equations. The equation is applied to the case of a well in a closed aquifer and analytical solutions are obtained to two different flow cases. It is shown that the radius of influence propagates with a finite velocity. The results show that the effect of fluid inertia could be of importance where transient flow in porous media is studied.List of symbols b Thickness of the aquifer, m - c ₀ Wave velocity, m/s - k Permeability of the porous medium, m² - n Porosity of the porous medium - p( ,t) Pressure, N/m² - Q Volume flux, m³/s - r Radial coordinate, m - r _w Radius of the well, m - s Transform variable - S Storativity of the aquifer - S _d(r, t) Drawdown, m - t Time, s - T Transmissivity of the aquifer, m²/s - ( ,t) Velocity of the fluid, m/s - Coordinate vector, m - z Vertical coordinate, m - Coefficient of compressibility, m²/N - Coefficient of fluid compressibility, m²/N - Relaxation time, s - (r, t) Hydraulic potential, m - Dynamic viscosity of the fluid, Ns/m² - Dimensionless radius - Density of the fluid, Ns²/m⁴ - (, ) Dimensionless drawdown - Dimensionless time - , x Dummy variables - ₀, ₁ Auxilary functions     

Stochastic response of structures with small geometric imperfections

Summary A probabilistic model of the geometric imperfections of a real structure is proposed, in order to provide a general theory of the stochastic response of structures in presence of small random deviations from the perfect scheme. The main statistical measures of the stochastic response are derived and an application to the study of a particular conservative elastic system is developed.

---

Sommario Si propone una teoria generale della risposta probabilistica di strutture, in presenza di piccole deviazioni aleatorie dei dati iniziali rispetto allo schema geometrico perfetto. Si deducono le principali proprietà statistiche della risposta della struttura a sollecitazioni esterne deterministiche, e si sviluppa una applicazione riguardante il comportamento aleatorio di un particolare sistema elastico conservativo.

---

List of symbols element of the sample space of events - _kn random variables modelling the structural imperfections - P(o) probability density of random variables - random imperfection of the unloaded structure - u additional displacement of the loaded structure - u_o deterministic fundamental solution for the perfect structure - difference between the additional displacement of the loaded structure and the deterministic fundamental solution for the perfect structure - V₁=u₁ buckling mode of the perfect structure - _i intrinsic coordinates of the structure - suitable measure of the magnitude of the random imperfections - scalar geometric variable representing the internal product - random imperfection divided by - single scalar variable denoting the magnitude of the prescribed loads - potential energy of the structure - potential energy of the perfect structure - difference between and - _c lowest critical load - _s real local maximum for the magnitude of the prescribed loads - _c divided by _S - E{} expected value of a random variable - ² variance of a random variable - , random variables defined by Eq. (21)     

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     

Saint-Venant's principle in linear two-dimensional elasticity for non-striplike domains

Let be a simply connected domain in the x ₁-x ₂ plane which lies within the strip 0<x ₂, is a simple closed piecewise smooth curve. Let l= [(x ₁, x ₂): (x ₁, x ₂) and x ₁>0], _l = [(x ₁ x ₂): (x ₁ ,x ₂) and x ₁>1>0].Suppose that a two-dimensional homogeneous isotropic elastic body occupies , that a self-equilibrated stress loading is applied to - l, and that l is stress-free. Knowles [2] and Flavin [6] showed that the elastic energy in _l decays exponentially with respect to l with an exponential decay constant of the form k/b, where k is a universal constant. It is shown here that a decay constant of the form c/ may be obtained where c is a universal constant and is a characteristic dimension of , which is more appropriate than b for general non-striplike domains. In addition, an appropriate decay theorem is obtained for coil-like domains.     

Thermal free convection from a solid sphere

Summary Thermal free convection from a sphere has been studied by melting solid benzene spheres in excess liquid benzene (Pr=8,3; 10⁸<Gr<10⁹). Overall heat transfer as well as local heat transfer were investigated. For the effect of cold liquid produced by the melting a correction has been applied. Results are compared with those obtained by other workers who used alternative experimental methods.Nomenclature coefficient of heat transfer - d characteristic length, here diameter of sphere - thermal conductivity - g acceleration of free fall - cubic expansion coefficient - T temperature difference between wall and fluid at infinity - kinematic viscosity - density - c specific heat capacity - a thermal diffusivity (=/c) - D diffusion coefficient - Nu dimensionless Nusselt number (=d/) - Nu* the analogous number for mass transfer (=kd/D) - mean value of Nusselt number - Gr dimensionless Grashof number (=gd ³T/ ²) - Gr* the analogous number for mass transfer (=gd ³x/ ²) - Pr dimensionless Prandtl number (=/a) - Sc dimensionless Schmidt number (=/D)     

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

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

Free convection on an arbitrarily inclined plate with uniform surface heat flux

This paper proposed a proper inclination parameter and transformation variables for the analysis of free convection from an inclined plate with uniform surface heat flux to fluids of any Prandtl number. Very accurate numerical results and a simple correlation equation are obtained for arbitrary inclination from the horizontal to the vertical and for 0.001 Pr. Maximum deviation between the correlated and calculated data is less than 1.2%.

---

Freie Konvektion an einer beliebig geneigten Platte mit erheblicher Wärmestromdichte an der Oberfläche

Zusammenfassung Für die Berechnung von freier Konvektion von Fluiden mit beliebiger Prandtl-Zahl an einer geneigten Platte mit einheitlicher Wärmestromdichte an der Oberfläche werden ein zweckmäßiger Neigungsparameter und Transformationsvariablen eingeführt. Sehr genaue numerische Ergebnisse und eine einfache Korrelationsgleichung wurden für beliebige Neigungen zwischen der Horizontalen und der Vertikalen und für 0.001Pr erhalten. Die größte Abweichung zwischen Korrelations- und berechneten Daten liegt bei weniger als 1.2%.

---

Nomenclature f reduced stream function - g gravitational acceleration - h local heat transfer coefficient - k thermal conductivity - Nu local Nusselt number - p static pressure difference - Pr Prandtl number - q _w wall heat flux - Ra* modified local Rayleigh number,g(q _w x/k)x ³/ - T fluid temperature - T temperature of ambient fluid - u velocity component inx-direction - v velocity component iny-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - (Ra* |sin|)^1/5/( Ra* cos)^1/6 - ( Ra* cos)^1/6+( Ra*|sin|)^1/5 - (y/x) - dimensionless temperature, (T–T )/(q _w x/k) - kinematic viscosity - [1+( Ra* cos)^1/6/( Ra*|sin|)^1/5]^–1 - density of fluid - Pr/(1+Pr) - _w wall shear stress - angle of inclination measured from the horizontal - stream function - dimensionless static pressure difference, p x ²/ ⁴     

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     

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     

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     

On Perturbative Expansions to the Stochastic Flow Problem

When analyzing stochastic steady flow, the hydraulic conductivity naturally appears logarithmically. Often the log conductivity is represented as the sum of an average plus a stochastic fluctuation. To make the problem tractable, the log conductivity fluctuation, f, about the mean log conductivity, lnK _G, is assumed to have finite variance, _f ². Historically, perturbation schemes have involved the assumption that _f ²<1. Here it is shown that _f may not be the most judicious choice of perturbation parameters for steady flow. Instead, we posit that the variance of the gradient of the conductivity fluctuation, _f ², is more appropriate hoice. By solving the problem withthis parameter and studying the solution, this conjecture can be refined and an even more appropriate perturbation parameter, , defined. Since the processes f and f can often be considered independent, further assumptions on f are necessary. In particular, when the two point correlation function for the conductivity is assumed to be exponential or Gaussian, it is possible to estimate the magnitude of f in terms of _f and various length scales. The ratio of the integral scale in the main direction of flow ( _x ) to the total domain length (L^*), _x ²=_x/L^*, plays an important role in the convergence of the perturbation scheme. For _x smaller than a critical value _c, _x < _c, the scheme's perturbation parameter is =_f/_x for one- dimensional flow, and =_f/_x ² for two-dimensional flow with mean flow in the x direction. For _x > _c, the parameter =_f/_x ³ may be thought as the perturbation parameter for two-dimensional flow. The shape of the log conductivity fluctuation two point correlation function, and boundary conditions influence the convergence of the perturbation scheme.