Analyzing the influence of compressibility on the rapid pressure–strain rate correlation in turbulent shear flows

The influence of compressibility on the rapid pressure–strain rate tensor is investigated using the Green’s function for the wave equation governing pressure fluctuations in compressible homogeneous shear flow. The solution for the Green’s function is obtained as a combination of parabolic cylinder functions; it is oscillatory with monotonically increasing frequency and decreasing amplitude at large times, and anisotropic in wave-vector space. The Green’s function depends explicitly on the turbulent Mach number M _t , given by the root mean square turbulent velocity fluctuations divided by the speed of sound, and the gradient Mach number M _g , which is the mean shear rate times the transverse integral scale of the turbulence divided by the speed of sound. Assuming a form for the temporal decorrelation of velocity fluctuations brought about by the turbulence, the rapid pressure–strain rate tensor is expressed exactly in terms of the energy (or Reynolds stress) spectrum tensor and the time integral of the Green’s function times a decaying exponential. A model for the energy spectrum tensor linear in Reynolds stress anisotropies and in mean shear is assumed for closure. The expression for the rapid pressure–strain correlation is evaluated using parameters applicable to a mixing layer and a boundary layer. It is found that for the same range of M _t there is a large reduction of the pressure–strain correlation in the mixing layer but not in the boundary layer. Implications for compressible turbulence modeling are also explored.      

Heat transfer in the flow of a viscoelastic fluid over a stretching surface

An analysis is made of heat transfer in the boundary layer of a viscoelastic fluid flowing over a stretching surface. The velocity of the surface varies linearly with the distance x from a fixed point and the surface is held at a uniform temperature T _w higher than the temperature T _∞ of the ambient fluid. An exact analytical solution for the temperature distribution is found by solving the energy equation after taking into account strain energy stored in the fluid (due to its elastic property) and viscous dissipation. It is shown that the temperature profiles are nonsimilar in marked contrast with the case when these profiles are found to be similar in the absence of viscous dissipation and strain energy. It is also found that temperature at a point increases due to the combined influence of these two effects in comparison with its corresponding value in the absence of these two effects. A novel result of this analysis is that for small values of x, heat flows from the surface to the fluid while for moderate and large values of x, heat flows from the fluid to the surface even when T _w >T _∞. Temperature distribution and the surface heat flux are determined for various values of the Prandtl number P, the elastic parameter K ₁ and the viscous dissipation parameter a. Numerical solutions are also obtained through a fourth-order accurate compact finite difference scheme. Received on 14 October 1997     

Comparison of energy release rate and minimum strain energy density for potential failure of composite with bi-material interface

An analytical method is developed to describe the fields of stress and displacement in a bi-material strip specimen with an edge interfacial crack. All of the basic governing equations, boundary conditions on crack surfaces and conditions of continuity along the interface are satisfied by the eigenfunction expansion method. The other boundary conditions are satisfied by the generalized variational principle. The stress intensity factors are calculated for determining the energy release rate and minimum strain energy density factor S_min that is used the strain energy density criterion. Problems with oscillatory singularity and contact zone are discussed. Not only the effects of bi-material modulus ratio, thickness ratio, Poisson's ratio and crack length to S_min, but also the influences of bi-material modulus ratio, thickness ratio to phase angle are presented. Among these parameters, particular situations where S_min become jeopardously high and lead to failure are discussed.     

Transient growth in the asymptotic suction boundary layer

In the present experimental setup, the transient disturbance growth in a spatially invariant boundary layer flow, i.e., the asymptotic suction boundary layer (ASBL), has been investigated. The choice of the ASBL brings along several advantages compared with an ordinary spatially growing boundary layer. A unique feature of the ASBL is that the Reynolds number (Re) can be varied without changing the boundary layer thickness, which in turn allows for parameter variations not possible to carry out in traditional boundary layer flows. A spanwise array of discrete surface roughness elements was mounted on the surface to trigger modes with different spanwise wavenumbers (β). It is concluded that for each mode there exists a threshold roughness Reynolds number (Re _k ), below which no significant transient growth is present. The experimental data suggests that this threshold Re _k is both a function of β and Re. An interesting result is that the energy growth curves respond differently to a change in Re _k when caused by a change in roughness height k, implying that Re remains constant, compared with a change in the free-stream velocity U_￥U_\infty, which also affects the Re. The scaling of the energy growth curves both in level and the downstream direction is treated and appropriate scalings are found. The result shows a complex non-linear receptivity mechanism. Optimal perturbation theory, which has failed to predict the energy evolution in growing boundary layers, is tested for the ASBL and shows that it may satisfactorily predict the evolution of all transiently growing modes that are triggered by the roughness elements.     

A bound on the strain energy for the traction problem in finite elasticity with localized non-zero surface data

For the boundary value problem in finite elasticity in which nonzero tractions are given on a connected subdomain of the boundary, the rest of the boundary is stress-free, and there are no body forces, a bound is obtained for the strain energy in terms of the L ₂ integral norm of the surface tractions with the constant involved depending only upon and the material constants.The result is obtained in the context of finite elasticity under the assumptions that the unstressed body occupies a convex domain and the displacement gradients are sufficiently small. In the context of the linear theory, the same result is obtained without these assumptions.

---

Zusammenfassung Wir betrachten ein Randwertproblem in der nichtlinearen Elastizität in dem nur ein zusammenhängendes Teilgebiet der Randfläche belastet ist, sonst aber keine Randbelastung oder Körperkräfte vorhanden sind. Eine Schranke für die Verzerrungsenergie wird mittels der L ₂ Integralnorme der Randbelastung hergeleitet, wobei die auftretende Konstante nur von dem Teilgebiet und von den Eigenschaften des Materials abhängig ist.Das Ergebnis gilt für die nichtlinearen Elastizität, unter den Vorraussetzungen dass das unbelastete Material ein konvexes Gebiet besetzt und dass die Verschiebungsgradiente hinreichend klein sind. Das gleiche Ergebnis gilt in der linearen Elastizität ohne diese Vorraussetzungen.

---

---

The second author was a visitor at Georgia Institute of Technology, School of Mathematics, at the time that the revised version was prepared.     

Asymptotic computation of the inhomogeneity energy in a body in an external stress field

We consider a problem on an ellipsoidal inhomogeneity in an infinitely extended homogeneous isotropic elastic medium. The inhomogeneity differs from the ambient body in the elastic moduli (Poisson’s ratio ν and shear modulus μ) and in that it has intrinsic strains. We use the equivalent inclusion method to write out expressions for the Helmholtz and Gibbs free energy of the inhomogeneity as quadratic forms in the intrinsic strains and strains at infinity. The general expressions for the coefficients of these quadratic forms are written out as three rank four tensors characterizing the contribution to the energy by the plastic strain (ɛ _p ²), by the strain at infinity (ɛ ₀²), and (only for the Gibbs energy) by the cross term ɛ ⁰ ɛ _p .     

Mixed mode fatigue crack growth and the strain energy density factor

The strain energy density factor S was first proposed by Sih for the prediction of the critical of the load and failure direction under monotonic, mixed mode loading condition. It seems a natural extension to apply the same concept to fatigue crack propagation. However, a close examination of the existing theory indicates that the Strain Energy Density Factor cannot logically account for the phenomena of the R-ratio effect and crack arrest. Thus, modification is necessary before the concept can be applied successfully for the prediction of mixed mode fatigue crack propagation.Based on the concept of hysteresis energy dissipation, an effective strain energy density factor range, ΔS_p,eff, is proposed for the correlation of fatigue crack growth data. ΔS_p,eff is consistent with the concept of crack closure. Experimental investigation indicates that it could predict the crack growth rates and trajectories.     

Non-self-similar crack growth in elastic-plastic finite thickness plate

This work is concerned with non-self-similar crack growth in medium strength metal plates while the loading step, plate thickness and material properties are altered. The three-dimensional elastic-plastic finite element stress analysis is combined with the strain energy density criterion for modeling the material damage process from crack initiation to final global instability including the intervening stage of slow crack growth. Both inelastic deformation and crack growth are accounted for each increment of loading such that the redistribution of stresses and strains are made for each new crack profile. Numerical results are obtained for the center cracked plate configuration under uniform extension with twenty-seven (27) different combinations of specimen thickness, loading step and material type. The fracture toughness S_c being related to K_1c for three different materials are predicted analytically from the corresponding uniaxial tensile test data. Effective strain energy density factor and half crack length are defined so that the results can be compared with their two-dimensional counterparts. Crack growth resistance curves (R-curves) are constructed by plotting as a function of . The condition is found to prevail during slow crack growth. Translation and/or rotation of the lines can yield results other than those calculated and serve a useful purpose for scaling component size and test time. The minimum thickness requirement for the ASTM valid K_1c test is also discussed in connection with predictions based on the strain energy density criterion. The corresponding K_1c for smaller specimens that exhibit moderate ductility and nonlinearity can also be obtained analytically. In such cases, the influence of loading step can be significant and should not be neglected. Notwithstanding the shortcomings of the theory of plasticity, the qualitative features of non-self-similar crack growth are predicted by the strain energy density criterion. Any refinements on the analytical modeling of the material damage process would only affect the results qualitatively, a subject that is left for future investigation.     

Damage accumulation and crack growth in bilinear materials with softening: Application of strain energy density theory

A pseudo-elastic damage-accumulation model is developed by application of the strain energy density theory. The three-point bending specimen is analyzed to illustrate the crack growth characteristics according to a linear elastic softening constitutive law that is typical of concrete materials. Damage accumulation is accounted for by the decrease of elastic modulus and fracture toughness. Both of these effects are assessed by means of the strain energy density functions in the elements around a slowly moving crack. The rate of change of the strain energy density factor S with crack growth as expressed by the relation dS/da = constant is shown to describe the failure behavior of concrete. Results are obtained for different loading steps that yield different slopes of lines in an S versus a (crack length) plot. The lines rotate about the common intersect in an anti-clockwise direction as the load steps are increased. The intersect shifts upward according to increase in the specimen size. In this way, the combined interaction of material properties, load steps and specimen geometry and size are easily analyzed in terms of the failure mode or behavior that can change from the very brittle to the ductile involving stable crack growth. An upper limit on specimen or structural size is established beyond which stable crack growth ceases to occur and failure corresponds to unstable crack propagation or catastrophic fracture. The parameters that control the failure mode are the threshold values of the strain energy density function (dW/dV)_c and the strain energy density factor S_c.     

A steady line heat source in a decagonal quasicrystalline half-space

The problem of a steady line heat source within a decagonal quasicrystalline half-space, in which the boundary is traction-free, is investigated in detail by applying the complex variable technique. The line heat source is infinitely long in the period direction, then the plane strain condition prevails. The four complex stress functions are derived. It is found that the determination of the complex stress functions is independent of the thermal boundary conditions. The stress fields induced by the heat source are explicitly given based on the complex stress functions. Furthermore, we consider, both qualitatively and quantitatively, the distribution of σ₁₁ on the boundary y = 0 of the half-space to demonstrate the influence of the phonon–phason coupling constant R on the stress field.     

Strain energy density criteria for dynamic fracture and dynamic crack branching

Dynamic extension of Sih's fracture criterion based on strain energy density factor, r_c (dW/dV), is used to analyze dynamic crack propagation and branching. Influence of the nonsingular components, which are known as the higher order terms (HOT) in the crack tip stress field, on the strain energy density distribution at a critical distance surrounding the crack tip moving at constant crack velocity is examined. This r_c (dW/dV) fracture criterion is then used to analyze available dynamic photoelastic results of crack branching and of engineering materials.     

Two-dimensional problem of anisotropic elastic body with a hyperbolic boundary

The general and simplified formula for anisotropic medium with a hyperbolic boundary subjected to pure bending M_o is provided in this paper. The stress and strain fields in medium are obtained. Based upon the above results, we analyse the hoop stress along the hyperbolic curve and the stress distributions on the planex₂=0 for aluminium (cubic crystal). When the boundary curve degenerates into an external crack three kinds of stress intensity factors (k₁, k₂, k₃) are obtained, and it is easily found that the first stress intensity factor k₁ is independent of the material constants.     

Integral bounds on the strain energy for the traction problem in finite elasticity

For the traction boundary value problem in finite elasticity, a bound is obtained for the total strain energy in terms of the L₂ integral norms of the surface tractions and body forces, under the assumptions that the unstressed body occupies a convex domain and the displacement gradients are sufficiently small.This is an extension of known results in linear (infinitesimal) elasticity into finite elasticity.     

The flow of a non-Newtonian liquid near a rotating disk

Synopsis The flow of non-Newtonian liquid near a rotating disk has been discussed by using second order stress strain velocity relations of classical hydrodynamics. It is found that the effect of cross-viscosity depends on a non-dimensional number R _c. The boundary layer thickness decreases and the dimensionless moment coefficient increases with the increase of R _c.     

On thermodynamic potentials in linear thermoelasticity

The four thermodynamic potentials, the internal energy u=u(ε_ij,s), the Helmholtz free energy f=f(ε_ij,T), the Gibbs energy g=g(σ_ij,T) and the enthalpy h=h(σ_ij,s) are derived, independently of each other, by using the Duhamel–Neumann extension of Hooke's law and an assumed linear dependence of the specific heat on temperature. A systematic procedure is then presented to express all thermodynamic potentials in terms of four possible pairs of independent state variables. This procedure circumvents a tedious transition from one potential to another, based on the formal change of variables, and inversions of the stress–strain and entropy–temperature relations. The general results are applied to uniaxial loading paths under isothermal, adiabatic, constant stress, and constant strain conditions. An interplay of adiabatic and isothermal elastic constants in the expressions for exchanged heat along certain thermodynamic paths is indicated.     

Prediction of mixed-mode fracture based on an elasto–plastic energy balance criterion

The maximum energy release rate criterion, i.e., G _max criterion, is commonly used for crack propagation analysis. This fracture criterion is based on the elastic macroscopic strength of materials. In the present investigation, however, the G _max criterion has been modified in order to accommodate the consideration of plastic strain energy. This modified criterion is extended to study the fatigue crack growth characteristics of mixed-mode cracks. To predict crack propagation due to fatigue loads, a new elasto–plastic energy model is presented. This new model includes the effects of material properties such as strain hardening exponent n, yield strength σ _y , and fracture toughness and stress intensity factor ranges. The results obtained are compared with those obtained using the commonly employed crack growth law and the experimental data.     

Two-scale micromechanics-based probabilistic modeling of domain switching in ferroelectric ceramics

A two-scale micromechanics model is developed in this paper to analyze domain switching in ferroelectric ceramics, using a probabilistic domain switching criterion based on energetic analysis. The microstructure of ferroelectric ceramics at two distinct length scales, domains and grains, has been carefully analyzed. The interaction at domain level is accounted for by energy minimization theory, while the fluctuation at grain level is analyzed using ellipsoidal two-point correlation function. The model has been implemented by Monte Carlo method, and applied to simulate the electric poling and mechanical depoling of Pb(Zr_xTi_1-x)O₃ (PZT) ceramics across morphotropic phase boundary (MPB). The drastically different switching characteristics of PZT ceramics across MPB has been captured, and good agreement with experiments has been observed. The effects of the transformation strains and spontaneous polarizations are highlighted, confirming the proposition of Li et al. [2005. Domain switching in polycrystalline ferroelectric ceramics. Nature Materials 4, 776–781] that the strain compatibility plays a dominant role in domain switching in ferroelectric ceramics.     

Influence of indenter geometry on half-plane with edge crack subjected to fretting condition

An edge crack is analyzed to study fretting failure. A flat punch with rounded corners and a half-plane are regarded as an indenter and a substrate, respectively. Plane strain condition is considered. Contact shear traction in the case of partial slip is evaluated numerically. It is assumed that an initial crack is extended to the point of minimum strain energy density in the half-plane from the trailing edge of contact. Dislocation density function method is used to evaluate K_I and K_II. The variations of K_I and K_II during crack growth are examined in the case of indentation by a punch with different ratio of the flat region (l) to the punch width (L). Sih's minimum strain energy density theory [1] is also applied to predict the propagation direction of the initial crack. The direction evaluated is similar to that found in the experiment. Stress intensity factor ranges (ΔK_I and ΔK_II) are examined during cyclic shear on the contact. For the design of contacting bodies, a suggestible geometry of punch for alleviating cracking failure is studied.     

On the Stress Tensor Jump at the Interface between a Polyatomic Gas and a Condensed Phase

For the components of the hydrodynamic stress tensor the boundary condition at the interface between a polyatomic gas and a condensed phase is obtained. The boundary-value problem is solved within the framework of the previously proposed kinetic model by the method of semi-spatial moments, taking into account the rotational degrees of freedom of the gas molecules. The gas-kinetic coefficient C _p entering into the boundary condition for the stress tensor components depends on the accommodation coefficient of the tangential momentum, q, the accommodation coefficients of the translational, _t , and rotational, _r , energy components, and the Prandtl number. This coefficient is calculated for several polyatomic gases.     

Analytical solutions for finite cylinders compressed between platens and the strain effect on the valence-band structure of Si1−xGex alloy

This paper presents an analytical solution for inhomogeneous strain and stress distributions within finite circular cylinders of Si_1−xGe_x alloy under compression test with end friction. The method follows Lekhnitskii’s stress function approach, but a new expression for the stress function is proposed so that all of the governing equations and boundary conditions are satisfied exactly. Numerical results show that the axial, radial, circumferential and shear strains are all inhomogeneous within finite cylinders, and local strain concentrations near two end surfaces were usually developed as long as friction exists between end surfaces and loading platens. Moreover, by using envelope-function method, the effect of strain on the valence-band structure of Si_1−xGe_x alloy is also studied. It was found that strain can induce band splitting, alteration of the shape of constant energy surfaces of the heavy-hole and the light-hole bands of Si_1−xGe_x alloy.