Nonlinear rheological behavior of a concentrated spherical silica suspension

Linear and nonlinear viscoelastic properties were examined for a 50 wt% suspension of spherical silica particles (with radius of 40 nm) in a viscous medium, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The effective volume fraction of the particles evaluated from zero-shear viscosities of the suspension and medium was 0.53. At a quiescent state the particles had a liquid-like, isotropic spatial distribution in the medium. Dynamic moduli G^* obtained for small oscillatory strain (in the linear viscoelastic regime) exhibited a relaxation process that reflected the equilibrium Brownian motion of those particles. In the stress relaxation experiments, the linear relaxation modulus G(t) was obtained for small step strain (0.2) while the nonlinear relaxation modulus G(t, ) characterizing strong stress damping behavior was obtained for large (>0.2). G(t, ) obeyed the time-strain separability at long time scales, and the damping function h() (–G(t, )/G(t)) was determined. Steady flow measurements revealed shear-thinning of the steady state viscosity () for small shear rates (< ^–1; = linear viscoelastic relaxation time) and shear-thickening for larger (>^–1). Corresponding changes were observed also for the viscosity growth and decay functions on start up and cessation of flow, + (t, ) and _– (t, ). In the shear-thinning regime, the and dependence of ₊(t,) and _–(t,) as well as the dependence of () were well described by a BKZ-type constitutive equation using the G(t) and h() data. On the other hand, this equation completely failed in describing the behavior in the shear-thickening regime. These applicabilities of the BKZ equation were utilized to discuss the shearthinning and shear-thickening mechanisms in relation to shear effects on the structure (spatial distribution) and motion of the suspended particles.Dedicated to the memory of Prof. Dale S. Parson     

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

Similar solutions and the asymptotics of filtration equations

In this paper we consider the asymptotic behavior of solutions of the quasilinear equation of filtration as t. We prove that similar solutions of the equation u _t = (u )_xx asymptotically represent solutions of the Cauchy problem for the full equation u _t = [(u)]_xx if (u) is close to u for small u.     

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)     

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.     

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.

     

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

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

The dynamic and steady shear properties of synovial fluid and of the components making up synovial fluid

Steady-shear and dynamic properties of a pooled sample of cattle synovial fluid have been measured using techniques developed for low viscosity fluids. The rheological properties of synovial fluid were found to exhibit typical viscoelastic behaviour and can be described by the Carreau type A rheological model. Typical model parameters for the fluid are given; these may be useful for the analysis of the complex flow problems of joint lubrication.The two major constituents, hyaluronic acid and proteins, have been successfully separated from the pooled sample of synovial fluid. The rheological properties of the hyaluronic acid and the recombined hyaluronic acid-protein solutions of both equal and half the concentration of the constituents found in the original synovial fluid have been measured. These properties, when compared to those of the original synovial fluid, show an undeniable contribution of proteins to the flow behaviour of synovial fluid in joints. The effect of protein was found to be more prominent in hyaluronic acid of half the normal concentration found in synovial fluid, thus providing a possible explanation for the differences in flow behaviour observed between synovial fluid from certain diseased joints compared to normal joint fluid.Nomenclature A Ratio of angular amplitude of torsion head to oscillation input signal - G Storage modulus - G Loss modulus - I Moment of inertia of upper platen — torsion head assembly - K Restoring constant of torsion bar - N ₁ First normal-stress difference - R Platen radius - S ⁽ⁱ⁾ Geometric factor in the dynamic property analysis - t ₁ Characteristic time parameter of the Carreau model - X, Y Carreau model parameters - Z () Reimann Zeta function of - Carreau model parameter - Shear rate - Apparent steady-shear viscosity - ^* Complex dynamic viscosity - Dynamic viscosity - Imaginary part of the complex dynamic viscosity - ₀ Zero-shear viscosity - ₀ Cone angle - Carreau model characteristic time - Density of fluid - Shear stress - Phase difference between torsion head and oscillation input signals - ₀ Zero-shear rate first normal-stress coefficient - Oscillatory frequency     

Material functions generated by the complex viscosity

Based on the complex viscosity model various steady-state and transient material functions have been completed. The model is investigated in terms of a corotational frame reference. Also, BKZ-type integral constitutive equations have been studied. Some relations between material functions have been derived. C ^–1 Finger tensor - F[], (F ^–1[]) Fourier (inverse) transform - rate of deformation tensor in corotating frame - h(I, II) Wagner's damping function - J (x) Bessel function - m parameter inh (I, II) - m(s) memory function - m _k, n_k integers (powers in complex viscosity model) - P principal value of the integral - parameter in the complex viscosity model - rate of deformation tensor - shear rates - [], [] incomplete gamma function - (a) gamma function - steady-shear viscosity - ^* complex viscosity - , real and imaginary parts of ^* - ₀ zero shear viscosity - ⁺, ₁ ⁺ stress growth functions - ^–, ₁ ^- stress relaxation functions - (s) relaxation modulus - ₁(s) primary normal-stress coefficient - ø(a, b; z) degenerate hypergeometric function - ₁, ₂ time constants (parameters of ^*) - frequency - extra stress tensor     

Pulsatile blood flow through arterioles

An analytical study was made to examine the effect of vascular deformability on the pulsatile blood flow in arterioles through the use of a suitable mathematical model. The blood in arterioles is assumed to consist of two layers — both Newtonian but with differing coefficients of viscosity. The flow characteristics of blood as well as the resistance to flow have been determined using the numerical computations of the resulting expressions. The applicability of the model is illustrated using numerical results based on the existing experimental data. r, z coordinate system - u, axial/longitudinal velocity component of blood - p pressure exerted by blood - _b density of blood - µ viscosity of blood - t time - , displacement components of the vessel wall - T _t0,T ₀ known initial stresses - density of the wall material - h thickness of the vessel wall - T _t,T stress components of the vessel - K _l,K _r components of the spring coefficient - C _l,C _r components of the friction coefficient - M _a additional mass of the mechanical model - r ₁ outer radius of the vessel - thickness of the plasma layer - r ₁ inner radius of the vessel - circular frequency of the forced oscillation - k wave number - E ₀,E _t, , _t material parameters for the arterial segment - µ _p viscosity of the plasma layer - Q total flux - Q _p flux across the plasma zone - Q _h flux across the core region - Q mean flow rate - resistance to flow - P pressure difference - l length of the segment of the vessel     

A noninvasive method for the measurement of flow-induced surface displacement of a compliant surface

A noninvasive optical method is described which allows the measurement of the vertical component of the instantaneous displacement of a surface at one or more points. The method has been used to study the motion of a passive compliant layer responding to the random forcing of a fully developed turbulent boundary layer. However, in principle, the measurement technique described here can be used equally well with any surface capable of scattering light and to which optical access can be gained. The technique relies on the use of electro-optic position-sensitive detectors; this type of transducer produces changes in current which are linearly proportional to the displacement of a spot of light imaged onto the active area of the detector. The system can resolve displacements as small as 2 m for a point 1.8 mm in diameter; the final output signal of the system is found to be linear for displacements up to 200 m, and the overall frequency response is from DC to greater than 1 kHz. As an example of the use of the system, results detailing measurements obtained at both one and two points simultaneously are presented.List of symbols C _t elastic transverse wave speed = (G/)^1/2 - d ⁺ spot diameter normalized by viscous length scale - G frequency average of G() - G() shear storage modulus - G() shear loss modulus - l. viscous length scale = v/u _* - N total number of sampled data values - r separation vector for 2-point measurements = (, ) - rms root-mean-square value - R momentum thickness Reynolds number = U _t8/v - t time - u (y) mean streamwise component of velocity in boundary layer - u _* friction velocity = (t _w/)^1/2 - U free-stream velocity - x, y, z longitudinal, normal and spanwise directions - y _o undisturbed surface position - vertical component of compliant surface displacement - ₉₉ boundary layer thickness for which u(y) = 0.99 U _t8 - _l viscous sublayer thickness 5 l _* - frequency average of G()/ - boundary layer momentum thicknes = - fluid dynamic viscosity - v fluid kinematic viscosity = / - , longitudinal, spanwise components of separation vector r - fluid density - time delay - _w wall shear stress     

General viscosity-composition relationship for dispersions,solutions and binary liquid systems

An equation for the viscosity of a mixture of two imaginary Newtonian liquids is derived. In the derivation the mathematical assumption is used that the effective activation energy for viscous flow of a binary liquid mixture is a linear combination of the reciprocals of the activation energy of the components. It contains two dependent fitting constants and has the same structure as the Mooney equation for dispersions of spherical solid particles, the Huggins equation for polymer solutions and is identical to an equation by Hoffmann and Rother, when written in the variables that the last authors used.As a consequence it can be shown that the viscosity of binary liquid mixtures, liquid resion solutions, dispersions of solid spherical particles and polymer solutions can be described very well by one and the same equation, up to the highest concentrations.It has further been found that the viscosity of dispersions of non-spherical particles, solutions of solids in organic solvents and solutions of electrolytes and non-electrolytes in water can also be described by this formula. The equation permits the construction of a straight line on which all liquids can be plotted.An algebraic analysis of the equation shows that each series of viscosity composition data can be placed in one of three rheological groups independent of the type of fraction that is used to characterize the composition.Seventy-four binary systems, covering a wide range of liquids have been used to show the applicability of the developed equation.It has been found that in most cases the data are best described by splitting them into two regions, each with its own set of dependent constants. General symbol for the fraction or concentration of the component with the higher viscosity determining the composition of a binary mixture [—] - _v Volume fraction of the component with the higher viscosity [—] - _w Weight fraction of the component with the higher viscosity [—] - _mw Molecular weight fraction of the component with the higher viscosity [—] - c Concentration of the component with the higher viscosity [g/cm³] - E ₂,E ₁,E Activation energy for viscous flow referring to the component with the higher viscosity, the lower viscosity and the viscosity of the binary mixtures, respectively [J] - ₂, ₁, Experimental parameter (with the dimension of energy) referring to the component with the higher viscosity, the component with the lower viscosity and to the binary mixtures, respectively [J] - ₁, ₂ Viscosity of the component with the lower and the higher viscosity, respectively [Pa · s] - Viscosity of a binary mixture [Pa · s] - [] The usual intrinsic viscosity of the component with the highest viscosity [cm³/g] - _r / ₁ [—] - _{sp r} – 1 [—] - [] -intrinsic viscosity [—] - [] _v Volume intrinsic viscosity [—] - [] _w Weight intrinsic viscosity [—] - [] _c Concentration intrinsic viscosity, identical to [] [cm³/g] - T _e Temperature at which the two liquids have the same viscosity [K] - _e Viscosity at temperatureT _e [Pa · s] - P ₁,P ₂ Density of the component with the lower and the higher viscosity, respectively - R Gas constant [J · Mol^–1 · K^–1]     

Fractal-time stochastic processes and dynamic functions of polymeric systems

Dynamic material functions of polymeric systems are calculated via a defect-diffusion model. The random motion of defects is modelled by a fractaltime stochastic process. It is shown that the dynamic functions of polymeric solutions can be approximated by the defect-diffusion process of the mixed type. The relaxation modulus of Kohlrausch type is obtained for a fractal-time defect-diffusion process, and it is shown that this modulus is capable of portraying the dynamic behavior of typical viscoelastic solutions.The Fourier transforms of the Kohlrausch function are calculated to obtain and. A three-parameter model for and is compared with the previous calculations. Experimental measurements for five polymer solutions are compared with model predictions. D rate of deformation tensor - G(t) mechanical relaxation modulus - H relaxation spectrum - I(t) flux of defects - P _n (s) probability of finding a walker ats aftern-steps - P generating function ofP _n (s) - s(t) fraction of surviving defects - , () gamma function (incomplete) - ₀ zero shear viscosity - ^* () complex viscosity - frequency - t _n n-th moment - F[] Fourier transform - f ^* (u) Laplace transform off(t) - , components of ^* - G _f, _f ^* fractional model - G ₃, ₃ ^* three parameter model - complex conjugate ofz - material time derivative ofD     

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

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

Power spectra of fluid velocities measured by laser Doppler velocimetry

The power spectrum and the correlation of the laser Doppler velocimeter velocity signal obtained by sampling and holding the velocity at each new Doppler burst are studied. Theory valid for low fluctuation intensity flows shows that the measured spectrum is filtered at the mean sample rate and that it contains a filtered white noise spectrum caused by the steps in the sample and hold signal. In the limit of high data density, the step noise vanishes and the sample and hold signal is statistically unbiased for any turbulence intensity.List of symbols A cross-section of the LDV measurement volume, m² - A empirical constant - B bandwidth of velocity spectrum, Hz - C concentration of particles that produce valid signals, number/m³ - d _m diameter of LDV measurement volume, m - f(₁, ₂ | u) probability density of t _i; and t _j given (t) for all t, Hz² - probability density for t _j-t_i, Hz - n (t, t) number of valid bursts in (t, t) = N + n - N (t, t) mean number of valid bursts in (t, t) - N _e mean number of particles in LDV measurement volume - valid signal arrival rate, Hz - mean valid signal arrival rate, Hz - R _uu time delayed autocorrelation of velocity, m²/s² - S _u power spectrum of velocity, m²/s²/Hz - t ₁, t ₂ times at which velocity is correlated, s - t _i, t _j arrival times of the bursts that immediately precede t ₁ and t ₂, respectively, s - t _ij t _j–t _i s - T averaging time for spectral estimator, s - T _u integral time scale of u (t), s - T Taylor's microscale for u (t), s - u velocity vector = U + u, m/s - u fluctuating component of velocity, m/s - U mean velocity, m/s - u _m sampled and held signal, m/s Greek symbols (t) noise signal, m/s - _m (t) sampled and held noise signal, m/s - bandwidth of spectral estimator window, radians/s - time between arrivals in pdf, s - Taylor's microscale of length = UT m - kinematic viscosity - ₁, ₂ arrival times in pdf, s - root mean square of noise signal, m/s - _u root mean square of u, m/s - delay time = t ₂ - t ₁ s - B duration of a Doppler burst, s - circular frequency, radians/s - _c low pass frequency of signal spectrum radians/s Other symbols ensemble average - conditional average - ^ estimate     

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.     

Hyperbolic phenomena in a strongly degenerate parabolic equation

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

Transport in ordered and disordered porous media III: Closure and comparison between theory and experiment

In this paper we examine the closure problem associated with the volume averaged form of the Stokes equations presented in Part II. For both ordered and disordered porous media, we make use of a spatially periodic model of a porous medium. Under these circumstances the closure problem, in terms of theclosure variables, is independent of the weighting functions used in the spatial smoothing process. Comparison between theory and experiment suggests that the geometrical characteristics of the unit cell dominate the calculated value of the Darcy's law permeability tensor, whereas the periodic conditions required for thelocal form of the closure problem play only a minor role.Roman Letters A interfacial area of the- interface contained within the macroscopic region, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A interfacial area of the- interface associated with the local closure problem, m² - A _p surface area of a particle, m² - b vector used to represent the pressure deviation, m^–1 - B ⁰ B+I, a second order tensor that maps v _m ontov - B second-order tensor used to represent the velocity deviation - d _p 6V _p/A_p, effective particle diameter, m - d a vector related to the pressure, m - D a second-order tensor related to the velocity, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor calculated on the basis of a spatially periodic model, m² - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p characteristic length for the volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - characteristic length (pore scale) for the-phase - _i i=1, 2, 3 lattice vectors, m - weighting function - m(-y) , convolution product weighting function - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - n unit normal vector pointing from the-phase toward the -phase - p pressure in the-phase, N/m² - p _m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p – p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function - r position vector, m - r position vector locating points in the-phase, m. - V averaging volume, m³ - B volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - v velocity vector in the-phase, m/s - v _m superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - v traditional superficial volume averaged velocity, m/s - v – v _m , spatial deviation velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the -phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - _m m * , weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

Transient flow towards a well in an aquifer including the effect of fluid inertia

Transient propagation of weak pressure perturbations in a homogeneous, isotropic, fluid saturated aquifer has been studied. A damped wave equation for the pressure in the aquifer is derived using the macroscopic, volume averaged, mass conservation and momentum equations. The equation is applied to the case of a well in a closed aquifer and analytical solutions are obtained to two different flow cases. It is shown that the radius of influence propagates with a finite velocity. The results show that the effect of fluid inertia could be of importance where transient flow in porous media is studied.List of symbols b Thickness of the aquifer, m - c ₀ Wave velocity, m/s - k Permeability of the porous medium, m² - n Porosity of the porous medium - p( ,t) Pressure, N/m² - Q Volume flux, m³/s - r Radial coordinate, m - r _w Radius of the well, m - s Transform variable - S Storativity of the aquifer - S _d(r, t) Drawdown, m - t Time, s - T Transmissivity of the aquifer, m²/s - ( ,t) Velocity of the fluid, m/s - Coordinate vector, m - z Vertical coordinate, m - Coefficient of compressibility, m²/N - Coefficient of fluid compressibility, m²/N - Relaxation time, s - (r, t) Hydraulic potential, m - Dynamic viscosity of the fluid, Ns/m² - Dimensionless radius - Density of the fluid, Ns²/m⁴ - (, ) Dimensionless drawdown - Dimensionless time - , x Dummy variables - ₀, ₁ Auxilary functions     

Theory of single well tests in a leaky aquifer. Analytical solution for non steady flow

The evaluation of a pump test or a slug test in a single well that completely penetrates a leaky aquifer does not yield a unique relation between the hydraulic properties of the aquifer, independent of the testing conditions. If the flow is transient, the drawdown is characterized by a single similarity parameter that does not distinguish between the storativity and the leakage factor. If the flow is quasi stationary, the drawdown is characterized by a single similarity parameter that does not distinguish between the transmissivity and the leakage factor. The general non steady solution, which is derived in closed form, is characterized bythree similarity parameters.Nomenclature a e 0.8905 = auxiliary parameter - b thickness of the aquifer - b _c thickness of the semipervious stratum - B() auxiliary function - f(s),g(s) auxiliary functions in the complex plane - F(t),G(t) auxiliary functions of time - h undisturbed level of the phreatic surface - K conductivity of the aquifer - K _c conductivity of the semipervious stratum - m ₀ leakage factor - m dimensionless leakage factor - N(s) auxiliary function in the complex plane - Q _w (t) discharge flux - Q steady discharge flux - Q ₀ constant discharge flux during limited time - Q(t) dimensionless discharge flux - r ₀ radius of the well - r radial coordinate - r dimensionless radial coordinate - s complex variable - s ₀ pole - S storativity of the aquifer - S _n n'th part of an integration contour - t time - t dimensionless time - T transmissivity of the aquifer - ,,,,, dimensionless parameters - Euler's number - dummy variable - ₁(), ₂() auxiliary functions - (r, t) drawdown - ₀(t) drawdown in the well - (r, t) dimensionless drawdown - ₀(t) dimensionless drawdown in the well