Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A ^* interfacial area of the- interface contained within a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

Stress analysis of the cylinder subjected to arbitrarily distributed axisymmetrical loads

Summary Stress analysis has been carried out for a finite cylinder subjected to arbitrarily distributed axisymmetrical surface loads. Direct stress _x in the axial direction is assumed to be of the form _x = ₀+r ₁ +r ₂ where ₀ to ₂ are functions of x. Using the equations of equilibrium and compatibility the other direct stresses and the shearing stress are expressed by ₁ and ₂. Fundamental equations governing ₁ and ₂ are introduced using the variational principle of complementary energy. From the results of the present analysis it is evident that the boundary conditions can be satisfied completely even for the case where the external forces are specified in complicated form, and that more accurate solutions can easily be obtained by introducing additional terms in _x.

---

Spannungsanalyse für den Zylinder unter axialsymmetrischer Last in beliebiger Verteilung

Übersicht Für einen endlichen Zylinder unter axialsymmetrischer Oberflächenlast in beliebiger Verteilung werden die Spannungen ermittelt. Die Normalspannung in Axialrichtung wird in der Form _x = ₀+r ₁ +r ₂ angesetzt mit ₀, ₁, ₂ als Funktionen von x. Mit Hilfe der Gleichgewichtsund Verträglichkeitsbedingungen werden die anderen Normalspannungen und die Schubspannung durch ₁ und ₂ ausgedrückt. Über das Variationsprinzip für die Komplementärenergie werden die grundlegenden Gleichungen für ₁ und ₂ eingeführt. Die Ergebnisse zeigen, daß die Randbedingungen selbst für komplizierte Belastungsarten vollständig erfüllbar sind und mit zusätzlichen Termen in _x mühelos noch genauere Lösungen bestimmt werden können.

     

A rotating field mhd hydrostatic thrust bearing

It is found that the load capacity of a magnetohydrodynamic thrust bearing with a rotating disk can be increased by rotating the axial magnetic field at a suitable speed in a direction opposite to that of the disk rotation. This method of improving the bearing performance is considered to be efficient if the Hartmann number is not too large. Thus for a given load, the size and weight of the magnet to be used in a thrust bearing with rotating field can be reduced considerably.Nomenclature a radius of plenum recess - b outside disk radius - B ₀ magnetic induction of applied axial magnetic field - hE ₀ ^1/2/a ^1/2, nondimensionalized electric field - E ₀ radial electric field at r=a - E _r radial electric field - h half of lubricant film thickness - M (B ₀ ² h ²/)^1/2, Hartmann number - P pressure - P _e pressure at r=b - P ₀ pressure at r=a - Q volume flow rate of lubricant - Q ₀ volume flow rate of a nonrotating bearing in the absence of applied magnetic field - r radial coordinate - u, v fluid velocity components in radial and circumferential directions, respectively - W load capacity of bearing - W ₀ load capacity of a nonrotating bearing in the absence of a magnetic field having a flow rate which the same bearing would have at Hartmann number M - z axial coordinate - azimuthal coordinate - coefficient of viscosity of lubricant - _e magnetic permeability - fluid density - electrical conductivity - angular velocity of rotating disk - _C critical disk velocity at which W=0 - _M angular velocity of axial magnetic field - optimum angular velocity of magnetic field On leave of absence from Department of Aero-Space Engineering, University of Notre Dame, Notre Dame (Ind.), U.S.A.     

Zur Variationsformulierung nichtlinearer Randwertprobleme

Übersicht Es werden verschiedene Bedingungen aufgestellt, die es erlauben, die durch die beiden (Systeme von) nichtlinearen DifferentialgleichungenA (u, ) = q, B (u, ) = und Randbedingungen zusammen mit den nichtlinearen algebraischen Relationenq = C(u, ), = D(u, ) beschriebene Aufgabe durch äquivalente Variationsprobleme zu ersetzen. Dabei zeigt sich ein enger Zusammenhang mit den in der Festkörpermechanik wohlbekannten Prinzipien der virtuellen Verschiebungen und der virtuellen Kräfte. Die auf systematischem Weg konstruierten Variationsfunktionale enthalten viele in der Physik bekannte Funktionale als Sonderfälle, insbesondere jene, die in der Elastomechanik nach Green, Castigliano, Hellinger, Reißner, Hu und Washizu benannt werden.

---

Summary In this paper there are established various conditions which allow a variational formulation of the problem described by the two (systems of) nonlinear differential equationsA(u, ) = q, B(u, ) = and boundary conditions together with the nonlinear algebraic relationsq = C(u, ), = D(u, ). Besides a close relationship is revealed to the principles of virtual displacements and virtual forces which are wellknown in solid mechanics. The systematically constructed variational functional contain many functionals in physics as special cases, mainly those of Green, Castigliano, Hellinger, Reißner, Hu and Washizu in elastomechanics.

     

Higher-order correction of effective permeability of heterogeneous isotropic formations of lognormal conductivity distribution

Steady flow of an incompressible fluid takes place in a porous formation of spatially variable hydraulic conductivityK. The latter is regarded as a lognormal stationary random space function and Y=ln(K/K _G ), whereK _G is the geometric mean ofK, is characterized by its variance ² and correlation scale I. Exact results are known for the effective conductivityK _eff in one- and two-dimensional flows. In contrast, only a first-order term in a perturbation expansion in ² has been derived exactly for the three-dimensional flow. A conjecture has been made in the past onK _eff for any ², but it was not yet proved exactly. This study derived the exact nonlinear correction, i.e. the termO(⁴) ofK _eff, which is found to be the one resulting from the conjecture, strengthening the confidence in it. It is also shown that the self-consistent approximation leads to the exact results for one-dimensional and two-dimensional flows, but underestimates the nonlinear correction ofK _eff for in the three-dimensional case.     

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.     

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

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

---

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

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

---

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

The 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     

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     

Axial dispersion of solids in rotary solid flow systems

Summary A simple unidirectional diffusion model is employed to analyze the axial dispersion of solid particles flowing through a rotary solid flow system, namely a rotary dryer. It is shown that the reciprocal of the Peclet number D/uL is uniquely correlated as a function of the dimensionless number F/dSN which characterizes the operating conditions of the rotary dryer.List of Symbols C concentration of tracer, mass/(length)³ - d diameter of rotary dryer, length - d _p diameter of solid particles, length - D longitudinal dispersion coefficient or axial mixing coefficient, (length)²/time - F volumetric flow rate of solid, (length)³/(length)² time - L length of rotary dryer, length - N rate of rotations of dryer, time^–1 - Q volume of tracer injected, based on bulk density of particles, (length)³ - S slope of the rotary dryer - u average flow velocity, based on effective flow volume of dryer, length/time - v volumetric flow rate, based on bulk density of particles, (length)³/time - V effective volume of rotary dryer, (length)³ - x distance from entrance of experimental section of dryer, length Greek letters time, measured from instant of introducing tracer into flowing material - 1– volumetric solid hold-up fraction - standard deviation - ² variance - _r relative standard deviation     

On the effects of a periodic, spanwise variation of suction upon asymptotic suction flow

Summary The effects of superposing streamwise vorticity, periodic in the lateral direction, upon two-dimensional asymptotic suction flow are analyzed. Such vorticity, generated by prescribing a spanwise variation in the suction velocity, is known to play an important role in unstable and turbulent boundary layers. The flow induced by the variation has been obtained for a freestream velocity which (i) is steady, (ii) oscillates periodically in time, (iii) changes impulsively from rest. For the oscillatory case it is shown that a frequency can exist which maximizes the induced, unsteady wall shear stress for a given spanwise period. For steady flow the heat transfer to, or from a wall at constant temperature has also been computed.Nomenclature (x, y, z) spatial coordinates - (u, v, w) corresponding components of velocity - (, , ) corresponding components of vorticity - t time - stream function for v and w - v _w mean wall suction velocity - nondimensional amplitude of variation in wall suction velocity - characteristic wavenumber for variation in direction of z - T temperature - P pressure - density - coefficient of kinematic viscosity - coefficient of thermal diffusivity - (/v _w)² - frequency of oscillation of freestream velocity - nondimensional amplitude of freestream oscillation - /v _w ² - z z - y –v _w y/ - v _w ² t/4 - /v _w - U ₀ characteristic freestream velocity - u/U ₀ - coefficient of viscosity - _w wall shear stress - Prandtl number (/) - q heat transfer to wall - T _w wall temperature - T (T _w–T)/(T _w–)     

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

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

Unterkühltes Sieden strömender Flüssigkeitsgemische

Zusammenfassung Es werden Meßergebnisse zum Wärmeübergang beim unterkühlten Sieden von Isopropanol/Wasser-Gemischen in einem senkrechten Ringspalt bei Aufwärtsströmung vorgestellt. Der Einfluß der Versuchsparameter Wärmestromdichte, Flüssigkeitsgeschwindigkeit, Flüssigkeitsunterkühlung und Flüssigkeitszusammensetzung auf den Wärmeübergang im Blasensiedebereich wird dargestellt. Die gemessenen Wärmeübergangskoeffizienten werden mit den Vorhersagen zweier Korrelationen für gesättigtes Sieden aus der Literatur verglichen. Die Übereinstimmung zwischen Meßwerten und Korrelationen ist zufriedenstellend.

---

Subcooled flow-boiling of mixtures

Experimental results for subcooled flow-boiling heat transfer to isopropanol/water-mixtures in a vertical annulus during upward flow are reported. The influence of heat flux, flow velocity subcooling and mixture composition on heat transfer in the nucleate boiling region is discussed. The measured subcooled boiling heat transfer coefficients are compared with the predictions of two correlations for saturated boiling suggested in the literature. It was found that these correlations agree well with the present measurements.

---

Formelzeichen A ₀ empirische Stoffgröße - D relative Abweichung - D _b Blasendurchmesser, m - g Erdbeschleunigung, m/s² - h _v Verdampfungsenthalpie, J/kg - p Druck, Pa - Wärmestromdichte, W/m² - s Abstand des Thermoelementes von der Heizstaboberfläche, m - T Temperatur, K - T _sub Unterkühlung, K - x Flüssigkeitskonzentration in Molprozent, mol% - y Gaskonzentration in Molprozent, mol% - Wärmeübergangskoeffizient nach Gleichung (4), W/m² K - ₀ Wärmeübergangskoeffizient nach Gleichung (2), W/m² K - Stoffübergangskoeffizient, m/s - Wärmeleitfähigkeit, W/m K - Dichte, kg/m³ - Oberflächenspannung, N/m - Randwinkel, Grad Indizes b Kern (bulk) - g Gas - id ideal - krit kritisch - l Flüssigkeit (liquid) - s Sättigung - th Thermoelement - w Wand - 1 leichter siedende Komponente (Isopropanol) - 2 schwerer siedende Komponente (Wasser) Herrn Prof. Dr. rer. nat. K. Bier zum 65. Geburtstag gewidmet     

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

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

Rheological properties of polyethylene oxide solutions

Summary Measurements were made on solutions of Polyethylene oxide (WSR-301) varying in concentration from .0511 g/dl to 4.014 g/dl, prepared from two samples of dry material of different ages (I, II), using aWeissenberg Rheogoniometer with cone-and-plate and parallel-plate geometries, and also using capillary viscometers. Steady shear data were obtained for eight decades of strain-rates (10^–3 <k < 10⁵ sec^–1), and oscillatory data for over four decades of frequency (10^–3 <f < 10¹ Hz). Results are presented for the shear-dependent viscosity,(k), normal stress differences, ₁(k), ₂(k), and the complex viscosity, ^*(f). It was found that characteristic fluid times obtained from continuum arguments correlated the experimental(k), (f) andG(f) data.Using the ^*(f) data, the stress relaxation function,(), was calculated, from which the second-order fluid coefficients and ₀ were obtained, and compared to the directly measured values.Evidence is given to show that the sign of ₁(k) varies both with concentration and strainrate.Using solutions prepared from sample II, correlations with the material properties of solutions of sample I were found which indicated the effect of aging on stored dry samples.

---

Zusammenfassung Es wurden Viskositätsmessungen an Polyäthylenoxid-(WSR-301)-Lösungen mit Konzentrationen zwischen 0,0511 g/dl und 4,014 g/dl ausgeführt, die aus zwei Proben von Trockenmaterial verschiedenen Alters (I, II) genommen waren. Verwendet wurde einWeissenberg-Rheogoniometer mit Kegel-Platte- und Parallel-Platten-Geometrie sowie verschiedene Kapillarviskosimeter. Werte für die stationäre Scherung wurden über acht Dekaden der Schergeschwindigkeit (10^–3 <k < 10⁵ s^–1) erhalten, solche für periodische Beanspruchung über mehr als vier Dekaden der Frequenz (10^–3 <f < 10 Hz). Es werden die Werte der Scherviskosität(k), der Normalspannungsdifferenzen ₁(k) und ₂(k) sowie der komplexen Viskosität ^*(f) mitgeteilt. Man findet, daß die experimentell ermittelten Werte von(k), (f) undG(f) mit Hilfe charakteristischer Zeitkonstanten, die man aus kontinuumsmechanischen Überlegungen gewinnt, korreliert werden können.Aus dem Verlauf von ^*(f) wurde die Spannungsrelaxationsfunktion() berechnet, woraus sich die Koeffizienten zweiter Ordnung und ₀ bestimmen lassen. Diese wurden mit den auf direkte Weise gewonnenen Werten verglichen. Es wird nachgewiesen, daß das Vorzeichen von ₁(k) sowohl bei der Veränderung der Konzentration als auch der Deformationsgeschwindigkeit wechselt.Durch Vergleich der an den Proben I und II erhaltenen Ergebnisse wird auf Alterungserscheinungen bei der trocken gelagerten Probe geschlossen.

---

---

With 19 figures and 1 table     

Local energy decay for the wave equation on the exterior of two balls or convex bodies

An algebraic rate of decay of local energy, nonuniform with respect to the initial data, is established for solutions of the Dirichlet and Neumann problems for the scalar wave equation defined on the exterior V³ of two balls or of two convex bodies. That is, for given initial data f(x)=u(x), 0 and g(x)= u _t (x, 0), if u solves u _tt in V with either u(x, t)=0 or u _n (x,t)+(x) u(x,t,)-0 ((x)0) on V, then there exists a constant T ₀, depending upon (f, g), such that the local energy (the energy in any compact set) of u at t=T is bounded from above by QE(0)T ^–1 for TT ₀, where E(0) is the total initial energy of u and Q is a positive constant, independent of u, that depends upon V.     

Zur Kapillarviskosimetrie verdünnter Lösungen

Zusammenfassung Es werden allgemeine Dimensionierungsregeln für Kapillarviskosimeter angegeben. Der Einfluß der Oberflächenspannung () der zu untersuchenden Flüssigkeit auf die mittlere hydrostatische Druckhöhe ist für den Fall berechnet, daß das Gefäß zwischen den beiden Meßmarken eine Kugel mit dem Volumen von 50 cm³ ist. Die berechnete mittlere Druckhöhe ist eine lineare Funktionvon /. Daher ist es im Fall desUbbelohde-Viskosimeters möglich, den Kapillarzug durch eine entsprechende Krümmung des hängenden Niveaus im Ausfluß der Kapillare zu kompensieren.

---

Summary General rules to dimension capillary viscosimeters are given. The effect of the surface tension () of the fluid under test on the mean hydrostatic head is calculated for a sphere of 50 cm³ between the measuring marks.The calculated mean hydrostatic head is a linear function of/. Therefore it is possible in the case of theUbbelohde-viscosimeter to compensate the capillary traction by appropriate curvature of the suspended level in the outflow of the capillary.

     

A generalized Norton-Hoff model and the Prandtl-Reuss law of plasticity

The aim of this article is to study the quasistatic evolution of a three-dimensional elastic-perfectly plastic solid which satisfies the Prandtl-Reuss law. The evolution of the field of stresses -which solves a time dependent variational inequality — and that of the field of displacements u, have been described in previous works [15], [26], [35], [36], [37] but it was not shown there that and u satisfy indeed the Prandtl-Reuss constitutive law. In this article we find and u in a class of functions which are sufficiently regular for the Prandtl-Reuss law to make sense and we prove that and u satisfy the constitutive law. This result is attained by considering the elastic-perfectly plastic model as the limit of a family of elastic-visco-plastic models like those of Norton and Hoff. The Norton-Hoff type models which we introduce depend on a viscosity parameter > 0; we study the perturbed models (i.e. > 0 fixed) and then we pass to the limit 0.Dedicated to James Serrin on the occasion of his 60^th Birthday     

Transport in ordered and disordered porous media V: Geometrical results for two-dimensional systems

In this paper we continue the geometrical studies of computer generated two-phase systems that were presented in Part IV. In order to reduce the computational time associated with the previous three-dimensional studies, the calculations presented in this work are restricted to two dimensions. This allows us to explore more thoroughly the influence of the size of the averaging volume and to learn something about the use of anon-representative region in the determination of averaged quantities.

Nomenclature

Roman Letters A interfacial area of the– interface associated with the local closure problem, m² - a _i i=1, 2, gaussian probability distribution used to locate the position of particles - l unit tensor - characteristic length for the-phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1,2,3 lattice vectors, m - m convolution product weighting function - m _V special convolution product weighting function associated with a unit cell - n _i i=1, 2 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume,V, m³ - x position of the centroid of an averaging area, m - x ₀ reference position of the centroid of an averaging area, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters V /V, volume average porosity - _{a i} standard deviation ofa _i - _r standard deviation ofr - intrinsic phase average of      

Diffusion of moisture in deformable porous media. Anisotropic fibrous materials

This paper presents a study on the deformation of anisotropic fibrous porous media subjected to moistening by water in the liquid phase. The deformation of the medium is studied by applying the concept of effective stress. Given the structure of the medium, the displacement of the solid matrix is not taken into account with respect to the displacement of the liquid phase. The transport equations are derived from the model proposed by Narasimhan. The transport coefficients and the relation between the variation in apparent density and effective stress are obtained by test measurements. A numerical model has been established and applied for studying drip moistening of mineral wool samples capable or incapable of deformation.Nomenclature D mass diffusion coefficient [L²t^–1] - e void fraction - g gravity acceleration [Lt^–2] - J mass transfer density [ML^–2t^–1] - K hydraulic conductivity [Lt^–1] - K _s hydraulic conductivity of the solid phase [Lt^–1] - K ^* hydraulic conductivity of the deformable porous medium [Lt^–1] - P pressure of moistening liquid [ML^–1 t^–2] - S degree of saturation - t time [t] - V speed [Lt^–1] - X horizontal coordinate [L] - Z vertical coordinate measured from the bottom of porous medium [L] - z z-coordinate [L] Greek Letters porosity - ₁ total hydric potential [L] - _g gas density [ML^–3] - ₁ liquid density [ML^–3] - ₀ apparent density [ML^–3] - _s density of the solid phase [ML^–3] - density of the moist porous medium [ML^–3] - external load [ML^–1t^–2] - effective stress [ML^–1t^–2] - bishop's parameter - matrix potential or capillary suction [L] Indices g gas - 1 moistening liquid - p direction perpendicular to fiber planes - s solid matrix - t direction parallel to fiber planes - v pore Exponent * movement of solid particles taken into account