Some aspects of dilute suspension theory

A theory proposed by the author as representative of the flow of a general suspension contains three interaction forces, f, S and N. For a quasi-concentrated suspension and for a dilute suspension, N and S, N are omitted, respectively. For the latter special case, we treat diffusion of a fluid through an elastic solid. For a quasi-concentrated suspension, we show that F and S depend on the gradient of the motion gradient. We demonstrate the existence of interesting phenomena: non-simple behavior, dissipative effects, generalized lift and drag forces.Presented at the second conference Recent Developments in Structured Continua, May 23 – 25, 1990, in Sherbrooke, Québec, Canada.     

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²     

Flow of a simple non-newtonian fluid past a sphere

Summary Creeping flow past a sphere is solved for a limiting case of fluid behaviour: an abrupt change in viscosity.List of Symbols d _ij Component of rate-of-deformation tensor - F _d Drag force exerted on sphere by fluid - G ^(d) Coefficients in expression for _ij in terms of d _ij - G _YOJK ^(d) Coefficients in power series representing G ^(d) - R Radius of sphere - r Spherical coordinate - V Velocity of fluid very far from sphere - v _i Component of the velocity vector - x Dimensionless radial distance, r/R - x _i Rectangular Cartesian coordinate - Dimensionless quantity defined by (26) - ^(d) Potential defined by (7) - Value of x denoting border between Regions 1 and 2 as a function of - ₁, ₂ Lower and upper limiting viscosities defined by (10) - Spherical coordinate - * Value of for which =1 - Value of denoting border between regions 1 and 2 as a function of x - Newtonian viscosity - _ij Component of the stress tensor - Spherical coordinate - ₁, ₂ Stream functions defined by (12) and (14) - Second and third invariants of the stress tensor and of the rate-of-deformation tensor, defined by (3)     

Stability criterion for steady-state burning of powders

A model of a powder with a variable surface temperature is considered. A stability criterion for the steady-state burning of powders at constant pressure is found on the assumption that all the processes in the gas phase and the reaction layer of the condensed phase are inertialess. It is shown that the stability region is determined by two parameters: K=(T₁ – t₀) ( ln m/T₀)p and r=(T₁/T₀)p (here T₁ is the surface temperature, T₁ is the initial temperature of the powder, m is the mass burning rate, p is pressure). Burning is always stable if < 1. If > 1, the burning process is stable only when r > ( – 1)²/ / + 1.The author expresses his appreciation to A. S. Kompaneits, O. I. Leipunskii, A. G. Istratov, V. B. Librovich, and S. S. Novikov for discussing his work.     

Shear Instability at the “Explosion Product–Metal” Interface for Sliding Detonation of an Explosive Charge

Periodic perturbations at the explosion product–metal interface were studied experimentally. Experiments were performed for both spherical and plane geometry. Critical conditions of wave formation (detonation velocity of an explosive charge D 6.9 mm/sec) are determined, and an explanation of this effect is given. It is found experimentally that a dynamic pulse causes intense plastic strains at the explosion products–metal interface, leading to thermal softening of the steel boundary layer. In this layer, Kelvin–Helmholtz instability occurs. Calculationanalytical estimates of the critical boundary unstable wavelength agree satisfactorily with experimental results.     

Whose Thoughts Are They, Anyway? Dimensionally Exploding Bion's “Double-Headed Arrow” into Coadapting, Transitional Space

The physics and biology that found psychoanalysis account for discontinuous experience only in the presence of nonmeasurable, metaphysical operators; these include the ego and its subsystems as well as biological experience inherited through Lamarckian principles. Complex, self-organizing systems, however, can link biology to experience without metaphysics. They can also account for psychoanalytically relevant behaviors without appealing to stable internal representations. These behaviors include what W. R. Bion called transformation in O and its corollary, the appearance of the selected fact. By dimensionally exploding the double-headed arrow that he used to link the states Ps and D in his model for thinking (Ps D), we can generate a space that is, at once, psychoanalytically imaginal and dynamically coadapting. Isomorphic to D. W. Winnicott's transitional space, it is self-organizing. It is describable according to dynamics formulated by W J. Freeman, S. Kauffman and C. Langton and it can generate instantaneous conscious contents by way of a selective process analogous to spatio-temporal binding. As a whole, this model supports a clinical stance advanced by D. W. Winnicott as play, within transitional space.     

Stability of elastic bodies under omnilateral compression

Conclusions We have analyzed here the stability of the equilibrium of a simply connected isotropic compressible body with the elastic potential of arbitrary form and under uniform omnilateral deformation. A survey has been given here of earlier results obtained by other authors. The basic celations have been stated in a general form covering the theory of finite subcritical strains and two variants of the theory of small subcritical strains. For the latter theory new relations have been rigorously derived from which perturbations of tracking surface loads can be calculated, on the basis of corresponding expressions in the theory of finite subcritical strains. It has been proven that the sufficient conditions for the applicability of the static method of analysis are satisfied when the same boundary conditions are given over the entire body surface, as well as in several cases of different boundary conditions given at different segments of the boundary surface. It has been shown in a general form, for the theory of finite subcritical strains and for two variants of the theory of small subcritical strains, that the equilibrium of an elastic body under omnilateral deformation is stable, if a tracking load, is given over the entire boundary surface. As an example of problems with different boundary conditions at different segments of the boundary surface, we have considered the conventional problem concerning the stability of a bar on hinge supports and under uniform omnilateral deformation. It has been rigorously proven that in this case the equilibrium is stable when tracking loads are given at the lateral surfaces and is unstable when dead loads are given at the lateral surfaces. These conclusions apply to the theory of finite subcritical strains as well as to the theory of small subcritical strains, and they represent the complete version pertaining to compressible bodies.Institute of Mechanics, Academy of Sciences of the Ukrainian SSR, Kiev. Translated from Prikladnaya Mekhanika, Vol. 12, No. 6, pp. 3–27, June, 1976.     

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

Flow stability of a Grad-model liquid layer on an inclined plane

A study is presented of the flow of stability of a Grad-model liquid layer [1, 2] flowing over an inclined plane under the influence of the gravity force.It is assumed that at every point of the considered material continuum, along with the conventional velocity vector v, there is defined an angular velocity vector , the internal moment stresses are negligibly small, and in the general case the force stress tensor _kj is asymmetric. The model is characterized by the usual Newtonian viscosity , the Newtonian rolling viscosity _r, and the relaxation time = J/4 _r, where J is a scalar constant of the medium with dimensions of moment of inertia per unit mass, is the density. It is assumed that the medium is incompressible, the coefficients , _r, J are constant [2].The exact solution of the equations of motion, corresponding to flow of a layer with a plane surface, coincides with the solution of the Navier-Stokes equations in the case of flow of a layer of Newtonian fluid. The equations for three-dimensional periodic disturbances differ considerably from the corresponding equations for the problem of the flow stability of a layer of a Newtonian medium. It is shown that the Squire theorem is valid for parallel flows of a Grad liquid.The flow stability of the layer with respect to long-wave disturbances is studied using the method of sequential approximations suggested in [3, 4].     

A similarity law for acute circular cones at large angles of attack in a supersonic flow

For thin bodies placed in a hypersonic flow at a small angle of attack the similarity law is known. From this law it follows that for various numbers M, angles of attack , and relative thicknesses the similarity conditions will be observed if in the flows under consideration the parameters M and / are the same. This similarity law is obtained with the assumption M 1, 1. But even for M=3 and 1/3 the results of solving the complete system of gasdynamic equations for affino-similar bodies is in a good agreement with the similarity law [1], In [2] it is shown that this similarity law is generalized for the case of a flow around a thin pointed body at large angles of attack. According to the similarity law, at large angles of attack the flows near bodies with an identical distribution of cross-sectional shapes will be similar if the parameters K₁= cotan and K₂=m sin for all cases have one and the same value. As the angle of attack decreases, the requirements of constancy of K₁ and K₂ become analogous to the conditions M=const, /=const.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 78–83, May–June, 1976.The authors thank V. V. Lunev for the useful discussions and valuable observations.     

Shear softening and thixotropic properties of wheat flour doughs in dynamic testing at high shear strain

Shear softening and thixotropic properties of wheat flour doughs are demonstrated in dynamic testing with a constant stress rheometer. This behaviour appears beyond the strictly linear domain (strain amplitude ₀ 0.2%),G,G and |^*| decreasing with ₀, the strain response to a sine stress wave yet retaining a sinusoidal shape. It is also shown thatG recovers progressively in function of rest time. In this domain, as well as in the strictly linear domain, the Cox-Merz rule did not apply but() and | ^*())| may be superimposed by using a shift factor, its value decreasing in the former domain when ₀ increases. Beyond a strain amplitude of about 10–20%, the strain response is progressively distorted and the shear softening effects become irreversible following rest.     

On the Ellipticity of the Middle Surface of a Shell and its Application to the Asymptotic Analysis of ‘Membrane’ Shells

We consider a surface S = (), where ² is a bounded, connected, open set with a smooth boundary and : ³ is a smooth map; let () denote the components of the two-dimensional linearized strain tensor of S and let 0 with length 0 > 0. We assume the the norm ,|| ()||0, in the space V0() = { H¹() × H¹() × L²(); = 0 on 0 } is equivalent to the usual product norm on this space. We then establish that this assumption implies that the surface S is uniformly elliptic and that we necessarily have 0 = .     

Propagation of a spherical wave in nonlinearly compressible and elastoplastic materials

The problem of spherical wave propagation in soil under the action of an intense uniformly decreasing load ₀(t) applied to the boundary of a cavity with radius r₀ is considered. Soil with a high stress level is modeled either by ideally nonlinearly compressible or elastoplastic material, taking account of linear irreversible unloading for the material. In contrast to [1–7], in order to describe material movement use is made of strain theory [8] with determining functions = (), _i=_i(_i), where , _i, , _i are the first and second invariants of strain and stress tensors. During material loading these functions are presented in the form of polynomials ()=(_i+₂¦¦), _ii)=(_i-_2i)_i, in which constant coefficients _i, _i=1, 2) are determined by experiment, taking account of the triaxial stressed state of soil. Solution of the problem is constructed by an analytically reversible method, with prescribed shape for the shock-wave (SW) surface in the form of a second-degree polynomial relating to time t and a numerical method of characteristics for a prescribed arbitrarily decreasing load _i(t). On the basis of the analytical equations obtained, calculations are carried out for material parameters (including loading profile) in a computer and stresses and mass velocity of plastic and elastoplastic materials are compared.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 95–100, July–August, 1986.The authors express their sincere thanks to Kh. A. Rakhmatulin for discussing the results of this work.     

A rotating hot-wire technique for spatial sampling of disturbed and manipulated duct flows

The documentation and control of flow disturbances downstream of various open inlet contractions was the primary focus with which to evaluate a spatial sampling technique. An X-wire probe was rotated about the center of a cylindrical test section at a radius equal to one-half that of the test section. This provided quasi-instantaneous multi-point measurements of the streamwise and azimuthal components of the velocity to investigate the temporal and spatial characteristics of the flowfield downstream of various contractions. The extent to which a particular contraction is effective in controlling ingested flow disturbances was investigated by artificially introducing disturbances upstream of the contractions. Spatial as well as temporal mappings of various quantities are presented for the streamwise and azimuthal components of the velocity. It was found that the control of upstream disturbances is highly dependent on the inlet contraction; for example, reduction of blade passing frequency noise in the ground testing of jet engines should be achieved with the proper choice of inlet configurations.List of symbols K _uv correlation coefficient= - P percentage of time that an azimuthal fluctuating velocity derivative dv/d is found - U streamwise velocity component U=U (, t) - V azimuthal or tangential velocity component due to flow and probe rotation V=V (, t) - mean value of streamwise velocity component - U _m resultant velocity from and - mean value of azimuthal velocity component induced by rotation - u fluctuating streamwise component of velocity u=u(, t) - v fluctuating azimuthal component of velocity v = v (, t) - u phase-averaged fluctuating streamwise component of velocity u=u(0) - v phase-averaged fluctuating azimuthal component of velocity v=v() - û average of phase-averaged fluctuating streamwise component of velocity (u()) over cases I-1, II-1 and III-1 û = û() - average of phase-averaged fluctuating azimuthal component of velocity (v()) over cases I-1, II-1 and III-1 - u fluctuating streamwise component of velocity corrected for non-uniformity of probe rotation and/or phase-related vibration u = u(0, t) - v fluctuating azimuthal component of velocity corrected for non-uniformity or probe rotation and/or phase-related vibration v=v (, t) - u ² rms value of corrected fluctuating streamwise component of velocity - rms value of corrected fluctuating azimuthal component of velocity - phase or azimuthal position of X-probe     

On surface waves in generalized thermoelasticity

Knowles' representation theorem for harmonically time-dependent free surface waves on a homogeneous, isotropic elastic half-space is extended to include harmonically time-dependent free processes for thermoelastic surface waves in generalized thermoelasticity of Lord and Shulman and of Green and Lindsay.r , , r , , .This work was done when author was unemployed.     

Measurements of slip at the wall during flow of high-density polyethylene through a rectangular conduit

A hot-film probe has been used to measure slip of high-density polyethylene flowing through a conduit with a rectangular cross section. A transition from no slip to a stick-slip condition has been observed and associated with irregular extrudate shape. Appreciable extrudate roughness was initiated at the same flow rate as that at which the relationship between Nusselt number and Péclet number for heat transfer from the probe departed from the behavior expected for a no-slip condition at the conduit wall. A ₁ constant defined by eq. (A3) - C dimensionless group used in eq. (7) - C _p heat capacity - D constant in eq. (13) - f u _s/u - f _lin defined by eq. (A6) - G storage modulus - G loss modulus - k thermal conductivity - L length of hot film in thex-direction - L _eff effective length of large probe found from eq. (A3) - Nu _L Nusselt number, defined for a lengthL by eq. (2) - (Nu _L)₀ value ofNu _L atPe = 0 (eq. (A 1)) - Pe Péclet number,uL/ - Pe ₀ Péclet number in slip flow, eq. (6) - Pe ₁ Péclet number in shear flow, eq. (4) - q _L average heat flux over hot film of lengthL - R _i resistances defined by figure 8 - R _AB correlation coefficient defined by eq. (14) for signalsA andB - T temperature - T _s temperature of probe surface - T ambient temperature - T T _s –T - u average velocity - u _s slip velocity - V _b voltage indicated in figure 8 - W probe dimension (figure 6) - x distance in flow direction (figure 1) - y distance perpendicular to flow direction (figure 1) - thermal diffusivity,k/C _p - wall shear rate - _5% thickness of lubricating layer during probe calibration for introduction of an error no greater than 5%; see Appendix I - ^* complex viscosity - density - time - _c critical shear stress, eq. (13) - _w wall shear stress - frequency characterizing extrudate distortion (figures 12 and 13), or frequency of oscillation during rheometric characterization (figures 18–20) - * quantity obtained from normalized Nusselt number, eq. (A1), or complex viscosity ^* - A actual (small) probe (see Appendix I) - M model (large) probe (see Appendix I)     

Investigation of stationary flow of viscous gas in a thin three-dimensional shock layer

A three-dimensional shock layer near the blunt surface of a fairly smooth body is analyzed asymptotically. Equations of the first approximation are obtained and justified in various cases of the limit 1, 0, ( – 1)^–1M ^-2 0. These equations are simplified for the flow near the stagnation point of a body with double curvature and near the blunt leading edge of a sweptback wing. The results of some calculations are given and compared with the results of [17, 18] in the case of axisymmetric flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 115–126, September–October, 1980.     

Analysis of turbulent boundary layer interaction with an external supersonic flow behind a step

An integral method of analyzing turbulent flow behind plane and axisymmetric steps is proposed, which will permit calculation of the pressure distribution, the displacement thickness, the momentum-loss thickness, and the friction in the zone of boundary layer interaction with an external ideal flow. The characteristics of an incompressible turbulent equilibrium boundary layer are used to analyze the flow behind the step, and the parameters of the compressible boundary layer flow are connected with the parameters of the incompressible boundary layer flow by using the Cowles-Crocco transformation.A large number of theoretical and experimental papers devoted to this topic can be mentioned. Let us consider just two [1, 2], which are similar to the method proposed herein, wherein the parameter distribution of the flow of a plane nearby turbulent wake is analyzed. The flow behind the body in these papers is separated into a zone of isobaric flow and a zone of boundary layer interaction with an external ideal flow. The jet boundary layer in the interaction zone is analyzed by the method of integral relations.The flow behind plane and axisymmetric steps is analyzed on the basis of a scheme of boundary layer interaction with an external ideal supersonic stream. The results of the analysis by the method proposed are compared with known experimental data.Notation x, y longitudinal and transverse coordinates - X, Y transformed longitudinal and transverse coordinates - , ^*, ^** boundary layer thickness, displacement thickness, momentum-loss thickness of a boundary layer - , ^*, ^** layer thickness, displacement thickness, momentum-loss thickness of an incompressible boundary layer - u, velocity and density of a compressible boundary layer - U, velocity and density of the incompressible boundary layer - , stream function of the compressible and incompressible boundary layers - , dynamic coefficient of viscosity of the compressible and incompressible boundary layers - r₁ radius of the base part of an axisymmetric body - r radius - R transformed radius - M Mach number - friction stress - p pressure - a speed of sound - s enthalpy - v Prandtl-Mayer angle - P Prandtl number - P_t turbulent Prandtl number - r₂ radius of the base sting - b step depth - =0 for plane flow - =1 for axisymmetric flow Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 33–40, May–June, 1971.In conclusion, the authors are grateful to M. Ya. Yudelovich and E. N. Bondarev for useful comments and discussions.     

A note on three-dimensional boundary layer growth

Summary This note is an extension of the work of Görtler²) on two-dimensional boundary layer growth to the three-dimensional case. The solutions of three-dimensional boundary layer equations are obtained by considering the potential flow of the body to be governed by the functions At U ₀(, ) and At U ₀(, ) where is any positive number.     

Die Viskosität reiner Flüssigkeiten in der Nähe der kritischen Temperatur

Zusammenfassung In einem Druck-Temperatur-Diagramm mit der Viskosität als Parameter sind die Isoviskosen reiner Flüssigkeiten in der Nähe der kritischen Temperatur in guter Näherung Geraden. Unter Berücksichtigung dieses phänomenologischen Sachverhaltes sowie einiger Gesetzmäßigkeiten am kritischen Punkt gelingt es, eine Gleichung für dieses Zustandsgebiet aufzustellen, in der nur zwei Konstanten den Meßwerten angepaßt werden müssen. Diese Gleichung wird auf Wasser und Para-Wasserstoff angewandt. Bei Wasser gibt die Gleichung im Bereich von 300 °C bis 374 °C, in dem bisher keine international anerkannte Interpolationsgleichung existiert, die Werte der internationalen Rahmentafel [1] innerhalb deren Toleranz wieder. Sie ist daher als Interpolationsgleichung für diesen Bereich geeignet.

---

The viscosity of pure liquids near the critical temperature

In a pressure-temperature-viscosity-diagram, the lines of constant viscosity of pure liquids near the critical temperature are straight to a good approximation. This phenomenological fact as well as some laws at the critical point lead to an equation for the viscosity in which only two constants have to be adjusted to the experimental values. The equation is evaluated for water and parahydrogen. The results for water are within the limits of the International Seeleton Tables in the region between 300 °C and 374 °C where no international equation for interpolation exists.

---

Formelzeichen T Temperatur - T^* Grenztemperatur für Bereich I - T' Grenztemperatur für die Restviskosität - T_kr kritische Temperatur - T_N normierte Temperatur; T_N=T/T_kr - A₀(T) Viskosität auf der Siedelinie - A₁(T) Steigung der Isothermen im ,p, T-Bild - p Druck - P_kr kritischer Druck - P_N normierter Druck: p_N = p/p_kr - Psätt Druck auf der Siedelinie bei der TemperaturT - Dichte - _kr kritische Dichte - _N normierte Dichte; _N=/kr - Viskosität - _kr kritische Viskosität - _kr normierte Viskosität; _N=/_kr - Restviskosität - (_N) Funktion in Gl. (3) - (_N) Punktion in Gl. (3) - Ri Riedel-Zahl nach Gl. (4) - spezifisches Volum