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     

The wedge subjected to tractions: a paradox resolved

The classical two-dimensional solution provided by Lévy for the stress distribution in an elastic wedge, loaded by a uniform pressure on one face, becomes infinite when the opening angle 2 of the wedge satisfies the equation tan 2 = 2. Such pathological behavior prompted the investigation in this paper of the stresses and displacements that are induced by tractions of O(r ^–) as r0. The key point is to choose an Airy stress function which generates stresses capable of accommodating unrestricted loading. Fortunately conditions can be derived which pre-determine the form of the necessary Airy stress function. The results show that inhomogeneous boundary conditions can induce stresses of O(r ^–), O(r ^– ln r), or O(r ^– ln² r) as r0, depending on which conditions are satisfied. The stress function used by Williams is sufficient only if the induced stress and displacement behavior is of the power type. The wedge loaded by uniform antisymmetric shear tractions is shown in this paper to exhibit stresses of O(ln r) as r0 for the half-plane or crack geometry. At the critical opening angle 2, uniform antisymmetric normal and symmetric shear tractions also induce the above type of stress singularity. No anticipating such stresses, Lévy used an insufficiently general Airy stress function that led to the observed pathological behavior at 2.     

Laminar flow in the entrance region of a porous-wall channel

Summary The effect of fluid injection at the walls of a two-dimensional channel on the development of flow in the entrance region of the channel has been investigated. The integral forms of the boundary layer equations for flow in the channel were set up for an injection velocity uniformly distributed along the channel walls.With an assumed polynomial of the n-th degree for the one-parameter velocity profile a solution of the above boundary layer equations was obtained by an iteration method. A closed form solution was also obtained for the case when a similar velocity profile was assumed. The agreement between the entrance region velocity profiles of the present analysis for an impermeable-walled channel and of Schlichting¹) and Bodoia and Osterle²) is found to be very good.The results of the analysis show that fluid injection at the channel walls increases the rate of the growth of the boundary layer thickness, and hence reduces considerably the entrance length required for a fully developed flow.Nomenclature h half channel thickness - L entrance length with wall-injection - L ₀ entrance length without wall-injection - p static pressure - p=p/U ₀ ² dimensionless pressure - Re=U ₀ h/ Reynolds number at inlet cross-section - u velocity in the x direction at any point in the channel - =u/U ₀ dimensionless velocity in the x direction at any point in the channel - U _av average velocity at a channel cross-section - U _c center line velocity - U ₀ inlet cross-section velocity - _c =U _c /U ₀ dimensionless center line velocity - v velocity in the y direction at any point in the channel - v ₀ constant injection velocity of fluid at the wall - v=v/v ₀ dimensionless velocity in the y direction at any point in the channel - x distance along the channel wall measured from the inlet cross-section - x=x/hRe dimensionless distance in the x direction - y distance perpendicular to the channel wall - y=y/h dimensionless distance in the y direction - thickness of the boundary layer - =/h dimensionless boundary layer thickness - =/ dimensionless distance within the boundary layer region - =v ₀ h/ injection parameter or injection Reynolds number - kinematic viscosity - 1+ie - mass density of the fluid - parameter defined in (14)     

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

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.     

The moving, plane heat source in an elastic medium

Summary This note presents an exact solution for the stress and displacement field in an unbounded and transversely constrained elastic medium resulting from the motion of a plane heat source travelling through the medium at constant speed in the direction normal to the source plane.Nomenclature mass density - diffusivity - thermal conductivity - Q heat emitted by plane heat source per unit time per unit area - speed of propagation of plane heat source - shear modulus - Poisson's ratio - T temperature - _x, _y, _z normal stress components - u _x, u_y, u_z displacement components - c speed of irrotational waves - t time - x, y, z Cartesian coordinates - =x–vt moving coordinate     

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 = ₂     

Optical co-ordinates in general relativity 2 — the problem of distance

Summary Let denote the congruence of null geodesics associated with a given optical observer inV ₄. We prove that determines a unique collection of vector fieldsM₍₎ ( =1, 2, 3) and ₍₀₎ overV ₄, satisfying a weak version of Killing's conditions.This allows a natural interpretation of these fields as the infinitesimal generators of spatial rotations and temporal translation relative to the given observer. We prove also that the definition of the fieldsM₍₎ and ₍₀₎ is mathematically equivalent to the choice of a distinguished affine parameter f along the curves of, playing the role of a retarded distance from the observer.The relation between f and other possible definitions of distance is discussed.

---

Sommario Sia la congruenza di geodetiche nulle associata ad un osservatore ottico assegnato nello spazio-tempoV ₄. Dimostriamo che determina un'unica collezione di campi vettorialiM₍₎ ( =1, 2, 3) e ₍₀₎ inV ₄ che soddisfano una versione in forma debole delle equazioni di Killing. Ciò suggerisce una naturale interpretazione di questi campi come generatori infinitesimi di rotazioni spaziali e traslazioni temporali relative all'osservatore assegnato. Dimostriamo anche che la definizione dei campiM₍₎, ₍₀₎ è matematicamente equivalente alla scelta di un parametro affine privilegiato f lungo le curve di, che gioca il ruolo di distanza ritardata dall'osservatore. Successivamente si esaminano i legami tra f ed altre possibili definizioni di distanza in grande.

---

---

Work performed in the sphere of activity of: Gruppo Nazionale per la Fisica Matematica del CNR.     

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

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.

     

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

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

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

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

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     

Diffusion in anisotropic porous media

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

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.     

Supersonic nonequilibrium flow past segmental bodies of a mixture simulating the atmosphere of venus

Calculations of the flow of the mixture 0.94 CO₂+0.05 N₂+0.01 Ar past the forward portion of segmentai bodies are presented. The temperature, pressure, and concentration distributions are given as a function of the pressure ahead of the shock wave and the body velocity. Analysis of the concentration distribution makes it possible to formulate a simplified model for the chemical reaction kinetics in the shock layer that reflects the primary flow characteristics. The density distributions are used to verify the validity of the binary similarity law throughout the shock layer region calculated.The flow of a CO₂+N₂+Ar gas mixture of varying composition past a spherical nose was examined in [1]. The basic flow properties in the shock layer were studied, particularly flow dependence on the free-stream CO₂ and N₂ concentration.New revised data on the properties of the Venusian atmosphere have appeared in the literature [2, 3] One is the dominant CO₂ concentration. This finding permits more rigorous formulation of the problem of blunt body motion in the Venus atmosphere, and attention can be concentrated on revising the CO₂ thermodynamic and kinetic properties that must be used in the calculation.The problem of supersonic nonequilibrium flow past a blunt body is solved within the framework of the problem formulation of [4].Notation V body velocity - shock wave standoff - universal gas constant - ratio of frozen specific heats - hRt/m enthalpy per unit mass undisturbed stream P pressure - density - T temperature - m molecular weight - c_p specific heat at constant pressure - (X) concentration of component X (number of particles in unit mass) - R body radius of curvature at the stagnation point - _j rate of j-th chemical reaction shock layer P V ² pressure - density - TT temperature - mm molecular weight Translated from Izv. AN SSSR. Mekhanika Zhidkosti i Gaza, Vol. 5, No. 2, pp. 67–72, March–April, 1970.The author thanks V. P. Stulov for guidance in this study.     

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.     

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     

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