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.     

Asymptotic theory for calculating heat fluxes near the corner of a body

A method is presented for calculating the distribution of the thermal fluxes, friction stresses, and pressure near the corner point of a body contour in whose vicinity the outer supersonic flow passes through an expansion wave. The method is based on a study of the asymptotic solutions of the Navier-Stokes equations as the Reynolds number R approaches infinity for the flow region in which the longitudinal gradients of the flow functions are large, invalidating conventional boundary layer theory. This problem was examined in part in [1], in which the distribution of the friction and pressure in a region with length on the order of a few thicknesses of the approaching boundary layer was obtained in the first approximation. The leading term of the expansion for the thermal flux to the surface of the body vanishes for a value of the Prandtl number equal to unity and for other values of the Prandtl number does not match directly with its value in the undisturbed boundary layer.The thermal-flux distribution is obtained for values of the Prandtl number approaching unity. For this purpose it was necessary to consider a more general double passage to the limit as 1 and 0 for a finite value of the parameter B=[(–1)/] [–ln ^1/4/]^1/4 characterizing the ratio of the effects of thermal conduction, viscous dissipation, and convection. The solution obtained previously [1] corresponds to the particular case B and therefore for actual values of R=10⁴–10⁶, ~ 0.7 overestimates considerably the effect of the dissipative term on heat transfer, although even in first approximation it describes the pressure distribution well and the friction distribution satisfactorily. For smooth matching of the solutions with the corresponding flow functions in the undisturbed boundary layer it was necessary to introduce a flow region with free interaction for the expansion flow. Equations and boundary conditions which describe the flow as a whole are presented. Examples are given of numerical calculations and comparison with experiment.     

Flow in porous media III: Deformable media

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

Stress analysis of the cylinder subjected to arbitrarily distributed axisymmetrical loads

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

---

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

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

     

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     

Mixed convection from a horizontal plate to fluids of any Prandtl number

Laminar mixed convection over a horizontal plate with uniform wall temperature or uniform wall heat flux is analyzed by introducing proper buoyancy parameters and transformation variables for fluids of any Prandtl number between 0.001 and 10,000. Both cases of buoyancy assisting and opposing flow conditions are investigated. For the buoyancy-assisting case, the obtained numerical results are very accurate over the entire range of mixed convection intensity from pure forced convection limit to pure free convection limit. For the buoyancy-opposing case, solutions are obtained from the forced convection limit to the point of breakdown.

---

Mischkonvektion an einer horizontalen Platte für Fluide mit beliebiger Prandtl-Zahl

Zusammenfassung Es wurde laminare Mischkonvektion an einer horizontalen Platte mit einheitlicher Wandtemperatur oder einheitlicher Wandwärmestromdichte bei Einführung zweckmäßiger Auftriebsparameter und Transformationsvariablen für Fluide mit beliebiger Prandtl-Zahl zwischen 0,001 und 10 000 untersucht. Es wurden die Fälle der Strömung entgegen und in Richtung der Auftriebskraft untersucht. Für den Fall der Strömung in Richtung der Auftriebskraft wurden sehr genaue numerische Ergebnisse für den gesamten Bereich der gemischten Konvektion von rein erzwungener Konvektion bis zu rein freier Konvektion erhalten. Für den Fall der Strömung entgegen der Auftriebsrichtung wurden Lösungen für erzwungene Konvektion bis zum Umkehrpunkt erhalten.

---

Nomenclature C _f local friction coefficient - f reduced stream function - g gravitational acceleration - Gr local Grashof number for UWT,g (T _w –T )x ³/ ² - Gr* local Grashof number for UHF,g q _w x ⁴/k ² - m =10 for UWT; and =6 for UHF - n =5 for UWT; and =3 for UHF - Nu local Nusselt number - p pressure - Pr Prandtl number,/ - q _w wall heat flux - Ra local Rayleigh number for UWT,Gr Pr - Ra* local Rayleigh number for UHF,Gr*Pr - Re local Reynolds number,u x/ - T fluid temperature - T _w wall temperature - T free-stream temperature - u velocity component inx-direction - u free-stream velocity - v velocity component iny-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - =0 for UWT; and =1 for UHF - buoyancy parameter, =( Ra)^1/5/( Re)^1/2 for UWT; and =( Ra*)^1/6/( Re)^1/2 for UHF - pseudo-similarity variable, (y/x) - dimensionless temperature, =(T–T )/(T _w –T ) for UWT; and =(T–T )/(q _w x/k) for UHF - =[( Re)^1/2+( Ra)^1/5] for UWT; and =[( Re)^1/2+( Ra*)^1/6] for UHF - dynamic viscosity - kinematic viscosity - /(1+) - dimensionless pressure - density - Pr/(1+Pr) - _w wall shear stress,(u/y) _y=0 - stream function - Pr/(1+Pr)^1/3     

Scalar Minimizers with Fractal Singular Sets

Lack of regularity of local minimizers for convex functionals with non-standard growth conditions is considered. It is shown that for every >0 there exists a function aC() such that the functional admits a local minimizer uW^1,p() whose set of non-Lebesgue points is a closed set with dim()>N–p–, and where 1<p<N<N+<q<+.     

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

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

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²     

Laser velocimetry measurements in a horizontal gas-solid pipe flow

This paper presents laser measurements of particle velocities in a horizontal turbulent two-phase pipe flow. A phase Doppler particle analyzer, (PDPA), was used to obtain particle size, velocity, and rms values of velocity fluctuations. The particulate phase consisted of glass spheres 50 m in diameter with the volume fraction of the suspension in the range _p=10^-4 to _p=10^-3. The results show that the turbulence increases with particle loading.List of symbols a particle diameter - C _va velocity diameter cross-correlation - d pipe diameter - Fr ² Froude number - g gravitational constant - p(a) Probability density of the particle diameter - Re pipe Reynolds number based on the friction velocity - T characteristic time scale of the energy containing eddies - T _L integral scale of the turbulence sampled along the particle path - u, U, u characteristic fluid velocities: fluctuating, mean and friction - v characteristic velocity of the paricle fluctuations - f expected value of any random variable f - f¦g expected value of f given a value of the random variable g - _p particle volume fraction - _p particle response time - absolute fluid viscosity - v kinematic fluid viscosity - _p, _f densities, particle and fluid - _a ² particle diameter variance - _va ² velocity variance due to the particle diameter variance - _vT ² total particle velocity variance - _vt ² particle velocity variance due to the response to the turbulent field     

Symmetry groups of attractors

For maps equivariant under the action of a finite group on ⁿ, the possible symmetries of fixed points are known and correspond to the isotropy subgroups. This paper investigates the possible symmetries of arbitrary, possibly chaotic, attractors and finds that the necessary conditions of Melbourne, Dellnitz & Golubitsky [15] are sufficient, at least for continuous maps.This result shows that the reflection hyperplanes are important in determining those groups which are admissible; more precisely, a subgroup of is admissible as the symmetry group of an attractor if there exists a with / cyclic such that fixes a connected component of the complement of the set of reflection hyperplanes of reflections in but not in . For finite reflection groups this condition on reduces to the condition that is an isotropy subgroup. Our results are illustrated for finite subgroups of O(3).     

Flow visualization of compressible vortex structures using density gradient techniques

Mathematical results are derived for the schlieren and shadowgraph contrast variation due to the refraction of light rays passing through two-dimensional compressible vortices with viscous cores. Both standard and small-disturbance solutions are obtained. It is shown that schlieren and shadowgraph produce substantially different contrast profiles. Further, the shadowgraph contrast variation is shown to be very sensitive to the vortex velocity profile and is also dependent on the location of the peak peripheral velocity (viscous core radius). The computed results are compared to actual contrast measurements made for rotor tip vortices using the shadowgraph flow visualization technique. The work helps to clarify the relationships between the observed contrast and the structure of vortical structures in density gradient based flow visualization experiments.Nomenclature a Unobstructed height of schlieren light source in cutoff plane, m - c Blade chord, m - f Focal length of schlieren focusing mirror, m - C _T Rotor thrust coefficient, T/( ² R ⁴) - I Image screen illumination, Lm/m ² - l Distance from vortex to shadowgraph screen, m - n _b Number of blades - p Pressure,N/m ² - p Ambient pressure, N/m ² - r, , z Cylindrical coordinate system - r _c Vortex core radius, m - Non-dimensional radial coordinate, (r/r _c ) - R Rotor radius, m - Tangential velocity, m/s - Specific heat ratio of air - Circulation (strength of vortex), m ²/s - Non-dimensional quantity, ² 8²p r _c ² - Refractive index of fluid medium - ₀ Refractive index of fluid medium at reference conditions - Gladstone-Dale constant, m ³/kg - Density, kg/m ³ - Density at ambient conditions, kg/m ³ - Non-dimensional density, (/ ) - Rotor solidity, (n _b c/ R) - Rotor rotational frequency, rad/s     

On the properties of compressible gas flow in a porous media

The flow of an adiabatic gas through a porous media is treated analytically for steady one- and two-dimensional flows. The effect on a compressible Darcy flow by inertia and Forchheimer terms is studied. Finally, wave solutions are found which exhibit a cut-off frequency and a phase shift between pressure and velocity of the gas, with the velocity lagging behind the pressure.Nomenclature A area of tube for one-dimensional flow - B drag coefficient associated with Forchheimer term - c speed of sound - M Mach number - p ^* gas pressure - p dimensionless gas pressure - s coordinate along the axis of tube - t ^* time variable - t dimensionless time variable - V^* gas velocity in the porous media - V dimensionless gas velocity Greek Letters ratio of specific heat capacities - phase angle between gas pressure and velocity for linear waves - parameter indicating the importance of the inertia term - viscosity - _p natural frequency of the porous media - ^* gas density - dimensionless gas density - parameter indicating the importance of the Forchheimer term - porosity of porous media - velocity potential - stream function     

Free convection heat transfer from an isothermal plate with arbitrary inclination

This paper presents a new formulation for the laminar free convection from an arbitrarily inclined isothermal plate to fluids of any Prandtl number between 0.001 and infinity. A novel inclination parameter is proposed such that all cases of the horizontal, inclined and vertical plates can be described by a single set of transformed equations. Moreover, the self-similar equations for the limiting cases of the horizontal and vertical plates are recovered from the transformed equations by setting=0 and=1, respectively. Heated upward-facing plates with positive and negative inclination angles are investigated. A very accurate correlation equation of the local Nusselt number is developed for arbitrary inclination angle and for 0.001 Pr .

---

Wärmeübertragung bei freier Konvektion an einer isothermen Platte mit beliebiger Neigung

Zusammenfasssung Diese Untersuchung stellt eine neue Formulierung der laminaren freien Konvektion von Flüssigkeiten mit einer Prandtl-Zahl zwischen 0,001 und unendlich an einer beliebig schräggestellten isothermen Platte dar. Ein neuer Neigungsparameter wird eingeführt, so daß alle Fälle der horizontalen, geneigten oder vertikalen Platte von einem einzigen Satz transformierter Gleichungen beschrieben werden können. Die unabhängigen Gleichungen für die beiden Fälle der horizontalen and vertikalen Platte wurden für=0 und=1 aus den transformierten Gleichungen wieder abgeleitet. Es wurden erwärmte aufwärtsgerichtete Platten mit positiven und negativen Neigungswinkeln untersucht. Eine sehr genaue Gleichung wurde für die lokale Nusselt-Zahl bei beliebigen Neigungswinkeln und für 0,001 Pr entwickelt.

---

Nomenclature C _p specific heat - f reduced stream function - g gravitational acceleration - Gr local Grashof number,g(T _w –T _w ) x³/v² - h local heat transfer coefficient - k thermal conductivity - n constant exponent - Nu local Nusselt number,hx/k - p pressure - Pr Prandtl number, v/ - Ra local Rayleigh number,g(T _w –T )J x³/v - T fluid temperature - T _w wall temperature - T temperature of ambient fluid - u velocity component in x-direction - v velocity component in y-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - (Ra¦sin¦)^1/4/( Ra cos()^1/5 - pseudo-similarity variable, (y/) - dimensionless temperature, (T–T )/(T _w–T ) - ( Ra cos)^1/5+(Rasin)^1/4 - v kinematic viscosity - 1/[1 +(Ra cos)^1/5/( Ra¦sin)^1/4] - density of fluid - Pr/(1+Pr) - _w wall shear stress - angle of plate inclination measured from the horizontal - stream function - dimensionless dynamic pressure     

Laminar source flow between parallel porous disks with equal suction and injection

An analysis is presented for laminar source flow between parallel stationary porous disks with suction at one of the disks and equal injection at the other. The solution is in the form of an infinite series expansion about the solution at infinite radius, and is valid for all suction and injection rates. Expressions for the velocity, pressure, and shear stress are presented and the effect of the cross flow is discussed.Nomenclature a distance between disks - A, B, ..., J functions of R _w only - F static pressure - p dimensionless static pressure, p(a ²/ ²) - Q volumetric flow rate of the source - r radial coordinate - r dimensionless radial coordinate, r/a - R radial coordinate of a point in the flow region - R dimensionless radial coordinate of a point in the flow region, R - Re source Reynolds number, Q/2a - R _w wall Reynolds number, Va/ - reduced Reynolds number, Re/r ² - critical Reynolds number - velocity component in radial direction - u dimensionless velocity component in radial direction, a/ - average radial velocity, Q/2a - u dimensionless average radial velocity, Re/r - ratio of radial velocity to average radial velocity, u/u - velocity component in axial direction - v dimensionless velocity component in axial direction, v - V magnitude of suction or injection velocity - z axial coordinate - z dimensionless axial coordinate, z a - viscosity - density - kinematic viscosity, / - shear stress at lower disk - shear stress at upper disk - ₀ dimensionless shear stress at lower disk, - ₁ dimensionless shear stress at upper disk, - dimensionless stream function     

On the closure problem for Darcy's law

In a previous derivation of Darcy's law, the closure problem was presented in terms of an integro-differential equation for a second-order tensor. In this paper, we show that the closure problem can be transformed to a set of Stokes-like equations and we compare solutions of these equations with experimental data. The computational advantages of the transformed closure problem are considerable.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the averaging volume, m² - B second-order tensor used to respresent the velocity deviation - b vector used to represent the pressure deviation, m^–1 - C second-order tensor related to the permeability tensor, m^–2 - D second-order tensor used to represent the velocity deviation, m² - d vector used to represent the pressure deviation, m - g gravity vector, m/s² - I unit tensor - K C ^–1,–D, Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - l _i i=1, 2, 3, lattice vectors, m - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of the-phase - p pressure in the-phase, N/m ² - p intrinsic phase average pressure, N/m² - p –p , spatial deviation of the pressure in the-phase, N/m² - r position vector locating points in the-phase, m - r ₀ radius of the averaging volume, m - t time, s - v velocity vector in the-phase, m/s - v intrinsic phase average velocity in the-phase, m/s - v phase average or Darcy velocity in the \-phase, m/s - v –v , spatial deviation of the velocity in the-phase m/s - V averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ Greek Letters V /V volume fraction of the-phase - mass density of the-phase, kg/m³ - viscosity of the-phase, Nt/m²     

Transport in ordered and disordered porous media IV: Computer generated porous media for three-dimensional systems

In the method of volume averaging, the difference between ordered and disordered porous media appears at two distinct points in the analysis, i.e. in the process of spatial smoothing and in the closure problem. In theclosure problem, the use of spatially periodic boundary conditions isconsistent with ordered porous media and the fields under consideration when the length-scale constraint,r ₀L is satisfied. For disordered porous media, spatially periodic boundary conditions are an approximation in need of further study.In theprocess of spatial smoothing, average quantities must be removed from area and volume integrals in order to extractlocal transport equations fromnonlocal equations. This leads to a series of geometrical integrals that need to be evaluated. In Part II we indicated that these integrals were constants for ordered porous media provided that the weighting function used in the averaging process contained thecellular average. We also indicated that these integrals were constrained by certain order of magnitude estimates for disordered porous media. In this paper we verify these characteristics of the geometrical integrals, and we examine their values for pseudo-periodic and uniformly random systems through the use of computer generated porous media.

Nomenclature

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 _i i=1, 2, 3 gaussian probability distribution used to locate the position of particles - I unit tensor - L general characteristic length for volume averaged quantities, m - L characteristic length for , m - L characteristic length for , m - characteristic length for the -phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1, 2, 3 lattice vectors, m - m convolution product weighting function - m _v special convolution product weighting function associated with the traditional volume average - n _i i=1, 2, 3 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - r position vector, m - r _m support of the weighting functionm, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume,, m³ - x positional vector locating the centroid of an averaging volume, m - x ₀ reference position vector associated with the centroid of an averaging volume, m - y position vector locating points 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 - /L, small parameter in the method of spatial homogenization - standard deviation ofa _i - _r standard deviation ofr - _r intrinsic phase average of      

Theoretical derivation of tangential velocity profiles in a flat vortex chamber-influence of turbulence and wall friction

Summary As part of a study on the hydrodynamics of a cyclone separator, a theoretical investigation of the flow pattern in a flat box cyclone (vortex chamber) has been carried out. Expressions have been derived for the tangential velocity profile as influenced by internal friction (eddy viscosity) and wall friction. The most important parameter controlling the tangential velocity profile is = –u ₀ R/(v+ ), where u ₀ is the radial velocity at the outer radius R of the cyclone, the kinematic liquid viscosity and is the kinematic eddy viscosity. For values of greater than about 10 the tangential velocity profile is nearly hyperbolic, for smaller than 1 the tangential velocity even decreases towards the centre. It is shown how and also the wall friction coefficient may be obtained from experimental velocity profiles with the aid of suitable graphs. Because of the close relation between eddy viscosity and eddy diffusion, measurements of velocity profiles in flat box cyclones will also provide information on the eddy motion of particles in a cyclone, a motion reducing its separation efficiency.List of symbols A cross-sectional area of cyclone inlet - h height of cyclone - p static pressure in cyclone - p static pressure difference in cyclone between two points on different radius - r radius in cyclone - r ₁ radius of cyclone outlet - R radius of cyclone circumference - u radial velocity in cyclone - u ₀ radial velocity at circumference of flat box cyclone - v tangential velocity - v ₀ tangential velocity at circumference of flat box cyclone - w axial velocity - z axial co-ordinate in cyclone - friction coefficient in flat box cyclone (for definition see § 5) - ₁ value of friction coefficient for ₁<< ₂ - ₂ value of friction coefficient for ₂<<1 - = - ₁ value of for ₁<< ₂ - ₂ value of for ₂<<1 - thickness of laminar boundary layer - =/h - turbulent kinematic viscosity - ratio of z to h - k ratio of height of cyclone to radius R of cyclone - parameter describing velocity profile in cyclone =–u ₀ R/(+) - kinematic viscosity of fluid - density of fluid - ratio of r to R - ₁ value of at outlet of cyclone - ₂ value of at inner radius of cyclone inlet - _w shear stress at cyclone wall - angular momentum in cyclone/angular momentum in cyclone inlet - ₁ value of at = ₁ - ₂ value of at = ₂     

Approximate solutions of heat and momentum transfer in laminar plane wall jet

In the present paper approximate solutions for the fluid and thermal boundary layers in an incompressible laminar plane wall jet with isothermal and adiabatic walls have been studied respectively, and comparisons with the known exact solutions have been made wherever possible. It is found that the present method is simple and straightforward, and gives results being in good agreement with the exact solutions. For moderate values of the Prandtl number the method may be used for calculating the heat transfer from an isothermal wall and temperature recovery factor for an adiabatic wall respectively.Nomenclature a* dimensionless temperature gradient at the wall - c _p specific heat at constant pressure - K momentum flux through a cross-section of the jet - Q volume flux through a cross-section of the jet - r* temperature recovery factor - T temperature of the fluid in the boundary layer - T _r adiabatic wall temperature - T temperature of the fluid at rest - u, v velocity components along and normal to the plane wall respectively - x, y rectangular coordinates along and normal to the plane wall respectively - z Greek symbols fluid boundary layer thickness - _t, _T thermal boundary layer thickness for an isothermal and an adiabatic wall respectively - dimensionless y-coordinate - dimensionless temperature difference (T–T )/T - coefficient of thermal conductivity - coefficient of viscosity - coefficient of kinematic viscosity - Prandtl number - _w shearing stress on the plane wall     

Relation between viscoelasticity and shear-thinning behaviour in liquids

Summary The shear-thinning behaviour of a liquid is represented in terms of a relaxation time, defined by the ratio ₀/G₀ of initial viscous and elastic constants. The relationship provides a very simple basis for the evaluation of andG ₀ from viscosity/shear data. Results are compared with relaxation times and moduli from primary normal-stress measurement, from stress relaxation and from direct measurement of recoverable shear strain. Good agreement is found but there is experimental evidence the recoverable shear strain _e is related to normal stressN ₁ and shear stress by _e = N₁/3, which does not agree with the theoretical prediction of eitherWeissenberg orLodge.

---

Zusammenfassung Das Scherentzähungsverhalten einer Flüssigkeit wird mittels einer Relaxationszeit beschrieben, die durch das Verhältnis der Anfangswerte von Viskosität und Elastizitätsmodul ₀/G₀ definiert ist. Diese Beziehung eröffnet eine einfache Methode zur Bestimmung von undG ₀ aus Scherviskositätsmessungen. Die damit erhaltenen Ergebnisse werden mit Relaxationszeiten und Moduln verglichen, die durch Messung der ersten Normalspannungsdifferenz, der Spannungsrelaxation und der Scherdehnungsrückstellung (recoverable shear strain) gewonnen worden sind. Es wird eine gute Übereinstimmung gefunden, zugleich aber wird der experimentelle Nachweis geführt, daß die Scherdehnungsrückstellung _e mit der ersten NormalspannungsdifferenzN ₁ und der Schubspannung durch die Beziehung _e = N₁/3 verknüpft ist, was sowohl zu der theoretischen Voraussage vonWeissenberg als auch zu derjenigen vonLodge im Widerspruch steht.

---

---

With 10 figures and 1 table