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

Elliptical notches in rectangular beams under pure bending are examined photoelastically. Stress-concentration factors due to a single elliptical notch are obtained for wide ranges of 2a/h andd/h, where 2a, d, andh are the width of notch, depth of notch, and depth of beam, respectively. In particular, the geometries of the optimum elliptical notches producing the least stress concentrations are obtained. Almost the whole elliptical boundary of these notches are stressed to the same peak, which indicates that these notches will probably produce the least stress concentrations among all notches, elliptical or nonelliptical. The graphs herein will enable the designers to find the stress-concentration factors of elliptical notches and to pick out the geometry of the optimum notch which will give the least stress concentration for any given values of 2a/h andd/h. Stresses and the stress-concentration factors at the bottom of the beam opposite the notch are also obtained. These stresses, though smaller in magnitude, are of an opposite sign to the peak stress at the notch. For brittle materials, a smaller tensile stress may be more critical than a large compressive stress; therefore, these stress-concentration factors are also given.     

Plane acoustic wave propagation through a composite of elastic and Kelvin–Voigt viscoelastic material layers

The problem of plane wave propagation through a plane composite layer of thickness h is considered. The composite consists of periodically repeated elastic and Kelvin–Voigt viscoelastic material layers, and all layers are either parallel or perpendicular to the incident wave front. Moreover, it is assumed that the thickness of each separate layer of the composite is much less than the acoustic wave length and the thickness h of the entire composite. We study the problem by using a homogenized model of the composite, which allows us to find the reflection and transmission factors and the variation in the sound intensity level as it propagates though the composite layer of thickness h.     

Transient response of bonded strain gages

Transient-response characteristics of bonded strain gages have been studied by measuring the elastic step wave produced in a steel bar. A quenched-steel bar with a circumferential notch is stretched statically along its axis until it fractures at the notch. Thus a new cross section is suddenly created in a bar under static tension. At this moment, a sharp step wave of zero stress is produced from the new section, and it travels in the bar at the velocity of sound. By measuring this step wave with strain gages, it is shown that the rise time of the gage itself is less than 0.5 μs+0.8L/c, whereL is the gage length andc is the velocity of longitudinal elastic wave in the bar.     

Mach Wave Effect on Laminar-Turbulent Transition in Supersonic Flow over a Flat Plate

The effect of a Mach wave (N wave) on laminar-turbulent transition induced by the first instability mode (Tollmien–Schlichting wave) in the flat-plate boundary layer is investigated on the basis of the numerical solution of Navier–Stokes equations at the freestream Mach number of 2.5. In accordance with the experiment, the N wave is generated by a two-dimensional roughness at the computation domain boundary corresponding to the side wall of the test section of a wind tunnel. It is shown that the disturbance induced by the backward front of the N wave in the boundary layer has no effect on the beginning of transition but displaces downstream the nonlinear stage of the first mode development. The disturbance induced by the forward front of the N wave displaces the beginning of transition upstream.     

Photographic method to sharpen and double isochromatic fringes

A photographic method is described which gives, upon the superposition of ordinary light- and darkfield isochromatic fringe patterns, a new fringe pattern (mixed field). This mixed-field pattern has fringes at the (N/4 and 3N/4 positions. Use of the mixed-field fringe pattern coupled with the oridinary light- and dark-field patterns permit fringes to be read at theN/4,N/2, 3N/4 andN positions, and thus represents a factor of 2 increase in the number of countable fringes. The method is illustrated on two plane-stress examples, a three-dimensional example, and a case employing birefringent coatings.     

Wave reflection at a free interface in an anisotropic pyroelectric medium with nonclassical thermoelasticity

In this paper, the well-established two-dimensional mathematical model for linear pyroelectric materials is employed to investigate the reflection of waves at the boundary between a vacuum and an elastic, transversely isotropic, pyroelectric material. A comparative study between the solutions of (a) classical thermoelasticity, (b) Cattaneo–Lord–Shulman theory and (c) Green–Lindsay theory equations, characterised by none, one and two relaxation times, respectively, is presented. Suitable boundary conditions are considered in order to determine the reflection coefficients when incident elasto–electro–thermal waves impinge the free interface. It is established that, in the quasi-electrostatic approximation, three different classes of waves: (1) two principally elastic waves, namely a quasi-longitudinal Primary (qP) wave and a quasi-transverse Secondary (qS) wave; and (2) a mainly thermal (qT) wave. The observed electrical effects are, on the other hand, a direct consequence of mechanical and thermal phenomena due to pyroelectric coupling. The computed reflection coefficients of plane qP waves are found to depend upon the angle of incidence, the elastic, electric and thermal parameters of the medium, as well as the thermal relaxation times. The special cases of normal and grazing incidence are also derived and discussed. Finally, the reflection coefficients are computed for cadmium selenide observing the influence of (1) the anisotropy of the material, (2) the electrical potential and (3) temperature variations and (4) the thermal relaxation times on the reflection coefficients.     

Characterization of Viscoelasticity Due to Shear Wave Propagation: a Comparison of Existing Methods Based on Computational Simulation and Experimental Data

Shear wave propagation causes microvibrations within a medium; measuring the wave attenuation coefficient, α, and phase velocity, c _s , the medium shear modulus, μ, and shear viscosity, η, are determined based on a viscoelastic model that includes both c _s and α. The present work compares the performances of nine processing methods, based on cross-correlation and quadrature demodulation, used to extract the motion waveform from a sequence of radio-frequency (RF) echo signals from the medium. Kalman filtering determined the amplitude and the phase of the extracted motion waveform. The comparisons were done with regard to computational simulation and experiments with a gel phantom. Estimates obtained for μ and η of the medium considered different conditions for the vibration amplitude and the signal-to-noise ratio (SNR) of the RF echo signals and the waveform extracted by means of single frequency and shear wave dispersion ultrasound vibration (SDUV) methods. According to the simulated results, the cross-correlation-based processing techniques are more precise and accurate in comparison to quadrature demodulation techniques. The results for c _s , α, μ and η of the phantom and those obtained under the same setup conditions for experimental and computational tests agree with each other. Comparing the estimates based on single frequency and SDUV techniques, they presented similar performances at high SNR of the RF echo signal. On the other hand, the former technique prevailed for low SNR.     

Streamwise Vortices and Velocity Streaks in a Locally Drag-Reduced Turbulent Boundary Layer

This work aims to understand the changes associated with the near-wall streaky structures in a turbulent boundary layer (TBL) where the local skin-friction drag is substantially reduced. The Reynolds number is R e _?? = 1000 based on the momentum thickness or R e _τ = 440 based on the friction velocity of the uncontrolled flow. The TBL is perturbed via a local surface oscillation produced by an array of spanwise-aligned piezo-ceramic (PZT) actuators and measurements are made in two orthogonal planes using particle image velocimetry (PIV). Data analyses are conducted using the vortex detection, streaky structure identification, spatial correlation and proper orthogonal decomposition (POD) techniques. It is found that the streaky structures are greatly modified in the near-wall region. Firstly, the near-wall streamwise vortices are increased in number and swirling strength but decreased in size, and are associated with greatly altered velocity correlations. Secondly, the velocity streaks grow in number and strength but contract in width and spacing, exhibiting a regular spatial arrangement. Other aspects of the streaky structures are also characterized; they include the spanwise gradient of the longitudinal fluctuating velocity and both streamwise and spanwise integral length scales. The POD analysis indicates that the turbulent kinetic energy of the streaky structures is reduced. When possible, our results are compared with those obtained by other control techniques such as a spanwise-wall oscillation, a spanwise oscillatory Lorentz force and a transverse traveling wave.     

Effects of the Reynolds number on a scale-similarity model of Lagrangian velocity correlations in isotropic turbulent flows

A scale-similarity model of a two-point two-time Lagrangian velocity correlation(LVC) was originally developed for the relative dispersion of tracer particles in isotropic turbulent flows(HE, G. W., JIN, G. D., and ZHAO, X. Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence. Physical Review E, 80, 066313(2009)). The model can be expressed as a two-point Eulerian space correlation and the dispersion velocity V. The dispersion velocity denotes the rate at which one moving particle departs from another fixed particle. This paper numerically validates the robustness of the scale-similarity model at high Taylor micro-scale Reynolds numbers up to 373, which are much higher than the original values(R_λ = 66, 102). The effect of the Reynolds number on the dispersion velocity in the scale-similarity model is carefully investigated. The results show that the scale-similarity model is more accurate at higher Reynolds numbers because the two-point Lagrangian velocity correlations with different initial spatial separations collapse into a universal form compared with a combination of the initial separation and the temporal separation via the dispersion velocity.Moreover, the dispersion velocity V normalized by the Kolmogorov velocity V_η≡η/τ_η in which η and τ_η are the Kolmogorov space and time scales, respectively, scales with the Reynolds number R_λ as V/V_η∝ R_λ~(1.39) obtained from the numerical data.     

Fracture in Microscale SU-8 Polymer Thin Films

In this work, the fracture of spin coated SU-8 epoxy thin films was investigated under mode I loading using in situ optical experiments on specimens with double edge notched tensile geometry. A method was developed to fabricate 3 μm thick SU-8 films with tapered Chevron type notches using a combination of electron beam and ultra-violet lithography techniques. Subsequently, through speckle patterning under tensile loading, the local deformation fields around the crack tip were extracted using digital image correlation. Since the notches were tapered through the thickness, a crack nucleated from them and grew stably until it spanned the entire thickness before propagating unstably leading to catastrophic failure. As SU-8 underwent brittle fracture with no evidence of a large process zone, the critical energy release rate, J _{I C} was computed from deformation fields, and was found to be 106.6 ± 12.03 J /m ². As the film thickness was small compared to lateral dimensions, assuming plane stress conditions, the critical stress intensity factor was calculated as 0.57 ± 0.03 MPa\(\sqrt {m}\). Furthermore, to assess the validity of the experimental method, a finite element simulation on the exact specimen geometry was conducted with experimentally evaluated far field displacement boundary conditions. The strain fields and J-integral value obtained from the simulation were in good agreement with the experimental results, implying the validity of the in situ experimental method proposed given the challenges of small scale specimens. Furthermore, using fractography and optical imaging it was confirmed that the unstable crack propagation started once the crack front reached full thickness, thereby providing sharp crack at the time of failure, which is necessary for brittle materials for valid fracture toughness experiment. It is expected that the proposed methods of specimen preparation and fracture experiments on microscale polymer thin films can be used on other materials.     

Operator method for solving the plane elasticity problem for a strip with periodic cuts

In the present paper, we use the conformal mapping z/c = ζ?2a sin ζ (a, c?const, ζ = u + iv) of the strip {|v| ≤ v ₀, |u| < ∞} onto the domain D, which is a strip with symmetric periodic cuts. For the domain D, in the orthogonal system of isometric coordinates u, v, we solve the plane elasticity problem. We seek the biharmonic function in the form F = C ψ ₀ + S ψ*₀ + x(C ψ ₁ ? S ψ ₂) + y(C ψ ₂ + S ψ ₁), where C(v) and S(v) are the operator functions described in [1] and ψ ₀(u), …, ψ ₂(u) are the desired functions. The boundary conditions for the function F posed for v = ±v ₀ are equivalent to two operator equations for ψ ₁(u) and ψ ₂(u) and to two ordinary differential equations of first order for ψ ₀(u) and ψ*₀(u) [2]. By finding the functions ψ _j (u) in the form of trigonometric series with indeterminate coefficients and by solving the operator equations, we obtain infinite systems of linear equations for the unknown coefficients. We present an efficient method for solving these systems, which is based on studying stable recursive relations. In the present paper, we give an example of analysis of a specific strip (a = 1/4, v ₀ = 1) loaded on the boundary v = v ₀ by a normal load of intensity p. We find the particular solutions corresponding to the extension of the strip by the longitudinal force X ^∞ and to the transverse and pure bending of the strip due to the transverse force Y ^∞ and the constant moment M ^∞, respectively. We also present the graphs of normal and tangential stresses in the transverse cross-section x = 0 and study the stress concentration effect near the cut bottom.     

Effect of tube diameter and capillary number on platelet margination and near-wall dynamics

I investigate the effect of tube diameter D and red blood cell capillary number Ca (i.e. the ratio of viscous to elastic forces) on platelet margination in blood flow at ≈37 % tube haematocrit. The system is modelled as three-dimensional suspension of deformable red blood cells and nearly rigid platelets using a combination of the lattice-Boltzmann, immersed boundary and finite element methods. Results of simulations during the dynamics before the steady state has been reached show that a non-diffusive radial platelet transport facilitates margination. This non-diffusive effect is important near the edge of the cell-free layer, but only for Ca > 0.2, when red blood cells are tank-treading. I also show that platelet trapping in the cell-free layer is reversible for Ca ≤ 0.2. Margination is essentially independent of Ca only for the smallest investigated tube diameter (D = 10 μm). Once platelets have reached the cell-free layer, they tend to slide rather than tumble. The tumbling rate is essentially independent of Ca but increases with D. Strong confinement suppresses tumbling due to the relatively small cell-free layer thickness at ≈ 37 % tube haematocrit.     

Surface waves under constrained deformation

The possibility of the existence of surface waves in a range of velocities greater than the velocity of transverse waves, but smaller than the velocity of longitudinal waves is shown. It turns out that, in the boundary value problem for an elastic half-space in this velocity range, there are the surface waves whose velocity is constant and equal to \(\sqrt 2 \) b, where b is the velocity of transverse waves. These waves as well as the Rayleigh surface waves have no dispersion. Their velocity is specified only by the elastic constants and density of the medium. It is also shown that the existence of such a velocity is possibly related to the velocity of surface waves that appear as unloading waves under constrained deformation.     

A simple criterion for finite time stability with application to impacted buckling of elastic columns

The existing theories of finite-time stability depend on a prescribed bound on initial disturbances and a prescribed threshold for allowable responses. It remains a challenge to identify the critical value of loading parameter for finite time instability observed in experiments without the need of specifying any prescribed threshold for allowable responses. Based on an energy balance analysis of a simple dynamic system, this paper proposes a general criterion for finite time stability which indicates that finite time stability of a linear dynamic system with constant coefficients during a given time interval [0, t _f ] is guaranteed provided the product of its maximum growth rate (determined by the maximum eigen-root p₁ >0) and the duration t _f does not exceed 2, i.e., p₁t _f <2. The proposed criterion (p₁t _f =2) is applied to several problems of impacted buckling of elastic columns: (i) an elastic column impacted by a striking mass, (ii) longitudinal impact of an elastic column on a rigid wall, and (iii) an elastic column compressed at a constant speed (“Hoff problem”), in which the time-varying axial force is replaced approximately by its average value over the time duration. Comparison of critical parameters predicted by the proposed criterion with available experimental and simulation data shows that the proposed criterion is in robust reasonable agreement with the known data, which suggests that the proposed simple criterion (p₁t _f =2) can be used to estimate critical parameters for finite time stability of dynamic systems governed by linear equations with constant coefficients.     

Effect of Flow Structures on Turbulence Statistics of Taylor-Couette Flow in the Torque Transition State

Direct numerical simulations of Taylor-Couette flow from Re= 8000 to 25000 have been conducted to investigate changes of turbulence statistics in the transition of the Reynolds number dependency of the mean torque near Re= 10000. The velocity fluctuations are decomposed into the contributions of the Taylor vortex and remaining turbulent fluctuations. Significant Reynolds number dependencies of these components are observed in the radial profiles of the Reynolds stress and the transmission of the mean torque. The contributions of Taylor vortex and turbulent components in the net amount of mean torque are evaluated. The Taylor vortex component is overtaken by the turbulent counterpart around Re= 15000 when they are defined as the azimuthally averaged component and the remnants. The results show that the torque transition can be explained by the competition between the contributions of azimuthally averaged Taylor vortex and the remaining turbulent components.     

Dual yielding in capillary suspensions

Rheological measurements were performed to examine the yielding behavior of capillary suspensions prepared by mixing cocoa powder as dispersed phase, vegetable oil as the continuous primary fluid, and water as the secondary fluid. Here, we investigated the yielding behavior of solid-fluid-fluid systems with varying particle volume fraction, ?, spanning the regime from a low volume fraction (? = 0.25) to a highly filled regime (? = 0.65) using dynamic oscillatory measurements. While for ? ≤ 0.4 with a fixed water volume fraction (? _w ) of 0.06 as the secondary fluid, capillary suspensions exhibited a single yield point due to rupturing of aqueous capillary bridges between the particles, while capillary suspensions with ? ≥ 0.45 showed a two-step yielding behavior. On plotting elastic stress (G^′ γ) as a function of applied strain (γ), two distinct peaks, indicating two yield stresses, were observed. Both the yield stresses and storage modulus at low strains were found to increase with ? following a power law dependence. With increasing ? _w (0 – 0.08) at a fixed ? = 0.65, the system shifted to a frustrated, jammed state with particles strongly held together shown by rapidly increasing first and second yield stresses. In particular, the first yield stress was found to increase with ? _w following a power law dependence, while the second yield stress was found to increase exponentially with ? _w . Transient steady shear tests were also performed. The single stress overshoot for ? ≤ 0.4 with ? _w = 0.06 reflected one-step yielding behavior. In contrast, for high ? (≥ 0.45) values with ? _w = 0.06, two stress overshoots were observed in agreement with the two-step yielding behavior shown in the dynamic oscillatory measurements. Experiments on the effect of resting time on microstructure recovery demonstrated that aggregates could reform after resting under quiescent conditions.     

Wall slip and multi-tier yielding in capillary suspensions

Impact of wall slip on the yield stress measurement is examined for capillary suspensions consisting of cocoa powder as the dispersed phase, vegetable oil as the continuous primary fluid, and water as the secondary fluid using smooth and serrated parallel plates. Using dynamic oscillatory measurements, we investigated the yielding behavior of this ternary solid-fluid-fluid system with varying particle volume fraction, ?, from 0.45 to 0.65 and varying water volume fraction, ?_w, from 0.02 to 0.08. Yield stress is defined as the maximum in the elastic stress (G^′γ), which is obtained by plotting the product of elastic modulus (G^′) and strain amplitude (γ) as a function of applied strain amplitude. With serrated plates, which offer minimal slippage, capillary suspensions with ? ≥?0.45 and a fixed ?_w =?0.06 showed a two-step yielding behavior as indicated by two peaks in the plots of elastic stress as a function of strain amplitude. On the other hand with smooth plates, the capillary suspensions showed strong evidence of wall slip as evident by the presence of three distinct peaks and lowered first yield stresses for all ? and ?_w. These results can be interpreted based on the fact that a particle-depleted layer, which is known to be responsible for slip, is present in the vicinity of the smooth surfaces. The slip layer presents itself as an additional “pseudo-microstructure” (characteristic length scale) besides the two microstructures, aqueous bridges and solid particle agglomerates, that may occur in the system. With serrated plates, both the yield stresses (σ₁σ₂) and storage moduli plateau at lower strain (before the first yield point) and at higher strain (before the second yield point) (G\(^{\prime }_{p1}\), G\(^{\prime }_{p2}\)) were found to increase with ? (at a fixed ?_w =?0.06) following power-law dependences. Similarly with increasing ?_w (0.02 – 0.08) at a fixed ? =?0.62, the system behaved as a solid-like material in a jammed state with particles strongly held together as manifested by rapidly increasing σ₁ and σ₂. The usage of smooth surfaces primarily affected σ₁ which was reflected by an approximately 70–90% decrement in the measured σ₁ for all values of ?. By contrast, σ₂ and G\(^{\prime }_{p2}\) were found to be unaffected as shown by close agreement of values obtained using serrated geometry due to vanishing slip layers at higher strain amplitudes.     

Direct Numerical Simulation and Theory of a Wall-Bounded Flow with Zero Skin Friction

We study turbulent plane Couette-Poiseuille (CP) flows in which the conditions (relative wall velocity ΔU _w ≡ 2U _w , pressure gradient dP/dx and viscosity ν) are adjusted to produce zero mean skin friction on one of the walls, denoted by APG for adverse pressure gradient. The other wall, FPG for favorable pressure gradient, provides the friction velocity u _τ , and h is the half-height of the channel. This leads to a one-parameter family of one-dimensional flows of varying Reynolds number Re ≡ U _w h/ν. We apply three codes, and cover three Reynolds numbers stepping by a factor of two each time. The agreement between codes is very good, and the Reynolds-number range is sizable. The theoretical questions revolve around Reynolds-number independence in both the core region (free of local viscous effects) and the two wall regions. The core region follows Townsend’s hypothesis of universal behavior for the velocity and shear stress, when they are normalized with u _τ and h; on the other hand universality is not observed for all the Reynolds stresses, any more than it is in Poiseuille flow or boundary layers. The FPG wall region obeys the classical law of the wall, again for velocity and shear stress. For the APG wall region, Stratford conjectured universal behavior when normalized with the pressure gradient, leading to a square-root law for the velocity. The literature, also covering other flows with zero skin friction, is ambiguous. Our results are very consistent with both of Stratford’s conjectures, suggesting that at least in this idealized flow turbulence theory is successful like it was for the classical logarithmic law of the wall. We appear to know the constants of the law within a 10% bracket. On the other hand, that again does not extend to Reynolds stresses other than the shear stress, but these stresses are passive in the momentum equation.     

Effects of Fuel Lewis Number on Localised Forced Ignition of Turbulent Mixing Layers

The influences of fuel Lewis number Le _F on localised forced ignition of inhomogeneous mixtures are analysed using three-dimensional compressible Direct Numerical Simulations (DNS) of turbulent mixing layers for Le _F  = 0.8, 1.0 and 1.2 and a range of different root-mean-square turbulent velocity fluctuation u′ values. For all Le _F cases a tribrachial flame has been observed in case of successful ignition. However, the lean premixed branch tends to merge with the diffusion flame on the stoichiometric mixture fraction isosurface at later stages of the flame evolution. It has been observed that the maximum values of temperature and reaction rate increase with decreasing Le _F during the period of external energy addition. Moreover, Le _F is found to have a significant effect on the behaviours of mean temperature and fuel reaction rate magnitude conditional on mixture fraction values. It is also found that reaction rate and mixture fraction gradient magnitude \(\vert \nabla \xi \vert \) are negatively correlated at the most reactive region for all values of Le _F explored. The probability of finding high values of \(\vert \nabla \xi \vert \) increases with increasing Le _F . For a given value of u′, the extent of burning decreases with increasing Le _F . A moderate increase in u′ gives rise to an increase in the extent of burning for Le _F  = 0.8 and 1.0, which starts to decrease with further increases in u′. For Le _F  = 1.2, the extent of burning decreases monotonically with increasing u′. The extent of edge flame propagation on the stoichiometric mixture fraction ξ = ξ _st isosurface is characterised by the probability of finding burned gas on this isosurface, which decreases with increasing u′ and Le _F . It has been found that it is easier to obtain self-sustained combustion following localised forced ignition in case of inhomogeneous mixtures than that in the case of homogeneous mixtures with the same energy input, energy deposition duration when the ignition centre is placed at the stoichiometric mixture. The difficultly to sustain combustion unaided by external energy addition in homogeneous mixture is particularly prevalent in the case of Le _F  = 1.2.     

Measurement of local forcing on a turbulent boundary layer using PIV

An experimental study was carried out to investigate the effect of periodic blowing and suction on a turbulent boundary layer. Particle image velocimetry (PIV) was used to probe the characteristics of the flow. Local forcing was introduced to the boundary layer via a sinusoidally-oscillating jet issuing from a thin spanwise slot. Three forcing frequencies (f⁺=0.44, 0.66 and 0.88) with a fixed forcing amplitude (A⁺=0.6) were employed at Re _θ =690. The effect of three different forcing angles (α=60°, 90° and l20°) was investigated under a fixed forcing frequency (f⁺=0.088). The PIV results showed that the wall-region velocity decreases on imposition of the local forcing. Inspection of the phase-averaged velocity profiles revealed that spanwise large-scale vortices are generated downstream of the slot and persist farther downstream. The highest reduction in skin friction was achieved at the highest forcing frequency (f⁺=0.088) and a forcing angle of α=120°. The spatial fraction of the vortices was examined to analyze the skin friction reduction.