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.     

Control of low-speed turbulent separated flow using jet vortex generators

A parametric study has been performed with jet vortex generators to determine their effectiveness in controlling flow separation associated with low-speed turbulent flow over a two-dimensional rearward-facing ramp. Results indicate that flow-separation control can be accomplished, with the level of control achieved being a function of jet speed, jet orientation (with respect to the free-stream direction), and jet location (distance from the separation region in the free-stream direction). Compared to slot blowing, jet vortex generators can provide an equivalent level of flow control over a larger spanwise region (for constant jet flow area and speed).Nomenclature C _p pressure coefficient, 2(P-P_)/V ² - C _Q total flow coefficient, Q/ v - D ₀ jet orifice diameter - Q total volumetric flow rate - R Reynolds number based on momentum thickness - u fluctuating velocity component in the free-stream (x) direction - V free-stream flow speed - VR ratio of jet speed to free-stream flow speed - x coordinate along the wall in the free-stream direction - jet inclination angle (angle between the jet axis and the wall) - jet azimuthal angle (angle between the jet axis and the free-stream direction in a horizontal plane) - boundary-layer thickness - momentum thickness - lateral distance between jet orifices A version of this paper was presented at the 12th Symposium on Turbulence, University of Missouri-Rolla, 24–26 Sept. 1990     

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.     

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.     

A note on the unsteady motion of a viscous conducting liquid between two porous concentric circular cylinders acted on by a radial magnetic field

In this note the author has investigated some problems of a flow of conducting liquid through two porous non-conducting infinite circular cylinders rotating with various angular velocities for some time in the presence of a radial magnetic field.It is assumed that the rate of suction at the inner cylinder is equal to the rate of injection at the outer cylinder. Furthermore the induced electric and magnetic fields are neglected.Nomenclature a, b radii of coaxial cylinders - h ₁, h ₂ constants - m, n constants - P hydrostatic pressure - t time - A, B functions of r and b - B ₀ magnetic induction vector B ₀= _e H - H magnetic field vector - H intensity of magnetic field - L, M functions of a and b - S Suction parameter - Coefficient of viscosity - _e magnetic permeability - Kinematic coefficient of viscosity - density of liquid - conductivity of liquid - Small time - Constant - ₁, ₂ angular vrlocities - , ₁, ₂ angular velocities.     

Thin-Walled Beams: the Case of the Rectangular Cross-Section

In this paper we present an asymptotic analysis of the three-dimensional problem for a thin linearly elastic cantilever =×(0,l) with rectangular cross-section of sides and ², as goes to zero. Under suitable assumptions on the given loads, we show that the three-dimensional problem converges in a variational sense to the classical one-dimensional model for extension, flexure and torsion of thin-walled beams. Mathematics Subject Classifications (2000) 474K20, 74B10, 49J45.     

Investigations of mean and turbulent flow characteristics of a two dimensional plane diffuser

This paper reports the investigation of mean and turbulent flow characteristics of a two-dimensional plane diffuser. Both experimental and theoretical details are considered. The experimental investigation consists of the measurement of mean velocity profiles, wall static pressure and turbulence stresses. Theoretical study involves the prediction of downstream velocity profiles and the distribution of turbulence kinetic energy using a well tested finite difference procedure. Two models, viz., Prandtl's mixing length hypothesis and k- model of turbulence, have been used and compared. The nondimensional static pressure distribution, the longitudinal pressure gradient, the pressure recovery coefficient, percentage recovery of static pressure, the variation of U _max/U _bar along the length of the diffuser and the blockage factor have been valuated from the predicted results and compared with the experimental data. Further, the predicted and the measured value of kinetic energy of turbulence have also been compared. It is seen that for the prediction of mean flow characteristics and to evaluate the performance of the diffuser, a simple turbulence model like Prandtl's mixing length hypothesis is quite adequate.List of symbols C ₁ , C ₂ ,C turbulence model constants - F _x body force - k kinetic energy of turbulence - l _m mixing length - L length of the diffuser - u, v, w rms value of the fluctuating velocity - u, v, w turbulent component of the velocity - mean velocity in the x direction - _A average velocity at inlet - U _bar average velocity in any cross section - U _max maximum velocity in any cross section - V mean velocity in the y direction - W local width of the diffuser at any cross section - x, y coordinates - dissipation rate of turbulence - _m eddy diffusivity - Von Karman constant - mixing length constant - _l laminar viscosity - _eff effective viscosity - v kinematic viscosity - density - _k effective Schmidt number for k - effective Schmidt number for - stream function - non dimensional stream function     

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

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.     

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

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

---

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

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

     

On nonlinear problems of mixed type: A qualitative theory using infinite-dimensional center manifolds

We consider the equation a(y)u_xx+div_y(b(y)_yu)+c(y)u=g(y, u) in the cylinder (–l,l)×, being elliptic where b(y)>0 and hyperbolic where b(y)<0. We construct self-adjoint realizations in L₂() of the operatorAu= (1/a) div_y(b_yu)+(c/a) in the case ofb changing sign. This leads to the abstract problem u_xx+Au=g(u), whereA has a spectrum extending to + as well as to –. For l= it is shown that all sufficiently small solutions lie on an infinite-dimensional center manifold and behave like those of a hyperbolic problem. Anx-independent cross-sectional integral E=E(u, u_x) is derived showing that all solutions on the center manifold remain bounded forx ±. For finitel, all small solutionsu are close to a solution on the center manifold such that u(x)-(x) Ce ^-(1-|x|) for allx, whereC and are independent ofu. Hence, the solutions are dominated by hyperbolic properties, except close to the terminal ends {±1}×, where boundary layers of elliptic type appear.     

Evaluation of the dynamic-mechanical properties of solids from a cooperative model for stress relaxation

For many solid materials the stress relaxation process obeys the universal relationF = – (d/d lnt)_max = (0.1 ± 0.01) ( ₀ — _i ), regardless of the structure of the material. Here denotes the stress,t the time, ₀ the initial stress of the experiment and _i the internal stress. A cooperative model accounting for the similarity in relaxation behaviour between different materials was developed earlier. Since this model has a spectral character, the concepts of linear viscoelasticity are used here to evaluate the corresponding prediction of the dynamic mechanical properties, i.e. the frequency dependence of the storageE () and lossE () moduli. Useful numerical approximations ofE () andE () are also evaluated. It is noted that the universal relation in stress relaxation had a counterpart in the frequency dependence ofE (). The theoretical prediction of the loss factor for high-density polyethylene is compared with experimental results. The agreement is good.     

Qualitative behavior of dissipative wave equations on bounded domains

The qualitative behavior of solutions of the mixed problem u_tt = u-a(x)u_t in IR x , u=0 on IR x , is studied in the case when a>0 and IRⁿ is bounded. Roughly speaking, if aa_min>0, then solutions decay at least as fast as exp t( –1/2a_min), with the possible exception of a finite dimensional set of smooth solutions whose existence is associated with a phenomenon of overdamping. If a_max is sufficiently small, depending on , then no overdamping occurs.Partially supported by NSF grant NSF GP 34260.This work was partially supported by the National Science Foundation under Grant No. GP 34260     

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.

     

About the Regularized Navier–Stokes Equations

The first goal of this paper is to study the large time behavior of solutions to the Cauchy problem for the 3-dimensional incompressible Navier–Stokes system. The Marcinkiewicz space L^3, is used to prove some asymptotic stability results for solutions with infinite energy. Next, this approach is applied to the analysis of two classical regularized Navier–Stokes systems. The first one was introduced by J. Leray and consists in mollifying the nonlinearity. The second one was proposed by J.-L. Lions, who added the artificial hyper-viscosity (–)^{/ 2}, > 2 to the model. It is shown in the present paper that, in the whole space, solutions to those modified models converge as t toward solutions of the original Navier–Stokes system.     

Unsteady convective diffusion in a rotating parallel plate channel

The paper presents an exact analysis of the dispersion of a passive contaminant in a viscous fluid flowing in a parallel plate channel driven by a uniform pressure gradient. The channel rotates about an axis perpendicular to its walls with a uniform angular velocity resulting in a secondary flow. Using a generalized dispersion model which is valid for all time, we evaluate the longitudinal dispersion coefficientsK _i (i=1, 2, ...) as functions of time. It is shown thatK ₁=0 andK ₃,K ₄, ... decay rapidly in comparison withK ₂. ButK ₂ decreases with increasing (the dimensionless rotation parameter) for values of upto approximately =2.2. ThereafterK ₂ increases with further increase in and its value gets saturated for large values of (say, 500) and does not change any further with increase in . A physical explanation of this anomalous behaviour ofK ₂ is given.

---

Instationäre konvektive Diffusion in einem rotierenden Parallelplattenkanal

Zusammenfassung In dieser Untersuchung wird eine exakte Analyse der Ausbreitung eines passiven Kontaminierungsstoffes in einer zähen Flüssigkeit gegeben, die, befördert durch einen gleichförmigen Druckgradienten, in einem Parallelplattenkanal strömt. Der Kanal rotiert mit gleichförmiger Winkelgeschwindigkeit um eine zu seinen Wänden senkrechte Achse, wodurch sich eine Sekundärströmung ausbildet. Unter Verwendung eines generalisierten, für alle Zeiten gültigen Dispersionsmodells werden die longitudinalen DispersionskoeffizientenK _i (i=1, 2, ...) als Funktionen der Zeit ermittelt. Es wird gezeigt, daßK ₁=0 gilt und dieK ₃,K ₄, ... gegenüberK ₂ schnell abnehmen.K ₂ nimmt ab, wenn , der dimensionslose Rotationsparameter, bis etwa zum Wert 2,2 ansteigt. Danach wächstK ₂ mit bis auf einem Endwert an, der etwa ab =500 erreicht wird. Dieses anomale Verhalten vonK ₂ findet eine physikalische Erklärung.

---

List of symbols C solute concentration - D molecular diffusivity - K _i longitudinal dispersion coefficients - 2L depth of the channel - P ₀ dimensionless pressure gradient along main flow - Pe Péclet number - q velocity vector - Q _x,Q _y mass flux along the main flow and the secondary flow directions - dimensionless average velocity along the main flow direction - (x, y, z) Cartesian co-ordinates Greek symbols dimensionless rotation parameter - the inclination of side walls withx-axis - kinematic viscosity - fluid density - dimensionless time - angular velocity of the channel - dimensionless distance along the main flow direction - dimensionless distance along the vertical direction - dimensionless solute concentration - integral of the dispersion coefficientK ₂() over a time interval     

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