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.     

Analysis of a method and a system to determine heat transfer and properties of fluids in the absence of temperature measurements and a modified Fourier Equation

The technique to determine by capacitance measurements heat transfer, thermal transport and dielectric properties of fluids introduced recently is now analyzed for a simple system of spherical geometry. The temperature distribution under programmed heat input to a fluid annulus between solid walls is computed by finite difference method for the determination of the capacitance time function of the arrangement. A system of heavy wall structure and heated long enough will produce a capacitance-time curve which is a function of thermal conductivity only. Thermal diffusivity is of influence in thin wall systems. The capacitance change of a heavy wall arrangement is related to the thermal conductivity of the test fluid by a modified Fourier equation. This equation describes the heat flow through the fluid layer but includes the thermal expansion of the solid walls. The change of geometry with T is therefore accounted for. For other multicomposite structures the Fourier equation must be further modified by including the thermal expansion of all materials of the structure and possibly also their compressibilities.

---

Zusammenfassung Die kürzlich eingeführte Methode der Bestimmung von Wärmeübergang, thermischen Transport und dielektrischen Größen mittels Kapazitäts-Zeit-Messung wird analysiert für ein einfaches kugeliges System. Die Temperaturverteilung in der Flüssigkeit im Kugelspalt zwischen zwei festen Körpern wird für konstante Wärmezufuhr von außen mittels der Differenzmethode bestimmt und daraus die Kapazitäts-Zeit-Funktion ermittelt. Es wird gezeigt, daß die Kapazitäts-Zeit-Kurve nur eine Funktion der Wärmeleitzahl ist für den Fall dickwandiger Anordnungen. Für dünnwandige Systeme wird sie auch abhängig von der Temperaturleitzahl. Es wird eine modifizierte Fourier-Gleichung eingeführt, die den Wärmetransport durch die Flüssigkeit beschreibt, dabei aber die Änderung der Geometrie der Schicht berücksichtigt, die sich wegen der thermischen Ausdehnung der festen Wände bei der Einstellung der Temperaturdifferenz ergibt. Für andere mehrschichtige Körper muß die Fourier-Gleichung weiterhin modifiziert werden durch Berücksichtigung der thermischen Ausdehnungskoeffizienten aller beteiligten Materialien und möglicherweise auch ihrer Kompressibilitäten.

---

Nomenclature A average cross-sectional area of fluid layer - A coefficient matrix - B matrix defined by Eq. (20) - B₀ geometric constant of fluid layer (A/L) at reference temperature - C capacitance of arrangement - C_i, C_r capacitance of layer of fluid i and reference fluid at temperature T - capacitances at reference temperature - C_H, c_l specific heats of outer and inner wall - FA...FE constants defined in Eqs. (13 ... 17) - L thickness of fluid layer - M_H, M_L mass of outer and inner wall - P power input to the system - R constant defined by Eq. (24) - T temperature - T_ref reference temperature - T (O, t), T (L, t) temperatures of outer and inner wall at time t - T _i ⁿ , T _i+0 ^n+m temperatures at location i and time n (m=number of t's; 0=number of x's) - T temperature difference across fluid layer - T apparent temperature difference - t_h, T_l temperature increases of outer and inner wall - T_max temperature change of system from one to another thermal equilibrium condition a thermal diffusivity - k, k_i, k_r thermal qonductivity of fluids and of fluid i and reference fluid - q heat flow through fluid layer - ^rh,^rl inner radius of outer wall and outer radius of inner wall - ^rOH,^rOL radii at reference temperature - t time - t time interval - x coordinate - ¯x vector of unknown T_i ⁿ⁺¹ - x length interval Greek symbols linear thermal expansion coefficient - _H, _L linear thermal expansion coefficient of materials of outer and inner wall - dielectric constant - _i, _ref dielectric constant of fluid i and reference fluid - ₀ permittivity of free space - multiplyer of conduction Eq. (7) in finite difference form - time needed to establish quasi-steady state conditions in the system heated by a constant power input In honor of Prof. Dr. E. Schmidt to his 80th Birthday.     

Finite difference analysis of plane Poiseuille and Couette flow developments

Summary The problem of flow development from an initially flat velocity profile in the plane Poiseuille and Couette flow geometry is investigated for a viscous fluid. The basic governing momentum and continuity equations are expressed in finite difference form and solved numerically on a high speed digital computer for a mesh network superimposed on the flow field. Results are obtained for the variations of velocity, pressure and resistance coefficient throughout the development region. A characteristic development length is defined and evaluated for both types of flow.Nomenclature h width of channel - L ratio of development length to channel width - p fluid pressure - p ₀ pressure at channel mouth - P dimensionless pressure, p/ ² - P ₀ dimensionless pressure at channel mouth - P pressure defect, P ₀–P - (P)₀ pressure defect neglecting inertia - Re Reynolds number, uh/ - u fluid velocity in x-direction - mean u velocity across channel - u ₀ wall velocity - U dimensionles u velocity u/ - U _c dimensionless centreline velocity - U ₀ dimensionless wall velocity - v fluid velocity in y-direction - V dimensionless v velocity, hv/ - x coordinate along channel - X dimensionless x-coordinate, x/h ² - y coordinate across channel - Y dimensionless y-coordinate, y/h - resistance coefficient, - ₀ resistance coefficient neglecting inertia - fluid density - fluid viscosity     

A solitary wave in a porous medium

The one-phase Darcy continuity equation, including the quadratic gradient term, is considered. The exact linearization of the equation is found by a functional transformation for an arbitrary spatial dimension in the limit case where the constant fluid compressibility is much more dominant than the constant compressibilities of the reservoir parameters.The equation permits a solution representing a localized wave travelling through a one-dimensional reservoir without changing its form. This is the actual long-time limit of the transient solution for a constant sandface-rate injection of a compressible fluid with a constant compressibility if the fluid is much more compressible than the matrix. A solitary wave solution is not possible for production.A fully developed solitary wave would appear only for very high pressure increases, but the first signs of the emerging solitary wave are detectable at the sandface for moderate pressure increases which can appear under physical reservoir conditions.Latin symbols a Dimensionless wave propagation velocity - A _N Sandface area (N = 0, 1, 2) - c ₁, c ₂ Sums of compressibilities - c _x Generic (generalized) compressibility - c Fluid compressibility - c _h Reservoir height (i.e. bulk volume) compressibility (N = 0, 1) - c _k , c , c Generalized compressibilities - D Spatial reservoir dimensionality (D = 1, 2, 3) - f Fractional change of p _n1 due to nonlinear effects - h Reservoir height (proportional to bulk volume for N = 0, 1) - Horizontal reservoir width (N = 0) - k Reservoir permeability - K _N Constant with dimension of pressure (N = 0, 1, 2) - n Sum index - N Integer variable (N = D – 1) - p Reservoir pressure - p^* Overburden pressure - p _D Dimensionless (scaled) version of p - p ₀ Initial pressure - q Volumetric flow rate referred to sandface - r Radial (or linear) spatial distance from center of well - r _w Well radius - r _e External reservoir radius (or length) from center of well - t Time variable - t _f Injection/production time corresponding to fraction f - T Cole-Hopf-transformed version of dimensionless pressure y - u Rescaled (dimensionless) version of v _D - v Darcy velocity - v _d Dimensionless (scaled) version of v - x Generic symbol in compressibility expression (also used for auxiliary function and for auxiliary variable) - y Rescaled (dimensionless) version of p _D - z Dimensionless (scaled) version of r Greek symbols Coefficient of inertial resistance - Variable in wave solution for y - p _n1 Absolute change in physical sandface pressure due to production or injection - _p Pressure change over (dimensionless) distance behind and far away from front - _r Physical distance at constant time corresponding to - Characteristic (dimensionless) width of solitary wave - Formation porosity - ₁, ₂ Integration constants - Dimensionless (scaled) length of finite reservoir - Fluid viscosity - Fluid density - Dimensionless (scaled) version of t - Wave solution for dimensionless pressure y - Integer variable (±1) distinguishing between production and injection     

New concepts in relative permeabilities at high capillary numbers for surfactant flooding

Flooding oil reservoirs with surfactant solutions can increase the amount of oil that can be recovered. Macroscopic modelling of the process requires relative permeabilities to be functions of saturation and capillary number. With only limited experimental data, relative permeabilities have usually been assumed to be linear functions of saturation at high capillary numbers. The experimental data is reviewed, some of which suggest that this assumption is not necessarily correct. The basis for the assumption is therefore reviewed and it is concluded that the linear model corresponds to microscopically segregated flow in the porous medium. Based on new but equally plausible complementary assumptions about the flow pattern, a mixed flow model is derived. These models are then shown to be limiting cases of a droplet model which represents the mixing scale within the porous medium and gives a physical basis for interpolating between the models. The models are based on physical concepts of flow in a porous medium and so the approach described here represents a significant improvement in the understanding of high capillary number flow. This is shown by the fact that fewer parameters are needed to describe experimental data.Notation A total cross-sectional area assigned to capillary bundle - A ⁽ⁱ⁾ physical cross-sectional area of tube i - c ⁽ⁱ⁾ ordered configurational label for droplets in tube i - c configuration label for tube i (order not considered) - D defined by Equation (26) - E(...) expectation value with respect to the trinomial distribution - S _r ⁽⁾ fractional flow of phase - k absolute permeability - k _r relative permeability of phase - k _r ⁰ endpoint relative permeability of phase - L capillary tube length in bundle model - m ⁽ⁱ⁾ number of droplets of phase a occupying tube i - n exponent for phase a in Equation (2) - N number of droplets in bundle model - N _c capillary number - p pressure - p(c') probability of configuration c - Q ⁽ⁱ⁾ total volume flow rate in tube i - S saturation of phase - S flowing saturation of phase - S _r residual saturation of phase - S _r ⁽⁾ saturations when fractional flow of phase is 1 in the case of varying residual saturations for three-phase flow ( ) - t _c residence time for droplet configuration c - v ⁽ⁱ⁾ total fluid velocity in bundle tube i - , phase label - p pressure differential across capillary bundle - ⁽ⁱ⁾ tube conductivity defined by Equation (7) - viscosity of phase - interfacial tension - gradient operator - ... average over tube droplet configurations     

Mean and turbulent velocity measurements of supersonic mixing layers

The behavior of supersonic mixing layers under three conditions has been examined by schlieren photography and laser Doppler velocimetry. In the schlieren photographs, some large-scale, repetitive patterns were observed within the mixing layer; however, these structures do not appear to dominate the mixing layer character under the present flow conditions. It was found that higher levels of secondary freestream turbulence did not increase the peak turbulence intensity observed within the mixing layer, but slightly increased the growth rate. Higher levels of freestream turbulence also reduced the axial distance required for development of the mean velocity. At higher convective Mach numbers, the mixing layer growth rate was found to be smaller than that of an incompressible mixing layer at the same velocity and freestream density ratio. The increase in convective Mach number also caused a decrease in the turbulence intensity ( _u/U).List of symbols a speed of sound - b total mixing layer thickness between U ₁ – 0.1 U and U ₂ + 0.1 U - f normalized third moment of u-velocity, f u³/(U)³ - g normalized triple product of u² , g u²/(U)³ - h normalized triple product of u ², h u ²/(U)³ - l _u axial distance for similarity in the mean velocity - l _u axial distance for similarity in the turbulence intensity - M Mach number - M _c convective Mach number (for ₁ = ₂), M _c (U ₁ – U ₂)/(a ₁ + a ₂) - P static pressure - r freestream velocity ratio, r U ₂/U ₁ - Re unit Reynolds number, Re U/ - s freestream density ratio, s ₂/₁ - T _t total temperature - u instantaneous streamwise velocity - u deviation of u-velocity, uu – U - U local mean streamwise velocity - U ₁ primary freestream velocity - U ₂ secondary freestream velocity - average of freestream velocities, (U ₁ + U ₂)/2 - U freestream velocity difference, U U ₁ – U ₂ - instantaneous transverse velocity - v deviation of -velocity, – V - V local mean transverse velocity - x streamwise coordinate - y transverse coordinate - y ₀ transverse location of the mixing layer centerline - ensemble average - ratio of specific heats - boundary layer thickness (y-location at 99.5% of free-stream velocity) - similarity coordinate, (y – y ₀)/b - compressible boundary layer momentum thickness - viscosity - density - standard deviation - dimensionless velocity, (U – U ₂)/U - 1 primary stream - 2 secondary stream A version of this paper was presented at the 11th Symposium on Turbulence, October 17–19, 1988, University of Missouri-Rolla     

Transport in ordered and disordered porous media I: The cellular average and the use of weighting functions

In this work we consider transport in ordered and disordered porous media using singlephase flow in rigid porous mediaas an example. We defineorder anddisorder in terms of geometrical integrals that arise naturally in the method of volume averaging, and we show that dependent variables for ordered media must generally be defined in terms of thecellular average. The cellular average can be constructed by means of a weighting function, thus transport processes in both ordered and disordered media can be treated with a single theory based on weighted averages. Part I provides some basic ideas associated with ordered and disordered media, weighted averages, and the theory of distributions. In Part II a generalized averaging procedure is presented and in Part III the closure problem is developed and the theory is compared with experiment. Parts IV and V provide some geometrical results for computer generated porous media.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² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor, m² - L general characteristic length for volume averaged quantities, m - characteristic length (pore scale) for the-phase - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - N unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p0 reference pressure in the-phase, N/m² - p traditional intrinsic volume averaged pressure, N/m² - r0 radius of a spherical averaging volume, m - r position vector, m - r position vector locating points in the-phase, m - averaging volume, m³ - V 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 traditional superficial volume averaged 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 - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

Chaos and integrability in a nonlinear wave equation

We consider the parametrized family of equations _tt ,u- _xx u-au+u ₂ ² u=O,x(0,L), with Dirichlet boundary conditions. This equation has finite-dimensional invariant manifolds of solutions. Studying the reduced equation to a four-dimensional manifold, we prove the existence of transversal homoclinic orbits to periodic solutions and of invariant sets with chaotic dynamics, provided that =2, 3, 4,.... For =1 we prove the existence of infinitely many first integrals pairwise in involution.     

Viscoelastic properties of semidilute poly(methyl methacrylate) solutions

Viscoelastic properties were examined for semidilute solutions of poly(methyl methacrylate) (PMMA) and polystyrene (PS) in chlorinated biphenyl. The number of entanglement per molecule, N, was evaluated from the plateau modulus, G _N . Two time constants, _s and ₁, respectively, characterizing the glass-to-rubber transition and terminal flow regions, were evaluated from the complex modulus and the relaxation modulus. A time constant _k supposedly characterizing the shrink of an extended chain, was evaluated from the relaxation modulus at finite strains. The ratios ₁/ _s and _k / _s were determined solely by N for each polymer species. The ratio ₁/ _s was proportional to N ^4.5 and N ^3.5 for PMMA and PS solutions, respectively. The ratio _k / _s was approximately proportional to N ^2.0 in accord with the prediction of the tube model theory, for either of the polymers. However, the values for PMMA were about four times as large as those for PS. The result is contrary to the expectation from the tube model theory that the viscoelasticity of a polymeric system, with given molecular weight and concentration, is determined if two material constants _s and G _N are known.     

Lower semicontinuity of the attractor for a singularly perturbed hyperbolic equation

For a smooth, bounded domain R, n 3, and a real, positive parameter, we consider the hyperbolic equationu _tt +u _t –u=–f(u) –g in with Dirichlet boundary conditions. Under certain conditions onf, this equation has a global attractorA inH ₀ ¹ () ×L ²(). For=0, the parabolic equation also has a global attractor which can be naturally embedded into a compact setA ₀ inH ₀ ¹ () ×L ²(). If all of the equilibrium points of the parabolic equation are hyperbolic, it is shown that the setsA are lower semicontinuous at=0. Moreover, we give an estimate of the symmetric distance betweenA ₀ andA .     

A visualization study of the interaction of a free vortex with the wake behind an airfoil

The two-dimensional interaction of a single vortex with a thin symmetrical airfoil and its vortex wake has been investigated in a low turbulence wind tunnel having velocity of about 2 m/s in the measuring section. The flow Reynolds number based on the airfoil chord length was 4.5 × 10³. The investigation was carried out using a smoke-wire visualization technique with some support of standard hot-wire measurements. The experiment has proved that under certain conditions the vortex-airfoil-wake interaction leads to the formation of new vortices from the part of the wake positioned closely to the vortex. After the formation, the vortices rotate in the direction opposite to that of the incident vortex.List of symbols c test airfoil chord - C vortex generator airfoil chord - TA test airfoil - TE test airfoil trailing edge - TE _G vortex generator airfoil trailing edge - t dimensionless time-interval measured from the vortex passage by the test airfoil trailing edge: gDt=(T-T- _TE)·U/c - T time-interval measured from the start of VGA rotation - U free stream velocity - U vortex induced velocity fluctuation - VGA vortex generator airfoil - y distance in which the vortex passes the test airfoil - Z vortex circulation coefficient: Z=/(U · c/2) - vortex generator airfoil inclination angle - vortex circulation - vortex strength: =/2     

On slip effect in free coating of non-Newtonian fluids

The failure of the current theories to predict the coating thickness of non-Newtonian fluids in free coating operations is shown to be a result of the effective slip at the moving rigid surface being coated. This slip phenomenon is a consequence of stress induced diffusion occurring in flow of structured liquids in non-homogeneous flow fields. Literature data have been analysed to substantiate the slip hypothesis proposed in this work. The experimentally observed coating thickness is shown to lie between an upper bound, which is estimated by a no-slip condition for homogeneous solution and a lower bound, which is estimated by using solvent properties. Some design considerations have been provided, which will serve as useful guidelines for estimating coating thickness in industrial practice.fa exponent in eq. (15) - b n/(4 –n)(n + 1) - Ca Capillary number - D diffusivity - De Deborah number - g acceleration due to gravity - G Goucher number - h thickness profile - h ₀ final coating thickness - K consistency index - L length available for diffusion - L _t tube length - n power-law index - P pressure drop - Q flow rate - R cylinder radius - R _t tube radius - t time available for diffusion - T ₀ dimensionless thickness without slip - T _s dimensionless thickness with slip - U _c theoretically calculated withdrawal velocity to match the film thickness - u _s slip velocity - U withdrawal velocity - U _w theoretically calculated withdrawal velocity based on solvent properties - U ^* effective withdrawal velocity - x distance in the direction of flow - y distance transverse to the flow direction - curvature coefficient - slip coefficient - curvature coefficient - rate of deformation tensor - u _s /U - relaxation time - density - surface tension - shear stress in tube flow - _w wall shear stress in tube flow - stress tensor - _w wall shear stress - T _s /T ₀ NCL-Communication No. 2818     

Zur Berechnung der Filmverdunstung von Kohlenwasserstoffen in einem Heißluftstrom

The results of an analytical approximation method to predict the film vaporization are compared with the predictions of a finite difference method of Hermitian type. The analytically estimated rate of vaporization of different hydrocarbons, which is the most important value for practical applications, deviates only a few percents from the numerically estimated value.

---

Zur Berechnung der Filmverdunstung von Kohlenwasserstoffen in einem Heißluftstrom

Zusammenfassung Es wird ein Näherungsverfahren zur Berechnung der Filmverdunstung dargestellt, bei dem eine vollständige Lösung der miteinander gekoppelten Grenzschichtgleichungen entfallt. Die nach dieser analytischen Methode ermittelte Verdunstung verschiedener Kohlenwasserstoffe wird mit Werten verglichen, die nach einem Differenzenverfahren vom Hermiteschen Typ berechnet wurden. Es zeigt sich, daß die analytisch berechnete Verdunstungsrate, die für praktische Anwendungen wichtigste Größe, nur wenige Prozent von dem numerisch ermittelten Wert abweicht.

---

Formelzeichen c _ew Konzentrationsdifferenz c_1e -c _1w - c _i Massenkonzentration der Komponentei - c_p, c_pi spezifische Wärmekapazität bei konstantem Druck des Gemisches — der Komponentei - D ₁₂ binärer Diffusionskoeffizient - f dimensionslose Stromfunktion - f dimensionslose Geschwindigkeit - g () allgemeine Funktion - m ₁ Massenstromdichte der Komponente 1 - m ^* dimensionslose Massenstromdichte, G1. (4.8) - M, M_i Molgewicht, — der Komponentei - P, P _i Druck, Partialdruck der Komponentei - Pr Prandtlzahl,C _p/ - q Wärmestromdichte - r ₁ Verdampfungswärme - R allgemeine Gaskonstante - Sc Schmidtzahl/D ₁₂ - T absolute Temperatur - u Geschwindigkeitskomponente inx-Richtung - v Geschwindigkeitskomponente iny-Richtung - x Längskoordinate - y Querkoordinate - z dimensionslose Konzentration - dimensionslose Funktion/ _{e e} - transformierte Koordinatey - dimensionslose Temperatur (T-T _w)/(T_e-T_w) - Wärmeleitfähigkeit des Gemisches - Zähigkeit des Gemisches - transformierte Koordinate - Dichte des Gemisches - Stromfunktion Indizes e am Außenrand der Grenzschicht - i Stoffi - w an der Filmoberfläche - 1, 2 Komponente 1, 2 - () Ableitung ()/ _n      

Measurements of thermally stratified pipe flow using image-processing techniques

The cross-correlation technique and Laser Induced Fluorescence (LIF) have been adopted to measure the time-dependent and two-dimensional velocity and temperature fields of a stably thermal-stratified pipe flow. One thousand instantaneous and simultaneous velocity and temperature maps were obtained at overall Richardson numberRi = 0 and 2.5, from which two-dimensional vorticity, Reynolds stress and turbulent heat flux vector were evaluated. The quasi-periodic inclined vortices (which connected to the crest) were revealed from successive instantaneous maps and temporal variation of vorticity and temperature. It has been recognized that these vortices are associated with the crest and valley in the roll-up motion.List of symbols A Fraction of the available light collected - C Concentration of fluorescence - D Pipe diameter - I Fluorescence intensity - L Sampling length along the incident beam - I ₀ Intensity of an excitation beam - I _c (T) Calibration curve between temperature and fluorescence intensity - I _ref Reference intensity of fluorescence radiation - Re _b Reynolds number based on bulk velocity,U _b D/v - Ri Overall Richardson number based on velocity difference,gDT/U ² - t Time - t Time interval between the reference and corresponding matrix - T Temperature - T ₁,T ₂ Temperature of lower and upper layer - T ^* Normalized temperature, (T–T ₁)/T - T _c (I) Inverse function of temperature as a function ofI _c - T _ref Reference temperature - T Temperature difference between upper and lower flow,T ₂ –T ₁ - U ₁ Velocity of lower stream - U ₂ Velocity of upper stream - U _b Bulk velocity - U _c Streamwise mean velocity atY/D=0 - U Streamwise velocity difference between upper and lower flow,U ₁–U ₂ - u, v, T Fluctuating component ofU, V, T - U, V Velocity component of X, Y direction - X Streamwise distance from the splitter plate - Y Transverse distance from the centerline of the pipe - Z Spanwise distance from the centerline of the pipe - Quantum yield - Absorptivity - vorticity calculated from a circulation - Kinematic viscosity - circulation     

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     

On the effects of uniform high suction on the rotationally symmetric flow of a conducting liquid near a stationary disc in the presence of a transverse magnetic field

In the present paper an attempt has been made to find out effects of uniform high suction in the presence of a transverse magnetic field, on the motion near a stationary plate when the fluid at a large distance above it rotates with a constant angular velocity. Series solutions for velocity components, displacement thickness and momentum thickness are obtained in the descending powers of the suction parameter a. The solutions obtained are valid for small values of the non-dimensional magnetic parameter m (= 4 _e ² H ₀ ² /) and large values of a (a2).Nomenclature a suction parameter - E electric field - E _r , E , E _z radial, azimuthal and axial components of electric field - F, G, H reduced radial, azimuthal and axial velocity components - H magnetic field - H _r , H , H _z radial, azimuthal and axial components of magnetic field - H ₀ uniform magnetic field - H* displacement thickness and momentum thickness ratio, */ - h induced magnetic field - h _r , h , h _z radial, azimuthal and axial components of induced magnetic field - J current density - m nondimensional magnetic parameter - p pressure - P reduced pressure - R Reynolds number - U ₀ representative velocity - V velocity - V _r , V , V _z radial, azimuthal and axial velocity components - w ₀ uniform suction through the disc. - density - electrical conductivity - kinematic viscosity - _e magnetic permeability - a parameter, (/)^1/2 z - a parameter, a - * displacement thickness - momentum thickness - angular velocity     

Heat transfer between gas-liquid spray stream flowing perpendicularly to the row of the cylinders

Experimental investigation and analysis of heat transfer process between a gas-liquid spray flow and the row of smooth cylinders placed in the surface perpendicular to the flow has been performed. Among others, there was taken into account in the analysis the phenomenon of droplets bouncing and omitting the cylinder as well as the phenomenon of the evaporation process from the liquid film surface.In the experiments test cylinders were used, which were placed between two other cylinders standing in the row.From the experiments and the analysis the conclusion can be drawn that the heat transfer coefficients values for a row of the cylinders are higher than for a single cylinder placed in the gasliquid spray flow.

---

Wärmeübergang an eine senkrecht anf eine Zylinderreihe auftreffende Gas-Flüssigkeits-Sprüh-Strömung

Zusammenfassung Es wurden Messungen und theoretische Analysen des Wärmeübergangs zwischen einer Gas-FlüssigkeitsSprüh-Strömung und den glatten Oberflächen einer Zylinderreihe durchgeführt, die senkrecht zum Sprühstrahl angeordnet waren. Dabei wurde in der Analyse unter anderem das Phänomen betrachtet, daß die Tropfen die Zylinderwand treffen und verfehlen können und daß sich ein Verdampfungsprozeß aus dem flüssigen Film an der Zylinderoberfläche einstellt.Gemessen wurde an einem zwischen zwei Randzylindern befindlichen Zylinder.Die Experimente und die Analyse gestatten die Schlußfolgerung, daß der Wärmeübergangskoeffizient für eine Zylinderreihe höher ist als für einen einzelnen Zylinder in der Sprühströmung.

---

Nomenclature a distance between axes of cylinders, m - c _l specific heat capacity of liquid, J/kg K - c _g specific heat capacity of gas, J/kg K - D cylinder diameter, m - g _l mass velocity of liquid, kg/m²s - ¯k average volume ratio of liquid entering film to amount of liquid directed at the cylinder in gas-liquid spray flow, dimensionless - k() local volume ratio of liquid entering film to amount of liquid directed at the cylinder in gas-liquid spray flow, dimensionless - L specific latent heat of vaporisation, J/kg - m mass fraction of water in gas-liquid spray flow, dimensionless - M constant in Eq. (9) - p pressure, Pa - p _g statical pressure of gas, Pa - p _w pressure of gas on the cylinder surface, Pa - p external pressure on the liquid film surface, Pa - r cylindrical coordinate, m - R radius of cylinder, m - T temperature, K, °C - T _l, t_l liquid temperature in the gas-liquid spray, K, °C - T _w,t_w temperature of cylinder surface, K, °C - T temperature of gas-liquid film interface, K - U liquid film velocity, m/s - w gas velocity on cylinder surface, m/s - w _g gas velocity in free stream, m/s - W _l liquid vapour mass ratio in free stream, dimensionless - W liquid vapour mass ratio at the edge of a liquid film, dimensionless - x coordinate, m - y coordinate, m - z complex variable, dimensionless - average heat transfer coefficient, W/m²K - local heat transfer coefficient, W/m² K - average heat transfer coefficient between cylinder surface and gas, W/m² K - _g, local heat transfer coefficient between cylinder surface and gas, W/m² K - mass transfer coefficient, kg/m²s - liquid film thickness, m - _lg dynamic diffusion coefficient of liquid vapour in gas, kg/m s - pressure distribution function on a cylinder surface - function defined by Eq. (3) - _l liquid dynamic viscosity, kg/m s - _g gas dynamic viscosity, kg/m s - cylindrical coordinate, rad, deg - _l thermal conductivity of liquid, W/m K - _g thermal conductivity of gas, W/m K - mass transfer driving force, dimensionless - _l density of liquid, kg/m³ - _g density of gas, kg/m³ - _w shear stress on the cylinder surface, N/m² - _w shear stress exerted by gas at the liquid film surface, N/m² - air relative humidity, dimensionless - T -T _w - _w =T _wT_l Dimensionless parameters I= enhancement factor of heat transfer - m ^*=M _l/M_g molar mass of liquid to the molar mass of gas ratio - Nu _g= D/ _g gas Nusselt number - Pr _g=c _{g g}/_g gas Prandtl number - Pr _l=c_lll liquid Prandtl number - ¯r=(r–R)/ dimensionless coordinate - Re _g=w_gD _g/_g gas Reynolds number - Re _g,max=w_g,max D _g/_g gas Reynolds number calculated for the maximal gas velocity between the cylinders - Sc=m ^* _g/_l–g Schmidt number =/R dimensionless film thickness     

Transport in ordered and disordered porous media II: Generalized volume averaging

In this paper we develop the averaged form of the Stokes equations in terms of weighting functions. The analysis clearly indicates at what point one must choose a media-specific weighting function in order to achieve spatially smoothed transport equations. The form of the weighting function that produces the cellular average is derived, and some important geometrical theorems are presented.Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - A _p surface area of a particle, m² - d _p 6V _p/A_p, effective particle diameter, m - g gravity vector, m/s² - I unit tensor - 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 general characteristic length for volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - l characteristic length (pore scale) for the-phase - l _i i=1, 2, 3 lattice vectors, m - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - 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 - m _D special convolution product weighting function for disordered media - m _M master convolution product weighting function for ordered and disordered media - 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 p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function, m - r position vector, m - r position vector locating points in the-phase, m - V averaging volume, m³ - V 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 - vm superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - V volume of the-phase contained in the averaging volume, m³ - V _p volume of a particle, m³ - v traditional superficial volume averaged velocity, m/s - v v p _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² - V /V, volume fraction of the-phase     

Heat and mass transfer in laminar flow between parallel porous plates

Summary The problem of heat transfer in a two-dimensional porous channel has been discussed by Terrill [6] for small suction at the walls. In [6] the heat transfer problem of a discontinuous change in wall temperature was solved. In the present paper the solution of Terrill for small suction at the walls is revised and the whole problem is extended to the cases of large suction and large injection at the walls. It is found that, for all values of the Reynolds number R, the limiting Nusselt number Nu increases with increasing R.Nomenclature stream function - 2h channel width - x, y distances measured parallel and perpendicular to the channel walls respectively - U velocity of fluid at x=0 - V constant velocity of fluid at the wall - =y/h nondimensional distance perpendicular to the channel walls - f() function defined in equation (1) - coefficient of kinematic viscosity - R=Vh/ suction Reynolds number - density of fluid - C _p specific heat at constant pressure - K thermal conductivity - T temperature - x=x ₀ position where temperature of walls changes - T ₀, T ₁ temperature of walls for x<x ₀, x>x ₀ respectively - = (T – T ₁)/T ₀ – T ₁) nondimensional temperature - =x/h nondimensional distance along channel - R ^* = Uh/v channel Reynolds number - Pr = C _p/K Prandtl number - _n eigenvalues - B _n() eigenfunctions - B _n ⁽⁰⁾ , () eigenfunctions for R=0 - B ₀ ⁽ⁱ⁾ , B ₀ ⁽ⁱⁱ⁾ ... change in eigenfunctions when R0 and small - K _n constants given by equation (13) - h heat transfer coefficient - Nu Nusselt number - _m mean temperature - C _n constants given by equation (18) - perturbation parameter - B _0i () perturbation approximations to B ₀() - Q = B ₀/ ₀ derivative of eigenfunction with respect to eigenvalue - z nondimensional distance perpendicular to the channel walls - F(z) function defined by (54)     

A class of hypo-elastic non-elastic materials and their thermodynamics

Thermodynamics is developed for a class of thermo-hypo-elastic materials. It is shown that materials of this class obey the laws of thermodynamics, but are not elastic.

Table of Symbols

Latin Letters A _ijkl tensor-valued function of t _ij appearing in hypo-elastic constitutive relation - B _ijkl another tensor-valued function. See equation (4.2) - B the square of - d _ij rate of deformation tensor - d _ij deviator of rate of deformation - f, k functions of pressure, p - g, h functions of the invariant - p pressure - q _i heat flux vector - s _ij stress deviator - _ij co-rotational derivative of stress deviator - t time - t ₁ t ₂ specific values of time - t _ij stress tensor - t _ij ⁰ a specific value of stress - T Temperature - T ₀ a specific value of temperature - u _i velocity - V(t) a material volume as a function of time, t - V ₀ a material volume at a reference configuration - W work (W = work done in a deformation—section 5) Sript Letters Specific internal energy - Specific Helmholtz free energy - G Specific Gibbs function Greek Letters an invariant of the stress deviator—see eq. (2.4) - _ij kroneker delta - (W = work done in a deformation—section 5) - specific entropy - hypo-elastic potential - hypo-elastic potential - mass density - ₀ mass density in a reference configuration - specific volume = 1/ - a function of p - _ijkl a constant tensor—see eq. (2.5) - –G/ - _ij rate of rotation tensor This work is dedicated to Jerald L. Ericksen, without whose influence it would not have been possible