Torsion of a hemisphere embedded in an elastic half-space

This paper deals with the stress distribution in a homogeneous isotropic elastic hemisphere embedded in a semi-infinite homogeneous isotropic elastic medium when a rigid circular disc is clamped to the plane face of the hemisphere and the stresses are caused by the rotation of the disc through an angle . The problem is reduced to the solution of a Fredholm integral equation of the second kind in the auxiliary function (t). An analytical expression for the torque T required to rotate the die through an angle is obtained in terms of (t). The Fredholm integral equation is solved numerically, and the numerical values of T are graphed.This work has been supported by the National Research Council of Canada through NRC-Research Grant No. A4177.     

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.     

Filament stretching equation in Lagrangian coordinates

The general integral of the very simple equation ^21/n/() was found to describe the cross sectional area of filaments of isothermal power law fluids while in transient stretching where is time and is the initial location of fluid molecules at time = 0 given as the distance from a reference point fixed in space. Any such stretching transient given as a solution of the above equation is physically realizable subject to the restrictions > 0 and/ < 0.     

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     

The effect of background particulates on the delayed transition of a heated 9:1 ellipsoid

An experimental study was done to quantify the effects of a variety of background particulates on the delayed laminar-turbulent transition of a thermally stabilized boundary layer in water. A Laser-Doppler Velocimeter system was used to measure the location of boundary layer transition on a 50 mm diameter, 9:1 fineness ratio ellipsoid. The ellipsoid had a 0.15 m RMS surface finish. Boundary layer transition locations were determined for length Reynolds numbers ranging from 3.0 × 10⁶ to 7.5 × 106. The ellipsoid was tested in three different heating conditions in water seeded with particles of four distinct size ranges. For each level of boundary layer heating, measurements of transition were made for clean water and subsequently, water seeded with 12.5 m, 38.9 m, 85.5 m and 123.2 m particles, alternately. The three surface heating conditions tested were no heating, T = 10°C and T = 15°C where T is the difference between the inlet model heating water temperature, T _i, and free stream water temperature, T . The effects of particle concentration were studied for 85.5 m and 123.2 m particulates.The results of the study can be summarized as follows. The 12.5 m and 38.9 m particles has no measurable effect on transition for any of the test conditions. However, transition was significantly affected by the 85.5 m and 123.2 m particles. Above a length Reynolds number of 4 × 10⁶ the boundary layer transition location moved forward on the body due to the effect of the 85.5 m particles for all heating conditions. The largest percentage changes in transition location from clean water, were observed for 85.5 m particles seeded water.Transition measurements made with varied concentrations of background particulates indicated that the effect of the 85.5 m particles on the transition of the model reached a plateau between 2.65 particulates/ml concentration and 4.2 particles/ml. Measurements made with 123.3 m particles at concentrations up to 0.3 part/ml indicated no similar plateau.     

Effect of Streamwise Streaky Structures on Turbulization of a Circular Jet

Experimental investigations of the influence of streamwise streaky structures on turbulization of a circular laminar jet are described. The qualitative characteristics of jet evolution are studied by smoke visualization of the flow pattern in the jet and by filming the transverse and longitudinal sections of the jet illuminated by the laser sheet with image stroboscopy. It is shown that the streaky structures can be generated directly at the nozzle exit, and their interaction with the Kelvin–Helmholtz ring vortices leads to emergence of azimuthal beams ( structures) by a mechanism similar to threedimensional distortion of the twodimensional Tollmien–Schlichting wave at the nonlinear stage of the classical transition in nearwall flows. The effect of the jetexhaustion velocity and acoustic action on jet turbulization is considered.     

Mathematics and Rheology

Summary The great advances made in Rheology during the last forty years owe much to the application of mathematics. But in some cases, there have been misunderstandings. As a result, rheologists have sometimes been unnecessarily restricted. Mathematics is limited only by self-consistency: Rheology deals with the physical world.In particular, three terms are discussed: infinity, zero and negative, all of which have at least two meanings which have been confused. Rheological phenomena cannot be explained merely by mathematical formulation. Unnecessary mathematics should be avoided but professional rheologists must not shirk mathematics.This paper was read in an abbreviated form, at aConference of the British Society of Rheology in April 1971 The Proceedings were not published.     

Interpolation formulas for maxwell viscosity of certain metals as a function of shear-strain intensity and temperature

The purpose of this study is the construction of interpolation formulas for the dependence of Maxwell viscosity, a quantity which is the reciprocal of shear-strain relaxation time , on shear-strain intensity and temperature for several metals: iron, aluminum, copper, and lead. This function was interpolated in various temperature and deformation velocity ranges in accordance with available experimental data for iron (0 10⁷ sec^–1, 200 ° T 1500 °); aluminum (0 10⁷ sec^–1, 300 ° T 900 °); copper (0 10⁵ sec^–1, 300 ° T 1300 °); lead (0 10⁶ sec^–1, 90 ° T 400 °); temperatures in °K.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 114–118, July–August, 1974.     

On laminar flow through a uniformly porous pipe

Numerous investigations ([1] and [4–9]) have been made of laminar flow in a uniformly porous circular pipe with constant suction or injection applied at the wall. The object of this paper is to give a complete analysis of the numerical and theoretical solutions of this problem. It is shown that two solutions exist for all values of injection as well as the dual solutions for suction which had been noted by previous investigators. Analytical solutions are derived for large suction and injection; for large suction a viscous layer occurs at the wall while for large injection one solution has a viscous layer at the centre of the channel and the other has no viscous layer anywhere. Approximate analytic solutions are also given for small values of suction and injection.

Nomenclature

General r distance measured radially - z distance measured along axis of pipe - u velocity component in direction of z increasing - v velocity component in direction of r increasing - p pressure - density - coefficient of kinematic viscosity - a radius of pipe - V velocity of suction at the wall - r ²/a ² - R wall or suction Reynolds number, Va/ - f() similarity function defined in (6) - u ₀() eigensolution - U(0) a velocity at z=0 - K an arbitrary constant - B _K Bernoulli numbers Particular Section 5 perturbation parameter, –2/R - ² a constant, –K - x / - g(x) f()/ Section 6 perturbation parameter, –R/2 - ² a constant, –K - g() f() - g _c ()=g() near centre of pipe - * point where g()=0 Section 7 2/R - ² K - t (1–)/ - w(t, ) [1–f(t)]/ - ₀, ₁ constants - g() f()– ₀ - ₀/ - ₀ a constant - * point where f()=0     

Second-Order Estimates for the Macroscopic Response and Loss of Ellipticity in Porous Rubbers at Large Deformations

This work presents the application of a recently proposed second-order homogenization method (Ponte Castañeda, 2002) to generate estimates for effective behavior and loss of ellipticity in hyperelastic porous materials with random microstructures that are subjected to finite deformations. The main concept behind the method is the introduction of an optimally selected linear thermoelastic comparison composite, which can then be used to convert available linear homogenization estimates into new estimates for the nonlinear hyperelastic composite. In this paper, explicit results are provided for the case where the matrix is taken to be isotropic and strongly elliptic. In spite of the strong ellipticity of the matrix phase, the homogenized second-order estimates for the overall behavior are found to lose ellipticity at sufficiently large compressive deformations corresponding to the possible development of shear band-type instabilities (Abeyaratne and Triantafyllidis, 1984). The reasons for this result have been linked to the evolution of the microstructure, which, under appropriate loading conditions, can induce geometric softening leading to overall loss of ellipticity. Furthermore, the second-order homogenization method has the merit that it recovers the exact evolution of the porosity under a finite-deformation history in the limit of incompressible behavior for the matrix. Mathematics Subject Classifications (2000) 49S05, 74B20, 74Q15, 74Q05.     

Calculation of the pressure on the side surface of bodies consisting of a spherical segment joined to an inverted cone in a supersonic stream

The results of investigations of inviscid flow over inverted cones with nose consisting of a spherical segment were published for the first time in Soviet literature in [1–4]. In the present paper, a numerical solution to this problem is obtained using the improved algorithms of [5, 6], which have proved themselves well in problems of exterior flow over surfaces with positive angles of inclination to the oncoming flow. It is shown that the Mach number 2 M , equilibrium and nonequilibrium physicochemical transformations in air (H = 60 km, V = 7.4 km/sec, R₀ = 1 m), and the angle of attack 0 40° influence the investigated pressure distributions. A comparison of the results of the calculations with drainage experiments for M = 6, = 0-25° confirms the extended region of applicability of the developed numerical methods. Also proposed is a simple correlation of the dependence on the Mach number in the range 1.5 M of the shape of the shock wave near a sphere in a stream of ideal gas with adiabatic exponent = 1.4.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 178–183, January–February, 1981.     

Closure of microcracks during unloading and the formation of reverse dislocations

Conclusion During unloading, local closure of the crack is followed by the emission of the first pair of reverse dislocations. The formation of an area of local contact during unloading occurs with a bang; the critical SIF of the external loads at this moment coincides with the value of the overall SIF of the field of the primary sessile dislocations. There is always a free gap between the crack tip and the area of local contact. External loads acting in the crack plane and not resulting in a stress concentration at the crack front play an important role in the emission of reverse dislocations (in particular, they may block their departure from the S-structure).Pacific Oceanological Institute, Academy of Sciences of the USSR, Vladivostok. Translated from Prikladnaya Mekhanika, Vol. 24, No. 7, pp. 3–18, July, 1988.     

Effect of Dynamic Prehistory on Aerodynamics of a Laminar Separated Flow in a Channel Behind a Rectangular Backward‐Facing Step

The effect of dynamic prehistory of the flow and the channelexpansion ratio on aerodynamics of a steady separated laminar flow behind a rectangular backwardfacing step located in a planeparallel channel is numerically studied. It is shown that the boundary layer upstream of the flow separation exerts a strong effect on flow characteristics behind the step. A decrease in the boundarylayer thickness in the cross section of the step leads to a decrease in the separationregion length, and an increase in the channelexpansion ratio with a fixed initial boundarylayer thickness and Reynolds number leads to an increase in the separationregion length.     

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     

Extended Korn's inequalities and the associated best possible constants

We are concerned with the coerciveness of the strain energy E(u) (in linear elasticity) associated with a displacement vector u on the Sobolev space H¹ () or its subspaces, a domain in ⁿ representing an isotropic elastic body—certain specific cases are called Korn's inequalities. Sufficient (and necessary) conditions on the Lamé moduli for E(·) to be coercive (or uniformly positive) on such spaces are given, and the associated best possible constants are obtained for some cases.     

EXACT SOLUTION TO THE PROBLEM OF FLOW OF SLIGHTLY COMPRESSIBLE FLUIDS IN A BOUNDED CONFINED “PRACTURE-PORE” MEDIUM AND ITS APPLICATION TO WELL TESTING

The problem of flow of slightly compressible fluids through a bounded confined fracture-pore medium is solved and studied thoroughly in this paper. Some essential natures of flow of elastic liquids through a medium with double porosity under the condition of neglecting the flow in matrix system were revealed and clarified further. The method of estimating all parameters commonly interested in a bounded confined fracture-pore medium reservoir through a series of flow tests in wells by use of the solution obtained is presented.     

Normal forms for random diffeomorphisms

Given a dynamical system (,, ,) and a random diffeomorphism (): ^{d d} with fixed point at x=0. The normal form problem is to construct a smooth near-identity nonlinear random coordinate transformation h() to make the random diffeomorphism ()=h()^–1() h() as simple as possible, preferably linear. The linearization D(, 0)=:A() generates a matrix cocycle for which the multiplicative ergodic theorem holds, providing us with stochastic analogues of eigenvalues (Lyapunov exponents) and eigenspaces. Now the development runs pretty much parallel to the deterministic one, the difference being that the appearance of turns all problems into infinite-dimensional ones. In particular, the range of the homological operator is in general not closed, making the conceptof-normal form necessary. The stochastic versions of resonance and averaging are developed. The case of simple Lyapunov spectrum is treated in detail.     

Transverse Decay of Solutions in an Elastic Cylinder

In a semiinfinite cylinder composed of anisotropic linearised elastic material, loaded on the base and clamped along the lateral surface, it is known that the solution as measured, for example, by the strainenergy flux through a plane crossection decays longitudinally at most exponentially with respect to the axial distance from the base. There is, however, also a transverse radial decay of the solution, again measured for example by the strainenergy, occurring from the region close to the cylinder's axis to the region near the lateral surface, where the energy vanishes.This problem is considered in the present paper which discusses a circular semiinfinite cylinder and derives an estimate for the strainenergy contained in a cylindrical annulus at a given distance from the base and of variable height, and whose outer surface coincides with the lateral surface of the cylinder. It is shown that the strainenergy decays at most algebraically to zero as the inner radius of the annulus increases to that of the cylinder.Sommario. E'noto che in un cilindro semiinfinito composto da materiale elastico lineare anisotropo, caricato sulla base ed incastrato lungo la superficie laterale, la soluzione elastica, misurata, per esempio, dal flusso di energia di deformazione attraverso una sezione trasversale piana, decade con legge al più esponenziale con la distanza dalla base. C'è tuttavia, anche un decadimento radiale della soluzione, misurato, per esempio, dall'energia di deformazione che passa dalla regione vicina all'asse del cilindro a quella vicin alla superficie laterale dove l'energia si annullaQuesto problema è qui studiato. Si discute in particolare un cilindro circolare semiinfinito e si deduce una stima per l'energia di deformazione contenuta in un anello cilindrico ad una distanza assegnata dalla base e di altezza variabile, e la cui superficie esterna coincide con la superficie laterale del cilindro. Si dimostra che l'energia di deformazione decade al più con legge algebrica a zero quando il raggio interno del cilindro si avvicina a quello esterno.     

A matrix theory of elastic-locking structures

Summary The theory developed exhibits the following peculiar features: structures are discretizised in finite elements, the constitutive laws piecewise linearized, the problem is split in a preliminary linear elastic solution and a corrective nonlinear subproblem; concepts and techniques of quadratic and linear programming theory are utilized. The main results are: for the analysis under given loads and dislocations, a pair of extremum theorems for locking stresses, corresponding to dual quadratic programs; for the limit analysis with respect to locking situations two already known theorems, which are here deduced from the solvability conditions of the above quadratic programs and formulated as dual linear programs. The extension of the results to imperfectly locking behavior is carried out. Some examples illustrate the solution techniques based on the theory expounded.

---

Sommario Caratteristiche peculiari della teoria sviluppata sono: discretizzazione della struttura in elementi finiti, linearizzazione a tratti delle leggi costitutive, suddivisione del problema in una preliminare soluzione elastica lineare ed in una fase nonlineare correttiva, utilizzazione di concetti e tecniche della teoria delle programmazioni quadratiche e lineari.Si ottengono, per l'analisi di continui a deformabilità limitata soggetti a dati carichi e distorsioni, due teoremi di estremo corrispondenti a programmi quadratici negli sforzi di bloccaggio; per l'analisi limite al bloccaggio del sistema, due teoremi già noti, che qui si deducono dalle condizioni di solubilità dei precedenti programmi quadratici e si presentano in veste di programmi lineari duali. I risultati ottenuti vengono generalizzati ai casi definibili a deformabilità imperfettamente limitata. Alcuni esempi illustrano le tecniche risolutive che sorgono dalla teoria svolta.

---

---

The study presented here forms part of a research program supported by the C.N.R.