On the continuum modelling of porous media containing fluid: a molecular viewpoint with particular attention to scale

Mass conservation and linear momentum balance relations for a porous body and any fluid therein, valid at any given length scale in excess of nearest-neighbour molecular separations, are established in terms of local weighted averages of molecular quantities. The mass density field for the porous body at a given scale is used to identify its boundary at this scale, and a porosity field is defined for any pair of distinct length scales. Specific care is paid to the interpretation of the stress tensor associated with each of the body and fluid at macroscopic scales, and of the force per unit volume each exerts on the other. Consequences for the usual microscopic and macroscopic viewpoints are explored.Nomenclature material system; Section 2.1. - porous body (example of a material system); Sections 2.1, 3.1, 4.1 - fluid body (example of a material system); Sections 2.1, 3.1, 4.1 - weighting function; Sections 2.1, 2.3 - ,h weighting function corresponding to spherical averaging regions of radius and boundary mollifying layer of thicknessh; Section 3.2 - Euclidean space; Section 2.1 - V space of all displacements between pairs of points in; Section 2.1 - mass density field corresponding to; (2.3)₁ - ^P , ^f mass density fields for , ; (4.1) - P momentum density field corresponding to; (2.3)₂ - v velocity field corresponding to; (2.4) - S _r (X) interior of sphere of radiusr with centre at pointx; (3.3) - boundary ofany region - region in which ^p > 0 with = ,h; (3.1) - subset of whose points lie at least+h from boundary of ; (3.4) - abbreviated versions of ; Section 3.2, Remark 4 - strict interior of ; (3.7) - analogues of for fluid system ; Section 3.2 - general version of corresponding to any choice of weighting function; (4.6) - interfacial region at scale; (3.8) - ₀ scale of nearest-neighbour separations in ; Section 3.2. Remark 1 - porosity field at scales ( ₁; ₂); (3.9) - pore space at scales ( ₁; ₂); (3.12)     

A technique for evaluation of three-dimensional behavior in turbulent boundary layers using computer augmented hydrogen bubble-wire flow visualization

A system is described which allows the recreation of the three-dimensional motion and deformation of a single hydrogen bubble time-line in time and space. By digitally interfacing dualview video sequences of a bubble time-line with a computer-aided display system, the Lagrangian motion of the bubble-line can be displayed in any viewing perspective desired. The u and v velocity history of the bubble-line can be rapidly established and displayed for any spanwise location on the recreated pattern. The application of the system to the study of turbulent boundary layer structure in the near-wall region is demonstrated.List of Symbols Reynolds number based on momentum thickness u /v - t⁺ nondimensional time - u shear velocity - u local streamwise velocity, x-direction - u ⁺ nondimensional streamwise velocity - v local normal velocity, -direction - x ⁺ nondimensional coordinate in streamwise direction - ⁺ nondimensional coordinate normal to wall - ₊ ^wire wire nondimensional location of hydrogen bubble-wire normal to wall - z ⁺ nondimensional spanwise coordinate - momentum thickness - v kinematic viscosity - _W wall shear stress     

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

Effects of couple-stresses on the dispersion of a soluble matter in a pipe flow of blood

Summary An analysis of the effects of couple-stresses on the effective Taylor diffusion coefficient has been carried out with the help of two non-dimensional parameters based on the concentration of suspensions and , a constant associated with the couple-stresses. It is observed that the concentration distribution increases with increasing or The effective Taylor diffusion coefficient falls rapidly with increasing when is negative.

---

Zusammenfassung Der Einfluß der Momentenspannungen auf den effektiven Taylorschen Diffusionskoeffizienten wird untersucht. Dabei treten zwei dimensionslose Parameter and auf: Der erste bezieht sich auf die Suspensionskonzentration, der zweite kennzeichnet die Momentenspannungen. Man findet, daß die Verteilungsgeschwindigkeit mit wachsendem oder zunimmt. Dagegen fällt der Taylorsche Diffusionskoeffizient bei wachsendem stark ab, wenn negativ ist.

---

a Tube radius - C Concentration - C _i Body moment vector - C ₀ Concentration at the axis of the tube - C _m Mean concentration - D Molecular diffusion coefficient - d _ij Symmetric part of velocity gradient - F Function of and characterising effective Taylor diffusion coefficient - f _i Body force vector - H A function of and - K ₂ Integration constant - K ^* Effective Taylor diffusion coefficient - k Radius of gyration of a unit cuboid with its sides normal to the spatial axes - I _n Modified Bessel's function ofnth order - L Length of the tube over which the concentration is spread - M Function ofH and - M _ij Couple stress tensor - P Function of - p Fluid pressure - Q Volume rate of the transport of the solute across a section of the tube - r Radial distance from the axis of the tube - T _ij Stress tensor - t Time coordinate - T _ij ^A Antisymmetric part of the stress tensor - u Relative fluid velocity - Average velocity - v _i Velocity vector - Fluid velocity at any point of the tube - v ₀ ⁿ Velocity of Newtonian fluid at the axis of the tube - _i Vorticity vector - x Axial coordinate - x ₁ Relative axial coordinate - z Non-Dimensional radial coordinate - Density - _ij Symmetric part of the stress tensor - µ Viscosity of the fluid - µ _ij Deviatoric part ofM _ij - , Constants associated with couple-stress With 3 figures     

Response of a viscous annular liquid layer in zero-gravity to different axial excitations of the rigid boundary plates

A cylindrical annular liquid layer between two plates and around a rigid center-core consisting of incompressible and viscous liquid is subjected to different axial excitations, such as one-sided, counter-directional and double-sided unequal excitations. The response of the free liquid surface, the velocity- and pressure-distribution has been determined.

---

Zusammenfassung Eine zylindrische Flüssigkeitsschicht bestehend aus inkompressibler und viskoser Flüssigkeit wurde verschiedenen harmonischen Anregungsformen ausgesetzt. Dabei wurden die Fälle einseitiger, doppelseitiger entgegengesetzter und ungleicher doppelseitiger Anregung mit Phase behandelt. Die Vergrößerungsfunktionen für die freie Flüssigkeitsoberfläche, für die Geschwindigkeits- und Druckverteilung wurden bestimmt.

---

List of symbols a radius of liquid layer - b radius of inner cylindrical core - (a–b) thickness of layer - e _r , e , k unit vectors in the radial, angular and axial direction resp. - h length of layer - I _m , K _m modified Bessel functions of first and second kind and order m - diameter ratio - p pressure - q 2na/h - q* na/h - r, , z cylindrical coordinates - complex frequency - S sa ²/ - t time - u, w velocity components in the radial- and axial direction - ₀ excitation amplitude - abbreviation - surface tension parameter - surface tension - dynamic viscosity - kinematic viscosity - density of liquid - free liquid surface elevation - dimensionless time - _rz shear stress - reduced forcing frequency - forcing frequency - stream function - _mn natural frequency of non-viscous liquid     

Instability of Buckley-Leverett flow in a heterogeneous medium

We study the simultaneous one-dimensional flow of water and oil in a heterogeneous medium modelled by the Buckley-Leverett equation. It is shown both by analytical solutions and by numerical experiments that this hyperbolic model is unstable in the following sense: Perturbations in physical parameters in a tiny region of the reservoir may lead to a totally different picture of the flow. This means that simulation results obtained by solving the hyperbolic Buckley-Leverett equation may be unreliable.Symbols and Notation f fractional flow function varying withs andx - value off outsideI - value off insideI - local approximation off around¯x - f ^–,f ⁺ values of - f _j ⁿ value off atS _j ⁿ andx _j - g acceleration due to gravity [ms^–2] - I interval containing a low permeable rock - k dimensionless absolute permeability - k ^* absolute permeability [m²] - k _c ^* characteristic absolute permeability [m²] - k _ro relative oil permeability - k _rw relative water permeability - L ^* characteristic length [m] - L ₁ the space of absolutely integrable functions - L the space of bounded functions - P _c dimensionless capillary pressure function - P _c ^* capillary pressure function [Pa] - P _c ^* characteristic pressure [Pa] - S similarity solution - S _j ⁿ numerical approximation tos(x_j, t_n) - S ₁, S₂,S ₃ constant values ofs - s water saturation - value ofs at - s ^L left state ofs (wrt. ) - s ^R right state ofs (wrt. ) - s s for a fixed value of in Section 3 - T value oft - t dimensionless time coordinate - t ^* time coordinate [s] - t _c ^* characteristic time [s] - t _n temporal grid point,t _n=n t - v ^* total filtration (Darcy) velocity [ms^–1] - W, , v dimensionless numbers defined by Equations (4), (5) and (6) - x dimensionless spatial coordinate [m] - x ^* spatial coordinate [m] - x _j spatial grid piont,x _j=j x - discontinuity curve in (x, t) space - right limiting value of¯x - left limiting value of¯x - angle between flow direction and horizontal direction - t temporal grid spacing - x spatial grid spacing - length ofI - parameter measuring the capillary effects - argument ofS - _o dimensionless dynamic oil viscosity - _w dimensionless dynamic water viscosity - _c ^* characteristic viscosity [kg m^–1s^–1] - _o ^* dynamic oil viscosity [kg m^–1s^–1] - _w ^* dynamic water viscosity [k gm^–1s^–1] - _o dimensionless density of oil - _w dimensionless density of water - _c ^* characteristic density [kgm^–3] - _o ^* density of oil [kgm^–3] - _w ^* density of water [kgm^–3] - porosity - dimensionless diffusion function varying withs andx - ^* dimensionless function varying with s andx ^* [kg^–1m³s] - _j ⁿ value of atS _j ⁿ andx _j This research has been supported by VISTA, a research cooperation between the Norwegian Academy of Science and Letters and Den norske stats oljeselskap a.s. (Statoil).     

Rheological properties of reactive polymer melts.

A new method for describing the rheological properties of reactive polymer melts, which was presented in an earlier paper, is developed in more detail. In particular, a detailed derivation of the equation of a first-order rheometrical flow surface is given and a procedure for determining parameters and functions occurring in this equation is proposed. The experimental verification of the presented approach was carried out using our data for polyamide-6.Notation E Dimensionless reduced viscosity, eq. (34) - E ₀ Newtonian asymptote of the function (36) - E power-law asymptote of the function (36) - E _{= 1} the value ofE at = 1 - k degradation reaction rate constant, s^–1 - k ₁ rate constant of function (t), eq. (26), s^–1 - k ₂ rate constant of function (t), eq. (29), s^–1 - K(t) residence-time-dependent consistency factor, eq. (22) - M _w weight-average molecular weight - M _x x-th moment of the molecular weight distribution - R gas constant - S _x M _x /M _w - t residence time in molten state, s - t _j thej-th value oft, s - T temperature, K - % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xd9vqpe0x% c9q8qqaqFn0dXdir-xcvk9pIe9q8qqaq-xir-f0-yqaqVeLsFr0-vr% 0-vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaaieGaceWFZo% Gbaiaaaaa!3B4E!\[\dot \gamma \] shear rate, s^–1 - _i thei-th value of , s^–1 - _r =1 the value of at = 1, s^–1 - ^* reduced shear rate, eq. (44), s^–1 - dimensionless reduced shear rate, eq. (35) - viscosity, Pa · s - shear-rate and residence-time dependent viscosity, Pa · s - zero-shear-rate degradation curve - degradation curve at - _{t0 (t)} zero-residence-time flow curve - Newtonian asymptote of the RFS - instantaneous flow curve - power-law asymptote of the RFS - _0,0 zero-shear-rate and zero-residence-time viscosity, Pa · s - _E=1 value of viscosity atE=1, Pa · s - ^* reduced viscosity, eq. (43), Pa · s - zero-residence-time rheological time constant, s - density, kg/m³ - (t),(t) residence time functions     

The flow regime in the turbulent near wake of a hemisphere

The work presented is a wind tunnel study of the near wake region behind a hemisphere immersed in three different turbulent boundary layers. In particular, the effect of different boundary layer profiles on the generation and distribution of near wake vorticity and on the mean recirculation region is examined. Visualization of the flow around a hemisphere has been undertaken, using models in a water channel, in order to obtain qualitative information concerning the wake structure.List of symbols C _p pressure coefficient, - D diameter of hemisphere - n vortex shedding frequency - p pressure on model surface - p ₀ static pressure - Re Reynolds number, - St Strouhal number, - U, V, W local mean velocity components - mean freestream velocity inX direction - U ^* shear velocity, - u, v, w velocity fluctuations inX, Y andZ directions - X Cartesian coordinate in longitudinal direction - Y Cartesian coordinate in lateral direction - Z Cartesian coordinate in direction perpendicular to the wall - it* boundary layer displacement thickness, - diameter of model surface roughness - elevation angleI - O boundary layer momentum thickness, - _w wall shearing stress - dynamic viscosity of fluid - density of fluid - streamfunction - _x longitudinal component of vorticity, - _y lateral component of vorticity, - _z vertical component of vorticity, This paper was presented at the Ninth symposium on turbulence, University of Missouri-Rolla, October 1–3, 1984     

Zur analytischen Beschreibung des Füllstoffeinflusses auf das Fließverhalten von Kunststoff-Schmelzen

From literature some representative equations have been compiled describing the influence of filler on the viscosity of polymer melts. By application of these on the experimental results obtained from GF-SAN it was found that the relative viscosity _R , i.e. the ratio of the viscosities of the filled and unfilled melt, shows a pronounced dependence on the shear rate but not on the shear stress. Defining _R with constant and not with constant (as it is usually done), an analytical approach is possible independent of Further the influence of pressure, temperature and filler content on the zero-shear viscosity of filled polymer melts may be expressed by a modified Arrhenius equation.

     

Velocity field in isothermal turbulent bubbly gas-liquid flow through a pipe

Velocity field was measured by laser Doppler velocimetry in isothermal, turbulent bubbly gas-liquid flow through a 26.6 mm inner diameter vertical pipe. The measurements were made about 33 diameters downstream from the pipe entrance, gas injection being just upstream of the entrance. The gas phase radial distribution at the measurement plane exhibited influence of the injection device in that higher gas fraction existed in the central region of the pipe. For comparison, velocity field was also measured in isothermal, turbulent single-phase liquid flow through the same pipe at the same axial plane. Measured were the radial distributions of liquid mean axial and radial velocities, axial and radial turbulent intensities, and axial Reynolds shear stress. The radial distributions of gas bubble mean axial velocity and axial velocity fluctuation intensity were also measured by LDV. A dualsensor fiberoptic probe was used at the same time to measure the radial distributions of gas fraction, bubble mean axial velocity and size slightly downstream of the LDV measurement plane.List of Symbols an average gas bubble diameter - f, f _TP friction factor, friction factor for gas-liquid flow - k _L liquid turbulent kinetic energy - , gas, liquid mass flow rate - R inner radius of pipe - r, {sitR}^* radial coordinate; nondimensional radial coordinate (=r/R) - Re _L liquid Reynolds number - U _G mean axial velocity of gas bubble - U _L mean axial velocity of liquid - U _LO mean axial velocity for flow at the total mass velocity with properties of the liquid phase - u _L ⁺ nondimensional mean axial velocity of liquid in wall coordinate - friction velocity - axial velocity fluctuation intensity of liquid - axial velocity fluctuation intensity of gas bubbles - V_L mean radial velocity of liquid - v _L radial velocity fluctuation intensity of liquid - (uv)_L single-point cross-correlation between axial and radial velocity fluctuations of liquid ( axial Reynolds shear stress) - T _in mean liquid temperature at test section inlet - x flow quality - y normal distance from wall - y ⁺ nondimensional normal distance from wall in wall coordinate (=yu/v_L) - _G gas phase residence time fraction - _L rate of dissipation in the liquid - _L Kolmogorov length scale in the liquid - _L liquid kinematic viscosity - _L characteristic turbulence length scale in the liquid - _G, _L density of gas, liquid - _m gas-liquid mixture density This work was partly supported by National Science Foundation, Thermal Transport and Thermal Processing Program, Chemical and Thermal Systems Division, under Grant No. CTS-9411898.     

Pore structures and transport properties of sandstone

To quantitatively analyze the macroscopic properties of the flow in porous media by means of the continuum approach, detailed information (velocity and pressure fields) on the microscopic scale is necessary. In this paper, the numerical solution for incompressible, Newtonian flow in a diverging-converging representative unit cell (RUC) is presented. A new solution procedure for the problem is introduced. A review of the accuracy of the computational method is given.Nomenclature A _ff ^* area of entrance and exit of RUC - A _fs ^* interfacial area between the fluid and solid phases - d throat diameter of RUC (m) - D pore diameter of RUC (m) - i, j unit vector for RUC - L ^* wave length of a unit cell - L _p pore length of RUC (m) - L _t throat length of RUC (m) - n unit outwardly directed vector for the fluid phase - p ^* fluid pressure - ^* cross-sectional mean pressure - _en ^* entrance cross-sectional mean pressure - Re _d Reynolds number - x ^*, r^* cylindrical coordinates - u ^*, v^* velocity - u _cl ^* centerline velocity - _d mean velocity at the throat of RUC (m/s) - _D mean velocity at the large segment of RUC (m/s) Greek viscosity coefficient (Ns/m²) - _p excess momentum loss factor defined in (4.1) - fluid density (kg/m³) - ^* stream function - ^* vorticity - dimensionless circulation defined in (2.7) Symbols - the mean value - * dimensionless quantities     

Solutions for steady groundwater flow in multi-layer aquifer systems

A method is shown that enables closed-form solutions to be calculated for steady flow in leaky multi-layer aquifer systems. The method requires, as a prerequisite, the numerical solution of a generalized eigenvalue problem. This eigenvalue problem always has real positive eigenvalues, and the eigenvectors are always orthogonal relative to the transmissivity matrix. Partial solutions are given for a number of examples, and a numerical example is used to show how a solution can be obtained for flow to a well when the well abstracts water from more than one aquifer.List of notation A area - , specified vector functions - piezometric head vector - h _i piezometric head in aquifer i - I ₀, I ₁ modified Bessel functions of the first kind - K ₀ , K ₁ modified Bessel functions of the second kind - l half width of an irrigated strip - vertical flux velocity vector - P _i vertical flux velocity in aquifer i - well flow rate vector - Q _i well flow rate in aquifer i - R leakance matrix - R _i leakance of aquitard i - R radius of an irrigated area - r distance between two points - r ₀ well screen radius - T transmissivity matrix - tt_i transmissivity of acquifer i - _n eigenvector - Euler's constant - _n eigenvalue - v arc length normal to a boundary     

Zum Mischen rheologisch inhomogener Stoffsysteme auf Einschneckenmaschinen

Zusammenfassung Das Mischen von Stoffen mit unterschiedlichen rheologischen Eigenschaften in Schneckenmaschinen ist in der Kunststoffauf- und -verarbeitung eine Standardaufgabe. Trotzdem gibt es hierfür kein zufriedenstellendes mathematisch-physikalisches Modell. Daher werden zunächst einfache Mischmodelle diskutiert. Auf der Basis dieser Modelle wird dann unter Berücksichtigung der Besonderheiten des Plastifizierextruderprozesses eine Mischgütebeziehung mathematisch formuliert. Die experimentelle Überprüfung erfolgt mit Hilfe der Grauwertanalyse extrudierter Zweistoffsysteme, bei denen ein Stoff mit Ruß eingefärbt war. Da der Mischprozeß hochgradig stochastisch ist, streuen die Meßergebnisse. Unter Berücksichtigung dieses Tatbestandes ist der theoretische Ansatz zufriedenstellend.

---

Mixing of polymer resins with different rheological properties is a usual demand in plastics processing using screw extruders. A mathematical model describing this processing problem sufficiently is not known, however. Therefore, simple mixing models will be discussed. Based on these, a concept for the calculation of mixing homogeneity will be presented, including the particular requirement of the plasticating screw process. An experimental investigation utilizes the grey-value analysis of extruded two-component materials, which in one phase is carbon-black filled. Considering the fact that the mixing process is highly random, the theoretical model leads to a good level of aggreement with the scattering measurement data.

---

b Schneckenkanalbreite - B Bandbreite der Grauwerte - c Konstante - mittlere Konzentration, bezogen auf die Grauwertbandbreite - h Höhe, Gangtiefe, Schneckenkanalhöhe - h ₀ Gangtiefe der Einzugszone - h ₁ Gangtiefe der Ausstoßzone - L Länge - gemittelte Schmelzebettlänge - n Exponent des Potenzfließgesetzes - s Standardabweichung der Grauwerte bezogen auf die Grauwertbandbreite - S Standardabweichung der Grauwerte - t Verweilzeit - t ₁ kürzeste Verweilzeit - mittlere Verweilzeit - ₀ Umfangsgeschwindigkeit - mittlere Geschwindigkeit - V Volumenstrom - w Dicke eines Kontrollelements - w Ausstreichdicke eines Kontrollelements - x Koordinate - Mittelwert der Grauwerte - y Koordinate - Scherdeformationswinkel - Scherdeformation - mittlere Scherdeformation - Schergeschwindigkeit - Viskosität - ₁ dimensionslose kürzeste Verweilzeit - dimensionsloser Volumenstrom - _LSM laminarer Schermischgrad - _{LSM, the} theoretischer laminarer Schermischgrad - _{LSM, exp} experimenteller laminarer Schermischgrad - ² Varianz der Verweilzeit im Schmelzebett - Schubspannung - Gangsteigungswinkel der Schnecke - ø Volumenanteil - dimensionslose Kennzahl     

Modeling of contaminant transport resulting from dissolution of nonaqueous phase liquid pools in saturated porous media

A mathematical model for transient contaminant transport resulting from the dissolution of a single component nonaqueous phase liquid (NAPL) pool in two-dimensional, saturated, homogeneous porous media was developed. An analytical solution was derived for a semi-infinite medium under local equilibrium conditions accounting for solvent decay. The solution was obtained by taking Laplace transforms to the equations with respect to time and Fourier transforms with respect to the longitudinal spatial coordinate. The analytical solution is given in terms of a single integral which is easily determined by numerical integration techniques. The model is applicable to both denser and lighter than water NAPL pools. The model successfully simulated responses of a 1,1,2-trichloroethane (TCA) pool at the bottom of a two-dimensional porous medium under controlled laboratory conditions.Notation a,a ₁ defined in (45a) and (45b), respectively - b defined in (45c) - b vector of true model parameters (n×1) - vector of estimated model parameters (n×1) - c liquid phase solute concentration (solute mass/liquid volume), M/L³ - c _s aqueous saturation concentration (solubility), M/L³ - C dimensionless liquid phase solute concentration, equal toc/c _s - molecular diffusion coefficient, L²/t - ^e effective molecular diffusion coefficient, equal to / ^*, L²/t - D _x longitudinal hydrodynamic dispersion coefficient, L²/t - D _z hydrodynamic dispersion coefficient in the vertical direction, L²/t - e random vector with zero mean (m×1) - erf[x] error function, equal to (2/ ^1/2) - f vector of fitting errors or residuals (m×1) - Fourier operator - _-1 Fourier inverse operator - g vector of model simulated data (m×1) - k mass transfer coefficient, L/t - average mass transfer coefficient, L/t - K _d partition or distribution coefficient (liquid volume/solids mass), L³/M - pool length, L - _o distance between the pool and the origin of the specified Cartesian coordinate system, L - Laplace operator - _-1 Laplace inverse operator - m number of observations - M Laplace/Fourier function defined in (38) - n number of model parameters - N Laplace/Fourier function defined in (39) - p defined in (46) - Pe _x Péclet number, equal toU _x /D _x - Pe _z Péclet number, equal toU _x /D _z - q defined in (47) - R retardation factor - s Laplace transform variable - S objective function - Sh local Sherwood number, equal tok/ _e - Sh _o overall Sherwood number, equal to l/ _e - t time,t - T dimensionless time, equal toU _x t/ - u dummy integration variable - u vector of independent variables - U _x average interstitial velocity, L/t - x spatial coordinate in the longitudinal direction, L - X dimensionless longitudinal length, equal to (x–)/ - y vector of observed data (m×1) - z spatial coordinate in the vertical direction, L - Z dimensionless vertical length, equal toz/ - Fourier transform variable - defined in (37) - defined in (50) - porosity (liquid volume/aquifer volume), L³/L³ - defined in (52a) and (52b), respectively - decay coefficient, t^–1 - dimensionless decay coefficient, equal to /U _x - bulk density of the solid matrix (solids mass/aquifer volume), M/L³ - dummy integration variable - ^* tortuosity     

Production of temperature fluctuations in grid turbulence: Wiskind's experiment revisited

Measurements have been made in nearly-isotropic grid turbulence on which is superimposed a linearly-varying transverse temperature distribution. The mean-square temperature fluctuations, , increase indefinitely with streamwise distance, in accordance with theoretical predictions, and consistent with an excess of production over dissipation some 50% greater than values recorded in previous experiments. This high level of production has the effect of reducing the ratio,r, of the time scales of the fluctuating velocity and temperature fields. The results have been used to estimate the coefficient,C, in Monin's return-to-isotropy model for the slow part of the pressure terms in the temperature-flux equations. An empirical expression by Shih and Lumley is consistent with the results of earlier experiments in whichr 1.5, C 3.0, but not with the present data where r 0.5, C 1.6. Monin's model is improved when it incorporates both time scales.List of symbols C coefficient in Monin model, Eq. (5) - M grid mesh length - m exponent in power law for temperature variance, x ^m - n turbulence-energy decay exponent,q ² x ^-n - p production rate of - p pressure - q ² - R microscale Reynolds number - r time-scale ratiot/t - T mean temperature - U mean velocity - mean-square velocity fluctuations (turbulent energy components) - turbulent temperature flux - x, y, z spatial coordinates - temperature gradient dT/dy - thermal diffusivity - dissipation rate ofq ²/2 - dissipation rate of - Taylor microscale (₂=5q₂/) - temperature microscale - _v temperature-flux correlation coefficient, /v - dimensionless distance from the grid,x/M     

The microscopic analysis of high forchheimer number flow in porous media

High Forchheimer number flow through a rigid porous medium is numerically analysed by means of the volumetric averaging concept. The microscopic flow mechanisms, which must be known in order to understand the macroscopic flow phenomena, are studied by utilising a periodic diverging-converging representative unit cell (RUC). The detailed information for the microscopic flow field, in association with the locally averaged momentum balance, makes it possible to quantitatively demonstrate that the microscopic inertial phenomenon, which leads to distorted velocity and pressure fields, is the fundamental reason for the onset of nonlinear (non-Darcy) effects as velocity increases. The hydrodynamic definitions for Darcy's law permeabilityk, the inertial coefficient and Forchheimer number Fo are obtained by applying the averaging theorem to the pore level Navier-Stokes equations. Finally, these macroscopic parameters are numerically calculated at various combinations of micro-geometry and flow rate, and graphically correlated with the relevant microscopic parameters.Nomenclature a _i body force acceleration (m/s²) - A viscous integral term defined in (4.6) - A _f area of entrance and exist of RUC (m²) - A _fs interfacial area between the fluid and solid phases (m²) - B pressure integral term defined in (4.4) - d throat diameter of RUC (m) - D pore diameter of RUC (m) - Fo Forchheimer number defined in (4.1) and (4.10) - g gravitational acceleration (m/s²) - i, j microscopic unit vector for RUC - k Darcy's law permeability (m²) - k _v velocity dependent permeability defined in (4.1) (m²) - L length of a unit cell (m) - L _p pore length of RUC (m) - L _t throat length of RUC (m) - n unit outwardly directed vector for the fluid phase - p microscopic fluid pressure (N/m²) - P macroscopic fluid pressure (N/m²) - _en mean pressure at entrance of RUC (N/m²) - _ex mean pressure at exit of RUC (N/m²) - r _i,r coordinate on the macroscopic scale (m) - Re _d Reynolds number defined in (4.5) - u _i,u microscopic velocity (m/s) - specific discharge (m/s) - _d mean velocity at the throat of RUC (m/s) - v microscopic velocity (m/s) - V _b representative elementary volume (REV) (m³) - V _f volume occupied by the fluid within REV (m³) - V _s volume occupied by the solid within REV (m³) - x _i,x coordinate on the microscopic scale (m) - X _i,X coordinate on the macroscopic scale (m) Greek the inertia coefficient (1/m) - viscosity coefficient (Ns/m²) - _i microscopic unit vector - areosity at the entrance and the exit cross-section of RUC - fluid density (kg/m³) - porosity - _f a general property of the fluid phase Symbols ^f intrinsic phase average - the fluctuating part of _f - the mean value of _f - _f ^* the dimensionless value of _f     

Unsteady flow in an annulus between two concentric rotating spheres

An analysis is presented for the unsteady laminar flow of an incompressible Newtonian fluid in an annulus between two concentric spheres rotating about a common axis of symmetry. A solution of the Navier-Stokes equations is obtained by employing an iterative technique. The solution is valid for small values of Reynolds numbers and acceleration parameters of the spheres. In applying the results of this analysis to a rotationally accelerating sphere, a virtual moment of intertia is introduced to account for the local inertia of the fluid.Nomenclature R _i radius of the inner sphere - R _o radius of the outer sphere - radial coordinate - r dimensionless radial coordinate, - meridional coordinate - azimuthal coordinate - time - t dimensionless time, - Re _i instantaneous Reynolds number of the inner sphere, _i R _k ² / - Re _o instantaneous Reynolds number of the outer sphere, _o R _o ² / - radial velocity component - V _r dimensionless radial velocity component, - meridional velocity component - V dimensionless meridional velocity component, - azimuthal velocity component - V dimensionless azimuthal velocity component, - viscous torque - T dimensionless viscous torque, - viscous torque at surface of inner sphere - T _i dimensionless viscous torque at surface of inner sphere, - viscous torque at surface of outer sphere - T _o dimensionless viscous torque at surface of outer sphere, - externally applied torque on inner sphere - T _p,i dimensionless applied torque on inner sphere, - moment of inertia of inner sphere - Z _i dimensionless moment of inertia of inner sphere, - virtual moment of inertia of inner sphere - Z _i,v dimensionless virtual moment of inertia of inner sphere, - virtual moment of inertia of outer sphere - _i instantaneous angular velocity of the inner sphere - _o instantaneous angular velocity of the outer sphere - density of fluid - viscosity of fluid - kinematic viscosity of fluid,/ - radius ratio,R _i/R _o - swirl function, - dimensionless swirl function, - stream function - dimensionless stream function, - _i acceleration parameter for the inner sphere, - _o acceleration parameter for the outer sphere, - shear stress - _r dimensionless shear stress,      

Stochastic response of structures with small geometric imperfections

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

---

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

---

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

Flow of second-order fluids over an enclosed rotating disc

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

A generalized rheological model of thixotropic materials

Summary A generalization of the rheological model of thixotropic materials, presented previously, was carried out. In the generalized rheological equation of state the yield stress depending on the structural parameter was introduced. In the generalized rate equation the difference in the destruction and recovery rates of the material structure was taken into account. A procedure leading to the determination of nine rheological parameters of the generalized model was worked out. The model was checked experimentally for a thixotropic paint.

---

Zusammenfassung Eine früher dargestellte Theorie thixotroper Stoffe wird verallgemeinert, wobei eine von dem Strukturparameter abhängige Fließspannung eingeführt wird. Weiterhin wird der Unterschied zwischen der Zerstörungs-und der Wiederaufbaugeschwindigkeit der Stoffstruktur berücksichtigt. Eine Methode zur Bestimmung der neun benötigten Stoffparameter wird ausgearbeitet. Das Modell wird am Beispiel einer thixotropen Farbe experimentell geprüft.

---

Notation a rheological parameter in eq. [26], s^–1 - A rheological parameter in eq. [16] - b rheological parameter in eq. [26] - c function in eq. [21] - averaged value of functionc in eq. [28] - c function in the rate equation [23], defined by eq. [21] - G function [1] defining material of the rate type - h function [2] determining the state of thixotropic fluid - k rheological parameter in the Herschel-Bulkley equation [17] or, in special case, in eq. [8], Ns ⁿ /m² - K function in eq. [18], Ns ^m /m² - m rheological parameter in eq. [18] or, in special case, in eq. [10] - n rheological parameter in the Herschel-Bulkley model [17] or, in special case, in eq. [8] - s rheological parameter in eq. [16] - t time, s - x arbitrary real variable - rheological parameter in eq. [9], s - shear rate, s^–1 - structural parameter, defined by eq. [2] - substantial derivative of structural parameter, s^–1 - _e function [6] describing the equilibrium curve in the coordinate system ( ) - ₀ initial value of structural parameter (att = 0) - natural time function of the thixotropic material, defined by eq. [22] - shear stress, N/m² - substantial derivative of shear stress, N/m² s - _e function describing equilibrium flow curve in the coordinate system ( ) - ₀ equilibrium yield stress, defined by eq. [12], N/m² - _y function of structural parameter describing the yield stress - function in eq. [11] Notation used in the algorithm:(Appendix) i,j,k integer - k _e (i) ordinal number of the experimental point at which the line of _i = const intersects the equilibrium flow curve - l _i number of the experiments of the type stepchange of the shear rate - l _j number of experimental points in one experiment of the type step-change of the shear rate - n _e number of experimental points on the equilibrium flow curve - n _k number of experimental points on the line of constant - n _y number of lines of constant - t(j) measured time interval (from the moment of the step-change of shear rate) - abscissa of the experimental point of ordinal numberk on the line of _i = const, in the coordinate system ( ) - abscissa of the experimental point of ordinal numberi on the equilibrium flow curve, in the coordinate system ( ) - shear rate at which the experiment of the type step-change of shear rate was carried out - _e (i) ordinate of the experimental point of ordinal numberi on the equilibrium flow curve, in the coordinate system ( ) - _y (i) value of yield stress at = _i - _s (i,j) experimental value of shear stress at constant value of shear rate (2i) for time intervalt(j) - (i,k) ordinate of the experimental point of ordinal numberk on the line of _i = const, in the coordinate system ( ) - ₀ the admissible value of the difference between the experimental and theoretical value of shear stress With 4 figures and 1 table