Renormalization Group Model of Large-Scale Turbulence in Porous Media

Renormalization group methods are used to develop a macroscopic (large-scale) turbulence model for incompressible flow in porous media. The model accounts for the large-distance and large-time behavior of velocity correlations generated by the momentum equation for a randomly stirred, incompressible flow. Utilizing the renormalization procedure, the transport equations for the large-scale modes and expressions for effective transport coefficients are obtained. Expressions for renormalized turbulent viscosity, which accounts for the ultraviolet subrange of the turbulent kinetic energy spectrum, are also obtained.     

Parametric stress-concentration factors in greenstick bone fracture

The reported work is a part of an ongoing research program concerned with structural analysis of fractured long bone and methods of internal fixation. The stress-concentration factors for equine metacarpus bones containing greenstick fractures and “through” fractures (surgically repaired) were determined for the compression, flexural and torsional modes of loading based on whole bone (unfractured) strengths. The greenstick type of fracture was simulated with saw cuts at the mid-span of the bone, and the parameters varied were depth of fracture and orientation of fracture. All specimens consisted of fresh dead bone which had been placed in a freezer within 4 hr after expiration. The maximum stress-concentration factors for the simulated greenstick fractures studied were about 3.4 for compression, 4.3 for torsion and 16 for flexure. The stressconcentration factors for fractured bones surgically repaired with commercial plates were about 3.0 for compression, 2.7 for torsion and 6.1 for flexure.     

Equivalent linearizations for practical hysteretic systems

Experimental testing of a friction damped base isolation system has indicated a need for a new model of friction damping and for an appropriate equivalent linearization technique. The model for the damping adopted is a combination of viscous damping, constant Coulomb friction and linear Coulomb friction.This model is incorporated into the equation of motion for a single-degree-of-freedom system and the exact solutions are given for free vibrations and for steady-state vibrations excited by a harmonic force. The exact solution is taken as a basis for an equivalent linearization technique that can be used in conjunction with conventional design spectra for a practical design of such a system.     

Singularities in 2D anisotropic potential problems in multi-material corners: Real variable approach

An analysis of singular solutions at corners consisting of several different homogeneous wedges is presented for anisotropic potential theory in plane. The concept of transfer matrix is applied for a singularity analysis first of single wedge problems and then of multi-material corner problems. Explicit forms of eigenequations for evaluation of singularity exponent in the case of multi-material corners are derived both for all combinations of homogeneous Neumann and Dirichlet boundary conditions at faces of open corners and for multi-material planes with singular interior points. Perfect transmission conditions at wedge interfaces are considered in both cases. It is proved that singularity exponents are real for open anisotropic multi-material corners, and a sufficient condition for the singularity exponents to be real for anisotropic multi-material planes is deduced. A case of a complex singularity exponent for an anisotropic multi-material plane is reported, apparently for the first time in potential theory. Simple expressions of eigenequations are presented first for open bi-material corners and bi-material planes and second for a crack terminating at a bi-material interface, as examples of application of the theory developed here. Analytical solutions of these eigenequations are presented for interface cracks with any combination of homogeneous boundary conditions along the interface crack faces, and also for a special case of a crack perpendicular to a bi-material interface. A numerical study of variation of the singularity exponent as a function of inclination of a crack terminating at a bi-material interface is presented.     

Green's function for SH-waves in a cylindrically monoclinic material

Green's function for SH-waves in a cylindrically monoclinic material is considered for impulsive and time-harmonic sources. Closed form expressions for the Green's function are derived for a few limited values of anisotropic parameters. A very interesting time development of the wave front shape is illustrated and the wave front singularity is discussed for the transient SH-wave. Contours of the displacement amplitude for the time-harmonic wave are also shown.     

The re-examination of the weakly nonlinear theory of hydrodynamic stability

The weakly nonlinear theory has been widely applied in the problem of hydrodynamicstability and also in other fields.However,although its application has been successful forsome problems,yet,for other problems,the results obtained are not satisfactory,especiallyfor problems like transition or the evolution of the vortex in the free shear flow,for whichthe goal of the theoretical investigation is not seeking for a steady state,but predicting anevolutional process.In this paper,we shall examine the reason for the unsuccessfulness andsuggest ways for its amendment.     

Calculation of the shock-wave adiabatics for some heterogeneous mixtures

A method has been proposed [1] for calculating shock-wave adiabatics for two-phase mixtures which does not use the additivity principle, and this is here extended to a medium consisting of n phases; the shock-wave adiabatic is derived and compared with experiment for a three-phase mixture; an alloy consisting of copper, zinc, and lead (brass). It is shown for a paraffin-tungsten mixture that the adiabatics for the individual phases deviate from ones for the corresponding solid materials. The known equations of state for aluminum and epoxide resin have been used to determine the shock-wave adiabatic for such a mixture. A comparison is made with the analogous calculation made from the additivity principle.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 130–137, July–August, 1973.I am indebted to V. N. Nikolaevskii for suggesting the topics and for useful discussion.     

Mechanical energy balances for a capillary viscometer

The entrance and exit flow processes for a cylindrical geometry are analyzed by writing macroscopic mechanical energy balances for a capillary viscometer. These equations can be used to compute the entrance and exit excess dissipation integrals from measured pressure differences if viscometric normal stress data are available for the material of interest. Upper and lower bounds are derived for these integrals for cases when high shear rate normal stress data are not available. The utilization of macroscopic mechanical energy balances in the interpretation of capillary viscometer results is illustrated using numerical solutions for a Maxwell fluid and experimental pressure drop data for high density polythylene.     

A sharp condition for existence of an inertial manifold

It is shown that a perturbation argument that guarantees persistence of inertial (invariant and exponentially attracting) manifolds for linear perturbations of linear evolution equations applies also when the perturbation is nonlinear. This gives a simple but sharp condition for existence of inertial manifolds for semi-linear parabolic as well as for some nonlinear hyperbolic equations. Fourier transform of the explicitly given equation for the tracking solution together with the Plancherel's theorem for Banach valued functions are used.     

On unsteady unidirectional flows of a second grade fluid

Some properties of unsteady unidirectional flows of a fluid of second grade are considered for flows produced by the sudden application of a constant pressure gradient or by the impulsive motion of one or two boundaries. Exact analytical solutions for these flows are obtained and the results are compared with those of a Newtonian fluid. It is found that the stress at the initial time on the stationary boundary for flows generated by the impulsive motion of a boundary is infinite for a Newtonian fluid and is finite for a second grade fluid. Furthermore, it is shown that initially the stress on the stationary boundary, for flows started from rest by sudden application of a constant pressure gradient is zero for a Newtonian fluid and is not zero for a fluid of second grade. The required time to attain the asymptotic value of a second grade fluid is longer than that for a Newtonian fluid. It should be mentioned that the expressions for the flow properties, such as velocity, obtained by the Laplace transform method are exactly the same as the ones obtained for the Couette and Poiseuille flows and those which are constructed by the Fourier method. The solution of the governing equation for flows such as the flow over a plane wall and the Couette flow is in a series form which is slowly convergent for small values of time. To overcome the difficulty in the calculation of the value of the velocity for small values of time, a practical method is given. The other property of unsteady flows of a second grade fluid is that the no-slip boundary condition is sufficient for unsteady flows, but it is not sufficient for steady flows so that an additional condition is needed. In order to discuss the properties of unsteady unidirectional flows of a second grade fluid, some illustrative examples are given.     

Development of co-current air–water flow in a vertical pipe

Measurements of the cross-sectional distribution of the gas fraction and bubble size distributions were conducted in a vertical pipe with an inner diameter of 51.2 mm and a length of about 3 m for air/water bubbly and slug flow regimes. The use of a wire-mesh sensor obtained a high resolution of the gas fraction data in space as well as in time. From this data, time averaged values for the two-dimensional gas fraction profiles were decomposed into a large number of bubble size classes. This allowed the extraction of the radial gas fraction profiles for a given range of bubble sizes as well as data for local bubble size distributions. The structure of the flow can be characterized by such data. The measurements were performed for up to 10 different inlet lengths and for about 100 combinations of gas and liquid volume flow rates. The data is very useful for the development and validation of meso-scale models to account for the forces acting on a bubble in a shear liquid flow and models for bubble coalescence and break-up. Such models are necessary for the validation of CFD codes for the simulation of bubbly flows.     

An integral representation of the surface-impedance tensor for incompressible elastic materials

A fundamental result in anisotropic elasticity and surface-wave theory is the integral representation for the surface-impedance tensor first derived by Barnett and Lothe in 1973. However, this representation is only valid for compressible materials but not valid for incompressible materials. In this paper the corresponding integral representation for the surface-impedance tensor valid for incompressible materials is derived and is used to establish the uniqueness of surface-wave speed and to obtain an expression for the tensor Green's function for the infinite space. Mathematics subject classifications (2000) 74B05, 74B15, 74B20, 74J15     

Indentation by nominally flat or conical indenters with rounded corners

Axisymmetric indentation of a flat surface is considered: specifically, the case of flat-ended indenter with rounded edges, and the case of a shallow cone with a rounded tip. Analytical solutions are obtained for the normal and sequential tangential loading, in both full or partial slip conditions (with the Cattaneofn9Mindlin approximation) , and for the complete interior stress field in all these conditions.Implications for strength of the contact are discussed with reference to metallic or brittle materials, with the intention to shed more light in particular to the understanding of common fretting fatigue or indentation testings with nominally flat or conical indenters. It is found that the strength of the contact, which is nominally zero for perfectly sharp flat or conical indenters, is well defined even for a small radius of curvature. This is particularly true for the flat indenter, for which the strength is even significantly higher than for the classical Hertzian indenter for a wide range of geometrical and loading conditions, rendering it very attractive for design purposes.     

Rotating beams and non-rotating beams with shared eigenpair for pinned-free boundary condition

In this paper we look for a rotating beam, with pinned-free boundary conditions, whose eigenpair (frequency and mode-shape) is same as that of a uniform non-rotating beam for a particular mode. It is seen that for any given mode, there exists a flexural stiffness function (FSF) for which the ith mode eigenpair of a rotating beam with uniform mass distribution, is identical to that of a corresponding non-rotating beam with same length and mass distribution. Inserting these derived FSF’s in a finite element code for a rotating pinned-free beam, the frequencies and mode shapes of a non-rotating pinned-free beam are obtained. For the first mode, a physically realistic equivalent rotating beam is possible, but for higher modes, the FSF has internal singularities. Strategies for addressing these singularities in the FSF for finite element analysis are provided. The proposed functions can be used as test functions for rotating beam codes and also for targeted destiffening of rotating beams.     

Stress-concentration factors for multiple semielliptical notches in beams under pure bending

Uniformly spaced semielliptical notches in rectangular beams under pure bending are examined photoelastically. Minimum stress-concentration factors, produced by multiple elliptical notching of beams, are obtained for wide ranges of notch width, semiminor elliptical axis, notch depth, notch pitch and depth of beam. In particular, the geometries of the optimum elliptical notch producing the least stress concentration are obtained for a practical range of parameters. Stress-concentration factors for beams with multiple semi-elliptical notches are compared to those for beams with single semielliptical notches and to those for beams with semicircular notches. The maximum reduction of stress-concentration factor for beams with multiple semielliptical notches is to approximately 37 percent of the stress-concentration factor for beams with single semielliptical notches. Within the range of parameters investigated, the stress-concentration factor for beams with multiple semielliptical notches was 15 to 37 percent less than that for multiple semicircular notches.     

Stochastic equations for continuum and determination of hydraulic drag coefficients for smooth flat plate and smooth round tube with taking into account intensity and scale of turbulent flow

The stochastic equations of continuum are used for determining the hydraulic drag coefficients. As a result, the formulas for the hydraulic drag coefficients dependent on the turbulence intensity and scale instead of only on the Reynolds number are proposed for the classic flows of an incompressible fluid along a smooth flat plate and a round smooth tube. It is shown that the new expressions for the classical drag coefficients, which depend only on the Reynolds number, should be obtained from these new general formulas if to use the well-known experimental data for the initial turbulence. It is found that the limitations of classical empirical and semiempirical formulas for the hydraulic drag coefficients and their deviation from the experimental data depend on different parameters of initial fluctuations in the flow for different experiments in a wide range of Reynolds numbers. On the basis of these new dependencies, it is possible to explain that the differences between the experimental results for the fixed Reynolds number are caused by the difference in the values of flow fluctuations for each experiment instead of only due to the systematic error in the processing of experiments. Accordingly, the obtained general dependencies for the smooth flat plate and the smooth round tube can serve as the basis for clarifying the results of experiments and the experimental formulas, which used for continuum flows in different devices.     

Long-wavelength dispersion of acoustic waves in transversely inhomogeneous anisotropic plates

Long-wavelength onset of the fundamental branches is described for a free anisotropic plate with arbitrary through-plate variation of material properties. Main attention is given to the flexural branch. Closed-form expressions for the leading-order dispersion coefficient of the velocity and displacement are derived for a generic case and exemplified for the various types of either continuous, or discrete, or periodic inhomogeneity and for the monoclinic symmetry. The relevance of the static averaging is examined in detail. The bounds for the slope of the flexural velocity branch are established. The upper fundamental branches are considered for the case when these are uncoupled inplane and shear horizontal ones.     

Corrections for spatial velocity derivatives in a turbulent shear flow

Possible corrections for measured spatial velocity derivatives have been inferred from a direct numerical simulation database for a fully developed turbulent channel flow. The magnitude of the correction depends much less on the distance from the wall for derivatives in the spanwise direction than for those in the wall-normal direction. Corrections based on local isotropy are better approximations for spanwise derivatives than for wall-normal derivatives.The support of the Australian Research Council is gratefully ackowledged.     

Inelastic behavior of thick-walled cylinders subjected to nonproportionate loading

Inelastic behavior of thick-walled cylinders subjected to nonproportionate loading was studied by the testing of specimens made of C1045 steel and of annealed copper. Several theories were reviewed. A closed-form solution proposed by Mendelson¹² was used to predict external strains for open-end and closed-end thick-walled cylinders. An incremental theory proposed by Chu¹³ was used to provide incremental solutions for open-end thick-walled cylinders, and for cylinders subjected to nonproportionate loading. Test data for open-end and closed-end thick-walled cylinders made of C1045 steel and of annealed copper were in excellent agreement with the incremental theory. Larger values were predicted by use of the closed-form solution for circumferential strains than actual test data for open-end thick-walled cylinders at large depth of yielding. For cylinders subjected to nonproportionate loading, excellent agreement was indicated between the incremental theory and the experiments for the plot of axial load vs. circumferential strain for specimens made of both metals. Agreement between the incremental theory prediction of axial strains for the specimens made of annealed copper and test data is quite satisfactory. Larger values were predicted by the incremental theory for axial strain than experimental data for specimens made of C1045 steel. The error was conservative.