A three-parameter model describing the behaviour of a viscoelastic liquid in a tangential annular flow

A three-parameter model describing the shear rate-shear stress relation of viscoelastic liquids and in which each parameter has a physical significance, is applied to a tangential annular flow in order to calculate the velocity profile and the shear rate distribution. Experiments were carried out with a 5000 wppm aqueous solution of polyacrylamide and different types of rheometers. In a shear-rate range of seven decades (5 10^–3 s^–1 < < 1.2 10⁵ s^–1) a good agreement is obtained between apparent viscosities calculated with our model and those measured with three different types of rheometers, i.e. Couette rheometers, a cone-and-plate rheogoniometer and a capillary tube rheometer. a physical quantity defined by:a = {1 – ( / ₀)}/ ₀ (Pa^–1) - C constant of integration (1) - r distancer from the center (m) - r ₁,r ₂ radius of the inner and outer cylinder (m) - v _r local tangential velocity at a distancer from the center (v _r = _r r) (m s^–1) - v ₂ local tangential velocity at a distancer ₂ from the center (m s^–1) - shear rate (s^–1) - local shear rate (s^–1) - ₁ wall shear rate at the inner cylinder (s^–1) - dynamic viscosity (Pa s) - _a apparent viscosity (_a = / ) (Pa s) - _a1 apparent viscosity at the inner cylinder (Pa s) - ₀ zero-shear viscosity (Pa s) - infinite-shear viscosity (Pa s) - shear stress (Pa) - _r local shear stress at a distancer from the center (Pa) - ₀ yield stress (Pa) - ₁, ₂ wall shear-stress at the inner and outer cylinder (Pa) - _r local angular velocity (s^–1) - ₂ angular velocity of the outer cylinder (s^–1)     

Two-phase flow in heterogeneous porous media III: Laboratory experiments for flow parallel to a stratified system

Two-phase flow in stratified porous media is a problem of central importance in the study of oil recovery processes. In general, these flows are parallel to the stratifications, and it is this type of flow that we have investigated experimentally and theoretically in this study. The experiments were performed with a two-layer model of a stratified porous medium. The individual strata were composed of Aerolith-10, an artificial: sintered porous medium, and Berea sandstone, a natural porous medium reputed to be relatively homogeneous. Waterflooding experiments were performed in which the saturation field was measured by gamma-ray absorption. Data were obtained at 150 points distributed evenly over a flow domain of 0.1 × 0.6 m. The slabs of Aerolith-10 and Berea sandstone were of equal thickness, i.e. 5 centimeters thick. An intensive experimental study was carried out in order to accurately characterize the individual strata; however, this effort was hampered by both local heterogeneities and large-scale heterogeneities.The theoretical analysis of the waterflooding experiments was based on the method of large-scale averaging and the large-scale closure problem. The latter provides a precise method of discussing the crossflow phenomena, and it illustrates exactly how the crossflow influences the theoretical prediction of the large-scale permeability tensor. The theoretical analysis was restricted to the quasi-static theory of Quintard and Whitaker (1988), however, the dynamic effects described in Part I (Quintard and Whitaker 1990a) are discussed in terms of their influence on the crossflow.Roman Letters A interfacial area between the -region and the -region contained within V, m² - a vector that maps onto , m - b vector that maps onto , m - b vector that maps onto , m - B second order tensor that maps onto , m² - C second order tensor that maps onto , m² - E energy of the gamma emitter, keV - f fractional flow of the -phase - g gravitational vector, m/s² - h characteristic length of the large-scale averaging volume, m - H height of the stratified porous medium , m - i unit base vector in the x-direction - K local volume-averaged single-phase permeability, m² - K - {K}, large-scale spatial deviation permeability - { K} large-scale volume-averaged single-phase permeability, m² - K ^* large-scale single-phase permeability, m² - K ^** equivalent large-scale single-phase permeability, m² - K local volume-averaged -phase permeability in the -region, m² - K local volume-averaged -phase permeability in the -region, m² - K - {K } , large-scale spatial deviation for the -phase permeability, m² - K ^* large-scale permeability for the -phase, m² - l thickness of the porous medium, m - l characteristic length for the -region, m - l characteristic length for the -region, m - L length of the experimental porous medium, m - characteristic length for large-scale averaged quantities, m - n outward unit normal vector for the -region - n outward unit normal vector for the -region - n unit normal vector pointing from the -region toward the -region (n = - n ) - N number of photons - p pressure in the -phase, N/m² - p ₀ reference pressure in the -phase, N/m² - local volume-averaged intrinsic phase average pressure in the -phase, N/m² - large-scale volume-averaged pressure of the -phase, N/m² - large-scale intrinsic phase average pressure in the capillary region of the -phase, N/m² - - , large-scale spatial deviation for the -phase pressure, N/m² - p_c , capillary pressure, N/m² - p _c capillary pressure in the -region, N/m² - p capillary pressure in the -region, N/m² - {p _c } ^c large-scale capillary pressure, N/m² - q -phase velocity at the entrance of the porous medium, m/s - q -phase velocity at the entrance of the porous medium, m/s - S_wi irreducible water saturation - S /, local volume-averaged saturation for the -phase - S _i initial saturation for the -phase - S _r residual saturation for the -phase - S ^* { }^*/}^*, large-scale average saturation for the -phase - S saturation for the -phase in the -region - S saturation for the -phase in the -region - t time, s - v -phase velocity vector, m/s - v local volume-averaged phase average velocity for the -phase, m/s - {v } large-scale averaged velocity for the -phase, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - V large-scale averaging volume, m³ - y position vector relative to the centroid of the large-scale averaging volume, m - {y}^c large-scale average of y over the capillary region, m Greek Letters local porosity - local porosity in the -region - local porosity in the -region - local volume fraction for the -phase - local volume fraction for the -phase in the -region - local volume fraction for the -phase in the -region - {}^* { }^*+{ }^*, large-scale spatial average volume fraction - { }^* large-scale spatial average volume fraction for the -phase - mass density of the -phase, kg/m³ - mass density of the -phase, kg/m³ - viscosity of the -phase, N s/m² - viscosity of the -phase, Ns/m² - V /V , volume fraction of the -region ( + =1) - V /V , volume fraction of the -region ( + =1) - attenuation coefficient to gamma-rays, m^-1 - -      

Some characteristic features of diffusion of a magnetic field into a moving conductor

The study of the diffusion of a magnetic field into a moving conductor is of interest in connection with the production of ultra-high-strength magnetic fields by rapid compression of conducting shells [1,2]. In [3,4] it is shown that when a magnetic field in a plane slit is compressed at constant velocity, the entire flux enters the conductor. In the present paper we formulate a general result concerning the conservation of the sum current in the cavity and conductor for arbitrary motion of the latter. We also consider a special case of conductor motion when the flux in the cavity remains constant despite the finite conductivity of the material bounding the magnetic field.Notation ₁, * flux which has diffused into the conductor - ₂ flux in the cavity - ₀ sum flux - r radius - r* cavity boundary - thickness of the skin layer - (r) delta function of r - t time - q intensity of the fluid sink - v velocity - flux which has diffused to a depth larger than r - x self-similar variable - dimensionless fraction of the flux which has diffused to a depth larger than r - * fraction of the flux which has diffused into the conductor - a conductivity - c electrodynamic constant - R_m magnetic Reynolds number - dimensionless parameter     

The theory of a vibrating-rod viscometer

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

Parameter determination for thermodynamical models of the pseudoelastic behaviour of shape memory alloy by resistivity measurements and infrared thermography

Two thermodynamical models of pseudoelastic behaviour of shape memory alloys have been formulated. The first corresponds to the ideal reversible case. The second takes into account the hysteresis loop characteristic of this shape memory alloys.Two totally independent techniques are used during a loading-unloading tensile test to determine the whole set of model parameters, namely resistivity and infrared thermography measurements. In the ideal case, there is no difficulty in identifying parameters.Infrared thermography measurements are well adapted for observing the phase transformation thermal effects.Notations 1 austenite 2 martensite - ₍₎ Macroscopic infinitesimal strain tensor of phase - ₍₂₎ ^f Traceless strain tensor associated with the formation of martensite phase - Macroscopic infiniesimal strain tensor - Macroscopic infinitesimal strain tensor deviator - _f Trace - Equivalent strain - ^pe Macroscopic pseudoelastic strain tensor - x Distortion due to parent (austenite =1)product (martensite =2) phase transformation (traceless symmetric second order tensor) - M Total mass of a system - M() Total mass of phase - V Total volume of a system - V() Total volume of phase - z=M(2)/M Weight fraction of martensite - 1-z=M(1)/M Weight fraction of austenite - u ₀ ^* () Specific internal energy of phase (=1,2) - s ₀ ^* () Specific internal entropy of phase - Specific configurational energy - Specific configurational entropy - ₀ ^f (T) Driving force for temperature-induced martensitic transformation at stress free state ( ₀ ^f T) = T ^*–Ts ^*) - Kirchhoff stress tensor - Kirchhoff stress tensor deviator - Equivalent stress - Cauchy stress tensor - Mass density - K Bulk moduli (K ₀=K) - L Elastic moduli tensor (order 4) - E Young modulus - Energetic shear (₀ = ) - Poisson coefficient - M _s ^o (M _F ^o ) Martensite start (finish) temperature at stress free state - A _s ^o (A _F ^o ) Austenite start (finish) temperature at stress free state - C _v Specific heat at constant volume - k Conductivity - Pseudoelastic strain obtained in tensile test after complete phase transformation (AM) (unidimensional test) - ₀ Thermal expansion tensor - r Resistivity - 1MPa 10⁶ N/m ² - ⁽⁾ Specific free energy of phase - _n Specific free energy at non equilibrium (R model) - _n ^eq Specific free energy at equilibrium (R model) - _n ^v Volumic part of ^eq - Specific free energy at non equilibrium (R _L model) - _conf Specific coherency energy (R _L model) - _c Specific free energy at constrained equilibria (R _L model) - _it (T) Coherency term (R _L model)     

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.     

The theory of a vibrating-rod densimeter

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

On flow and stability of a Newtonian fluid past a rotating plane

An exact solution is given for the steady flow of a Newtonian fluid occupying the halfspace past the plane z=0 uniformly rotating about a fixed normal axis (Oz). This solution is obtained in a velocity field of the form considered by Berker [2] and can be deduced as a limiting case, as h+, of the solution to the problem relative to the strip 0zh imposing at z=h either the adherence boundary conditions or the free surface conditions. Furthermore, the stability of this flow, subject to periodic disturbances of finite amplitude, is studied using the energy method and the result is compared with those corresponding to stability of flows in the strip 0zh.

---

Sommario In questa nota si mostra che-oltre alla calssica soluzione di von Karman [1] — esiste, per opportuni valori del gradiente di pressione all'infinito, una soluzione esatta per il moto stazionario di un fluido Newtoniano posto nel semispazio limitato dal piano z=0 uniformemente rotante attorno ad un asse ad esso perpendicolare (Oz). Tale soluzione, ottenuta sulla scia del lavoro di Berker [2], si può dedurre anche come limite, per h+, della soluzione del problema relativo alla striscia 0zh quando sul piano z=h si assegnano o le condizioni di aderenza o le condizioni di frontiera libera. Si studia poi la stabilità di tale moto rispetto a perturbazioni spazialmente periodiche di ampiezza finita col metodo dell'energia e si confronta il risultato ottenuto con quelli relativi alla stabilità dei moti nella striscia 0zh.

     

A distributed electrical analog for transient diffusion-techniques of model fabrication and calibration

Simulation of transient two-dimensional diffusion by means of a distributed electrical analog is discussed. After present techniques of analog model construction and calibration are reviewed, an improved calibration technique is presented and a convenient method of analog fabrication, not previously reported, is described. The proposed new method allows complete access to any point on the analog model during a test. Frequency response and step response measurements indicate that an adequate simulation is provided by this particular type of analog model.

---

Zusammenfassung Die Simulierung eines nichtstationären, zweidimensionalen Diffusionsvorganges mittels eines kontinuierlichen, elektro-thermischen Analogapparates wird besprochen. Eine Übersicht der gegenwärtigen Methoden für die Konstruktion und Kalibrierung von elektro-thermischen Analogmodellen wird gegeben. Ein verbessertes Verfahren für die Kalibrierung und eine handliche Fertigungsmethode von Analogmodellen, die noch nicht in der Literatur beschrieben wurden, werden dargestellt. Das vorgeschlagene neue Verfahren gestattet vollständigen Zugang zu jedem Punkt im Analogmodell während des Experiments. Meßbeobachtungen des periodischen Frequenzverhaltens und des nichtstationären Verhaltens zeigen, daß dieses spezielle Analogmodell eine ausreichende Simulierung des Diffusionsvorganges gestattet.

---

Nomenclature A ₀ Potential amplitude atx=0 (see Fig. 5) - A _x Potential amplitude at locationx (see Fig. 5) - A _L Potential amplitude atx=L (see Fig. 5) - C Capacitance per unit area - j Frequency - L Characteristic length - R Resistance per square - t Time - x Coordinate - X Dimensionless distance,x/L - y Coordinate - Y Dimensionless distance,y/L - Diffusivity - _x Phase angle at locationx (see Fig. 5) - _L Phase angle atx=L (see Fig. 5) - _d Thickness of dielectric sheet - _r Thickness of resistance sheet - Dielectric constant - Resistivity - Dimensionless time,t/L ² - Potential - ₁ Reference potential - ₂ Reference potential - Dimensionless potential, ( – ₁)/( ₂ – ₁) - Angular frequency, 2f - Dimensionless frequency,L ²/ The investigation was performed while the first author was Visiting Associate Professor at Purdue University during 1967/68.     

Heat transfer from laminar flow in ducts with elliptic cross-section

Summary The cooling of a hot fluid in laminar Newtonian flow through cooled elliptic tubes has been calculated theoretically. Numerical data have been computed for the two values 1.25 and 4 of the axial ratio of the elliptic cross-section . For =1.25 the influence of non-zero thermal resistance between outmost fluid layer and isothermal surroundings has also been investigated. Special attention has been given to the distribution of heat flux around the perimeter; when increases the flux varies more with the position at the circumference. This positional dependence becomes less pronounced, however, as the (position-independent) thermal resistance of the wall increases.Flattening of the conduit, while maintaining its cross-sectional area constant, improves the cooling. Comparison with rectangular pipes shows that this improvement is not as marked with elliptic as with rectangular pipes.Nomenclature A _k =A _{m, n} coefficients of expansion (6) - a, b half-axes of ellipse, b<a - a _p =a _{r, s} coefficients of representation (V) - D hydraulic diameter, = 4S/P; S = cross-sectional area, P = perimeter - D _e equivalent diameter, according to (13) - n coordinate (outward) normal to the tube wall - T temperature of fluid - T _i temperature of fluid at the inlet - T _s temperature of surroundings - v ₀ mean velocity of fluid - v _z longitudinal velocity of fluid - x, y carthesian coordinates coinciding with axes of ellipse - z coordinate in flow direction - , dimensionless half-axes of ellipse, =a/D and =b/D - _t heat transfer coefficient from fluid at bulk temperature to surroundings; equation (11) - _w heat transfer coefficient at the wall; equation (3) - axial ratio of ellipse, = a/b = / - , , , dimensionless coordinates; =x/D, =y/D, =z/D, =n/D - dimensionless temperature, = (T–T _s)/(T _i–T _s) - ₀ cup-mixing mean value of ; equation (10) - thermal conductivity of fluid - _m,n = _k eigenvalue - c volumetric heat capacity of fluid - _{m, n} = _k = _k eigenfunction; equations (6) and (I) - Nu total Nusselt number, = _t D/ - Nusselt number at large distance from the inlet - Nu _w wall Nusselt number, = _w D/, based on _w - Pé Péclet number, = ₀ Dc/     

Finite difference solution of unsteady natural convection boundary layer flow over an inclined plate with variable surface temperature

The unsteady natural convection boundary layer flow over a semi-infinite inclined plate is considered with the wall temperatureT _w ,(x) (=T +ax ⁿ )varying as the power of the axial coordinate. The governing equations are solved by an implicit finite difference scheme of Crank-Nicolson type. Numerical results are obtained for different values of Prandtl number, Grashof number and n for different angles of inclination. The steadystate velocity and temperature profiles, local and average skin frictions and Nusselt numbers are shown graphically. The effects of the angle of inclination and exponent n on velocity and temperature profiles, skin friction and Nusselt number have been discussed. The velocity, temperature and Nusselt number of the present study are compared with the available results and a good agreement is found to exist between the two.

---

Differenzlösung für nichtstationäre natürliche Grenzschicht-Konvektionsströmung an einer geneigten Platte mit veränderlicher Oberflächentemperatur

Zusammenfassung Die nichtstationäre natürliche Grenzschicht-Konvektionsströmung an einer halbunendlichen geneigten Platte wird unter Zugrundelegung der GesetzmäßigkeitT _w ,(x) (=T +ax ⁿ für die Wandtemperatur als Funktion der Achsialkoordinate untersucht, und zwar mit Hilfe eines impliziten Differenzverfahrens vom Crank-Nicolson Typ und bei Variation der Prandtl- und Grashof-Zahlen, des Exponenten n und des Neigungswinkels. Graphisch dargestellt sind die Geschwindigkeits-und Temperaturprofile im stationären Zustand, die örtlichen und gemittelten Reibungsbeiwerte und der Nusselt-Zahlen. Der Einfluß des Neigungswinkels und des Exponenten n auf diese Größe wird diskutiert. Im Vergleich mit den Ergebnissen aus anderen Arbeiten konnte gute Übereinstimmung festgestellt werden.

---

Nomenclature g acceleration due to gravity - Gr _L Grashof number at x=L - L length of the plate - n exponent in the power law variation of the wall temperature - Nu _x local Nusselt number - Nu _X dimensionless local Nusselt number - average Nusselt number - dimensionless average Nusselt number - p pressure - Pr Prandtl number - t time - t dimensionless time - T temperature - Tw temperature on the plate - T dimensionless temperature - u x-velocity component - U dimensionlessX-velocity component - v y-velocity component - V dimensionlessY-velocity component - x spatial coordinate along the plate - X dimensionless spatial coordinate along the plate - y spatial coordinate normal to the plate - Y dimensionless spatial coordinate normal to the plate Greek symbols thermal diffusivity - ß volumetric coefficient of thermal expansion - t dimensionless time-step - X dimensionless finite difference grid spacing in theX-direction - Y dimensionless finite difference grid spacing in theY-direction - angle of inclination of plate with horizontal - kinematic viscosity - density - _x local skin friction - _X dimensionless local skin friction - average skin friction - dimensionless average skin friction     

Effects of MHD free convection and mass transfer on the flow past an oscillating infinite coaxial vertical circular cylinder

The effects of MHD free convection and mass transfer are taken into account on the flow past oscillating infinite coaxial vertical circular cylinder. The analytical expressions for velocity, temperature and concentration of the fluid are obtained by using perturbation technique.

---

Einwirkungen von freier MHD-Konvektion und Stoffübertragung auf eine Strömung nach einem schwingenden unendlichen koaxialen vertikalen Zylinder

Zusammenfassung Die Einwirkungen der freien MHD-Konvektion und Stoffübertragung auf eine Strömung nach einem schwingenden, unendlichen, koaxialen, vertikalen Zylinder wurden untersucht. Die analytischen Ausdrücke der Geschwindigkeit, Temperatur und Fluidkonzentration sind durch die Perturbationstechnik erhalten worden.

---

Nomenclature C _p Specific heat at constant temperature - C the species concentration near the circular cylinder - C _w the species concentration of the circular cylinder - C the species concentration of the fluid at infinite - * dimensionless species concentration - D chemical molecular diffusivity - g acceleration due to gravity - Gr Grashof number - Gm modified Grashof number - K thermal conductivity - Pr Prandtl number - r _a ,r _b radius of inner and outer cylinder - a, b dimensionless inner and outer radius - r,r coordinate and dimensionless coordinate normal to the circular cylinder - Sc Schmidt number - t time - t dimensionless time - T temperature of the fluid near the circular cylinder - T _w temperature of the circular cylinder - T temperature of the fluid at infinite - u velocity of the fluid - u dimensionless velocity of the fluid - U ₀ reference velocity - z,z coordinate and dimensionless coordinate along the circular cylinder - coefficient of volume expansion - * coefficient of thermal expansion with concentration - dimensionless temperature - H ₀ magnetic field intensity - coefficient of viscosity - _e permeability (magnetic) - kinematic viscosity - electric conductivity - density - M Hartmann number - dimensionless skin-friction - frequency - dimensionless frequency     

Swirling turbulent flow through a curved pipe

An experimental study of swirling turbulent flow through a curved bend and its downstream tangent has been carried out. This study reports on the recovery from swirl and bend curvature and relies on measurements obtained in the downstream tangent and data reported in Part 1 to assess the recovery. Unlike the nonswirling flow case, the present measurements show that the cross-stream secondary flow is dominated by the decay of the solid-body rotation and the total wall shear stress measured at the inner and outer bend (furthest away from the bend center of curvature) is approximately equal. The shear distribution is fairly uniform, even at 1 D downstream of the bend exit. At 49D downstream of the bend exit, the mean axial velocity has recovered to its measured profile at 18D upstream of the bend entrance. Furthermore, the mean tangential velocity is close to zero everywhere and the turbulent shear and normal stresses take another 15D to approximately approach their stationary straight pipe values. Therefore, complete flow recovery from swirl and bend curvature takes a total length of about 85D from the bend entrance. This compares with a recovery length of about 78D for bend curvature alone. The recovery length is substantially shorter than that measured previously in swirling flow through straight pipes and is a consequence of the angular momentum decreasing by approximately 74% across the curved bend. Consequently, the effect of bend curvature is to accelerate swirl decay in a pipe flow.List of symbols C _f total skin friction coefficient, = 2 _w / w ₀ ² - D pipe diameter, = 7.62 cm - De Dean number, = ^1/2 Re = 13,874 - M angular momentum - N _s swirl number, = D/2 W ₀ = 1 - r radial coordinate - R mean bend radius of curvature, = 49.5 cm - Re pipe Reynolds number, = DW ₀ /v= 50,000 - S axial coordinate along the upstream (measured negative) and downstream (measured positive) tangent - U, V, W mean velocities along the radial, tangential and axial directions, respectively - u, v, w mean fluctuating velocities along the radial, tangential and axial directions, respectively - u, v, w root mean square normal stress along the radial, tangential and axial directions, respectively - W ₀ mean bulk velocity, 10 m/s - w total wall friction velocity, = _w / - (w ) _s total wall friction velocity measured as S/D = -18 - turbulent shear stresses - pipe-to-bend radius ratio, = D/2R = 0.077 - axial coordinate measured from bend entrance - fluid kinetic viscosity - fluid density - _w total wall shear stress - azimuthal coordinate measured zero from pipe horizontal diameter near outer bend - angular speed of the rotating section     

Anwendung des Übertragungsverfahrens bei stark abklingenden Lösungsfunktionen

Übersicht Bei stark abklingenden Funktionen wird die Übertragungsmatrix U() aufgespalten in die Anteilc U ₁() e^– und U ₂() e. Der zweite Term spielt am Rand = 0 keinc Rolle. Die unbekannten Anfangswerte sind über die Matrix U ₁(0) an die bekannten gebunden und eindeutig bestimmbar.

---

Summary For strongly decaying solution functions the transfer matrix U() is splitted into the parts U ₁() e^– and U ₂() e. The second term does not influence at the boundary = 0. The unknown initial values are related by the matrix U ₁(0) to the known values and they can be uniquely determined.

     

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     

Hyperbolic phenomena in a strongly degenerate parabolic equation

We consider the equation u _t =((u) (u _x )) _x , where >0 and where is a strictly increasing function with lim _s = <. We solve the associated Cauchy problem for an increasing initial function, and discuss to what extent the solution behaves qualitatively like solutions of the first-order conservation law u _t = ((u)) _x . Equations of this type arise, for example, in the theory of phase transitions where the corresponding free-energy functional has a linear growth rate with respect to the gradient.     

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     

Flow of second-order fluids over an enclosed rotating disc

Summary This paper is devoted to a study of the flow of a second-order fluid (flowing with a small mass rate of symmetrical radial outflow m, taken negative for a net radial inflow) over a finite rotating disc enclosed within a coaxial cylinderical casing. The effects of the second-order terms are observed to depend upon two dimensionless parameters ₁ and ₂. Maximum values ₁ and ₂ of the dimensionless radial distances at which there is no recirculation, for the cases of net radial outflow (m>0) and net radial inflow (m<0) respectively, decrease with an increase in the second-order effects [represented by T(=₁+₂)]. The velocities at ₁ and ₂ as well as at some other fixed radii have been calculated for different T and the associated phenomena of no-recirculation/recirculation discussed. The change in flow phenomena due to a reversal of the direction of net radial flow has also been studied. The moment on the rotating disc increases with T.Nomenclature , , z coordinates in a cylindrical polar system - z ₀ distance between rotor and stator (gap length) - =/z ₀, dimensionless radial distance - =z/z ₀, dimensionless axial distance - _s = _s/z₀, dimensionless disc radius - V =(u, v, w), velocity vector - dimensionless velocity components - uniform angular velocity of the rotor - , p fluid density and pressure - P =p/(₂ z ₀₂ ² , dimensionless pressure - ₁, ₂, ₃ kinematic coefficients of Newtonian viscosity, elastico-viscosity and cross-viscosity respectively - ₁, _{2 2}/z ₀ ² , resp. ₃/z ₀ ² , dimensionless parameters representing the ratio of second-order and inertial effects - m = , mass rate of symmetrical radial outflow - l a number associated with induced circulatory flow - Rm =m/(z ₀₁), Reynolds number of radial outflow - R _l =l/(z ₀₁), Reynolds number of induced circulatory flow - Rz =z ₀ ² /₁, Reynolds number based on the gap - ₁, ₂ maximum radii at which there is no recirculation for the cases Rm>0 and Rm<0 respectively - ₁(T), ₂(T) ₁ and ₂ for different T - U _1(T) ⁽⁺⁾ = dimensionless radial velocity, Rm>0 - V _1(T) ⁽⁺⁾ = , dimensionless transverse velocity, Rm>0 - U _2(T) ^(–) = , dimensionless radial velocity, Rm=–Rn<0, m=–n - V _2(T) ^(–) = , dimensionless transverse velocity, Rm<0 - C _m moment coefficient     

Positively invariant regions for a problem in phase transitions

Positively invariant regions for the system v _t + p(W) _x = V _xx , W _t –V _x = W _xx are constructed where p < 0, w < , w > , p(w) = 0, w , > 0. Such a choice of p is motivated by the Maxwell construction for a van der Waals fluid. The method of an analysis is a modification of earlier ideas of Chueh, Conley, & Smoller [1]. The results given here provide independent L bounds on the solution (w, v).Dedicated to Professor James Serrin on the occasion of his sixtieth birthday     

On natural convection from a vertical plate with a prescribed surface heat flux in porous media

This paper presents a theoretical and numerical investigation of the natural convection boundary-layer along a vertical surface, which is embedded in a porous medium, when the surface heat flux varies as (1 +x ²)), where is a constant andx is the distance along the surface. It is shown that for > -1/2 the solution develops from a similarity solution which is valid for small values ofx to one which is valid for large values ofx. However, when -1/2 no similarity solutions exist for large values ofx and it is found that there are two cases to consider, namely < -1/2 and = -1/2. The wall temperature and the velocity at large distances along the plate are determined for a range of values of .Notation g Gravitational acceleration - k Thermal conductivity of the saturated porous medium - K Permeability of the porous medium - l Typical streamwise length - q _w Uniform heat flux on the wall - Ra Rayleigh number, =gK(q _w /k)l/(v) - T Temperature - Too Temperature far from the plate - u, v Components of seepage velocity in the x and y directions - x, y Cartesian coordinates - Thermal diffusivity of the fluid saturated porous medium - The coefficient of thermal expansion - An undetermined constant - Porosity of the porous medium - Similarity variable, =y(1+x ) ^/3/x ^1/3 - A preassigned constant - Kinematic viscosity - Nondimensional temperature, =(T – T )Ra^1/3 k/qw - Similarity variable, = =y(log_e x)^1/3/x ^2/3 - Similarity variable, =y/x ^2/3 - Stream function